Changeset 879 for branches/gui/eml/stage_separators/reboiler.mso

Timestamp:

Nov 9, 2009, 2:11:31 PM (13 years ago)

Author:

gerson bicca

Message:

updated dynamic models

File:

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

branches/gui/eml/stage_separators/reboiler.mso

@@ -18 +18 @@
 *--------------------------------------------------------------------*#
 
-using "streams";
+using "tank";
 
 Model thermosyphon
@@ -373 +373 @@
         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}");
-        
         Levelpercent_Initial            as positive     (Brief="Initial liquid height in Percent", Default = 0.70);
         Initial_Temperature                     as temperature  (Brief="Initial Temperature of Reboiler");
         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
+
+        Geometry                as VesselVolume (Brief="Vessel Geometry", Symbol=" ");
 
 in      InletLiquid     as stream                       (Brief="Liquid inlet stream", PosX=0.17, PosY=1, Symbol="_{in}^{Liquid}");
@@ -417 +397 @@
         vV                      as volume_mol   (Brief="Vapour Molar volume", Protected = true);
         rhoV            as dens_mass    (Brief="Vapour Density", Protected = true, Symbol="\rho");
-        Level           as length               (Brief="Level of liquid phase", Protected = true);
         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
 
@@ -438 +412 @@
 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;
@@ -520 +422 @@
 
 "Energy Holdup"
-        E = ML*OutletLiquid.h + MV*OutletVapour.h - OutletLiquid.P*Vtotal;
+        E = ML*OutletLiquid.h + MV*OutletVapour.h - OutletLiquid.P*Geometry.Vtotal;
 
 "Mol Fraction Normalisation"
@@ -550 +452 @@
 
 "Geometry Constraint"
-        Vtotal = ML*vL + MV*vV;
+        Geometry.Vtotal = ML*vL + MV*vV;
 
 "Liquid Level"
-        ML * vL = Vfilled;
+        ML * vL = Geometry.Vfilled;
 
 "Temperature Indicator"
@@ -562 +464 @@
 
 "Level indicator"
-        LI*Vtotal= Vfilled;
+        LI*Geometry.Vtotal= Geometry.Vfilled;
 
 end