Changeset 869 for branches/gui

Timestamp:

Oct 29, 2009, 11:01:32 AM (14 years ago)

Author:

gerson bicca

Message:

added volume sizing for dynamic reboiler model

File:

• branches/gui/eml/stage_separators/reboiler.mso (modified) (8 diffs)

branches/gui/eml/stage_separators/reboiler.mso

@@ -343 +343 @@
 * thermodynamics equilibrium;
 * no liquid entrainment in the vapour stream.
-        
+
+== SET ==
+*Orientation: vessel position - vertical or horizontal;
+*Heads (bottom and top heads are identical)
+**elliptical: 2:1 elliptical heads (25% of vessel diameter);
+**hemispherical: hemispherical heads (50% of vessel diameter);
+*Diameter: Vessel diameter;
+*Lenght: Side length of the cylinder shell;
+
 == SPECIFY ==
 * the InletLiquid stream;
@@ -364 +372 @@
         outer PP                as Plugin       (Brief = "External Physical Properties", Type="PP");
         outer NComp     as Integer      (Brief="Number of Components");
-        Across                  as area         (Brief="Cross Section Area of reboiler");
-        V                               as volume       (Brief="Total volume of reboiler");
+        
+        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}");
         
         Levelpercent_Initial            as positive     (Brief="Initial liquid height in Percent", Default = 0.70);
@@ -371 +391 @@
         Initial_Composition(NComp)      as positive     (Brief="Initial Liquid Composition",Lower=1E-6);
 
+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;
+        
 VARIABLES
 
@@ -392 +420 @@
         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 reboiler",Protected=true);
+
+        
 INITIAL
 
@@ -405 +438 @@
 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
+
 "Component Molar Balance"
         diff(M)= InletLiquid.F*InletLiquid.z    - OutletLiquid.F*OutletLiquid.z - OutletVapour.F*OutletVapour.z;
@@ -415 +520 @@
 
 "Energy Holdup"
-        E = ML*OutletLiquid.h + MV*OutletVapour.h - OutletLiquid.P*V;
+        E = ML*OutletLiquid.h + MV*OutletVapour.h - OutletLiquid.P*Vtotal;
 
 "Mol Fraction Normalisation"
@@ -445 +550 @@
 
 "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"
@@ -457 +562 @@
 
 "Level indicator"
-        LI*V = Level*Across;
-
-end
+        LI*Vtotal= Vfilled;
+
+end