source: mso/eml/heat_exchangers/HeatExchangerDetailed.mso @ 100

Last change on this file since 100 was 100, checked in by gerson bicca, 16 years ago

updated

  • Property svn:eol-style set to native
  • Property svn:keywords set to Id
File size: 70.7 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* Author: Gerson Balbueno Bicca
17* $Id: HeatExchangerDetailed.mso 100 2007-01-09 14:15:56Z bicca $
18*------------------------------------------------------------------*#
19
20using "HEX_Engine";
21
22#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
23#       Basic Model for Detailed Shell and Tubes Heat Exchangers
24#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
25
26Model HeatExchangerDetailed_Basic
27       
28PARAMETERS
29
30ext PP          as CalcObject   (Brief="External Physical Properties");
31ext NComp       as Integer      (Brief="Number of Components");
32        HE              as CalcObject   (Brief="STHE Calculations",File="heatex");
33        M(NComp)        as molweight    (Brief="Component Mol Weight");
34       
35VARIABLES
36
37in  Inlet               as Inlet_Main_Stream;   
38out Outlet              as Outlet_Main_Stream; 
39        Properties      as Main_Properties;             
40        Details                 as Details_Main;
41        Tubes                   as Tube_Side_Main;     
42        Shell                   as Shell_Side_Main;
43        Resistances     as Main_Resistances;
44        Baffles         as Baffles_Main;
45
46SET
47        M   = PP.MolecularWeight();
48
49EQUATIONS
50
51#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
52#                                               Properties                                                                      #
53#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
54
55"Hot Stream Average Temperature"
56        Properties.Hot.Average.T = 0.5*Inlet.Hot.T + 0.5*Outlet.Hot.T;
57       
58"Cold Stream Average Temperature"
59        Properties.Cold.Average.T = 0.5*Inlet.Cold.T + 0.5*Outlet.Cold.T;
60       
61"Hot Stream Average Pressure"
62        Properties.Hot.Average.P = 0.5*Inlet.Hot.P+0.5*Outlet.Hot.P;
63       
64"Cold Stream Average Pressure"
65        Properties.Cold.Average.P = 0.5*Inlet.Cold.P+0.5*Outlet.Cold.P;
66       
67"Hot Stream Average Molecular Weight"
68        Properties.Hot.Average.Mw = sum(M*Inlet.Hot.z);
69
70"Cold Stream Average Molecular Weight"
71        Properties.Cold.Average.Mw = sum(M*Inlet.Cold.z);
72       
73
74if Inlet.Cold.v equal 0
75       
76        then   
77       
78"Cold Stream Average Heat Capacity"
79        Properties.Cold.Average.Cp      =       PP.LiquidCp(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);
80
81"Cold Stream Inlet Heat Capacity"
82        Properties.Cold.Inlet.Cp        =       PP.LiquidCp(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);
83
84"Cold Stream Outlet Heat Capacity"
85        Properties.Cold.Outlet.Cp       =       PP.LiquidCp(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);
86
87"Cold Stream Average Mass Density"
88        Properties.Cold.Average.rho =   PP.LiquidDensity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);
89
90"Cold Stream Inlet Mass Density"
91        Properties.Cold.Inlet.rho       =       PP.LiquidDensity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);
92
93"Cold Stream Outlet Mass Density"
94        Properties.Cold.Outlet.rho      =       PP.LiquidDensity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);
95
96"Cold Stream Average Viscosity"
97        Properties.Cold.Average.Mu      =       PP.LiquidViscosity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);
98
99"Cold Stream inlet Viscosity"
100        Properties.Cold.Inlet.Mu        =       PP.LiquidViscosity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);
101       
102"Cold Stream Outlet Viscosity"
103        Properties.Cold.Outlet.Mu       =       PP.LiquidViscosity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);
104
105"Cold Stream Average Conductivity"
106        Properties.Cold.Average.K       =       PP.LiquidThermalConductivity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);
107
108"Cold Stream Inlet Conductivity"       
109        Properties.Cold.Inlet.K         =       PP.LiquidThermalConductivity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);
110
111"Cold Stream Outlet Conductivity"
112        Properties.Cold.Outlet.K        =       PP.LiquidThermalConductivity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);
113
114"Cold Stream Heat Capacity at Wall Temperature"
115        Properties.Cold.Wall.Cp         =       PP.LiquidCp(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z);
116       
117"Cold Stream Viscosity at Wall Temperature"
118        Properties.Cold.Wall.Mu         =       PP.LiquidViscosity(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z);
119
120"Cold Stream Conductivity at Wall Temperature"
121        Properties.Cold.Wall.K          =       PP.LiquidThermalConductivity(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z);
122
123
124        else
125
126"Cold Stream Average Heat Capacity"
127        Properties.Cold.Average.Cp      =       PP.VapourCp(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);
128
129"Cold Stream Inlet Heat Capacity"       
130        Properties.Cold.Inlet.Cp        =       PP.VapourCp(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);
131
132"Cold Stream Outlet Heat Capacity"     
133        Properties.Cold.Outlet.Cp       =       PP.VapourCp(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);
134
135"Cold Stream Average Mass Density"
136        Properties.Cold.Average.rho =   PP.VapourDensity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);
137
138"Cold Stream Inlet Mass Density"
139        Properties.Cold.Inlet.rho       =       PP.VapourDensity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);
140
141"Cold Stream Outlet Mass Density"       
142        Properties.Cold.Outlet.rho      =       PP.VapourDensity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);
143
144"Cold Stream Average Viscosity "
145        Properties.Cold.Average.Mu      =       PP.VapourViscosity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);
146
147"Cold Stream Inlet Viscosity " 
148        Properties.Cold.Inlet.Mu        =       PP.VapourViscosity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);
149
150"Cold Stream Outlet Viscosity "
151        Properties.Cold.Outlet.Mu       =       PP.VapourViscosity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);
152
153"Cold Stream Average Conductivity "
154        Properties.Cold.Average.K       =       PP.VapourThermalConductivity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);
155
156"Cold Stream Inlet Conductivity "
157        Properties.Cold.Inlet.K         =       PP.VapourThermalConductivity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);
158
159"Cold Stream Outlet Conductivity "
160        Properties.Cold.Outlet.K        =       PP.VapourThermalConductivity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);
161       
162"Cold Stream Heat Capacity at Wall Temperature"
163        Properties.Cold.Wall.Cp         =       PP.VapourCp(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z);
164
165
166"Cold Stream Viscosity at Wall Temperature"
167        Properties.Cold.Wall.Mu         =       PP.VapourViscosity(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z);
168
169"Cold Stream Conductivity at Wall Temperature"
170        Properties.Cold.Wall.K          =       PP.VapourThermalConductivity(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z);
171       
172       
173       
174end
175
176if Inlet.Hot.v equal 0
177
178        then
179
180"Hot Stream Average Heat Capacity"
181        Properties.Hot.Average.Cp       =               PP.LiquidCp(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z);
182
183"Hot Stream Inlet Heat Capacity"
184        Properties.Hot.Inlet.Cp         =               PP.LiquidCp(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);
185
186"Hot Stream Outlet Heat Capacity"
187        Properties.Hot.Outlet.Cp        =               PP.LiquidCp(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z);
188
189"Hot Stream Average Mass Density"
190        Properties.Hot.Average.rho      =               PP.LiquidDensity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z);
191
192"Hot Stream Inlet Mass Density"
193        Properties.Hot.Inlet.rho        =               PP.LiquidDensity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);
194
195"Hot Stream Outlet Mass Density"       
196        Properties.Hot.Outlet.rho       =               PP.LiquidDensity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z);
197
198"Hot Stream Average Viscosity"
199        Properties.Hot.Average.Mu       =               PP.LiquidViscosity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z);     
200
201"Hot Stream Inlet Viscosity"
202        Properties.Hot.Inlet.Mu         =               PP.LiquidViscosity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);       
203
204"Hot Stream Outlet Viscosity"
205        Properties.Hot.Outlet.Mu        =               PP.LiquidViscosity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z);     
206
207"Hot Stream Average Conductivity"
208        Properties.Hot.Average.K        =               PP.LiquidThermalConductivity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z);   
209
210"Hot Stream Inlet Conductivity"
211        Properties.Hot.Inlet.K          =               PP.LiquidThermalConductivity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);     
212
213"Hot Stream Outlet Conductivity"
214        Properties.Hot.Outlet.K         =               PP.LiquidThermalConductivity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z);   
215
216"Hot Stream Heat Capacity at Wall Temperature"
217        Properties.Hot.Wall.Cp          =               PP.LiquidCp(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z);
218
219"Hot Stream Viscosity  at Wall Temperature"
220        Properties.Hot.Wall.Mu          =               PP.LiquidViscosity(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z);     
221
222"Hot Stream Conductivity at Wall Temperature"
223        Properties.Hot.Wall.K           =               PP.LiquidThermalConductivity(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z);   
224       
225
226        else
227
228"Hot Stream Average Heat Capacity"
229        Properties.Hot.Average.Cp       =               PP.VapourCp(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z);
230
231"Hot Stream Inlet Heat Capacity"
232        Properties.Hot.Inlet.Cp         =               PP.VapourCp(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);
233
234"Hot Stream Outlet Heat Capacity"
235        Properties.Hot.Outlet.Cp        =               PP.VapourCp(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z);
236
237"Hot Stream Average Mass Density"
238        Properties.Hot.Average.rho      =               PP.VapourDensity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z);
239
240"Hot Stream Inlet Mass Density"
241        Properties.Hot.Inlet.rho        =               PP.VapourDensity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);
242
243"Hot Stream Outlet Mass Density"
244        Properties.Hot.Outlet.rho       =               PP.VapourDensity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z);
245
246"Hot Stream Average Viscosity"
247        Properties.Hot.Average.Mu       =               PP.VapourViscosity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z);
248
249"Hot Stream Inlet Viscosity"
250        Properties.Hot.Inlet.Mu         =               PP.VapourViscosity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);
251
252"Hot Stream Outlet Viscosity"
253        Properties.Hot.Outlet.Mu        =               PP.VapourViscosity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z);
254
255"Hot Stream Average Conductivity"
256        Properties.Hot.Average.K        =               PP.VapourThermalConductivity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z);   
257
258"Hot Stream Inlet Conductivity"
259        Properties.Hot.Inlet.K          =               PP.VapourThermalConductivity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);     
260       
261"Hot Stream Outlet Conductivity"
262        Properties.Hot.Outlet.K         =               PP.VapourThermalConductivity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z);   
263
264"Hot Stream Heat Capacity at Wall Temperature"
265        Properties.Hot.Wall.Cp          =               PP.VapourCp(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z);
266
267"Hot Stream Viscosity at Wall Temperature"
268        Properties.Hot.Wall.Mu          =               PP.VapourViscosity(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z);
269
270"Hot Stream Conductivity at Wall Temperature"
271        Properties.Hot.Wall.K           =               PP.VapourThermalConductivity(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z);   
272
273
274end
275
276#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
277#                                       Thermal Details                                                                 #
278#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
279
280"Hot Stream Heat Capacity"
281        Details.Ch =Inlet.Hot.F*Properties.Hot.Average.Cp;
282       
283"Cold Stream Heat Capacity"
284        Details.Cc =Inlet.Cold.F*Properties.Cold.Average.Cp;
285       
286"Minimum Heat Capacity"
287        Details.Cmin  = min([Details.Ch,Details.Cc]);
288       
289"Maximum Heat Capacity"
290        Details.Cmax  = max([Details.Ch,Details.Cc]);
291
292"Heat Capacity Ratio"
293        Details.Cr    = Details.Cmin/Details.Cmax;
294
295#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
296#                                       Energy Balance
297#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
298
299"Energy Balance Hot Stream"
300        Details.Q = Inlet.Hot.F*(Inlet.Hot.h-Outlet.Hot.h);
301
302"Energy Balance Cold Stream"
303        Details.Q =-Inlet.Cold.F*(Inlet.Cold.h-Outlet.Cold.h);
304
305#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
306#                                       Material Balance
307#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
308"Flow Mass Inlet Cold Stream"
309        Properties.Cold.Inlet.Fw        =  sum(M*Inlet.Cold.z)*Inlet.Cold.F;
310
311"Flow Mass Outlet Cold Stream"
312        Properties.Cold.Outlet.Fw       =  sum(M*Outlet.Cold.z)*Outlet.Cold.F;
313
314"Flow Mass Inlet Hot Stream"
315        Properties.Hot.Inlet.Fw         =  sum(M*Inlet.Hot.z)*Inlet.Hot.F;
316
317"Flow Mass Outlet Hot Stream"   
318        Properties.Hot.Outlet.Fw        =  sum(M*Outlet.Hot.z)*Outlet.Hot.F;
319
320"Molar Balance Hot Stream"
321        Inlet.Hot.F  = Outlet.Hot.F;
322       
323"Molar Balance Cold Stream"
324        Inlet.Cold.F = Outlet.Cold.F;
325
326#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
327#                                       Constraints
328#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
329"Hot Stream Molar Fraction Constraint"
330        Outlet.Hot.z=Inlet.Hot.z;
331       
332"Cold Stream Molar Fraction Constraint"
333        Outlet.Cold.z=Inlet.Cold.z;
334       
335"No Phase Change In Cold Stream"
336        Inlet.Cold.v=Outlet.Cold.v;
337
338"No Phase Change In Hot Stream"
339        Inlet.Hot.v=Outlet.Hot.v;
340       
341end
342
343Model Heatex_Detailed           as HeatExchangerDetailed_Basic
344
345PARAMETERS
346
347#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
348#                               Basic  Parameters                                                                       #
349#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
350ext PP          as CalcObject   (Brief="External Physical Properties");
351side        as Integer          (Brief="Fluid Alocation Flag",Lower=0,Upper=1);
352Pi                      as constant     (Brief="Pi Number",Default=3.14159265);
353#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
354#                               Shell Geometrical Parameters                                            #
355#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
356Tpass                   as Integer              (Brief="Number of Tube Passes",Lower=1);
357Nss                             as Integer              (Brief="Number of Sealing Strips pairs",Lower=1);
358Dishell                 as length               (Brief="Inside Shell Diameter",Lower=10e-6);
359Donozzle_Shell  as length               (Brief="Shell Outlet Nozzle Diameter",Lower=10e-6);
360Dinozzle_Shell  as length               (Brief="Shell Inlet Nozzle Diameter",Lower=10e-6);
361Hinozzle_Shell  as length               (Brief="Height Under Shell Inlet Nozzle",Lower=10e-6);
362Honozzle_Shell  as length               (Brief="Height Under Shell Outlet Nozzle",Lower=10e-6);
363Lcf                     as length               (Brief="Bundle-to-Shell Clearance",Lower=10e-8);
364
365#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
366#                               Tubes Geometrical Parameters                                            #
367#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
368Ntt                       as Integer            (Brief="Total Number of Tubes in Shell",Default=100,Lower=1);
369Pattern       as Integer                (Brief="Tube Layout Characteristic Angle",Lower=30);
370Ltube             as length                     (Brief="Effective Tube Length",Lower=0.1);
371pitch             as length                     (Brief="Tube Pitch",Lower=1e-8);
372Kwall             as conductivity       (Brief="Tube Wall Material Thermal Conductivity");
373Dotube            as length                     (Brief="Tube Outside Diameter",Lower=10e-6);
374Ditube            as length                     (Brief="Tube Inside Diameter",Lower=10e-6);
375Donozzle_Tube as length                 (Brief="Tube Outlet Nozzle Diameter",Lower=10e-6);
376Dinozzle_Tube as length                 (Brief="Tube Inlett Nozzle Diameter",Lower=10e-6);
377
378#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
379#                               Baffles Geometrical Parameters                                          #
380#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
381Bc     as Integer (Brief="Baffle Cut",Default=25,Lower=25);
382Nb     as Real    (Brief="Number of Baffles",Lower=1);
383Lcd    as length  (Brief="Baffle-to-Shell Clearance",Lower=10e-8);
384Ltd        as length  (Brief="Tube-to-Bafflehole Clearance",Lower=10e-8);
385
386SET
387
388#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
389#             Set Parameters for heatex Calculation                                     #
390#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
391        Pi                              = 3.14159265;
392        HE.Tpass                        = Tpass;
393        HE.Nss                          = Nss;
394        HE.Ntt                          = Ntt;
395        HE.Pattern                      = Pattern;
396        HE.Bc                           = Bc;
397        HE.Donozzle_Shell       = Donozzle_Shell;
398        HE.Dinozzle_Shell       = Dinozzle_Shell;
399        HE.Honozzle_Shell       = Honozzle_Shell;
400        HE.Hinozzle_Shell       = Hinozzle_Shell;
401        HE.Donozzle_Tube        = Donozzle_Tube;
402        HE.Dinozzle_Tube        = Dinozzle_Tube;
403        HE.Nb                       = Nb;
404        HE.Dishell                      = Dishell;             
405        HE.Lcf                          = Lcf; 
406        HE.Ltube                        = Ltube;       
407        HE.pitch                        = pitch;               
408        HE.Dotube                       = Dotube;       
409        HE.Ditube                       = Ditube;       
410        HE.Lcd                          = Lcd;
411        HE.Ltd                          = Ltd;
412        side                            = HE.FluidAlocation();
413       
414EQUATIONS
415
416#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
417#                               Heat Transfer Correction Factors                                        #
418#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
419"Ji Factor"
420        Shell.HeatTransfer.Ji = HE.JiFactor(Shell.HeatTransfer.Re);
421
422"Jc Factor"
423        Shell.HeatTransfer.Jc = HE.JcFactor();
424       
425"Jl Factor"
426        Shell.HeatTransfer.Jl = HE.JlFactor(Shell.HeatTransfer.Sm);
427
428"Jb Factor"
429        Shell.HeatTransfer.Jb = HE.JbFactor(Shell.HeatTransfer.Re,Baffles.Ls,Shell.HeatTransfer.Sm);
430
431"Jr Factor"     
432        Shell.HeatTransfer.Jr = HE.JrFactor(Shell.HeatTransfer.Re);
433       
434"Total J Factor"
435        Shell.HeatTransfer.Jtotal =     Shell.HeatTransfer.Jc*Shell.HeatTransfer.Jl*Shell.HeatTransfer.Jb*Shell.HeatTransfer.Jr*Shell.HeatTransfer.Js;
436
437#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
438#                                        Nozzles rho-V^2                                                                #
439#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
440
441if side equal 1
442       
443        then
444
445"Shell Side inlet Nozzle rho-V^2"
446        Shell.PressureDrop.RVsquare_in = Properties.Hot.Inlet.rho*(Shell.PressureDrop.Vnozzle_in)^2;
447
448"Shell Side Outlet Nozzle rho-V^2"
449        Shell.PressureDrop.RVsquare_out = Properties.Hot.Outlet.rho*(Shell.PressureDrop.Vnozzle_out)^2;
450       
451        else
452       
453"Shell Side inlet Nozzle rho-V^2"
454        Shell.PressureDrop.RVsquare_in = Properties.Cold.Inlet.rho*(Shell.PressureDrop.Vnozzle_in)^2;
455
456"Shell Side Outlet Nozzle rho-V^2"
457        Shell.PressureDrop.RVsquare_out = Properties.Cold.Outlet.rho*(Shell.PressureDrop.Vnozzle_out)^2;
458
459end
460
461#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
462#                                        Phi correction                                                                 #
463#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
464
465if side equal 1
466       
467        then
468
469"Shell Side Phi correction"
470        Shell.HeatTransfer.Phi  = HE.PhiCorrection(Properties.Hot.Average.Mu,Properties.Hot.Wall.Mu);
471       
472"Tube Side Phi correction"
473        Tubes.HeatTransfer.Phi  = HE.PhiCorrection(Properties.Cold.Average.Mu,Properties.Cold.Wall.Mu);
474       
475        else
476       
477"Shell Side Phi correction"
478        Shell.HeatTransfer.Phi  = HE.PhiCorrection(Properties.Cold.Average.Mu,Properties.Cold.Wall.Mu);
479       
480"Tube Side Phi correction"
481        Tubes.HeatTransfer.Phi  = HE.PhiCorrection(Properties.Hot.Average.Mu,Properties.Hot.Wall.Mu);
482
483end
484
485#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
486#                                                                                                                                       #
487#                                                                                                                                       #
488#                                                                                                                                       #
489#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
490
491if side equal 1
492       
493        then
494       
495"Hot Wall Temperature"
496        Properties.Hot.Wall.Twall       = (Properties.Hot.Average.T+Properties.Cold.Average.T)/2;
497
498"ColdWall Temperature"
499        Properties.Cold.Wall.Twall  = (Properties.Hot.Average.T+Properties.Cold.Average.T)/2;
500
501"Tube Side Velocity"
502        Tubes.HeatTransfer.Vtube        = Properties.Cold.Inlet.Fw*Tpass/((Pi*Ditube*Ditube/4)*Properties.Cold.Average.rho*Ntt);
503
504"Tube Side Reynolds Number"
505        Tubes.HeatTransfer.Re           = (Properties.Cold.Average.rho*Tubes.HeatTransfer.Vtube*Ditube)/Properties.Cold.Average.Mu;
506       
507"Tube Side Prandtl Number"
508        Tubes.HeatTransfer.PR           = ((Properties.Cold.Average.Cp/Properties.Cold.Average.Mw)*Properties.Cold.Average.Mu)/Properties.Cold.Average.K;
509
510"Tube Side Prandtl Number at Wall Temperature"
511        Tubes.HeatTransfer.PRw          = ((Properties.Cold.Wall.Cp/Properties.Cold.Average.Mw)*Properties.Cold.Wall.Mu)/Properties.Cold.Wall.K;
512
513"Tube Side Film Coefficient"
514        Tubes.HeatTransfer.htube        = HE.TubeFilmCoeff(Tubes.HeatTransfer.Re,Tubes.HeatTransfer.PR,Properties.Cold.Average.K)*Tubes.HeatTransfer.Phi;
515
516"Shell Side Prandtl Number"
517        Shell.HeatTransfer.PR           = ((Properties.Hot.Average.Cp/Properties.Hot.Average.Mw)*Properties.Hot.Average.Mu)/Properties.Hot.Average.K;
518
519"Shell Side Prandtl Number at Wall Temperature"
520        Shell.HeatTransfer.PRw          = ((Properties.Hot.Wall.Cp/Properties.Hot.Average.Mw)*Properties.Hot.Wall.Mu)/Properties.Hot.Wall.K;
521
522"Tube Side Pressure Drop"
523        Tubes.PressureDrop.PdTube               = HE.DeltaPtube(Tubes.HeatTransfer.Re,Properties.Cold.Average.rho,Tubes.HeatTransfer.Vtube,Tubes.HeatTransfer.Phi);
524
525"Pressure Drop Tube Side Inlet Nozzle"
526        Tubes.PressureDrop.Pdnozzle_in  = HE.DeltaPtubeNozzlein(Properties.Cold.Inlet.rho,Properties.Cold.Inlet.Fw);
527
528"Velocity Tube Side Inlet Nozzle"
529        Tubes.PressureDrop.Vnozzle_in   = HE.TubeVelocityNozzlein(Properties.Cold.Inlet.rho,Properties.Cold.Inlet.Fw);
530
531"Pressure Drop Tube Side Outlet Nozzle"
532        Tubes.PressureDrop.Pdnozzle_out = HE.DeltaPtubeNozzleout(Properties.Cold.Outlet.rho,Properties.Cold.Inlet.Fw);
533
534"Velocity Tube Side Outlet Nozzle"
535        Tubes.PressureDrop.Vnozzle_out  = HE.TubeVelocityNozzleout(Properties.Cold.Outlet.rho,Properties.Cold.Inlet.Fw);
536
537"Shell Pressure Drop Inlet Nozzle"
538        Shell.PressureDrop.Pdnozzle_in  = HE.DeltaPshellNozzleIn(Properties.Hot.Inlet.rho,Properties.Hot.Inlet.Fw);
539
540"Velocity Shell Side Inlet Nozzle"
541        Shell.PressureDrop.Vnozzle_in   = HE.ShellVelocityNozzleIn(Properties.Hot.Inlet.rho,Properties.Hot.Inlet.Fw);
542
543"Shell Pressure Drop Outlet Nozzle"
544        Shell.PressureDrop.Pdnozzle_out = HE.DeltaPshellNozzleOut(Properties.Hot.Outlet.rho,Properties.Hot.Inlet.Fw);
545
546"Velocity Shell Side Outlet Nozzle"
547        Shell.PressureDrop.Vnozzle_out  =HE.ShellVelocityNozzleOut(Properties.Hot.Outlet.rho,Properties.Hot.Inlet.Fw);
548
549"Pressure Drop Hot Stream"
550        Outlet.Hot.P  = Inlet.Hot.P- Shell.PressureDrop.Pdtotal;       
551       
552"Pressure Drop Cold Stream"
553        Outlet.Cold.P  = Inlet.Cold.P - Tubes.PressureDrop.Pdtotal;
554
555
556        else
557
558"Hot Wall Temperature"
559        Properties.Hot.Wall.Twall  = (Properties.Hot.Average.T+Properties.Cold.Average.T)/2;
560
561"Cold Wall Temperature"
562        Properties.Cold.Wall.Twall  = (Properties.Hot.Average.T+Properties.Cold.Average.T)/2;
563
564"Tube Side Velocity"
565        Tubes.HeatTransfer.Vtube  = Properties.Hot.Inlet.Fw*Tpass/((Pi*Ditube*Ditube/4)*Properties.Hot.Average.rho*Ntt);
566
567"Tube Side Reynolds Number"
568        Tubes.HeatTransfer.Re=(Properties.Hot.Average.rho*Tubes.HeatTransfer.Vtube*Ditube)/Properties.Hot.Average.Mu;
569
570"Tube Side Prandtl Number"
571        Tubes.HeatTransfer.PR = ((Properties.Hot.Average.Cp/Properties.Hot.Average.Mw)*Properties.Hot.Average.Mu)/Properties.Hot.Average.K;
572
573"Tube Side Prandtl Number at Wall"
574        Tubes.HeatTransfer.PRw = ((Properties.Hot.Wall.Cp/Properties.Hot.Average.Mw)*Properties.Hot.Wall.Mu)/Properties.Hot.Wall.K;
575
576"Tube Side Film Coefficient"
577        Tubes.HeatTransfer.htube= HE.TubeFilmCoeff(Tubes.HeatTransfer.Re,Tubes.HeatTransfer.PR,Properties.Hot.Average.K)*Tubes.HeatTransfer.Phi;
578       
579"Shell Side Prandtl Number"
580        Shell.HeatTransfer.PR = ((Properties.Cold.Average.Cp/Properties.Cold.Average.Mw)*Properties.Cold.Average.Mu)/Properties.Cold.Average.K;
581
582"Shell Side Prandtl Number at Wall"
583        Shell.HeatTransfer.PRw=((Properties.Cold.Wall.Cp/Properties.Cold.Average.Mw)*Properties.Cold.Wall.Mu)/Properties.Cold.Wall.K;
584
585"Tube Side Pressure Drop"
586        Tubes.PressureDrop.PdTube               = HE.DeltaPtube(Tubes.HeatTransfer.Re,Properties.Hot.Average.rho,Tubes.HeatTransfer.Vtube,Tubes.HeatTransfer.Phi);
587
588"Pressure Drop Tube Side Inlet Nozzle"
589        Tubes.PressureDrop.Pdnozzle_in  = HE.DeltaPtubeNozzlein(Properties.Hot.Inlet.rho,Properties.Hot.Inlet.Fw);
590
591"Velocity Tube Side Inlet Nozzle"
592        Tubes.PressureDrop.Vnozzle_in   = HE.TubeVelocityNozzlein(Properties.Hot.Inlet.rho,Properties.Hot.Inlet.Fw);
593
594"Pressure Drop Tube Side Outlet Nozzle"
595        Tubes.PressureDrop.Pdnozzle_out = HE.DeltaPtubeNozzleout(Properties.Hot.Outlet.rho,Properties.Hot.Inlet.Fw);   
596
597"Velocity Tube Side Outlet Nozzle"
598        Tubes.PressureDrop.Vnozzle_out  = HE.TubeVelocityNozzleout(Properties.Hot.Outlet.rho,Properties.Hot.Inlet.Fw); 
599
600"Shell Pressure Drop Inlet Nozzle"
601        Shell.PressureDrop.Pdnozzle_in  = HE.DeltaPshellNozzleIn(Properties.Cold.Inlet.rho,Properties.Cold.Inlet.Fw);
602
603"Velocity Shell Side Inlet Nozzle"
604        Shell.PressureDrop.Vnozzle_in   = HE.ShellVelocityNozzleIn(Properties.Cold.Inlet.rho,Properties.Cold.Inlet.Fw);
605
606"Shell Pressure Drop Outlet Nozzle"
607        Shell.PressureDrop.Pdnozzle_out = HE.DeltaPshellNozzleOut(Properties.Cold.Outlet.rho,Properties.Cold.Inlet.Fw);
608
609"Velocity Shell Side Outlet Nozzle"
610        Shell.PressureDrop.Vnozzle_out  = HE.ShellVelocityNozzleOut(Properties.Cold.Outlet.rho,Properties.Cold.Inlet.Fw);
611
612"Pressure Drop Hot Stream"
613        Outlet.Hot.P  = Inlet.Hot.P- Tubes.PressureDrop.Pdtotal;       
614       
615"Pressure Drop Cold Stream"
616        Outlet.Cold.P  = Inlet.Cold.P - Shell.PressureDrop.Pdtotal;
617
618end
619
620#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
621#                                               Resistances                                                                     #
622#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
623
624"Tube Resistance"       
625        Resistances.Rtube*(Tubes.HeatTransfer.htube*Ditube) = Dotube;
626       
627"Wall Resistance"
628        Resistances.Rwall=Dotube*ln(Dotube/Ditube)/(2*Kwall);
629       
630"Shell Resistance"
631        Resistances.Rshell*(Shell.HeatTransfer.hshell)=1;
632
633#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
634#                                                                                                                                       #
635#                                                                                                                                       #
636#                                                                                                                                       #
637#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
638
639"Shell Side Cross Flow Area"
640        Shell.HeatTransfer.Sm = HE.CrossFlowArea(Baffles.Ls);
641       
642"Overall Heat Transfer Coefficient Dirty"
643        Details.Ud=1/(Dotube/(Tubes.HeatTransfer.htube*Ditube)+Resistances.Rfo+Resistances.Rfi*(Dotube/Ditube)+(Dotube*ln(Dotube/Ditube)/(2*Kwall))+(1/(Shell.HeatTransfer.hshell)));
644
645"Overall Heat Transfer Coefficient Clean"
646        (1/Details.Ud)=(1/Details.Uc)+Resistances.Rfo+Resistances.Rfi*(Dotube/Ditube);
647
648"Exchange Surface Area"
649        Details.A=Pi*Dotube*Ntt*Ltube;
650
651"Baffles Spacing"
652        Ltube = Baffles.Lsi+Baffles.Lso+Baffles.Ls*(Nb-1);
653
654end
655
656Model Heatex_Detailed_NTU               as Heatex_Detailed
657       
658VARIABLES
659
660Eft                      as positive (Brief="Effectiveness",Default=0.05,Lower=1e-8);
661       
662EQUATIONS
663"Duty"
664        Details.Q       = Eft*Details.Cmin*(Inlet.Hot.T-Inlet.Cold.T); 
665       
666end
667
668Model Heatex_Detailed_LMTD              as Heatex_Detailed
669       
670VARIABLES
671
672DT0     as temp_delta   (Brief="Temperature Difference at Inlet",Lower=1);
673DTL             as temp_delta   (Brief="Temperature Difference at Outlet",Lower=1);
674LMTD    as temp_delta   (Brief="Logarithmic Mean Temperature Difference",Lower=5);
675Fc              as positive             (Brief="LMTD Correction Factor",Lower=0.4);
676MTD             as temp_delta   (Brief="Mean Temperature Difference",Lower=5);
677
678EQUATIONS
679#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
680#                       Log Mean Temperature Difference
681#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
682if abs(DT0 - DTL) > 0.05*max(abs([DT0,DTL]))
683       
684        then
685"Log Mean Temperature Difference"
686        LMTD= (DT0-DTL)/ln(DT0/DTL);
687
688        else
689       
690if DT0*DTL equal 0
691       
692        then
693"Log Mean Temperature Difference"
694        LMTD = 0.5*(DT0+DTL);
695       
696        else
697"Log Mean Temperature Difference"
698        LMTD = 0.5*(DT0+DTL)*(1-(DT0-DTL)^2/(DT0*DTL)*(1+(DT0-DTL)^2/(DT0*DTL)/2)/12);
699       
700end
701       
702end
703
704"Exchange Surface Area"
705        Details.Q   = Details.Ud*Details.A*MTD;
706       
707"Mean Temperature Difference"   
708        MTD = Fc*LMTD;
709       
710end
711
712Model E_Shell_NTU_Det                   as Heatex_Detailed_NTU
713#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
714#       Shell and Tubes Heat Exchanger with 1 shell pass - LMTD Method
715#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
716EQUATIONS
717
718"TEMA E Shell Effectiveness"
719#       Eft = 2*(1+Details.Cr+sqrt(1+Details.Cr^2)*((1+exp(-Details.NTU*sqrt(1+Details.Cr^2)))/(1-exp(-Details.NTU*sqrt(1+Details.Cr^2)))) )^-1;
720        Eft = HE.EshellEffectiveness(Details.Cr,Details.NTU);
721
722"Js Factor"     
723        Shell.HeatTransfer.Js = HE.JsFactor(Shell.HeatTransfer.Re,Baffles.Lsi,Baffles.Lso,Baffles.Ls);
724
725if side equal 1
726       
727        then
728       
729"Shell Side Reynolds Number"
730        Shell.HeatTransfer.Re=(Dotube*Properties.Hot.Inlet.Fw/Shell.HeatTransfer.Sm)/Properties.Hot.Average.Mu;
731
732"Shell Heat Transfer Coefficient"
733        Shell.HeatTransfer.hshell=Shell.HeatTransfer.Ji*(Properties.Hot.Average.Cp/Properties.Hot.Average.Mw)*(Properties.Hot.Inlet.Fw/Shell.HeatTransfer.Sm)*(Shell.HeatTransfer.PR^(-2/3))*Shell.HeatTransfer.Jtotal*Shell.HeatTransfer.Phi;
734
735"Shell Pressure Drop Cross Flow"
736        Shell.PressureDrop.PdCross              = HE.DeltaPcross(Shell.HeatTransfer.Re,Baffles.Ls,Baffles.Lso,Baffles.Lsi,Properties.Hot.Inlet.Fw,Shell.HeatTransfer.Phi,Properties.Hot.Average.rho);
737
738"Shell Pressure Baffle Window"
739        Shell.PressureDrop.Pdwindow             = HE.DeltaPwindow(Properties.Hot.Inlet.Fw,Shell.HeatTransfer.Sm,Properties.Hot.Average.rho,Properties.Hot.Average.Mu,Baffles.Ls);
740
741"Shell Pressure End Zones"
742        Shell.PressureDrop.PdEndZones  = HE.DeltaPendZones(Shell.HeatTransfer.Re,Baffles.Ls,Baffles.Lso,Baffles.Lsi,Properties.Hot.Inlet.Fw,Shell.HeatTransfer.Phi,Properties.Hot.Average.rho);
743
744
745        else
746
747"Shell Side Reynolds Number"
748        Shell.HeatTransfer.Re=(Dotube*Properties.Cold.Inlet.Fw/Shell.HeatTransfer.Sm)/Properties.Cold.Average.Mu;
749
750"Shell Heat Transfer Coefficient"
751        Shell.HeatTransfer.hshell=Shell.HeatTransfer.Ji*(Properties.Cold.Average.Cp/Properties.Cold.Average.Mw)*(Properties.Cold.Inlet.Fw/Shell.HeatTransfer.Sm)*(Shell.HeatTransfer.PR^(-2/3))*Shell.HeatTransfer.Jtotal*Shell.HeatTransfer.Phi;
752
753"Shell Pressure Drop Cross Flow"
754        Shell.PressureDrop.PdCross      = HE.DeltaPcross(Shell.HeatTransfer.Re,Baffles.Ls,Baffles.Lso,Baffles.Lsi,Properties.Cold.Inlet.Fw,Shell.HeatTransfer.Phi,Properties.Cold.Average.rho);
755
756"Shell Pressure Baffle Window"
757        Shell.PressureDrop.Pdwindow     = HE.DeltaPwindow(Properties.Cold.Inlet.Fw,Shell.HeatTransfer.Sm,Properties.Cold.Average.rho,Properties.Cold.Average.Mu,Baffles.Ls);
758
759"Shell Pressure End Zones"
760        Shell.PressureDrop.PdEndZones  = HE.DeltaPendZones(Shell.HeatTransfer.Re,Baffles.Ls,Baffles.Lso,Baffles.Lsi,Properties.Cold.Inlet.Fw,Shell.HeatTransfer.Phi,Properties.Cold.Average.rho);
761
762
763end
764
765end
766
767Model F_Shell_NTU_Det                   as Heatex_Detailed_NTU
768#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
769#       Shell and Tubes Heat Exchanger with 1 shell pass - LMTD Method
770#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
771VARIABLES
772
773Eft1    as positive (Brief="Effectiveness Correction",Lower=0.01,Default=0.5);
774
775EQUATIONS
776
777"Effectiveness Correction"
778        Eft1 = 2*(1+Details.Cr+sqrt(1+Details.Cr^2)*((1+exp(-Details.NTU*sqrt(1+Details.Cr^2)))/(1-exp(-Details.NTU*sqrt(1+Details.Cr^2)))) )^-1;
779
780"TEMA F Shell Effectiveness"
781        Eft = ( ((1-Eft1*Details.Cr)/(1-Eft1))^2 -1  )*( ((1-Eft1*Details.Cr)/(1-Eft1))^2 - Details.Cr )^-1;
782
783"Js Factor"     
784        Shell.HeatTransfer.Js = HE.JsFactor(Shell.HeatTransfer.Re,Baffles.Lsi,Baffles.Lso,Baffles.Ls);
785
786if side equal 1
787       
788        then
789       
790"Shell Side Reynolds Number"
791        Shell.HeatTransfer.Re=(Dotube*Properties.Hot.Inlet.Fw/Shell.HeatTransfer.Sm)/Properties.Hot.Average.Mu;
792
793"Shell Heat Transfer Coefficient"
794        Shell.HeatTransfer.hshell= Shell.HeatTransfer.Ji*(Properties.Hot.Average.Cp/Properties.Hot.Average.Mw)*(Properties.Hot.Inlet.Fw/Shell.HeatTransfer.Sm)*(Shell.HeatTransfer.PR^(-2/3))*Shell.HeatTransfer.Jtotal*Shell.HeatTransfer.Phi;
795
796"Shell Pressure Drop Cross Flow"
797        Shell.PressureDrop.PdCross              = HE.DeltaPcross(Shell.HeatTransfer.Re,Baffles.Ls,Baffles.Lso,Baffles.Lsi,Properties.Hot.Inlet.Fw,Shell.HeatTransfer.Phi,Properties.Hot.Average.rho);
798
799"Shell Pressure Baffle Window"
800        Shell.PressureDrop.Pdwindow             = HE.DeltaPwindow(Properties.Hot.Inlet.Fw,Shell.HeatTransfer.Sm,Properties.Hot.Average.rho,Properties.Hot.Average.Mu,Baffles.Ls);
801
802"Shell Pressure End Zones"
803        Shell.PressureDrop.PdEndZones  = HE.DeltaPendZones(Shell.HeatTransfer.Re,Baffles.Ls,Baffles.Lso,Baffles.Lsi,Properties.Hot.Inlet.Fw,Shell.HeatTransfer.Phi,Properties.Hot.Average.rho);
804
805
806        else
807
808"Shell Side Reynolds Number"
809        Shell.HeatTransfer.Re=(Dotube*Properties.Cold.Inlet.Fw/Shell.HeatTransfer.Sm)/Properties.Cold.Average.Mu;
810
811"Shell Heat Transfer Coefficient"
812        Shell.HeatTransfer.hshell=Shell.HeatTransfer.Ji*(Properties.Cold.Average.Cp/Properties.Cold.Average.Mw)*(Properties.Cold.Inlet.Fw/Shell.HeatTransfer.Sm)*(Shell.HeatTransfer.PR^(-2/3))*Shell.HeatTransfer.Jtotal*Shell.HeatTransfer.Phi;
813
814"Shell Pressure Drop Cross Flow"
815        Shell.PressureDrop.PdCross      = HE.DeltaPcross(Shell.HeatTransfer.Re,Baffles.Ls,Baffles.Lso,Baffles.Lsi,Properties.Cold.Inlet.Fw,Shell.HeatTransfer.Phi,Properties.Cold.Average.rho);
816
817"Shell Pressure Baffle Window"
818        Shell.PressureDrop.Pdwindow     = HE.DeltaPwindow(Properties.Cold.Inlet.Fw,Shell.HeatTransfer.Sm,Properties.Cold.Average.rho,Properties.Cold.Average.Mu,Baffles.Ls);
819
820"Shell Pressure End Zones"
821        Shell.PressureDrop.PdEndZones  = HE.DeltaPendZones(Shell.HeatTransfer.Re,Baffles.Ls,Baffles.Lso,Baffles.Lsi,Properties.Cold.Inlet.Fw,Shell.HeatTransfer.Phi,Properties.Cold.Average.rho);
822
823
824end
825
826
827end
828
829Model Multipass_NTU_Det         
830#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
831#       Shell and Tubes Heat Exchanger In Series with 1 shell pass - LMTD Method
832#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
833PARAMETERS
834
835        Nshell          as Integer              (Brief="N Shell in Series",Default=2);
836ext PP                  as CalcObject   (Brief="External Physical Properties");
837        HE                      as CalcObject   (Brief="heatex Calculations", File="heatex");
838        side        as Integer          (Brief="Fluid Alocation",Lower=0,Upper=1);
839        Pi                      as constant     (Brief="Pi Number",Default=3.14159265);
840
841#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
842#       Shell Geometrical Parameters
843#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
844Tpass                   as Integer              (Brief="Number of Tube Passes",Lower=1);
845Nss                             as Integer              (Brief="Number of Sealing Strips pairs",Lower=1);
846Dishell                 as length               (Brief="Inside Shell Diameter",Lower=10e-6);
847Donozzle_Shell  as length               (Brief="Shell Outlet Nozzle Diameter",Lower=10e-6);
848Dinozzle_Shell  as length               (Brief="Shell Inlet Nozzle Diameter",Lower=10e-6);
849Hinozzle_Shell  as length               (Brief="Height Under Shell Inlet Nozzle",Lower=10e-6);
850Honozzle_Shell  as length               (Brief="Height Under Shell Outlet Nozzle",Lower=10e-6);
851Lcf                     as length               (Brief="Bundle-to-Shell Clearance",Lower=10e-8);
852#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
853#       Tubes Geometrical Parameters
854#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
855Ntt                             as Integer              (Brief="Total Number of Tubes in Shell",Default=100,Lower=1);
856Pattern         as Integer              (Brief="Tube Layout Characteristic Angle",Lower=30);
857Ltube                   as length               (Brief="Effective Tube Length",Lower=0.1);
858pitch                   as length               (Brief="Tube Pitch",Lower=1e-8);
859Kwall                   as conductivity (Brief="Tube Wall Material Thermal Conductivity");
860Dotube                  as length               (Brief="Tube Outside Diameter",Lower=10e-6);
861Ditube                  as length               (Brief="Tube Inside Diameter",Lower=10e-6);
862Donozzle_Tube   as length               (Brief="Tube Outlet Nozzle Diameter",Lower=10e-6);
863Dinozzle_Tube   as length               (Brief="Tube Inlett Nozzle Diameter",Lower=10e-6);
864#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
865#       Baffles Geometrical Parameters
866#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
867Bc                      as Integer              (Brief="Baffle Cut",Default=25,Lower=25);
868Nb              as Real                 (Brief="Number of Baffles",Lower=1);
869Lcd                     as length               (Brief="Baffle-to-Shell Clearance",Lower=10e-8);
870Ltd                             as length               (Brief="Tube-to-Bafflehole Clearance",Lower=10e-8);
871#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
872VARIABLES
873
874Eft(Nshell)       as positive (Brief="Effectiveness",Default=0.05,Lower=1e-8);
875Unity(Nshell) as HeatExchangerDetailed_Basic;# "Shell in Series"
876
877CONNECTIONS
878
879Unity([1:Nshell-1]).Outlet.Hot  to Unity([2:Nshell]).Inlet.Hot;
880Unity([2:Nshell]).Outlet.Cold   to Unity([1:Nshell-1]).Inlet.Cold;
881
882EQUATIONS
883
884for i in [1:Nshell]
885       
886"Shell Side Cross Flow Area"
887        Unity(i).Shell.HeatTransfer.Sm = HE.CrossFlowArea(Unity(i).Baffles.Ls);
888
889"Ji Factor"
890        Unity(i).Shell.HeatTransfer.Ji =HE.JiFactor(Unity(i).Shell.HeatTransfer.Re);
891
892"Jc Factor"
893        Unity(i).Shell.HeatTransfer.Jc = HE.JcFactor();
894       
895"Jl Factor"
896        Unity(i).Shell.HeatTransfer.Jl = HE.JlFactor(Unity(i).Shell.HeatTransfer.Sm);
897
898"Jb Factor"
899        Unity(i).Shell.HeatTransfer.Jb = HE.JbFactor(Unity(i).Shell.HeatTransfer.Re,Unity(i).Baffles.Ls,Unity(i).Shell.HeatTransfer.Sm);
900
901"Jr Factor"     
902        Unity(i).Shell.HeatTransfer.Jr = HE.JrFactor(Unity(i).Shell.HeatTransfer.Re);
903       
904"Total J Factor"
905        Unity(i).Shell.HeatTransfer.Jtotal=Unity(i).Shell.HeatTransfer.Jc*Unity(i).Shell.HeatTransfer.Jl*Unity(i).Shell.HeatTransfer.Jb*Unity(i).Shell.HeatTransfer.Jr*Unity(i).Shell.HeatTransfer.Js;
906
907#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
908#  Fluid Alocation
909#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
910
911if side equal 1 then
912       
913"Shell Side Phi correction for viscosity"
914        Unity(i).Shell.HeatTransfer.Phi = HE.PhiCorrection(Unity(i).Properties.Hot.Average.Mu,Unity(i).Properties.Hot.Wall.Mu);
915       
916"Tube Side Phi correction for viscosity"
917        Unity(i).Tubes.HeatTransfer.Phi  = HE.PhiCorrection(Unity(i).Properties.Cold.Average.Mu,Unity(i).Properties.Cold.Wall.Mu);
918       
919        else
920       
921"Shell Side Phi correction for viscosity"
922        Unity(i).Shell.HeatTransfer.Phi = HE.PhiCorrection(Unity(i).Properties.Cold.Average.Mu,Unity(i).Properties.Cold.Wall.Mu);
923       
924"Tube Side Phi correction for viscosity"
925        Unity(i).Tubes.HeatTransfer.Phi  = HE.PhiCorrection(Unity(i).Properties.Hot.Average.Mu,Unity(i).Properties.Hot.Average.Mu);
926
927end
928
929
930if side equal 1
931       
932        then
933       
934"Hot Wall Temperature"
935        Unity(i).Properties.Hot.Wall.Twall  = (Unity(i).Properties.Hot.Average.T+Unity(i).Properties.Cold.Average.T)/2;
936
937"Cold Wall Temperature"
938        Unity(i).Properties.Cold.Wall.Twall  = (Unity(i).Properties.Hot.Average.T+Unity(i).Properties.Cold.Average.T)/2;
939
940"Tube Side Velocity"
941        Unity(i).Tubes.HeatTransfer.Vtube  = Unity(i).Properties.Cold.Inlet.Fw*Tpass/((Pi*Ditube*Ditube/4)*Unity(i).Properties.Cold.Average.rho*Ntt);
942
943"Tube Side Reynolds Number"
944        Unity(i).Tubes.HeatTransfer.Re=(Unity(i).Properties.Cold.Average.rho*Unity(i).Tubes.HeatTransfer.Vtube*Ditube)/Unity(i).Properties.Cold.Average.Mu;
945
946"Tube Side Prandtl Number"
947        Unity(i).Tubes.HeatTransfer.PR = ((Unity(i).Properties.Cold.Average.Cp/Unity(i).Properties.Cold.Average.Mw)*Unity(i).Properties.Cold.Average.Mu)/Unity(i).Properties.Cold.Average.K;
948
949"Tube Side Prandtl Number at Wall Temperature"
950        Unity(i).Tubes.HeatTransfer.PRw = ((Unity(i).Properties.Cold.Wall.Cp/Unity(i).Properties.Cold.Average.Mw)*Unity(i).Properties.Cold.Wall.Mu)/Unity(i).Properties.Cold.Wall.K;
951
952"Tube Side Film Coefficient"
953        Unity(i).Tubes.HeatTransfer.htube= HE.TubeFilmCoeff(Unity(i).Tubes.HeatTransfer.Re,Unity(i).Tubes.HeatTransfer.PR,Unity(i).Properties.Cold.Average.K)*Unity(i).Tubes.HeatTransfer.Phi;
954       
955"Shell Side Prandtl Number"
956        Unity(i).Shell.HeatTransfer.PR=((Unity(i).Properties.Hot.Average.Cp/Unity(i).Properties.Hot.Average.Mw)*Unity(i).Properties.Hot.Average.Mu)/Unity(i).Properties.Hot.Average.K; 
957
958"Shell Side Prandtl Number at Wall Temperature"
959        Unity(i).Shell.HeatTransfer.PRw = ((Unity(i).Properties.Hot.Wall.Cp/Unity(i).Properties.Hot.Average.Mw)*Unity(i).Properties.Hot.Wall.Mu)/Unity(i).Properties.Hot.Wall.K;       
960
961"Tube Side Pressure Drop"
962        Unity(i).Tubes.PressureDrop.PdTube              = HE.DeltaPtube(Unity(i).Tubes.HeatTransfer.Re,Unity(i).Properties.Cold.Average.rho,Unity(i).Tubes.HeatTransfer.Vtube,Unity(i).Tubes.HeatTransfer.Phi);
963
964"Pressure Drop Tube Side Inlet Nozzle"
965        Unity(i).Tubes.PressureDrop.Pdnozzle_in = HE.DeltaPtubeNozzlein(Unity(i).Properties.Cold.Inlet.rho,Unity(i).Properties.Cold.Inlet.Fw);
966
967"Velocity Tube Side Inlet Nozzle"
968        Unity(i).Tubes.PressureDrop.Vnozzle_in  = HE.TubeVelocityNozzlein(Unity(i).Properties.Cold.Inlet.rho,Unity(i).Properties.Cold.Inlet.Fw);
969
970"Pressure Drop Tube Side Outlet Nozzle"
971        Unity(i).Tubes.PressureDrop.Pdnozzle_out        = HE.DeltaPtubeNozzleout(Unity(i).Properties.Cold.Outlet.rho,Unity(i).Properties.Cold.Inlet.Fw);
972
973"Velocity Tube Side Outlet Nozzle"
974        Unity(i).Tubes.PressureDrop.Vnozzle_out = HE.TubeVelocityNozzleout(Unity(i).Properties.Cold.Outlet.rho,Unity(i).Properties.Cold.Inlet.Fw);
975
976"Shell Pressure Drop Inlet Nozzle"
977        Unity(i).Shell.PressureDrop.Pdnozzle_in = HE.DeltaPshellNozzleIn(Unity(i).Properties.Hot.Inlet.rho,Unity(i).Properties.Hot.Inlet.Fw);
978
979"Velocity Shell Side Inlet Nozzle"
980        Unity(i).Shell.PressureDrop.Vnozzle_in  = HE.ShellVelocityNozzleIn(Unity(i).Properties.Hot.Inlet.rho,Unity(i).Properties.Hot.Inlet.Fw);
981
982"Shell Pressure Drop Outlet Nozzle"
983        Unity(i).Shell.PressureDrop.Pdnozzle_out        =HE.DeltaPshellNozzleOut(Unity(i).Properties.Hot.Outlet.rho,Unity(i).Properties.Hot.Inlet.Fw);
984
985"Velocity Shell Side Outlet Nozzle"
986        Unity(i).Shell.PressureDrop.Vnozzle_out =HE.ShellVelocityNozzleOut(Unity(i).Properties.Hot.Outlet.rho,Unity(i).Properties.Hot.Inlet.Fw);
987
988"Pressure Drop Hot Stream"
989        Unity(i).Outlet.Hot.P  = Unity(i).Inlet.Hot.P- Unity(i).Shell.PressureDrop.Pdtotal;     
990       
991"Pressure Drop Cold Stream"
992        Unity(i).Outlet.Cold.P  = Unity(i).Inlet.Cold.P - Unity(i).Tubes.PressureDrop.Pdtotal;
993
994
995        else
996
997"Hot Wall Temperature"
998        Unity(i).Properties.Hot.Wall.Twall  = (Unity(i).Properties.Hot.Average.T+Unity(i).Properties.Cold.Average.T)/2;
999
1000"Cold Wall Temperature"
1001        Unity(i).Properties.Cold.Wall.Twall  = (Unity(i).Properties.Hot.Average.T+Unity(i).Properties.Cold.Average.T)/2;
1002
1003"Tube Side Velocity"
1004        Unity(i).Tubes.HeatTransfer.Vtube  = Unity(i).Properties.Hot.Inlet.Fw*Tpass/((Pi*Ditube*Ditube/4)*Unity(i).Properties.Hot.Average.rho*Ntt);
1005
1006"Tube Side Reynolds Number"
1007        Unity(i).Tubes.HeatTransfer.Re          =       (Unity(i).Properties.Hot.Average.rho*Unity(i).Tubes.HeatTransfer.Vtube*Ditube)/Unity(i).Properties.Hot.Average.Mu;
1008       
1009"Tube Side Prandtl Number"
1010        Unity(i).Tubes.HeatTransfer.PR          =       ((Unity(i).Properties.Hot.Average.Cp/Unity(i).Properties.Hot.Average.Mw)*Unity(i).Properties.Hot.Average.Mu)/Unity(i).Properties.Hot.Average.K;
1011       
1012"Tube Side Prandtl Number at Wall"
1013        Unity(i).Tubes.HeatTransfer.PRw     =  ((Unity(i).Properties.Hot.Wall.Cp/Unity(i).Properties.Hot.Average.Mw)*Unity(i).Properties.Hot.Wall.Mu)/Unity(i).Properties.Hot.Wall.K;
1014
1015"Tube Side Film Coefficient"
1016        Unity(i).Tubes.HeatTransfer.htube       =       HE.TubeFilmCoeff(Unity(i).Tubes.HeatTransfer.Re,Unity(i).Tubes.HeatTransfer.PR,Unity(i).Properties.Hot.Average.K)*Unity(i).Tubes.HeatTransfer.Phi;
1017       
1018"Shell Side Prandtl Number"
1019        Unity(i).Shell.HeatTransfer.PR          =       ((Unity(i).Properties.Cold.Average.Cp/Unity(i).Properties.Cold.Average.Mw)*Unity(i).Properties.Cold.Average.Mu)/Unity(i).Properties.Cold.Average.K;
1020       
1021"Shell Side Prandtl Number at Wall"
1022        Unity(i).Shell.HeatTransfer.PRw     =   ((Unity(i).Properties.Cold.Wall.Cp/Unity(i).Properties.Cold.Average.Mw)*Unity(i).Properties.Cold.Wall.Mu)/Unity(i).Properties.Cold.Wall.K;
1023
1024"Tube Side Pressure Drop"
1025        Unity(i).Tubes.PressureDrop.PdTube      =       HE.DeltaPtube(Unity(i).Tubes.HeatTransfer.Re,Unity(i).Properties.Hot.Average.rho,Unity(i).Tubes.HeatTransfer.Vtube,Unity(i).Tubes.HeatTransfer.Phi);
1026
1027"Pressure Drop Tube Side Inlet Nozzle"
1028        Unity(i).Tubes.PressureDrop.Pdnozzle_in = HE.DeltaPtubeNozzlein(Unity(i).Properties.Hot.Inlet.rho,Unity(i).Properties.Hot.Inlet.Fw);
1029
1030"Velocity Tube Side Inlet Nozzle"
1031        Unity(i).Tubes.PressureDrop.Vnozzle_in  = HE.TubeVelocityNozzlein(Unity(i).Properties.Hot.Inlet.rho,Unity(i).Properties.Hot.Inlet.Fw);
1032
1033"Pressure Drop Tube Side Outlet Nozzle"
1034        Unity(i).Tubes.PressureDrop.Pdnozzle_out        = HE.DeltaPtubeNozzleout(Unity(i).Properties.Hot.Outlet.rho,Unity(i).Properties.Hot.Inlet.Fw); 
1035
1036"Velocity Tube Side Outlet Nozzle"
1037        Unity(i).Tubes.PressureDrop.Vnozzle_out = HE.TubeVelocityNozzleout(Unity(i).Properties.Hot.Outlet.rho,Unity(i).Properties.Hot.Inlet.Fw);       
1038
1039"Shell Pressure Drop Inlet Nozzle"
1040        Unity(i).Shell.PressureDrop.Pdnozzle_in = HE.DeltaPshellNozzleIn(Unity(i).Properties.Cold.Inlet.rho,Unity(i).Properties.Cold.Inlet.Fw);
1041
1042"Velocity Shell Side Inlet Nozzle"
1043        Unity(i).Shell.PressureDrop.Vnozzle_in  = HE.ShellVelocityNozzleIn(Unity(i).Properties.Cold.Inlet.rho,Unity(i).Properties.Cold.Inlet.Fw);
1044
1045"Shell Pressure Drop Outlet Nozzle"
1046        Unity(i).Shell.PressureDrop.Pdnozzle_out = HE.DeltaPshellNozzleOut(Unity(i).Properties.Cold.Outlet.rho,Unity(i).Properties.Cold.Inlet.Fw);
1047
1048"Velocity Shell Side Outlet Nozzle"
1049        Unity(i).Shell.PressureDrop.Vnozzle_out  = HE.ShellVelocityNozzleOut(Unity(i).Properties.Cold.Outlet.rho,Unity(i).Properties.Cold.Inlet.Fw);
1050
1051"Pressure Drop Hot Stream"
1052        Unity(i).Outlet.Hot.P  = Unity(i).Inlet.Hot.P- Unity(i).Tubes.PressureDrop.Pdtotal;     
1053       
1054"Pressure Drop Cold Stream"
1055        Unity(i).Outlet.Cold.P  = Unity(i).Inlet.Cold.P - Unity(i).Shell.PressureDrop.Pdtotal;
1056
1057
1058end
1059
1060if side equal 1
1061       
1062        then
1063
1064"Shell Side inlet Nozzle rho-V^2"
1065        Unity(i).Shell.PressureDrop.RVsquare_in = Unity(i).Properties.Hot.Inlet.rho*(Unity(i).Shell.PressureDrop.Vnozzle_in)^2;
1066
1067"Shell Side Outlet Nozzle rho-V^2"
1068        Unity(i).Shell.PressureDrop.RVsquare_out = Unity(i).Properties.Hot.Outlet.rho*(Unity(i).Shell.PressureDrop.Vnozzle_out)^2;
1069       
1070        else
1071       
1072"Shell Side inlet Nozzle rho-V^2"
1073        Unity(i).Shell.PressureDrop.RVsquare_in = Unity(i).Properties.Cold.Inlet.rho*(Unity(i).Shell.PressureDrop.Vnozzle_in)^2;
1074
1075"Shell Side Outlet Nozzle rho-V^2"
1076        Unity(i).Shell.PressureDrop.RVsquare_out = Unity(i).Properties.Cold.Outlet.rho*(Unity(i).Shell.PressureDrop.Vnozzle_out)^2;
1077
1078end
1079
1080"Tube Resistance"       
1081        Unity(i).Resistances.Rtube*(Unity(i).Tubes.HeatTransfer.htube*Ditube) = Dotube;
1082       
1083"Wall Resistance"
1084        Unity(i).Resistances.Rwall=Dotube*ln(Dotube/Ditube)/(2*Kwall);
1085       
1086"Shell Resistance"
1087        Unity(i).Resistances.Rshell*(Unity(i).Shell.HeatTransfer.hshell)=1;
1088       
1089"Overall Heat Transfer Coefficient Clean"
1090        Unity(i).Details.Uc=1/(Dotube/(Unity(i).Tubes.HeatTransfer.htube*Ditube)+(Dotube*ln(Dotube/Ditube)/(2*Kwall))+(1/(Unity(i).Shell.HeatTransfer.hshell)));
1091
1092"Overall Heat Transfer Coefficient Dirty"
1093        (1/Unity(i).Details.Ud)=(1/Unity(i).Details.Uc)+Unity(i).Resistances.Rfo+Unity(i).Resistances.Rfi*(Dotube/Ditube);
1094
1095"Exchange Surface Area"
1096        Unity(i).Details.A=Pi*Dotube*Ntt*Ltube;
1097
1098"Baffles Spacing"
1099        Ltube = Unity(i).Baffles.Lsi+Unity(i).Baffles.Lso+Unity(i).Baffles.Ls*(Nb-1);
1100
1101"TEMA E Shell Effectiveness"
1102        Eft(i) = 2*(1+Unity(i).Details.Cr+sqrt(1+Unity(i).Details.Cr^2)*((1+exp(-Unity(i).Details.NTU*sqrt(1+Unity(i).Details.Cr^2)))/(1-exp(-Unity(i).Details.NTU*sqrt(1+Unity(i).Details.Cr^2)))) )^-1;
1103
1104"Js Factor"     
1105        Unity(i).Shell.HeatTransfer.Js = HE.JsFactor(Unity(i).Shell.HeatTransfer.Re,Unity(i).Baffles.Lsi,Unity(i).Baffles.Lso,Unity(i).Baffles.Ls);
1106
1107if side equal 1
1108       
1109        then
1110       
1111"Shell Side Reynolds Number"
1112        Unity(i).Shell.HeatTransfer.Re=(Dotube*Unity(i).Properties.Hot.Inlet.Fw/Unity(i).Shell.HeatTransfer.Sm)/Unity(i).Properties.Hot.Average.Mu;
1113
1114"Shell Heat Transfer Coefficient"
1115        Unity(i).Shell.HeatTransfer.hshell=Unity(i).Shell.HeatTransfer.Ji*(Unity(i).Properties.Hot.Average.Cp/Unity(i).Properties.Hot.Average.Mw)*(Unity(i).Properties.Hot.Inlet.Fw/Unity(i).Shell.HeatTransfer.Sm)*(Unity(i).Shell.HeatTransfer.PR^(-2/3))*Unity(i).Shell.HeatTransfer.Jtotal*Unity(i).Shell.HeatTransfer.Phi;
1116
1117"Shell Pressure Drop Cross Flow"
1118        Unity(i).Shell.PressureDrop.PdCross             = HE.DeltaPcross(Unity(i).Shell.HeatTransfer.Re,Unity(i).Baffles.Ls,Unity(i).Baffles.Lso,Unity(i).Baffles.Lsi,Unity(i).Properties.Hot.Inlet.Fw,Unity(i).Shell.HeatTransfer.Phi,Unity(i).Properties.Hot.Average.rho);
1119
1120"Shell Pressure Baffle Window"
1121        Unity(i).Shell.PressureDrop.Pdwindow            = HE.DeltaPwindow(Unity(i).Properties.Hot.Inlet.Fw,Unity(i).Shell.HeatTransfer.Sm,Unity(i).Properties.Hot.Average.rho,Unity(i).Properties.Hot.Average.Mu,Unity(i).Baffles.Ls);
1122
1123"Shell Pressure End Zones"
1124        Unity(i).Shell.PressureDrop.PdEndZones  = HE.DeltaPendZones(Unity(i).Shell.HeatTransfer.Re,Unity(i).Baffles.Ls,Unity(i).Baffles.Lso,Unity(i).Baffles.Lsi,Unity(i).Properties.Hot.Inlet.Fw,Unity(i).Shell.HeatTransfer.Phi,Unity(i).Properties.Hot.Average.rho);
1125
1126
1127        else
1128
1129"Shell Side Reynolds Number"
1130        Unity(i).Shell.HeatTransfer.Re=(Dotube*Unity(i).Properties.Cold.Inlet.Fw/Unity(i).Shell.HeatTransfer.Sm)/Unity(i).Properties.Cold.Average.Mu;
1131
1132"Shell Heat Transfer Coefficient"
1133        Unity(i).Shell.HeatTransfer.hshell=Unity(i).Shell.HeatTransfer.Ji*(Unity(i).Properties.Cold.Average.Cp/Unity(i).Properties.Cold.Average.Mw)*(Unity(i).Properties.Cold.Inlet.Fw/Unity(i).Shell.HeatTransfer.Sm)*(Unity(i).Shell.HeatTransfer.PR^(-2/3))*Unity(i).Shell.HeatTransfer.Jtotal*Unity(i).Shell.HeatTransfer.Phi;
1134
1135"Shell Pressure Drop Cross Flow"
1136        Unity(i).Shell.PressureDrop.PdCross      = HE.DeltaPcross(Unity(i).Shell.HeatTransfer.Re,Unity(i).Baffles.Ls,Unity(i).Baffles.Lso,Unity(i).Baffles.Lsi,Unity(i).Properties.Cold.Inlet.Fw,Unity(i).Shell.HeatTransfer.Phi,Unity(i).Properties.Cold.Average.rho);
1137
1138"Shell Pressure Baffle Window"
1139        Unity(i).Shell.PressureDrop.Pdwindow     = HE.DeltaPwindow(Unity(i).Properties.Cold.Inlet.Fw,Unity(i).Shell.HeatTransfer.Sm,Unity(i).Properties.Cold.Average.rho,Unity(i).Properties.Cold.Average.Mu,Unity(i).Baffles.Ls);
1140
1141"Shell Pressure End Zones"
1142        Unity(i).Shell.PressureDrop.PdEndZones  = HE.DeltaPendZones(Unity(i).Shell.HeatTransfer.Re,Unity(i).Baffles.Ls,Unity(i).Baffles.Lso,Unity(i).Baffles.Lsi,Unity(i).Properties.Cold.Inlet.Fw,Unity(i).Shell.HeatTransfer.Phi,Unity(i).Properties.Cold.Average.rho);
1143
1144
1145end
1146
1147"Duty"
1148        Unity(i).Details.Q      = Eft(i)*Unity(i).Details.Cmin*(Unity(i).Inlet.Hot.T-Unity(i).Inlet.Cold.T);
1149
1150end
1151
1152        SET
1153#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
1154#  Set Parameters for heatex Calculation
1155#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
1156        Pi                              = 3.14159265;
1157        HE.Tpass                        = Tpass;
1158        HE.Nss                          = Nss;
1159        HE.Ntt                          = Ntt;
1160        HE.Pattern                      = Pattern;
1161        HE.Bc                           = Bc;
1162        HE.Donozzle_Shell       = Donozzle_Shell;
1163        HE.Dinozzle_Shell       = Dinozzle_Shell;
1164        HE.Honozzle_Shell       = Honozzle_Shell;
1165        HE.Hinozzle_Shell       = Hinozzle_Shell;
1166        HE.Donozzle_Tube        = Donozzle_Tube;
1167        HE.Dinozzle_Tube        = Dinozzle_Tube;
1168        HE.Nb                       = Nb;
1169        HE.Dishell                      = Dishell;             
1170        HE.Lcf                          = Lcf; 
1171        HE.Ltube                        = Ltube;       
1172        HE.pitch                        = pitch;               
1173        HE.Dotube                       = Dotube;       
1174        HE.Ditube                       = Ditube;       
1175        HE.Lcd                          = Lcd;
1176        HE.Ltd                          = Ltd; 
1177        side                            = HE.FluidAlocation();
1178
1179end
1180
1181Model E_Shell_LMTD_Det                  as Heatex_Detailed_LMTD
1182#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
1183#       Shell and Tubes Heat Exchanger with 1 shell pass - LMTD Method
1184#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
1185EQUATIONS
1186
1187"LMTD Correction Factor"
1188        Fc = HE.EshellCorrectionFactor(Inlet.Hot.T,Outlet.Hot.T,Inlet.Cold.T,Outlet.Cold.T);
1189
1190"Temperature Difference at Inlet"
1191        DT0 = Inlet.Hot.T - Outlet.Cold.T;
1192
1193"Temperature Difference at Outlet"
1194        DTL = Outlet.Hot.T - Inlet.Cold.T;
1195       
1196"Js Factor"     
1197        Shell.HeatTransfer.Js = HE.JsFactor(Shell.HeatTransfer.Re,Baffles.Lsi,Baffles.Lso,Baffles.Ls);
1198
1199if side equal 1
1200       
1201        then
1202       
1203"Shell Side Reynolds Number"
1204        Shell.HeatTransfer.Re=(Dotube*Properties.Hot.Inlet.Fw/Shell.HeatTransfer.Sm)/Properties.Hot.Average.Mu;
1205
1206"Shell Heat Transfer Coefficient"
1207        Shell.HeatTransfer.hshell=Shell.HeatTransfer.Ji*(Properties.Hot.Average.Cp/Properties.Hot.Average.Mw)*(Properties.Hot.Inlet.Fw/Shell.HeatTransfer.Sm)*(Shell.HeatTransfer.PR^(-2/3))*Shell.HeatTransfer.Jtotal*Shell.HeatTransfer.Phi;
1208
1209"Shell Pressure Drop Cross Flow"
1210        Shell.PressureDrop.PdCross              = HE.DeltaPcross(Shell.HeatTransfer.Re,Baffles.Ls,Baffles.Lso,Baffles.Lsi,Properties.Hot.Inlet.Fw,Shell.HeatTransfer.Phi,Properties.Hot.Average.rho);
1211
1212"Shell Pressure Baffle Window"
1213        Shell.PressureDrop.Pdwindow             = HE.DeltaPwindow(Properties.Hot.Inlet.Fw,Shell.HeatTransfer.Sm,Properties.Hot.Average.rho,Properties.Hot.Average.Mu,Baffles.Ls);
1214
1215"Shell Pressure End Zones"
1216        Shell.PressureDrop.PdEndZones  = HE.DeltaPendZones(Shell.HeatTransfer.Re,Baffles.Ls,Baffles.Lso,Baffles.Lsi,Properties.Hot.Inlet.Fw,Shell.HeatTransfer.Phi,Properties.Hot.Average.rho);
1217
1218
1219        else
1220
1221"Shell Side Reynolds Number"
1222        Shell.HeatTransfer.Re=(Dotube*Properties.Cold.Inlet.Fw/Shell.HeatTransfer.Sm)/Properties.Cold.Average.Mu;
1223
1224"Shell Heat Transfer Coefficient"
1225        Shell.HeatTransfer.hshell=Shell.HeatTransfer.Ji*(Properties.Cold.Average.Cp/Properties.Cold.Average.Mw)*(Properties.Cold.Inlet.Fw/Shell.HeatTransfer.Sm)*(Shell.HeatTransfer.PR^(-2/3))*Shell.HeatTransfer.Jtotal*Shell.HeatTransfer.Phi;
1226
1227"Shell Pressure Drop Cross Flow"
1228        Shell.PressureDrop.PdCross      = HE.DeltaPcross(Shell.HeatTransfer.Re,Baffles.Ls,Baffles.Lso,Baffles.Lsi,Properties.Cold.Inlet.Fw,Shell.HeatTransfer.Phi,Properties.Cold.Average.rho);
1229
1230"Shell Pressure Baffle Window"
1231        Shell.PressureDrop.Pdwindow     = HE.DeltaPwindow(Properties.Cold.Inlet.Fw,Shell.HeatTransfer.Sm,Properties.Cold.Average.rho,Properties.Cold.Average.Mu,Baffles.Ls);
1232
1233"Shell Pressure End Zones"
1234        Shell.PressureDrop.PdEndZones  = HE.DeltaPendZones(Shell.HeatTransfer.Re,Baffles.Ls,Baffles.Lso,Baffles.Lsi,Properties.Cold.Inlet.Fw,Shell.HeatTransfer.Phi,Properties.Cold.Average.rho);
1235
1236
1237end
1238
1239
1240
1241end
1242
1243Model F_Shell_LMTD_Det                  as Heatex_Detailed_LMTD
1244#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
1245#       Shell and Tubes Heat Exchanger with 2 shell pass - LMTD Method
1246#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
1247
1248EQUATIONS
1249
1250"LMTD Correction Factor"
1251        Fc = HE.FshellCorrectionFactor(Inlet.Hot.T,Outlet.Hot.T,Inlet.Cold.T,Outlet.Cold.T);
1252
1253"Temperature Difference at Inlet"
1254        DT0 = Inlet.Hot.T - Outlet.Cold.T;
1255
1256"Temperature Difference at Outlet"
1257        DTL = Outlet.Hot.T - Inlet.Cold.T;
1258       
1259"Js Factor"     
1260        Shell.HeatTransfer.Js = HE.JsFactor(Shell.HeatTransfer.Re,Baffles.Lsi,Baffles.Lso,Baffles.Ls);
1261
1262if side equal 1
1263       
1264        then
1265       
1266"Shell Side Reynolds Number"
1267        Shell.HeatTransfer.Re=(Dotube*Properties.Hot.Inlet.Fw/Shell.HeatTransfer.Sm)/Properties.Hot.Average.Mu;
1268
1269"Shell Heat Transfer Coefficient"
1270        Shell.HeatTransfer.hshell=Shell.HeatTransfer.Ji*(Properties.Hot.Average.Cp/Properties.Hot.Average.Mw)*(Properties.Hot.Inlet.Fw/Shell.HeatTransfer.Sm)*(Shell.HeatTransfer.PR^(-2/3))*Shell.HeatTransfer.Jtotal*Shell.HeatTransfer.Phi;
1271
1272"Shell Pressure Drop Cross Flow"
1273        Shell.PressureDrop.PdCross              = HE.DeltaPcross(Shell.HeatTransfer.Re,Baffles.Ls,Baffles.Lso,Baffles.Lsi,Properties.Hot.Inlet.Fw,Shell.HeatTransfer.Phi,Properties.Hot.Average.rho);
1274
1275"Shell Pressure Baffle Window"
1276        Shell.PressureDrop.Pdwindow             = HE.DeltaPwindow(Properties.Hot.Inlet.Fw,Shell.HeatTransfer.Sm,Properties.Hot.Average.rho,Properties.Hot.Average.Mu,Baffles.Ls);
1277
1278"Shell Pressure End Zones"
1279        Shell.PressureDrop.PdEndZones  = HE.DeltaPendZones(Shell.HeatTransfer.Re,Baffles.Ls,Baffles.Lso,Baffles.Lsi,Properties.Hot.Inlet.Fw,Shell.HeatTransfer.Phi,Properties.Hot.Average.rho);
1280
1281
1282        else
1283
1284"Shell Side Reynolds Number"
1285        Shell.HeatTransfer.Re=(Dotube*Properties.Cold.Inlet.Fw/Shell.HeatTransfer.Sm)/Properties.Cold.Average.Mu;
1286
1287"Shell Heat Transfer Coefficient"
1288        Shell.HeatTransfer.hshell=Shell.HeatTransfer.Ji*(Properties.Cold.Average.Cp/Properties.Cold.Average.Mw)*(Properties.Cold.Inlet.Fw/Shell.HeatTransfer.Sm)*(Shell.HeatTransfer.PR^(-2/3))*Shell.HeatTransfer.Jtotal*Shell.HeatTransfer.Phi;
1289
1290"Shell Pressure Drop Cross Flow"
1291        Shell.PressureDrop.PdCross      = HE.DeltaPcross(Shell.HeatTransfer.Re,Baffles.Ls,Baffles.Lso,Baffles.Lsi,Properties.Cold.Inlet.Fw,Shell.HeatTransfer.Phi,Properties.Cold.Average.rho);
1292
1293"Shell Pressure Baffle Window"
1294        Shell.PressureDrop.Pdwindow     = HE.DeltaPwindow(Properties.Cold.Inlet.Fw,Shell.HeatTransfer.Sm,Properties.Cold.Average.rho,Properties.Cold.Average.Mu,Baffles.Ls);
1295
1296"Shell Pressure End Zones"
1297        Shell.PressureDrop.PdEndZones  = HE.DeltaPendZones(Shell.HeatTransfer.Re,Baffles.Ls,Baffles.Lso,Baffles.Lsi,Properties.Cold.Inlet.Fw,Shell.HeatTransfer.Phi,Properties.Cold.Average.rho);
1298
1299
1300end
1301
1302end
1303
1304Model Multipass_LMTD_Det         
1305#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
1306#       Shell and Tubes Heat Exchanger In Series with 1 shell pass - LMTD Method
1307#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
1308PARAMETERS
1309
1310Nshell          as Integer              (Brief="N Shell in Series",Default=2);
1311        HE              as CalcObject   (Brief="STHE Calculations",File="heatex");
1312ext PP          as CalcObject   (Brief="External Physical Properties");
1313side        as Integer          (Brief="Fluid Alocation",Lower=0,Upper=1);
1314Pi                      as constant     (Brief="Pi Number",Default=3.14159265);
1315
1316#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
1317#       Shell Geometrical Parameters
1318#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
1319Tpass                   as Integer              (Brief="Number of Tube Passes",Lower=1);
1320Nss                             as Integer              (Brief="Number of Sealing Strips pairs",Lower=1);
1321Dishell                 as length               (Brief="Inside Shell Diameter",Lower=10e-6);
1322Donozzle_Shell  as length               (Brief="Shell Outlet Nozzle Diameter",Lower=10e-6);
1323Dinozzle_Shell  as length               (Brief="Shell Inlet Nozzle Diameter",Lower=10e-6);
1324Hinozzle_Shell  as length               (Brief="Height Under Shell Inlet Nozzle",Lower=10e-6);
1325Honozzle_Shell  as length               (Brief="Height Under Shell Outlet Nozzle",Lower=10e-6);
1326Lcf                     as length               (Brief="Bundle-to-Shell Clearance",Lower=10e-8);
1327#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
1328#       Tubes Geometrical Parameters
1329#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
1330Ntt                             as Integer              (Brief="Total Number of Tubes in Shell",Default=100,Lower=1);
1331Pattern         as Integer              (Brief="Tube Layout Characteristic Angle",Lower=30);
1332Ltube                   as length               (Brief="Effective Tube Length",Lower=0.1);
1333pitch                   as length               (Brief="Tube Pitch",Lower=1e-8);
1334Kwall                   as conductivity (Brief="Tube Wall Material Thermal Conductivity");
1335Dotube                  as length               (Brief="Tube Outside Diameter",Lower=10e-6);
1336Ditube                  as length               (Brief="Tube Inside Diameter",Lower=10e-6);
1337Donozzle_Tube   as length               (Brief="Tube Outlet Nozzle Diameter",Lower=10e-6);
1338Dinozzle_Tube   as length               (Brief="Tube Inlett Nozzle Diameter",Lower=10e-6);
1339#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
1340#       Baffles Geometrical Parameters
1341#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
1342Bc                      as Integer              (Brief="Baffle Cut",Default=25,Lower=25);
1343Nb              as Real                 (Brief="Number of Baffles",Lower=1);
1344Lcd                     as length               (Brief="Baffle-to-Shell Clearance",Lower=10e-8);
1345Ltd                             as length               (Brief="Tube-to-Bafflehole Clearance",Lower=10e-8);
1346#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
1347VARIABLES
1348
1349LMTD(Nshell)            as temp_delta   (Brief="Logarithmic Mean Temperature Difference",Lower=5);
1350Fc(Nshell)                      as positive             (Brief="LMTD Correction Factor",Lower=0.5);
1351MTD(Nshell)                     as temp_delta   (Brief="Mean Temperature Difference",Lower=5);
1352Unity(Nshell)           as HeatExchangerDetailed_Basic;
1353
1354CONNECTIONS
1355
1356Unity([1:Nshell-1]).Outlet.Hot  to Unity([2:Nshell]).Inlet.Hot;
1357Unity([2:Nshell]).Outlet.Cold to Unity([1:Nshell-1]).Inlet.Cold;
1358
1359EQUATIONS
1360
1361for i in [1:Nshell]
1362       
1363"Shell Side Cross Flow Area"
1364        Unity(i).Shell.HeatTransfer.Sm = HE.CrossFlowArea(Unity(i).Baffles.Ls);
1365
1366"Ji Factor"
1367        Unity(i).Shell.HeatTransfer.Ji = HE.JiFactor(Unity(i).Shell.HeatTransfer.Re);
1368
1369"Jc Factor"
1370        Unity(i).Shell.HeatTransfer.Jc = HE.JcFactor();
1371       
1372"Jl Factor"
1373        Unity(i).Shell.HeatTransfer.Jl = HE.JlFactor(Unity(i).Shell.HeatTransfer.Sm);
1374
1375"Jb Factor"
1376        Unity(i).Shell.HeatTransfer.Jb = HE.JbFactor(Unity(i).Shell.HeatTransfer.Re,Unity(i).Baffles.Ls,Unity(i).Shell.HeatTransfer.Sm);
1377
1378"Jr Factor"     
1379        Unity(i).Shell.HeatTransfer.Jr = HE.JrFactor(Unity(i).Shell.HeatTransfer.Re);
1380       
1381"Total J Factor"
1382        Unity(i).Shell.HeatTransfer.Jtotal=Unity(i).Shell.HeatTransfer.Jc*Unity(i).Shell.HeatTransfer.Jl*Unity(i).Shell.HeatTransfer.Jb*Unity(i).Shell.HeatTransfer.Jr*Unity(i).Shell.HeatTransfer.Js;
1383
1384#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
1385#  Fluid Alocation
1386#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
1387
1388if side equal 1 then
1389       
1390       
1391"Shell Side Phi correction for viscosity"
1392        Unity(i).Shell.HeatTransfer.Phi = HE.PhiCorrection(Unity(i).Properties.Hot.Average.Mu,Unity(i).Properties.Hot.Wall.Mu);
1393       
1394"Tube Side Phi correction for viscosity"
1395        Unity(i).Tubes.HeatTransfer.Phi  = HE.PhiCorrection(Unity(i).Properties.Cold.Average.Mu,Unity(i).Properties.Cold.Wall.Mu);
1396       
1397        else
1398       
1399       
1400"Shell Side Phi correction for viscosity"
1401        Unity(i).Shell.HeatTransfer.Phi = HE.PhiCorrection(Unity(i).Properties.Cold.Average.Mu,Unity(i).Properties.Cold.Wall.Mu);
1402       
1403"Tube Side Phi correction for viscosity"
1404        Unity(i).Tubes.HeatTransfer.Phi  = HE.PhiCorrection(Unity(i).Properties.Hot.Average.Mu,Unity(i).Properties.Hot.Wall.Mu);
1405
1406end
1407
1408
1409if side equal 1
1410       
1411        then
1412       
1413"Cold Wall Temperature"
1414        Unity(i).Properties.Cold.Wall.Twall  = (Unity(i).Properties.Hot.Average.T+Unity(i).Properties.Cold.Average.T)/2;
1415
1416"Hot Wall Temperature" 
1417        Unity(i).Properties.Hot.Wall.Twall  = (Unity(i).Properties.Hot.Average.T+Unity(i).Properties.Cold.Average.T)/2;
1418
1419"Tube Side Velocity"
1420        Unity(i).Tubes.HeatTransfer.Vtube  = Unity(i).Properties.Cold.Inlet.Fw*Tpass/((Pi*Ditube*Ditube/4)*Unity(i).Properties.Cold.Average.rho*Ntt);
1421
1422"Tube Side Reynolds Number"
1423        Unity(i).Tubes.HeatTransfer.Re=(Unity(i).Properties.Cold.Average.rho*Unity(i).Tubes.HeatTransfer.Vtube*Ditube)/Unity(i).Properties.Cold.Average.Mu;
1424
1425"Tube Side Prandtl Number"
1426        Unity(i).Tubes.HeatTransfer.PR = ((Unity(i).Properties.Cold.Average.Cp/Unity(i).Properties.Cold.Average.Mw)*Unity(i).Properties.Cold.Average.Mu)/Unity(i).Properties.Cold.Average.K;
1427
1428"Tube Side Prandtl Number at Wall"
1429        Unity(i).Tubes.HeatTransfer.PRw = ((Unity(i).Properties.Cold.Wall.Cp/Unity(i).Properties.Cold.Average.Mw)*Unity(i).Properties.Cold.Wall.Mu)/Unity(i).Properties.Cold.Wall.K;
1430
1431"Tube Side Film Coefficient"
1432        Unity(i).Tubes.HeatTransfer.htube= HE.TubeFilmCoeff(Unity(i).Tubes.HeatTransfer.Re,Unity(i).Tubes.HeatTransfer.PR,Unity(i).Properties.Cold.Average.K)*Unity(i).Tubes.HeatTransfer.Phi;
1433       
1434"Shell Side Prandtl Number"
1435        Unity(i).Shell.HeatTransfer.PR=((Unity(i).Properties.Hot.Average.Cp/Unity(i).Properties.Hot.Average.Mw)*Unity(i).Properties.Hot.Average.Mu)/Unity(i).Properties.Hot.Average.K; 
1436
1437"Shell Side Prandtl Number at Wall"
1438        Unity(i).Shell.HeatTransfer.PRw=((Unity(i).Properties.Hot.Wall.Cp/Unity(i).Properties.Hot.Average.Mw)*Unity(i).Properties.Hot.Wall.Mu)/Unity(i).Properties.Hot.Wall.K; 
1439
1440"Tube Side Pressure Drop"
1441        Unity(i).Tubes.PressureDrop.PdTube              = HE.DeltaPtube(Unity(i).Tubes.HeatTransfer.Re,Unity(i).Properties.Cold.Average.rho,Unity(i).Tubes.HeatTransfer.Vtube,Unity(i).Tubes.HeatTransfer.Phi);
1442
1443"Pressure Drop Tube Side Inlet Nozzle"
1444        Unity(i).Tubes.PressureDrop.Pdnozzle_in = HE.DeltaPtubeNozzlein(Unity(i).Properties.Cold.Inlet.rho,Unity(i).Properties.Cold.Inlet.Fw);
1445
1446"Velocity Tube Side Inlet Nozzle"
1447        Unity(i).Tubes.PressureDrop.Vnozzle_in  = HE.TubeVelocityNozzlein(Unity(i).Properties.Cold.Inlet.rho,Unity(i).Properties.Cold.Inlet.Fw);
1448
1449"Pressure Drop Tube Side Outlet Nozzle"
1450        Unity(i).Tubes.PressureDrop.Pdnozzle_out        = HE.DeltaPtubeNozzleout(Unity(i).Properties.Cold.Outlet.rho,Unity(i).Properties.Cold.Inlet.Fw);
1451
1452"Velocity Tube Side Outlet Nozzle"
1453        Unity(i).Tubes.PressureDrop.Vnozzle_out = HE.TubeVelocityNozzleout(Unity(i).Properties.Cold.Outlet.rho,Unity(i).Properties.Cold.Inlet.Fw);
1454
1455"Shell Pressure Drop Inlet Nozzle"
1456        Unity(i).Shell.PressureDrop.Pdnozzle_in = HE.DeltaPshellNozzleIn(Unity(i).Properties.Hot.Inlet.rho,Unity(i).Properties.Hot.Inlet.Fw);
1457
1458"Velocity Shell Side Inlet Nozzle"
1459        Unity(i).Shell.PressureDrop.Vnozzle_in  = HE.ShellVelocityNozzleIn(Unity(i).Properties.Hot.Inlet.rho,Unity(i).Properties.Hot.Inlet.Fw);
1460
1461"Shell Pressure Drop Outlet Nozzle"
1462        Unity(i).Shell.PressureDrop.Pdnozzle_out        =HE.DeltaPshellNozzleOut(Unity(i).Properties.Hot.Outlet.rho,Unity(i).Properties.Hot.Inlet.Fw);
1463
1464"Velocity Shell Side Outlet Nozzle"
1465        Unity(i).Shell.PressureDrop.Vnozzle_out =HE.ShellVelocityNozzleOut(Unity(i).Properties.Hot.Outlet.rho,Unity(i).Properties.Hot.Inlet.Fw);
1466
1467"Pressure Drop Hot Stream"
1468        Unity(i).Outlet.Hot.P  = Unity(i).Inlet.Hot.P- Unity(i).Shell.PressureDrop.Pdtotal;     
1469       
1470"Pressure Drop Cold Stream"
1471        Unity(i).Outlet.Cold.P  = Unity(i).Inlet.Cold.P - Unity(i).Tubes.PressureDrop.Pdtotal;
1472
1473
1474        else
1475
1476"Cold Wall Temperature"
1477        Unity(i).Properties.Cold.Wall.Twall  = (Unity(i).Properties.Hot.Average.T+Unity(i).Properties.Cold.Average.T)/2;
1478
1479"Hot Wall Temperature"
1480        Unity(i).Properties.Hot.Wall.Twall  = (Unity(i).Properties.Hot.Average.T+Unity(i).Properties.Cold.Average.T)/2;
1481
1482"Tube Side Velocity"
1483        Unity(i).Tubes.HeatTransfer.Vtube  = Unity(i).Properties.Hot.Inlet.Fw*Tpass/((Pi*Ditube*Ditube/4)*Unity(i).Properties.Hot.Average.rho*Ntt);
1484
1485"Tube Side Reynolds Number"
1486        Unity(i).Tubes.HeatTransfer.Re=(Unity(i).Properties.Hot.Average.rho*Unity(i).Tubes.HeatTransfer.Vtube*Ditube)/Unity(i).Properties.Hot.Average.Mu;
1487       
1488"Tube Side Prandtl Number"
1489        Unity(i).Tubes.HeatTransfer.PR=((Unity(i).Properties.Hot.Average.Cp/Unity(i).Properties.Hot.Average.Mw)*Unity(i).Properties.Hot.Average.Mu)/Unity(i).Properties.Hot.Average.K;
1490
1491"Tube Side Prandtl Number at Wall"
1492        Unity(i).Tubes.HeatTransfer.PRw=((Unity(i).Properties.Hot.Wall.Cp/Unity(i).Properties.Hot.Average.Mw)*Unity(i).Properties.Hot.Wall.Mu)/Unity(i).Properties.Hot.Wall.K;
1493
1494"Tube Side Film Coefficient"
1495        Unity(i).Tubes.HeatTransfer.htube= HE.TubeFilmCoeff(Unity(i).Tubes.HeatTransfer.Re,Unity(i).Tubes.HeatTransfer.PR,Unity(i).Properties.Hot.Average.K)*Unity(i).Tubes.HeatTransfer.Phi;
1496       
1497"Shell Side Prandtl Number"
1498        Unity(i).Shell.HeatTransfer.PR=((Unity(i).Properties.Cold.Average.Cp/Unity(i).Properties.Cold.Average.Mw)*Unity(i).Properties.Cold.Average.Mu)/Unity(i).Properties.Cold.Average.K;
1499
1500"Shell Side Prandtl Number at Wall"
1501        Unity(i).Shell.HeatTransfer.PRw=((Unity(i).Properties.Cold.Wall.Cp/Unity(i).Properties.Cold.Average.Mw)*Unity(i).Properties.Cold.Wall.Mu)/Unity(i).Properties.Cold.Wall.K;
1502
1503"Tube Side Pressure Drop"
1504        Unity(i).Tubes.PressureDrop.PdTube              = HE.DeltaPtube(Unity(i).Tubes.HeatTransfer.Re,Unity(i).Properties.Hot.Average.rho,Unity(i).Tubes.HeatTransfer.Vtube,Unity(i).Tubes.HeatTransfer.Phi);
1505
1506"Pressure Drop Tube Side Inlet Nozzle"
1507        Unity(i).Tubes.PressureDrop.Pdnozzle_in = HE.DeltaPtubeNozzlein(Unity(i).Properties.Hot.Inlet.rho,Unity(i).Properties.Hot.Inlet.Fw);
1508
1509"Velocity Tube Side Inlet Nozzle"
1510        Unity(i).Tubes.PressureDrop.Vnozzle_in  = HE.TubeVelocityNozzlein(Unity(i).Properties.Hot.Inlet.rho,Unity(i).Properties.Hot.Inlet.Fw);
1511
1512"Pressure Drop Tube Side Outlet Nozzle"
1513        Unity(i).Tubes.PressureDrop.Pdnozzle_out        = HE.DeltaPtubeNozzleout(Unity(i).Properties.Hot.Outlet.rho,Unity(i).Properties.Hot.Inlet.Fw); 
1514
1515"Velocity Tube Side Outlet Nozzle"
1516        Unity(i).Tubes.PressureDrop.Vnozzle_out = HE.TubeVelocityNozzleout(Unity(i).Properties.Hot.Outlet.rho,Unity(i).Properties.Hot.Inlet.Fw);       
1517
1518"Shell Pressure Drop Inlet Nozzle"
1519        Unity(i).Shell.PressureDrop.Pdnozzle_in = HE.DeltaPshellNozzleIn(Unity(i).Properties.Cold.Inlet.rho,Unity(i).Properties.Cold.Inlet.Fw);
1520
1521"Velocity Shell Side Inlet Nozzle"
1522        Unity(i).Shell.PressureDrop.Vnozzle_in  = HE.ShellVelocityNozzleIn(Unity(i).Properties.Cold.Inlet.rho,Unity(i).Properties.Cold.Inlet.Fw);
1523
1524"Shell Pressure Drop Outlet Nozzle"
1525        Unity(i).Shell.PressureDrop.Pdnozzle_out = HE.DeltaPshellNozzleOut(Unity(i).Properties.Cold.Outlet.rho,Unity(i).Properties.Cold.Inlet.Fw);
1526
1527"Velocity Shell Side Outlet Nozzle"
1528        Unity(i).Shell.PressureDrop.Vnozzle_out  = HE.ShellVelocityNozzleOut(Unity(i).Properties.Cold.Outlet.rho,Unity(i).Properties.Cold.Inlet.Fw);
1529
1530"Pressure Drop Hot Stream"
1531        Unity(i).Outlet.Hot.P  = Unity(i).Inlet.Hot.P- Unity(i).Tubes.PressureDrop.Pdtotal;     
1532       
1533"Pressure Drop Cold Stream"
1534        Unity(i).Outlet.Cold.P  = Unity(i).Inlet.Cold.P - Unity(i).Shell.PressureDrop.Pdtotal;
1535
1536
1537end
1538
1539if side equal 1
1540       
1541        then
1542
1543"Shell Side inlet Nozzle rho-V^2"
1544        Unity(i).Shell.PressureDrop.RVsquare_in = Unity(i).Properties.Hot.Inlet.rho*(Unity(i).Shell.PressureDrop.Vnozzle_in)^2;
1545
1546"Shell Side Outlet Nozzle rho-V^2"
1547        Unity(i).Shell.PressureDrop.RVsquare_out = Unity(i).Properties.Hot.Outlet.rho*(Unity(i).Shell.PressureDrop.Vnozzle_out)^2;
1548       
1549        else
1550       
1551"Shell Side inlet Nozzle rho-V^2"
1552        Unity(i).Shell.PressureDrop.RVsquare_in = Unity(i).Properties.Cold.Inlet.rho*(Unity(i).Shell.PressureDrop.Vnozzle_in)^2;
1553
1554"Shell Side Outlet Nozzle rho-V^2"
1555        Unity(i).Shell.PressureDrop.RVsquare_out = Unity(i).Properties.Cold.Outlet.rho*(Unity(i).Shell.PressureDrop.Vnozzle_out)^2;
1556
1557end
1558
1559"Tube Resistance"       
1560        Unity(i).Resistances.Rtube*(Unity(i).Tubes.HeatTransfer.htube*Ditube) = Dotube;
1561       
1562"Wall Resistance"
1563        Unity(i).Resistances.Rwall=Dotube*ln(Dotube/Ditube)/(2*Kwall);
1564       
1565"Shell Resistance"
1566        Unity(i).Resistances.Rshell*(Unity(i).Shell.HeatTransfer.hshell)=1;
1567       
1568"Overall Heat Transfer Coefficient Clean"
1569        Unity(i).Details.Uc*(Dotube/(Unity(i).Tubes.HeatTransfer.htube*Ditube)+(Dotube*ln(Dotube/Ditube)/(2*Kwall))+(1/(Unity(i).Shell.HeatTransfer.hshell)))=1;
1570
1571"Overall Heat Transfer Coefficient Dirty"
1572        (1/Unity(i).Details.Ud)=(1/Unity(i).Details.Uc)+Unity(i).Resistances.Rfo+Unity(i).Resistances.Rfi*(Dotube/Ditube);
1573
1574"Exchange Surface Area"
1575        Unity(i).Details.A=Pi*Dotube*Ntt*Ltube;
1576
1577"Baffles Spacing"
1578        Ltube = Unity(i).Baffles.Lsi+Unity(i).Baffles.Lso+Unity(i).Baffles.Ls*(Nb-1);
1579       
1580"LMTD Correction Factor"
1581        Fc(i) = HE.EshellCorrectionFactor(Unity(i).Inlet.Hot.T,Unity(i).Outlet.Hot.T,Unity(i).Inlet.Cold.T,Unity(i).Outlet.Cold.T);
1582
1583"Counter Flow LMTD"
1584        LMTD(i) = HE.CounterLMTD(Unity(i).Inlet.Hot.T,Unity(i).Outlet.Hot.T,Unity(i).Inlet.Cold.T,Unity(i).Outlet.Cold.T);
1585
1586"Js Factor"     
1587        Unity(i).Shell.HeatTransfer.Js = HE.JsFactor(Unity(i).Shell.HeatTransfer.Re,Unity(i).Baffles.Lsi,Unity(i).Baffles.Lso,Unity(i).Baffles.Ls);
1588
1589if side equal 1
1590       
1591        then
1592       
1593"Shell Side Reynolds Number"
1594        Unity(i).Shell.HeatTransfer.Re=(Dotube*Unity(i).Properties.Hot.Inlet.Fw/Unity(i).Shell.HeatTransfer.Sm)/Unity(i).Properties.Hot.Average.Mu;
1595
1596"Shell Heat Transfer Coefficient"
1597        Unity(i).Shell.HeatTransfer.hshell=Unity(i).Shell.HeatTransfer.Ji*(Unity(i).Properties.Hot.Average.Cp/Unity(i).Properties.Hot.Average.Mw)*(Unity(i).Properties.Hot.Inlet.Fw/Unity(i).Shell.HeatTransfer.Sm)*(Unity(i).Shell.HeatTransfer.PR^(-2/3))*Unity(i).Shell.HeatTransfer.Jtotal*Unity(i).Shell.HeatTransfer.Phi;
1598
1599"Shell Pressure Drop Cross Flow"
1600        Unity(i).Shell.PressureDrop.PdCross             = HE.DeltaPcross(Unity(i).Shell.HeatTransfer.Re,Unity(i).Baffles.Ls,Unity(i).Baffles.Lso,Unity(i).Baffles.Lsi,Unity(i).Properties.Hot.Inlet.Fw,Unity(i).Shell.HeatTransfer.Phi,Unity(i).Properties.Hot.Average.rho);
1601
1602"Shell Pressure Baffle Window"
1603        Unity(i).Shell.PressureDrop.Pdwindow            = HE.DeltaPwindow(Unity(i).Properties.Hot.Inlet.Fw,Unity(i).Shell.HeatTransfer.Sm,Unity(i).Properties.Hot.Average.rho,Unity(i).Properties.Hot.Average.Mu,Unity(i).Baffles.Ls);
1604
1605"Shell Pressure End Zones"
1606        Unity(i).Shell.PressureDrop.PdEndZones  = HE.DeltaPendZones(Unity(i).Shell.HeatTransfer.Re,Unity(i).Baffles.Ls,Unity(i).Baffles.Lso,Unity(i).Baffles.Lsi,Unity(i).Properties.Hot.Inlet.Fw,Unity(i).Shell.HeatTransfer.Phi,Unity(i).Properties.Hot.Average.rho);
1607
1608
1609        else
1610
1611"Shell Side Reynolds Number"
1612        Unity(i).Shell.HeatTransfer.Re=(Dotube*Unity(i).Properties.Cold.Inlet.Fw/Unity(i).Shell.HeatTransfer.Sm)/Unity(i).Properties.Cold.Average.Mu;
1613
1614"Shell Heat Transfer Coefficient"
1615        Unity(i).Shell.HeatTransfer.hshell=Unity(i).Shell.HeatTransfer.Ji*(Unity(i).Properties.Cold.Average.Cp/Unity(i).Properties.Cold.Average.Mw)*(Unity(i).Properties.Cold.Inlet.Fw/Unity(i).Shell.HeatTransfer.Sm)*(Unity(i).Shell.HeatTransfer.PR^(-2/3))*Unity(i).Shell.HeatTransfer.Jtotal*Unity(i).Shell.HeatTransfer.Phi;
1616
1617"Shell Pressure Drop Cross Flow"
1618        Unity(i).Shell.PressureDrop.PdCross      = HE.DeltaPcross(Unity(i).Shell.HeatTransfer.Re,Unity(i).Baffles.Ls,Unity(i).Baffles.Lso,Unity(i).Baffles.Lsi,Unity(i).Properties.Cold.Inlet.Fw,Unity(i).Shell.HeatTransfer.Phi,Unity(i).Properties.Cold.Average.rho);
1619
1620"Shell Pressure Baffle Window"
1621        Unity(i).Shell.PressureDrop.Pdwindow     = HE.DeltaPwindow(Unity(i).Properties.Cold.Inlet.Fw,Unity(i).Shell.HeatTransfer.Sm,Unity(i).Properties.Cold.Average.rho,Unity(i).Properties.Cold.Average.Mu,Unity(i).Baffles.Ls);
1622
1623"Shell Pressure End Zones"
1624        Unity(i).Shell.PressureDrop.PdEndZones  = HE.DeltaPendZones(Unity(i).Shell.HeatTransfer.Re,Unity(i).Baffles.Ls,Unity(i).Baffles.Lso,Unity(i).Baffles.Lsi,Unity(i).Properties.Cold.Inlet.Fw,Unity(i).Shell.HeatTransfer.Phi,Unity(i).Properties.Cold.Average.rho);
1625
1626
1627end
1628
1629"Exchange Surface Area"
1630        Unity(i).Details.Q   = Unity(i).Details.Ud*Pi*Dotube*Ntt*Ltube*Fc(i)*LMTD(i);
1631       
1632"Mean Temperature Difference"   
1633        MTD(i)   = Fc(i)*LMTD(i);
1634
1635end
1636
1637        SET
1638#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
1639#  Set Parameters for heatex Calculation
1640#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
1641        Pi                              = 3.14159265;
1642        HE.Tpass                        = Tpass;
1643        HE.Nss                          = Nss;
1644        HE.Ntt                          = Ntt;
1645        HE.Pattern                      = Pattern;
1646        HE.Bc                           = Bc;
1647        HE.Donozzle_Shell       = Donozzle_Shell;
1648        HE.Dinozzle_Shell       = Dinozzle_Shell;
1649        HE.Honozzle_Shell       = Honozzle_Shell;
1650        HE.Hinozzle_Shell       = Hinozzle_Shell;
1651        HE.Donozzle_Tube        = Donozzle_Tube;
1652        HE.Dinozzle_Tube        = Dinozzle_Tube;
1653        HE.Nb                       = Nb;
1654        HE.Dishell                      = Dishell;             
1655        HE.Lcf                          = Lcf; 
1656        HE.Ltube                        = Ltube;       
1657        HE.pitch                        = pitch;               
1658        HE.Dotube                       = Dotube;       
1659        HE.Ditube                       = Ditube;       
1660        HE.Lcd                          = Lcd;
1661        HE.Ltd                          = Ltd; 
1662        side                            = HE.FluidAlocation();
1663
1664end
1665
Note: See TracBrowser for help on using the repository browser.