Changeset 492
- Timestamp:
- Apr 4, 2008, 12:57:01 PM (16 years ago)
- Location:
- trunk/eml/heat_exchangers
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/eml/heat_exchangers/HEX_Engine.mso
r490 r492 76 76 VARIABLES 77 77 78 Inlet as Properties_In_Out (Brief="Properties at Inlet Stream", Symbol = "^{in}");79 Average as Properties_Average (Brief="Properties at Average Temperature", Symbol = " ^{avg}");80 Outlet as Properties_In_Out (Brief="Properties at Outlet Stream", Symbol = "^{out}");81 Wall as Properties_Wall (Brief="Properties at Wall Temperature", Symbol = "^{wall}");78 Inlet as Properties_In_Out (Brief="Properties at Inlet Stream", Symbol = "_{in}"); 79 Average as Properties_Average (Brief="Properties at Average Temperature", Symbol = "_{avg}"); 80 Outlet as Properties_In_Out (Brief="Properties at Outlet Stream", Symbol = "_{out}"); 81 Wall as Properties_Wall (Brief="Properties at Wall Temperature", Symbol = "_{wall}"); 82 82 83 83 end … … 87 87 ATTRIBUTES 88 88 Pallete = false; 89 Brief = "to be documented"; 90 Info = 91 "to be documented"; 92 93 VARIABLES 94 PdTube as press_delta (Brief="Tube Pressure Drop",Default=0.01, Lower=1e-10,DisplayUnit='kPa'); 95 Pdtotal as press_delta (Brief="Total Pressure Drop",Default=0.01, Lower=1e-10,DisplayUnit='kPa'); 96 Pdnozzle_in as press_delta (Brief="Inlet Nozzle Pressure Drop",Default=0.01, Lower=0,DisplayUnit='kPa'); 97 Pdnozzle_out as press_delta (Brief="Outlet Nozzle Pressure Drop",Default=0.01, Lower=0,DisplayUnit='kPa'); 98 Vnozzle_in as velocity (Brief="Inlet Nozzle Velocity",Default=1, Upper=1e5, Lower=0); 99 Vnozzle_out as velocity (Brief="Outlet Nozzle Velocity",Default=1, Upper=1e5, Lower=0); 100 fi as fricfactor (Brief="Friction Factor", Default=0.05, Lower=1e-10, Upper=2000); 89 Brief = "Pressure drop and velocities in the tube side section of a shell and tube heat exchanger."; 90 91 VARIABLES 92 TubeFriction as press_delta (Brief="Tube Pressure Drop due to friction", Symbol = "\Delta P_{tube}", Default=0.01, Lower=1E-10,DisplayUnit='kPa'); 93 InletNozzle as press_delta (Brief="Inlet Nozzle Pressure Drop", Symbol = "\Delta P_{Nozzle\_In}", Default=0.01, Lower=0,DisplayUnit='kPa'); 94 OutletNozzle as press_delta (Brief="Outlet Nozzle Pressure Drop", Symbol = "\Delta P_{Nozzle\_Out}", Default=0.01, Lower=0,DisplayUnit='kPa'); 95 Total as press_delta (Brief="Total Pressure Drop", Symbol = "\Delta P_{total}", Default=0.01, Lower=1E-10,DisplayUnit='kPa'); 96 Vnozzle_in as velocity (Brief="Inlet Nozzle Velocity", Symbol = "V_{Nozzle\_In}", Default=1, Upper=1e5, Lower=0); 97 Vnozzle_out as velocity (Brief="Outlet Nozzle Velocity", Symbol = "V_{Nozzle\_Out}", Default=1, Upper=1E5, Lower=0); 98 FricFactor as fricfactor (Brief="Friction Factor", Symbol = "f_i", Default=0.05, Lower=1e-10, Upper=2000); 101 99 102 100 EQUATIONS 103 101 "Total Pressure Drop" 104 Pdtotal = PdTube + Pdnozzle_in + Pdnozzle_out;102 Total = TubeFriction + InletNozzle + OutletNozzle; 105 103 106 104 end … … 129 127 ATTRIBUTES 130 128 Pallete = false; 131 Brief = "to be documented"; 132 Info = 133 "to be documented"; 134 135 VARIABLES 136 Pideal as press_delta (Brief="Ideal Pressure Drop",Default=0.01, Lower=0,DisplayUnit='kPa'); 137 PdCross as press_delta (Brief="Cross Flow Pressure Drop",Default=0.01, Lower=0,DisplayUnit='kPa'); 138 PdEndZones as press_delta (Brief="End Zones Pressure Drop",Default=0.01, Lower=0,DisplayUnit='kPa'); 139 Pdwindow as press_delta (Brief="Window Pressure Drop",Default=0.01, Lower=1e-10,DisplayUnit='kPa'); 140 Pdtotal as press_delta (Brief="Total Pressure Drop",Default=0.01, Lower=0,DisplayUnit='kPa'); 141 Pdnozzle_in as press_delta (Brief="Inlet Nozzle Pressure Drop",Default=0.01, Lower=0,DisplayUnit='kPa'); 142 Pdnozzle_out as press_delta (Brief="Outlet Nozzle Pressure Drop",Default=0.01, Lower=0,DisplayUnit='kPa'); 143 fi as fricfactor (Brief="Friction Factor", Default=0.05, Lower=1e-10, Upper=2000); 144 Vnozzle_in as velocity (Brief="Inlet Nozzle Velocity",Default=1, Upper=1e5, Lower=0); 145 Vnozzle_out as velocity (Brief="Outlet Nozzle Velocity",Default=1, Upper=1e5, Lower=0); 129 Brief = "Pressure drop and velocities in the shell side section of a shell and tube heat exchanger."; 130 131 VARIABLES 132 Ideal as press_delta (Brief="Ideal Pressure Drop", Symbol = "\Delta P_{ideal}", Default=0.01, Lower=0,DisplayUnit='kPa'); 133 CrossFlow as press_delta (Brief="Cross Flow Pressure Drop", Symbol = "\Delta P_{CrossFlow}", Default=0.01, Lower=0,DisplayUnit='kPa'); 134 EndZones as press_delta (Brief="End Zones Pressure Drop", Symbol = "\Delta P_{EndZones}", Default=0.01, Lower=0,DisplayUnit='kPa'); 135 Window as press_delta (Brief="Window Pressure Drop", Symbol = "\Delta P_{Window}", Default=0.01, Lower=1e-10,DisplayUnit='kPa'); 136 InletNozzle as press_delta (Brief="Inlet Nozzle Pressure Drop", Symbol = "\Delta P_{Nozzle\_In}", Default=0.01, Lower=0,DisplayUnit='kPa'); 137 OutletNozzle as press_delta (Brief="Outlet Nozzle Pressure Drop", Symbol = "\Delta P_{Nozzle\_Out}", Default=0.01, Lower=0,DisplayUnit='kPa'); 138 Total as press_delta (Brief="Total Pressure Drop", Symbol = "\Delta P_{Total}", Default=0.01, Lower=0,DisplayUnit='kPa'); 139 FricFactor as fricfactor (Brief="Friction Factor", Symbol = "f_i", Default=0.05, Lower=1e-10, Upper=2000); 140 Vnozzle_in as velocity (Brief="Inlet Nozzle Velocity", Symbol = "V_{Nozzle\_In}", Default=1, Upper=1e5, Lower=0); 141 Vnozzle_out as velocity (Brief="Outlet Nozzle Velocity", Symbol = "V_{Nozzle\_Out}", Default=1, Upper=1e5, Lower=0); 146 142 RVsquare_out as positive (Brief = "Outlet Nozzle rho-V^2", Default=1, Upper=1e6, Unit = 'kg/s^2/m'); 147 143 RVsquare_in as positive (Brief = "Inlet Nozzle rho-V^2", Default=1, Upper=1e6, Unit = 'kg/s^2/m'); … … 149 145 EQUATIONS 150 146 "Shell Side Total Pressure Drop" 151 Pdtotal = PdCross + PdEndZones + Pdnozzle_in + Pdnozzle_out + Pdwindow;147 Total = CrossFlow+ EndZones + InletNozzle + OutletNozzle + Window; 152 148 153 149 end … … 303 299 VARIABLES 304 300 305 PressureDrop as Tube_Pdrop (Brief="Tube Side Pressure Drop" );306 HeatTransfer as Tube_Heat_Transfer (Brief="Tube Side Heat Transfer" );307 Properties as Physical_Properties (Brief="Tube Side Properties" );301 PressureDrop as Tube_Pdrop (Brief="Tube Side Pressure Drop", Symbol = " "); 302 HeatTransfer as Tube_Heat_Transfer (Brief="Tube Side Heat Transfer", Symbol = " "); 303 Properties as Physical_Properties (Brief="Tube Side Properties", Symbol = " "); 308 304 309 305 end … … 324 320 VARIABLES 325 321 326 PressureDrop as Shell_Pdrop (Brief="Shell Side Pressure Drop" );327 HeatTransfer as Shell_Heat_Transfer (Brief= "Shell Side Heat Transfer" );328 Properties as Physical_Properties (Brief="ShellSide Properties" );322 PressureDrop as Shell_Pdrop (Brief="Shell Side Pressure Drop", Symbol = " "); 323 HeatTransfer as Shell_Heat_Transfer (Brief= "Shell Side Heat Transfer", Symbol = " "); 324 Properties as Physical_Properties (Brief="ShellSide Properties", Symbol = " "); 329 325 330 326 end -
trunk/eml/heat_exchangers/HeatExchangerDetailed.mso
r490 r492 24 24 Brief = "Basic Model for Detailed Shell and Tube Heat Exchanger."; 25 25 Info = 26 "to be documented."; 26 "to be documented. 27 28 == Assumptions == 29 * to be documented 30 31 == Specify == 32 * to be documented 33 34 == Setting Parameters == 35 * to be documented 36 37 == References == 38 [1] E.A.D. Saunders, Heat Exchangers: Selection, Design and 39 Construction, Longman, Harlow, 1988. 40 41 [2] Taborek, J., Shell-and-tube heat exchangers, in Heat Exchanger Design Handbook, Vol. 3 42 Hemisphere Publishing Corp., New York, 1988. 43 44 [3] Bell, K. J., Mueller, A. C., Wolverine Engineering Data Book II. Wolverine Tube, Inc., <www.wlv.com>, 2001. 45 46 [4] Fakheri, A. , Alternative approach for determining log mean temperature difference correction factor 47 and number of shells of shell and tube heat exchangers, Journal of Enhanced Heat Transfer, v. 10, p. 407- 420, 2003. 48 49 [5] Gnielinski, V., Forced convection in ducts, in Heat Exchanger Design Handbook, Vol. 2 50 Hemisphere Publishing Corp., New York, 1988."; 27 51 28 52 PARAMETERS … … 30 54 HotSide as Switcher (Brief="Hot Side in the Exchanger",Valid=["shell","tubes"],Default="shell"); 31 55 ShellType as Switcher (Brief="TEMA Designation",Valid=["Eshell","Fshell"],Default="Eshell"); 32 33 Pattern as Switcher (Brief="Tube Layout Characteristic Angle",Valid=["Triangle","Rotated Square","Square"],Default="Triangle"); 56 Pattern as Switcher (Brief="Tube Layout Characteristic Angle",Valid=["Triangle","Rotated Square","Square"],Default="Triangle"); 34 57 35 58 VARIABLES 36 59 37 Tubes as Tube_Side_Main (Brief="Tube Side Section" );38 Shell as Shell_Side_Main (Brief="Shell Side Section" );60 Tubes as Tube_Side_Main (Brief="Tube Side Section" , Symbol="^{tube}"); 61 Shell as Shell_Side_Main (Brief="Shell Side Section" , Symbol="^{shell}"); 39 62 Baffles as Baffles_Main (Brief="Baffle Section", Symbol=" "); 40 63 Clearances as Clearances_Main (Brief="Diametral Clearances", Symbol=" "); 41 64 42 in InletTube as stream (Brief="Inlet Tube Stream", PosX=0.08, PosY=0, Symbol="_{in Tube}");43 out OutletTube as streamPH (Brief="Outlet Tube Stream", PosX=0.08, PosY=1, Symbol="_{out Tube}");44 in InletShell as stream (Brief="Inlet Shell Stream", PosX=0.2237, PosY=0, Symbol="_{in Shell}");45 out OutletShell as streamPH (Brief="Outlet Shell Stream", PosX=0.8237, PosY=1, Symbol="_{out Shell}");46 Details as Details_Main (Brief="Details in Heat Exchanger" );65 in InletTube as stream (Brief="Inlet Tube Stream", PosX=0.08, PosY=0, Symbol="_{in }^{tube}"); 66 out OutletTube as streamPH (Brief="Outlet Tube Stream", PosX=0.08, PosY=1, Symbol="_{out }^{tube}"); 67 in InletShell as stream (Brief="Inlet Shell Stream", PosX=0.2237, PosY=0, Symbol="_{in }^{shell}"); 68 out OutletShell as streamPH (Brief="Outlet Shell Stream", PosX=0.8237, PosY=1, Symbol="_{out }^{shell}"); 69 Details as Details_Main (Brief="Details in Heat Exchanger", Symbol = " "); 47 70 48 71 #+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++# … … 57 80 Rspd as Real (Brief = "Pressure Drop Correction Factor for Unequal Baffle Spacing", Hidden = true, Lower=1e-3); 58 81 mw as Real (Brief = "Mass Velocity in Window Zone", Hidden = true, Unit='kg/m^2/s'); 59 Ji as constant (Brief="Shell Side Ji Factor", Hidden = true, Default=0.05);60 Jr as positive (Brief="Shell Side Jr Factor", Hidden = true, Lower=10e-6);61 Jl as positive (Brief="Shell Side Jl Factor", Hidden = true, Lower=10e-6);62 Jb as positive (Brief="Shell Side Jb Factor", Hidden = true, Lower=10e-6);63 Jc as positive (Brief="Shell Side Jc Factor", Hidden = true, Lower=10e-6);64 Js as positive (Brief="Shell Side Js Factor", Hidden = true, Lower=10e-6);65 Jtotal as positive (Brief="Shell Side Jtotal Factor", Hidden = true, Lower=10e-6);66 Sm as area (Brief="Shell Side Cross Flow Area", Hidden = true, Default=0.05,Lower=10e-6);82 Ji as constant (Brief="Shell Side Ji Factor", Symbol ="J_i", Hidden = true, Default=0.05); 83 Jr as positive (Brief="Shell Side Jr Factor", Symbol ="J_r", Hidden = true, Lower=10e-6); 84 Jl as positive (Brief="Shell Side Jl Factor", Symbol ="J_l", Hidden = true, Lower=10e-6); 85 Jb as positive (Brief="Shell Side Jb Factor", Symbol ="J_b", Hidden = true, Lower=10e-6); 86 Jc as positive (Brief="Shell Side Jc Factor", Symbol ="J_c", Hidden = true, Lower=10e-6); 87 Js as positive (Brief="Shell Side Js Factor", Symbol ="J_s", Hidden = true, Lower=10e-6); 88 Jtotal as positive (Brief="Shell Side Jtotal Factor", Symbol ="J_{total}", Hidden = true, Lower=10e-6); 89 Sm as area (Brief="Shell Side Cross Flow Area", Symbol ="S_m", Hidden = true, Default=0.05,Lower=10e-6); 67 90 68 91 PARAMETERS … … 83 106 # Auxiliar Parameters - Must be hidden 84 107 #+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++# 85 Pi as constant (Brief="Pi Number", Hidden = true, Default=3.14159265 );86 Aonozzle_Shell as area (Brief="Shell Outlet Nozzle Area", Hidden = true, Lower=1E-6 );87 Ainozzle_Shell as area (Brief="Shell Inlet Nozzle Area", Hidden = true, Lower=1E-6 );88 Aeonozzle_Shell as area (Brief="Shell Outlet Escape Area Under Nozzle", Hidden = true, Lower=1E-6 );89 Aeinozzle_Shell as area (Brief="Shell Inlet Escape Area Under Nozzle",Hidden = true, Lower=1E-6 );90 Aonozzle_Tube as area (Brief="Tube Outlet Nozzle Area", Hidden = true, Lower=1E-6 );91 Ainozzle_Tube as area (Brief="Tube Inlet Nozzle Area", Hidden = true, Lower=1E-6 );92 Kinlet_Tube as positive (Brief="Tube Inlet Nozzle Pressure Loss Coeff", Hidden = true, Default=1.1 );93 Koutlet_Tube as positive (Brief="Tube Outlet Nozzle Pressure Loss Coeff",Hidden = true, Default=0.7 );108 Pi as constant (Brief="Pi Number", Hidden = true, Default=3.14159265, Symbol = "\pi"); 109 Aonozzle_Shell as area (Brief="Shell Outlet Nozzle Area", Hidden = true, Lower=1E-6 , Symbol ="A_{nozzle\_out }^{shell}"); 110 Ainozzle_Shell as area (Brief="Shell Inlet Nozzle Area", Hidden = true, Lower=1E-6 , Symbol ="A_{nozzle\_in }^{shell}"); 111 Aeonozzle_Shell as area (Brief="Shell Outlet Escape Area Under Nozzle", Hidden = true, Lower=1E-6 , Symbol ="Aescape_{nozzle\_out }^{shell}"); 112 Aeinozzle_Shell as area (Brief="Shell Inlet Escape Area Under Nozzle",Hidden = true, Lower=1E-6 , Symbol ="Aescape_{nozzle\_in }^{shell}"); 113 Aonozzle_Tube as area (Brief="Tube Outlet Nozzle Area", Hidden = true, Lower=1E-6 , Symbol ="A_{nozzle\_out }^{tube}"); 114 Ainozzle_Tube as area (Brief="Tube Inlet Nozzle Area", Hidden = true, Lower=1E-6 , Symbol ="A_{nozzle\_in }^{tube}"); 115 Kinlet_Tube as positive (Brief="Tube Inlet Nozzle Pressure Loss Coeff", Hidden = true, Default=1.1, Symbol ="K_{in }^{tube}"); 116 Koutlet_Tube as positive (Brief="Tube Outlet Nozzle Pressure Loss Coeff",Hidden = true, Default=0.7, Symbol ="K_{out }^{tube}"); 94 117 Ods as Real (Brief="Baffle cut angle in degrees", Symbol = "\theta _{ds}", Hidden = true); 95 118 Octl as Real (Brief="Baffle cut angle relative to the centerline in degrees", Symbol = "\theta _{ctl}", Hidden = true); … … 159 182 160 183 "Tube Stream Average Temperature" 161 Tubes.Properties.Average.T = 0.5* OutletTube.T + 0.5*OutletTube.T;184 Tubes.Properties.Average.T = 0.5*InletTube.T + 0.5*OutletTube.T; 162 185 163 186 "Shell Stream Average Pressure" … … 165 188 166 189 "Tube Stream Average Pressure" 167 Tubes.Properties.Average.P = 0.5* OutletTube.P+0.5*OutletTube.P;190 Tubes.Properties.Average.P = 0.5*InletTube.P+0.5*OutletTube.P; 168 191 169 192 "Shell Stream Average Molecular Weight" … … 342 365 343 366 "Friction Factor for Pressure Drop in Laminar Flow" 344 Tubes.PressureDrop. fi= 16/Tubes.HeatTransfer.Re;367 Tubes.PressureDrop.FricFactor = 16/Tubes.HeatTransfer.Re; 345 368 346 369 switch LaminarCorrelation … … 366 389 367 390 "Friction Factor for Pressure Drop in Transition Flow" 368 Tubes.PressureDrop. fi= 0.0122;391 Tubes.PressureDrop.FricFactor = 0.0122; 369 392 370 393 switch TransitionCorrelation … … 391 414 392 415 "Friction Factor for Pressure Drop in Turbulent Flow" 393 Tubes.PressureDrop. fi= 0.0035 + 0.264*Tubes.HeatTransfer.Re^(-0.42);416 Tubes.PressureDrop.FricFactor = 0.0035 + 0.264*Tubes.HeatTransfer.Re^(-0.42); 394 417 395 418 switch TurbulentCorrelation … … 434 457 435 458 "Ideal Shell Side Pressure Drop" 436 Shell.PressureDrop. Pideal= 2*Shell.PressureDrop.fi*(Shell.ShellID*(1-0.02*Baffles.BaffleCut)/(0.866*Tubes.TubePitch))*(Shell.Properties.Inlet.Fw/Sm)^2/(Shell.Properties.Average.rho*Shell.HeatTransfer.Phi);459 Shell.PressureDrop.Ideal= 2*Shell.PressureDrop.FricFactor*(Shell.ShellID*(1-0.02*Baffles.BaffleCut)/(0.866*Tubes.TubePitch))*(Shell.Properties.Inlet.Fw/Sm)^2/(Shell.Properties.Average.rho*Shell.HeatTransfer.Phi); 437 460 438 461 "Shell Pressure End Zones" 439 Shell.PressureDrop. PdEndZones = Shell.PressureDrop.Pideal*(1+ (Ncw/(Shell.ShellID*(1-0.02*Baffles.BaffleCut)/(0.866*Tubes.TubePitch))))*Rb*Rspd;462 Shell.PressureDrop.EndZones = Shell.PressureDrop.Ideal*(1+ (Ncw/(Shell.ShellID*(1-0.02*Baffles.BaffleCut)/(0.866*Tubes.TubePitch))))*Rb*Rspd; 440 463 441 464 switch ShellRange … … 447 470 448 471 "Shell Side Pressure Drop Friction Factor" 449 Shell.PressureDrop. fi=48*((1.33*Tubes.TubeOD/Tubes.TubePitch)^b)*Shell.HeatTransfer.Re^-1;472 Shell.PressureDrop.FricFactor=48*((1.33*Tubes.TubeOD/Tubes.TubePitch)^b)*Shell.HeatTransfer.Re^-1; 450 473 451 474 when Shell.HeatTransfer.Re > 10 switchto "range2"; … … 457 480 458 481 "Shell Side Pressure Drop Friction Factor" 459 Shell.PressureDrop. fi=45.10*((1.33*Tubes.TubeOD/Tubes.TubePitch)^b)*Shell.HeatTransfer.Re^-0.973;482 Shell.PressureDrop.FricFactor=45.10*((1.33*Tubes.TubeOD/Tubes.TubePitch)^b)*Shell.HeatTransfer.Re^-0.973; 460 483 461 484 when Shell.HeatTransfer.Re > 100 switchto "range3"; … … 467 490 468 491 "Shell Side Pressure Drop Friction Factor" 469 Shell.PressureDrop. fi=4.570*((1.33*Tubes.TubeOD/Tubes.TubePitch)^b)*Shell.HeatTransfer.Re^-0.476;492 Shell.PressureDrop.FricFactor=4.570*((1.33*Tubes.TubeOD/Tubes.TubePitch)^b)*Shell.HeatTransfer.Re^-0.476; 470 493 471 494 when Shell.HeatTransfer.Re > 1000 switchto "range4"; … … 477 500 478 501 "Shell Side Pressure Drop Friction Factor" 479 Shell.PressureDrop. fi=0.486*((1.33*Tubes.TubeOD/Tubes.TubePitch)^b)*Shell.HeatTransfer.Re^-0.152;502 Shell.PressureDrop.FricFactor=0.486*((1.33*Tubes.TubeOD/Tubes.TubePitch)^b)*Shell.HeatTransfer.Re^-0.152; 480 503 481 504 when Shell.HeatTransfer.Re > 10000 switchto "range5"; … … 487 510 488 511 "Shell Side Pressure Drop Friction Factor" 489 Shell.PressureDrop. fi=0.372*((1.33*Tubes.TubeOD/Tubes.TubePitch)^b)*Shell.HeatTransfer.Re^-0.123;512 Shell.PressureDrop.FricFactor=0.372*((1.33*Tubes.TubeOD/Tubes.TubePitch)^b)*Shell.HeatTransfer.Re^-0.123; 490 513 491 514 when Shell.HeatTransfer.Re < 10000 switchto "range4"; … … 514 537 515 538 "Ideal Shell Side Pressure Drop" 516 Shell.PressureDrop. Pideal= 2*Shell.PressureDrop.fi*(Shell.ShellID*(1-0.02*Baffles.BaffleCut)/(0.707*Tubes.TubePitch))*(Shell.Properties.Inlet.Fw/Sm)^2/(Shell.Properties.Average.rho*Shell.HeatTransfer.Phi);539 Shell.PressureDrop.Ideal= 2*Shell.PressureDrop.FricFactor*(Shell.ShellID*(1-0.02*Baffles.BaffleCut)/(0.707*Tubes.TubePitch))*(Shell.Properties.Inlet.Fw/Sm)^2/(Shell.Properties.Average.rho*Shell.HeatTransfer.Phi); 517 540 518 541 "Shell Pressure End Zones" 519 Shell.PressureDrop. PdEndZones = Shell.PressureDrop.Pideal*(1+ (Ncw/(Shell.ShellID*(1-0.02*Baffles.BaffleCut)/(0.707*Tubes.TubePitch))))*Rb*Rspd;542 Shell.PressureDrop.EndZones = Shell.PressureDrop.Ideal*(1+ (Ncw/(Shell.ShellID*(1-0.02*Baffles.BaffleCut)/(0.707*Tubes.TubePitch))))*Rb*Rspd; 520 543 521 544 switch ShellRange … … 527 550 528 551 "Shell Side Pressure Drop Friction Factor" 529 Shell.PressureDrop. fi=32*((1.33*Tubes.TubeOD/Tubes.TubePitch)^b)*Shell.HeatTransfer.Re^-1;552 Shell.PressureDrop.FricFactor=32*((1.33*Tubes.TubeOD/Tubes.TubePitch)^b)*Shell.HeatTransfer.Re^-1; 530 553 531 554 when Shell.HeatTransfer.Re > 10 switchto "range2"; … … 537 560 538 561 "Shell Side Pressure Drop Friction Factor" 539 Shell.PressureDrop. fi=26.20*((1.33*Tubes.TubeOD/Tubes.TubePitch)^b)*Shell.HeatTransfer.Re^-0.913;562 Shell.PressureDrop.FricFactor=26.20*((1.33*Tubes.TubeOD/Tubes.TubePitch)^b)*Shell.HeatTransfer.Re^-0.913; 540 563 541 564 when Shell.HeatTransfer.Re > 100 switchto "range3"; … … 547 570 548 571 "Shell Side Pressure Drop Friction Factor" 549 Shell.PressureDrop. fi=3.50*((1.33*Tubes.TubeOD/Tubes.TubePitch)^b)*Shell.HeatTransfer.Re^-0.476;572 Shell.PressureDrop.FricFactor=3.50*((1.33*Tubes.TubeOD/Tubes.TubePitch)^b)*Shell.HeatTransfer.Re^-0.476; 550 573 551 574 when Shell.HeatTransfer.Re > 1000 switchto "range4"; … … 557 580 558 581 "Shell Side Pressure Drop Friction Factor" 559 Shell.PressureDrop. fi=0.333*((1.33*Tubes.TubeOD/Tubes.TubePitch)^b)*Shell.HeatTransfer.Re^-0.136;582 Shell.PressureDrop.FricFactor=0.333*((1.33*Tubes.TubeOD/Tubes.TubePitch)^b)*Shell.HeatTransfer.Re^-0.136; 560 583 561 584 when Shell.HeatTransfer.Re > 10000 switchto "range5"; … … 567 590 568 591 "Shell Side Pressure Drop Friction Factor" 569 Shell.PressureDrop. fi=0.303*((1.33*Tubes.TubeOD/Tubes.TubePitch)^b)*Shell.HeatTransfer.Re^-0.126;592 Shell.PressureDrop.FricFactor=0.303*((1.33*Tubes.TubeOD/Tubes.TubePitch)^b)*Shell.HeatTransfer.Re^-0.126; 570 593 571 594 when Shell.HeatTransfer.Re < 10000 switchto "range4"; … … 594 617 595 618 "Ideal Shell Side Pressure Drop" 596 Shell.PressureDrop. Pideal= 2*Shell.PressureDrop.fi*(Shell.ShellID*(1-0.02*Baffles.BaffleCut)/Tubes.TubePitch)*(Shell.Properties.Inlet.Fw/Sm)^2/(Shell.Properties.Average.rho*Shell.HeatTransfer.Phi);619 Shell.PressureDrop.Ideal= 2*Shell.PressureDrop.FricFactor*(Shell.ShellID*(1-0.02*Baffles.BaffleCut)/Tubes.TubePitch)*(Shell.Properties.Inlet.Fw/Sm)^2/(Shell.Properties.Average.rho*Shell.HeatTransfer.Phi); 597 620 598 621 "Shell Pressure End Zones" 599 Shell.PressureDrop. PdEndZones = Shell.PressureDrop.Pideal*(1+ (Ncw/(Shell.ShellID*(1-0.02*Baffles.BaffleCut)/Tubes.TubePitch)))*Rb*Rspd;622 Shell.PressureDrop.EndZones = Shell.PressureDrop.Ideal*(1+ (Ncw/(Shell.ShellID*(1-0.02*Baffles.BaffleCut)/Tubes.TubePitch)))*Rb*Rspd; 600 623 601 624 switch ShellRange … … 607 630 608 631 "Shell Side Pressure Drop Friction Factor" 609 Shell.PressureDrop. fi=35*((1.33*Tubes.TubeOD/Tubes.TubePitch)^b)*Shell.HeatTransfer.Re^-1;632 Shell.PressureDrop.FricFactor=35*((1.33*Tubes.TubeOD/Tubes.TubePitch)^b)*Shell.HeatTransfer.Re^-1; 610 633 611 634 when Shell.HeatTransfer.Re > 10 switchto "range2"; … … 617 640 618 641 "Shell Side Pressure Drop Friction Factor" 619 Shell.PressureDrop. fi=32.10*((1.33*Tubes.TubeOD/Tubes.TubePitch)^b)*Shell.HeatTransfer.Re^-0.963;642 Shell.PressureDrop.FricFactor=32.10*((1.33*Tubes.TubeOD/Tubes.TubePitch)^b)*Shell.HeatTransfer.Re^-0.963; 620 643 621 644 when Shell.HeatTransfer.Re > 100 switchto "range3"; … … 627 650 628 651 "Shell Side Pressure Drop Friction Factor" 629 Shell.PressureDrop. fi=6.090*((1.33*Tubes.TubeOD/Tubes.TubePitch)^b)*Shell.HeatTransfer.Re^-0.602;652 Shell.PressureDrop.FricFactor=6.090*((1.33*Tubes.TubeOD/Tubes.TubePitch)^b)*Shell.HeatTransfer.Re^-0.602; 630 653 631 654 when Shell.HeatTransfer.Re > 1000 switchto "range4"; … … 637 660 638 661 "Shell Side Pressure Drop Friction Factor" 639 Shell.PressureDrop. fi=0.0815*((1.33*Tubes.TubeOD/Tubes.TubePitch)^b)*Shell.HeatTransfer.Re^0.022;662 Shell.PressureDrop.FricFactor=0.0815*((1.33*Tubes.TubeOD/Tubes.TubePitch)^b)*Shell.HeatTransfer.Re^0.022; 640 663 641 664 when Shell.HeatTransfer.Re > 10000 switchto "range5"; … … 647 670 648 671 "Shell Side Pressure Drop Friction Factor" 649 Shell.PressureDrop. fi=0.391*((1.33*Tubes.TubeOD/Tubes.TubePitch)^b)*Shell.HeatTransfer.Re^-0.148;672 Shell.PressureDrop.FricFactor=0.391*((1.33*Tubes.TubeOD/Tubes.TubePitch)^b)*Shell.HeatTransfer.Re^-0.148; 650 673 651 674 when Shell.HeatTransfer.Re < 10000 switchto "range4"; … … 675 698 676 699 "Shell Pressure Drop Baffle Window" 677 Shell.PressureDrop. Pdwindow = Baffles.NumberOfBaffles*((26/Shell.Properties.Average.rho)*mw*Shell.Properties.Average.Mu*(Ncw/(Tubes.TubePitch-Tubes.TubeOD)+ Baffles.Central_Spacing/(Dw*Dw))+ 0.5*mw*mw/Shell.Properties.Average.rho)*exp(-1.33*(1+Rs)*((Scd + Std)/Sm)^(-0.15*(1+Rs) + 0.8));700 Shell.PressureDrop.Window = Baffles.NumberOfBaffles*((26/Shell.Properties.Average.rho)*mw*Shell.Properties.Average.Mu*(Ncw/(Tubes.TubePitch-Tubes.TubeOD)+ Baffles.Central_Spacing/(Dw*Dw))+ 0.5*mw*mw/Shell.Properties.Average.rho)*exp(-1.33*(1+Rs)*((Scd + Std)/Sm)^(-0.15*(1+Rs) + 0.8)); 678 701 679 702 when Shell.HeatTransfer.Re > 20 switchto "laminar"; … … 697 720 698 721 "Shell Pressure Drop Baffle Window" 699 Shell.PressureDrop. Pdwindow = Baffles.NumberOfBaffles*((26/Shell.Properties.Average.rho)*mw*Shell.Properties.Average.Mu*(Ncw/(Tubes.TubePitch-Tubes.TubeOD)+ Baffles.Central_Spacing/(Dw*Dw))+ 0.5*mw*mw/Shell.Properties.Average.rho)*exp(-1.33*(1+Rs)*((Scd + Std)/Sm)^(-0.15*(1+Rs) + 0.8));722 Shell.PressureDrop.Window = Baffles.NumberOfBaffles*((26/Shell.Properties.Average.rho)*mw*Shell.Properties.Average.Mu*(Ncw/(Tubes.TubePitch-Tubes.TubeOD)+ Baffles.Central_Spacing/(Dw*Dw))+ 0.5*mw*mw/Shell.Properties.Average.rho)*exp(-1.33*(1+Rs)*((Scd + Std)/Sm)^(-0.15*(1+Rs) + 0.8)); 700 723 701 724 when Shell.HeatTransfer.Re < 20 switchto "deep laminar"; … … 720 743 721 744 "Shell Pressure Drop Baffle Window" 722 Shell.PressureDrop. Pdwindow = Baffles.NumberOfBaffles*((2+0.6*Ncw)*0.5*mw*mw/Shell.Properties.Average.rho)*exp(-1.33*(1+Rs)*((Scd + Std)/Sm)^(-0.15*(1+Rs) + 0.8));745 Shell.PressureDrop.Window = Baffles.NumberOfBaffles*((2+0.6*Ncw)*0.5*mw*mw/Shell.Properties.Average.rho)*exp(-1.33*(1+Rs)*((Scd + Std)/Sm)^(-0.15*(1+Rs) + 0.8)); 723 746 724 747 when Shell.HeatTransfer.Re < 100 switchto "laminar"; … … 727 750 728 751 "Shell Pressure Drop Cross Flow" 729 Shell.PressureDrop. PdCross = Shell.PressureDrop.Pideal*Rb*(Baffles.NumberOfBaffles-1)*exp(-1.33*(1+Rs)*((Scd + Std)/Sm)^(-0.15*(1+Rs) + 0.8));752 Shell.PressureDrop.CrossFlow= Shell.PressureDrop.Ideal*Rb*(Baffles.NumberOfBaffles-1)*exp(-1.33*(1+Rs)*((Scd + Std)/Sm)^(-0.15*(1+Rs) + 0.8)); 730 753 731 754 "Shell Side Phi correction" … … 742 765 743 766 "Tube Side Pressure Drop" 744 Tubes.PressureDrop. PdTube = 2*Tubes.PressureDrop.fi*Tubes.TubeLength*Tubes.Properties.Average.rho*(Tubes.HeatTransfer.Vtube^2)*Tubes.Tubepasses/(Tubes.TubeID*Tubes.HeatTransfer.Phi);767 Tubes.PressureDrop.TubeFriction = 2*Tubes.PressureDrop.FricFactor*Tubes.TubeLength*Tubes.Properties.Average.rho*(Tubes.HeatTransfer.Vtube^2)*Tubes.Tubepasses/(Tubes.TubeID*Tubes.HeatTransfer.Phi); 745 768 746 769 "Pressure Drop Tube Side Inlet Nozzle" 747 Tubes.PressureDrop. Pdnozzle_in= 0.5*Kinlet_Tube*Tubes.Properties.Inlet.rho*Tubes.PressureDrop.Vnozzle_in^2;770 Tubes.PressureDrop.InletNozzle = 0.5*Kinlet_Tube*Tubes.Properties.Inlet.rho*Tubes.PressureDrop.Vnozzle_in^2; 748 771 749 772 "Velocity Tube Side Inlet Nozzle" … … 751 774 752 775 "Pressure Drop Tube Side Outlet Nozzle" 753 Tubes.PressureDrop. Pdnozzle_out= 0.5*Koutlet_Tube*Tubes.Properties.Outlet.rho*Tubes.PressureDrop.Vnozzle_out^2;776 Tubes.PressureDrop.OutletNozzle = 0.5*Koutlet_Tube*Tubes.Properties.Outlet.rho*Tubes.PressureDrop.Vnozzle_out^2; 754 777 755 778 "Velocity Tube Side Outlet Nozzle" … … 757 780 758 781 "Shell Pressure Drop Inlet Nozzle" 759 Shell.PressureDrop. Pdnozzle_in= (0.5*Shell.Properties.Inlet.Fw^2/Shell.Properties.Inlet.rho)*((1/Ainozzle_Shell^2)+(1/Aeinozzle_Shell^2));782 Shell.PressureDrop.InletNozzle = (0.5*Shell.Properties.Inlet.Fw^2/Shell.Properties.Inlet.rho)*((1/Ainozzle_Shell^2)+(1/Aeinozzle_Shell^2)); 760 783 761 784 "Velocity Shell Side Inlet Nozzle" … … 763 786 764 787 "Shell Pressure Drop Outlet Nozzle" 765 Shell.PressureDrop. Pdnozzle_out= (0.5*Shell.Properties.Outlet.Fw^2/Shell.Properties.Outlet.rho)*((1/Ainozzle_Shell^2)+(1/Aeinozzle_Shell^2));788 Shell.PressureDrop.OutletNozzle = (0.5*Shell.Properties.Outlet.Fw^2/Shell.Properties.Outlet.rho)*((1/Ainozzle_Shell^2)+(1/Aeinozzle_Shell^2)); 766 789 767 790 "Velocity Shell Side Outlet Nozzle" … … 769 792 770 793 "Pressure Drop Shell Stream" 771 OutletShell.P = InletShell.P - Shell.PressureDrop. Pdtotal;794 OutletShell.P = InletShell.P - Shell.PressureDrop.Total; 772 795 773 796 "Pressure Drop Tube Stream" 774 OutletTube.P = InletTube.P - Tubes.PressureDrop. Pdtotal;797 OutletTube.P = InletTube.P - Tubes.PressureDrop.Total; 775 798 776 799 "Shell Wall Temperature"
Note: See TracChangeset
for help on using the changeset viewer.