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

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

REMOVED SOME VARIABLES

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