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

Timestamp:

Jul 18, 2009, 5:08:06 PM (14 years ago)

Author:

gerson bicca

Message:

updated tank model (vertical/horizontal orientation with heads flat,elliptical or hemispherical)

File:

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

branches/gui/eml/stage_separators/tank.mso

@@ -12 +12 @@
 * EMSO is distributed under the therms of the ALSOC LICENSE as
 * available at http://www.enq.ufrgs.br/alsoc.
-*
-*-------------------------------------------------------------------
-* Model of tanks
-*-------------------------------------------------------------------- 
-*       Streams:
-*               * an inlet stream
-*               * an outlet stream
-*
-*       Specify:
-*               * the Inlet stream
-*               * the Outlet flow
-*               * the tank Q
-*
-*       Initial:
-*               * the tank temperature (OutletL.T)
-*               * the tank level (h)
-*               * (NoComps - 1) Outlet compositions
 *----------------------------------------------------------------------
 * Author: Paula B. Staudt
@@ -124 +107 @@
 end
 
-Model RefluxDrum
+Model TankVL
 
 ATTRIBUTES
         Pallete         = true;
-        Icon            = "icon/RefluxDrum"; 
-        Brief           = "Model of a Reflux Drum with a lain cylinder geometry.";
+        Icon            = "icon/TankVL"; 
+        Brief           = "Model of a Tank With Thermodynamic Equilibrium.";
+        Info            =
+"== ASSUMPTIONS ==
+* perfect mixing of both phases;
+* thermodynamics equilibrium.
+
+== 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);
+**flat: flat heads (0% of vessel diameter);
+*Diameter: Vessel diameter;
+*Lenght: Side length of the cylinder shell;
+        
+== SPECIFY ==
+* the Inlet stream;
+* the outlet flows: OutletVapour.F and OutletLiquid.F;
+* the InletQ (the model requires an energy stream, also you can use a controller for setting the heat duty using the heat_flow model).
+
+== OPTIONAL ==
+* the TankVL model has three control ports
+** TI OutletLiquid Temperature Indicator;
+** PI OutletLiquid Pressure Indicator;
+** LI Level Indicator;
+
+== INITIAL CONDITIONS ==
+* Initial_Temperature :  the Tank temperature (OutletLiquid.T);
+* Initial_Level : the Tank liquid level (Level);
+* Initial_Composition : (NoComps) OutletLiquid compositions.
+";
+        
+PARAMETERS
+outer PP                as Plugin       (Brief = "External Physical Properties", Type="PP");
+outer NComp     as Integer      (Brief = "Number of components", Lower = 1); 
+
+        Mw(NComp)               as molweight    (Brief="Mol Weight", Hidden=true);
+        pi                      as positive             (Brief="Pi value", Default=3.141593,Hidden=true, Symbol="\pi");
+        
+        Orientation     as Switcher     (Valid=["vertical","horizontal"],Default="vertical");
+        Heads                   as Switcher     (Valid=["elliptical","hemispherical","flat"],Default="flat");
+        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);
+        Temperature_Initial                             as temperature  (Brief="Initial Liquid Temperature", Default = 330);
+        Composition_Initial(NComp)              as fraction             (Brief="Initial Composition", Default = 0.10);
+
+SET
+
+        Mw=PP.MolecularWeight();
+
+        Vhead_elliptical        = (pi*Diameter^3)/12;
+        Vhead_hemispherical = (pi*Diameter^3)/6;
+        Vcylinder = 0.25*(pi*Diameter^2)*Lenght;
+        radius = 0.5*Diameter;
+
+VARIABLES
+
+in      Inlet                   as stream                       (Brief="Feed Stream", PosX=0.22, PosY=0, Symbol="_{in}");
+out     OutletLiquid    as liquid_stream        (Brief="Liquid outlet stream", PosX=0.43, PosY=1, Symbol="_{out}^{Liquid}");
+out     OutletVapour    as vapour_stream        (Brief="Vapour outlet stream", PosX=0.68, PosY=0, Symbol="_{out}^{Vapour}");
+in      InletQ                  as power                        (Brief="Heat Duty", PosX=0.735, PosY=1, Protected =true,Symbol="Q_{in}");
+
+        Vtotal                  as volume                       (Brief="Vessel total volume",Protected=true, Symbol="V_{total}");
+        Vfilled                 as volume                       (Brief="Vessel volume content",Protected=true, Symbol="V_{filled}");
+
+        TotalHoldup(NComp)              as mol  (Brief="Molar Holdup in the Vessel", Protected=true);
+        LiquidHoldup                    as mol  (Brief="Molar liquid holdup", Protected=true);
+        VapourHoldup                    as mol  (Brief="Molar vapour holdup", Protected=true);
+        
+        E                       as energy               (Brief="Total Energy Holdup in the Vessel", Protected=true);
+        vL                      as volume_mol   (Brief="Liquid Molar Volume", Protected=true);
+        vV                      as volume_mol   (Brief="Vapour Molar volume", Protected=true);
+        Level           as length               (Brief="liquid height", Protected=true);
+        Across          as area                 (Brief="Vessel cylinder shell Cross section area", Hidden=true, Symbol="A_{cross}");
+        vfrac           as positive     (Brief="Vapourization fraction", Symbol="\phi", Protected=true);
+        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);
+
+out     TI as control_signal    (Brief="Temperature Indicator", PosX=0.525, PosY=0, Protected=true);
+out     PI as control_signal    (Brief="Pressure Indicator", PosX=0.368, PosY=0, Protected=true);
+out     LI as control_signal    (Brief="Level Indicator", PosX=1, PosY=0.6, Protected=true);
+
+INITIAL
+
+"Initial level Percent"
+        LI = Levelpercent_Initial;
+        
+"Initial Outlet Liquid Temperature"
+        OutletLiquid.T = Temperature_Initial;
+        
+"Initial Outlet Liquid Composition Normalized"
+        OutletLiquid.z(1:NComp - 1) = Composition_Initial(1:NComp - 1)/sum(Composition_Initial);
+
+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
+
+case "flat":
+
+"Vessel Total Volume"
+        Vtotal = Vcylinder;
+
+"Vessel Filled Volume"
+        Vfilled = Across*Level;
+
+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;
+
+case "flat":
+
+"Vessel Total Volume"
+        Vtotal = Vcylinder;
+
+"Vessel Filled Volume"
+        Vfilled = Across*Lenght;
+
+end
+
+end
+
+"Component Molar Balance"
+        diff(TotalHoldup)=Inlet.F*Inlet.z - OutletLiquid.F*OutletLiquid.z - OutletVapour.F*OutletVapour.z;
+        
+"Energy Balance"
+        diff(E) = Inlet.F*Inlet.h - OutletLiquid.F*OutletLiquid.h - OutletVapour.F*OutletVapour.h + InletQ;
+        
+"Molar Holdup"
+        TotalHoldup = LiquidHoldup*OutletLiquid.z + VapourHoldup*OutletVapour.z; 
+        
+"Energy Holdup"
+        E = LiquidHoldup*OutletLiquid.h + VapourHoldup*OutletVapour.h - OutletLiquid.P*Vtotal;
+        
+"Mol fraction normalisation"
+        sum(OutletLiquid.z)=1.0;
+
+"Mol fraction normalisation"
+        sum(OutletLiquid.z)=sum(OutletVapour.z);
+
+"Vaporization Fraction"
+        OutletVapour.F = Inlet.F * vfrac;
+
+"Liquid Volume"
+        vL = PP.LiquidVolume(OutletLiquid.T, OutletLiquid.P, OutletLiquid.z);
+
+"Vapour Volume"
+        vV = PP.VapourVolume(OutletVapour.T, OutletVapour.P, OutletVapour.z);
+        
+"Chemical Equilibrium"
+        PP.LiquidFugacityCoefficient(OutletLiquid.T, OutletLiquid.P, OutletLiquid.z)*OutletLiquid.z = 
+                PP.VapourFugacityCoefficient(OutletVapour.T, OutletVapour.P, OutletVapour.z)*OutletVapour.z;
+        
+"Thermal Equilibrium"
+        OutletVapour.T = OutletLiquid.T;
+        
+"Mechanical Equilibrium"
+        OutletVapour.P = OutletLiquid.P;
+
+"Pressure Drop"
+        OutletLiquid.P  = Inlet.P - Pdrop;
+
+"Pressure Ratio"
+        OutletLiquid.P = Inlet.P * Pratio;
+
+"Geometry Constraint"
+        Vtotal = LiquidHoldup * vL + VapourHoldup * vV;
+
+"Temperature indicator"
+        TI * 'K' = OutletLiquid.T;
+
+"Pressure indicator"
+        PI * 'atm' = OutletLiquid.P;
+
+"Level indicator"
+        LI*Vtotal= Vfilled;
+
+"Liquid Level"
+        LiquidHoldup * vL = Vfilled;
+
+end
+
+Model tank_cylindrical
+        ATTRIBUTES
+        Pallete         = true;
+        Icon            = "icon/TankHorizontal"; 
+        Brief           = "Model of a tank with a lain cylinder geometry.";
         Info            =
 "== Specify ==
@@ -157 +394 @@
 
 SET
-        radius = 0.5*Diameter;
+        radius = Diameter/2;
 
 VARIABLES
@@ -216 +453 @@
 end
 
-Model tank_cylindrical
-        ATTRIBUTES
-        Pallete         = true;
-        Icon            = "icon/TankHorizontal"; 
-        Brief           = "Model of a tank with a lain cylinder geometry.";
-        Info            =
-"== Specify ==
-* the Inlet stream;
-* the outlet flow;
-* the InletQ (requires an energy source).
-
-== Initial Conditions ==
-* the tank initial temperature;
-* the tank initial level;
-* the tank initial composition.
-";
-        
-PARAMETERS
-        outer PP                as Plugin               (Brief = "External Physical Properties", Type="PP");
-        outer NComp     as Integer              (Brief = "Number of Components");
-        
-        pi                      as positive     (Brief="Pi value", Default=3.141593,Hidden=true);
-        eps                     as positive     (Brief="small number",Default=1E-8,Hidden=true);
-        Diameter                as length       (Brief="Tank internal Diameter",Default=1.5);
-        radius                  as length       (Brief="Tank radius",Hidden=true);
-        L                               as length       (Brief="Tank length",Default=5);
-        
-        Initial_Level                           as length               (Brief="Initial Level of the Tank",Default=1);
-        Initial_Temperature                     as temperature  (Brief="Initial Temperature of Liquid",Default=300);
-        Initial_Composition(NComp)      as positive     (Brief="Initial Liquid Composition",Lower=1E-8,Default=0.1);
-
-SET
-        radius = Diameter/2;
-
-VARIABLES
-in      Inlet           as stream                       (Brief="Inlet stream", PosX=0.1825, PosY=0, Symbol="_{in}");
-out     Outlet          as liquid_stream        (Brief="Outlet liquid stream", PosX=1, PosY=1, Symbol="_{out}");
-in      InletQ          as power                        (Brief="Rate of heat supply", PosX=1, PosY=0.6160, Symbol="_{in}",Protected=true); 
-        Level           as length                       (Brief="Tank level",Protected=true);
-        Vtotal          as volume                       (Brief="Tank total volume",Protected=true);
-        Vfilled         as volume                       (Brief="Tank volume content",Protected=true);
-        Across          as area                         (Brief="Tank cross section area", Default=2,Protected=true);
-        E                       as energy                       (Brief="Total Energy Holdup on tank",Protected=true);
-        vL                      as volume_mol           (Brief="Liquid Molar Volume",Protected=true);
-        M(NComp)        as mol                          (Brief="Molar Holdup in the tank",Protected=true);
-
-INITIAL
-
-"Initial Level" 
-        Level = Initial_Level;
-
-"Initial Liquid Temperature"    
-        Outlet.T = Initial_Temperature;
-
-"Initial Liquid Composition"    
-        Outlet.z(1:NComp-1) = Initial_Composition(1:NComp-1)/sum(Initial_Composition);
-
-EQUATIONS
-
-"Tank total volume"
-        Vtotal = (0.25*pi*(Diameter)^2)*L;
-        
-"Mass balance"
-        diff(M) = Inlet.F*Inlet.z - Outlet.F*Outlet.z;
-        
-"Energy balance"
-        diff(E) = Inlet.F*Inlet.h - Outlet.F*Outlet.h + InletQ;
-
-"Energy Holdup"
-        E = sum(M)*Outlet.h;
-        
-"Mechanical Equilibrium"
-        Inlet.P = Outlet.P;
-        
-"Liquid Volume"
-        vL = PP.LiquidVolume(Outlet.T, Outlet.P, Outlet.z);
-        
-"Composition"
-        M = Outlet.z*sum(M);
-        
-"Cylindrical Area"
-        Across = radius^2 * (asin(1) - asin((radius-Level)/radius) ) + (Level-radius)*sqrt(Level*(2*radius - Level)+eps*'m^2');
-
-"Level of liquid phase"
-        L*Across = sum(M)*vL;
-
-"Volume Filled of liquid phase"
-        Vfilled = L*Across;
-
-end
-
 Model tank_feed
         ATTRIBUTES