source: trunk/sample/heat_exchangers/Sample_DoublePipe_convergence.mso @ 494

Last change on this file since 494 was 473, checked in by gerson bicca, 16 years ago

Double Pipe: added sample to explain how to solve convergence problems
File size: 9.2 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 a Double Pipe Heat Exchanger - LMTD Method
17*--------------------------------------------------------------------
18*
19* This sample file needs VRTherm (www.vrtech.com.br) to run.
20*
21*----------------------------------------------------------------------
22* Author: Gerson B. Bicca
23* $Id: Sample_DoublePipe_convergence.mso $
24*--------------------------------------------------------------------*#
25using "heat_exchangers/DoublePipe";
26
27
28#* Note: This sample file explain how to obtain a solution using a double pipe simplified model at conditions with better convergence characteristics
29
30step 1:
31        - use the option below to obtain a good approximation of the final solution:
32
33        CalculationApproach = "Simplified";             setting the simplified model
34        Qestimated = 70*'kW';                                                   estimate heat duty for the simplified model (a good value to heat duty)
35
36step 2:
37        - Once you got convergence, use the EMSO graphical interface and save the result to a file (eg. result.rlt)
38
39step 3:
40        - CalculationApproach = "Full";                                 setting the full model
41        - guessFile="result.rlt"                                                                use the option guessFile
42
43*#
44
45FlowSheet TESTE_1
46
47PARAMETERS
48        PP                      as Plugin (Type ="PP", LiquidModel = "PR", VapourModel = "PR",
49                                                                        Components = ["aniline","benzene"]);
50        NComp           as Integer;
51       
52DEVICES
53
54        exchanger       as DoublePipe_LMTD;
55        InletHot                as simple_source;
56        InletCold       as simple_source;
57        OutletHot               as simple_sink;
58        OutletCold      as simple_sink;
59       
60       
61CONNECTIONS
62
63        InletCold.Outlet                                to exchanger.InletInner;
64        InletHot.Outlet                         to exchanger.InletOuter;
65        exchanger.OutletOuter   to OutletHot.Inlet;
66        exchanger.OutletInner   to OutletCold.Inlet;
67       
68SET
69
70NComp                   = PP.NumberOfComponents;
71
72#exchanger.CalculationApproach = "Simplified"; 
73exchanger.Qestimated = 30*'kW';
74
75exchanger.CalculationApproach = "Full";
76
77#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
78#   Options
79#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
80
81exchanger.FlowDirection                                                 = "counter";
82exchanger.InnerTurbulentCorrelation     = "SiederTate";
83exchanger.InnerLaminarCorrelation       = "Schlunder";
84
85exchanger.OuterTurbulentCorrelation     = "SiederTate";
86exchanger.OuterLaminarCorrelation       = "Schlunder";
87       
88#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
89# Double Pipe Geometrical Parameters and Alocation
90#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
91
92exchanger.DoInner               =       42.16*'mm';
93exchanger.DiInner               =       35.05*'mm';
94exchanger.DiOuter               =       52.50*'mm';
95exchanger.Kwall                         = 0.057 *'kW/m/K';
96exchanger.Lpipe                 =       9*'m';
97
98#exchanger.Donozzle_Inner       =       36.63*'mm';
99#exchanger.Dinozzle_Inner       =       36.63*'mm';
100
101#exchanger.Donozzle_Outer       =       42.72*'mm';
102#exchanger.Dinozzle_Outer       =       42.72*'mm';
103#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
104#       Fouling
105#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
106exchanger.Rfi = 0.00018*'m^2*K/W';
107exchanger.Rfo = 0*'m^2*K/W';
108
109SPECIFY
110#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
111#               Hot Stream
112#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
113        #InletHot.Fw                    = 4396           * 'kg/h';
114        InletHot.Outlet.F                       = 47.2038        * 'kmol/h';
115        InletHot.Outlet.T       = 338.71   * 'K';
116        InletHot.Outlet.P       = 413.685 * 'kPa';
117        InletHot.Outlet.z       = [1,0];
118       
119#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
120#       Cold Stream
121#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
122        #InletCold.Fw                   = 4536          * 'kg/h';
123        InletCold.Outlet.F      = 58.0694       * 'kmol/h';
124        InletCold.Outlet.T      = 288.71        * 'K';
125        InletCold.Outlet.P      = 413.685 * 'kPa';
126    InletCold.Outlet.z  = [0,1];
127       
128OPTIONS
129
130Dynamic  = false;
131GuessFile = "TESTE_1";
132
133end
134
135FlowSheet TESTE_2
136
137PARAMETERS
138        PP                      as Plugin (Type ="PP", LiquidModel = "PR", VapourModel = "PR",
139                                                                        Components = ["n-hexane","benzene"]);
140        NComp           as Integer;
141       
142DEVICES
143
144        exchanger       as DoublePipe_LMTD;
145        InletHot                as simple_source;
146        InletCold       as simple_source;
147        OutletHot               as simple_sink;
148        OutletCold      as simple_sink;
149       
150       
151CONNECTIONS
152
153        InletCold.Outlet                                to exchanger.InletInner;
154        InletHot.Outlet                         to exchanger.InletOuter;
155        exchanger.OutletOuter   to OutletHot.Inlet;
156        exchanger.OutletInner   to OutletCold.Inlet;
157       
158SET
159
160NComp                   = PP.NumberOfComponents;
161
162#exchanger.CalculationApproach = "Simplified";
163exchanger.Qestimated = 40*'kW';
164
165exchanger.CalculationApproach = "Full";
166
167#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
168#   Options
169#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
170
171exchanger.FlowDirection                                                 = "counter";
172exchanger.InnerTurbulentCorrelation     = "SiederTate";
173exchanger.InnerLaminarCorrelation       = "Schlunder";
174
175exchanger.OuterTurbulentCorrelation     = "SiederTate";
176exchanger.OuterLaminarCorrelation       = "Schlunder";
177       
178#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
179# Double Pipe Geometrical Parameters and Alocation
180#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
181
182exchanger.DoInner               =       42.16*'mm';
183exchanger.DiInner               =       35.05*'mm';
184exchanger.DiOuter               =       52.50*'mm';
185exchanger.Kwall                         = 0.057 *'kW/m/K';
186exchanger.Lpipe                 =       9*'m';
187
188#exchanger.Donozzle_Inner       =       36.63*'mm';
189#exchanger.Dinozzle_Inner       =       36.63*'mm';
190
191#exchanger.Donozzle_Outer       =       42.72*'mm';
192#exchanger.Dinozzle_Outer       =       42.72*'mm';
193#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
194#       Fouling
195#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
196exchanger.Rfi = 0.00018*'m^2*K/W';
197exchanger.Rfo = 0*'m^2*K/W';
198
199SPECIFY
200#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
201#               Hot Stream
202#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
203        #InletHot.Fw                    = 4396           * 'kg/h';
204        InletHot.Outlet.F                       = 47.2038        * 'kmol/h';
205        InletHot.Outlet.T       = 338.71   * 'K';
206        InletHot.Outlet.P       = 413.685 * 'kPa';
207        InletHot.Outlet.z       = [1,0];
208       
209#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
210#       Cold Stream
211#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
212        #InletCold.Fw                   = 4536          * 'kg/h';
213        InletCold.Outlet.F      = 58.0694       * 'kmol/h';
214        InletCold.Outlet.T      = 288.71        * 'K';
215        InletCold.Outlet.P      = 413.685 * 'kPa';
216    InletCold.Outlet.z  = [0,1];
217       
218OPTIONS
219
220Dynamic  = false;
221GuessFile = "TESTE_2";
222
223end
224
225FlowSheet TESTE_3
226
227PARAMETERS
228        PP                      as Plugin (Type ="PP", LiquidModel = "PR", VapourModel = "PR",
229                                                                        Components = [ "aniline", "toluene"]);
230        NComp           as Integer;
231       
232DEVICES
233
234        exchanger       as DoublePipe_LMTD;
235        InletHot                as simple_source;
236        InletCold       as simple_source;
237        OutletHot               as simple_sink;
238        OutletCold      as simple_sink;
239       
240       
241CONNECTIONS
242
243        InletCold.Outlet                                to exchanger.InletInner;
244        InletHot.Outlet                         to exchanger.InletOuter;
245        exchanger.OutletOuter   to OutletHot.Inlet;
246        exchanger.OutletInner   to OutletCold.Inlet;
247       
248SET
249
250NComp                   = PP.NumberOfComponents;
251
252#exchanger.CalculationApproach = "Simplified";
253exchanger.Qestimated = 30*'kW';
254
255exchanger.CalculationApproach = "Full";
256
257#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
258#   Options
259#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
260
261exchanger.FlowDirection                                                 = "counter";
262exchanger.InnerTurbulentCorrelation     = "SiederTate";
263exchanger.InnerLaminarCorrelation       = "Schlunder";
264
265exchanger.OuterTurbulentCorrelation     = "SiederTate";
266exchanger.OuterLaminarCorrelation       = "Schlunder";
267       
268#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
269# Double Pipe Geometrical Parameters and Alocation
270#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
271
272exchanger.DoInner               =       42.16*'mm';
273exchanger.DiInner               =       35.05*'mm';
274exchanger.DiOuter               =       52.50*'mm';
275exchanger.Kwall                         = 0.057 *'kW/m/K';
276exchanger.Lpipe                 =       9*'m';
277
278#exchanger.Donozzle_Inner       =       36.63*'mm';
279#exchanger.Dinozzle_Inner       =       36.63*'mm';
280
281#exchanger.Donozzle_Outer       =       42.72*'mm';
282#exchanger.Dinozzle_Outer       =       42.72*'mm';
283#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
284#       Fouling
285#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
286exchanger.Rfi = 0.00018*'m^2*K/W';
287exchanger.Rfo = 0*'m^2*K/W';
288
289SPECIFY
290#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
291#               Hot Stream
292#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
293        #InletHot.Fw                    = 4396           * 'kg/h';
294        InletHot.Outlet.F                       = 47.2038        * 'kmol/h';
295        InletHot.Outlet.T       = 413   * 'K';
296        InletHot.Outlet.P       = 413.685 * 'kPa';
297        InletHot.Outlet.z       = [1,0];
298       
299#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
300#       Cold Stream
301#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
302        #InletCold.Fw                   = 4536          * 'kg/h';
303        InletCold.Outlet.F      = 58.0694       * 'kmol/h';
304        InletCold.Outlet.T      = 330           * 'K';
305        InletCold.Outlet.P      = 413.685 * 'kPa';
306    InletCold.Outlet.z  = [0,1];
307       
308OPTIONS
309
310Dynamic  = false;
311GuessFile = "TESTE_3";
312end
Note: See TracBrowser for help on using the repository browser.