Changeset 820 for branches/gui

Timestamp:

Aug 10, 2009, 4:55:23 PM (14 years ago)

Author:

mamuller

Message:

Added the static head contribution to the liquid outlet pressure.

File:

• branches/gui/eml/stage_separators/tank.mso (modified) (5 diffs)

branches/gui/eml/stage_separators/tank.mso

@@ -583 +583 @@
         Pallete         = true;
         Icon            = "icon/SumpTank"; 
-        Brief           = "Model of a Tank With Thermodynamic Equilibrium.";
+        Brief           = "Model of a tank with 2 material inlet streams and with thermodynamic equilibrium.";
         Info            =
 "== ASSUMPTIONS ==
@@ -621 +621 @@
         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);
         
         Heads                   as Switcher     (Valid=["elliptical","hemispherical","flat"],Default="flat");
@@ -639 +640 @@
         Mw=PP.MolecularWeight();
 
+        g = 9.81 * 'm/(s^2)';
         Vhead_elliptical        = (pi*Diameter^3)/12;
         Vhead_hemispherical = (pi*Diameter^3)/6;
@@ -666 +668 @@
         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     LI as control_signal    (Brief="Level Indicator", PosX=1, PosY=0.7, Protected=true);
 out     TI as control_signal    (Brief="Temperature Indicator", PosX=1, PosY=0.6, Protected=true);
@@ -750 +755 @@
 
 "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"
-        #OutletLiquid.P  = min([InletLiquid.P,InletVapour.P]) - Pdrop;
-        OutletLiquid.P  = InletLiquid.P - Pdrop;
+        Pdrop = min([InletLiquid.P,InletVapour.P]) - OutletLiquid.P;
+        #OutletLiquid.P  = InletLiquid.P - Pdrop;
 
 "Geometry Constraint"