source: branches/newlanguage/eml/heat_exchangers/HEX_Engine.mso @ 156

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

updated double pipe model

  • Property svn:eol-style set to native
  • Property svn:keywords set to Id
File size: 13.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*                       Heat Exchangers Abstract Models
17*--------------------------------------------------------------------
18*       - Inlet_Main_Stream     : Inlet Streams
19*--------------------------------------------------------------------
20*                                       - Hot   : Inlet Hot Stream
21*                                       - Cold  : Inlet Cold Stream
22*--------------------------------------------------------------------
23*       - Outlet_Main_Stream    : Outlet Streams
24*--------------------------------------------------------------------
25*                                       - Hot   : Outlet Hot Stream
26*                                       - Cold  : Outlet Cold Stream
27*--------------------------------------------------------------------
28*       - Main_Properties               : Physical Properties for Hot and Cold Side
29*--------------------------------------------------------------------
30*               Physical_Properties
31*                       Properties_In_Out       :       Inlet/Outlet Physical Properties
32*                       Properties_Average      :       Average Physical Properties
33*                       Properties_Wall         :       Physical Properties at Wall Temperature
34*--------------------------------------------------------------------
35*       - Tube_Side_Main        : Tube Side Main Variables
36*----------------------------------------------------------------------
37*                       Tube_Pdrop                      : Tube Side Pressure Drop
38*                       Tube_Heat_Transfer      : Tube Side Heat Transfer
39*----------------------------------------------------------------------
40*       - Shell_Side_Main       : Shell Side Main Variables
41*----------------------------------------------------------------------
42*                       Shell_Pdrop                     : Shell Side Pressure Drop
43*                       Shell_Heat_Transfer     : Shell Side Heat Transfer
44*----------------------------------------------------------------------
45*       - Baffles_Main  : Baffles Spacing
46*----------------------------------------------------------------------
47*       - Main_Resistances      : Thermal Resistances
48*----------------------------------------------------------------------
49*       - Details_Main  : Heat Exchanger Thermal Details
50*----------------------------------------------------------------------
51*       - Main_Pdrop : Heat Exchanger Pressure Drop (Simplified)
52*               Basic_Pdrop     : Hot and Cold Side Pressure Drop
53*----------------------------------------------------------------------
54*       - Main_DoublePipe : Double Pipe Heat Exchanger Block
55*               DoublePipe_HeatTransfer
56*               DoublePipe_PressureDrop
57*----------------------------------------------------------------------
58* Author: Gerson Balbueno Bicca
59* $Id: HEX_Engine.mso 156 2007-02-07 13:40:14Z bicca $
60*--------------------------------------------------------------------*#
61
62using "streams";
63
64Model Inlet_Main_Stream
65       
66ATTRIBUTES
67        Pallete = false;
68        Brief = "write some information";
69        Info =
70        "write some information";
71       
72VARIABLES
73Hot     as stream       (Brief="Inlet Hot Stream");
74Cold    as stream       (Brief="Inlet Cold Stream");
75
76end
77
78Model Outlet_Main_Stream
79       
80ATTRIBUTES
81        Pallete = false;
82        Brief = "write some information";
83        Info =
84        "write some information";
85       
86VARIABLES
87# Must be streamPH
88Hot     as streamPH (Brief="Outlet Hot Stream");
89Cold    as streamPH (Brief="Outlet Cold Stream");
90
91end
92
93#=====================================================================
94# Heat Exchangers Physical Properties
95#=====================================================================
96 
97Model Properties_Average
98       
99ATTRIBUTES
100        Pallete = false;
101        Brief = "write some information";
102        Info =
103        "write some information";
104       
105VARIABLES
106Mw              as molweight            (Brief="Average Mol Weight",Default=75, Lower=1, Upper=1e8);
107T                       as temperature          (Brief="Average  Temperature",Lower=50);
108P                       as pressure                     (Brief="Average  Pressure",Default=1, Lower=1e-10, Upper=30);
109rho             as dens_mass            (Brief="Stream Density" ,Default=1000, Lower=1e-3, Upper=5e5);
110Mu              as viscosity            (Brief="Stream Viscosity",Lower=0.0001);
111Cp                      as cp_mol                       (Brief="Stream Molar Heat Capacity", Upper=1e10);
112K                       as conductivity         (Brief="Stream Thermal Conductivity", Default=1.0, Lower=1e-5, Upper=500);
113
114end
115
116Model Properties_In_Out
117       
118ATTRIBUTES
119        Pallete = false;
120        Brief = "write some information";
121        Info =
122        "write some information";
123       
124VARIABLES
125Fw              as flow_mass            (Brief="Stream Mass Flow");
126rho             as dens_mass            (Brief="Stream Density" ,Default=1000, Lower=1e-3, Upper=5e5);
127Mu              as viscosity            (Brief="Stream Viscosity",Default=1, Lower=1e-5, Upper=1e5);
128Cp                      as cp_mol                       (Brief="Stream Molar Heat Capacity", Upper=1e10);
129K                       as conductivity         (Brief="Stream Thermal Conductivity", Default=1.0, Lower=1e-5, Upper=500);
130
131end
132
133Model Properties_Wall
134       
135ATTRIBUTES
136        Pallete = false;
137        Brief = "write some information";
138        Info =
139        "write some information";
140       
141VARIABLES
142Mu              as viscosity            (Brief="Stream Viscosity",Default=1, Lower=1e-5, Upper=1e5);
143Cp                      as cp_mol                       (Brief="Stream Molar Heat Capacity", Upper=1e10);
144K                       as conductivity         (Brief="Stream Thermal Conductivity", Default=1.0, Lower=1e-5, Upper=500);
145Twall           as temperature  (Brief="Wall Temperature",Lower=50);
146
147end
148
149Model Physical_Properties
150       
151ATTRIBUTES
152        Pallete = false;
153        Brief = "write some information";
154        Info =
155        "write some information";
156       
157VARIABLES
158Inlet           as Properties_In_Out    (Brief="Properties at Inlet Stream");
159Average         as Properties_Average   (Brief="Properties at Average Temperature");
160Outlet          as Properties_In_Out    (Brief="Properties at Outlet Stream");
161Wall                    as Properties_Wall              (Brief="Properties at Wall Temperature");
162
163end
164
165Model Main_Properties
166ATTRIBUTES
167        Pallete = false;
168        Brief = "write some information";
169        Info =
170        "write some information";
171       
172VARIABLES
173
174Hot  as Physical_Properties  (Brief="Hot Stream");
175Cold as Physical_Properties  (Brief="Cold Stream");
176
177end
178
179Model Tube_Pdrop
180       
181ATTRIBUTES
182        Pallete = false;
183        Brief = "write some information";
184        Info =
185        "write some information";
186       
187VARIABLES
188PdTube                  as press_delta  (Brief="Tube Pressure Drop",Default=0.01, Lower=1e-10);
189Pdtotal                         as press_delta  (Brief="Total Pressure Drop",Default=0.01, Lower=1e-10);
190Pdnozzle_in     as press_delta  (Brief="Inlet Nozzle Pressure Drop",Default=0.01, Lower=0);
191Pdnozzle_out    as press_delta  (Brief="Outlet Nozzle Pressure Drop",Default=0.01, Lower=0);
192Vnozzle_in      as velocity             (Brief="Inlet Nozzle Velocity",Default=1, Upper=1e5, Lower=0);
193Vnozzle_out     as velocity             (Brief="Outlet Nozzle Velocity",Default=1, Upper=1e5, Lower=0);
194
195EQUATIONS
196"Total Pressure Drop"
197        Pdtotal = PdTube + Pdnozzle_in + Pdnozzle_out;
198
199end
200
201Model Tube_Heat_Transfer
202       
203ATTRIBUTES
204        Pallete = false;
205        Brief = "write some information";
206        Info =
207        "write some information";
208       
209VARIABLES
210Re              as positive                             (Brief="Tube Side Reynolds Number",Default=1000,Lower=1);
211htube   as heat_trans_coeff (Brief="Tube Side Film Coefficient",Default=1,Lower=1e-12, Upper=1e6);
212PR              as positive                             (Brief="Tube Side Prandtl Number",Default=0.5,Lower=1e-8);
213PRw     as positive                             (Brief="Tube Side Prandtl Number at Wall Temperature",Default=0.5,Lower=1e-8);
214Phi             as positive                             (Brief="Phi Correction",Default=1,Lower=1e-3);
215Vtube   as velocity                             (Brief="Tube Side Velocity",Lower=1e-8);
216
217end
218
219Model Shell_Pdrop
220       
221ATTRIBUTES
222        Pallete = false;
223        Brief = "write some information";
224        Info =
225        "write some information";
226       
227VARIABLES
228PdCross                         as press_delta          (Brief="Cross Flow Pressure Drop",Default=0.01, Lower=0);
229PdEndZones              as press_delta          (Brief="End Zones Pressure Drop",Default=0.01, Lower=0);
230Pdwindow                        as press_delta          (Brief="Window Pressure Drop",Default=0.01, Lower=1e-10);
231Pdtotal                                 as press_delta          (Brief="Total Pressure Drop",Default=0.01, Lower=0);
232Pdnozzle_in             as press_delta          (Brief="Inlet Nozzle Pressure Drop",Default=0.01, Lower=0);
233Pdnozzle_out            as press_delta          (Brief="Outlet Nozzle Pressure Drop",Default=0.01, Lower=0);
234Vnozzle_in              as velocity                     (Brief="Inlet Nozzle Velocity",Default=1, Upper=1e5, Lower=0);
235Vnozzle_out             as velocity                     (Brief="Outlet Nozzle Velocity",Default=1, Upper=1e5, Lower=0);
236RVsquare_out            as positive                     (Brief = "Outlet Nozzle rho-V^2", Default=1, Upper=1e6, Unit = "kg/s^2/m");
237RVsquare_in             as positive                     (Brief = "Inlet Nozzle rho-V^2", Default=1, Upper=1e6, Unit = "kg/s^2/m");
238
239EQUATIONS
240"Shell Side Total Pressure Drop"
241        Pdtotal = PdCross + PdEndZones + Pdnozzle_in + Pdnozzle_out + Pdwindow;
242
243end
244
245Model Shell_Heat_Transfer
246       
247ATTRIBUTES
248        Pallete = false;
249        Brief = "write some information";
250        Info =
251        "write some information";
252       
253VARIABLES
254Re              as positive                             (Brief="Shell Side Reynolds Number",Default=100,Lower=1);
255PR              as positive                             (Brief="Shell Side Prandtl Number",Default=0.7,Lower=1e-6);
256PRw             as positive                             (Brief="Shell Side Prandtl Number",Default=0.5,Lower=1e-8);
257hshell          as heat_trans_coeff     (Brief="Shell Side Film Coefficient",Default=1,Lower=1e-12, Upper=1e6);
258Phi     as positive                             (Brief="Phi Correction",Default=1,Lower=1e-3);
259Ji                      as constant                     (Brief="Shell Side Ji Factor",Default=0.05);
260Jr                      as positive                             (Brief="Shell Side Jr Factor",Lower=10e-6);
261Jl                      as positive                             (Brief="Shell Side Jl Factor",Lower=10e-6);
262Jb                      as positive                             (Brief="Shell Side Jb Factor",Lower=10e-6);
263Jc                      as positive                             (Brief="Shell Side Jc Factor",Lower=10e-6);
264Js                      as positive                             (Brief="Shell Side Js Factor",Lower=10e-6);
265Jtotal          as positive                             (Brief="Shell Side Jtotal Factor",Lower=10e-6);
266Sm              as area                                 (Brief="Shell Side Cross Flow Area",Default=0.05,Lower=10e-6);
267
268end
269
270Model Baffles_Main
271       
272ATTRIBUTES
273        Pallete = false;
274        Brief = "write some information";
275        Info =
276        "write some information";
277       
278VARIABLES
279Ls              as length               (Brief="Central Baffle Spacing",Lower=1e-8);
280Lsi             as length               (Brief="Inlet Baffle Spacing",Lower=1e-8);
281Lso             as length               (Brief="Outlet Baffle Spacing",Lower=1e-8);
282
283end
284
285Model Main_Resistances
286       
287ATTRIBUTES
288        Pallete = false;
289        Brief = "write some information";
290        Info =
291        "write some information";
292       
293VARIABLES
294
295Rtube   as positive     (Brief="Tube Resistance",Unit='m^2*K/kW',Lower=1e-6);
296Rwall           as positive     (Brief="Wall Resistance",Unit='m^2*K/kW',Lower=1e-6);
297Rshell          as positive     (Brief="Shell Resistance",Unit='m^2*K/kW',Lower=1e-6);
298Rfi                     as positive     (Brief="Inside Fouling Resistance",Unit='m^2*K/kW',Default=1e-6,Lower=0);
299Rfo             as positive     (Brief="Outside Fouling Resistance",Unit='m^2*K/kW',Default=1e-6,Lower=0);
300
301end
302
303Model Details_Main
304       
305ATTRIBUTES
306        Pallete = false;
307        Brief = "write some information";
308        Info =
309        "write some information";
310       
311VARIABLES
312A               as area                                 (Brief="Exchange Surface Area");
313Q               as power                                (Brief="Heat Transfer", Default=7000, Lower=1e-6, Upper=1e10);
314Uc      as heat_trans_coeff (Brief="Overall Heat Transfer Coefficient Clean",Default=1,Lower=1e-6,Upper=1e10);
315Ud      as heat_trans_coeff (Brief="Overall Heat Transfer Coefficient Dirty",Default=1,Lower=1e-6,Upper=1e10);
316Ch      as positive                             (Brief="Hot Stream Heat Capacity",Lower=1e-3,Default=1e3,Unit='W/K');
317Cc      as positive                             (Brief="Cold Stream Heat Capacity",Lower=1e-3,Default=1e3,Unit='W/K');
318Cr      as positive                             (Brief="Heat Capacity Ratio",Default=0.5,Lower=1e-6);
319Cmin  as positive                               (Brief="Minimum Heat Capacity",Lower=1e-10,Default=1e3,Unit='W/K');
320Cmax as positive                                (Brief="Maximum Heat Capacity",Lower=1e-10,Default=1e3,Unit='W/K');
321NTU     as positive                             (Brief="Number of Units Transference",Default=0.05,Lower=1e-10);
322
323EQUATIONS
324"Number of Units Transference"
325        NTU*Cmin = Ud*A;
326       
327end
328
329Model Tube_Side_Main
330
331ATTRIBUTES
332        Pallete = false;
333        Brief = "write some information";
334        Info =
335        "write some information";
336       
337VARIABLES
338PressureDrop    as Tube_Pdrop                   (Brief="Tube Side Pressure Drop");
339HeatTransfer    as Tube_Heat_Transfer (Brief="Tube Side Heat Transfer");
340
341end
342
343Model Shell_Side_Main
344       
345ATTRIBUTES
346        Pallete = false;
347        Brief = "write some information";
348        Info =
349        "write some information";
350       
351VARIABLES
352PressureDrop    as Shell_Pdrop                           (Brief="Shell Side Pressure Drop");
353HeatTransfer    as Shell_Heat_Transfer  (Brief= "Shell Side Heat Transfer");
354
355end
356
357Model Basic_Pdrop
358       
359ATTRIBUTES
360        Pallete = false;
361        Brief = "write some information";
362        Info =
363        "write some information";
364       
365VARIABLES
366Pdrop   as press_delta  (Brief="Pressure Drop",Default=0.01, Lower=0);
367FPdrop  as Real                 (Brief="Pressure Drop : Fraction of Inlet",Lower=0,Upper=0.8);
368
369end
370
371Model Main_Pdrop
372       
373ATTRIBUTES
374        Pallete = false;
375        Brief = "write some information";
376        Info =
377        "write some information";
378       
379VARIABLES
380Hot     as Basic_Pdrop (Brief="Hot Stream");
381Cold    as Basic_Pdrop (Brief="Cold Stream");
382
383end     
384
385Model DoublePipe_HeatTransfer
386       
387ATTRIBUTES
388        Pallete = false;
389        Brief = "write some information";
390        Info =
391        "write some information";
392       
393VARIABLES
394As                      as area                                         (Brief="Cross Sectional Area for Flow",Default=0.05,Lower=1e-8);
395Dh              as length                                       (Brief="Hydraulic Diameter of Pipe for Heat Transfer",Lower=1e-8);
396Re              as positive                             (Brief="Reynolds Number",Default=100,Lower=1);
397hcoeff  as heat_trans_coeff (Brief="Film Coefficient",Default=1,Lower=1e-12, Upper=1e6);
398fi              as fricfactor                   (Brief="Friction Factor", Default=0.05, Lower=1e-10, Upper=2000);
399Nu              as positive                             (Brief="Nusselt Number",Default=0.5,Lower=1e-8);
400PR              as positive                             (Brief="Prandtl Number",Default=0.5,Lower=1e-8);
401Phi             as positive                             (Brief="Phi Correction",Default=1,Lower=1e-3);
402Vmean   as velocity                             (Brief="Tube Velocity",Lower=1e-8);
403
404end
405
406Model DoublePipe_PressureDrop
407       
408ATTRIBUTES
409        Pallete = false;
410        Brief = "write some information";
411        Info =
412        "write some information";
413       
414VARIABLES
415
416Dh      as length               (Brief="Hydraulic Diameter of Pipe for Pressure Drop",Lower=1e-6);
417Pdrop   as press_delta  (Brief="Pressure Drop",Default=0.01, Lower=1e-10);
418fi      as fricfactor   (Brief="Friction Factor", Default=0.05, Lower=1e-10, Upper=2000);
419Re      as positive             (Brief="Reynolds Number",Default=100,Lower=1);
420
421end     
422
423Model Main_DoublePipe
424       
425ATTRIBUTES
426        Pallete = false;
427        Brief = "write some information";
428        Info =
429        "write some information";
430       
431VARIABLES
432HeatTransfer as DoublePipe_HeatTransfer (Brief="Double Pipe Heat Transfer");
433PressureDrop as DoublePipe_PressureDrop (Brief="Double Pipe Pressure Drop");
434
435end
Note: See TracBrowser for help on using the repository browser.