Changeset 438 for trunk/eml

Timestamp:

Dec 20, 2007, 7:06:05 PM (16 years ago)

Author:

gerson bicca

Message:

added nozzles pressure drop for double pipe heat exchanger model

File:

• trunk/eml/heat_exchangers/DoublePipeIncr.mso (modified) (8 diffs)

trunk/eml/heat_exchangers/DoublePipeIncr.mso

@@ -160 +160 @@
 fi(N)                                           as fricfactor           (Brief="Incremental Friction Factor", Default=0.05, Lower=1e-10, Upper=2000);
 Re(N)                                           as positive                     (Brief="Incremental Reynolds Number",Default=100,Lower=1);
-
+Pdnozzle_in             as press_delta          (Brief="Inlet Nozzle Pressure Drop",Default=0.01, Lower=0,DisplayUnit='kPa');
+Pdnozzle_out            as press_delta          (Brief="Outlet Nozzle Pressure Drop",Default=0.01, Lower=0,DisplayUnit='kPa');
+Vnozzle_in              as velocity                     (Brief="Inlet Nozzle Velocity",Default=1, Upper=1e5, Lower=0);
+Vnozzle_out             as velocity                     (Brief="Outlet Nozzle Velocity",Default=1, Upper=1e5, Lower=0);
+RVsquare_in             as positive                     (Brief = "Inlet Nozzle rho-V^2", Default=1, Upper=1e6, Unit = 'kg/s^2/m');
+RVsquare_out            as positive                     (Brief = "Outlet Nozzle rho-V^2", Default=1, Upper=1e6, Unit = 'kg/s^2/m');
 end     
 
@@ -220 +225 @@
         Rfo                     as positive                     (Brief="Outside Fouling Resistance",Unit='m^2*K/kW',Default=1e-6,Lower=0);
 
+Donozzle_Inner          as length               (Brief="Inner Side Outlet Nozzle Diameter",Default = 0.036,Lower=10e-6);
+Dinozzle_Inner          as length               (Brief="Inner Side Inlet Nozzle Diameter",Default = 0.036,Lower=10e-6);
+
+Donozzle_Outer          as length               (Brief="Outer Side Outlet Nozzle Diameter",Default = 0.036,Lower=10e-6);
+Dinozzle_Outer  as length               (Brief="Outer Side Inlet Nozzle Diameter",Default = 0.036,Lower=10e-6);
+
+InnerKinlet             as positive                     (Brief="Inner Side Inlet Nozzle Pressure Loss Coeff",Default=1.1);
+InnerKoutlet    as positive                     (Brief="Inner Side Outlet Nozzle Pressure Loss Coeff",Default=0.7);
+
+OuterKinlet             as positive                     (Brief="Outer Side Inlet Nozzle Pressure Loss Coeff",Default=1.1);
+OuterKoutlet    as positive                     (Brief="Outer Side Outlet Nozzle Pressure Loss Coeff",Default=0.7);
+
 VARIABLES
 
@@ -840 +857 @@
 
 "Total Pressure Drop Outer Stream"
-        Outer.PressureDrop.Pdrop  = Outer.PressureDrop.Pd_fric(Npoints);
+        Outer.PressureDrop.Pdrop  = Outer.PressureDrop.Pd_fric(Npoints)+Outer.PressureDrop.Pdnozzle_in+Outer.PressureDrop.Pdnozzle_out;
 
 "Outer Pipe Pressure Drop for friction"
@@ -850 +867 @@
 for i in [1:N]
    
-"Outer Pipe Local Pressure"
+"Outer Pipe Local Pressure"# FIXME: NOZZLE PRESSURE DROP MUST BE ADDED
         Outer.PressureDrop.Plocal(i+1) =        Outer.PressureDrop.Plocal(1) - Outer.PressureDrop.Pd_fric(i+1);
 
@@ -858 +875 @@
 
 "Total Pressure Drop Outer Stream"
-        Outer.PressureDrop.Pdrop  = Outer.PressureDrop.Pd_fric(1);
+        Outer.PressureDrop.Pdrop  = Outer.PressureDrop.Pd_fric(1)+Outer.PressureDrop.Pdnozzle_in+Outer.PressureDrop.Pdnozzle_out;
+
 
 for i in [1:N]
@@ -872 +890 @@
 for i in [1:N]
    
-"Outer Pipe Local Pressure"
+"Outer Pipe Local Pressure"# FIXME: NOZZLE PRESSURE DROP MUST BE ADDED
         Outer.PressureDrop.Plocal(i) =  Outer.PressureDrop.Plocal(Npoints) - Outer.PressureDrop.Pd_fric(i+1);
 
@@ -880 +898 @@
 
 "Total Pressure Drop Inner Stream"
-        Inner.PressureDrop.Pdrop  = Inner.PressureDrop.Pd_fric(Npoints);
+        Inner.PressureDrop.Pdrop  = Inner.PressureDrop.Pd_fric(Npoints)+Inner.PressureDrop.Pdnozzle_in+Inner.PressureDrop.Pdnozzle_out;
         
 "Inner Pipe Pressure Drop for friction"
@@ -890 +908 @@
 for i in [1:N]
 
-"Inner Pipe Local Pressure"
+"Inner Pipe Local Pressure"# FIXME: NOZZLE PRESSURE DROP MUST BE ADDED
         Inner.PressureDrop.Plocal(i+1) =        Inner.PressureDrop.Plocal(1) - Inner.PressureDrop.Pd_fric(i+1);
 
 end
 
-end
+
+"Velocity Inner Side Inlet Nozzle"
+        Inner.PressureDrop.Vnozzle_in   = Inner.Properties.Inlet.Fw/(Inner.Properties.Inlet.rho*(0.25*Pi*Dinozzle_Inner^2));
+
+"Velocity Inner Side Outlet Nozzle"
+        Inner.PressureDrop.Vnozzle_out  = Inner.Properties.Outlet.Fw/(Inner.Properties.Outlet.rho*(0.25*Pi*Donozzle_Inner^2));
+
+"Velocity Outer Side Inlet Nozzle"
+        Outer.PressureDrop.Vnozzle_in   = Outer.Properties.Inlet.Fw/(Outer.Properties.Inlet.rho*(0.25*Pi*Dinozzle_Outer^2));
+
+"Velocity Outer Side Outlet Nozzle"
+        Outer.PressureDrop.Vnozzle_out  = Outer.Properties.Outlet.Fw/(Outer.Properties.Outlet.rho*(0.25*Pi*Donozzle_Outer^2));
+
+"Pressure Drop Inner Side Inlet Nozzle"
+        Inner.PressureDrop.Pdnozzle_in  = 0.5*InnerKinlet*Inner.Properties.Inlet.rho*Inner.PressureDrop.Vnozzle_in^2;
+
+"Pressure Drop Inner Side Outlet Nozzle"
+        Inner.PressureDrop.Pdnozzle_out = 0.5*InnerKoutlet*Inner.Properties.Outlet.rho*Inner.PressureDrop.Vnozzle_out^2;
+
+"Pressure Drop Outer Side Inlet Nozzle"
+        Outer.PressureDrop.Pdnozzle_in  = 0.5*OuterKinlet*Outer.Properties.Inlet.rho*Outer.PressureDrop.Vnozzle_in^2;
+
+"Pressure Drop Outer Side Outlet Nozzle"
+        Outer.PressureDrop.Pdnozzle_out = 0.5*OuterKoutlet*Outer.Properties.Outlet.rho*Outer.PressureDrop.Vnozzle_out^2;
+
+"Inner Side Inlet Nozzle rho-V^2"
+        Inner.PressureDrop.RVsquare_in = Inner.Properties.Inlet.rho*(Inner.PressureDrop.Vnozzle_in)^2;
+
+"Inner Side Outlet Nozzle rho-V^2"
+        Inner.PressureDrop.RVsquare_out = Inner.Properties.Outlet.rho*(Inner.PressureDrop.Vnozzle_out)^2;
+
+"Outer Side Inlet Nozzle rho-V^2"
+        Outer.PressureDrop.RVsquare_in = Outer.Properties.Inlet.rho*(Outer.PressureDrop.Vnozzle_in)^2;
+
+"Outer Side Outlet Nozzle rho-V^2"
+        Outer.PressureDrop.RVsquare_out = Outer.Properties.Outlet.rho*(Outer.PressureDrop.Vnozzle_out)^2;
+
+end