Ignore:
Timestamp:
Dec 20, 2007, 7:06:05 PM (15 years ago)
Author:
gerson bicca
Message:

added nozzles pressure drop for double pipe heat exchanger model

File:
1 edited

Legend:

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

    r437 r438  
    160160fi(N)                                           as fricfactor           (Brief="Incremental Friction Factor", Default=0.05, Lower=1e-10, Upper=2000);
    161161Re(N)                                           as positive                     (Brief="Incremental Reynolds Number",Default=100,Lower=1);
    162 
     162Pdnozzle_in             as press_delta          (Brief="Inlet Nozzle Pressure Drop",Default=0.01, Lower=0,DisplayUnit='kPa');
     163Pdnozzle_out            as press_delta          (Brief="Outlet Nozzle Pressure Drop",Default=0.01, Lower=0,DisplayUnit='kPa');
     164Vnozzle_in              as velocity                     (Brief="Inlet Nozzle Velocity",Default=1, Upper=1e5, Lower=0);
     165Vnozzle_out             as velocity                     (Brief="Outlet Nozzle Velocity",Default=1, Upper=1e5, Lower=0);
     166RVsquare_in             as positive                     (Brief = "Inlet Nozzle rho-V^2", Default=1, Upper=1e6, Unit = 'kg/s^2/m');
     167RVsquare_out            as positive                     (Brief = "Outlet Nozzle rho-V^2", Default=1, Upper=1e6, Unit = 'kg/s^2/m');
    163168end     
    164169
     
    220225        Rfo                     as positive                     (Brief="Outside Fouling Resistance",Unit='m^2*K/kW',Default=1e-6,Lower=0);
    221226
     227Donozzle_Inner          as length               (Brief="Inner Side Outlet Nozzle Diameter",Default = 0.036,Lower=10e-6);
     228Dinozzle_Inner          as length               (Brief="Inner Side Inlet Nozzle Diameter",Default = 0.036,Lower=10e-6);
     229
     230Donozzle_Outer          as length               (Brief="Outer Side Outlet Nozzle Diameter",Default = 0.036,Lower=10e-6);
     231Dinozzle_Outer  as length               (Brief="Outer Side Inlet Nozzle Diameter",Default = 0.036,Lower=10e-6);
     232
     233InnerKinlet             as positive                     (Brief="Inner Side Inlet Nozzle Pressure Loss Coeff",Default=1.1);
     234InnerKoutlet    as positive                     (Brief="Inner Side Outlet Nozzle Pressure Loss Coeff",Default=0.7);
     235
     236OuterKinlet             as positive                     (Brief="Outer Side Inlet Nozzle Pressure Loss Coeff",Default=1.1);
     237OuterKoutlet    as positive                     (Brief="Outer Side Outlet Nozzle Pressure Loss Coeff",Default=0.7);
     238
    222239VARIABLES
    223240
     
    840857
    841858"Total Pressure Drop Outer Stream"
    842         Outer.PressureDrop.Pdrop  = Outer.PressureDrop.Pd_fric(Npoints);
     859        Outer.PressureDrop.Pdrop  = Outer.PressureDrop.Pd_fric(Npoints)+Outer.PressureDrop.Pdnozzle_in+Outer.PressureDrop.Pdnozzle_out;
    843860
    844861"Outer Pipe Pressure Drop for friction"
     
    850867for i in [1:N]
    851868   
    852 "Outer Pipe Local Pressure"
     869"Outer Pipe Local Pressure"# FIXME: NOZZLE PRESSURE DROP MUST BE ADDED
    853870        Outer.PressureDrop.Plocal(i+1) =        Outer.PressureDrop.Plocal(1) - Outer.PressureDrop.Pd_fric(i+1);
    854871
     
    858875
    859876"Total Pressure Drop Outer Stream"
    860         Outer.PressureDrop.Pdrop  = Outer.PressureDrop.Pd_fric(1);
     877        Outer.PressureDrop.Pdrop  = Outer.PressureDrop.Pd_fric(1)+Outer.PressureDrop.Pdnozzle_in+Outer.PressureDrop.Pdnozzle_out;
     878
    861879
    862880for i in [1:N]
     
    872890for i in [1:N]
    873891   
    874 "Outer Pipe Local Pressure"
     892"Outer Pipe Local Pressure"# FIXME: NOZZLE PRESSURE DROP MUST BE ADDED
    875893        Outer.PressureDrop.Plocal(i) =  Outer.PressureDrop.Plocal(Npoints) - Outer.PressureDrop.Pd_fric(i+1);
    876894
     
    880898
    881899"Total Pressure Drop Inner Stream"
    882         Inner.PressureDrop.Pdrop  = Inner.PressureDrop.Pd_fric(Npoints);
     900        Inner.PressureDrop.Pdrop  = Inner.PressureDrop.Pd_fric(Npoints)+Inner.PressureDrop.Pdnozzle_in+Inner.PressureDrop.Pdnozzle_out;
    883901       
    884902"Inner Pipe Pressure Drop for friction"
     
    890908for i in [1:N]
    891909
    892 "Inner Pipe Local Pressure"
     910"Inner Pipe Local Pressure"# FIXME: NOZZLE PRESSURE DROP MUST BE ADDED
    893911        Inner.PressureDrop.Plocal(i+1) =        Inner.PressureDrop.Plocal(1) - Inner.PressureDrop.Pd_fric(i+1);
    894912
    895913end
    896914
    897 end
     915
     916"Velocity Inner Side Inlet Nozzle"
     917        Inner.PressureDrop.Vnozzle_in   = Inner.Properties.Inlet.Fw/(Inner.Properties.Inlet.rho*(0.25*Pi*Dinozzle_Inner^2));
     918
     919"Velocity Inner Side Outlet Nozzle"
     920        Inner.PressureDrop.Vnozzle_out  = Inner.Properties.Outlet.Fw/(Inner.Properties.Outlet.rho*(0.25*Pi*Donozzle_Inner^2));
     921
     922"Velocity Outer Side Inlet Nozzle"
     923        Outer.PressureDrop.Vnozzle_in   = Outer.Properties.Inlet.Fw/(Outer.Properties.Inlet.rho*(0.25*Pi*Dinozzle_Outer^2));
     924
     925"Velocity Outer Side Outlet Nozzle"
     926        Outer.PressureDrop.Vnozzle_out  = Outer.Properties.Outlet.Fw/(Outer.Properties.Outlet.rho*(0.25*Pi*Donozzle_Outer^2));
     927
     928"Pressure Drop Inner Side Inlet Nozzle"
     929        Inner.PressureDrop.Pdnozzle_in  = 0.5*InnerKinlet*Inner.Properties.Inlet.rho*Inner.PressureDrop.Vnozzle_in^2;
     930
     931"Pressure Drop Inner Side Outlet Nozzle"
     932        Inner.PressureDrop.Pdnozzle_out = 0.5*InnerKoutlet*Inner.Properties.Outlet.rho*Inner.PressureDrop.Vnozzle_out^2;
     933
     934"Pressure Drop Outer Side Inlet Nozzle"
     935        Outer.PressureDrop.Pdnozzle_in  = 0.5*OuterKinlet*Outer.Properties.Inlet.rho*Outer.PressureDrop.Vnozzle_in^2;
     936
     937"Pressure Drop Outer Side Outlet Nozzle"
     938        Outer.PressureDrop.Pdnozzle_out = 0.5*OuterKoutlet*Outer.Properties.Outlet.rho*Outer.PressureDrop.Vnozzle_out^2;
     939
     940"Inner Side Inlet Nozzle rho-V^2"
     941        Inner.PressureDrop.RVsquare_in = Inner.Properties.Inlet.rho*(Inner.PressureDrop.Vnozzle_in)^2;
     942
     943"Inner Side Outlet Nozzle rho-V^2"
     944        Inner.PressureDrop.RVsquare_out = Inner.Properties.Outlet.rho*(Inner.PressureDrop.Vnozzle_out)^2;
     945
     946"Outer Side Inlet Nozzle rho-V^2"
     947        Outer.PressureDrop.RVsquare_in = Outer.Properties.Inlet.rho*(Outer.PressureDrop.Vnozzle_in)^2;
     948
     949"Outer Side Outlet Nozzle rho-V^2"
     950        Outer.PressureDrop.RVsquare_out = Outer.Properties.Outlet.rho*(Outer.PressureDrop.Vnozzle_out)^2;
     951
     952end
Note: See TracChangeset for help on using the changeset viewer.