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

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

Adding symbol for latex.

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