Changeset 523 for trunk/eml/heat_exchangers/PHE.mso
- Timestamp:
- May 23, 2008, 4:44:08 PM (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/eml/heat_exchangers/PHE.mso
r420 r523 18 18 using "HEX_Engine"; 19 19 20 Model PHE_PressureDrop 21 22 ATTRIBUTES 23 Pallete = false; 24 Brief = "to be documented"; 25 Info = 26 "to be documented"; 27 28 VARIABLES 29 30 DPchannel as press_delta (Brief="Channel Pressure Drop",Default=0.01, Lower=1E10,DisplayUnit='kPa', Symbol ="\Delta P^{channel}"); 31 DPports as press_delta (Brief="Ports Pressure Drop",Default=0.01, Lower=1E-10,DisplayUnit='kPa', Symbol ="\Delta P^{ports}"); 32 Pdrop as press_delta (Brief="Total Pressure Drop",Default=0.01, Lower=1E-10,DisplayUnit='kPa', Symbol ="\Delta P"); 33 fi as fricfactor (Brief="Friction Factor", Default=0.05, Lower=1E-10, Upper=2000); 34 Vchannel as velocity (Brief="Stream Velocity in Channel",Lower=1E-8, Symbol ="V^{channel}"); 35 Vports as velocity (Brief="Stream Velocity in Ports",Lower=1E-8, Symbol ="V^{ports}"); 36 Npassage as positive (Brief="Number of Channels per Pass", Symbol ="N^{passage}"); 37 38 end 39 40 Model PHE_HeatTransfer 41 42 ATTRIBUTES 43 Pallete = false; 44 Brief = "to be documented"; 45 Info = 46 "to be documented"; 47 48 VARIABLES 49 50 Re as positive (Brief="Reynolds Number",Default=100,Lower=1); 51 PR as positive (Brief="Prandtl Number",Default=0.5,Lower=1e-8); 52 NTU as positive (Brief="Number of Units Transference",Default=0.05,Lower=1E-10); 53 WCp as positive (Brief="Stream Heat Capacity",Lower=1E-3,Default=1E3,Unit='W/K'); 54 hcoeff as heat_trans_coeff (Brief="Film Coefficient",Default=1,Lower=1E-12, Upper=1E6); 55 Gchannel as flux_mass (Brief ="Channel Mass Flux", Default=1, Lower=1E-6, Symbol ="G^{channel}"); 56 Gports as flux_mass (Brief ="Ports Mass Flux", Default=1, Lower=1E-6, Symbol ="G^{ports}"); 57 Phi as positive (Brief="Viscosity Correction",Default=1,Lower=1E-6, Symbol="\phi"); 58 59 end 60 61 Model Main_PHE 62 63 ATTRIBUTES 64 Pallete = false; 65 Brief = "to be documented"; 66 Info = 67 "to be documented"; 68 69 VARIABLES 70 71 HeatTransfer as PHE_HeatTransfer (Brief="PHE Heat Transfer", Symbol = " "); 72 PressureDrop as PHE_PressureDrop (Brief="PHE Pressure Drop", Symbol = " "); 73 Properties as Physical_Properties (Brief="PHE Properties", Symbol = " "); 74 75 end 76 77 Model Thermal_PHE 78 79 ATTRIBUTES 80 Pallete = false; 81 Brief = "to be documented"; 82 Info = 83 "to be documented"; 84 85 VARIABLES 86 Cr as positive (Brief="Heat Capacity Ratio",Default=0.5,Lower=1E-6); 87 Cmin as positive (Brief="Minimum Heat Capacity",Lower=1E-10,Default=1E3,Unit='W/K'); 88 Cmax as positive (Brief="Maximum Heat Capacity",Lower=1E-10,Default=1E3,Unit='W/K'); 89 NTU as positive (Brief="Number of Units Transference",Default=0.05,Lower=1E-10); 90 Eft as positive (Brief="Effectiveness",Default=0.5,Lower=0.1,Upper=1.1, Symbol = "\varepsilon"); 91 Q as power (Brief="Heat Transfer", Default=7000, Lower=1E-6, Upper=1E10); 92 Uc as heat_trans_coeff (Brief="Overall Heat Transfer Coefficient Clean",Default=1,Lower=1E-6,Upper=1E10); 93 Ud as heat_trans_coeff (Brief="Overall Heat Transfer Coefficient Dirty",Default=1,Lower=1E-6,Upper=1E10); 94 95 end 96 97 Model PHE_Geometry 98 99 ATTRIBUTES 100 Pallete = false; 101 Brief = "Parameters for a gasketed plate heat exchanger."; 102 103 PARAMETERS 104 105 outer PP as Plugin (Brief="External Physical Properties", Type="PP"); 106 outer NComp as Integer (Brief="Number of Chemical Components",Hidden=true); 107 108 Pi as constant (Brief="Pi Number",Default=3.14159265, Hidden=true,Symbol = "\pi"); 109 N1 as Integer (Brief="Auxiliar Constant", Hidden=true,Default = 15); 110 N2 as Integer (Brief="Auxiliar Constant",Hidden=true,Default = 14); 111 Kp1(N1) as constant (Brief="First constant in Kumar calculation for Pressure Drop", Hidden=true); 112 Kp2(N1) as constant (Brief="Second constant in Kumar calculation for Pressure Drop", Hidden=true); 113 Kc1(N2) as constant (Brief="First constant in Kumar calculation for Heat Transfer", Hidden=true); 114 Kc2(N2) as constant (Brief="Second constant Kumar calculation for Heat Transfer", Hidden=true); 115 M(NComp) as molweight (Brief="Component Mol Weight", Hidden=true); 116 117 118 Lv as length (Brief="Vertical Ports Distance",Lower=0.1); 119 Nplates as Integer (Brief="Total Number of Plates in The Whole Heat Exchanger",Default=25, Symbol ="N_{plates}"); 120 NpassHot as Integer (Brief="Number of Passes for Hot Side", Symbol ="Npasshot"); 121 NpassCold as Integer (Brief="Number of Passes for Cold Side", Symbol ="Npasscold"); 122 Dports as length (Brief="Ports Diameter",Lower=1e-6, Symbol ="D_{ports}"); 123 Lw as length (Brief="Plate Width",Lower=0.1); 124 pitch as length (Brief="Plate Pitch",Lower=0.1); 125 pt as length (Brief="Plate Thickness",Lower=0.1); 126 Kwall as conductivity (Brief="Plate Thermal Conductivity",Default=1.0, Symbol ="K_{wall}"); 127 Rfh as positive (Brief="Hot Side Fouling Resistance",Unit='m^2*K/kW',Default=1e-6,Lower=0); 128 Rfc as positive (Brief="Cold Side Fouling Resistance",Unit='m^2*K/kW',Default=1e-6,Lower=0); 129 PhiFactor as Real (Brief="Enlargement Factor",Lower=1e-6, Symbol ="\phi"); 130 131 Atotal as area (Brief="Total Effective Area",Lower=1e-6, Symbol ="A_{total}", Protected=true); 132 Aports as area (Brief="Port Opening Area of Plate",Lower=1e-6, Symbol ="A_{ports}", Protected=true); 133 Achannel as area (Brief="Cross-Sectional Area for Channel Flow",Lower=1e-6, Symbol ="A_{channel}", Protected=true); 134 Dh as length (Brief="Equivalent Diameter of Channel",Lower=1e-6, Protected=true); 135 Depth as length (Brief="Corrugation Depth",Lower=1e-6, Protected=true); 136 Nchannels as Integer (Brief="Total Number of Channels in The Whole Heat Exchanger", Protected=true); 137 Lp as length (Brief="Plate Vertical Distance between Port Centers",Lower=0.1, Protected=true); 138 Lpack as length (Brief="Compact Plate Pack Length",Lower=0.1, Protected=true); 139 Lh as length (Brief="Plate Horizontal Distance between Port Centers",Lower=0.1, Protected=true); 140 141 SET 142 143 #"Vector Length of constants for Kumar's calculating Pressure Drop" 144 N1 = 15; 145 146 #"Vector Length of constants for Kumar's calculating Heat Transfer" 147 N2 = 14; 148 149 #"First constant for Kumar's calculating Pressure Drop" 150 Kp1 = [50,19.40,2.990,47,18.290,1.441,34,11.250,0.772,24,3.240,0.760,24,2.80,0.639]; 151 152 #"Second constant for Kumar's calculating Pressure Drop" 153 Kp2 = [1,0.589,0.183,1,0.652,0.206,1,0.631,0.161,1,0.457,0.215,1,0.451,0.213]; 154 155 #"First constant for Kumar's calculating Heat Transfer" 156 Kc1 = [0.718,0.348,0.718,0.400,0.300,0.630,0.291,0.130,0.562,0.306,0.108,0.562,0.331,0.087]; 157 158 #"Second constant for Kumar's calculating Heat Transfer" 159 Kc2 = [0.349,0.663,0.349,0.598,0.663,0.333,0.591,0.732,0.326,0.529,0.703,0.326,0.503,0.718]; 160 161 #"Component Molecular Weight" 162 M = PP.MolecularWeight(); 163 164 #"Pi Number" 165 Pi = 3.14159265; 166 167 #"Plate Vertical Distance between Port Centers" 168 Lp = Lv - Dports; 169 170 #"Corrugation Depth" 171 Depth=pitch-pt; 172 173 #"Plate Horizontal Distance between Port Centers" 174 Lh=Lw-Dports; 175 176 #"Hydraulic Diameter" 177 Dh=2*Depth/PhiFactor; 178 179 #"Ports Area" 180 Aports=0.25*Pi*Dports*Dports; 181 182 #"Channel Area" 183 Achannel=Depth*Lw; 184 185 #"Pack Length" 186 Lpack=Depth*(Nplates-1)+Nplates*pt; 187 188 #"Total Number of Channels" 189 Nchannels = Nplates -1; 190 191 #"Exchange Surface Area" 192 Atotal =(Nplates-2)*Lw*Lp*PhiFactor; 193 194 end 195 20 196 Model PHE 21 197 … … 23 199 Icon = "icon/phe"; 24 200 Pallete = true; 25 Brief = "Shortcut model for plate and Frame heat exchanger.";201 Brief = "Shortcut model for Plate and Frame heat exchanger."; 26 202 Info = 27 203 "Model of a gasketed plate heat exchanger. … … 69 245 70 246 outer PP as Plugin (Brief="External Physical Properties", Type="PP"); 71 outer NComp as Integer (Brief="Number of Chemical Components"); 72 Pi as constant (Brief="Pi Number",Default=3.14159265, Symbol = "\pi"); 73 N1 as Integer (Brief="Auxiliar Constant",Default = 15); 74 N2 as Integer (Brief="Auxiliar Constant",Default = 14); 75 Kp1(N1) as constant (Brief="First constant in Kumar calculation for Pressure Drop"); 76 Kp2(N1) as constant (Brief="Second constant in Kumar calculation for Pressure Drop"); 77 Kc1(N2) as constant (Brief="First constant in Kumar calculation for Heat Transfer"); 78 Kc2(N2) as constant (Brief="Second constant Kumar calculation for Heat Transfer"); 79 M(NComp) as molweight (Brief="Component Mol Weight"); 80 81 ChevronAngle as Switcher (Brief="Chevron Corrugation Inclination Angle in Degrees ",Valid=["A30_Deg","A45_Deg","A50_Deg","A60_Deg","A65_Deg"],Default="A30_Deg"); 82 SideOne as Switcher (Brief="Fluid Alocation in the Side I - (The odd channels)",Valid=["hot","cold"],Default="hot"); 83 Nchannels as Integer (Brief="Total Number of Channels in The Whole Heat Exchanger"); 84 Nplates as Integer (Brief="Total Number of Plates in The Whole Heat Exchanger",Default=25, Symbol ="N_{plates}"); 85 NpassHot as Integer (Brief="Number of Passes for Hot Side", Symbol ="Npasshot"); 86 NpassCold as Integer (Brief="Number of Passes for Cold Side", Symbol ="Npasscold"); 87 Dports as length (Brief="Ports Diameter",Lower=1e-6, Symbol ="D_{ports}"); 88 Atotal as area (Brief="Total Effective Area",Lower=1e-6, Symbol ="A_{total}"); 89 Aports as area (Brief="Port Opening Area of Plate",Lower=1e-6, Symbol ="A_{ports}"); 90 Achannel as area (Brief="Cross-Sectional Area for Channel Flow",Lower=1e-6, Symbol ="A_{channel}"); 91 Dh as length (Brief="Equivalent Diameter of Channel",Lower=1e-6); 92 Depth as length (Brief="Corrugation Depth",Lower=1e-6); 93 PhiFactor as Real (Brief="Enlargement Factor",Lower=1e-6, Symbol ="\phi"); 94 Lp as length (Brief="Plate Vertical Distance between Port Centers",Lower=0.1); 95 Lpack as length (Brief="Compact Plate Pack Length",Lower=0.1); 96 Lv as length (Brief="Vertical Ports Distance",Lower=0.1); 97 Lh as length (Brief="Plate Horizontal Distance between Port Centers",Lower=0.1); 98 Lw as length (Brief="Plate Width",Lower=0.1); 99 pitch as length (Brief="Plate Pitch",Lower=0.1); 100 pt as length (Brief="Plate Thickness",Lower=0.1); 101 Kwall as conductivity (Brief="Plate Thermal Conductivity",Default=1.0, Symbol ="K_{wall}"); 102 Rfh as positive (Brief="Hot Side Fouling Resistance",Unit='m^2*K/kW',Default=1e-6,Lower=0); 103 Rfc as positive (Brief="Cold Side Fouling Resistance",Unit='m^2*K/kW',Default=1e-6,Lower=0); 247 outer NComp as Integer (Brief="Number of Chemical Components"); 248 249 ChevronAngle as Switcher (Brief="Chevron Corrugation Inclination Angle in Degrees ",Valid=["A30_Deg","A45_Deg","A50_Deg","A60_Deg","A65_Deg"],Default="A30_Deg"); 250 SideOne as Switcher (Brief="Fluid Alocation in the Side I - (The odd channels)",Valid=["hot","cold"],Default="hot"); 104 251 105 252 VARIABLES 106 253 107 in InletHot as stream (Brief="Inlet Hot Stream", PosX=0, PosY=0.75, Symbol="^{inHot}"); 108 in InletCold as stream (Brief="Inlet Cold Stream", PosX=0, PosY=0.25, Symbol="^{inCold}"); 109 out OutletHot as streamPH (Brief="Outlet Hot Stream", PosX=1, PosY=0.25, Symbol="^{outHot}"); 110 out OutletCold as streamPH (Brief="Outlet Cold Stream", PosX=1, PosY=0.75, Symbol="^{outCold}"); 111 112 HotSide as Main_PHE (Brief="Plate Heat Exchanger Hot Side", Symbol="_{hot}"); 113 ColdSide as Main_PHE (Brief="Plate Heat Exchanger Cold Side", Symbol="_{cold}"); 114 Thermal as Thermal_PHE (Brief="Thermal Results", Symbol = " "); 115 116 SET 117 #"Vector Length of constants for Kumar's calculating Pressure Drop" 118 N1 = 15; 119 120 #"Vector Length of constants for Kumar's calculating Heat Transfer" 121 N2 = 14; 122 123 #"First constant for Kumar's calculating Pressure Drop" 124 Kp1 = [50,19.40,2.990,47,18.290,1.441,34,11.250,0.772,24,3.240,0.760,24,2.80,0.639]; 125 126 #"Second constant for Kumar's calculating Pressure Drop" 127 Kp2 = [1,0.589,0.183,1,0.652,0.206,1,0.631,0.161,1,0.457,0.215,1,0.451,0.213]; 128 129 #"First constant for Kumar's calculating Heat Transfer" 130 Kc1 = [0.718,0.348,0.718,0.400,0.300,0.630,0.291,0.130,0.562,0.306,0.108,0.562,0.331,0.087]; 131 132 #"Second constant for Kumar's calculating Heat Transfer" 133 Kc2 = [0.349,0.663,0.349,0.598,0.663,0.333,0.591,0.732,0.326,0.529,0.703,0.326,0.503,0.718]; 134 135 #"Component Molecular Weight" 136 M = PP.MolecularWeight(); 137 138 #"Pi Number" 139 Pi = 3.14159265; 140 141 #"Plate Vertical Distance between Port Centers" 142 Lp = Lv - Dports; 143 144 #"Corrugation Depth" 145 Depth=pitch-pt; 146 147 #"Plate Horizontal Distance between Port Centers" 148 Lh=Lw-Dports; 149 150 #"Hydraulic Diameter" 151 Dh=2*Depth/PhiFactor; 152 153 #"Ports Area" 154 Aports=0.25*Pi*Dports*Dports; 155 156 #"Channel Area" 157 Achannel=Depth*Lw; 158 159 #"Pack Length" 160 Lpack=Depth*(Nplates-1)+Nplates*pt; 161 162 #"Total Number of Channels" 163 Nchannels = Nplates -1; 164 165 #"Exchange Surface Area" 166 Atotal =(Nplates-2)*Lw*Lp*PhiFactor; 167 254 Geometry as PHE_Geometry (Brief="Plate Heat Exchanger Geometrical Parameters", Symbol=" "); 255 in InletHot as stream (Brief="Inlet Hot Stream", PosX=0, PosY=0.75, Symbol="^{inHot}"); 256 in InletCold as stream (Brief="Inlet Cold Stream", PosX=0, PosY=0.25, Symbol="^{inCold}"); 257 out OutletHot as streamPH (Brief="Outlet Hot Stream", PosX=1, PosY=0.25, Symbol="^{outHot}"); 258 out OutletCold as streamPH (Brief="Outlet Cold Stream", PosX=1, PosY=0.75, Symbol="^{outCold}"); 259 260 261 HotSide as Main_PHE (Brief="Plate Heat Exchanger Hot Side", Symbol="_{hot}"); 262 ColdSide as Main_PHE (Brief="Plate Heat Exchanger Cold Side", Symbol="_{cold}"); 263 Thermal as Thermal_PHE (Brief="Thermal Results", Symbol = " "); 264 168 265 EQUATIONS 169 266 … … 187 284 188 285 "Hot Stream Average Molecular Weight" 189 HotSide.Properties.Average.Mw = sum( M*InletHot.z);286 HotSide.Properties.Average.Mw = sum(Geometry.M*InletHot.z); 190 287 191 288 "Cold Stream Average Molecular Weight" 192 ColdSide.Properties.Average.Mw = sum( M*InletCold.z);289 ColdSide.Properties.Average.Mw = sum(Geometry.M*InletCold.z); 193 290 194 291 if InletCold.v equal 0 … … 300 397 301 398 "Flow Mass Inlet Cold Stream" 302 ColdSide.Properties.Inlet.Fw = sum( M*InletCold.z)*InletCold.F;399 ColdSide.Properties.Inlet.Fw = sum(Geometry.M*InletCold.z)*InletCold.F; 303 400 304 401 "Flow Mass Outlet Cold Stream" 305 ColdSide.Properties.Outlet.Fw = sum( M*OutletCold.z)*OutletCold.F;402 ColdSide.Properties.Outlet.Fw = sum(Geometry.M*OutletCold.z)*OutletCold.F; 306 403 307 404 "Flow Mass Inlet Hot Stream" 308 HotSide.Properties.Inlet.Fw = sum( M*InletHot.z)*InletHot.F;405 HotSide.Properties.Inlet.Fw = sum(Geometry.M*InletHot.z)*InletHot.F; 309 406 310 407 "Flow Mass Outlet Hot Stream" 311 HotSide.Properties.Outlet.Fw = sum( M*OutletHot.z)*OutletHot.F;408 HotSide.Properties.Outlet.Fw = sum(Geometry.M*OutletHot.z)*OutletHot.F; 312 409 313 410 "Molar Balance Hot Stream" … … 328 425 329 426 "Total Number of Passages Cold Side" 330 ColdSide.PressureDrop.Npassage = (2* Nchannels+1+(-1)^(Nchannels+1))/(4*NpassCold);427 ColdSide.PressureDrop.Npassage = (2*Geometry.Nchannels+1+(-1)^(Geometry.Nchannels+1))/(4*Geometry.NpassCold); 331 428 332 429 "Total Number of Passages Hot Side" 333 HotSide.PressureDrop.Npassage = (2* Nchannels-1+(-1)^(Nchannels))/(4*NpassHot);430 HotSide.PressureDrop.Npassage = (2*Geometry.Nchannels-1+(-1)^(Geometry.Nchannels))/(4*Geometry.NpassHot); 334 431 335 432 case "hot": 336 433 337 434 "Total Number of Passages Cold Side" 338 HotSide.PressureDrop.Npassage = (2* Nchannels+1+(-1)^(Nchannels+1))/(4*NpassHot);435 HotSide.PressureDrop.Npassage = (2*Geometry.Nchannels+1+(-1)^(Geometry.Nchannels+1))/(4*Geometry.NpassHot); 339 436 340 437 "Total Number of Passages Hot Side" 341 ColdSide.PressureDrop.Npassage = (2* Nchannels-1+(-1)^(Nchannels))/(4*NpassCold);438 ColdSide.PressureDrop.Npassage = (2*Geometry.Nchannels-1+(-1)^(Geometry.Nchannels))/(4*Geometry.NpassCold); 342 439 343 440 end 344 441 345 442 "Hot Stream Mass Flux in the Channel" 346 HotSide.HeatTransfer.Gchannel=HotSide.Properties.Inlet.Fw/(HotSide.PressureDrop.Npassage* Achannel);443 HotSide.HeatTransfer.Gchannel=HotSide.Properties.Inlet.Fw/(HotSide.PressureDrop.Npassage*Geometry.Achannel); 347 444 348 445 "Hot Stream Mass Flux in the Ports" 349 HotSide.HeatTransfer.Gports=HotSide.Properties.Inlet.Fw/ Aports;446 HotSide.HeatTransfer.Gports=HotSide.Properties.Inlet.Fw/Geometry.Aports; 350 447 351 448 "Cold Stream Mass Flux in the Ports" 352 ColdSide.HeatTransfer.Gports=ColdSide.Properties.Inlet.Fw/ Aports;449 ColdSide.HeatTransfer.Gports=ColdSide.Properties.Inlet.Fw/Geometry.Aports; 353 450 354 451 "Cold Stream Mass Flux in the Channel" 355 ColdSide.HeatTransfer.Gchannel=ColdSide.Properties.Inlet.Fw/(ColdSide.PressureDrop.Npassage* Achannel);452 ColdSide.HeatTransfer.Gchannel=ColdSide.Properties.Inlet.Fw/(ColdSide.PressureDrop.Npassage*Geometry.Achannel); 356 453 357 454 "Hot Stream Pressure Drop in Ports" 358 HotSide.PressureDrop.DPports =1.5* NpassHot*HotSide.HeatTransfer.Gports^2/(2*HotSide.Properties.Average.rho);455 HotSide.PressureDrop.DPports =1.5*Geometry.NpassHot*HotSide.HeatTransfer.Gports^2/(2*HotSide.Properties.Average.rho); 359 456 360 457 "Cold Stream Pressure Drop in Ports" 361 ColdSide.PressureDrop.DPports =1.5* NpassCold*ColdSide.HeatTransfer.Gports^2/(2*ColdSide.Properties.Average.rho);458 ColdSide.PressureDrop.DPports =1.5*Geometry.NpassCold*ColdSide.HeatTransfer.Gports^2/(2*ColdSide.Properties.Average.rho); 362 459 363 460 "Hot Stream Pressure Drop in Channels" 364 HotSide.PressureDrop.DPchannel =2*HotSide.PressureDrop.fi* NpassHot*Lv*HotSide.HeatTransfer.Gchannel^2/(HotSide.Properties.Average.rho*Dh*HotSide.HeatTransfer.Phi^0.17);461 HotSide.PressureDrop.DPchannel =2*HotSide.PressureDrop.fi*Geometry.NpassHot*Geometry.Lv*HotSide.HeatTransfer.Gchannel^2/(HotSide.Properties.Average.rho*Geometry.Dh*HotSide.HeatTransfer.Phi^0.17); 365 462 366 463 "Cold Stream Pressure Drop in Channels" 367 ColdSide.PressureDrop.DPchannel =2*ColdSide.PressureDrop.fi* NpassCold*Lv*ColdSide.HeatTransfer.Gchannel^2/(ColdSide.Properties.Average.rho*Dh*ColdSide.HeatTransfer.Phi^0.17);464 ColdSide.PressureDrop.DPchannel =2*ColdSide.PressureDrop.fi*Geometry.NpassCold*Geometry.Lv*ColdSide.HeatTransfer.Gchannel^2/(ColdSide.Properties.Average.rho*Geometry.Dh*ColdSide.HeatTransfer.Phi^0.17); 368 465 369 466 "Hot Stream Total Pressure Drop" … … 379 476 if HotSide.HeatTransfer.Re < 10 380 477 then 381 HotSide.PressureDrop.fi = Kp1(1)/HotSide.HeatTransfer.Re^Kp2(1);382 ColdSide.PressureDrop.fi = Kp1(1)/ColdSide.HeatTransfer.Re^Kp2(1);478 HotSide.PressureDrop.fi = Geometry.Kp1(1)/HotSide.HeatTransfer.Re^Geometry.Kp2(1); 479 ColdSide.PressureDrop.fi = Geometry.Kp1(1)/ColdSide.HeatTransfer.Re^Geometry.Kp2(1); 383 480 else 384 481 if HotSide.HeatTransfer.Re < 100 385 482 then 386 HotSide.PressureDrop.fi = Kp1(2)/HotSide.HeatTransfer.Re^Kp2(2);387 ColdSide.PressureDrop.fi = Kp1(2)/ColdSide.HeatTransfer.Re^Kp2(2);388 else 389 HotSide.PressureDrop.fi = Kp1(3)/HotSide.HeatTransfer.Re^Kp2(3);390 ColdSide.PressureDrop.fi = Kp1(3)/ColdSide.HeatTransfer.Re^Kp2(3);483 HotSide.PressureDrop.fi = Geometry.Kp1(2)/HotSide.HeatTransfer.Re^Geometry.Kp2(2); 484 ColdSide.PressureDrop.fi = Geometry.Kp1(2)/ColdSide.HeatTransfer.Re^Geometry.Kp2(2); 485 else 486 HotSide.PressureDrop.fi = Geometry.Kp1(3)/HotSide.HeatTransfer.Re^Geometry.Kp2(3); 487 ColdSide.PressureDrop.fi = Geometry.Kp1(3)/ColdSide.HeatTransfer.Re^Geometry.Kp2(3); 391 488 end 392 489 … … 397 494 if HotSide.HeatTransfer.Re < 15 398 495 then 399 HotSide.PressureDrop.fi = Kp1(4)/HotSide.HeatTransfer.Re^Kp2(4);400 ColdSide.PressureDrop.fi = Kp1(4)/ColdSide.HeatTransfer.Re^Kp2(4);496 HotSide.PressureDrop.fi = Geometry.Kp1(4)/HotSide.HeatTransfer.Re^Geometry.Kp2(4); 497 ColdSide.PressureDrop.fi = Geometry.Kp1(4)/ColdSide.HeatTransfer.Re^Geometry.Kp2(4); 401 498 else 402 499 if HotSide.HeatTransfer.Re < 300 403 500 then 404 HotSide.PressureDrop.fi = Kp1(5)/HotSide.HeatTransfer.Re^Kp2(5);405 ColdSide.PressureDrop.fi = Kp1(5)/ColdSide.HeatTransfer.Re^Kp2(5);406 else 407 HotSide.PressureDrop.fi = Kp1(6)/HotSide.HeatTransfer.Re^Kp2(6);408 ColdSide.PressureDrop.fi = Kp1(6)/ColdSide.HeatTransfer.Re^Kp2(6);501 HotSide.PressureDrop.fi = Geometry.Kp1(5)/HotSide.HeatTransfer.Re^Geometry.Kp2(5); 502 ColdSide.PressureDrop.fi = Geometry.Kp1(5)/ColdSide.HeatTransfer.Re^Geometry.Kp2(5); 503 else 504 HotSide.PressureDrop.fi = Geometry.Kp1(6)/HotSide.HeatTransfer.Re^Geometry.Kp2(6); 505 ColdSide.PressureDrop.fi = Geometry.Kp1(6)/ColdSide.HeatTransfer.Re^Geometry.Kp2(6); 409 506 end 410 507 … … 415 512 if HotSide.HeatTransfer.Re < 20 416 513 then 417 HotSide.PressureDrop.fi = Kp1(7)/HotSide.HeatTransfer.Re^Kp2(7);418 ColdSide.PressureDrop.fi = Kp1(7)/ColdSide.HeatTransfer.Re^Kp2(7);514 HotSide.PressureDrop.fi = Geometry.Kp1(7)/HotSide.HeatTransfer.Re^Geometry.Kp2(7); 515 ColdSide.PressureDrop.fi = Geometry.Kp1(7)/ColdSide.HeatTransfer.Re^Geometry.Kp2(7); 419 516 else 420 517 if HotSide.HeatTransfer.Re < 300 421 518 then 422 HotSide.PressureDrop.fi = Kp1(8)/HotSide.HeatTransfer.Re^Kp2(8);423 ColdSide.PressureDrop.fi = Kp1(8)/ColdSide.HeatTransfer.Re^Kp2(8);424 else 425 HotSide.PressureDrop.fi = Kp1(9)/HotSide.HeatTransfer.Re^Kp2(9);426 ColdSide.PressureDrop.fi = Kp1(9)/ColdSide.HeatTransfer.Re^Kp2(9);519 HotSide.PressureDrop.fi = Geometry.Kp1(8)/HotSide.HeatTransfer.Re^Geometry.Kp2(8); 520 ColdSide.PressureDrop.fi = Geometry.Kp1(8)/ColdSide.HeatTransfer.Re^Geometry.Kp2(8); 521 else 522 HotSide.PressureDrop.fi = Geometry.Kp1(9)/HotSide.HeatTransfer.Re^Geometry.Kp2(9); 523 ColdSide.PressureDrop.fi = Geometry.Kp1(9)/ColdSide.HeatTransfer.Re^Geometry.Kp2(9); 427 524 end 428 525 … … 433 530 if HotSide.HeatTransfer.Re < 40 434 531 then 435 HotSide.PressureDrop.fi = Kp1(10)/HotSide.HeatTransfer.Re^Kp2(10);436 ColdSide.PressureDrop.fi = Kp1(10)/ColdSide.HeatTransfer.Re^Kp2(10);532 HotSide.PressureDrop.fi = Geometry.Kp1(10)/HotSide.HeatTransfer.Re^Geometry.Kp2(10); 533 ColdSide.PressureDrop.fi = Geometry.Kp1(10)/ColdSide.HeatTransfer.Re^Geometry.Kp2(10); 437 534 else 438 535 if HotSide.HeatTransfer.Re < 400 439 536 then 440 HotSide.PressureDrop.fi = Kp1(11)/HotSide.HeatTransfer.Re^Kp2(11);441 ColdSide.PressureDrop.fi = Kp1(11)/ColdSide.HeatTransfer.Re^Kp2(11);442 else 443 HotSide.PressureDrop.fi = Kp1(12)/HotSide.HeatTransfer.Re^Kp2(12);444 ColdSide.PressureDrop.fi = Kp1(12)/ColdSide.HeatTransfer.Re^Kp2(12);537 HotSide.PressureDrop.fi = Geometry.Kp1(11)/HotSide.HeatTransfer.Re^Geometry.Kp2(11); 538 ColdSide.PressureDrop.fi = Geometry.Kp1(11)/ColdSide.HeatTransfer.Re^Geometry.Kp2(11); 539 else 540 HotSide.PressureDrop.fi = Geometry.Kp1(12)/HotSide.HeatTransfer.Re^Geometry.Kp2(12); 541 ColdSide.PressureDrop.fi = Geometry.Kp1(12)/ColdSide.HeatTransfer.Re^Geometry.Kp2(12); 445 542 end 446 543 … … 451 548 if HotSide.HeatTransfer.Re < 50 452 549 then 453 HotSide.PressureDrop.fi = Kp1(13)/HotSide.HeatTransfer.Re^Kp2(13);454 ColdSide.PressureDrop.fi = Kp1(13)/ColdSide.HeatTransfer.Re^Kp2(13);550 HotSide.PressureDrop.fi = Geometry.Kp1(13)/HotSide.HeatTransfer.Re^Geometry.Kp2(13); 551 ColdSide.PressureDrop.fi = Geometry.Kp1(13)/ColdSide.HeatTransfer.Re^Geometry.Kp2(13); 455 552 else 456 553 if HotSide.HeatTransfer.Re < 500 457 554 then 458 HotSide.PressureDrop.fi = Kp1(14)/HotSide.HeatTransfer.Re^Kp2(14);459 ColdSide.PressureDrop.fi = Kp1(14)/ColdSide.HeatTransfer.Re^Kp2(14);460 else 461 HotSide.PressureDrop.fi = Kp1(15)/HotSide.HeatTransfer.Re^Kp2(15);462 ColdSide.PressureDrop.fi = Kp1(15)/ColdSide.HeatTransfer.Re^Kp2(15);555 HotSide.PressureDrop.fi = Geometry.Kp1(14)/HotSide.HeatTransfer.Re^Geometry.Kp2(14); 556 ColdSide.PressureDrop.fi = Geometry.Kp1(14)/ColdSide.HeatTransfer.Re^Geometry.Kp2(14); 557 else 558 HotSide.PressureDrop.fi = Geometry.Kp1(15)/HotSide.HeatTransfer.Re^Geometry.Kp2(15); 559 ColdSide.PressureDrop.fi = Geometry.Kp1(15)/ColdSide.HeatTransfer.Re^Geometry.Kp2(15); 463 560 end 464 561 … … 473 570 if HotSide.HeatTransfer.Re < 10 474 571 then 475 HotSide.HeatTransfer.hcoeff=(HotSide.Properties.Average.K*HotSide.HeatTransfer.PR^(1/3)*HotSide.HeatTransfer.Phi^0.17* Kc1(1)*HotSide.HeatTransfer.Re^Kc2(1))/Dh;476 ColdSide.HeatTransfer.hcoeff =(ColdSide.Properties.Average.K*ColdSide.HeatTransfer.PR^(1/3)*ColdSide.HeatTransfer.Phi^0.17* Kc1(1)*ColdSide.HeatTransfer.Re^Kc2(1))/Dh;477 else 478 HotSide.HeatTransfer.hcoeff=(HotSide.Properties.Average.K*HotSide.HeatTransfer.PR^(1/3)*HotSide.HeatTransfer.Phi^0.17* Kc1(2)*HotSide.HeatTransfer.Re^Kc2(2))/Dh;479 ColdSide.HeatTransfer.hcoeff =(ColdSide.Properties.Average.K*ColdSide.HeatTransfer.PR^(1/3)*ColdSide.HeatTransfer.Phi^0.17* Kc1(2)*ColdSide.HeatTransfer.Re^Kc2(2))/Dh;572 HotSide.HeatTransfer.hcoeff=(HotSide.Properties.Average.K*HotSide.HeatTransfer.PR^(1/3)*HotSide.HeatTransfer.Phi^0.17*Geometry.Kc1(1)*HotSide.HeatTransfer.Re^Geometry.Kc2(1))/Geometry.Dh; 573 ColdSide.HeatTransfer.hcoeff =(ColdSide.Properties.Average.K*ColdSide.HeatTransfer.PR^(1/3)*ColdSide.HeatTransfer.Phi^0.17*Geometry.Kc1(1)*ColdSide.HeatTransfer.Re^Geometry.Kc2(1))/Geometry.Dh; 574 else 575 HotSide.HeatTransfer.hcoeff=(HotSide.Properties.Average.K*HotSide.HeatTransfer.PR^(1/3)*HotSide.HeatTransfer.Phi^0.17*Geometry.Kc1(2)*HotSide.HeatTransfer.Re^Geometry.Kc2(2))/Geometry.Dh; 576 ColdSide.HeatTransfer.hcoeff =(ColdSide.Properties.Average.K*ColdSide.HeatTransfer.PR^(1/3)*ColdSide.HeatTransfer.Phi^0.17*Geometry.Kc1(2)*ColdSide.HeatTransfer.Re^Geometry.Kc2(2))/Geometry.Dh; 480 577 end 481 578 … … 484 581 if HotSide.HeatTransfer.Re < 10 485 582 then 486 HotSide.HeatTransfer.hcoeff=(HotSide.Properties.Average.K*HotSide.HeatTransfer.PR^(1/3)*HotSide.HeatTransfer.Phi^0.17* Kc1(3)*HotSide.HeatTransfer.Re^Kc2(3))/Dh;487 ColdSide.HeatTransfer.hcoeff =(ColdSide.Properties.Average.K*ColdSide.HeatTransfer.PR^(1/3)*ColdSide.HeatTransfer.Phi^0.17* Kc1(3)*ColdSide.HeatTransfer.Re^Kc2(3))/Dh;583 HotSide.HeatTransfer.hcoeff=(HotSide.Properties.Average.K*HotSide.HeatTransfer.PR^(1/3)*HotSide.HeatTransfer.Phi^0.17*Geometry.Kc1(3)*HotSide.HeatTransfer.Re^Geometry.Kc2(3))/Geometry.Dh; 584 ColdSide.HeatTransfer.hcoeff =(ColdSide.Properties.Average.K*ColdSide.HeatTransfer.PR^(1/3)*ColdSide.HeatTransfer.Phi^0.17*Geometry.Kc1(3)*ColdSide.HeatTransfer.Re^Geometry.Kc2(3))/Geometry.Dh; 488 585 else 489 586 if HotSide.HeatTransfer.Re < 100 490 587 then 491 HotSide.HeatTransfer.hcoeff=(HotSide.Properties.Average.K*HotSide.HeatTransfer.PR^(1/3)*HotSide.HeatTransfer.Phi^0.17* Kc1(4)*HotSide.HeatTransfer.Re^Kc2(4))/Dh;492 ColdSide.HeatTransfer.hcoeff =(ColdSide.Properties.Average.K*ColdSide.HeatTransfer.PR^(1/3)*ColdSide.HeatTransfer.Phi^0.17* Kc1(4)*ColdSide.HeatTransfer.Re^Kc2(4))/Dh;493 else 494 HotSide.HeatTransfer.hcoeff=(HotSide.Properties.Average.K*HotSide.HeatTransfer.PR^(1/3)*HotSide.HeatTransfer.Phi^0.17* Kc1(5)*HotSide.HeatTransfer.Re^Kc2(5))/Dh;495 ColdSide.HeatTransfer.hcoeff =(ColdSide.Properties.Average.K*ColdSide.HeatTransfer.PR^(1/3)*ColdSide.HeatTransfer.Phi^0.17* Kc1(5)*ColdSide.HeatTransfer.Re^Kc2(5))/Dh;588 HotSide.HeatTransfer.hcoeff=(HotSide.Properties.Average.K*HotSide.HeatTransfer.PR^(1/3)*HotSide.HeatTransfer.Phi^0.17*Geometry.Kc1(4)*HotSide.HeatTransfer.Re^Geometry.Kc2(4))/Geometry.Dh; 589 ColdSide.HeatTransfer.hcoeff =(ColdSide.Properties.Average.K*ColdSide.HeatTransfer.PR^(1/3)*ColdSide.HeatTransfer.Phi^0.17*Geometry.Kc1(4)*ColdSide.HeatTransfer.Re^Geometry.Kc2(4))/Geometry.Dh; 590 else 591 HotSide.HeatTransfer.hcoeff=(HotSide.Properties.Average.K*HotSide.HeatTransfer.PR^(1/3)*HotSide.HeatTransfer.Phi^0.17*Geometry.Kc1(5)*HotSide.HeatTransfer.Re^Geometry.Kc2(5))/Geometry.Dh; 592 ColdSide.HeatTransfer.hcoeff =(ColdSide.Properties.Average.K*ColdSide.HeatTransfer.PR^(1/3)*ColdSide.HeatTransfer.Phi^0.17*Geometry.Kc1(5)*ColdSide.HeatTransfer.Re^Geometry.Kc2(5))/Geometry.Dh; 496 593 end 497 594 end … … 501 598 if HotSide.HeatTransfer.Re < 20 502 599 then 503 HotSide.HeatTransfer.hcoeff=(HotSide.Properties.Average.K*HotSide.HeatTransfer.PR^(1/3)*HotSide.HeatTransfer.Phi^0.17* Kc1(6)*HotSide.HeatTransfer.Re^Kc2(6))/Dh;504 ColdSide.HeatTransfer.hcoeff =(ColdSide.Properties.Average.K*ColdSide.HeatTransfer.PR^(1/3)*ColdSide.HeatTransfer.Phi^0.17* Kc1(6)*ColdSide.HeatTransfer.Re^Kc2(6))/Dh;600 HotSide.HeatTransfer.hcoeff=(HotSide.Properties.Average.K*HotSide.HeatTransfer.PR^(1/3)*HotSide.HeatTransfer.Phi^0.17*Geometry.Kc1(6)*HotSide.HeatTransfer.Re^Geometry.Kc2(6))/Geometry.Dh; 601 ColdSide.HeatTransfer.hcoeff =(ColdSide.Properties.Average.K*ColdSide.HeatTransfer.PR^(1/3)*ColdSide.HeatTransfer.Phi^0.17*Geometry.Kc1(6)*ColdSide.HeatTransfer.Re^Geometry.Kc2(6))/Geometry.Dh; 505 602 else 506 603 if HotSide.HeatTransfer.Re < 300 507 604 then 508 HotSide.HeatTransfer.hcoeff=(HotSide.Properties.Average.K*HotSide.HeatTransfer.PR^(1/3)*HotSide.HeatTransfer.Phi^0.17* Kc1(7)*HotSide.HeatTransfer.Re^Kc2(7))/Dh;509 ColdSide.HeatTransfer.hcoeff =(ColdSide.Properties.Average.K*ColdSide.HeatTransfer.PR^(1/3)*ColdSide.HeatTransfer.Phi^0.17* Kc1(7)*ColdSide.HeatTransfer.Re^Kc2(7))/Dh;510 else 511 HotSide.HeatTransfer.hcoeff=(HotSide.Properties.Average.K*HotSide.HeatTransfer.PR^(1/3)*HotSide.HeatTransfer.Phi^0.17* Kc1(8)*HotSide.HeatTransfer.Re^Kc2(8))/Dh;512 ColdSide.HeatTransfer.hcoeff =(ColdSide.Properties.Average.K*ColdSide.HeatTransfer.PR^(1/3)*ColdSide.HeatTransfer.Phi^0.17* Kc1(8)*ColdSide.HeatTransfer.Re^Kc2(8))/Dh;605 HotSide.HeatTransfer.hcoeff=(HotSide.Properties.Average.K*HotSide.HeatTransfer.PR^(1/3)*HotSide.HeatTransfer.Phi^0.17*Geometry.Kc1(7)*HotSide.HeatTransfer.Re^Geometry.Kc2(7))/Geometry.Dh; 606 ColdSide.HeatTransfer.hcoeff =(ColdSide.Properties.Average.K*ColdSide.HeatTransfer.PR^(1/3)*ColdSide.HeatTransfer.Phi^0.17*Geometry.Kc1(7)*ColdSide.HeatTransfer.Re^Geometry.Kc2(7))/Geometry.Dh; 607 else 608 HotSide.HeatTransfer.hcoeff=(HotSide.Properties.Average.K*HotSide.HeatTransfer.PR^(1/3)*HotSide.HeatTransfer.Phi^0.17*Geometry.Kc1(8)*HotSide.HeatTransfer.Re^Geometry.Kc2(8))/Geometry.Dh; 609 ColdSide.HeatTransfer.hcoeff =(ColdSide.Properties.Average.K*ColdSide.HeatTransfer.PR^(1/3)*ColdSide.HeatTransfer.Phi^0.17*Geometry.Kc1(8)*ColdSide.HeatTransfer.Re^Geometry.Kc2(8))/Geometry.Dh; 513 610 end 514 611 end … … 518 615 if HotSide.HeatTransfer.Re < 20 519 616 then 520 HotSide.HeatTransfer.hcoeff=(HotSide.Properties.Average.K*HotSide.HeatTransfer.PR^(1/3)*HotSide.HeatTransfer.Phi^0.17* Kc1(9)*HotSide.HeatTransfer.Re^Kc2(9))/Dh;521 ColdSide.HeatTransfer.hcoeff =(ColdSide.Properties.Average.K*ColdSide.HeatTransfer.PR^(1/3)*ColdSide.HeatTransfer.Phi^0.17* Kc1(9)*ColdSide.HeatTransfer.Re^Kc2(9))/Dh;617 HotSide.HeatTransfer.hcoeff=(HotSide.Properties.Average.K*HotSide.HeatTransfer.PR^(1/3)*HotSide.HeatTransfer.Phi^0.17*Geometry.Kc1(9)*HotSide.HeatTransfer.Re^Geometry.Kc2(9))/Geometry.Dh; 618 ColdSide.HeatTransfer.hcoeff =(ColdSide.Properties.Average.K*ColdSide.HeatTransfer.PR^(1/3)*ColdSide.HeatTransfer.Phi^0.17*Geometry.Kc1(9)*ColdSide.HeatTransfer.Re^Geometry.Kc2(9))/Geometry.Dh; 522 619 else 523 620 if HotSide.HeatTransfer.Re < 400 524 621 then 525 HotSide.HeatTransfer.hcoeff=(HotSide.Properties.Average.K*HotSide.HeatTransfer.PR^(1/3)*HotSide.HeatTransfer.Phi^0.17* Kc1(10)*HotSide.HeatTransfer.Re^Kc2(10))/Dh;526 ColdSide.HeatTransfer.hcoeff =(ColdSide.Properties.Average.K*ColdSide.HeatTransfer.PR^(1/3)*ColdSide.HeatTransfer.Phi^0.17* Kc1(10)*ColdSide.HeatTransfer.Re^Kc2(10))/Dh;527 else 528 HotSide.HeatTransfer.hcoeff=(HotSide.Properties.Average.K*HotSide.HeatTransfer.PR^(1/3)*HotSide.HeatTransfer.Phi^0.17* Kc1(11)*HotSide.HeatTransfer.Re^Kc2(11))/Dh;529 ColdSide.HeatTransfer.hcoeff =(ColdSide.Properties.Average.K*ColdSide.HeatTransfer.PR^(1/3)*ColdSide.HeatTransfer.Phi^0.17* Kc1(11)*ColdSide.HeatTransfer.Re^Kc2(11))/Dh;622 HotSide.HeatTransfer.hcoeff=(HotSide.Properties.Average.K*HotSide.HeatTransfer.PR^(1/3)*HotSide.HeatTransfer.Phi^0.17*Geometry.Kc1(10)*HotSide.HeatTransfer.Re^Geometry.Kc2(10))/Geometry.Dh; 623 ColdSide.HeatTransfer.hcoeff =(ColdSide.Properties.Average.K*ColdSide.HeatTransfer.PR^(1/3)*ColdSide.HeatTransfer.Phi^0.17*Geometry.Kc1(10)*ColdSide.HeatTransfer.Re^Geometry.Kc2(10))/Geometry.Dh; 624 else 625 HotSide.HeatTransfer.hcoeff=(HotSide.Properties.Average.K*HotSide.HeatTransfer.PR^(1/3)*HotSide.HeatTransfer.Phi^0.17*Geometry.Kc1(11)*HotSide.HeatTransfer.Re^Geometry.Kc2(11))/Geometry.Dh; 626 ColdSide.HeatTransfer.hcoeff =(ColdSide.Properties.Average.K*ColdSide.HeatTransfer.PR^(1/3)*ColdSide.HeatTransfer.Phi^0.17*Geometry.Kc1(11)*ColdSide.HeatTransfer.Re^Geometry.Kc2(11))/Geometry.Dh; 530 627 end 531 628 end … … 535 632 if HotSide.HeatTransfer.Re < 20 536 633 then 537 HotSide.HeatTransfer.hcoeff=(HotSide.Properties.Average.K*HotSide.HeatTransfer.PR^(1/3)*HotSide.HeatTransfer.Phi^0.17* Kc1(12)*HotSide.HeatTransfer.Re^Kc2(12))/Dh;538 ColdSide.HeatTransfer.hcoeff =(ColdSide.Properties.Average.K*ColdSide.HeatTransfer.PR^(1/3)*ColdSide.HeatTransfer.Phi^0.17* Kc1(12)*ColdSide.HeatTransfer.Re^Kc2(12))/Dh;634 HotSide.HeatTransfer.hcoeff=(HotSide.Properties.Average.K*HotSide.HeatTransfer.PR^(1/3)*HotSide.HeatTransfer.Phi^0.17*Geometry.Kc1(12)*HotSide.HeatTransfer.Re^Geometry.Kc2(12))/Geometry.Dh; 635 ColdSide.HeatTransfer.hcoeff =(ColdSide.Properties.Average.K*ColdSide.HeatTransfer.PR^(1/3)*ColdSide.HeatTransfer.Phi^0.17*Geometry.Kc1(12)*ColdSide.HeatTransfer.Re^Geometry.Kc2(12))/Geometry.Dh; 539 636 else 540 637 if HotSide.HeatTransfer.Re < 500 541 638 then 542 HotSide.HeatTransfer.hcoeff=(HotSide.Properties.Average.K*HotSide.HeatTransfer.PR^(1/3)*HotSide.HeatTransfer.Phi^0.17* Kc1(13)*HotSide.HeatTransfer.Re^Kc2(13))/Dh;543 ColdSide.HeatTransfer.hcoeff =(ColdSide.Properties.Average.K*ColdSide.HeatTransfer.PR^(1/3)*ColdSide.HeatTransfer.Phi^0.17* Kc1(13)*ColdSide.HeatTransfer.Re^Kc2(13))/Dh;544 else 545 HotSide.HeatTransfer.hcoeff=(HotSide.Properties.Average.K*HotSide.HeatTransfer.PR^(1/3)*HotSide.HeatTransfer.Phi^0.17* Kc1(14)*HotSide.HeatTransfer.Re^Kc2(14))/Dh;546 ColdSide.HeatTransfer.hcoeff =(ColdSide.Properties.Average.K*ColdSide.HeatTransfer.PR^(1/3)*ColdSide.HeatTransfer.Phi^0.17* Kc1(14)*ColdSide.HeatTransfer.Re^Kc2(14))/Dh;639 HotSide.HeatTransfer.hcoeff=(HotSide.Properties.Average.K*HotSide.HeatTransfer.PR^(1/3)*HotSide.HeatTransfer.Phi^0.17*Geometry.Kc1(13)*HotSide.HeatTransfer.Re^Geometry.Kc2(13))/Geometry.Dh; 640 ColdSide.HeatTransfer.hcoeff =(ColdSide.Properties.Average.K*ColdSide.HeatTransfer.PR^(1/3)*ColdSide.HeatTransfer.Phi^0.17*Geometry.Kc1(13)*ColdSide.HeatTransfer.Re^Geometry.Kc2(13))/Geometry.Dh; 641 else 642 HotSide.HeatTransfer.hcoeff=(HotSide.Properties.Average.K*HotSide.HeatTransfer.PR^(1/3)*HotSide.HeatTransfer.Phi^0.17*Geometry.Kc1(14)*HotSide.HeatTransfer.Re^Geometry.Kc2(14))/Geometry.Dh; 643 ColdSide.HeatTransfer.hcoeff =(ColdSide.Properties.Average.K*ColdSide.HeatTransfer.PR^(1/3)*ColdSide.HeatTransfer.Phi^0.17*Geometry.Kc1(14)*ColdSide.HeatTransfer.Re^Geometry.Kc2(14))/Geometry.Dh; 547 644 end 548 645 end … … 557 654 558 655 "Hot Stream Velocity in Ports" 559 HotSide.PressureDrop.Vports =HotSide.Properties.Inlet.Fw/( Aports*HotSide.Properties.Inlet.rho);656 HotSide.PressureDrop.Vports =HotSide.Properties.Inlet.Fw/(Geometry.Aports*HotSide.Properties.Inlet.rho); 560 657 561 658 "Cold Stream Velocity in Ports" 562 ColdSide.PressureDrop.Vports =ColdSide.Properties.Inlet.Fw/( Aports*ColdSide.Properties.Inlet.rho);659 ColdSide.PressureDrop.Vports =ColdSide.Properties.Inlet.Fw/(Geometry.Aports*ColdSide.Properties.Inlet.rho); 563 660 564 661 "Hot Stream Reynolds Number" 565 HotSide.HeatTransfer.Re = Dh*HotSide.HeatTransfer.Gchannel/HotSide.Properties.Average.Mu;662 HotSide.HeatTransfer.Re =Geometry.Dh*HotSide.HeatTransfer.Gchannel/HotSide.Properties.Average.Mu; 566 663 567 664 "Cold Stream Reynolds Number" 568 ColdSide.HeatTransfer.Re = Dh*ColdSide.HeatTransfer.Gchannel/ColdSide.Properties.Average.Mu;665 ColdSide.HeatTransfer.Re =Geometry.Dh*ColdSide.HeatTransfer.Gchannel/ColdSide.Properties.Average.Mu; 569 666 570 667 "Hot Stream Prandtl Number" … … 587 684 588 685 "Overall Heat Transfer Coefficient Clean" 589 Thermal.Uc/HotSide.HeatTransfer.hcoeff +Thermal.Uc* pt/Kwall+Thermal.Uc/ColdSide.HeatTransfer.hcoeff=1;686 Thermal.Uc/HotSide.HeatTransfer.hcoeff +Thermal.Uc*Geometry.pt/Geometry.Kwall+Thermal.Uc/ColdSide.HeatTransfer.hcoeff=1; 590 687 591 688 "Overall Heat Transfer Coefficient Dirty" 592 Thermal.Ud*(1/HotSide.HeatTransfer.hcoeff + pt/Kwall+1/ColdSide.HeatTransfer.hcoeff + Rfc +Rfh)=1;689 Thermal.Ud*(1/HotSide.HeatTransfer.hcoeff +Geometry.pt/Geometry.Kwall+1/ColdSide.HeatTransfer.hcoeff + Geometry.Rfc + Geometry.Rfh)=1; 593 690 594 691 "Duty" … … 614 711 615 712 "Number of Units Transference for Hot Side" 616 HotSide.HeatTransfer.NTU*HotSide.HeatTransfer.WCp = Thermal.Ud* Atotal;713 HotSide.HeatTransfer.NTU*HotSide.HeatTransfer.WCp = Thermal.Ud*Geometry.Atotal; 617 714 618 715 "Number of Units Transference for Cold Side" 619 ColdSide.HeatTransfer.NTU*ColdSide.HeatTransfer.WCp = Thermal.Ud* Atotal;716 ColdSide.HeatTransfer.NTU*ColdSide.HeatTransfer.WCp = Thermal.Ud*Geometry.Atotal; 620 717 621 718 if Thermal.Eft >= 1 #To be Fixed: Effectiveness in true counter flow !
Note: See TracChangeset
for help on using the changeset viewer.