source: mso/sample/stage_separators/sample_column.mso @ 46

Last change on this file since 46 was 46, checked in by Paula Bettio Staudt, 16 years ago

Fixed some sample files

  • Property svn:eol-style set to native
  • Property svn:keywords set to Id
File size: 5.0 KB
Line 
1using "stage_separators/column";
2
3# column section with 2 trays
4FlowSheet SectionColumn_Test_with2tray
5        PARAMETERS
6        PP      as CalcObject(Brief="Physical Properties",File="vrpp");
7        NComp   as Integer;
8
9        SET
10        PP.Components = [ "isobutane", "n-pentane"];
11        PP.LiquidModel = "PR";
12        PP.VapourModel = "PR";
13        NComp = PP.NumberOfComponents;
14       
15        DEVICES
16        sec as Section_Column;
17        feed as stream_therm;
18        reb as stream_therm;
19        cond as stream_therm;
20        zero as stream;
21       
22        CONNECTIONS
23        feed to sec.trays(2).Inlet;
24        zero to sec.trays(1).Inlet;
25        reb to sec.trays(2).InletV;
26        cond to sec.trays(1).InletL;
27       
28        SPECIFY
29        feed.F = 113.4 * "kmol/h";
30        feed.T = 291 * "K";
31        feed.P = 168.3 * "kPa";
32        feed.z = [0.5, 0.5];
33        feed.v = 0;
34
35        zero.F = 0 * "kmol/h";
36        zero.T = 300 * "K";
37        zero.P = 1 * "atm";
38        zero.z = [0.5, 0.5];
39        zero.v = 0;
40        zero.h = 0 * "J/mol";
41       
42        cond.F = 68 * "kmol/h";
43        cond.P = 150 * "kPa";
44        cond.T = 281.75 * "K";
45        cond.v = 0.0;
46        cond.z = [0.6664, 0.3336];
47
48        reb.F = 153 * "kmol/h";
49        reb.P = 185 * "kPa";
50        reb.T = 328.12 * "K";
51        reb.z = [0.001848, 0.9982];
52        reb.v = 1.0;
53       
54        sec.trays.Emv = 1;
55
56        SET
57        sec.NTrays = 2;
58        #COLUMN
59        sec.trays.V = 4 * "ft^3";
60        sec.trays.Ah = 0.394 * "ft^2";
61        sec.trays.lw = 20.94 * "in";
62        sec.trays.hw = 0.125 * "ft";
63        sec.trays.Q = 0 * "kW";
64        sec.trays.beta = 0.6;
65        sec.trays.alfa = 4;
66        sec.trays.Ap = 3.94 * "ft^2";
67
68        INITIAL
69        sec.trays.OutletL.T = 290 *"K";
70        sec.trays.Level = 0.9 * sec.trays.hw;
71        sec.trays.OutletL.z(1) = 0.5;
72       
73        OPTIONS
74        #relativeAccuracy = 1e-3;
75        NLASolver = "sundials";
76        time = [0:10:1000];
77end
78
79# column section with 8 trays
80FlowSheet SectionColumn_Test_with8tray
81        PARAMETERS
82        PP      as CalcObject(Brief="Physical Properties",File="vrpp");
83        NComp   as Integer;
84
85        SET
86        PP.Components = [ "isobutane", "n-pentane"];
87        PP.LiquidModel = "PR";
88        PP.VapourModel = "PR";
89        PP.Derivatives = 0;
90        NComp = PP.NumberOfComponents;
91       
92        DEVICES
93        sec as Section_Column;
94        feed as stream_therm;
95        reb as stream_therm;
96        cond as stream_therm;
97        zero as stream;
98       
99        CONNECTIONS
100        feed to sec.trays(5).Inlet;
101       
102        zero to sec.trays([1:4]).Inlet;
103        zero to sec.trays([6:8]).Inlet;
104       
105        reb to sec.trays(8).InletV;
106        cond to sec.trays(1).InletL;
107       
108        SPECIFY
109        feed.F = 113.4 * "kmol/h";
110        feed.T = 291 * "K";
111        feed.P = 168.3 * "kPa";
112        feed.z = [0.5, 0.5];
113        feed.v = 0;
114
115        zero.F = 0 * "kmol/h";
116        zero.T = 300 * "K";
117        zero.P = 1 * "atm";
118        zero.z = [0.5, 0.5];
119        zero.v = 0;
120        zero.h = 0 * "J/mol";
121       
122        cond.F = 68 * "kmol/h";
123        cond.P = 150 * "kPa";
124        cond.T = 281.75 * "K";
125        cond.v = 0.0;
126        cond.z = [0.6664, 0.3336];
127
128        reb.F = 153 * "kmol/h";
129        reb.P = 185 * "kPa";
130        reb.T = 328.12 * "K";
131        reb.z = [0.001848, 0.9982];
132        reb.v = 1.0;
133       
134        sec.trays.Emv = 1;
135
136        SET
137        sec.NTrays = 8;
138        #COLUMN
139        sec.trays.V = 4 * "ft^3";
140        sec.trays.Ah = 0.394 * "ft^2";
141        sec.trays.lw = 20.94 * "in";
142        sec.trays.hw = 0.125 * "ft";
143        sec.trays.Q = 0 * "kW";
144        sec.trays.beta = 0.6;
145        sec.trays.alfa = 4;
146        sec.trays.Ap = 3.94 * "ft^2";
147
148        INITIAL
149        sec.trays.OutletL.T = [290:(330-290)/(sec.NTrays-1):330] *"K";
150        sec.trays.Level = 0.3 * sec.trays.hw;
151        sec.trays.OutletL.z(1) = 0.5;
152       
153        OPTIONS
154        relativeAccuracy = 1e-5;
155        time = [0:1:100];
156end
157
158
159FlowSheet Distillation_kettle_cond_Test
160        PARAMETERS
161        PP      as CalcObject(Brief="Physical Properties",File="vrpp");
162        NComp   as Integer;
163       
164        VARIABLES
165        Qc as heat_rate (Brief="Heat rate removed from condenser");
166        Qr as heat_rate (Brief="Heat rate supplied to reboiler");
167       
168        SET
169        PP.Components = [ "isobutane", "n-pentane", "propylene",
170                "benzene", "isobutene" ];
171        PP.LiquidModel = "PR";
172        PP.VapourModel = "PR";
173        PP.Derivatives = 1;
174        NComp = PP.NumberOfComponents;
175
176        DEVICES
177        col as Distillation_kettle_cond;
178        feed as streamTP;
179        zero as stream;
180       
181        CONNECTIONS
182        feed to col.trays(5).Inlet;
183        zero to col.reb.Inlet;
184        zero to col.trays([1:4]).Inlet;
185        zero to col.trays([6:col.NTrays]).Inlet;
186        Qc to col.cond.Q;
187        Qr to col.reb.Q;
188       
189        SPECIFY
190        feed.F = 113.4 * "kmol/h";
191        feed.T = 291 * "K";
192        feed.P = 168.3 * "kPa";
193        feed.z = 1/NComp;
194       
195        zero.F = 0 * "kmol/h";
196        zero.T = 300 * "K";
197        zero.P = 1 * "atm";
198        zero.z = 1/NComp;
199        zero.v = 0;
200        zero.h = 0 * "J/mol";
201       
202        col.sptop.Outlet2.F = 85 * "kmol/h";
203        col.reb.OutletL.F = 28.4 * "kmol/h";
204        col.sptop.frac = 0.444445;
205        col.cond.OutletV.F = 0 * "kmol/h";
206        Qr = 3.7743e6 * "kJ/h";
207        Qc = -3.71e6 * "kJ/h";
208        col.pump1.dP = 16 * "kPa";
209        col.trays.Emv = 1;
210       
211        SET
212        col.NTrays = 8;
213        col.cond.V = 2 * "m^3";
214        col.cond.Across = 1 * "m^2";
215        col.trays.V = 4 * "ft^3";
216        col.trays.Ah = 0.394 * "ft^2";
217        col.trays.lw = 20.94 * "in";
218        col.trays.hw = 0.125 * "ft";
219        col.trays.Q = 0 * "kW";
220        col.trays.beta = 0.6;
221        col.trays.alfa = 4;
222        col.trays.Ap = 3.94 * "ft^2";
223        col.reb.V = 2 * "m^3";
224        col.reb.Across = 1 * "m^2";
225       
226        INITIAL
227        # condenser
228        col.cond.OutletL.T = 260 *"K";
229        col.cond.Level = 1 * "m";
230        col.cond.OutletL.z([1:4]) = [0.65, 0.05, 0.01, 0.01];
231
232        # reboiler
233        col.reb.OutletL.T = 330 *"K";
234        col.reb.Level = 1 * "m";
235        col.reb.OutletL.z([1:4]) = [0.1, 0.7, 0.01, 0.01];
236
237        # column trays
238        col.trays.OutletL.T = [290:(330-290)/(col.NTrays-1):330] * "K";
239        col.trays.Level = 1.2 * col.trays.hw;
240        col.trays.OutletL.z([1:4]) = [0.5, 0.05, 0.01, 0.01];
241
242        OPTIONS
243        relativeAccuracy = 1e-3;
244        time = [0:0.01:1, 2:50];
245end
Note: See TracBrowser for help on using the repository browser.