source: branches/gui/sample/sensitivity/sense_flash.mso @ 804

Last change on this file since 804 was 804, checked in by edsoncv, 13 years ago

Added Case study (dynamic and Steady State case) examples.
Added Sensitivity (only Steady State).

File size: 3.3 KB
RevLine 
[804]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 633 2008-09-24 00:26:30Z bicca $
24*--------------------------------------------------------------------*#
25using "stage_separators/flash";
26
27# exemplo baseado nos dados do artigo do Gani 1985.
28FlowSheet flash_Test
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        VARIABLES
40        Q       as energy_source (Brief="Heat supplied");
41       
42        SET
43        NComp = PP.NumberOfComponents;
44       
45        DEVICES
46        fl as flash;
47        s1 as source;
48       
49        CONNECTIONS
50        s1.Outlet to fl.Inlet;
51        Q.OutletQ to fl.InletQ;
52       
53        EQUATIONS
54        fl.OutletL.F = 400*sqrt(fl.Level/'m') * 'kmol/h';
55       
56        SPECIFY
57        s1.F = 496.3 * 'kmol/h';
58        s1.T = 338 * 'K';
59        s1.P = 507.1 * 'kPa';
60        s1.Composition = [0.2379,0.3082,0.09958,0.1373,0.08872,0.1283];
61
62        fl.OutletV.F = 68.5 * 'kmol/h';
63        Q.OutletQ.Q = 0 * 'kJ/h';
64       
65        SET
66        fl.V = 50 * 'm^3';
67        fl.diameter = 2 * 'm';
68        fl.orientation = "vertical";
69#       fl.orientation = "horizontal";
70
71        INITIAL
72        fl.OutletL.T = 338 * 'K';
73        fl.Level = 0.4 * 'm';
74
75        fl.OutletL.z(1) = 0.1;
76        fl.OutletL.z(2) = 0.1;
77        fl.OutletL.z(3) = 0.1;
78        fl.OutletL.z(4) = 0.1;
79        fl.OutletL.z(5) = 0.1;
80       
81        OPTIONS
82        Dynamic = true;
83        TimeStep = 0.1;
84        TimeEnd = 5;
85        TimeUnit = 'h';
86end
87
88# exemplo baseado nos dados do artigo do Gani 1985.
89FlowSheet flashSteady_Test
90        PARAMETERS
91        PP as Plugin(Brief="Physical Properties",
92                Type="PP",
93                Components = ["1,3-butadiene", "isobutene", "n-pentane",
94                        "1-pentene", "1-hexene", "benzene"],
95                LiquidModel = "PR",
96                VapourModel = "PR"
97        );
98        NComp as Integer;
99
100        VARIABLES
101        Q       as energy_source (Brief="Heat supplied");
102       
103        SET
104        NComp = PP.NumberOfComponents;
105       
106        DEVICES
107        fl as flash_steady;
108        s1 as source;
109       
110        CONNECTIONS
111        s1.Outlet to fl.Inlet;
112        Q.OutletQ to fl.InletQ;
113       
114        SPECIFY
115        s1.F = 496.3 * 'kmol/h';
116        s1.T = 338 * 'K';
117        s1.P = 507.1 * 'kPa';
118        s1.Composition = [0.2379,0.3082,0.09959,0.1373,0.08872,0.1283];
119       
120        fl.OutletL.P = 2.5 * 'atm';
121
122        #Q.OutletQ.Q = 0 * 'kJ/h';
123        fl.OutletL.T = 315.06 * 'K';
124       
125        OPTIONS
126        Dynamic = false;       
127end
128
129
130Sensitivity Sense_Flash as flashSteady_Test
131
132        VARY
133
134        s1.F 'kmol/h';
135
136        RESPONSE
137
138        fl.OutletV.F'kmol/h';
139
140end
141
142# Case study to check the responses
143
144CaseStudy caseFlash as flashSteady_Test
145
146        VARY
147
148        s1.F = [496.3, 496.31]*'kmol/h';
149
150        RESPONSE
151
152
153        fl.OutletV.F'kmol/h';
154
155end
Note: See TracBrowser for help on using the repository browser.