source: trunk/sample/costs/sample_HeatExchangerSimplified_cost.mso @ 488

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

Adding the cost models in the EML.
File size: 7.6 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 models of HeatExchangerSimplified_cost
17*--------------------------------------------------------------------
18*
19* This sample file needs VRTherm DEMO (www.vrtech.com.br) to run.
20*
21*----------------------------------------------------------------------
22* Author: Núbia do Carmo Ferreira
23* $Id: Sample_Simplified.mso 145 2007-01-30 18:20:00Z rafael $
24*--------------------------------------------------------------------*#
25
26using "costs/HeatExchangerSimplified_cost";
27
28
29FlowSheet sample_HeatExchanger_LMTD_cost
30
31DEVICES
32
33        exchanger       as HeatExchanger_LMTD_cost;
34        streamhot       as source;
35        streamcold      as source;
36
37
38CONNECTIONS
39
40        streamhot.Outlet        to exchanger.InletHot;
41        streamcold.Outlet       to exchanger.InletCold;
42       
43PARAMETERS
44        PP                      as Plugin       (Type="PP",LiquidModel    = "RK",VapourModel   = "RK",
45                                                                                Components      = ["water"]);
46        NComp           as Integer;
47
48SET
49
50        NComp                           = PP.NumberOfComponents;
51       
52        exchanger.FlowDirection      = "counter";
53
54        # costs
55        exchanger.Material = "Stainless steel 316";
56       
57        exchanger.Cost(1,:) = [8.202,0.01506,0.06811];
58        exchanger.Cost(2,:) = [-0.9003,0.0906,0];
59        exchanger.Cost(3,:) = [1.35,0,0];
60        exchanger.Cost(4,:) = [-0.7844,0.0830,0];
61        exchanger.Cost(5,:) = [0.8955,0.04981,0];
62        exchanger.Cost(6,:) = [1.2002,0.07140,0];
63        exchanger.Cost(7,:) = [1.4272,0.12088,0];
64        exchanger.Cost(8,:) = [0.8608,0.23296,0];
65       
66
67SPECIFY
68
69        exchanger.Method.Fc                                      = 1;
70       
71        exchanger.Details.A                                      = 139.35*'m^2';
72        exchanger.Details.Ud                                     = 0.75*'W/(m^2*K)';
73        exchanger.Details.Uc                                     = 0.75*'W/(m^2*K)';
74       
75        exchanger.HotSide.PressureDrop.Pdrop  = 0*'atm';
76        exchanger.ColdSide.PressureDrop.Pdrop = 0*'atm';
77       
78        streamhot.Outlet.F                                      = 1000    * 'kmol/h';
79        streamhot.Outlet.T                                      = 450     * 'K';
80        streamhot.Outlet.P                                      = 1.5     * 'atm';
81        streamhot.Outlet.z                                      = [1];
82
83        streamcold.Outlet.F                                     = 500   * 'kmol/h';
84        streamcold.Outlet.P                                     = 1     * 'atm';
85        streamcold.Outlet.T                                     = 315   * 'K';
86    streamcold.Outlet.z                                         = [1];
87       
88       
89OPTIONS
90
91 Dynamic = false;
92
93end
94
95
96FlowSheet sample_HeatExchanger_NTU_cost
97       
98DEVICES
99
100        exchanger       as HeatExchanger_NTU_cost;
101        streamhot       as source;
102        streamcold      as source;
103
104CONNECTIONS
105
106        streamhot.Outlet        to exchanger.InletHot;
107        streamcold.Outlet       to exchanger.InletCold;
108       
109PARAMETERS
110        PP                      as Plugin       (Type="PP",LiquidModel    = "RK",VapourModel   = "RK",
111                                                                                Components      = ["water"]);
112        NComp           as Integer;
113       
114SET
115
116
117        NComp                           = PP.NumberOfComponents;
118       
119        exchanger.FlowDirection      = "counter";
120
121        # costs
122        exchanger.Material = "Stainless steel 316";
123       
124        exchanger.Cost(1,:) = [8.202,0.01506,0.06811];
125        exchanger.Cost(2,:) = [-0.9003,0.0906,0];
126        exchanger.Cost(3,:) = [1.35,0,0];
127        exchanger.Cost(4,:) = [-0.7844,0.0830,0];
128        exchanger.Cost(5,:) = [0.8955,0.04981,0];
129        exchanger.Cost(6,:) = [1.2002,0.07140,0];
130        exchanger.Cost(7,:) = [1.4272,0.12088,0];
131        exchanger.Cost(8,:) = [0.8608,0.23296,0];
132
133SPECIFY
134
135        exchanger.Details.A                                      = 65.031*'m^2';
136        exchanger.Details.Ud                                     = 0.75*'W/(m^2*K)';
137        exchanger.Details.Uc                                     = 0.75*'W/(m^2*K)';
138       
139        exchanger.HotSide.PressureDrop.Pdrop  = 0*'atm';
140        exchanger.ColdSide.PressureDrop.Pdrop = 0*'atm';
141       
142        streamhot.Outlet.F                                      = 1000    * 'kmol/h';
143        streamhot.Outlet.T                                      = 450     * 'K';
144        streamhot.Outlet.P                                      = 1.5     * 'atm';
145        streamhot.Outlet.z                                      = [1];
146
147        streamcold.Outlet.F                                     = 500   * 'kmol/h';
148        streamcold.Outlet.P                                     = 1     * 'atm';
149        streamcold.Outlet.T                                     = 315   * 'K';
150    streamcold.Outlet.z                                         = [1];
151       
152OPTIONS
153
154Dynamic = false;
155GuessFile = "GuessSimplified";
156
157end
158
159FlowSheet sample_Shell_and_tubes_LMTD_cost
160
161PARAMETERS
162
163                PP              as Plugin (Brief="Physical Properties",Type="PP",LiquidModel = "PR",VapourModel = "PR",
164                                                                Components      = ["water","n-butane", "benzene" ]);
165        NComp   as Integer;
166       
167DEVICES
168
169        exchanger        as Shell_and_tubes_LMTD_cost;
170        streamhot        as source;
171        streamcold   as source;
172        OutCold      as sink;
173        OutHot       as sink;
174
175SET
176
177        NComp                           = PP.NumberOfComponents;
178
179        exchanger.LMTDcorrection = "Bowmann";
180        exchanger.ShellType             = "Eshell";
181       
182        # costs
183        exchanger.Material = "Stainless steel 316";
184       
185        exchanger.Cost(1,:) = [8.202,0.01506,0.06811];
186        exchanger.Cost(2,:) = [-0.9003,0.0906,0];
187        exchanger.Cost(3,:) = [1.35,0,0];
188        exchanger.Cost(4,:) = [-0.7844,0.0830,0];
189        exchanger.Cost(5,:) = [0.8955,0.04981,0];
190        exchanger.Cost(6,:) = [1.2002,0.07140,0];
191        exchanger.Cost(7,:) = [1.4272,0.12088,0];
192        exchanger.Cost(8,:) = [0.8608,0.23296,0];
193
194
195CONNECTIONS
196
197        streamhot.Outlet        to exchanger.InletHot;
198        streamcold.Outlet       to exchanger.InletCold;
199
200        exchanger.OutletHot  to OutHot.Inlet;
201        exchanger.OutletCold to OutCold.Inlet;
202
203SPECIFY
204
205exchanger.Details.Ud                                            = 210*'W/(m^2*K)';
206exchanger.Details.Uc                                            = 210*'W/(m^2*K)';
207
208exchanger.HotSide.PressureDrop.Pdrop    = 0.1*'kPa';
209exchanger.ColdSide.PressureDrop.Pdrop           = 0.1*'kPa';
210
211streamhot.Outlet.F              = 36    * 'kmol/h';
212streamhot.Outlet.T              = 363.15   * 'K';
213streamhot.Outlet.P              = 120   * 'kPa';
214streamhot.Outlet.z              = [1,0,0];
215
216streamcold.Outlet.F             = 10 * 'kmol/h';
217streamcold.Outlet.P             = 120 * 'kPa';
218streamcold.Outlet.T             = 303.15 * 'K';
219streamcold.Outlet.z             = [0.1,0.7, 0.2];
220
221exchanger.Details.A     = 15*'m^2';
222       
223OPTIONS
224
225 Dynamic        = false;
226 GuessFile   = "GuessShellandTubes";
227
228end
229
230
231FlowSheet sample_Shell_and_tubes_NTU_cost
232
233PARAMETERS
234
235        PP              as Plugin (Brief="Physical Properties",Type="PP",LiquidModel = "PR",VapourModel = "PR",
236                                                                Components      = ["water","n-butane", "benzene" ]);
237        NComp   as Integer;
238       
239DEVICES
240
241        exchanger        as Shell_and_tubes_NTU_cost;
242        streamhot        as source;
243        streamcold  as source;
244        OutCold      as sink;
245        OutHot       as sink;
246
247SET
248
249        NComp                           = PP.NumberOfComponents;
250
251        exchanger.ShellType             = "Eshell";
252
253        # costs
254        exchanger.Material = "Stainless steel 316";
255       
256        exchanger.Cost(1,:) = [8.202,0.01506,0.06811];
257        exchanger.Cost(2,:) = [-0.9003,0.0906,0];
258        exchanger.Cost(3,:) = [1.35,0,0];
259        exchanger.Cost(4,:) = [-0.7844,0.0830,0];
260        exchanger.Cost(5,:) = [0.8955,0.04981,0];
261        exchanger.Cost(6,:) = [1.2002,0.07140,0];
262        exchanger.Cost(7,:) = [1.4272,0.12088,0];
263        exchanger.Cost(8,:) = [0.8608,0.23296,0];
264
265CONNECTIONS
266
267        streamhot.Outlet        to exchanger.InletHot;
268        streamcold.Outlet       to exchanger.InletCold;
269
270        exchanger.OutletHot  to OutHot.Inlet;
271        exchanger.OutletCold to OutCold.Inlet;
272
273SPECIFY
274
275exchanger.Details.Ud                                            = 210*'W/(m^2*K)';
276exchanger.Details.Uc                                            = 210*'W/(m^2*K)';
277
278exchanger.HotSide.PressureDrop.Pdrop    = 0.1*'kPa';
279exchanger.ColdSide.PressureDrop.Pdrop           = 0.1*'kPa';
280
281streamhot.Outlet.F              = 36    * 'kmol/h';
282streamhot.Outlet.T              = 363.15   * 'K';
283streamhot.Outlet.P              = 120   * 'kPa';
284streamhot.Outlet.z              = [1,0,0];
285
286streamcold.Outlet.F             = 10 * 'kmol/h';
287streamcold.Outlet.P             = 120 * 'kPa';
288streamcold.Outlet.T             = 303.15 * 'K';
289streamcold.Outlet.z             = [0.1,0.7, 0.2];
290
291exchanger.Details.A     = 15*'m^2';
292       
293OPTIONS
294
295 Dynamic        = false;
296
297end
Note: See TracBrowser for help on using the repository browser.