source: branches/eml_v9/sample/costs/sample_flash_cost.mso @ 963

Last change on this file since 963 was 650, checked in by Argimiro Resende Secchi, 14 years ago

Fixing convergence problems!

File size: 2.6 KB
Line 
1#*-------------------------------------------------------------------
2* EMSO Model Library (EML) Copyright (C) 2004 - 2007 ALSOC.
3*
4* This LIBRARY is free software; you can distribute it and/or modify
5* it under the therms of the ALSOC FREE LICENSE as available at
6* http://www.enq.ufrgs.br/alsoc.
7*
8* EMSO Copyright (C) 2004 - 2007 ALSOC, original code
9* from http://www.rps.eng.br Copyright (C) 2002-2004.
10* All rights reserved.
11*
12* EMSO is distributed under the therms of the ALSOC LICENSE as
13* available at http://www.enq.ufrgs.br/alsoc.
14*
15*--------------------------------------------------------------------
16* Sample file for model flash and flashSteady
17*--------------------------------------------------------------------
18*
19* This sample file needs VRTherm (www.vrtech.com.br) to run.
20*
21*----------------------------------------------------------------------
22* Author: Núbia do Carmo Ferreira
23* $Id: sample_flash.mso 247 2007-04-24 13:44:18Z rafael $
24*--------------------------------------------------------------------*#
25using "costs/flash_cost";
26
27FlowSheet sample_flash_cost
28       
29        PARAMETERS
30        PP as Plugin(Brief="Physical Properties",
31                Type="PP",
32                Components = ["1,3-butadiene", "isobutene", "n-pentane",
33                        "1-pentene", "1-hexene", "benzene"],
34                LiquidModel = "PR",
35                VapourModel = "PR"
36        );
37        NComp as Integer;
38
39               
40        SET
41        NComp = PP.NumberOfComponents;
42       
43        DEVICES
44        fl as flash_cost;
45        s1 as source;
46        Q as energy_source;
47       
48        CONNECTIONS
49        s1.Outlet to fl.Inlet;
50        Q.OutletQ to fl.InletQ;
51       
52        EQUATIONS
53        fl.OutletL.F = 400*sqrt(fl.Level/'m') * 'kmol/h';
54       
55        SPECIFY
56        s1.F = 496.3 * 'kmol/h';
57        s1.T = 338 * 'K';
58        s1.P = 507.1 * 'kPa';
59        s1.Composition = [0.2379,0.3082,0.09958,0.1373,0.08872,0.1283];
60
61        fl.OutletV.F = 68.5 * 'kmol/h';
62        Q.OutletQ.Q = 0 * 'kJ/h';
63       
64        SET
65        fl.V = 50 * 'm^3';
66        fl.diameter = 2 * 'm';
67        fl.orientation = "vertical";
68#       fl.orientation = "horizontal";
69       
70#costs 
71        fl.flash_length = 6 * 'm';
72        fl.Div = 2 * 'm';
73        fl.Ts = 20 * 'mm';
74        fl.dens_mass_material = 7830 * 'kg/m^3';
75        fl.Dih = 2 * 'm';
76        fl.Material = "Stainless steell 304";
77       
78        fl.Cost(1,:) = [8.114,-0.16449,0.04333];
79        fl.Cost(2,:) = [8.600,-0.21651,0.04576];
80        fl.Cost(3,:) = [1288.3,0,0];
81        fl.Cost(4,:) = [1017.0,0,0];
82        fl.Cost(5,:) = [0.8116,0,0];
83        fl.Cost(6,:) = [1.7,0,0];#For "Stainless steell 304"
84
85        INITIAL
86        fl.OutletL.T = 338 * 'K';
87        fl.Level = 0.4 * 'm';
88
89        fl.OutletL.z(1) = 0.1;
90        fl.OutletL.z(2) = 0.1;
91        fl.OutletL.z(3) = 0.1;
92        fl.OutletL.z(4) = 0.1;
93        fl.OutletL.z(5) = 0.1;
94       
95        OPTIONS
96        Dynamic = true;
97        TimeStep = 0.1;
98        TimeEnd = 5;
99        TimeUnit = 'h';
100end
Note: See TracBrowser for help on using the repository browser.