source: mso/sample/stage_separators/sample_flash.mso @ 54

Last change on this file since 54 was 54, checked in by Argimiro Resende Secchi, 16 years ago

Fix inlet composition.

  • Property svn:eol-style set to native
  • Property svn:keywords set to Id
File size: 1.8 KB
Line 
1using "stage_separators/flash";
2
3# exemplo baseado nos dados do artigo do Gani 1985.
4FlowSheet flash_Test
5        PARAMETERS
6        PP      as CalcObject(Brief="Physical Properties",File="vrpp");
7        NComp   as Integer;
8
9        VARIABLES
10        Q       as heat_rate (Brief="Heat supplied");
11       
12        SET
13        PP.Components = ["1,3-butadiene", "isobutene", "n-pentane", "1-pentene", "1-hexene", "benzene"];
14        PP.LiquidModel = "PR";
15        PP.VapourModel = "PR";
16        NComp = PP.NumberOfComponents;
17       
18        DEVICES
19        fl as flash;
20        s1 as streamTP;
21       
22        CONNECTIONS
23        s1 to fl.Inlet;
24        Q to fl.Q;
25       
26        EQUATIONS
27        fl.OutletL.F = 400*sqrt(fl.Level/"m") * "kmol/h";
28       
29        SPECIFY
30        s1.F = 496.3 * "kmol/h";
31        s1.T = 338 * "K";
32        s1.P = 507.1 * "kPa";
33        #s1.v = 0.1380;
34        s1.z = [0.2379,0.3082,0.09958,0.1373,0.08872,0.1283];
35
36        fl.OutletV.F = 68.5 * "kmol/h";
37        Q = 1 * "kJ/h";
38       
39        SET
40        fl.V = 2000 * "m^3";
41        fl.Across = 1 * "m^2";
42       
43        INITIAL
44        fl.OutletL.T = 330 *"K";
45        fl.Level = 1 * "m";
46        fl.OutletL.z(1) = 0.1;
47        fl.OutletL.z(2) = 0.1;
48        fl.OutletL.z(3) = 0.1;
49        fl.OutletL.z(4) = 0.1;
50        fl.OutletL.z(5) = 0.1;
51       
52        OPTIONS
53        relativeAccuracy = 1e-7;
54        absoluteAccuracy = 1e-9;
55        #mode="steady";
56        time = [0:0.1:20]*"h";
57end
58
59# exemplo baseado nos dados do artigo do Gani 1985.
60FlowSheet flashSteady_Test
61        PARAMETERS
62        PP      as CalcObject(Brief="Physical Properties",File="vrpp");
63        NComp   as Integer;
64
65        VARIABLES
66        Q       as heat_rate (Brief="Heat supplied");
67       
68        SET
69        PP.Components = ["1,3-butadiene", "isobutene", "n-pentane", "1-pentene", "1-hexene", "benzene"];
70        PP.LiquidModel = "PR";
71        PP.VapourModel = "PR";
72        NComp = PP.NumberOfComponents;
73       
74        DEVICES
75        fl as flash_Steady;
76        s1 as streamTP;
77       
78        CONNECTIONS
79        s1 to fl.Inlet;
80        Q to fl.Q;
81       
82       
83        SPECIFY
84        s1.F = 496.3 * "kmol/h";
85        s1.T = 338 * "K";
86        s1.P = 507.1 * "kPa";
87        #s1.v = 0.1380;
88        s1.z = [0.2379,0.3082,0.09959,0.1373,0.08872,0.1283];
89       
90        fl.OutletL.P = 2.5 * "atm";
91
92        #Q = 0 * "kJ/h";
93        fl.OutletL.T = 315.06 * "K";
94       
95        OPTIONS
96        mode = "steady";
97end
Note: See TracBrowser for help on using the repository browser.