- Timestamp:
- Nov 27, 2006, 6:37:59 PM (16 years ago)
- Location:
- mso
- Files:
-
- 4 added
- 15 edited
Legend:
- Unmodified
- Added
- Removed
-
mso/eml/heat_exchangers
-
Property
svn:ignore
set to
*.dll
-
Property
svn:ignore
set to
-
mso/eml/heat_exchangers/HEX_Engine.mso
r45 r68 38 38 * Basic_Pdrop : Hot and Cold Side Pressure Drop 39 39 *---------------------------------------------------------------------- 40 * - Main_DoublePipe : Double Pipe Heat Exchanger Block 41 * DoublePipe_HeatTransfer 42 * DoublePipe_PressureDrop 43 *---------------------------------------------------------------------- 40 44 * Author: Gerson Balbueno Bicca 41 45 * $Id$ … … 277 281 Cold as Basic_Pdrop;# Cold Stream 278 282 end 283 284 Model DoublePipe_HeatTransfer 285 #===================================================================== 286 # Heat Transfer 287 #===================================================================== 288 PARAMETERS 289 As as area (Brief="Cross Sectional Area for Flow",Default=0.05,Lower=1e-8); 290 Dh as length (Brief="Hydraulic Diameter of Pipe for Heat Transfer",Lower=1e-8); 291 292 VARIABLES 293 Re as positive (Brief="Reynolds Number",Default=100,Lower=1); 294 hcoeff as heat_trans_coeff (Brief="Film Coefficient",Default=1,Lower=1e-12, Upper=1e6); 295 PR as positive (Brief="Prandtl Number",Default=0.5,Lower=1e-8); 296 Phi as positive (Brief="Phi Correction",Default=1,Lower=1e-3); 297 Vmean as velocity (Brief="Tube Velocity",Lower=1e-8); 298 end 299 300 Model DoublePipe_PressureDrop 301 #===================================================================== 302 # Pressure Drop 303 #===================================================================== 304 PARAMETERS 305 Dh as length (Brief="Hydraulic Diameter of Pipe for Pressure Drop",Lower=1e-6); 306 307 VARIABLES 308 Pdrop as press_delta (Brief="Pressure Drop",Default=0.01, Lower=1e-10); 309 fi as fricfactor (Brief="Friction Factor", Default=0.05, Lower=1e-10, Upper=2000); 310 Re as positive (Brief="Reynolds Number",Default=100,Lower=1); 311 end 312 313 Model Main_DoublePipe 314 VARIABLES 315 HeatTransfer as DoublePipe_HeatTransfer; 316 PressureDrop as DoublePipe_PressureDrop; 317 end -
mso/eml/heat_exchangers/HeatExchangerDetailed.mso
r45 r68 165 165 then 166 166 167 "H eat Capacity Hot Stream"167 "Hot Stream Average Heat Capacity" 168 168 Properties.Hot.Average.Cp = PP.LiquidCp(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z); 169 170 "Hot Stream Inlet Heat Capacity" 169 171 Properties.Hot.Inlet.Cp = PP.LiquidCp(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z); 172 173 "Hot Stream Outlet Heat Capacity" 170 174 Properties.Hot.Outlet.Cp = PP.LiquidCp(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z); 171 175 172 " Mass Density Hot Stream"176 "Hot Stream Average Mass Density" 173 177 Properties.Hot.Average.rho = PP.LiquidDensity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z); 178 179 "Hot Stream Inlet Mass Density" 174 180 Properties.Hot.Inlet.rho = PP.LiquidDensity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z); 181 182 "Hot Stream Outlet Mass Density" 175 183 Properties.Hot.Outlet.rho = PP.LiquidDensity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z); 176 184 177 " Viscosity Hot Stream"185 "Hot Stream Average Viscosity" 178 186 Properties.Hot.Average.Mu = PP.LiquidViscosity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z); 187 188 "Hot Stream Inlet Viscosity" 179 189 Properties.Hot.Inlet.Mu = PP.LiquidViscosity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z); 190 191 "Hot Stream Outlet Viscosity" 180 192 Properties.Hot.Outlet.Mu = PP.LiquidViscosity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z); 181 193 182 "Conductivity Hot Stream" 183 Properties.Hot.Average.K = PP.LiquidThermalConductivity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z); 184 Properties.Hot.Inlet.K = PP.LiquidThermalConductivity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z); 185 Properties.Hot.Outlet.K = PP.LiquidThermalConductivity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z); 186 187 "Heat Capacity Hot Stream" 188 Properties.Hot.Wall.Cp = PP.LiquidCp(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z); 189 190 "Viscosity Hot Stream" 191 Properties.Hot.Wall.Mu = PP.LiquidViscosity(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z); 192 193 "Conductivity Hot Stream" 194 Properties.Hot.Wall.K = PP.LiquidThermalConductivity(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z); 195 196 197 else 198 199 "Heat Capacity Hot Stream" 194 "Hot Stream Average Conductivity" 195 Properties.Hot.Average.K = PP.LiquidThermalConductivity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z); 196 197 "Hot Stream Inlet Conductivity" 198 Properties.Hot.Inlet.K = PP.LiquidThermalConductivity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z); 199 200 "Hot Stream Outlet Conductivity" 201 Properties.Hot.Outlet.K = PP.LiquidThermalConductivity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z); 202 203 "Hot Stream Heat Capacity at Wall Temperature" 204 Properties.Hot.Wall.Cp = PP.LiquidCp(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z); 205 206 "Hot Stream Viscosity at Wall Temperature" 207 Properties.Hot.Wall.Mu = PP.LiquidViscosity(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z); 208 209 "Hot Stream Conductivity at Wall Temperature" 210 Properties.Hot.Wall.K = PP.LiquidThermalConductivity(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z); 211 212 213 else 214 215 "Hot Stream Average Heat Capacity" 200 216 Properties.Hot.Average.Cp = PP.VapourCp(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z); 217 218 "Hot Stream Inlet Heat Capacity" 201 219 Properties.Hot.Inlet.Cp = PP.VapourCp(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z); 220 221 "Hot Stream Outlet Heat Capacity" 202 222 Properties.Hot.Outlet.Cp = PP.VapourCp(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z); 203 223 204 " Mass Density Hot Stream"224 "Hot Stream Average Mass Density" 205 225 Properties.Hot.Average.rho = PP.VapourDensity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z); 226 227 "Hot Stream Inlet Mass Density" 206 228 Properties.Hot.Inlet.rho = PP.VapourDensity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z); 229 230 "Hot Stream Outlet Mass Density" 207 231 Properties.Hot.Outlet.rho = PP.VapourDensity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z); 208 232 209 " Viscosity Hot Stream"233 "Hot Stream Average Viscosity" 210 234 Properties.Hot.Average.Mu = PP.VapourViscosity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z); 235 236 "Hot Stream Inlet Viscosity" 211 237 Properties.Hot.Inlet.Mu = PP.VapourViscosity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z); 238 239 "Hot Stream Outlet Viscosity" 212 240 Properties.Hot.Outlet.Mu = PP.VapourViscosity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z); 213 241 214 "Conductivity Hot Stream" 215 Properties.Hot.Average.K = PP.VapourThermalConductivity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z); 216 Properties.Hot.Inlet.K = PP.VapourThermalConductivity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z); 217 Properties.Hot.Outlet.K = PP.VapourThermalConductivity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z); 218 219 "Heat Capacity Hot Stream" 220 Properties.Hot.Wall.Cp = PP.VapourCp(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z); 221 222 "Viscosity Hot Stream" 223 Properties.Hot.Wall.Mu = PP.VapourViscosity(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z); 224 225 "Conductivity Hot Stream" 226 Properties.Hot.Wall.K = PP.VapourThermalConductivity(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z); 242 "Hot Stream Average Conductivity" 243 Properties.Hot.Average.K = PP.VapourThermalConductivity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z); 244 245 "Hot Stream Inlet Conductivity" 246 Properties.Hot.Inlet.K = PP.VapourThermalConductivity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z); 247 248 "Hot Stream Outlet Conductivity" 249 Properties.Hot.Outlet.K = PP.VapourThermalConductivity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z); 250 251 "Hot Stream Heat Capacity at Wall Temperature" 252 Properties.Hot.Wall.Cp = PP.VapourCp(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z); 253 254 "Hot Stream Viscosity at Wall Temperature" 255 Properties.Hot.Wall.Mu = PP.VapourViscosity(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z); 256 257 "Hot Stream Conductivity at Wall Temperature" 258 Properties.Hot.Wall.K = PP.VapourThermalConductivity(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z); 227 259 228 260 … … 423 455 424 456 "Shell Side Phi correction" 425 Shell.HeatTransfer.Phi = HE.PhiCorrection(Properties.Hot.Average.Mu,Properties.Hot.Wall.Mu);457 Shell.HeatTransfer.Phi = HE.PhiCorrection(Properties.Hot.Average.Mu,Properties.Hot.Wall.Mu); 426 458 427 459 "Tube Side Phi correction" … … 431 463 432 464 "Shell Side Phi correction" 433 Shell.HeatTransfer.Phi = HE.PhiCorrection(Properties.Cold.Average.Mu,Properties.Cold.Wall.Mu);465 Shell.HeatTransfer.Phi = HE.PhiCorrection(Properties.Cold.Average.Mu,Properties.Cold.Wall.Mu); 434 466 435 467 "Tube Side Phi correction" … … 449 481 450 482 "Hot Wall Temperature" 451 Properties.Hot.Wall.Twall = (Properties.Hot.Average.T+Properties.Cold.Average.T)/2;483 Properties.Hot.Wall.Twall = (Properties.Hot.Average.T+Properties.Cold.Average.T)/2; 452 484 453 485 "ColdWall Temperature" … … 455 487 456 488 "Tube Side Velocity" 457 Tubes.HeatTransfer.Vtube = Properties.Cold.Inlet.Fw*Tpass/((Pi*Ditube*Ditube/4)*Properties.Cold.Average.rho*Ntt);489 Tubes.HeatTransfer.Vtube = Properties.Cold.Inlet.Fw*Tpass/((Pi*Ditube*Ditube/4)*Properties.Cold.Average.rho*Ntt); 458 490 459 491 "Tube Side Reynolds Number" 460 Tubes.HeatTransfer.Re =(Properties.Cold.Average.rho*Tubes.HeatTransfer.Vtube*Ditube)/Properties.Cold.Average.Mu;492 Tubes.HeatTransfer.Re = (Properties.Cold.Average.rho*Tubes.HeatTransfer.Vtube*Ditube)/Properties.Cold.Average.Mu; 461 493 462 494 "Tube Side Prandtl Number" 463 Tubes.HeatTransfer.PR = ((Properties.Cold.Average.Cp/Properties.Cold.Average.Mw)*Properties.Cold.Average.Mu)/Properties.Cold.Average.K;495 Tubes.HeatTransfer.PR = ((Properties.Cold.Average.Cp/Properties.Cold.Average.Mw)*Properties.Cold.Average.Mu)/Properties.Cold.Average.K; 464 496 465 497 "Tube Side Prandtl Number at Wall Temperature" 466 Tubes.HeatTransfer.PRw = ((Properties.Cold.Wall.Cp/Properties.Cold.Average.Mw)*Properties.Cold.Wall.Mu)/Properties.Cold.Wall.K;498 Tubes.HeatTransfer.PRw = ((Properties.Cold.Wall.Cp/Properties.Cold.Average.Mw)*Properties.Cold.Wall.Mu)/Properties.Cold.Wall.K; 467 499 468 500 "Tube Side Film Coefficient" 469 Tubes.HeatTransfer.htube = HE.TubeFilmCoeff(Tubes.HeatTransfer.Re,Tubes.HeatTransfer.PR,Properties.Cold.Average.K)*Tubes.HeatTransfer.Phi;501 Tubes.HeatTransfer.htube = HE.TubeFilmCoeff(Tubes.HeatTransfer.Re,Tubes.HeatTransfer.PR,Properties.Cold.Average.K)*Tubes.HeatTransfer.Phi; 470 502 471 503 "Shell Side Prandtl Number" 472 Shell.HeatTransfer.PR = ((Properties.Hot.Average.Cp/Properties.Hot.Average.Mw)*Properties.Hot.Average.Mu)/Properties.Hot.Average.K;504 Shell.HeatTransfer.PR = ((Properties.Hot.Average.Cp/Properties.Hot.Average.Mw)*Properties.Hot.Average.Mu)/Properties.Hot.Average.K; 473 505 474 506 "Shell Side Prandtl Number at Wall Temperature" 475 Shell.HeatTransfer.PRw = ((Properties.Hot.Wall.Cp/Properties.Hot.Average.Mw)*Properties.Hot.Wall.Mu)/Properties.Hot.Wall.K;507 Shell.HeatTransfer.PRw = ((Properties.Hot.Wall.Cp/Properties.Hot.Average.Mw)*Properties.Hot.Wall.Mu)/Properties.Hot.Wall.K; 476 508 477 509 "Tube Side Pressure Drop" … … 497 529 498 530 "Shell Pressure Drop Outlet Nozzle" 499 Shell.PressureDrop.Pdnozzle_out = HE.DeltaPshellNozzleOut(Properties.Hot.Outlet.rho,Properties.Hot.Inlet.Fw);531 Shell.PressureDrop.Pdnozzle_out = HE.DeltaPshellNozzleOut(Properties.Hot.Outlet.rho,Properties.Hot.Inlet.Fw); 500 532 501 533 "Velocity Shell Side Outlet Nozzle" … … 621 653 622 654 VARIABLES 655 623 656 DT0 as temp_delta (Brief="Temperature Difference at Inlet",Lower=1); 624 657 DTL as temp_delta (Brief="Temperature Difference at Outlet",Lower=1); … … 1295 1328 VARIABLES 1296 1329 1297 LMTD(Nshell) as temp_delta (Brief="Logarithmic Mean Temperature Difference",Lower= 10);1298 Fc(Nshell) as positive (Brief="LMTD Correction Factor",Lower=0. 75);1299 MTD(Nshell) as temp_delta (Brief="Mean Temperature Difference" );1330 LMTD(Nshell) as temp_delta (Brief="Logarithmic Mean Temperature Difference",Lower=5); 1331 Fc(Nshell) as positive (Brief="LMTD Correction Factor",Lower=0.5); 1332 MTD(Nshell) as temp_delta (Brief="Mean Temperature Difference",Lower=5); 1300 1333 Unity(Nshell) as HeatExchangerDetailed_Basic; 1301 1334 … … 1515 1548 1516 1549 "Overall Heat Transfer Coefficient" 1517 Unity(i).Details.U =1/(Dotube/(Unity(i).Tubes.HeatTransfer.htube*Ditube)+(Dotube*ln(Dotube/Ditube)/(2*Kwall))+(1/(Unity(i).Shell.HeatTransfer.hshell)));1550 Unity(i).Details.U*(Dotube/(Unity(i).Tubes.HeatTransfer.htube*Ditube)+(Dotube*ln(Dotube/Ditube)/(2*Kwall))+(1/(Unity(i).Shell.HeatTransfer.hshell)))=1; 1518 1551 1519 1552 "Exchange Surface Area" … … 1573 1606 1574 1607 "Exchange Surface Area" 1575 Unity(i).Details.Q = Unity(i).Details.U*Unity(i).Details.A*Fc(i)*LMTD(i); 1608 # Unity(i).Details.Q = Unity(i).Details.U*Unity(i).Details.A*Fc(i)*LMTD(i); 1609 Unity(i).Details.Q = Unity(i).Details.U*Pi*Dotube*Ntt*Ltube*Fc(i)*LMTD(i); 1576 1610 1577 1611 "Mean Temperature Difference" -
mso/eml/heat_exchangers/HeatExchangerDiscretized.mso
r45 r68 60 60 then 61 61 62 " Heat Capacity Cold Stream"62 "Cold Stream Average Heat Capacity" 63 63 Properties.Cold.Average.Cp = PP.LiquidCp(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z); 64 65 "Cold Stream Inlet Heat Capacity" 64 66 Properties.Cold.Inlet.Cp = PP.LiquidCp(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z); 67 68 "Cold Stream Outlet Heat Capacity" 65 69 Properties.Cold.Outlet.Cp = PP.LiquidCp(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z); 66 70 67 " Mass Density Cold Stream"71 "Cold Stream Average Mass Density" 68 72 Properties.Cold.Average.rho = PP.LiquidDensity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z); 73 74 "Cold Stream Inlet Mass Density" 69 75 Properties.Cold.Inlet.rho = PP.LiquidDensity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z); 76 77 "Cold Stream Outlet Mass Density" 70 78 Properties.Cold.Outlet.rho = PP.LiquidDensity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z); 71 79 72 "Viscosity Cold Stream" 73 Properties.Cold.Average.Mu = PP.LiquidViscosity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z); 74 Properties.Cold.Inlet.Mu = PP.LiquidViscosity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z); 75 Properties.Cold.Outlet.Mu = PP.LiquidViscosity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z); 76 77 "Conductivity Cold Stream" 78 Properties.Cold.Average.K = PP.LiquidThermalConductivity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z); 79 Properties.Cold.Inlet.K = PP.LiquidThermalConductivity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z); 80 Properties.Cold.Outlet.K = PP.LiquidThermalConductivity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z); 81 82 "Heat Capacity Cold Stream" 80 "Cold Stream Average Viscosity" 81 Properties.Cold.Average.Mu = PP.LiquidViscosity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z); 82 83 "Cold Stream inlet Viscosity" 84 Properties.Cold.Inlet.Mu = PP.LiquidViscosity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z); 85 86 "Cold Stream Outlet Viscosity" 87 Properties.Cold.Outlet.Mu = PP.LiquidViscosity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z); 88 89 "Cold Stream Average Conductivity" 90 Properties.Cold.Average.K = PP.LiquidThermalConductivity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z); 91 92 "Cold Stream Inlet Conductivity" 93 Properties.Cold.Inlet.K = PP.LiquidThermalConductivity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z); 94 95 "Cold Stream Outlet Conductivity" 96 Properties.Cold.Outlet.K = PP.LiquidThermalConductivity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z); 97 98 "Cold Stream Heat Capacity at Wall Temperature" 83 99 Properties.Cold.Wall.Cp = PP.LiquidCp(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z); 84 100 85 " Viscosity Cold Stream"86 Properties.Cold.Wall.Mu = PP.LiquidViscosity(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z);87 88 "Co nductivity Cold Stream"89 Properties.Cold.Wall.K = PP.LiquidThermalConductivity(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z);101 "Cold Stream Viscosity at Wall Temperature" 102 Properties.Cold.Wall.Mu = PP.LiquidViscosity(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z); 103 104 "Cold Stream Conductivity at Wall Temperature" 105 Properties.Cold.Wall.K = PP.LiquidThermalConductivity(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z); 90 106 91 107 92 108 else 93 109 94 " Heat Capacity Cold Stream"110 "Cold Stream Average Heat Capacity" 95 111 Properties.Cold.Average.Cp = PP.VapourCp(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z); 112 113 "Cold Stream Inlet Heat Capacity" 96 114 Properties.Cold.Inlet.Cp = PP.VapourCp(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z); 115 116 "Cold Stream Outlet Heat Capacity" 97 117 Properties.Cold.Outlet.Cp = PP.VapourCp(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z); 98 118 99 "Mass Density Cold Stream" 100 Properties.Cold.Average.rho = PP.VapourDensity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z); 119 "Cold Stream Average Mass Density" 120 Properties.Cold.Average.rho = PP.VapourDensity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z); 121 122 "Cold Stream Inlet Mass Density" 101 123 Properties.Cold.Inlet.rho = PP.VapourDensity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z); 124 125 "Cold Stream Outlet Mass Density" 102 126 Properties.Cold.Outlet.rho = PP.VapourDensity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z); 103 127 104 "Viscosity Cold Stream" 105 Properties.Cold.Average.Mu = PP.VapourViscosity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z); 106 Properties.Cold.Inlet.Mu = PP.VapourViscosity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z); 107 Properties.Cold.Outlet.Mu = PP.VapourViscosity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z); 108 109 "Conductivity Cold Stream" 110 Properties.Cold.Average.K = PP.VapourThermalConductivity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z); 111 Properties.Cold.Inlet.K = PP.VapourThermalConductivity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z); 112 Properties.Cold.Outlet.K = PP.VapourThermalConductivity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z); 113 114 "Heat Capacity Cold Stream" 128 "Cold Stream Average Viscosity " 129 Properties.Cold.Average.Mu = PP.VapourViscosity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z); 130 131 "Cold Stream Inlet Viscosity " 132 Properties.Cold.Inlet.Mu = PP.VapourViscosity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z); 133 134 "Cold Stream Outlet Viscosity " 135 Properties.Cold.Outlet.Mu = PP.VapourViscosity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z); 136 137 "Cold Stream Average Conductivity " 138 Properties.Cold.Average.K = PP.VapourThermalConductivity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z); 139 140 "Cold Stream Inlet Conductivity " 141 Properties.Cold.Inlet.K = PP.VapourThermalConductivity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z); 142 143 "Cold Stream Outlet Conductivity " 144 Properties.Cold.Outlet.K = PP.VapourThermalConductivity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z); 145 146 "Cold Stream Heat Capacity at Wall Temperature" 115 147 Properties.Cold.Wall.Cp = PP.VapourCp(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z); 116 148 117 149 118 "Viscosity Cold Stream" 119 Properties.Cold.Wall.Mu = PP.VapourViscosity(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z); 120 121 "Conductivity Cold Stream" 122 Properties.Cold.Wall.K = PP.VapourThermalConductivity(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z); 123 124 125 126 end 150 "Cold Stream Viscosity at Wall Temperature" 151 Properties.Cold.Wall.Mu = PP.VapourViscosity(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z); 152 153 "Cold Stream Conductivity at Wall Temperature" 154 Properties.Cold.Wall.K = PP.VapourThermalConductivity(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z); 155 156 157 158 end 159 127 160 128 161 if Inlet.Hot.v equal 0 … … 130 163 then 131 164 132 "H eat Capacity Hot Stream"165 "Hot Stream Average Heat Capacity" 133 166 Properties.Hot.Average.Cp = PP.LiquidCp(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z); 167 168 "Hot Stream Inlet Heat Capacity" 134 169 Properties.Hot.Inlet.Cp = PP.LiquidCp(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z); 170 171 "Hot Stream Outlet Heat Capacity" 135 172 Properties.Hot.Outlet.Cp = PP.LiquidCp(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z); 136 173 137 " Mass Density Hot Stream"174 "Hot Stream Average Mass Density" 138 175 Properties.Hot.Average.rho = PP.LiquidDensity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z); 176 177 "Hot Stream Inlet Mass Density" 139 178 Properties.Hot.Inlet.rho = PP.LiquidDensity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z); 179 180 "Hot Stream Outlet Mass Density" 140 181 Properties.Hot.Outlet.rho = PP.LiquidDensity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z); 141 182 142 " Viscosity Hot Stream"183 "Hot Stream Average Viscosity" 143 184 Properties.Hot.Average.Mu = PP.LiquidViscosity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z); 185 186 "Hot Stream Inlet Viscosity" 144 187 Properties.Hot.Inlet.Mu = PP.LiquidViscosity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z); 188 189 "Hot Stream Outlet Viscosity" 145 190 Properties.Hot.Outlet.Mu = PP.LiquidViscosity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z); 146 191 147 " Conductivity Hot Stream"192 "Hot Stream Average Conductivity" 148 193 Properties.Hot.Average.K = PP.LiquidThermalConductivity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z); 149 Properties.Hot.Inlet.K = PP.LiquidThermalConductivity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z); 194 195 "Hot Stream Inlet Conductivity" 196 Properties.Hot.Inlet.K = PP.LiquidThermalConductivity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z); 197 198 "Hot Stream Outlet Conductivity" 150 199 Properties.Hot.Outlet.K = PP.LiquidThermalConductivity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z); 151 200 152 "H eat Capacity Hot Stream"153 Properties.Hot.Wall.Cp = PP.LiquidCp(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z);154 155 " Viscosity Hot Stream"156 Properties.Hot.Wall.Mu = PP.LiquidViscosity(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z);157 158 " Conductivity Hot Stream"159 Properties.Hot.Wall.K = PP.LiquidThermalConductivity(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z);201 "Hot Stream Heat Capacity at Wall Temperature" 202 Properties.Hot.Wall.Cp = PP.LiquidCp(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z); 203 204 "Hot Stream Viscosity at Wall Temperature" 205 Properties.Hot.Wall.Mu = PP.LiquidViscosity(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z); 206 207 "Hot Stream Conductivity at Wall Temperature" 208 Properties.Hot.Wall.K = PP.LiquidThermalConductivity(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z); 160 209 161 210 162 211 else 163 212 164 "H eat Capacity Hot Stream"213 "Hot Stream Average Heat Capacity" 165 214 Properties.Hot.Average.Cp = PP.VapourCp(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z); 215 216 "Hot Stream Inlet Heat Capacity" 166 217 Properties.Hot.Inlet.Cp = PP.VapourCp(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z); 218 219 "Hot Stream Outlet Heat Capacity" 167 220 Properties.Hot.Outlet.Cp = PP.VapourCp(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z); 168 221 169 " Mass Density Hot Stream"222 "Hot Stream Average Mass Density" 170 223 Properties.Hot.Average.rho = PP.VapourDensity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z); 224 225 "Hot Stream Inlet Mass Density" 171 226 Properties.Hot.Inlet.rho = PP.VapourDensity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z); 227 228 "Hot Stream Outlet Mass Density" 172 229 Properties.Hot.Outlet.rho = PP.VapourDensity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z); 173 230 174 " Viscosity Hot Stream"231 "Hot Stream Average Viscosity" 175 232 Properties.Hot.Average.Mu = PP.VapourViscosity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z); 233 234 "Hot Stream Inlet Viscosity" 176 235 Properties.Hot.Inlet.Mu = PP.VapourViscosity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z); 236 237 "Hot Stream Outlet Viscosity" 177 238 Properties.Hot.Outlet.Mu = PP.VapourViscosity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z); 178 239 179 " Conductivity Hot Stream"240 "Hot Stream Average Conductivity" 180 241 Properties.Hot.Average.K = PP.VapourThermalConductivity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z); 181 Properties.Hot.Inlet.K = PP.VapourThermalConductivity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z); 242 243 "Hot Stream Inlet Conductivity" 244 Properties.Hot.Inlet.K = PP.VapourThermalConductivity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z); 245 246 "Hot Stream Outlet Conductivity" 182 247 Properties.Hot.Outlet.K = PP.VapourThermalConductivity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z); 183 248 184 "H eat Capacity Hot Stream"185 Properties.Hot.Wall.Cp = PP.VapourCp(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z);186 187 " Viscosity Hot Stream"188 Properties.Hot.Wall.Mu = PP.VapourViscosity(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z);189 190 " Conductivity Hot Stream"191 Properties.Hot.Wall.K = PP.VapourThermalConductivity(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z);249 "Hot Stream Heat Capacity at Wall Temperature" 250 Properties.Hot.Wall.Cp = PP.VapourCp(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z); 251 252 "Hot Stream Viscosity at Wall Temperature" 253 Properties.Hot.Wall.Mu = PP.VapourViscosity(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z); 254 255 "Hot Stream Conductivity at Wall Temperature" 256 Properties.Hot.Wall.K = PP.VapourThermalConductivity(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z); 192 257 193 258 … … 273 338 "TEMA E Shell Effectiveness" 274 339 Eft = HE.EshellEffectiveness(Details.Cr,Details.NTU); 275 340 276 341 end 277 342 … … 292 357 MTD = Fc*LMTD; 293 358 294 #"LMTD Correction Factor"295 #Fc = HE.EshellCorrectionFactor(Inlet.Hot.T,Outlet.Hot.T,Inlet.Cold.T,Outlet.Cold.T);359 "LMTD Correction Factor" 360 Fc = HE.EshellCorrectionFactor(Inlet.Hot.T,Outlet.Hot.T,Inlet.Cold.T,Outlet.Cold.T); 296 361 297 362 "Temperature Difference at Inlet" … … 745 810 746 811 "Overall Heat Transfer Coefficient" 747 Unity(i).Details.U = 1/(Dotube/(Unity(i).Tubes.HeatTransfer.htube*Ditube)+(Dotube*ln(Dotube/Ditube)/(2*Kwall))+(1/(Unity(i).Shell.HeatTransfer.hshell)));812 Unity(i).Details.U*(Dotube/(Unity(i).Tubes.HeatTransfer.htube*Ditube)+(Dotube*ln(Dotube/Ditube)/(2*Kwall))+(1/(Unity(i).Shell.HeatTransfer.hshell)))=1; 748 813 749 814 "Exchange Surface Area" … … 955 1020 Unity(Nb+1).Shell.HeatTransfer.Sm = HE.CrossFlowArea(Unity(Nb+1).Baffles.Lso); 956 1021 957 "LMTD Correction Factor"958 Unity.Fc = HE.EshellCorrectionFactor(Unity(1).Inlet.Hot.T,Unity(Nb+1).Outlet.Hot.T,Unity(Nb+1).Inlet.Cold.T,Unity(1).Outlet.Cold.T);1022 #"LMTD Correction Factor" 1023 # Unity.Fc = HE.EshellCorrectionFactor(Unity(1).Inlet.Hot.T,Unity(Nb+1).Outlet.Hot.T,Unity(Nb+1).Inlet.Cold.T,Unity(1).Outlet.Cold.T); 959 1024 960 1025 if side equal 1 … … 1234 1299 1235 1300 "Overall Heat Transfer Coefficient" 1236 Unity(i).Details.U = 1/(Dotube/(Unity(i).Tubes.HeatTransfer.htube*Ditube)+(Dotube*ln(Dotube/Ditube)/(2*Kwall))+(1/(Unity(i).Shell.HeatTransfer.hshell)));1301 Unity(i).Details.U*(Dotube/(Unity(i).Tubes.HeatTransfer.htube*Ditube)+(Dotube*ln(Dotube/Ditube)/(2*Kwall))+(1/(Unity(i).Shell.HeatTransfer.hshell)))=1; 1237 1302 1238 1303 "Exchange Surface Area" -
mso/eml/heat_exchangers/HeatExchangerSimplified.mso
r45 r68 53 53 54 54 55 56 55 if Inlet.Cold.v equal 0 56 57 57 then 58 "Heat Capacity Cold Stream" 59 Properties.Cold.Average.Cp = PP.LiquidCp(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z); 60 Properties.Cold.Inlet.Cp = PP.LiquidCp(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z); 61 Properties.Cold.Outlet.Cp = PP.LiquidCp(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z); 62 63 "Mass Density Cold Stream" 64 Properties.Cold.Average.rho = PP.LiquidDensity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z); 65 Properties.Cold.Inlet.rho = PP.LiquidDensity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z); 66 Properties.Cold.Outlet.rho = PP.LiquidDensity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z); 67 68 "Viscosity Cold Stream" 69 Properties.Cold.Average.Mu = PP.LiquidViscosity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z); 70 Properties.Cold.Inlet.Mu = PP.LiquidViscosity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z); 71 Properties.Cold.Outlet.Mu = PP.LiquidViscosity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z); 72 73 "Conductivity Cold Stream" 74 Properties.Cold.Average.K = PP.LiquidThermalConductivity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z); 75 Properties.Cold.Inlet.K = PP.LiquidThermalConductivity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z); 76 Properties.Cold.Outlet.K = PP.LiquidThermalConductivity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z); 77 78 "Heat Capacity Cold Stream" 79 Properties.Cold.Wall.Cp = PP.LiquidCp(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z); 80 81 "Viscosity Cold Stream" 82 Properties.Cold.Wall.Mu = PP.LiquidViscosity(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z); 83 84 "Conductivity Cold Stream" 85 Properties.Cold.Wall.K = PP.LiquidThermalConductivity(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z); 86 87 88 else 89 90 "Heat Capacity Cold Stream" 58 59 "Cold Stream Average Heat Capacity" 60 Properties.Cold.Average.Cp = PP.LiquidCp(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z); 61 62 "Cold Stream Inlet Heat Capacity" 63 Properties.Cold.Inlet.Cp = PP.LiquidCp(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z); 64 65 "Cold Stream Outlet Heat Capacity" 66 Properties.Cold.Outlet.Cp = PP.LiquidCp(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z); 67 68 "Cold Stream Average Mass Density" 69 Properties.Cold.Average.rho = PP.LiquidDensity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z); 70 71 "Cold Stream Inlet Mass Density" 72 Properties.Cold.Inlet.rho = PP.LiquidDensity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z); 73 74 "Cold Stream Outlet Mass Density" 75 Properties.Cold.Outlet.rho = PP.LiquidDensity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z); 76 77 "Cold Stream Average Viscosity" 78 Properties.Cold.Average.Mu = PP.LiquidViscosity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z); 79 80 "Cold Stream inlet Viscosity" 81 Properties.Cold.Inlet.Mu = PP.LiquidViscosity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z); 82 83 "Cold Stream Outlet Viscosity" 84 Properties.Cold.Outlet.Mu = PP.LiquidViscosity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z); 85 86 "Cold Stream Average Conductivity" 87 Properties.Cold.Average.K = PP.LiquidThermalConductivity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z); 88 89 "Cold Stream Inlet Conductivity" 90 Properties.Cold.Inlet.K = PP.LiquidThermalConductivity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z); 91 92 "Cold Stream Outlet Conductivity" 93 Properties.Cold.Outlet.K = PP.LiquidThermalConductivity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z); 94 95 "Cold Stream Heat Capacity at Wall Temperature" 96 Properties.Cold.Wall.Cp = PP.LiquidCp(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z); 97 98 "Cold Stream Viscosity at Wall Temperature" 99 Properties.Cold.Wall.Mu = PP.LiquidViscosity(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z); 100 101 "Cold Stream Conductivity at Wall Temperature" 102 Properties.Cold.Wall.K = PP.LiquidThermalConductivity(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z); 103 104 105 else 106 107 "Cold Stream Average Heat Capacity" 91 108 Properties.Cold.Average.Cp = PP.VapourCp(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z); 109 110 "Cold Stream Inlet Heat Capacity" 92 111 Properties.Cold.Inlet.Cp = PP.VapourCp(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z); 112 113 "Cold Stream Outlet Heat Capacity" 93 114 Properties.Cold.Outlet.Cp = PP.VapourCp(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z); 94 115 95 "Mass Density Cold Stream" 96 Properties.Cold.Average.rho = PP.VapourDensity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z); 97 Properties.Cold.Inlet.rho = PP.VapourDensity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z); 98 Properties.Cold.Outlet.rho = PP.VapourDensity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z); 99 100 "Viscosity Cold Stream" 101 Properties.Cold.Average.Mu = PP.VapourViscosity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z); 102 Properties.Cold.Inlet.Mu = PP.VapourViscosity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z); 103 Properties.Cold.Outlet.Mu = PP.VapourViscosity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z); 104 105 "Conductivity Cold Stream" 106 Properties.Cold.Average.K = PP.VapourThermalConductivity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z); 107 Properties.Cold.Inlet.K = PP.VapourThermalConductivity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z); 108 Properties.Cold.Outlet.K = PP.VapourThermalConductivity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z); 109 110 "Heat Capacity Cold Stream" 111 Properties.Cold.Wall.Cp = PP.VapourCp(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z); 112 113 114 "Viscosity Cold Stream" 115 Properties.Cold.Wall.Mu = PP.VapourViscosity(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z); 116 117 "Conductivity Cold Stream" 118 Properties.Cold.Wall.K = PP.VapourThermalConductivity(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z); 116 "Cold Stream Average Mass Density" 117 Properties.Cold.Average.rho = PP.VapourDensity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z); 118 119 "Cold Stream Inlet Mass Density" 120 Properties.Cold.Inlet.rho = PP.VapourDensity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z); 121 122 "Cold Stream Outlet Mass Density" 123 Properties.Cold.Outlet.rho = PP.VapourDensity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z); 124 125 "Cold Stream Average Viscosity " 126 Properties.Cold.Average.Mu = PP.VapourViscosity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z); 127 128 "Cold Stream Inlet Viscosity " 129 Properties.Cold.Inlet.Mu = PP.VapourViscosity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z); 130 131 "Cold Stream Outlet Viscosity " 132 Properties.Cold.Outlet.Mu = PP.VapourViscosity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z); 133 134 "Cold Stream Average Conductivity " 135 Properties.Cold.Average.K = PP.VapourThermalConductivity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z); 136 137 "Cold Stream Inlet Conductivity " 138 Properties.Cold.Inlet.K = PP.VapourThermalConductivity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z); 139 140 "Cold Stream Outlet Conductivity " 141 Properties.Cold.Outlet.K = PP.VapourThermalConductivity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z); 142 143 "Cold Stream Heat Capacity at Wall Temperature" 144 Properties.Cold.Wall.Cp = PP.VapourCp(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z); 145 146 147 "Cold Stream Viscosity at Wall Temperature" 148 Properties.Cold.Wall.Mu = PP.VapourViscosity(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z); 149 150 "Cold Stream Conductivity at Wall Temperature" 151 Properties.Cold.Wall.K = PP.VapourThermalConductivity(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z); 119 152 120 153 … … 126 159 then 127 160 128 "H eat Capacity Hot Stream"161 "Hot Stream Average Heat Capacity" 129 162 Properties.Hot.Average.Cp = PP.LiquidCp(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z); 163 164 "Hot Stream Inlet Heat Capacity" 130 165 Properties.Hot.Inlet.Cp = PP.LiquidCp(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z); 166 167 "Hot Stream Outlet Heat Capacity" 131 168 Properties.Hot.Outlet.Cp = PP.LiquidCp(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z); 132 169 133 " Mass Density Hot Stream"170 "Hot Stream Average Mass Density" 134 171 Properties.Hot.Average.rho = PP.LiquidDensity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z); 172 173 "Hot Stream Inlet Mass Density" 135 174 Properties.Hot.Inlet.rho = PP.LiquidDensity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z); 175 176 "Hot Stream Outlet Mass Density" 136 177 Properties.Hot.Outlet.rho = PP.LiquidDensity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z); 137 178 138 " Viscosity Hot Stream"179 "Hot Stream Average Viscosity" 139 180 Properties.Hot.Average.Mu = PP.LiquidViscosity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z); 181 182 "Hot Stream Inlet Viscosity" 140 183 Properties.Hot.Inlet.Mu = PP.LiquidViscosity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z); 184 185 "Hot Stream Outlet Viscosity" 141 186 Properties.Hot.Outlet.Mu = PP.LiquidViscosity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z); 142 187 143 " Conductivity Hot Stream"188 "Hot Stream Average Conductivity" 144 189 Properties.Hot.Average.K = PP.LiquidThermalConductivity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z); 190 191 "Hot Stream Inlet Conductivity" 145 192 Properties.Hot.Inlet.K = PP.LiquidThermalConductivity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z); 193 194 "Hot Stream Outlet Conductivity" 146 195 Properties.Hot.Outlet.K = PP.LiquidThermalConductivity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z); 147 196 148 "H eat Capacity Hot Stream"197 "Hot Stream Heat Capacity at Wall Temperature" 149 198 Properties.Hot.Wall.Cp = PP.LiquidCp(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z); 150 199 151 " Viscosity Hot Stream"200 "Hot Stream Viscosity at Wall Temperature" 152 201 Properties.Hot.Wall.Mu = PP.LiquidViscosity(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z); 153 202 154 " Conductivity Hot Stream"203 "Hot Stream Conductivity at Wall Temperature" 155 204 Properties.Hot.Wall.K = PP.LiquidThermalConductivity(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z); 156 205 … … 158 207 else 159 208 160 "H eat Capacity Hot Stream"209 "Hot Stream Average Heat Capacity" 161 210 Properties.Hot.Average.Cp = PP.VapourCp(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z); 211 212 "Hot Stream Inlet Heat Capacity" 162 213 Properties.Hot.Inlet.Cp = PP.VapourCp(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z); 214 215 "Hot Stream Outlet Heat Capacity" 163 216 Properties.Hot.Outlet.Cp = PP.VapourCp(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z); 164 217 165 " Mass Density Hot Stream"218 "Hot Stream Average Mass Density" 166 219 Properties.Hot.Average.rho = PP.VapourDensity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z); 220 221 "Hot Stream Inlet Mass Density" 167 222 Properties.Hot.Inlet.rho = PP.VapourDensity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z); 223 224 "Hot Stream Outlet Mass Density" 168 225 Properties.Hot.Outlet.rho = PP.VapourDensity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z); 169 226 170 " Viscosity Hot Stream"227 "Hot Stream Average Viscosity" 171 228 Properties.Hot.Average.Mu = PP.VapourViscosity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z); 229 230 "Hot Stream Inlet Viscosity" 172 231 Properties.Hot.Inlet.Mu = PP.VapourViscosity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z); 232 233 "Hot Stream Outlet Viscosity" 173 234 Properties.Hot.Outlet.Mu = PP.VapourViscosity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z); 174 235 175 " Conductivity Hot Stream"236 "Hot Stream Average Conductivity" 176 237 Properties.Hot.Average.K = PP.VapourThermalConductivity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z); 238 239 "Hot Stream Inlet Conductivity" 177 240 Properties.Hot.Inlet.K = PP.VapourThermalConductivity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z); 241 242 "Hot Stream Outlet Conductivity" 178 243 Properties.Hot.Outlet.K = PP.VapourThermalConductivity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z); 179 244 180 "H eat Capacity Hot Stream"245 "Hot Stream Heat Capacity at Wall Temperature" 181 246 Properties.Hot.Wall.Cp = PP.VapourCp(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z); 182 247 183 " Viscosity Hot Stream"248 "Hot Stream Viscosity at Wall Temperature" 184 249 Properties.Hot.Wall.Mu = PP.VapourViscosity(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z); 185 250 186 " Conductivity Hot Stream"251 "Hot Stream Conductivity at Wall Temperature" 187 252 Properties.Hot.Wall.K = PP.VapourThermalConductivity(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z); 188 253 189 254 190 255 end 256 191 257 192 258 #===================================================================== … … 417 483 418 484 Model HeatExchanger_NTU as Heatex_Basic_NTU 419 420 EQUATIONS 485 486 PARAMETERS 487 488 Side as Integer (Brief="Flow Direction",Lower=0,Upper=1); 489 490 SET 491 492 Side = HE.FlowDir(); # Return Flow Direction 493 494 EQUATIONS 495 496 if Details.Cr equal 0 497 498 then 421 499 "Effectiveness" 422 Eft=HE.Effectiveness(Details.Cr,Details.NTU); 500 Eft = 1-exp(-Details.NTU); 501 502 else 503 504 if Side equal 0 505 506 then 507 "Effectiveness in Cocurrent Flow" 508 Eft = (1-exp(-Details.NTU*(1+Details.Cr)))/(1+Details.Cr); 509 510 else 511 512 if Details.Cr equal 1 513 514 then 515 "Effectiveness in Counter Flow" 516 Eft = Details.NTU/(1+Details.NTU); 517 518 else 519 "Effectiveness in Counter Flow" 520 Eft*(1-Details.Cr*exp(-Details.NTU*(1-Details.Cr))) = (1-exp(-Details.NTU*(1-Details.Cr))); 521 522 end 523 524 end 525 526 527 end 528 423 529 424 530 end -
mso/sample/heat_exchangers/Eshell_Detailed_LMTD.mso
r45 r68 2 2 #=============================================================== 3 3 # Heat Exchanger TEMA E Shell - LMTD Method 4 # converge após várias tentativas5 4 #=============================================================== 6 5 -
mso/sample/heat_exchangers/Eshell_Detailed_NTU.mso
r45 r68 35 35 # Shell Geometrical Parameters 36 36 #===================================================================== 37 exchanger.Tpass = 2;37 exchanger.Tpass = 4; 38 38 exchanger.Dishell = 0.75 *"m"; 39 39 exchanger.Lcf = 0.043 *"m"; … … 53 53 exchanger.Ditube = 0.013395 *"m"; 54 54 exchanger.Dotube = 0.015875 *"m"; 55 exchanger.Kwall = 0. 57 *"kW/m/K";55 exchanger.Kwall = 0.057 *"kW/m/K"; 56 56 exchanger.Donozzle_Tube = 0.203 *"m"; 57 57 exchanger.Dinozzle_Tube = 0.203 *"m"; -
mso/sample/heat_exchangers/Eshell_Discretized_LMTD.mso
r45 r68 100 100 # Simulation Options 101 101 #============================================ 102 mode 103 102 mode = "steady"; 103 guessFile = "E_Shell_Discretized_NTU"; 104 104 end -
mso/sample/heat_exchangers/Eshell_Discretized_NTU.mso
r45 r68 97 97 #============================================ 98 98 mode = "steady"; 99 99 guessFile = "E_Shell_Discretized_NTU"; 100 100 end -
mso/sample/heat_exchangers/Multipass_Detailed.mso
r45 r68 99 99 #============================================ 100 100 mode = "steady"; 101 101 guessFile = "Multipass_LMTD"; 102 102 end -
mso/sample/heat_exchangers/NTU_Method.mso
r45 r68 11 11 streamcold as streamTP; 12 12 13 14 15 13 CONNECTIONS 16 14 … … 20 18 PARAMETERS 21 19 22 PP as CalcObject (File="vrpp.dll"); 23 HE as CalcObject (File="heatex.dll"); 20 PP as CalcObject (File="vrpp"); 24 21 NComp as Integer; 25 22 … … 33 30 PP.Components = ["water"]; 34 31 NComp = PP.NumberOfComponents; 35 HE.HotSide = "Shell";32 exchanger.HE.HotSide = "Shell"; 36 33 37 #HE.TurbulentFlow = "Petukhov";38 34 exchanger.HE.TurbulentFlow = "Petukhov"; 35 # exchanger.HE.TurbulentFlow = "SiederTate"; 39 36 40 37 #===================================================================== -
mso/sample/heat_exchangers/sampleEshell_LMTD.mso
r45 r68 23 23 24 24 exchanger.HE.LMTDcorrection = "Bowmann"; 25 # HE.LMTDcorrection = "Fakeri";25 # exchanger.HE.LMTDcorrection = "Fakeri"; 26 26 27 27 -
mso/sample/heat_exchangers/sampleNTU.mso
r45 r68 9 9 NComp as Integer; 10 10 11 DEVICES 11 12 12 13 DEVICES14 13 exchanger as HeatExchanger_NTU; 15 14 streamhot as streamTP; … … 20 19 PP.LiquidModel = "PR"; 21 20 PP.VapourModel = "PR"; 22 PP.Components = ["water","n-butane", "benzene", "n-octane" ]; 23 NComp = PP.NumberOfComponents; 24 exchanger.HE.FlowDirection = "Cocurrent"; 21 PP.Components = ["water","n-butane", "benzene", "n-octane" ]; 22 NComp = PP.NumberOfComponents; 23 24 exchanger.HE.FlowDirection = "Cocurrent"; 25 25 26 26 CONNECTIONS … … 32 32 33 33 exchanger.Details.U = 210*"W/(m^2*K)"; 34 exchanger.PressureDrop.Hot.Pdrop = 0. 1*"kPa";34 exchanger.PressureDrop.Hot.Pdrop = 0.2*"kPa"; 35 35 exchanger.PressureDrop.Cold.Pdrop = 0.2*"kPa"; 36 36 37 streamhot.F = 20 * "kmol/h"; 37 38 streamhot.T = 450 * "K"; 38 39 streamhot.P = 120 * "kPa"; 39 40 streamhot.z = [1,0,0,0]; 41 40 42 streamcold.F = 10 * "kmol/h"; 41 43 streamcold.P = 120 * "kPa"; 42 44 streamcold.T = 300 * "K"; 43 45 streamcold.z = [0,0.5, 0.1, 0.4]; 44 exchanger.Outlet.Cold.T = 340*"K"; 45 46 47 exchanger.Details.A = 20*"m^2"; 48 46 49 OPTIONS 47 mode = "steady"; 48 end 50 51 mode = "steady"; 52 53 end -
mso/sample/heat_exchangers/samples1.mso
r45 r68 25 25 SET 26 26 27 PP.LiquidModel = " RK";28 PP.VapourModel = " RK";27 PP.LiquidModel = "PR"; 28 PP.VapourModel = "PR"; 29 29 PP.Components = ["water"]; 30 30 NComp = PP.NumberOfComponents;
Note: See TracChangeset
for help on using the changeset viewer.