Changeset 823 for branches/gui/eml/stage_separators/flash.mso

Timestamp:

Aug 11, 2009, 11:12:13 AM (13 years ago)

Author:

mamuller

Message:

Added the static head contribution to the OutletLiquid? Pressure for the Dynamic Flash Vessel model.

File:

• branches/gui/eml/stage_separators/flash.mso (modified) (4 diffs)

branches/gui/eml/stage_separators/flash.mso

@@ -62 +62 @@
         Mw(NComp)               as molweight    (Brief="Mol Weight", Hidden=true);
         pi                      as positive             (Brief="Pi value", Default=3.141593,Hidden=true, Symbol="\pi");
+        g                               as acceleration         (Brief="Gravity Acceleration",Default=9.81,Hidden=true);
         
         Orientation     as Switcher     (Valid=["vertical","horizontal"],Default="vertical");
@@ -85 +86 @@
         Vcylinder = 0.25*(pi*Diameter^2)*Lenght;
         radius = 0.5*Diameter;
+        g = 9.81 * 'm/(s^2)';
 
 VARIABLES
@@ -109 +111 @@
         Pdrop           as press_delta  (Brief = "Pressure Drop", DisplayUnit = 'kPa', Symbol ="\Delta P", Protected=true);
 
+        Peq                     as pressure             (Brief="Equilibrium pressure on the liquid surface", Protected=true, DisplayUnit = 'kPa', Symbol="\Delta P_ª{eq}");
+        Pstatic         as pressure             (Brief="Static head at the bottom of the tank", Protected = true, DisplayUnit = 'kPa', Symbol="P_{static}^{Liquid}");
+
 out     TI as control_signal    (Brief="Temperature Indicator", PosX=1, PosY=0.39, Protected=true);
 out     PI as control_signal    (Brief="Pressure Indicator", PosX=1, PosY=0.21, Protected=true);
@@ -236 +241 @@
 
 "Liquid Volume"
-        vL = PP.LiquidVolume(OutletLiquid.T, OutletLiquid.P, OutletLiquid.z);
+        vL = PP.LiquidVolume(OutletLiquid.T, Peq, OutletLiquid.z);
 
 "Vapour Volume"
-        vV = PP.VapourVolume(OutletVapour.T, OutletVapour.P, OutletVapour.z);
+        vV = PP.VapourVolume(OutletVapour.T, Peq, OutletVapour.z);
         
 "Chemical Equilibrium"
-        PP.LiquidFugacityCoefficient(OutletLiquid.T, OutletLiquid.P, OutletLiquid.z)*OutletLiquid.z = 
-                PP.VapourFugacityCoefficient(OutletVapour.T, OutletVapour.P, OutletVapour.z)*OutletVapour.z;
+        PP.LiquidFugacityCoefficient(OutletLiquid.T, Peq, OutletLiquid.z)*OutletLiquid.z = 
+                PP.VapourFugacityCoefficient(OutletVapour.T, Peq, OutletVapour.z)*OutletVapour.z;
         
 "Thermal Equilibrium"
         OutletVapour.T = OutletLiquid.T;
-        
-"Mechanical Equilibrium"
-        OutletVapour.P = OutletLiquid.P;
+
+"Mechanical Equilibrium for the Vapour Phase"
+        OutletVapour.P = Peq;
+        
+"Static Head"   
+        Pstatic = PP.LiquidDensity(OutletLiquid.T, Peq, OutletLiquid.z) * g * Level;
+
+"Mechanical Equilibrium for the Liquid Phase"
+        OutletLiquid.P = Peq + Pstatic; 
 
 "Pressure Drop"