Changeset 231

Timestamp:

Mar 26, 2007, 11:12:04 AM (17 years ago)

Author:

gerson bicca

Message:

updated simplified heat exchangers models to incorporate phase change

Location:

branches/newlanguage

Files:

• eml/heat_exchangers/HeatExchangerSimplified.mso (modified) (6 diffs)
• sample/heat_exchangers/GuessSimplified.rlt (added)
• sample/heat_exchangers/Sample_Simplified.mso (modified) (1 diff)
• sample/heat_exchangers/sample_ShellandTubes.mso (modified) (1 diff)

branches/newlanguage/eml/heat_exchangers/HeatExchangerSimplified.mso

@@ -21 +21 @@
 
 Model HeatExchangerSimplified_Basic 
-        
+
 ATTRIBUTES
         Pallete         = false;
@@ -41 +41 @@
 out     OutletCold  as streamPH         (Brief="Outlet Cold Stream");
 
+        xh(NComp)       as fraction(Brief = "Liquid Molar Fraction in Hot Side");
+        yh(NComp) as fraction(Brief = "Vapour Molar Fraction in Hot Side");
+        vh                      as fraction(Brief = "Vapour Molar Fraction in Hot Side");
+        
+        xc(NComp)       as fraction(Brief = "Liquid Molar Fraction in Cold Side");
+        yc(NComp) as fraction(Brief = "Vapour Molar Fraction in Cold Side");
+        vc                      as fraction(Brief = "Vapour Molar Fraction in Cold Side");
+
         Details         as Details_Main                         (Brief="Heat Exchanger Details");
         HotSide         as Main_Simplified              (Brief="Heat Exchanger Hot Side");
@@ -52 +60 @@
 EQUATIONS
 
+"Flash Calculation in Hot Side"
+        [vh, xh, yh] = PP.Flash(InletHot.T, InletHot.P, InletHot.z);
+
+"Flash Calculation in Cold Side"
+        [vc, xc, yc] = PP.Flash(InletCold.T, InletCold.P, InletCold.z);
+
 "Hot Stream Average Temperature"
         HotSide.Properties.Average.T = 0.5*InletHot.T + 0.5*OutletHot.T;
@@ -76 +90 @@
         ColdSide.Properties.Average.Mw = sum(M*InletCold.z);
 
-if InletCold.v equal 0
-        
-        then    
-        
 "Cold Stream Average Heat Capacity"
-        ColdSide.Properties.Average.Cp  =       PP.LiquidCp(ColdSide.Properties.Average.T,ColdSide.Properties.Average.P,InletCold.z);
+        ColdSide.Properties.Average.Cp  =       (1-InletCold.v)*PP.LiquidCp(ColdSide.Properties.Average.T,ColdSide.Properties.Average.P,xc) +
+                InletCold.v*PP.VapourCp(ColdSide.Properties.Average.T,ColdSide.Properties.Average.P,yc);
 
 "Cold Stream Inlet Heat Capacity"
-        ColdSide.Properties.Inlet.Cp    =       PP.LiquidCp(InletCold.T,InletCold.P,InletCold.z);
+        ColdSide.Properties.Inlet.Cp    =       (1-InletCold.v)*PP.LiquidCp(InletCold.T,InletCold.P,xc)+
+                InletCold.v*PP.VapourCp(InletCold.T,InletCold.P,yc);
 
 "Cold Stream Outlet Heat Capacity"
-        ColdSide.Properties.Outlet.Cp   =       PP.LiquidCp(OutletCold.T,OutletCold.P,OutletCold.z);
+        ColdSide.Properties.Outlet.Cp   =       (1-OutletCold.v)*PP.LiquidCp(OutletCold.T,OutletCold.P,OutletCold.x)+
+                OutletCold.v*PP.VapourCp(OutletCold.T,OutletCold.P,OutletCold.y);
 
 "Cold Stream Average Mass Density"
-        ColdSide.Properties.Average.rho =       PP.LiquidDensity(ColdSide.Properties.Average.T,ColdSide.Properties.Average.P,InletCold.z);
+        ColdSide.Properties.Average.rho =       (1-InletCold.v)*PP.LiquidDensity(ColdSide.Properties.Average.T,ColdSide.Properties.Average.P,xc)+
+                InletCold.v*PP.VapourDensity(ColdSide.Properties.Average.T,ColdSide.Properties.Average.P,yc);
 
 "Cold Stream Inlet Mass Density"
-        ColdSide.Properties.Inlet.rho   =       PP.LiquidDensity(InletCold.T,InletCold.P,InletCold.z);
+        ColdSide.Properties.Inlet.rho   =       (1-InletCold.v)*PP.LiquidDensity(InletCold.T,InletCold.P,xc)+
+                InletCold.v*PP.VapourDensity(InletCold.T,InletCold.P,yc);
 
 "Cold Stream Outlet Mass Density"
-        ColdSide.Properties.Outlet.rho  =       PP.LiquidDensity(OutletCold.T,OutletCold.P,OutletCold.z);
+        ColdSide.Properties.Outlet.rho  =       (1-OutletCold.v)*PP.LiquidDensity(OutletCold.T,OutletCold.P,OutletCold.x)+
+                OutletCold.v*PP.VapourDensity(OutletCold.T,OutletCold.P,OutletCold.y);
 
 "Cold Stream Average Viscosity"
-        ColdSide.Properties.Average.Mu  =       PP.LiquidViscosity(ColdSide.Properties.Average.T,ColdSide.Properties.Average.P,InletCold.z);
+        ColdSide.Properties.Average.Mu  =       (1-InletCold.v)*PP.LiquidViscosity(ColdSide.Properties.Average.T,ColdSide.Properties.Average.P,xc)+
+                InletCold.v*PP.VapourViscosity(ColdSide.Properties.Average.T,ColdSide.Properties.Average.P,yc);
 
 "Cold Stream inlet Viscosity"
-        ColdSide.Properties.Inlet.Mu    =       PP.LiquidViscosity(InletCold.T,InletCold.P,InletCold.z);
-        
+        ColdSide.Properties.Inlet.Mu    =       (1-InletCold.v)*PP.LiquidViscosity(InletCold.T,InletCold.P,xc)+
+                InletCold.v*PP.VapourViscosity(InletCold.T,InletCold.P,yc);
+
 "Cold Stream Outlet Viscosity"
-        ColdSide.Properties.Outlet.Mu   =       PP.LiquidViscosity(OutletCold.T,OutletCold.P,OutletCold.z);
+        ColdSide.Properties.Outlet.Mu   =       (1-OutletCold.v)*PP.LiquidViscosity(OutletCold.T,OutletCold.P,OutletCold.x)+
+                OutletCold.v*PP.VapourViscosity(OutletCold.T,OutletCold.P,OutletCold.y);
 
 "Cold Stream Average Conductivity"
-        ColdSide.Properties.Average.K   =       PP.LiquidThermalConductivity(ColdSide.Properties.Average.T,ColdSide.Properties.Average.P,InletCold.z);
+        ColdSide.Properties.Average.K   =       (1-InletCold.v)*PP.LiquidThermalConductivity(ColdSide.Properties.Average.T,ColdSide.Properties.Average.P,xc)+
+                InletCold.v*PP.VapourThermalConductivity(ColdSide.Properties.Average.T,ColdSide.Properties.Average.P,yc);
 
 "Cold Stream Inlet Conductivity"        
-        ColdSide.Properties.Inlet.K     =       PP.LiquidThermalConductivity(InletCold.T,InletCold.P,InletCold.z);
+        ColdSide.Properties.Inlet.K     =       (1-InletCold.v)*PP.LiquidThermalConductivity(InletCold.T,InletCold.P,xc)+
+                InletCold.v*PP.VapourThermalConductivity(InletCold.T,InletCold.P,yc);
 
 "Cold Stream Outlet Conductivity"
-        ColdSide.Properties.Outlet.K    =       PP.LiquidThermalConductivity(OutletCold.T,OutletCold.P,OutletCold.z);
-
+        ColdSide.Properties.Outlet.K    =       (1-OutletCold.v)*PP.LiquidThermalConductivity(OutletCold.T,OutletCold.P,OutletCold.x)+
+                OutletCold.v*PP.VapourThermalConductivity(OutletCold.T,OutletCold.P,OutletCold.y);
+        
 "Cold Stream Viscosity at Wall Temperature"
-        ColdSide.Properties.Wall.Mu     =       PP.LiquidViscosity(ColdSide.Properties.Wall.Twall,ColdSide.Properties.Average.P,InletCold.z);
-
-
-        else
-
-"Cold Stream Average Heat Capacity"
-        ColdSide.Properties.Average.Cp  =       PP.VapourCp(ColdSide.Properties.Average.T,ColdSide.Properties.Average.P,InletCold.z);
-
-"Cold Stream Inlet Heat Capacity"       
-        ColdSide.Properties.Inlet.Cp    =       PP.VapourCp(InletCold.T,InletCold.P,InletCold.z);
-
-"Cold Stream Outlet Heat Capacity"      
-        ColdSide.Properties.Outlet.Cp   =       PP.VapourCp(OutletCold.T,OutletCold.P,OutletCold.z);
-
-"Cold Stream Average Mass Density"
-        ColdSide.Properties.Average.rho =       PP.VapourDensity(ColdSide.Properties.Average.T,ColdSide.Properties.Average.P,InletCold.z);
-
-"Cold Stream Inlet Mass Density"
-        ColdSide.Properties.Inlet.rho   =       PP.VapourDensity(InletCold.T,InletCold.P,InletCold.z);
-
-"Cold Stream Outlet Mass Density"       
-        ColdSide.Properties.Outlet.rho  =       PP.VapourDensity(OutletCold.T,OutletCold.P,OutletCold.z);
-
-"Cold Stream Average Viscosity "
-        ColdSide.Properties.Average.Mu  =       PP.VapourViscosity(ColdSide.Properties.Average.T,ColdSide.Properties.Average.P,InletCold.z);
-
-"Cold Stream Inlet Viscosity "  
-        ColdSide.Properties.Inlet.Mu    =       PP.VapourViscosity(InletCold.T,InletCold.P,InletCold.z);
-
-"Cold Stream Outlet Viscosity " 
-        ColdSide.Properties.Outlet.Mu   =       PP.VapourViscosity(OutletCold.T,OutletCold.P,OutletCold.z);
-
-"Cold Stream Average Conductivity "
-        ColdSide.Properties.Average.K   =       PP.VapourThermalConductivity(ColdSide.Properties.Average.T,ColdSide.Properties.Average.P,InletCold.z);
-
-"Cold Stream Inlet Conductivity "
-        ColdSide.Properties.Inlet.K     =       PP.VapourThermalConductivity(InletCold.T,InletCold.P,InletCold.z);
-
-"Cold Stream Outlet Conductivity "
-        ColdSide.Properties.Outlet.K    =       PP.VapourThermalConductivity(OutletCold.T,OutletCold.P,OutletCold.z);
-        
-"Cold Stream Viscosity at Wall Temperature"
-        ColdSide.Properties.Wall.Mu     =       PP.VapourViscosity(ColdSide.Properties.Wall.Twall,ColdSide.Properties.Average.P,InletCold.z);
-
-end
-
-if InletHot.v equal 0
-
-        then
+        ColdSide.Properties.Wall.Mu     =       (1-InletCold.v)*PP.LiquidViscosity(ColdSide.Properties.Wall.Twall,ColdSide.Properties.Average.P,xc)+
+                InletCold.v*PP.VapourViscosity(ColdSide.Properties.Wall.Twall,ColdSide.Properties.Average.P,yc);
 
 "Hot Stream Average Heat Capacity"
-        HotSide.Properties.Average.Cp   =               PP.LiquidCp(HotSide.Properties.Average.T,HotSide.Properties.Average.P,InletHot.z);
+        HotSide.Properties.Average.Cp   =       (1-InletHot.v)*PP.LiquidCp(HotSide.Properties.Average.T,HotSide.Properties.Average.P,xc) +
+                InletHot.v*PP.VapourCp(HotSide.Properties.Average.T,HotSide.Properties.Average.P,yc);
 
 "Hot Stream Inlet Heat Capacity"
-        HotSide.Properties.Inlet.Cp     =               PP.LiquidCp(InletHot.T,InletHot.P,InletHot.z);
+        HotSide.Properties.Inlet.Cp     =       (1-InletHot.v)*PP.LiquidCp(InletHot.T,InletHot.P,xc)+
+                InletHot.v*PP.VapourCp(InletHot.T,InletHot.P,yc);
 
 "Hot Stream Outlet Heat Capacity"
-        HotSide.Properties.Outlet.Cp    =               PP.LiquidCp(OutletHot.T,OutletHot.P,OutletHot.z);
+        HotSide.Properties.Outlet.Cp    =       (1-OutletHot.v)*PP.LiquidCp(OutletHot.T,OutletHot.P,OutletHot.x)+
+                OutletHot.v*PP.VapourCp(OutletHot.T,OutletHot.P,OutletHot.y);
 
 "Hot Stream Average Mass Density"
-        HotSide.Properties.Average.rho  =               PP.LiquidDensity(HotSide.Properties.Average.T,HotSide.Properties.Average.P,InletHot.z);
-
-"Hot Stream Inlet Mass Density" 
-        HotSide.Properties.Inlet.rho    =               PP.LiquidDensity(InletHot.T,InletHot.P,InletHot.z);
-
-"Hot Stream Outlet Mass Density"        
-        HotSide.Properties.Outlet.rho   =               PP.LiquidDensity(OutletHot.T,OutletHot.P,OutletHot.z);
+        HotSide.Properties.Average.rho =        (1-InletHot.v)*PP.LiquidDensity(HotSide.Properties.Average.T,HotSide.Properties.Average.P,xc)+
+                InletHot.v*PP.VapourDensity(HotSide.Properties.Average.T,HotSide.Properties.Average.P,yc);
+
+"Hot Stream Inlet Mass Density"
+        HotSide.Properties.Inlet.rho    =       (1-InletHot.v)*PP.LiquidDensity(InletHot.T,InletHot.P,xc)+
+                InletHot.v*PP.VapourDensity(InletHot.T,InletHot.P,yc);
+
+"Hot Stream Outlet Mass Density"
+        HotSide.Properties.Outlet.rho   =       (1-OutletHot.v)*PP.LiquidDensity(OutletHot.T,OutletHot.P,OutletHot.x)+
+                OutletHot.v*PP.VapourDensity(OutletHot.T,OutletHot.P,OutletHot.y);
 
 "Hot Stream Average Viscosity"
-        HotSide.Properties.Average.Mu   =               PP.LiquidViscosity(HotSide.Properties.Average.T,HotSide.Properties.Average.P,InletHot.z);       
-
-"Hot Stream Inlet Viscosity"
-        HotSide.Properties.Inlet.Mu     =               PP.LiquidViscosity(InletHot.T,InletHot.P,InletHot.z);   
+        HotSide.Properties.Average.Mu   =       (1-InletHot.v)*PP.LiquidViscosity(HotSide.Properties.Average.T,HotSide.Properties.Average.P,xc)+
+                InletHot.v*PP.VapourViscosity(HotSide.Properties.Average.T,HotSide.Properties.Average.P,yc);
+
+"Hot Stream inlet Viscosity"
+        HotSide.Properties.Inlet.Mu     =       (1-InletHot.v)*PP.LiquidViscosity(InletHot.T,InletHot.P,xc)+
+                InletHot.v*PP.VapourViscosity(InletHot.T,InletHot.P,yc);
 
 "Hot Stream Outlet Viscosity"
-        HotSide.Properties.Outlet.Mu    =               PP.LiquidViscosity(OutletHot.T,OutletHot.P,OutletHot.z);        
+        HotSide.Properties.Outlet.Mu    =       (1-OutletHot.v)*PP.LiquidViscosity(OutletHot.T,OutletHot.P,OutletHot.x)+
+                OutletHot.v*PP.VapourViscosity(OutletHot.T,OutletHot.P,OutletHot.y);
 
 "Hot Stream Average Conductivity"
-        HotSide.Properties.Average.K    =               PP.LiquidThermalConductivity(HotSide.Properties.Average.T,HotSide.Properties.Average.P,InletHot.z);     
-
-"Hot Stream Inlet Conductivity"
-        HotSide.Properties.Inlet.K              =               PP.LiquidThermalConductivity(InletHot.T,InletHot.P,InletHot.z); 
+        HotSide.Properties.Average.K    =       (1-InletHot.v)*PP.LiquidThermalConductivity(HotSide.Properties.Average.T,HotSide.Properties.Average.P,xc)+
+                InletHot.v*PP.VapourThermalConductivity(HotSide.Properties.Average.T,HotSide.Properties.Average.P,yc);
+
+"Hot Stream Inlet Conductivity" 
+        HotSide.Properties.Inlet.K      =       (1-InletHot.v)*PP.LiquidThermalConductivity(InletHot.T,InletHot.P,xc)+
+                InletHot.v*PP.VapourThermalConductivity(InletHot.T,InletHot.P,yc);
 
 "Hot Stream Outlet Conductivity"
-        HotSide.Properties.Outlet.K     =               PP.LiquidThermalConductivity(OutletHot.T,OutletHot.P,OutletHot.z);      
-
-"Hot Stream Viscosity  at Wall Temperature"
-        HotSide.Properties.Wall.Mu              =               PP.LiquidViscosity(HotSide.Properties.Wall.Twall,HotSide.Properties.Average.P,InletHot.z);      
-
-        else
-
-"Hot Stream Average Heat Capacity"
-        HotSide.Properties.Average.Cp   =               PP.VapourCp(HotSide.Properties.Average.T,HotSide.Properties.Average.P,InletHot.z);
-
-"Hot Stream Inlet Heat Capacity"
-        HotSide.Properties.Inlet.Cp     =               PP.VapourCp(InletHot.T,InletHot.P,InletHot.z);
-
-"Hot Stream Outlet Heat Capacity"
-        HotSide.Properties.Outlet.Cp    =               PP.VapourCp(OutletHot.T,OutletHot.P,OutletHot.z);
-
-"Hot Stream Average Mass Density"
-        HotSide.Properties.Average.rho  =               PP.VapourDensity(HotSide.Properties.Average.T,HotSide.Properties.Average.P,InletHot.z);
-
-"Hot Stream Inlet Mass Density" 
-        HotSide.Properties.Inlet.rho    =               PP.VapourDensity(InletHot.T,InletHot.P,InletHot.z);
-
-"Hot Stream Outlet Mass Density"
-        HotSide.Properties.Outlet.rho   =               PP.VapourDensity(OutletHot.T,OutletHot.P,OutletHot.z);
-
-"Hot Stream Average Viscosity"
-        HotSide.Properties.Average.Mu   =               PP.VapourViscosity(HotSide.Properties.Average.T,HotSide.Properties.Average.P,InletHot.z);
-
-"Hot Stream Inlet Viscosity"
-        HotSide.Properties.Inlet.Mu     =               PP.VapourViscosity(InletHot.T,InletHot.P,InletHot.z);
-
-"Hot Stream Outlet Viscosity"
-        HotSide.Properties.Outlet.Mu    =               PP.VapourViscosity(OutletHot.T,OutletHot.P,OutletHot.z);
-
-"Hot Stream Average Conductivity"
-        HotSide.Properties.Average.K    =               PP.VapourThermalConductivity(HotSide.Properties.Average.T,HotSide.Properties.Average.P,InletHot.z);     
-
-"Hot Stream Inlet Conductivity"
-        HotSide.Properties.Inlet.K              =               PP.VapourThermalConductivity(InletHot.T,InletHot.P,InletHot.z); 
-        
-"Hot Stream Outlet Conductivity"
-        HotSide.Properties.Outlet.K     =               PP.VapourThermalConductivity(OutletHot.T,OutletHot.P,OutletHot.z);      
-
+        HotSide.Properties.Outlet.K     =       (1-OutletHot.v)*PP.LiquidThermalConductivity(OutletHot.T,OutletHot.P,OutletHot.x)+
+                OutletHot.v*PP.VapourThermalConductivity(OutletHot.T,OutletHot.P,OutletHot.y);
+        
 "Hot Stream Viscosity at Wall Temperature"
-        HotSide.Properties.Wall.Mu              =               PP.VapourViscosity(HotSide.Properties.Wall.Twall,HotSide.Properties.Average.P,InletHot.z);
-
-end
+        HotSide.Properties.Wall.Mu      =       (1-InletHot.v)*PP.LiquidViscosity(HotSide.Properties.Wall.Twall,HotSide.Properties.Average.P,xc)+
+                InletHot.v*PP.VapourViscosity(HotSide.Properties.Wall.Twall,HotSide.Properties.Average.P,yc);
 
 "Energy Balance Hot Stream"
@@ -307 +252 @@
 VARIABLES
 
-Method as LMTD_Basic    (Brief="LMTD Method of Calculation");
+        Method as LMTD_Basic    (Brief="LMTD Method of Calculation");
 
 EQUATIONS
@@ -324 +269 @@
         Method.DTL = OutletHot.T - OutletCold.T;
 
-case "counter":
+        case "counter":
         
 "Temperature Difference at Inlet"

branches/newlanguage/sample/heat_exchangers/Sample_Simplified.mso

@@ -72 +72 @@
 OPTIONS
 
- Dynamic                          = false;
+ Dynamic = false;
+GuessFile = "GuessSimplified";
 
 end

branches/newlanguage/sample/heat_exchangers/sample_ShellandTubes.mso

@@ -136 +136 @@
 
  Dynamic        = false;
- GuessFile   = "GuessShellandTubes";
+ #GuessFile   = "GuessShellandTubes";
 
 end