source: branches/newlanguage/sample/stage_separators/sample_flash.mso @ 117

Last change on this file since 117 was 117, checked in by Rafael de Pelegrini Soares, 16 years ago

Updated streams and flash models to the new language being proposed

  • Property svn:eol-style set to native
  • Property svn:keywords set to Id
File size: 2.9 KB
RevLine 
[86]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: Paula B. Staudt
23* $Id: sample_flash.mso 117 2007-01-15 18:43:56Z rafael $
24*--------------------------------------------------------------------*#
[1]25using "stage_separators/flash";
26
27# exemplo baseado nos dados do artigo do Gani 1985.
28FlowSheet flash_Test
29        PARAMETERS
[117]30        PP as Plugin(Brief="Physical Properties",File="vrpp");
31        NComp as Integer;
[1]32
33        VARIABLES
34        Q       as heat_rate (Brief="Heat supplied");
35       
36        SET
37        PP.Components = ["1,3-butadiene", "isobutene", "n-pentane", "1-pentene", "1-hexene", "benzene"];
38        PP.LiquidModel = "PR";
39        PP.VapourModel = "PR";
40        NComp = PP.NumberOfComponents;
41       
42        DEVICES
43        fl as flash;
[117]44        s1 as source;
[1]45       
46        CONNECTIONS
47        s1 to fl.Inlet;
48        Q to fl.Q;
49       
[54]50        EQUATIONS
51        fl.OutletL.F = 400*sqrt(fl.Level/"m") * "kmol/h";
52       
[1]53        SPECIFY
[117]54        s1.Outlet.F = 496.3 * "kmol/h";
55        s1.Outlet.T = 338 * "K";
56        s1.Outlet.P = 507.1 * "kPa";
[35]57        #s1.v = 0.1380;
[117]58        s1.Outlet.z = [0.2379,0.3082,0.09958,0.1373,0.08872,0.1283];
[1]59
60        fl.OutletV.F = 68.5 * "kmol/h";
61        Q = 1 * "kJ/h";
62       
63        SET
64        fl.V = 2000 * "m^3";
65        fl.Across = 1 * "m^2";
66       
67        INITIAL
68        fl.OutletL.T = 330 *"K";
69        fl.Level = 1 * "m";
70        fl.OutletL.z(1) = 0.1;
71        fl.OutletL.z(2) = 0.1;
72        fl.OutletL.z(3) = 0.1;
73        fl.OutletL.z(4) = 0.1;
74        fl.OutletL.z(5) = 0.1;
75       
76        OPTIONS
[117]77        RelativeAccuracy = 1e-7;
78        AbsoluteAccuracy = 1e-9;
79        Dynamic = true;
80        TimeStep = 0.1;
81        TimeEnd = 20;
82        TimeUnit = "h";
[1]83end
84
85# exemplo baseado nos dados do artigo do Gani 1985.
86FlowSheet flashSteady_Test
87        PARAMETERS
[117]88        PP as Plugin(Brief="Physical Properties",File="vrpp");
89        NComp as Integer;
[1]90
91        VARIABLES
92        Q       as heat_rate (Brief="Heat supplied");
93       
94        SET
95        PP.Components = ["1,3-butadiene", "isobutene", "n-pentane", "1-pentene", "1-hexene", "benzene"];
96        PP.LiquidModel = "PR";
97        PP.VapourModel = "PR";
98        NComp = PP.NumberOfComponents;
99       
100        DEVICES
[117]101        fl as flash_steady;
102        s1 as source;
[1]103       
104        CONNECTIONS
105        s1 to fl.Inlet;
106        Q to fl.Q;
107       
108        SPECIFY
[117]109        s1.Outlet.F = 496.3 * "kmol/h";
110        s1.Outlet.T = 338 * "K";
111        s1.Outlet.P = 507.1 * "kPa";
112        #s1.Outlet.v = 0.1380;
113        s1.Outlet.z = [0.2379,0.3082,0.09959,0.1373,0.08872,0.1283];
[1]114       
115        fl.OutletL.P = 2.5 * "atm";
116
117        #Q = 0 * "kJ/h";
118        fl.OutletL.T = 315.06 * "K";
119       
120        OPTIONS
[117]121        Dynamic = false;
[1]122end
Note: See TracBrowser for help on using the repository browser.