Changeset 101 for mso/eml

Timestamp:

Jan 9, 2007, 1:59:59 PM (17 years ago)

Author:

Rafael de Pelegrini Soares

Message:

Updated flashPH model

File:

• mso/eml/stage_separators/flashPH.mso (modified) (3 diffs)

mso/eml/stage_separators/flashPH.mso

@@ -48 +48 @@
 in      Q as heat_rate (Brief="Rate of heat supply"); 
         vfrac as fraction(Brief="Real vaporization fraction");
-        vsat as Real(Lower=-1, Upper=2, Brief="Vaporization fraction if saturated");
+        vsat as Real(Lower=-5, Upper=5, Brief="Vaporization fraction if saturated");
         Tsat as temperature(Lower=173, Upper=473, Brief="Temperature if saturated");
         xsat(NComp) as Real(Lower=-5, Upper=5, Brief="Liquid composition if saturated");
@@ -58 +58 @@
         EQUATIONS
         "Chemical equilibrium"
-        PP.LiquidFugacityCoefficient(Tsat, OutletL.P, OutletL.z)*xsat = 
-                PP.VapourFugacityCoefficient(Tsat, OutletV.P, OutletV.z)*ysat;
+        PP.LiquidFugacityCoefficient(Tsat, OutletL.P, xsat)*xsat = 
+                PP.VapourFugacityCoefficient(Tsat, OutletV.P, ysat)*ysat;
 
         "Global Molar Balance"
@@ -71 +71 @@
         "Energy Balance if saturated"
         Inlet.F*Inlet.h  + Q =
-                Inlet.F*(1-vsat)*PP.LiquidEnthalpy(Tsat, OutletL.P, OutletL.z) +
-                Inlet.F*vsat*PP.VapourEnthalpy(Tsat, OutletV.P, OutletV.z);
+                Inlet.F*(1-vsat)*PP.LiquidEnthalpy(Tsat, OutletL.P, xsat) +
+                Inlet.F*vsat*PP.VapourEnthalpy(Tsat, OutletV.P, ysat);
 
         "Real Energy Balance"