source: trunk/eml/heat_exchangers/HEX_Engine.mso @ 261

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

documenting PHE heat exchanger

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