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

Timestamp:

Jul 18, 2009, 10:17:25 PM (14 years ago)

Author:

gerson bicca

Message:

updated tank model

File:

• branches/gui/eml/stage_separators/tank.mso (modified) (1 diff)

branches/gui/eml/stage_separators/tank.mso

@@ -362 +362 @@
 end
 
-Model tank_cylindrical
-        ATTRIBUTES
+Model TankL
+
+ATTRIBUTES
         Pallete         = true;
-        Icon            = "icon/TankHorizontal"; 
-        Brief           = "Model of a tank with a lain cylinder geometry.";
+        Icon            = "icon/TankL"; 
+        Brief           = "Model of a Tank.";
         Info            =
-"== Specify ==
+"== ASSUMPTIONS ==
+* liquid phase only;
+
+== 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 flow;
-* the InletQ (requires an energy source).
-
-== Initial Conditions ==
-* the tank initial temperature;
-* the tank initial level;
-* the tank initial composition.
+* the 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 TankL 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");
-        
-        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);
+outer PP                as Plugin       (Brief = "External Physical Properties", Type="PP");
+outer NComp     as Integer      (Brief = "Number of components", Lower = 1); 
+
+        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
-        radius = Diameter/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
-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);
+
+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}");
+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);
+        
+        E                       as energy               (Brief="Total Energy Holdup in the Vessel", Protected=true);
+        vL                      as volume_mol   (Brief="Liquid 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}");
+
+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" 
-        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);
+"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
 
-"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;
+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;
+        
+"Energy Balance"
+        diff(E) = Inlet.F*Inlet.h - OutletLiquid.F*OutletLiquid.h + InletQ;
 
 "Energy Holdup"
-        E = sum(M)*Outlet.h;
-        
+        E = sum(TotalHoldup)*OutletLiquid.h;
+
 "Mechanical Equilibrium"
-        Inlet.P = Outlet.P;
-        
+        Inlet.P = OutletLiquid.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;
+        vL = PP.LiquidVolume(OutletLiquid.T, OutletLiquid.P, OutletLiquid.z);
+
+"Molar Holdup"
+        TotalHoldup = OutletLiquid.z*sum(TotalHoldup); 
+        
+"Liquid Level"
+        Vfilled = sum(TotalHoldup) * vL;
+        
+"Temperature indicator"
+        TI * 'K' = OutletLiquid.T;
+
+"Pressure indicator"
+        PI * 'atm' = OutletLiquid.P;
+
+"Level indicator"
+        LI*Vtotal= Vfilled;
 
 end