Changeset 872 for branches/gui/eml

Timestamp:

Oct 29, 2009, 12:23:40 PM (13 years ago)

Author:

gerson bicca

Message:

added volume sizing for dynamic condenser model

File:

• branches/gui/eml/stage_separators/condenser.mso (modified) (7 diffs)

branches/gui/eml/stage_separators/condenser.mso

@@ -303 +303 @@
         outer NComp     as Integer (Brief="Number of Components");
         
+        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");
+        Heads                   as Switcher (Valid=["elliptical","hemispherical"],Default="elliptical");
+
+        Diameter                as length       (Brief="Vessel diameter", Symbol="D_{i}");
+        Lenght                  as length       (Brief="Side length of the cylinder shell", Symbol="L_{vessel}");
+        
+        Vhead_elliptical                as volume       (Brief="Elliptical Head Total Volume",Hidden=true, Symbol="V_{head}^{elliptical}");
+        Vhead_hemispherical     as volume       (Brief="Hemispherical Head Total Volume",Hidden=true, Symbol="V_{head}^{hemispherical}");
+        Vcylinder                               as volume       (Brief="Cylinder Total Volume",Hidden=true, Symbol="V_{cylinder}");
+        radius                                  as length       (Brief="Vessel radius",Hidden=true, Symbol="R_{cylinder}");
+        
         Mw(NComp)       as molweight    (Brief = "Component Mol Weight",Hidden=true);
         low_flow        as flow_mol     (Brief = "Low Flow",Default = 1E-6, Hidden=true);
@@ -310 +324 @@
         VapourFlow      as Switcher     (Brief="Vapour Flow", Valid = ["on", "off"], Default = "on",Hidden=true);
 
-        V               as volume       (Brief="Condenser total volume");
-        Across  as area         (Brief="Cross Section Area of condenser");
         Kfactor as positive (Brief="K factor for pressure drop", Lower = 1E-8, Default = 1E-3);
         
@@ -339 +351 @@
         Pdrop           as press_delta  (Brief = "Pressure Drop", DisplayUnit = 'kPa', Symbol ="\Delta P", Protected=true);
 
+        Vtotal                  as volume       (Brief="Vessel total volume",Protected=true, Symbol="V_{total}");
+        Vfilled                 as volume       (Brief="Vessel volume content",Protected=true, Symbol="V_{filled}");
+        Across                  as area         (Brief="Cross Section Area of condenser",Protected=true);
+
 SET
+        pi = 3.141593;
+        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;
 
         Mw   = PP.MolecularWeight();
@@ -359 +381 @@
 EQUATIONS
 
+switch Orientation
+
+case "vertical":
+
+"Vessel Cross Section Area"
+        Across = 0.25*(pi*Diameter^2);
+
+switch Heads
+
+case "elliptical":
+
+"Vessel Total Volume"
+        Vtotal = Vhead_elliptical +     Vcylinder;
+
+if Level < 0.25*Diameter then
+
+"Vessel Filled Volume"
+        Vfilled = 0.25*pi*(((Diameter*Level)/(0.25*Diameter))^2)*(0.25*Diameter-Level/3);
+
+else
+
+"Vessel Filled Volume"
+        Vfilled = 0.25*pi*(Diameter^2)*(Level - 0.25*Diameter/3);
+
+end
+
+case "hemispherical":
+
+"Vessel Total Volume"
+        Vtotal = Vhead_hemispherical + Vcylinder;
+
+if Level < 0.5*Diameter then
+
+"Vessel Filled Volume"
+        Vfilled = 0.25*pi*(Level^2)*(2*Diameter-4*Level/3);
+
+else
+
+"Vessel Filled Volume"
+        Vfilled = 0.25*pi*((2/3)*((0.5*Diameter)^3) - (0.25*(Diameter)^3) + Level*Diameter^2);
+
+end
+
+end
+
+case "horizontal":
+
+"Vessel Cross Section Area"
+        Across = (radius^2)*acos((radius-Level)/radius)-(radius-Level)*sqrt((2*radius*Level-Level^2));
+
+switch Heads
+
+case "elliptical":
+
+"Vessel Total Volume"
+        Vtotal = Vhead_elliptical +     Vcylinder;
+
+"Vessel Filled Volume"
+        Vfilled = 0.5236*Level^2*(1.5*Diameter-Level) + Across*Lenght;
+
+case "hemispherical":
+
+"Vessel Total Volume"
+        Vtotal = Vhead_hemispherical + Vcylinder;
+
+"Vessel Filled Volume"
+        Vfilled = 1.0472*Level^2*(1.5*Diameter-Level) + Across*Lenght;
+
+end
+
+end
+
 switch VapourFlow
 
@@ -383 +477 @@
         
 "Energy Holdup"
-        E = ML*OutletLiquid.h + MV*OutletVapour.h - OutletVapour.P*V;
+        E = ML*OutletLiquid.h + MV*OutletVapour.h - OutletVapour.P*Vtotal;
         
 "Mol fraction normalisation"
@@ -414 +508 @@
 
 "Geometry Constraint"
-        V = ML*vL + MV*vV;
-
-"Level of liquid phase"
-        Level = ML*vL/Across;
+        Vtotal = ML*vL + MV*vV;
+
+"Liquid Level"
+        ML * vL = Vfilled;
 
 "Temperature indicator"
@@ -426 +520 @@
 
 "Level indicator"
-        LI*V = Level*Across;
-        
-end
+        LI*Vtotal= Vfilled;
+        
+end