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

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

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