Changeset 197 for branches/newlanguage/eml/heat_exchangers/DoublePipe.mso
 Timestamp:
 Mar 8, 2007, 11:31:57 AM (16 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

branches/newlanguage/eml/heat_exchangers/DoublePipe.mso
r176 r197 69 69 SET 70 70 71 #"Component Molecular Weight" 71 72 M = PP.MolecularWeight(); 72 73 74 #"Pi Number" 73 75 Pi = 3.14159265; 74 76 … … 93 95 EQUATIONS 94 96 95 "Outer Stream Average Temperature"97 "Outer Stream Average Temperature" 96 98 Outer.Properties.Average.T = 0.5*InletOuter.T + 0.5*OutletOuter.T; 97 99 … … 121 123 then 122 124 123 " Heat Capacity Inner Stream"125 "Average Heat Capacity Inner Stream" 124 126 Inner.Properties.Average.Cp = PP.LiquidCp(Inner.Properties.Average.T,Inner.Properties.Average.P,InletInner.z); 125 Inner.Properties.Inlet.Cp = PP.LiquidCp(InletInner.T,InletInner.P,InletInner.z); 127 128 "Inlet Heat Capacity Inner Stream" 129 Inner.Properties.Inlet.Cp = PP.LiquidCp(InletInner.T,InletInner.P,InletInner.z); 130 131 "Outlet Heat Capacity Inner Stream" 126 132 Inner.Properties.Outlet.Cp = PP.LiquidCp(OutletInner.T,OutletInner.P,OutletInner.z); 127 133 128 " Mass Density Inner Stream"134 "Average Mass Density Inner Stream" 129 135 Inner.Properties.Average.rho = PP.LiquidDensity(Inner.Properties.Average.T,Inner.Properties.Average.P,InletInner.z); 136 137 "Inlet Mass Density Inner Stream" 130 138 Inner.Properties.Inlet.rho = PP.LiquidDensity(InletInner.T,InletInner.P,InletInner.z); 131 Inner.Properties.Outlet.rho = PP.LiquidDensity(OutletInner.T,OutletInner.P,OutletInner.z); 132 133 "Viscosity Inner Stream" 139 140 "Outlet Mass Density Inner Stream" 141 Inner.Properties.Outlet.rho = PP.LiquidDensity(OutletInner.T,OutletInner.P,OutletInner.z); 142 143 "Average Viscosity Inner Stream" 134 144 Inner.Properties.Average.Mu = PP.LiquidViscosity(Inner.Properties.Average.T,Inner.Properties.Average.P,InletInner.z); 145 146 "Inlet Viscosity Inner Stream" 135 147 Inner.Properties.Inlet.Mu = PP.LiquidViscosity(InletInner.T,InletInner.P,InletInner.z); 136 Inner.Properties.Outlet.Mu = PP.LiquidViscosity(OutletInner.T,OutletInner.P,OutletInner.z); 137 138 "Conductivity Inner Stream" 148 149 "Outlet Viscosity Inner Stream" 150 Inner.Properties.Outlet.Mu = PP.LiquidViscosity(OutletInner.T,OutletInner.P,OutletInner.z); 151 152 "Average Conductivity Inner Stream" 139 153 Inner.Properties.Average.K = PP.LiquidThermalConductivity(Inner.Properties.Average.T,Inner.Properties.Average.P,InletInner.z); 154 155 "Inlet Conductivity Inner Stream" 140 156 Inner.Properties.Inlet.K = PP.LiquidThermalConductivity(InletInner.T,InletInner.P,InletInner.z); 157 158 "Outlet Conductivity Inner Stream" 141 159 Inner.Properties.Outlet.K = PP.LiquidThermalConductivity(OutletInner.T,OutletInner.P,OutletInner.z); 142 160 143 "Viscosity Cold Stream"161 "Viscosity Inner Stream at wall temperature" 144 162 Inner.Properties.Wall.Mu = PP.LiquidViscosity(Inner.Properties.Wall.Twall,Inner.Properties.Average.P,InletInner.z); 145 163 146 147 164 else 148 165 149 " Heat Capacity ColdStream"166 "Average Heat Capacity InnerStream" 150 167 Inner.Properties.Average.Cp = PP.VapourCp(Inner.Properties.Average.T,Inner.Properties.Average.P,InletInner.z); 151 Inner.Properties.Inlet.Cp = PP.VapourCp(InletInner.T,InletInner.P,InletInner.z); 168 169 "Inlet Heat Capacity Inner Stream" 170 Inner.Properties.Inlet.Cp = PP.VapourCp(InletInner.T,InletInner.P,InletInner.z); 171 172 "Outlet Heat Capacity Inner Stream" 152 173 Inner.Properties.Outlet.Cp = PP.VapourCp(OutletInner.T,OutletInner.P,OutletInner.z); 153 174 154 " Mass Density ColdStream"175 "Average Mass Density Inner Stream" 155 176 Inner.Properties.Average.rho = PP.VapourDensity(Inner.Properties.Average.T,Inner.Properties.Average.P,InletInner.z); 177 178 "Inlet Mass Density Inner Stream" 156 179 Inner.Properties.Inlet.rho = PP.VapourDensity(InletInner.T,InletInner.P,InletInner.z); 157 Inner.Properties.Outlet.rho = PP.VapourDensity(OutletInner.T,OutletInner.P,OutletInner.z); 158 159 "Viscosity Cold Stream" 160 Inner.Properties.Average.Mu = PP.VapourViscosity(Inner.Properties.Average.T,Inner.Properties.Average.P,InletInner.z); 180 181 "Outlet Mass Density Inner Stream" 182 Inner.Properties.Outlet.rho = PP.VapourDensity(OutletInner.T,OutletInner.P,OutletInner.z); 183 184 "Average Viscosity Inner Stream" 185 Inner.Properties.Average.Mu = PP.VapourViscosity(Inner.Properties.Average.T,Inner.Properties.Average.P,InletInner.z); 186 187 "Inlet Viscosity Inner Stream" 161 188 Inner.Properties.Inlet.Mu = PP.VapourViscosity(InletInner.T,InletInner.P,InletInner.z); 162 Inner.Properties.Outlet.Mu = PP.VapourViscosity(OutletInner.T,OutletInner.P,OutletInner.z); 163 164 "Conductivity Cold Stream" 189 190 "Outlet Viscosity Inner Stream" 191 Inner.Properties.Outlet.Mu = PP.VapourViscosity(OutletInner.T,OutletInner.P,OutletInner.z); 192 193 "Average Conductivity Inner Stream" 165 194 Inner.Properties.Average.K = PP.VapourThermalConductivity(Inner.Properties.Average.T,Inner.Properties.Average.P,InletInner.z); 195 196 "Inlet Conductivity Inner Stream" 166 197 Inner.Properties.Inlet.K = PP.VapourThermalConductivity(InletInner.T,InletInner.P,InletInner.z); 198 199 "Outlet Conductivity Inner Stream" 167 200 Inner.Properties.Outlet.K = PP.VapourThermalConductivity(OutletInner.T,OutletInner.P,OutletInner.z); 168 201 169 "Viscosity Cold Stream"202 "Viscosity Inner Stream at wall temperature" 170 203 Inner.Properties.Wall.Mu = PP.VapourViscosity(Inner.Properties.Wall.Twall,Inner.Properties.Average.P,InletInner.z); 171 204 … … 176 209 then 177 210 178 " Heat Capacity HotStream"211 "Average Heat Capacity Outer Stream" 179 212 Outer.Properties.Average.Cp = PP.LiquidCp(Outer.Properties.Average.T,Outer.Properties.Average.P,InletOuter.z); 180 Outer.Properties.Inlet.Cp = PP.LiquidCp(InletOuter.T,InletOuter.P,InletOuter.z); 213 214 "Inlet Heat Capacity Outer Stream" 215 Outer.Properties.Inlet.Cp = PP.LiquidCp(InletOuter.T,InletOuter.P,InletOuter.z); 216 217 "Outlet Heat Capacity Outer Stream" 181 218 Outer.Properties.Outlet.Cp = PP.LiquidCp(OutletOuter.T,OutletOuter.P,OutletOuter.z); 182 219 183 "Mass Density Hot Stream" 184 Outer.Properties.Average.rho = PP.LiquidDensity(Outer.Properties.Average.T,Outer.Properties.Average.P,InletOuter.z); 185 Outer.Properties.Inlet.rho = PP.LiquidDensity(InletOuter.T,InletOuter.P,InletOuter.z); 220 "Average Mass Density Outer Stream" 221 Outer.Properties.Average.rho = PP.LiquidDensity(Outer.Properties.Average.T,Outer.Properties.Average.P,InletOuter.z); 222 223 "Inlet Mass Density Outer Stream" 224 Outer.Properties.Inlet.rho = PP.LiquidDensity(InletOuter.T,InletOuter.P,InletOuter.z); 225 226 "Outlet Mass Density Outer Stream" 186 227 Outer.Properties.Outlet.rho = PP.LiquidDensity(OutletOuter.T,OutletOuter.P,OutletOuter.z); 187 228 188 " Viscosity HotStream"229 "Average Viscosity Outer Stream" 189 230 Outer.Properties.Average.Mu = PP.LiquidViscosity(Outer.Properties.Average.T,Outer.Properties.Average.P,InletOuter.z); 231 232 "Inlet Viscosity Outer Stream" 190 233 Outer.Properties.Inlet.Mu = PP.LiquidViscosity(InletOuter.T,InletOuter.P,InletOuter.z); 234 235 "Outlet Viscosity Outer Stream" 191 236 Outer.Properties.Outlet.Mu = PP.LiquidViscosity(OutletOuter.T,OutletOuter.P,OutletOuter.z); 192 237 193 " Conductivity HotStream"238 "Average Conductivity Outer Stream" 194 239 Outer.Properties.Average.K = PP.LiquidThermalConductivity(Outer.Properties.Average.T,Outer.Properties.Average.P,InletOuter.z); 240 241 "Inlet Conductivity Outer Stream" 195 242 Outer.Properties.Inlet.K = PP.LiquidThermalConductivity(InletOuter.T,InletOuter.P,InletOuter.z); 243 244 "Outlet Conductivity Outer Stream" 196 245 Outer.Properties.Outlet.K = PP.LiquidThermalConductivity(OutletOuter.T,OutletOuter.P,OutletOuter.z); 197 246 198 "Viscosity Hot Stream"247 "Viscosity Outer Stream at wall temperature" 199 248 Outer.Properties.Wall.Mu = PP.LiquidViscosity(Outer.Properties.Wall.Twall,Outer.Properties.Average.P,InletOuter.z); 200 249 … … 202 251 else 203 252 204 " Heat Capacity HotStream"253 "Average Heat Capacity Outer Stream" 205 254 Outer.Properties.Average.Cp = PP.VapourCp(Outer.Properties.Average.T,Outer.Properties.Average.P,InletOuter.z); 206 Outer.Properties.Inlet.Cp = PP.VapourCp(InletOuter.T,InletOuter.P,InletOuter.z); 255 256 "Inlet Heat Capacity Outer Stream" 257 Outer.Properties.Inlet.Cp = PP.VapourCp(InletOuter.T,InletOuter.P,InletOuter.z); 258 259 "Outlet Heat Capacity Outer Stream" 207 260 Outer.Properties.Outlet.Cp = PP.VapourCp(OutletOuter.T,OutletOuter.P,OutletOuter.z); 208 261 209 "Mass Density Hot Stream" 210 Outer.Properties.Average.rho = PP.VapourDensity(Outer.Properties.Average.T,Outer.Properties.Average.P,InletOuter.z); 211 Outer.Properties.Inlet.rho = PP.VapourDensity(InletOuter.T,InletOuter.P,InletOuter.z); 262 "Average Mass Density Outer Stream" 263 Outer.Properties.Average.rho = PP.VapourDensity(Outer.Properties.Average.T,Outer.Properties.Average.P,InletOuter.z); 264 265 "Inlet Mass Density Outer Stream" 266 Outer.Properties.Inlet.rho = PP.VapourDensity(InletOuter.T,InletOuter.P,InletOuter.z); 267 268 "Outlet Mass Density Outer Stream" 212 269 Outer.Properties.Outlet.rho = PP.VapourDensity(OutletOuter.T,OutletOuter.P,OutletOuter.z); 213 270 214 " Viscosity HotStream"271 "Average Viscosity Outer Stream" 215 272 Outer.Properties.Average.Mu = PP.VapourViscosity(Outer.Properties.Average.T,Outer.Properties.Average.P,InletOuter.z); 273 274 "Inlet Viscosity Outer Stream" 216 275 Outer.Properties.Inlet.Mu = PP.VapourViscosity(InletOuter.T,InletOuter.P,InletOuter.z); 276 277 "Outlet Viscosity Outer Stream" 217 278 Outer.Properties.Outlet.Mu = PP.VapourViscosity(OutletOuter.T,OutletOuter.P,OutletOuter.z); 218 279 219 " Conductivity HotStream"280 "Average Conductivity Outer Stream" 220 281 Outer.Properties.Average.K = PP.VapourThermalConductivity(Outer.Properties.Average.T,Outer.Properties.Average.P,InletOuter.z); 282 283 "Inlet Conductivity Outer Stream" 221 284 Outer.Properties.Inlet.K = PP.VapourThermalConductivity(InletOuter.T,InletOuter.P,InletOuter.z); 285 286 "Outlet Conductivity Outer Stream" 222 287 Outer.Properties.Outlet.K = PP.VapourThermalConductivity(OutletOuter.T,OutletOuter.P,OutletOuter.z); 223 288 224 "Viscosity Hot Stream"289 "Viscosity Outer Stream at wall temperature" 225 290 Outer.Properties.Wall.Mu = PP.VapourViscosity(Outer.Properties.Wall.Twall,Outer.Properties.Average.P,InletOuter.z); 226 291 … … 251 316 end 252 317 253 "Flow Mass Inlet ColdStream"318 "Flow Mass Inlet Inner Stream" 254 319 Inner.Properties.Inlet.Fw = sum(M*InletInner.z)*InletInner.F; 255 320 256 "Flow Mass Outlet ColdStream"321 "Flow Mass Outlet Inner Stream" 257 322 Inner.Properties.Outlet.Fw = sum(M*OutletInner.z)*OutletInner.F; 258 323 259 "Flow Mass Inlet HotStream"324 "Flow Mass Inlet Outer Stream" 260 325 Outer.Properties.Inlet.Fw = sum(M*InletOuter.z)*InletOuter.F; 261 326 262 "Flow Mass Outlet HotStream"327 "Flow Mass Outlet Outer Stream" 263 328 Outer.Properties.Outlet.Fw = sum(M*OutletOuter.z)*OutletOuter.F; 264 329 265 "Molar Balance HotStream"330 "Molar Balance Outer Stream" 266 331 OutletOuter.F = InletOuter.F; 267 332 268 "Molar Balance ColdStream"333 "Molar Balance Inner Stream" 269 334 OutletInner.F = InletInner.F; 270 335 271 " HotStream Molar Fraction Constraint"336 "Outer Stream Molar Fraction Constraint" 272 337 OutletOuter.z=InletOuter.z; 273 338 274 " ColdStream Molar Fraction Constraint"339 "InnerStream Molar Fraction Constraint" 275 340 OutletInner.z=InletInner.z; 276 341 … … 340 405 341 406 case "Hausen": 407 408 "Nusselt Number" 342 409 Inner.HeatTransfer.Nu = 3.665 + ((0.19*((DiInner/Lpipe)*Inner.HeatTransfer.Re*Inner.HeatTransfer.PR)^0.8)/(1+0.117*((DiInner/Lpipe)*Inner.HeatTransfer.Re*Inner.HeatTransfer.PR)^0.467)); 343 410 344 411 case "Schlunder": 412 413 "Nusselt Number" 345 414 Inner.HeatTransfer.Nu = (49.027896+4.173281*Inner.HeatTransfer.Re*Inner.HeatTransfer.PR*(DiInner/Lpipe))^(1/3); 346 415 … … 351 420 case "transition": 352 421 422 "Inner Side Friction Factor for Heat Transfer  transition Flow" 353 423 Inner.HeatTransfer.fi = 1/(0.79*ln(Inner.HeatTransfer.Re)1.64)^2; 354 424 … … 356 426 357 427 case "Gnielinski": 428 429 "Nusselt Number" 358 430 Inner.HeatTransfer.Nu*(1+(12.7*sqrt(0.125*Inner.HeatTransfer.fi)*((Inner.HeatTransfer.PR)^(2/3) 1))) = 0.125*Inner.HeatTransfer.fi*(Inner.HeatTransfer.Re1000)*Inner.HeatTransfer.PR; 359 431 360 432 case "ESDU": 433 434 "Nusselt Number" 361 435 Inner.HeatTransfer.Nu =1;#to be implemented 362 436 … … 371 445 372 446 case "Petukhov": 447 448 "Inner Side Friction Factor for Heat Transfer  turbulent Flow" 373 449 Inner.HeatTransfer.fi = 1/(1.82*log(Inner.HeatTransfer.Re)1.64)^2; 450 451 "Nusselt Number" 374 452 Inner.HeatTransfer.Nu*(1.07+(12.7*sqrt(0.125*Inner.HeatTransfer.fi)*((Inner.HeatTransfer.PR)^(2/3) 1))) = 0.125*Inner.HeatTransfer.fi*Inner.HeatTransfer.Re*Inner.HeatTransfer.PR; 375 453 376 454 case "SiederTate": 455 456 "Nusselt Number" 377 457 Inner.HeatTransfer.Nu = 0.027*(Inner.HeatTransfer.PR)^(1/3)*(Inner.HeatTransfer.Re)^(4/5); 458 459 "Inner Side Friction Factor for Heat Transfer  turbulent Flow" 378 460 Inner.HeatTransfer.fi = 1/(1.82*log(Inner.HeatTransfer.Re)1.64)^2; 379 461 … … 388 470 case "laminar": 389 471 472 "Outer Side Friction Factor for Heat Transfer  laminar Flow" 390 473 Outer.HeatTransfer.fi = 1/(0.79*ln(Outer.HeatTransfer.Re)1.64)^2; 474 391 475 switch OuterLaminarCorrelation 392 476 393 477 case "Hausen": 478 479 "Nusselt Number" 394 480 Outer.HeatTransfer.Nu = 3.665 + ((0.19*((Outer.HeatTransfer.Dh/Lpipe)*Outer.HeatTransfer.Re*Outer.HeatTransfer.PR)^0.8)/(1+0.117*((Outer.HeatTransfer.Dh/Lpipe)*Outer.HeatTransfer.Re*Outer.HeatTransfer.PR)^0.467)); 395 481 396 482 case "Schlunder": 483 484 "Nusselt Number" 397 485 Outer.HeatTransfer.Nu = (49.027896+4.173281*Outer.HeatTransfer.Re*Outer.HeatTransfer.PR*(Outer.HeatTransfer.Dh/Lpipe))^(1/3); 398 486 … … 406 494 407 495 case "Gnielinski": 496 497 "Outer Side Friction Factor for Heat Transfer  transition Flow" 408 498 Outer.HeatTransfer.fi = 1/(0.79*ln(Outer.HeatTransfer.Re)1.64)^2; 499 500 "Nusselt Number" 409 501 Outer.HeatTransfer.Nu*(1+(12.7*sqrt(0.125*Outer.HeatTransfer.fi)*((Outer.HeatTransfer.PR)^(2/3) 1))) = 0.125*Outer.HeatTransfer.fi*(Outer.HeatTransfer.Re1000)*Outer.HeatTransfer.PR; 410 502 411 503 case "ESDU": 504 505 "Nusselt Number" 412 506 Outer.HeatTransfer.Nu =1;#to be implemented 507 508 "Outer Side Friction Factor for Heat Transfer  transition Flow" 413 509 Outer.HeatTransfer.fi = 1/(0.79*ln(Outer.HeatTransfer.Re)1.64)^2; 414 510 … … 423 519 424 520 case "Petukhov": 521 522 "Outer Side Friction Factor for Heat Transfer  turbulent Flow" 425 523 Outer.HeatTransfer.fi = 1/(1.82*log(Outer.HeatTransfer.Re)1.64)^2; 524 525 "Nusselt Number" 426 526 Outer.HeatTransfer.Nu*(1.07+(12.7*sqrt(0.125*Outer.HeatTransfer.fi)*((Outer.HeatTransfer.PR)^(2/3) 1))) = 0.125*Outer.HeatTransfer.fi*Outer.HeatTransfer.Re*Outer.HeatTransfer.PR; 427 527 428 528 case "SiederTate": 529 530 "Nusselt Number" 429 531 Outer.HeatTransfer.Nu = 0.027*(Outer.HeatTransfer.PR)^(1/3)*(Outer.HeatTransfer.Re)^(4/5); 532 533 "Outer Side Friction Factor for Heat Transfer  turbulent Flow" 430 534 Outer.HeatTransfer.fi = 1/(1.82*log(Outer.HeatTransfer.Re)1.64)^2; 431 535 … … 442 546 Outer.HeatTransfer.hcoeff= (Outer.HeatTransfer.Nu*Outer.Properties.Average.K/Outer.HeatTransfer.Dh)*Outer.HeatTransfer.Phi; 443 547 444 "Pressure Drop HotStream"548 "Pressure Drop Outer Stream" 445 549 OutletOuter.P = InletOuter.P  Outer.PressureDrop.Pdrop; 446 550 447 "Pressure Drop ColdStream"551 "Pressure Drop Inner Stream" 448 552 OutletInner.P = InletInner.P  Inner.PressureDrop.Pdrop; 449 553 450 "Outer Pipe Pressure Drop "554 "Outer Pipe Pressure Drop for friction" 451 555 Outer.PressureDrop.Pdrop = (2*Outer.PressureDrop.fi*Lpipe*Outer.Properties.Average.rho*Outer.HeatTransfer.Vmean^2)/(Outer.PressureDrop.Dh*Outer.HeatTransfer.Phi); 452 556 453 "Inner Pipe Pressure Drop "557 "Inner Pipe Pressure Drop for friction" 454 558 Inner.PressureDrop.Pdrop = (2*Inner.PressureDrop.fi*Lpipe*Inner.Properties.Average.rho*Inner.HeatTransfer.Vmean^2)/(DiInner*Inner.HeatTransfer.Phi); 455 559 … … 507 611 VARIABLES 508 612 509 Method as NTU_Basic ;613 Method as NTU_Basic (Brief="NTU Method of Calculation"); 510 614 511 615 EQUATIONS … … 546 650 547 651 then 652 548 653 "Effectiveness in Counter Flow" 549 654 Method.Eft = Method.NTU/(1+Method.NTU); 550 655 551 656 else 657 552 658 "Effectiveness in Counter Flow" 553 659 Method.Eft = (1exp(Method.NTU*(1Method.Cr)))/(1Method.Cr*exp(Method.NTU*(1Method.Cr))); … … 606 712 VARIABLES 607 713 608 Method as LMTD_Basic ;714 Method as LMTD_Basic (Brief="LMTD Method of Calculation"); 609 715 610 716 EQUATIONS
Note: See TracChangeset
for help on using the changeset viewer.