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

Last change on this file since 548 was 526, checked in by gerson bicca, 15 years ago

updated double pipe sample
File size: 9.2 KB
RevLine 
[526]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.Geometry.DoInner              =       42.16*'mm';
93exchanger.Geometry.DiInner              =       35.05*'mm';
94exchanger.Geometry.DiOuter              =       52.50*'mm';
95exchanger.Geometry.Kwall                        = 0.057 *'kW/m/K';
96exchanger.Geometry.Lpipe                        =       9*'m';
97
98#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
99#       Fouling
100#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
101exchanger.Geometry.Rfi = 0.00018*'m^2*K/W';
102exchanger.Geometry.Rfo = 0*'m^2*K/W';
103
104SPECIFY
105#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
106#               Hot Stream
107#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
108        #InletHot.Fw                    = 4396           * 'kg/h';
109        InletHot.Outlet.F                       = 47.2038        * 'kmol/h';
110        InletHot.Outlet.T       = 338.71   * 'K';
111        InletHot.Outlet.P       = 413.685 * 'kPa';
112        InletHot.Outlet.z       = [1,0];
113       
114#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
115#       Cold Stream
116#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
117        #InletCold.Fw                   = 4536          * 'kg/h';
118        InletCold.Outlet.F      = 58.0694       * 'kmol/h';
119        InletCold.Outlet.T      = 288.71        * 'K';
120        InletCold.Outlet.P      = 413.685 * 'kPa';
121    InletCold.Outlet.z  = [0,1];
122       
123OPTIONS
124
125Dynamic  = false;
126GuessFile = "TESTE_1";
127
128end
129
130FlowSheet TESTE_2
131
132PARAMETERS
133        PP                      as Plugin (Type ="PP", LiquidModel = "PR", VapourModel = "PR",
134                                                                        Components = ["n-hexane","benzene"]);
135        NComp           as Integer;
136       
137DEVICES
138
139        exchanger       as DoublePipe_LMTD;
140        InletHot                as simple_source;
141        InletCold       as simple_source;
142        OutletHot               as simple_sink;
143        OutletCold      as simple_sink;
144       
145       
146CONNECTIONS
147
148        InletCold.Outlet                                to exchanger.InletInner;
149        InletHot.Outlet                         to exchanger.InletOuter;
150        exchanger.OutletOuter   to OutletHot.Inlet;
151        exchanger.OutletInner   to OutletCold.Inlet;
152       
153SET
154
155NComp                   = PP.NumberOfComponents;
156
157#exchanger.CalculationApproach = "Simplified";
158exchanger.Qestimated = 40*'kW';
159
160exchanger.CalculationApproach = "Full";
161
162#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
163#   Options
164#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
165
166exchanger.FlowDirection                                                 = "counter";
167exchanger.InnerTurbulentCorrelation     = "SiederTate";
168exchanger.InnerLaminarCorrelation       = "Schlunder";
169
170exchanger.OuterTurbulentCorrelation     = "SiederTate";
171exchanger.OuterLaminarCorrelation       = "Schlunder";
172       
173#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
174# Double Pipe Geometrical Parameters and Alocation
175#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
176
177exchanger.Geometry.DoInner              =       42.16*'mm';
178exchanger.Geometry.DiInner              =       35.05*'mm';
179exchanger.Geometry.DiOuter              =       52.50*'mm';
180exchanger.Geometry.Kwall                        = 0.057 *'kW/m/K';
181exchanger.Geometry.Lpipe                        =       9*'m';
182
183#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
184#       Fouling
185#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
186exchanger.Geometry.Rfi = 0.00018*'m^2*K/W';
187exchanger.Geometry.Rfo = 0*'m^2*K/W';
188
189SPECIFY
190#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
191#               Hot Stream
192#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
193        #InletHot.Fw                    = 4396           * 'kg/h';
194        InletHot.Outlet.F                       = 47.2038        * 'kmol/h';
195        InletHot.Outlet.T       = 338.71   * 'K';
196        InletHot.Outlet.P       = 413.685 * 'kPa';
197        InletHot.Outlet.z       = [1,0];
198       
199#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
200#       Cold Stream
201#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
202        #InletCold.Fw                   = 4536          * 'kg/h';
203        InletCold.Outlet.F      = 58.0694       * 'kmol/h';
204        InletCold.Outlet.T      = 288.71        * 'K';
205        InletCold.Outlet.P      = 413.685 * 'kPa';
206    InletCold.Outlet.z  = [0,1];
207       
208OPTIONS
209
210Dynamic  = false;
211GuessFile = "TESTE_2";
212
213end
214
215FlowSheet TESTE_3
216
217PARAMETERS
218        PP                      as Plugin (Type ="PP", LiquidModel = "PR", VapourModel = "PR",
219                                                                        Components = [ "aniline", "toluene"]);
220        NComp           as Integer;
221       
222DEVICES
223
224        exchanger       as DoublePipe_LMTD;
225        InletHot                as simple_source;
226        InletCold       as simple_source;
227        OutletHot               as simple_sink;
228        OutletCold      as simple_sink;
229       
230       
231CONNECTIONS
232
233        InletCold.Outlet                                to exchanger.InletInner;
234        InletHot.Outlet                         to exchanger.InletOuter;
235        exchanger.OutletOuter   to OutletHot.Inlet;
236        exchanger.OutletInner   to OutletCold.Inlet;
237       
238SET
239
240NComp                   = PP.NumberOfComponents;
241
242#exchanger.CalculationApproach = "Simplified";
243exchanger.Qestimated = 30*'kW';
244
245exchanger.CalculationApproach = "Full";
246
247#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
248#   Options
249#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
250
251exchanger.FlowDirection                                                 = "counter";
252exchanger.InnerTurbulentCorrelation     = "SiederTate";
253exchanger.InnerLaminarCorrelation       = "Schlunder";
254
255exchanger.OuterTurbulentCorrelation     = "SiederTate";
256exchanger.OuterLaminarCorrelation       = "Schlunder";
257       
258#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
259# Double Pipe Geometrical Parameters and Alocation
260#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
261
262exchanger.Geometry.DoInner              =       42.16*'mm';
263exchanger.Geometry.DiInner              =       35.05*'mm';
264exchanger.Geometry.DiOuter              =       52.50*'mm';
265exchanger.Geometry.Kwall                        = 0.057 *'kW/m/K';
266exchanger.Geometry.Lpipe                        =       9*'m';
267
268#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
269#       Fouling
270#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
271exchanger.Geometry.Rfi = 0.00018*'m^2*K/W';
272exchanger.Geometry.Rfo = 0*'m^2*K/W';
273
274SPECIFY
275#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
276#               Hot Stream
277#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
278        #InletHot.Fw                    = 4396           * 'kg/h';
279        InletHot.Outlet.F                       = 47.2038        * 'kmol/h';
280        InletHot.Outlet.T       = 413   * 'K';
281        InletHot.Outlet.P       = 413.685 * 'kPa';
282        InletHot.Outlet.z       = [1,0];
283       
284#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
285#       Cold Stream
286#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
287        #InletCold.Fw                   = 4536          * 'kg/h';
288        InletCold.Outlet.F      = 58.0694       * 'kmol/h';
289        InletCold.Outlet.T      = 330           * 'K';
290        InletCold.Outlet.P      = 413.685 * 'kPa';
291    InletCold.Outlet.z  = [0,1];
292       
293OPTIONS
294
295Dynamic  = false;
296GuessFile = "TESTE_3";
297end
Note: See TracBrowser for help on using the repository browser.