Changeset 782

Timestamp:

Jun 30, 2009, 10:53:56 AM (14 years ago)

Author:

gerson bicca

Message:

updated tank model

File:

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

branches/gui/eml/stage_separators/tank.mso

@@ -125 +125 @@
 * the Inlet stream;
 * the outlet flow;
-* the tank Q.
+* the InletQ (requires an energy source).
 
 == Initial Conditions ==
-* the tank initial temperature (OutletL.T);
-* the tank initial level (Level);
-* (NoComps - 1) OutletL (OR OutletV) compositions.
+* 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;
-        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");
+        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
-        Vtotal = (0.25*pi*(2*radius)^2)*L;
+        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}"); 
-        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");
-        E                               as energy                               (Brief="Total Energy Holdup on tank");
-        vL                      as volume_mol           (Brief="Liquid Molar Volume");
+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
 
-        Level                                   = Initial_Level;
-        Outlet.T                                = Initial_Temperature;
-        Outlet.z(1:NComp-1)     = Initial_Composition(1:NComp-1)/sum(Initial_Composition);
+"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;
@@ -186 +197 @@
         
 "Cylindrical Area"
-        Across = radius^2 * (asin(1) - asin((radius-Level)/radius) ) + (Level-radius)*sqrt(Level*(2*radius - Level)+1E-8);
+        Across = radius^2 * (asin(1) - asin((radius-Level)/radius) ) + (Level-radius)*sqrt(Level*(2*radius - Level)+eps*'m^2');
 
 "Level of liquid phase"
@@ -194 +205 @@
         Vfilled = L*Across;
 
-end
-
-Model tank_simplified
-        ATTRIBUTES
-        Pallete         = true;
-        Icon            = "icon/Tank"; 
-        Brief           = "Model of a simplified tank.";
-        Info            =
-"== Specify ==
-* the Inlet flow rate;
-
-== Initial Conditions ==
-* the tank initial level (Level);
-";
-
-        PARAMETERS
-        k as Real (Brief="Valve Constant", Unit = 'm^2.5/h', Default=4);
-        A as area (Brief="Tank area", Default=2);
-
-        VARIABLES
-        Level as length(Brief="Tank level");
-in      Fin  as flow_vol(Brief="Input flow", PosX=0.3037, PosY=0);
-out     Fout as flow_vol(Brief="Output flow", PosX=1, PosY=1);
-
-        EQUATIONS
-        "Mass balance"
-        diff(A*Level) = Fin - Fout;
-
-        "Valve equation"
-        Fout = k*sqrt(Level);           
 end