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

Timestamp:

Jul 15, 2009, 5:28:09 PM (13 years ago)

Author:

gerson bicca

Message:

updates

File:

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

branches/gui/eml/stage_separators/tank.mso

@@ -124 +124 @@
 end
 
-#*----------------------------------------------------------
-*
-*Model of a tank with a lain cylinder geometry
-*
-*---------------------------------------------------------*#
-Model tank_cylindrical
-        ATTRIBUTES
+Model RefluxDrum
+
+ATTRIBUTES
         Pallete         = true;
-        Icon            = "icon/TankHorizontal"; 
-        Brief           = "Model of a tank with a lain cylinder geometry.";
+        Icon            = "icon/RefluxDrum"; 
+        Brief           = "Model of a Reflux Drum with a lain cylinder geometry.";
         Info            =
 "== Specify ==
@@ -161 +157 @@
 
 SET
-        radius = Diameter/2;
+        radius = 0.5*Diameter;
 
 VARIABLES
@@ -220 +216 @@
 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