Changeset 831 for branches/gui/eml/stage_separators/tank.mso

Timestamp:

Aug 14, 2009, 5:20:34 PM (13 years ago)

Author:

gerson bicca

Message:

updates

File:

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

branches/gui/eml/stage_separators/tank.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");
@@ -73 +75 @@
         radius                                  as length               (Brief="Vessel radius",Hidden=true, Symbol="R_{cylinder}");
         
+        low_flow        as flow_mol     (Brief = "Low Flow",Default = 1E-6, Hidden=true);
+        zero_flow       as flow_mol     (Brief = "No Flow",Default = 0, Hidden=true);
+        KfConst         as Real                 (Brief="Constant for K factor pressure drop", Unit= 'mol/(s*(Pa^0.5))',Default = 1, Hidden=true);
+        Kfactor as positive (Brief="K factor for pressure drop", Lower = 1E-8, Default = 2);
+        
+        NormalFlow      as Switcher     (Brief="Normal Flow", Valid = ["on", "off"], Default = "on",Hidden=true);
+
         Levelpercent_Initial                    as positive     (Brief="Initial liquid height in Percent", Default = 0.70);
         Temperature_Initial                             as temperature  (Brief="Initial Liquid Temperature", Default = 330);
@@ -80 +89 @@
 
         Mw=PP.MolecularWeight();
-
+        g = 9.81 * 'm/(s^2)';
         Vhead_elliptical        = (pi*Diameter^3)/12;
         Vhead_hemispherical = (pi*Diameter^3)/6;
         Vcylinder = 0.25*(pi*Diameter^2)*Lenght;
         radius = 0.5*Diameter;
+        
+        low_flow = 1E-6 * 'kmol/h';
+        zero_flow = 0 * 'kmol/h';
+        KfConst = 1*'mol/(s*(Pa^0.5))';
 
 VARIABLES
@@ -105 +118 @@
         Level           as length               (Brief="liquid height", Protected=true);
         Across          as area                 (Brief="Vessel cylinder shell Cross section area", Hidden=true, Symbol="A_{cross}");
-        Pratio          as positive             (Brief = "Pressure Ratio", Symbol ="P_{ratio}", Protected=true);        
         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, Symbol="\Delta P_{eq}");
+        Pstatic         as pressure             (Brief="Static head at the bottom of the tank", Protected = true, Symbol="P_{static}^{Liquid}");
 
 out     TI as control_signal    (Brief="Temperature Indicator", PosX=0.525, PosY=0, Protected=true);
@@ -213 +227 @@
 end
 
+switch NormalFlow
+
+case "on":
+        Inlet.F = Kfactor *sqrt(Pdrop)*KfConst;
+
+        when Inlet.F < low_flow switchto "off";
+
+case "off":
+        Inlet.F = zero_flow;
+
+        when Inlet.P > OutletVapour.P switchto "on";
+
+end
+
 "Component Molar Balance"
         diff(TotalHoldup)=Inlet.F*Inlet.z - OutletLiquid.F*OutletLiquid.z - OutletVapour.F*OutletVapour.z;
@@ -232 +260 @@
 
 "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"
@@ -245 +273 @@
         
 "Mechanical Equilibrium"
-        OutletVapour.P = OutletLiquid.P;
+        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  = Inlet.P - Pdrop;
-
-"Pressure Ratio"
-        OutletLiquid.P = Inlet.P * Pratio;
+        Pdrop = Inlet.P - OutletVapour.P;
 
 "Geometry Constraint"
@@ -323 +354 @@
         radius                                  as length               (Brief="Vessel radius",Hidden=true, Symbol="R_{cylinder}");
         
+        low_flow        as flow_mol     (Brief = "Low Flow",Default = 1E-6, Hidden=true);
+        zero_flow       as flow_mol     (Brief = "No Flow",Default = 0, Hidden=true);
+        KfConst         as Real                 (Brief="Constant for K factor pressure drop", Unit= 'mol/(s*(Pa^0.5))',Default = 1, Hidden=true);
+        Kfactor as positive (Brief="K factor for pressure drop", Lower = 1E-8, Default = 2);
+        
+        NormalFlow      as Switcher     (Brief="Normal Flow", Valid = ["on", "off"], Default = "on",Hidden=true);
+
+
         Levelpercent_Initial                    as positive     (Brief="Initial liquid height in Percent", Default = 0.70);
         Temperature_Initial                             as temperature  (Brief="Initial Liquid Temperature", Default = 330);
@@ -333 +372 @@
         Vcylinder = 0.25*(pi*Diameter^2)*Lenght;
         radius = 0.5*Diameter;
+        
+        low_flow = 1E-6 * 'kmol/h';
+        zero_flow = 0 * 'kmol/h';
+        KfConst = 1*'mol/(s*(Pa^0.5))';
 
 VARIABLES
@@ -457 +500 @@
 end
 
+#*
+switch NormalFlow
+
+case "on":
+        Inlet.F = Kfactor *sqrt(Pstatic)*KfConst;
+
+        when Inlet.F < low_flow switchto "off";
+
+case "off":
+        Inlet.F = zero_flow;
+
+        when Inlet.P > OutletLiquid.P switchto "on";
+
+end
+*#
 "Component Molar Balance"
         diff(TotalHoldup)=Inlet.F*Inlet.z - OutletLiquid.F*OutletLiquid.z;