Changeset 492 for trunk/eml


Ignore:
Timestamp:
Apr 4, 2008, 12:57:01 PM (14 years ago)
Author:
gerson bicca
Message:

more updates in heat exchangers

Location:
trunk/eml/heat_exchangers
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/eml/heat_exchangers/HEX_Engine.mso

    r490 r492  
    7676VARIABLES
    7777
    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}");
     78Inlet           as Properties_In_Out    (Brief="Properties at Inlet Stream", Symbol = "_{in}");
     79Average         as Properties_Average   (Brief="Properties at Average Temperature", Symbol = "_{avg}");
     80Outlet          as Properties_In_Out    (Brief="Properties at Outlet Stream", Symbol = "_{out}");
     81Wall                    as Properties_Wall              (Brief="Properties at Wall Temperature", Symbol = "_{wall}");
    8282
    8383end
     
    8787ATTRIBUTES
    8888        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       
     91VARIABLES
     92TubeFriction    as press_delta  (Brief="Tube Pressure Drop due to friction", Symbol = "\Delta P_{tube}", Default=0.01, Lower=1E-10,DisplayUnit='kPa');
     93InletNozzle     as press_delta  (Brief="Inlet Nozzle Pressure Drop", Symbol = "\Delta P_{Nozzle\_In}", Default=0.01, Lower=0,DisplayUnit='kPa');
     94OutletNozzle    as press_delta  (Brief="Outlet Nozzle Pressure Drop", Symbol = "\Delta P_{Nozzle\_Out}", Default=0.01, Lower=0,DisplayUnit='kPa');
     95Total                           as press_delta  (Brief="Total Pressure Drop", Symbol = "\Delta P_{total}", Default=0.01, Lower=1E-10,DisplayUnit='kPa');
     96Vnozzle_in              as velocity             (Brief="Inlet Nozzle Velocity", Symbol = "V_{Nozzle\_In}", Default=1, Upper=1e5, Lower=0);
     97Vnozzle_out     as velocity             (Brief="Outlet Nozzle Velocity", Symbol = "V_{Nozzle\_Out}", Default=1, Upper=1E5, Lower=0);
     98FricFactor              as fricfactor           (Brief="Friction Factor", Symbol = "f_i", Default=0.05, Lower=1e-10, Upper=2000);
    10199
    102100EQUATIONS
    103101"Total Pressure Drop"
    104         Pdtotal = PdTube + Pdnozzle_in + Pdnozzle_out;
     102        Total = TubeFriction + InletNozzle + OutletNozzle;
    105103
    106104end
     
    129127ATTRIBUTES
    130128        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       
     131VARIABLES
     132Ideal                           as press_delta          (Brief="Ideal Pressure Drop", Symbol = "\Delta P_{ideal}", Default=0.01, Lower=0,DisplayUnit='kPa');
     133CrossFlow                       as press_delta          (Brief="Cross Flow Pressure Drop", Symbol = "\Delta P_{CrossFlow}", Default=0.01, Lower=0,DisplayUnit='kPa');
     134EndZones                        as press_delta          (Brief="End Zones Pressure Drop", Symbol = "\Delta P_{EndZones}", Default=0.01, Lower=0,DisplayUnit='kPa');
     135Window                                  as press_delta          (Brief="Window Pressure Drop", Symbol = "\Delta P_{Window}", Default=0.01, Lower=1e-10,DisplayUnit='kPa');
     136InletNozzle             as press_delta          (Brief="Inlet Nozzle Pressure Drop", Symbol = "\Delta P_{Nozzle\_In}", Default=0.01, Lower=0,DisplayUnit='kPa');
     137OutletNozzle            as press_delta          (Brief="Outlet Nozzle Pressure Drop", Symbol = "\Delta P_{Nozzle\_Out}", Default=0.01, Lower=0,DisplayUnit='kPa');
     138Total                                   as press_delta          (Brief="Total Pressure Drop", Symbol = "\Delta P_{Total}", Default=0.01, Lower=0,DisplayUnit='kPa');
     139FricFactor              as fricfactor                   (Brief="Friction Factor",  Symbol = "f_i", Default=0.05, Lower=1e-10, Upper=2000);
     140Vnozzle_in                      as velocity                     (Brief="Inlet Nozzle Velocity", Symbol = "V_{Nozzle\_In}", Default=1, Upper=1e5, Lower=0);
     141Vnozzle_out             as velocity                     (Brief="Outlet Nozzle Velocity", Symbol = "V_{Nozzle\_Out}", Default=1, Upper=1e5, Lower=0);
    146142RVsquare_out            as positive                     (Brief = "Outlet Nozzle rho-V^2", Default=1, Upper=1e6, Unit = 'kg/s^2/m');
    147143RVsquare_in             as positive                     (Brief = "Inlet Nozzle rho-V^2", Default=1, Upper=1e6, Unit = 'kg/s^2/m');
     
    149145EQUATIONS
    150146"Shell Side Total Pressure Drop"
    151         Pdtotal = PdCross + PdEndZones + Pdnozzle_in + Pdnozzle_out + Pdwindow;
     147        Total = CrossFlow+ EndZones + InletNozzle + OutletNozzle + Window;
    152148
    153149end
     
    303299VARIABLES
    304300
    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");
     301PressureDrop    as Tube_Pdrop                           (Brief="Tube Side Pressure Drop", Symbol = " ");
     302HeatTransfer    as Tube_Heat_Transfer   (Brief="Tube Side Heat Transfer", Symbol = " ");
     303Properties              as Physical_Properties          (Brief="Tube Side Properties", Symbol = " ");
    308304
    309305end
     
    324320VARIABLES
    325321
    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");
     322PressureDrop    as Shell_Pdrop                                  (Brief="Shell Side Pressure Drop", Symbol = " ");
     323HeatTransfer    as Shell_Heat_Transfer          (Brief= "Shell Side Heat Transfer", Symbol = " ");
     324Properties              as Physical_Properties          (Brief="ShellSide Properties", Symbol = " ");
    329325
    330326end
  • trunk/eml/heat_exchangers/HeatExchangerDetailed.mso

    r490 r492  
    2424        Brief  = "Basic Model for Detailed Shell and Tube Heat Exchanger.";
    2525        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.";
    2751
    2852PARAMETERS
     
    3054HotSide                 as Switcher     (Brief="Hot Side in the Exchanger",Valid=["shell","tubes"],Default="shell");
    3155ShellType       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");
     56Pattern         as Switcher     (Brief="Tube Layout Characteristic Angle",Valid=["Triangle","Rotated Square","Square"],Default="Triangle");
    3457
    3558VARIABLES
    3659
    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}");       
    3962        Baffles                 as Baffles_Main                         (Brief="Baffle Section", Symbol=" ");
    4063        Clearances      as Clearances_Main      (Brief="Diametral Clearances", Symbol=" ");
    4164       
    42 in  InletTube           as stream                                       (Brief="Inlet Tube Stream", PosX=0.08, PosY=0, Symbol="_{inTube}");
    43 out OutletTube  as streamPH                             (Brief="Outlet Tube Stream", PosX=0.08, PosY=1, Symbol="_{outTube}");
    44 in  InletShell          as stream                                       (Brief="Inlet Shell Stream", PosX=0.2237, PosY=0, Symbol="_{inShell}");
    45 out OutletShell as streamPH                             (Brief="Outlet Shell Stream", PosX=0.8237, PosY=1, Symbol="_{outShell}");
    46         Details                         as Details_Main                         (Brief="Details in Heat Exchanger");
     65in  InletTube           as stream                                       (Brief="Inlet Tube Stream", PosX=0.08, PosY=0, Symbol="_{in }^{tube}");
     66out OutletTube  as streamPH                             (Brief="Outlet Tube Stream", PosX=0.08, PosY=1, Symbol="_{out }^{tube}");
     67in  InletShell          as stream                                       (Brief="Inlet Shell Stream", PosX=0.2237, PosY=0, Symbol="_{in }^{shell}");
     68out 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 = " ");
    4770       
    4871#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
     
    5780        Rspd    as Real                                 (Brief = "Pressure Drop Correction Factor for Unequal Baffle Spacing", Hidden = true, Lower=1e-3);
    5881        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);
    6790
    6891PARAMETERS
     
    83106#                               Auxiliar Parameters     - Must be hidden                                                       
    84107#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
    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);
     108Pi                                                      as constant   (Brief="Pi Number", Hidden = true, Default=3.14159265, Symbol = "\pi");
     109Aonozzle_Shell          as area                 (Brief="Shell Outlet Nozzle Area", Hidden = true, Lower=1E-6 ,  Symbol ="A_{nozzle\_out }^{shell}");
     110Ainozzle_Shell          as area                 (Brief="Shell Inlet Nozzle Area", Hidden = true, Lower=1E-6 ,  Symbol ="A_{nozzle\_in }^{shell}");
     111Aeonozzle_Shell         as area                 (Brief="Shell Outlet Escape Area Under Nozzle", Hidden = true, Lower=1E-6 ,  Symbol ="Aescape_{nozzle\_out }^{shell}");
     112Aeinozzle_Shell         as area                 (Brief="Shell Inlet Escape Area Under Nozzle",Hidden = true, Lower=1E-6 ,  Symbol ="Aescape_{nozzle\_in }^{shell}");
     113Aonozzle_Tube   as area                 (Brief="Tube Outlet Nozzle Area", Hidden = true, Lower=1E-6 ,  Symbol ="A_{nozzle\_out }^{tube}");
     114Ainozzle_Tube           as area                 (Brief="Tube Inlet Nozzle Area", Hidden = true, Lower=1E-6 ,  Symbol ="A_{nozzle\_in }^{tube}");
     115Kinlet_Tube             as positive     (Brief="Tube Inlet Nozzle Pressure Loss Coeff", Hidden = true, Default=1.1, Symbol ="K_{in }^{tube}");
     116Koutlet_Tube            as positive     (Brief="Tube Outlet Nozzle Pressure Loss Coeff",Hidden = true, Default=0.7, Symbol ="K_{out }^{tube}");
    94117Ods                                             as Real                 (Brief="Baffle cut angle in degrees", Symbol = "\theta _{ds}", Hidden = true);
    95118Octl                                    as Real                 (Brief="Baffle cut angle relative to the centerline in degrees", Symbol = "\theta _{ctl}", Hidden = true);
     
    159182       
    160183"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;
    162185       
    163186"Shell Stream Average Pressure"
     
    165188       
    166189"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;
    168191       
    169192"Shell Stream Average Molecular Weight"
     
    342365       
    343366"Friction Factor for Pressure Drop in Laminar Flow"
    344         Tubes.PressureDrop.fi   = 16/Tubes.HeatTransfer.Re;
     367        Tubes.PressureDrop.FricFactor  = 16/Tubes.HeatTransfer.Re;
    345368       
    346369switch LaminarCorrelation
     
    366389       
    367390"Friction Factor for Pressure Drop in Transition Flow"
    368         Tubes.PressureDrop.fi   = 0.0122;
     391        Tubes.PressureDrop.FricFactor  = 0.0122;
    369392       
    370393switch TransitionCorrelation
     
    391414
    392415"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);
    394417
    395418switch TurbulentCorrelation
     
    434457       
    435458"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);
    437460
    438461"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;
    440463
    441464switch ShellRange
     
    447470       
    448471"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;
    450473       
    451474        when Shell.HeatTransfer.Re > 10 switchto "range2";
     
    457480       
    458481"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;
    460483       
    461484        when Shell.HeatTransfer.Re > 100 switchto "range3";
     
    467490       
    468491"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;
    470493       
    471494        when Shell.HeatTransfer.Re > 1000 switchto "range4";
     
    477500       
    478501"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;
    480503
    481504        when Shell.HeatTransfer.Re > 10000 switchto "range5";
     
    487510
    488511"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;
    490513       
    491514        when Shell.HeatTransfer.Re < 10000 switchto "range4";
     
    514537
    515538"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);
    517540
    518541"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;
    520543
    521544switch ShellRange
     
    527550
    528551"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;
    530553       
    531554        when Shell.HeatTransfer.Re > 10 switchto "range2";
     
    537560
    538561"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;
    540563       
    541564        when Shell.HeatTransfer.Re > 100 switchto "range3";
     
    547570
    548571"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;
    550573       
    551574        when Shell.HeatTransfer.Re > 1000 switchto "range4";
     
    557580       
    558581"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;
    560583
    561584        when Shell.HeatTransfer.Re > 10000 switchto "range5";
     
    567590       
    568591"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;
    570593       
    571594        when Shell.HeatTransfer.Re < 10000 switchto "range4";
     
    594617
    595618"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);
    597620
    598621"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;
    600623
    601624switch ShellRange
     
    607630       
    608631"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;
    610633
    611634        when Shell.HeatTransfer.Re > 10 switchto "range2";
     
    617640
    618641"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;
    620643
    621644        when Shell.HeatTransfer.Re > 100 switchto "range3";
     
    627650
    628651"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;
    630653
    631654        when Shell.HeatTransfer.Re > 1000 switchto "range4";
     
    637660
    638661"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;
    640663
    641664        when Shell.HeatTransfer.Re > 10000 switchto "range5";
     
    647670
    648671"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;
    650673
    651674        when Shell.HeatTransfer.Re < 10000 switchto "range4";
     
    675698
    676699"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));
    678701
    679702        when Shell.HeatTransfer.Re > 20 switchto "laminar";
     
    697720
    698721"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));
    700723
    701724        when Shell.HeatTransfer.Re < 20 switchto "deep laminar";
     
    720743
    721744"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));
    723746       
    724747        when Shell.HeatTransfer.Re < 100 switchto "laminar";
     
    727750
    728751"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));
    730753
    731754"Shell Side Phi correction"
     
    742765
    743766"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);
    745768
    746769"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;
    748771
    749772"Velocity Tube Side Inlet Nozzle"
     
    751774
    752775"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;
    754777
    755778"Velocity Tube Side Outlet Nozzle"
     
    757780
    758781"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));
    760783
    761784"Velocity Shell Side Inlet Nozzle"
     
    763786
    764787"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));
    766789
    767790"Velocity Shell Side Outlet Nozzle"
     
    769792
    770793"Pressure Drop Shell Stream"
    771         OutletShell.P  = InletShell.P - Shell.PressureDrop.Pdtotal;     
     794        OutletShell.P  = InletShell.P - Shell.PressureDrop.Total;       
    772795       
    773796"Pressure Drop Tube Stream"
    774         OutletTube.P  = InletTube.P - Tubes.PressureDrop.Pdtotal;
     797        OutletTube.P  = InletTube.P - Tubes.PressureDrop.Total;
    775798
    776799"Shell Wall Temperature"
Note: See TracChangeset for help on using the changeset viewer.