Ignore:
Timestamp:
Jan 12, 2007, 4:44:02 PM (16 years ago)
Author:
gerson bicca
Message:

modified some equations

File:
1 edited

Legend:

Unmodified
Added
Removed

  • mso/eml/heat_exchangers/HeatExchangerDiscretized.mso

    r100 r110  
    454454Donozzle_Shell  as length               (Brief="Shell Outlet Nozzle Diameter",Lower=10e-6);
    455455Dinozzle_Shell  as length               (Brief="Shell Inlet Nozzle Diameter",Lower=10e-6);
     456Aonozzle_Shell  as area                 (Brief="Shell Outlet Nozzle Area",Lower=10e-6);
     457Ainozzle_Shell  as area                 (Brief="Shell Inlet Nozzle Area",Lower=10e-6);
     458Aeonozzle_Shell as area                 (Brief="Shell Outlet Escape Area Under Nozzle",Lower=10e-6);
     459Aeinozzle_Shell as area                 (Brief="Shell Inlet Escape Area Under Nozzle",Lower=10e-6);
    456460Hinozzle_Shell  as length               (Brief="Height Under Shell Inlet Nozzle",Lower=10e-6);
    457461Honozzle_Shell  as length               (Brief="Height Under Shell Outlet Nozzle",Lower=10e-6);
    458462Lcf                     as length               (Brief="Bundle-to-Shell Clearance",Lower=10e-8);
    459 #+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
    460 #       Tubes Geometrical Parameters
    461 #+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
    462 Ntt                             as Integer              (Brief="Total Number of Tubes in Shell",Default=100,Lower=1);
    463 Pattern         as Integer              (Brief="Tube Layout Characteristic Angle",Lower=30);
    464 Ltube                   as length               (Brief="Effective Tube Length",Lower=0.1);
    465 pitch                   as length               (Brief="Tube Pitch",Lower=1e-8);
    466 Kwall                   as conductivity (Brief="Tube Wall Material Thermal Conductivity");
    467 Dotube                  as length               (Brief="Tube Outside Diameter",Lower=10e-6);
    468 Ditube                  as length               (Brief="Tube Inside Diameter",Lower=10e-6);
    469 Donozzle_Tube   as length               (Brief="Tube Outlet Nozzle Diameter",Lower=10e-6);
    470 Dinozzle_Tube   as length               (Brief="Tube Inlett Nozzle Diameter",Lower=10e-6);
     463
     464#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
     465#                               Tubes Geometrical Parameters                                            #
     466#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
     467Ntt                       as Integer            (Brief="Total Number of Tubes in Shell",Default=100,Lower=1);
     468Pattern       as Integer                (Brief="Tube Layout Characteristic Angle",Lower=30);
     469Ltube             as length                     (Brief="Effective Tube Length",Lower=0.1);
     470pitch             as length                     (Brief="Tube Pitch",Lower=1e-8);
     471Kwall             as conductivity       (Brief="Tube Wall Material Thermal Conductivity");
     472Dotube            as length                     (Brief="Tube Outside Diameter",Lower=10e-6);
     473Ditube            as length                     (Brief="Tube Inside Diameter",Lower=10e-6);
     474Donozzle_Tube as length                 (Brief="Tube Outlet Nozzle Diameter",Lower=10e-6);
     475Dinozzle_Tube as length                 (Brief="Tube Inlet Nozzle Diameter",Lower=10e-6);
     476Aonozzle_Tube as area                   (Brief="Tube Outlet Nozzle Area",Lower=10e-6);
     477Ainozzle_Tube as area                   (Brief="Tube Inlet Nozzle Area",Lower=10e-6);
     478Kinlet_Tube   as positive               (Brief="Tube Inlet Nozzle Pressure Loss Coeff",Default=1.1);
     479Koutlet_Tube  as positive               (Brief="Tube Outlet Nozzle Pressure Loss Coeff",Default=0.7);
     480
    471481#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
    472482#       Baffles Geometrical Parameters
    … …  
    572582
    573583"Pressure Drop Tube Side Inlet Nozzle"
    574         Unity(1).Tubes.PressureDrop.Pdnozzle_in = HE.DeltaPtubeNozzlein(Unity(1).Properties.Cold.Inlet.rho,Unity(1).Properties.Cold.Inlet.Fw);
     584        Unity(1).Tubes.PressureDrop.Pdnozzle_in = 0.5*Kinlet_Tube*Unity(1).Properties.Cold.Inlet.rho*Unity(1).Tubes.PressureDrop.Vnozzle_in^2;
    575585
    576586"Velocity Tube Side Inlet Nozzle"
    577         Unity(1).Tubes.PressureDrop.Vnozzle_in  = HE.TubeVelocityNozzlein(Unity(1).Properties.Cold.Inlet.rho,Unity(1).Properties.Cold.Inlet.Fw);
     587        Unity(1).Tubes.PressureDrop.Vnozzle_in  = Unity(1).Properties.Cold.Inlet.Fw/(Unity(1).Properties.Cold.Inlet.rho*Ainozzle_Tube);
    578588
    579589"Pressure Drop Tube Side Outlet Nozzle"
    580         Unity(Nb+1).Tubes.PressureDrop.Pdnozzle_out     = HE.DeltaPtubeNozzleout(Unity(Nb+1).Properties.Cold.Outlet.rho,Unity(Nb+1).Properties.Cold.Inlet.Fw);
     590        Unity(Nb+1).Tubes.PressureDrop.Pdnozzle_out     = 0.5*Koutlet_Tube*Unity(Nb+1).Properties.Cold.Outlet.rho*Unity(Nb+1).Tubes.PressureDrop.Vnozzle_out^2;
    581591
    582592"Velocity Tube Side Outlet Nozzle"
    583         Unity(Nb+1).Tubes.PressureDrop.Vnozzle_out      = HE.TubeVelocityNozzleout(Unity(Nb+1).Properties.Cold.Outlet.rho,Unity(Nb+1).Properties.Cold.Inlet.Fw);
     593        Unity(Nb+1).Tubes.PressureDrop.Vnozzle_out      = Unity(Nb+1).Properties.Cold.Inlet.Fw/(Unity(Nb+1).Properties.Cold.Outlet.rho*Aonozzle_Tube);
    584594
    585595"Shell Pressure Drop Inlet Nozzle"
    586         Unity(1).Shell.PressureDrop.Pdnozzle_in = HE.DeltaPshellNozzleIn(Unity(1).Properties.Hot.Inlet.rho,Unity(1).Properties.Hot.Inlet.Fw);
     596        Unity(1).Shell.PressureDrop.Pdnozzle_in = (0.5*Unity(1).Properties.Hot.Inlet.Fw^2/Unity(1).Properties.Hot.Inlet.rho)*((1/Ainozzle_Shell^2)+(1/Aeinozzle_Shell^2));
    587597
    588598"Velocity Shell Side Inlet Nozzle"
    589         Unity(1).Shell.PressureDrop.Vnozzle_in  = HE.ShellVelocityNozzleIn(Unity(1).Properties.Hot.Inlet.rho,Unity(1).Properties.Hot.Inlet.Fw);
     599        Unity(1).Shell.PressureDrop.Vnozzle_in  = Unity(1).Properties.Hot.Inlet.Fw/(Unity(1).Properties.Hot.Inlet.rho*Ainozzle_Shell);
    590600
    591601"Shell Pressure Drop Outlet Nozzle"
    592         Unity(Nb+1).Shell.PressureDrop.Pdnozzle_out     =HE.DeltaPshellNozzleOut(Unity(Nb+1).Properties.Hot.Outlet.rho,Unity(Nb+1).Properties.Hot.Inlet.Fw);
     602        Unity(Nb+1).Shell.PressureDrop.Pdnozzle_out     = (0.5*Unity(Nb+1).Properties.Hot.Outlet.Fw^2/Unity(Nb+1).Properties.Hot.Outlet.rho)*((1/Ainozzle_Shell^2)+(1/Aeinozzle_Shell^2));
    593603
    594604"Velocity Shell Side Outlet Nozzle"
    595         Unity(Nb+1).Shell.PressureDrop.Vnozzle_out      =HE.ShellVelocityNozzleOut(Unity(Nb+1).Properties.Hot.Outlet.rho,Unity(Nb+1).Properties.Hot.Inlet.Fw);
    596        
     605        Unity(Nb+1).Shell.PressureDrop.Vnozzle_out      = Unity(Nb+1).Properties.Hot.Outlet.Fw/(Unity(Nb+1).Properties.Hot.Outlet.rho*Aonozzle_Shell);
     606
    597607
    598608        else
    … …  
    612622        Unity(Nb+1).Baffles.Lsi,Unity(Nb+1).Properties.Cold.Inlet.Fw,Unity(Nb+1).Shell.HeatTransfer.Phi,Unity(Nb+1).Properties.Cold.Average.rho);
    613623
    614 
    615624"Pressure Drop Tube Side Inlet Nozzle"
    616         Unity(1).Tubes.PressureDrop.Pdnozzle_in = HE.DeltaPtubeNozzlein(Unity(1).Properties.Hot.Inlet.rho,Unity(1).Properties.Hot.Inlet.Fw);
    617 
     625        Unity(1).Tubes.PressureDrop.Pdnozzle_in = 0.5*Kinlet_Tube*Unity(1).Properties.Hot.Inlet.rho*Unity(1).Tubes.PressureDrop.Vnozzle_in^2;
     626       
    618627"Velocity Tube Side Inlet Nozzle"
    619         Unity(1).Tubes.PressureDrop.Vnozzle_in  = HE.TubeVelocityNozzlein(Unity(1).Properties.Hot.Inlet.rho,Unity(1).Properties.Hot.Inlet.Fw);
    620 
     628        Unity(1).Tubes.PressureDrop.Vnozzle_in  = Unity(1).Properties.Hot.Inlet.Fw/(Unity(1).Properties.Hot.Inlet.rho*Ainozzle_Tube);
     629       
    621630"Pressure Drop Tube Side Outlet Nozzle"
    622         Unity(Nb+1).Tubes.PressureDrop.Pdnozzle_out     = HE.DeltaPtubeNozzleout(Unity(Nb+1).Properties.Hot.Outlet.rho,Unity(Nb+1).Properties.Hot.Inlet.Fw);   
    623 
     631        Unity(Nb+1).Tubes.PressureDrop.Pdnozzle_out     = 0.5*Koutlet_Tube*Unity(Nb+1).Properties.Hot.Outlet.rho*Unity(Nb+1).Tubes.PressureDrop.Vnozzle_out^2;
     632       
    624633"Velocity Tube Side Outlet Nozzle"
    625         Unity(Nb+1).Tubes.PressureDrop.Vnozzle_out      = HE.TubeVelocityNozzleout(Unity(Nb+1).Properties.Hot.Outlet.rho,Unity(Nb+1).Properties.Hot.Inlet.Fw); 
     634        Unity(Nb+1).Tubes.PressureDrop.Vnozzle_out      = Unity(Nb+1).Properties.Hot.Inlet.Fw/(Unity(Nb+1).Properties.Hot.Outlet.rho*Aonozzle_Tube);   
    626635
    627636"Shell Pressure Drop Inlet Nozzle"
    628         Unity(1).Shell.PressureDrop.Pdnozzle_in = HE.DeltaPshellNozzleIn(Unity(1).Properties.Cold.Inlet.rho,Unity(1).Properties.Cold.Inlet.Fw);
     637        Unity(1).Shell.PressureDrop.Pdnozzle_in = (0.5*Unity(1).Properties.Cold.Inlet.Fw^2/Unity(1).Properties.Cold.Inlet.rho)*((1/Ainozzle_Shell^2)+(1/Aeinozzle_Shell^2));
    629638
    630639"Velocity Shell Side Inlet Nozzle"
    631         Unity(1).Shell.PressureDrop.Vnozzle_in  = HE.ShellVelocityNozzleIn(Unity(1).Properties.Cold.Inlet.rho,Unity(1).Properties.Cold.Inlet.Fw);
     640        Unity(1).Shell.PressureDrop.Vnozzle_in  = Unity(1).Properties.Cold.Inlet.Fw/(Unity(1).Properties.Cold.Inlet.rho*Ainozzle_Shell);
    632641
    633642"Shell Pressure Drop Outlet Nozzle"
    634         Unity(Nb+1).Shell.PressureDrop.Pdnozzle_out = HE.DeltaPshellNozzleOut(Unity(Nb+1).Properties.Cold.Outlet.rho,Unity(Nb+1).Properties.Cold.Inlet.Fw);
     643        Unity(Nb+1).Shell.PressureDrop.Pdnozzle_out = (0.5*Unity(Nb+1).Properties.Cold.Outlet.Fw^2/Unity(Nb+1).Properties.Cold.Outlet.rho)*((1/Ainozzle_Shell^2)+(1/Aeinozzle_Shell^2));
    635644
    636645"Velocity Shell Side Outlet Nozzle"
    637         Unity(Nb+1).Shell.PressureDrop.Vnozzle_out  = HE.ShellVelocityNozzleOut(Unity(Nb+1).Properties.Cold.Outlet.rho,Unity(Nb+1).Properties.Cold.Inlet.Fw);
    638 
     646        Unity(Nb+1).Shell.PressureDrop.Vnozzle_out  = Unity(Nb+1).Properties.Cold.Outlet.Fw/(Unity(Nb+1).Properties.Cold.Outlet.rho*Ainozzle_Shell);
    639647
    640648end
    … …  
    922930        HE.Ltd                          = Ltd;
    923931        side                            = HE.FluidAlocation();
     932       
     933#"Tube Side Inlet Nozzle Area"
     934        Ainozzle_Tube = (Pi*Dinozzle_Tube*Dinozzle_Tube)/4;
     935       
     936#"Tube Side Outlet Nozzle Area"
     937        Aonozzle_Tube = (Pi*Donozzle_Tube*Donozzle_Tube)/4;
     938       
     939#"Tube Inlet Nozzle Pressure Loss Coeff"
     940        Kinlet_Tube   = 1.1;
     941
     942#"Tube Outlet Nozzle Pressure Loss Coeff"
     943        Koutlet_Tube  = 0.7;
     944       
     945#"Shell Outlet Nozzle Area"
     946        Aonozzle_Shell  = (Pi*Donozzle_Shell*Donozzle_Shell)/4;
     947       
     948#"Shell Inlet Nozzle Area"
     949        Ainozzle_Shell  = (Pi*Dinozzle_Shell*Dinozzle_Shell)/4;
     950       
     951#"Shell Outlet Escape Area Under Nozzle"
     952        Aeonozzle_Shell = Pi*Donozzle_Shell*Honozzle_Shell + 0.6*Aonozzle_Shell*(1-Dotube/pitch);
     953       
     954#"Shell Inlet Escape Area Under Nozzle"
     955        Aeinozzle_Shell = Pi*Dinozzle_Shell*Hinozzle_Shell + 0.6*Ainozzle_Shell*(1-Dotube/pitch);
     956       
     957       
    924958end
    925959
    … …  
    943977Donozzle_Shell  as length               (Brief="Shell Outlet Nozzle Diameter",Lower=10e-6);
    944978Dinozzle_Shell  as length               (Brief="Shell Inlet Nozzle Diameter",Lower=10e-6);
     979Aonozzle_Shell  as area                 (Brief="Shell Outlet Nozzle Area",Lower=10e-6);
     980Ainozzle_Shell  as area                 (Brief="Shell Inlet Nozzle Area",Lower=10e-6);
     981Aeonozzle_Shell as area                 (Brief="Shell Outlet Escape Area Under Nozzle",Lower=10e-6);
     982Aeinozzle_Shell as area                 (Brief="Shell Inlet Escape Area Under Nozzle",Lower=10e-6);
    945983Hinozzle_Shell  as length               (Brief="Height Under Shell Inlet Nozzle",Lower=10e-6);
    946984Honozzle_Shell  as length               (Brief="Height Under Shell Outlet Nozzle",Lower=10e-6);
    947985Lcf                     as length               (Brief="Bundle-to-Shell Clearance",Lower=10e-8);
    948 #+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
    949 #       Tubes Geometrical Parameters
    950 #+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
    951 Ntt                             as Integer              (Brief="Total Number of Tubes in Shell",Default=100,Lower=1);
    952 Pattern         as Integer              (Brief="Tube Layout Characteristic Angle",Lower=30);
    953 Ltube                   as length               (Brief="Effective Tube Length",Lower=0.1);
    954 pitch                   as length               (Brief="Tube Pitch",Lower=1e-8);
    955 Kwall                   as conductivity (Brief="Tube Wall Material Thermal Conductivity");
    956 Dotube                  as length               (Brief="Tube Outside Diameter",Lower=10e-6);
    957 Ditube                  as length               (Brief="Tube Inside Diameter",Lower=10e-6);
    958 Donozzle_Tube   as length               (Brief="Tube Outlet Nozzle Diameter",Lower=10e-6);
    959 Dinozzle_Tube   as length               (Brief="Tube Inlett Nozzle Diameter",Lower=10e-6);
     986
     987#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
     988#                               Tubes Geometrical Parameters                                            #
     989#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
     990Ntt                       as Integer            (Brief="Total Number of Tubes in Shell",Default=100,Lower=1);
     991Pattern       as Integer                (Brief="Tube Layout Characteristic Angle",Lower=30);
     992Ltube             as length                     (Brief="Effective Tube Length",Lower=0.1);
     993pitch             as length                     (Brief="Tube Pitch",Lower=1e-8);
     994Kwall             as conductivity       (Brief="Tube Wall Material Thermal Conductivity");
     995Dotube            as length                     (Brief="Tube Outside Diameter",Lower=10e-6);
     996Ditube            as length                     (Brief="Tube Inside Diameter",Lower=10e-6);
     997Donozzle_Tube as length                 (Brief="Tube Outlet Nozzle Diameter",Lower=10e-6);
     998Dinozzle_Tube as length                 (Brief="Tube Inlet Nozzle Diameter",Lower=10e-6);
     999Aonozzle_Tube as area                   (Brief="Tube Outlet Nozzle Area",Lower=10e-6);
     1000Ainozzle_Tube as area                   (Brief="Tube Inlet Nozzle Area",Lower=10e-6);
     1001Kinlet_Tube   as positive               (Brief="Tube Inlet Nozzle Pressure Loss Coeff",Default=1.1);
     1002Koutlet_Tube  as positive               (Brief="Tube Outlet Nozzle Pressure Loss Coeff",Default=0.7);
     1003
    9601004#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
    9611005#       Baffles Geometrical Parameters
    … …  
    10441088       
    10451089        then
    1046 
     1090       
    10471091"Shell Side inlet Nozzle rho-V^2"
    10481092        Unity(1).Shell.PressureDrop.RVsquare_in = Unity(1).Properties.Hot.Inlet.rho*(Unity(1).Shell.PressureDrop.Vnozzle_in)^2;
    … …  
    10601104
    10611105"Pressure Drop Tube Side Inlet Nozzle"
    1062         Unity(1).Tubes.PressureDrop.Pdnozzle_in = HE.DeltaPtubeNozzlein(Unity(1).Properties.Cold.Inlet.rho,Unity(1).Properties.Cold.Inlet.Fw);
     1106        Unity(1).Tubes.PressureDrop.Pdnozzle_in = 0.5*Kinlet_Tube*Unity(1).Properties.Cold.Inlet.rho*Unity(1).Tubes.PressureDrop.Vnozzle_in^2;
    10631107
    10641108"Velocity Tube Side Inlet Nozzle"
    1065         Unity(1).Tubes.PressureDrop.Vnozzle_in  = HE.TubeVelocityNozzlein(Unity(1).Properties.Cold.Inlet.rho,Unity(1).Properties.Cold.Inlet.Fw);
     1109        Unity(1).Tubes.PressureDrop.Vnozzle_in  = Unity(1).Properties.Cold.Inlet.Fw/(Unity(1).Properties.Cold.Inlet.rho*Ainozzle_Tube);
    10661110
    10671111"Pressure Drop Tube Side Outlet Nozzle"
    1068         Unity(Nb+1).Tubes.PressureDrop.Pdnozzle_out     = HE.DeltaPtubeNozzleout(Unity(Nb+1).Properties.Cold.Outlet.rho,Unity(Nb+1).Properties.Cold.Inlet.Fw);
     1112        Unity(Nb+1).Tubes.PressureDrop.Pdnozzle_out     = 0.5*Koutlet_Tube*Unity(Nb+1).Properties.Cold.Outlet.rho*Unity(Nb+1).Tubes.PressureDrop.Vnozzle_out^2;
    10691113
    10701114"Velocity Tube Side Outlet Nozzle"
    1071         Unity(Nb+1).Tubes.PressureDrop.Vnozzle_out      = HE.TubeVelocityNozzleout(Unity(Nb+1).Properties.Cold.Outlet.rho,Unity(Nb+1).Properties.Cold.Inlet.Fw);
     1115        Unity(Nb+1).Tubes.PressureDrop.Vnozzle_out      = Unity(Nb+1).Properties.Cold.Inlet.Fw/(Unity(Nb+1).Properties.Cold.Outlet.rho*Aonozzle_Tube);
    10721116
    10731117"Shell Pressure Drop Inlet Nozzle"
    1074         Unity(1).Shell.PressureDrop.Pdnozzle_in = HE.DeltaPshellNozzleIn(Unity(1).Properties.Hot.Inlet.rho,Unity(1).Properties.Hot.Inlet.Fw);
     1118        Unity(1).Shell.PressureDrop.Pdnozzle_in = (0.5*Unity(1).Properties.Hot.Inlet.Fw^2/Unity(1).Properties.Hot.Inlet.rho)*((1/Ainozzle_Shell^2)+(1/Aeinozzle_Shell^2));
    10751119
    10761120"Velocity Shell Side Inlet Nozzle"
    1077         Unity(1).Shell.PressureDrop.Vnozzle_in  = HE.ShellVelocityNozzleIn(Unity(1).Properties.Hot.Inlet.rho,Unity(1).Properties.Hot.Inlet.Fw);
     1121        Unity(1).Shell.PressureDrop.Vnozzle_in  = Unity(1).Properties.Hot.Inlet.Fw/(Unity(1).Properties.Hot.Inlet.rho*Ainozzle_Shell);
    10781122
    10791123"Shell Pressure Drop Outlet Nozzle"
    1080         Unity(Nb+1).Shell.PressureDrop.Pdnozzle_out     =HE.DeltaPshellNozzleOut(Unity(Nb+1).Properties.Hot.Outlet.rho,Unity(Nb+1).Properties.Hot.Inlet.Fw);
     1124        Unity(Nb+1).Shell.PressureDrop.Pdnozzle_out     = (0.5*Unity(Nb+1).Properties.Hot.Outlet.Fw^2/Unity(Nb+1).Properties.Hot.Outlet.rho)*((1/Ainozzle_Shell^2)+(1/Aeinozzle_Shell^2));
    10811125
    10821126"Velocity Shell Side Outlet Nozzle"
    1083         Unity(Nb+1).Shell.PressureDrop.Vnozzle_out      =HE.ShellVelocityNozzleOut(Unity(Nb+1).Properties.Hot.Outlet.rho,Unity(Nb+1).Properties.Hot.Inlet.Fw);
    1084        
     1127        Unity(Nb+1).Shell.PressureDrop.Vnozzle_out      = Unity(Nb+1).Properties.Hot.Outlet.Fw/(Unity(Nb+1).Properties.Hot.Outlet.rho*Aonozzle_Shell);
     1128
    10851129
    10861130        else
    1087 
     1131       
    10881132"Shell Side inlet Nozzle rho-V^2"
    10891133        Unity(1).Shell.PressureDrop.RVsquare_in = Unity(1).Properties.Cold.Inlet.rho*(Unity(1).Shell.PressureDrop.Vnozzle_in)^2;
    … …  
    11001144        Unity(Nb+1).Baffles.Lsi,Unity(Nb+1).Properties.Cold.Inlet.Fw,Unity(Nb+1).Shell.HeatTransfer.Phi,Unity(Nb+1).Properties.Cold.Average.rho);
    11011145
    1102 
    11031146"Pressure Drop Tube Side Inlet Nozzle"
    1104         Unity(1).Tubes.PressureDrop.Pdnozzle_in = HE.DeltaPtubeNozzlein(Unity(1).Properties.Hot.Inlet.rho,Unity(1).Properties.Hot.Inlet.Fw);
    1105 
     1147        Unity(1).Tubes.PressureDrop.Pdnozzle_in = 0.5*Kinlet_Tube*Unity(1).Properties.Hot.Inlet.rho*Unity(1).Tubes.PressureDrop.Vnozzle_in^2;
     1148       
    11061149"Velocity Tube Side Inlet Nozzle"
    1107         Unity(1).Tubes.PressureDrop.Vnozzle_in  = HE.TubeVelocityNozzlein(Unity(1).Properties.Hot.Inlet.rho,Unity(1).Properties.Hot.Inlet.Fw);
    1108 
     1150        Unity(1).Tubes.PressureDrop.Vnozzle_in  = Unity(1).Properties.Hot.Inlet.Fw/(Unity(1).Properties.Hot.Inlet.rho*Ainozzle_Tube);
     1151       
    11091152"Pressure Drop Tube Side Outlet Nozzle"
    1110         Unity(Nb+1).Tubes.PressureDrop.Pdnozzle_out     = HE.DeltaPtubeNozzleout(Unity(Nb+1).Properties.Hot.Outlet.rho,Unity(Nb+1).Properties.Hot.Inlet.Fw);   
    1111 
     1153        Unity(Nb+1).Tubes.PressureDrop.Pdnozzle_out     = 0.5*Koutlet_Tube*Unity(Nb+1).Properties.Hot.Outlet.rho*Unity(Nb+1).Tubes.PressureDrop.Vnozzle_out^2;
     1154       
    11121155"Velocity Tube Side Outlet Nozzle"
    1113         Unity(Nb+1).Tubes.PressureDrop.Vnozzle_out      = HE.TubeVelocityNozzleout(Unity(Nb+1).Properties.Hot.Outlet.rho,Unity(Nb+1).Properties.Hot.Inlet.Fw); 
     1156        Unity(Nb+1).Tubes.PressureDrop.Vnozzle_out      = Unity(Nb+1).Properties.Hot.Inlet.Fw/(Unity(Nb+1).Properties.Hot.Outlet.rho*Aonozzle_Tube);   
    11141157
    11151158"Shell Pressure Drop Inlet Nozzle"
    1116         Unity(1).Shell.PressureDrop.Pdnozzle_in = HE.DeltaPshellNozzleIn(Unity(1).Properties.Cold.Inlet.rho,Unity(1).Properties.Cold.Inlet.Fw);
     1159        Unity(1).Shell.PressureDrop.Pdnozzle_in = (0.5*Unity(1).Properties.Cold.Inlet.Fw^2/Unity(1).Properties.Cold.Inlet.rho)*((1/Ainozzle_Shell^2)+(1/Aeinozzle_Shell^2));
    11171160
    11181161"Velocity Shell Side Inlet Nozzle"
    1119         Unity(1).Shell.PressureDrop.Vnozzle_in  = HE.ShellVelocityNozzleIn(Unity(1).Properties.Cold.Inlet.rho,Unity(1).Properties.Cold.Inlet.Fw);
     1162        Unity(1).Shell.PressureDrop.Vnozzle_in  = Unity(1).Properties.Cold.Inlet.Fw/(Unity(1).Properties.Cold.Inlet.rho*Ainozzle_Shell);
    11201163
    11211164"Shell Pressure Drop Outlet Nozzle"
    1122         Unity(Nb+1).Shell.PressureDrop.Pdnozzle_out = HE.DeltaPshellNozzleOut(Unity(Nb+1).Properties.Cold.Outlet.rho,Unity(Nb+1).Properties.Cold.Inlet.Fw);
     1165        Unity(Nb+1).Shell.PressureDrop.Pdnozzle_out = (0.5*Unity(Nb+1).Properties.Cold.Outlet.Fw^2/Unity(Nb+1).Properties.Cold.Outlet.rho)*((1/Ainozzle_Shell^2)+(1/Aeinozzle_Shell^2));
    11231166
    11241167"Velocity Shell Side Outlet Nozzle"
    1125         Unity(Nb+1).Shell.PressureDrop.Vnozzle_out  = HE.ShellVelocityNozzleOut(Unity(Nb+1).Properties.Cold.Outlet.rho,Unity(Nb+1).Properties.Cold.Inlet.Fw);
     1168        Unity(Nb+1).Shell.PressureDrop.Vnozzle_out  = Unity(Nb+1).Properties.Cold.Outlet.Fw/(Unity(Nb+1).Properties.Cold.Outlet.rho*Ainozzle_Shell);
    11261169
    11271170
    … …  
    14121455        HE.Ltd                          = Ltd;
    14131456        side                            = HE.FluidAlocation();
    1414 end
     1457       
     1458#"Tube Side Inlet Nozzle Area"
     1459        Ainozzle_Tube = (Pi*Dinozzle_Tube*Dinozzle_Tube)/4;
     1460       
     1461#"Tube Side Outlet Nozzle Area"
     1462        Aonozzle_Tube = (Pi*Donozzle_Tube*Donozzle_Tube)/4;
     1463       
     1464#"Tube Inlet Nozzle Pressure Loss Coeff"
     1465        Kinlet_Tube   = 1.1;
     1466
     1467#"Tube Outlet Nozzle Pressure Loss Coeff"
     1468        Koutlet_Tube  = 0.7;
     1469       
     1470#"Shell Outlet Nozzle Area"
     1471        Aonozzle_Shell  = (Pi*Donozzle_Shell*Donozzle_Shell)/4;
     1472       
     1473#"Shell Inlet Nozzle Area"
     1474        Ainozzle_Shell  = (Pi*Dinozzle_Shell*Dinozzle_Shell)/4;
     1475       
     1476#"Shell Outlet Escape Area Under Nozzle"
     1477        Aeonozzle_Shell = Pi*Donozzle_Shell*Honozzle_Shell + 0.6*Aonozzle_Shell*(1-Dotube/pitch);
     1478       
     1479#"Shell Inlet Escape Area Under Nozzle"
     1480        Aeinozzle_Shell = Pi*Dinozzle_Shell*Hinozzle_Shell + 0.6*Ainozzle_Shell*(1-Dotube/pitch);
     1481       
     1482end
Note: See TracChangeset for help on using the changeset viewer.