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

Timestamp:

Jun 25, 2009, 6:05:45 PM (13 years ago)

Author:

gerson bicca

Message:

updated column/tank model

File:

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

branches/gui/eml/stage_separators/tank.mso

@@ -57 +57 @@
         outer NComp     as Integer;
         Across                          as area                         (Brief="Tank cross section area", Default=2);
+        L                                       as length                       (Brief="Tank length");
+        Vtotal                          as volume       (Brief="Tank total volume");
         
         Initial_Level                           as length               (Brief="Initial Level of the Tank");
         Initial_Temperature                     as temperature  (Brief="Initial Temperature of Liquid");
         Initial_Composition(NComp)      as fraction     (Brief="Initial Liquid Composition");
+
+SET
+        Vtotal =Across*L;
         
 VARIABLES
@@ -66 +71 @@
 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.7859, Symbol="_{in}"); 
+        Vfilled         as volume                               (Brief="Tank volume content");
         Level                   as length                               (Brief="Tank level");
         M(NComp)        as mol                                  (Brief="Molar Holdup in the tank");
@@ -100 +106 @@
         Level = sum(M)*vL/Across;
 
+"Volume Filled of liquid phase"
+        Vfilled = Level*Across;
+        
 end
 
@@ -127 +136 @@
         outer PP                        as Plugin               (Brief = "External Physical Properties", Type="PP");
         outer NComp     as Integer;
+        pi                                      as positive     (Brief="Pi value", Default=3.141593,Hidden=true);
         radius                          as length               (Brief="Tank radius");
         L                                               as length               (Brief="Tank length");
+        Vtotal                          as volume       (Brief="Tank total volume");
         
         Initial_Level                   as length                       (Brief="Initial Level of the Tank");
         Initial_Temperature             as temperature                  (Brief="Initial Temperature of Liquid");
         Initial_Composition(NComp)      as fraction                             (Brief="Initial Liquid Composition");
-        
+
+SET
+        Vtotal = (0.25*pi*(2*radius)^2)*L;
+
 VARIABLES
 in      Inlet           as stream                               (Brief = "Inlet stream", PosX=0.1825, PosY=0, Symbol="_{in}");
@@ -139 +153 @@
 in      InletQ          as power                                (Brief="Rate of heat supply", PosX=1, PosY=0.6160, Symbol="_{in}"); 
         Level           as length                               (Brief="Tank level");
+        Vfilled         as volume                               (Brief="Tank volume content");
         Across          as area                                         (Brief="Tank cross section area", Default=2);
         M(NComp) as mol                                         (Brief="Molar Holdup in the tank");
@@ -171 +186 @@
         
 "Cylindrical Area"
-        Across = radius^2 * (asin(1) - asin((radius-Level)/radius) ) + (Level-radius)*sqrt(Level*(2*radius - Level));
+        Across = radius^2 * (asin(1) - asin((radius-Level)/radius) ) + (Level-radius)*sqrt(Level*(2*radius - Level)+1E-8);
 
 "Level of liquid phase"
         L*Across = sum(M)*vL;
+
+"Volume Filled of liquid phase"
+        Vfilled = L*Across;
 
 end