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 | |
---|
20 | using "HEX_Engine"; |
---|
21 | |
---|
22 | #+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++# |
---|
23 | # Basic Model for Detailed Shell and Tubes Heat Exchangers |
---|
24 | #+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++# |
---|
25 | |
---|
26 | Model HeatExchangerDetailed_Basic |
---|
27 | |
---|
28 | PARAMETERS |
---|
29 | |
---|
30 | ext PP as CalcObject (Brief="External Physical Properties"); |
---|
31 | ext 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 | |
---|
35 | VARIABLES |
---|
36 | |
---|
37 | in Inlet as Inlet_Main_Stream; |
---|
38 | out 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 | |
---|
46 | SET |
---|
47 | M = PP.MolecularWeight(); |
---|
48 | |
---|
49 | EQUATIONS |
---|
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 | |
---|
74 | if 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 | |
---|
174 | end |
---|
175 | |
---|
176 | if 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 | |
---|
274 | end |
---|
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 | |
---|
341 | end |
---|
342 | |
---|
343 | Model Heatex_Detailed as HeatExchangerDetailed_Basic |
---|
344 | |
---|
345 | PARAMETERS |
---|
346 | |
---|
347 | #+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++# |
---|
348 | # Basic Parameters # |
---|
349 | #+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++# |
---|
350 | ext PP as CalcObject (Brief="External Physical Properties"); |
---|
351 | side as Integer (Brief="Fluid Alocation Flag",Lower=0,Upper=1); |
---|
352 | Pi as constant (Brief="Pi Number",Default=3.14159265); |
---|
353 | #+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++# |
---|
354 | # Shell Geometrical Parameters # |
---|
355 | #+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++# |
---|
356 | Tpass as Integer (Brief="Number of Tube Passes",Lower=1); |
---|
357 | Nss as Integer (Brief="Number of Sealing Strips pairs",Lower=1); |
---|
358 | Dishell as length (Brief="Inside Shell Diameter",Lower=10e-6); |
---|
359 | Donozzle_Shell as length (Brief="Shell Outlet Nozzle Diameter",Lower=10e-6); |
---|
360 | Dinozzle_Shell as length (Brief="Shell Inlet Nozzle Diameter",Lower=10e-6); |
---|
361 | Hinozzle_Shell as length (Brief="Height Under Shell Inlet Nozzle",Lower=10e-6); |
---|
362 | Honozzle_Shell as length (Brief="Height Under Shell Outlet Nozzle",Lower=10e-6); |
---|
363 | Lcf as length (Brief="Bundle-to-Shell Clearance",Lower=10e-8); |
---|
364 | |
---|
365 | #+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++# |
---|
366 | # Tubes Geometrical Parameters # |
---|
367 | #+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++# |
---|
368 | Ntt as Integer (Brief="Total Number of Tubes in Shell",Default=100,Lower=1); |
---|
369 | Pattern as Integer (Brief="Tube Layout Characteristic Angle",Lower=30); |
---|
370 | Ltube as length (Brief="Effective Tube Length",Lower=0.1); |
---|
371 | pitch as length (Brief="Tube Pitch",Lower=1e-8); |
---|
372 | Kwall as conductivity (Brief="Tube Wall Material Thermal Conductivity"); |
---|
373 | Dotube as length (Brief="Tube Outside Diameter",Lower=10e-6); |
---|
374 | Ditube as length (Brief="Tube Inside Diameter",Lower=10e-6); |
---|
375 | Donozzle_Tube as length (Brief="Tube Outlet Nozzle Diameter",Lower=10e-6); |
---|
376 | Dinozzle_Tube as length (Brief="Tube Inlett Nozzle Diameter",Lower=10e-6); |
---|
377 | |
---|
378 | #+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++# |
---|
379 | # Baffles Geometrical Parameters # |
---|
380 | #+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++# |
---|
381 | Bc as Integer (Brief="Baffle Cut",Default=25,Lower=25); |
---|
382 | Nb as Real (Brief="Number of Baffles",Lower=1); |
---|
383 | Lcd as length (Brief="Baffle-to-Shell Clearance",Lower=10e-8); |
---|
384 | Ltd as length (Brief="Tube-to-Bafflehole Clearance",Lower=10e-8); |
---|
385 | |
---|
386 | SET |
---|
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 | |
---|
414 | EQUATIONS |
---|
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 | |
---|
441 | if 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 | |
---|
459 | end |
---|
460 | |
---|
461 | #+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++# |
---|
462 | # Phi correction # |
---|
463 | #+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++# |
---|
464 | |
---|
465 | if 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 | |
---|
483 | end |
---|
484 | |
---|
485 | #+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++# |
---|
486 | # # |
---|
487 | # # |
---|
488 | # # |
---|
489 | #+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++# |
---|
490 | |
---|
491 | if 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 | |
---|
618 | end |
---|
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 | |
---|
654 | end |
---|
655 | |
---|
656 | Model Heatex_Detailed_NTU as Heatex_Detailed |
---|
657 | |
---|
658 | VARIABLES |
---|
659 | |
---|
660 | Eft as positive (Brief="Effectiveness",Default=0.05,Lower=1e-8); |
---|
661 | |
---|
662 | EQUATIONS |
---|
663 | "Duty" |
---|
664 | Details.Q = Eft*Details.Cmin*(Inlet.Hot.T-Inlet.Cold.T); |
---|
665 | |
---|
666 | end |
---|
667 | |
---|
668 | Model Heatex_Detailed_LMTD as Heatex_Detailed |
---|
669 | |
---|
670 | VARIABLES |
---|
671 | |
---|
672 | DT0 as temp_delta (Brief="Temperature Difference at Inlet",Lower=1); |
---|
673 | DTL as temp_delta (Brief="Temperature Difference at Outlet",Lower=1); |
---|
674 | LMTD as temp_delta (Brief="Logarithmic Mean Temperature Difference",Lower=5); |
---|
675 | Fc as positive (Brief="LMTD Correction Factor",Lower=0.4); |
---|
676 | MTD as temp_delta (Brief="Mean Temperature Difference",Lower=5); |
---|
677 | |
---|
678 | EQUATIONS |
---|
679 | #+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++# |
---|
680 | # Log Mean Temperature Difference |
---|
681 | #+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++# |
---|
682 | if 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 | |
---|
690 | if 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 | |
---|
700 | end |
---|
701 | |
---|
702 | end |
---|
703 | |
---|
704 | "Exchange Surface Area" |
---|
705 | Details.Q = Details.Ud*Details.A*MTD; |
---|
706 | |
---|
707 | "Mean Temperature Difference" |
---|
708 | MTD = Fc*LMTD; |
---|
709 | |
---|
710 | end |
---|
711 | |
---|
712 | Model E_Shell_NTU_Det as Heatex_Detailed_NTU |
---|
713 | #+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++# |
---|
714 | # Shell and Tubes Heat Exchanger with 1 shell pass - LMTD Method |
---|
715 | #+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++# |
---|
716 | EQUATIONS |
---|
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 | |
---|
725 | if 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 | |
---|
763 | end |
---|
764 | |
---|
765 | end |
---|
766 | |
---|
767 | Model F_Shell_NTU_Det as Heatex_Detailed_NTU |
---|
768 | #+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++# |
---|
769 | # Shell and Tubes Heat Exchanger with 1 shell pass - LMTD Method |
---|
770 | #+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++# |
---|
771 | VARIABLES |
---|
772 | |
---|
773 | Eft1 as positive (Brief="Effectiveness Correction",Lower=0.01,Default=0.5); |
---|
774 | |
---|
775 | EQUATIONS |
---|
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 | |
---|
786 | if 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 | |
---|
824 | end |
---|
825 | |
---|
826 | |
---|
827 | end |
---|
828 | |
---|
829 | Model Multipass_NTU_Det |
---|
830 | #+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++# |
---|
831 | # Shell and Tubes Heat Exchanger In Series with 1 shell pass - LMTD Method |
---|
832 | #+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++# |
---|
833 | PARAMETERS |
---|
834 | |
---|
835 | Nshell as Integer (Brief="N Shell in Series",Default=2); |
---|
836 | ext 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 | #+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++# |
---|
844 | Tpass as Integer (Brief="Number of Tube Passes",Lower=1); |
---|
845 | Nss as Integer (Brief="Number of Sealing Strips pairs",Lower=1); |
---|
846 | Dishell as length (Brief="Inside Shell Diameter",Lower=10e-6); |
---|
847 | Donozzle_Shell as length (Brief="Shell Outlet Nozzle Diameter",Lower=10e-6); |
---|
848 | Dinozzle_Shell as length (Brief="Shell Inlet Nozzle Diameter",Lower=10e-6); |
---|
849 | Hinozzle_Shell as length (Brief="Height Under Shell Inlet Nozzle",Lower=10e-6); |
---|
850 | Honozzle_Shell as length (Brief="Height Under Shell Outlet Nozzle",Lower=10e-6); |
---|
851 | Lcf as length (Brief="Bundle-to-Shell Clearance",Lower=10e-8); |
---|
852 | #+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++# |
---|
853 | # Tubes Geometrical Parameters |
---|
854 | #+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++# |
---|
855 | Ntt as Integer (Brief="Total Number of Tubes in Shell",Default=100,Lower=1); |
---|
856 | Pattern as Integer (Brief="Tube Layout Characteristic Angle",Lower=30); |
---|
857 | Ltube as length (Brief="Effective Tube Length",Lower=0.1); |
---|
858 | pitch as length (Brief="Tube Pitch",Lower=1e-8); |
---|
859 | Kwall as conductivity (Brief="Tube Wall Material Thermal Conductivity"); |
---|
860 | Dotube as length (Brief="Tube Outside Diameter",Lower=10e-6); |
---|
861 | Ditube as length (Brief="Tube Inside Diameter",Lower=10e-6); |
---|
862 | Donozzle_Tube as length (Brief="Tube Outlet Nozzle Diameter",Lower=10e-6); |
---|
863 | Dinozzle_Tube as length (Brief="Tube Inlett Nozzle Diameter",Lower=10e-6); |
---|
864 | #+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++# |
---|
865 | # Baffles Geometrical Parameters |
---|
866 | #+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++# |
---|
867 | Bc as Integer (Brief="Baffle Cut",Default=25,Lower=25); |
---|
868 | Nb as Real (Brief="Number of Baffles",Lower=1); |
---|
869 | Lcd as length (Brief="Baffle-to-Shell Clearance",Lower=10e-8); |
---|
870 | Ltd as length (Brief="Tube-to-Bafflehole Clearance",Lower=10e-8); |
---|
871 | #+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++# |
---|
872 | VARIABLES |
---|
873 | |
---|
874 | Eft(Nshell) as positive (Brief="Effectiveness",Default=0.05,Lower=1e-8); |
---|
875 | Unity(Nshell) as HeatExchangerDetailed_Basic;# "Shell in Series" |
---|
876 | |
---|
877 | CONNECTIONS |
---|
878 | |
---|
879 | Unity([1:Nshell-1]).Outlet.Hot to Unity([2:Nshell]).Inlet.Hot; |
---|
880 | Unity([2:Nshell]).Outlet.Cold to Unity([1:Nshell-1]).Inlet.Cold; |
---|
881 | |
---|
882 | EQUATIONS |
---|
883 | |
---|
884 | for 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 | |
---|
911 | if 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 | |
---|
927 | end |
---|
928 | |
---|
929 | |
---|
930 | if 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 | |
---|
1058 | end |
---|
1059 | |
---|
1060 | if 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 | |
---|
1078 | end |
---|
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 | |
---|
1107 | if 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 | |
---|
1145 | end |
---|
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 | |
---|
1150 | end |
---|
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 | |
---|
1179 | end |
---|
1180 | |
---|
1181 | Model E_Shell_LMTD_Det as Heatex_Detailed_LMTD |
---|
1182 | #+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++# |
---|
1183 | # Shell and Tubes Heat Exchanger with 1 shell pass - LMTD Method |
---|
1184 | #+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++# |
---|
1185 | EQUATIONS |
---|
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 | |
---|
1199 | if 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 | |
---|
1237 | end |
---|
1238 | |
---|
1239 | |
---|
1240 | |
---|
1241 | end |
---|
1242 | |
---|
1243 | Model F_Shell_LMTD_Det as Heatex_Detailed_LMTD |
---|
1244 | #+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++# |
---|
1245 | # Shell and Tubes Heat Exchanger with 2 shell pass - LMTD Method |
---|
1246 | #+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++# |
---|
1247 | |
---|
1248 | EQUATIONS |
---|
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 | |
---|
1262 | if 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 | |
---|
1300 | end |
---|
1301 | |
---|
1302 | end |
---|
1303 | |
---|
1304 | Model Multipass_LMTD_Det |
---|
1305 | #+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++# |
---|
1306 | # Shell and Tubes Heat Exchanger In Series with 1 shell pass - LMTD Method |
---|
1307 | #+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++# |
---|
1308 | PARAMETERS |
---|
1309 | |
---|
1310 | Nshell as Integer (Brief="N Shell in Series",Default=2); |
---|
1311 | HE as CalcObject (Brief="STHE Calculations",File="heatex"); |
---|
1312 | ext PP as CalcObject (Brief="External Physical Properties"); |
---|
1313 | side as Integer (Brief="Fluid Alocation",Lower=0,Upper=1); |
---|
1314 | Pi as constant (Brief="Pi Number",Default=3.14159265); |
---|
1315 | |
---|
1316 | #+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++# |
---|
1317 | # Shell Geometrical Parameters |
---|
1318 | #+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++# |
---|
1319 | Tpass as Integer (Brief="Number of Tube Passes",Lower=1); |
---|
1320 | Nss as Integer (Brief="Number of Sealing Strips pairs",Lower=1); |
---|
1321 | Dishell as length (Brief="Inside Shell Diameter",Lower=10e-6); |
---|
1322 | Donozzle_Shell as length (Brief="Shell Outlet Nozzle Diameter",Lower=10e-6); |
---|
1323 | Dinozzle_Shell as length (Brief="Shell Inlet Nozzle Diameter",Lower=10e-6); |
---|
1324 | Hinozzle_Shell as length (Brief="Height Under Shell Inlet Nozzle",Lower=10e-6); |
---|
1325 | Honozzle_Shell as length (Brief="Height Under Shell Outlet Nozzle",Lower=10e-6); |
---|
1326 | Lcf as length (Brief="Bundle-to-Shell Clearance",Lower=10e-8); |
---|
1327 | #+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++# |
---|
1328 | # Tubes Geometrical Parameters |
---|
1329 | #+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++# |
---|
1330 | Ntt as Integer (Brief="Total Number of Tubes in Shell",Default=100,Lower=1); |
---|
1331 | Pattern as Integer (Brief="Tube Layout Characteristic Angle",Lower=30); |
---|
1332 | Ltube as length (Brief="Effective Tube Length",Lower=0.1); |
---|
1333 | pitch as length (Brief="Tube Pitch",Lower=1e-8); |
---|
1334 | Kwall as conductivity (Brief="Tube Wall Material Thermal Conductivity"); |
---|
1335 | Dotube as length (Brief="Tube Outside Diameter",Lower=10e-6); |
---|
1336 | Ditube as length (Brief="Tube Inside Diameter",Lower=10e-6); |
---|
1337 | Donozzle_Tube as length (Brief="Tube Outlet Nozzle Diameter",Lower=10e-6); |
---|
1338 | Dinozzle_Tube as length (Brief="Tube Inlett Nozzle Diameter",Lower=10e-6); |
---|
1339 | #+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++# |
---|
1340 | # Baffles Geometrical Parameters |
---|
1341 | #+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++# |
---|
1342 | Bc as Integer (Brief="Baffle Cut",Default=25,Lower=25); |
---|
1343 | Nb as Real (Brief="Number of Baffles",Lower=1); |
---|
1344 | Lcd as length (Brief="Baffle-to-Shell Clearance",Lower=10e-8); |
---|
1345 | Ltd as length (Brief="Tube-to-Bafflehole Clearance",Lower=10e-8); |
---|
1346 | #+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++# |
---|
1347 | VARIABLES |
---|
1348 | |
---|
1349 | LMTD(Nshell) as temp_delta (Brief="Logarithmic Mean Temperature Difference",Lower=5); |
---|
1350 | Fc(Nshell) as positive (Brief="LMTD Correction Factor",Lower=0.5); |
---|
1351 | MTD(Nshell) as temp_delta (Brief="Mean Temperature Difference",Lower=5); |
---|
1352 | Unity(Nshell) as HeatExchangerDetailed_Basic; |
---|
1353 | |
---|
1354 | CONNECTIONS |
---|
1355 | |
---|
1356 | Unity([1:Nshell-1]).Outlet.Hot to Unity([2:Nshell]).Inlet.Hot; |
---|
1357 | Unity([2:Nshell]).Outlet.Cold to Unity([1:Nshell-1]).Inlet.Cold; |
---|
1358 | |
---|
1359 | EQUATIONS |
---|
1360 | |
---|
1361 | for 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 | |
---|
1388 | if 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 | |
---|
1406 | end |
---|
1407 | |
---|
1408 | |
---|
1409 | if 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 | |
---|
1537 | end |
---|
1538 | |
---|
1539 | if 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 | |
---|
1557 | end |
---|
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 | |
---|
1589 | if 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 | |
---|
1627 | end |
---|
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 | |
---|
1635 | end |
---|
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 | |
---|
1664 | end |
---|
1665 | |
---|