Changeset 850 for branches/gui/eml/stage_separators/tray.mso

Timestamp:

Oct 8, 2009, 4:38:37 PM (13 years ago)

Author:

gerson bicca

Message:

packed column revised

File:

• branches/gui/eml/stage_separators/tray.mso (modified) (5 diffs)

branches/gui/eml/stage_separators/tray.mso

@@ -298 +298 @@
 ";      
         
-        
 PARAMETERS
 
 outer PP as Plugin(Brief = "External Physical Properties", Type="PP");
 outer NComp as Integer;
-#*      
-outer ColumnDiameter            as length               (Brief="Column diameter");
-outer AreaPerPackingVol as Real (Brief="Surface area per packing volume", Unit='m^2/m^3');
-outer VoidFraction as fraction (Brief="Void fraction of packing, m^3/m^3");
-outer hs as length (Brief="Height of the packing stage");
-outer ResistanceCoeff as positive (Brief="Resistance coefficient on the liquid load", Default = 0.6);
-*#
-        ColumnDiameter          as length               (Brief="Column diameter");
-        AreaPerPackingVol as Real (Brief="Surface area per packing volume", Unit='m^2/m^3');
-        VoidFraction as fraction (Brief="Void fraction of packing, m^3/m^3");
-        hs as length (Brief="Height of the packing stage");
-        ResistanceCoeff as positive (Brief="Resistance coefficient on the liquid load", Default = 0.6);
-
-        V as volume (Brief="Total volume of the tray");
-        HeatSupply      as heat_rate    (Brief="Rate of heat supply"); 
-        Gconst          as acceleration (Brief="Gravity Acceleration",Default=9.81,Hidden=true);
+
         Mw(NComp)       as molweight    (Brief = "Component Mol Weight");
-        VapourFlow as Switcher (Valid=["on", "off"], Default = "on");
-        
         
 VARIABLES
 
         Inlet                   as stream                               (Brief="Feed stream", PosX=0, PosY=0.4932, Symbol="_{in}");
-
-in      InletLiquid     as stream                               (Brief="Inlet liquid stream", PosX=0.5195, PosY=0, Symbol="_{inL}");
-in      InletVapour     as stream                               (Brief="Inlet vapour stream", PosX=0.4994, PosY=1, Symbol="_{inV}");
+in      InletLiquid                     as stream                               (Brief="Inlet liquid stream", PosX=0.5195, PosY=0, Symbol="_{inL}");
+in      InletVapour                     as stream                               (Brief="Inlet vapour stream", PosX=0.4994, PosY=1, Symbol="_{inV}");
 out     OutletLiquid    as liquid_stream        (Brief="Outlet liquid stream", PosX=0.8277, PosY=1, Symbol="_{outL}");
 out     OutletVapour    as vapour_stream        (Brief="Outlet vapour stream", PosX=0.8043, PosY=0, Symbol="_{outV}");
@@ -344 +325 @@
         rhoV as dens_mass;
         
-        deltaP as pressure (Lower = -10);
-        
-        uL as velocity (Brief="volume flow rate of liquid, m^3/m^2/s", Lower=0, Upper=100);
-        uV as velocity (Brief="volume flow rate of vapor, m^3/m^2/s", Lower=-0, Upper=100);
+        deltaP as pressure;
+        
+        uL as velocity (Brief="volume flow rate of liquid, m^3/m^2/s", Lower=-10, Upper=100);
+        uV as velocity (Brief="volume flow rate of vapor, m^3/m^2/s", Lower=-10, Upper=100);
         #dp as length (Brief="Particle diameter", Default=1e-3, Lower=0, Upper=10);
         #invK as positive (Brief="Wall factor", Default=1, Upper=10);
         #Rev as Real (Brief="Reynolds number of the vapor stream", Default=4000);
-        Al as area (Brief="Area occupied by the liquid", Default=0.001, Upper=10);
-        hl as positive (Brief="Column holdup", Unit='m^3/m^3', Default=0.04,Upper=1);
-
+        Al as area (Brief="Area occupied by the liquid", Default=0.001, Upper=1);
+        hl as positive (Brief="Column holdup", Unit='m^3/m^3', Default=0.01,Upper=10);
 
 SET
         Mw = PP.MolecularWeight();
 
-
 EQUATIONS
 
@@ -364 +343 @@
         diff(M)=Inlet.F*Inlet.z + InletLiquid.F*InletLiquid.z + InletVapour.F*InletVapour.z- OutletLiquid.F*OutletLiquid.z - OutletVapour.F*OutletVapour.z;
 
-"Energy Balance"
-    diff(E) = ( Inlet.F*Inlet.h + InletLiquid.F*InletLiquid.h + InletVapour.F*InletVapour.h - OutletLiquid.F*OutletLiquid.h - OutletVapour.F*OutletVapour.h + HeatSupply);
-
 "Molar Holdup"
         M = ML*OutletLiquid.z + MV*OutletVapour.z;
         
-"Energy Holdup"
-    E = ML*OutletLiquid.h + MV*OutletVapour.h - OutletLiquid.P*V;       
-        
 "Mol fraction normalisation"
         sum(OutletLiquid.z)= 1.0;
-        sum(OutletLiquid.z)=sum(OutletVapour.z);
+        sum(OutletLiquid.z)= sum(OutletVapour.z);
         
 "Liquid Volume"
@@ -392 +365 @@
         OutletLiquid.P = OutletVapour.P;
         
-"Geometry Constraint"
-    V*VoidFraction = ML*vL + MV*vV;
-
 "Liquid Density"
         rhoL = PP.LiquidDensity(OutletLiquid.T, OutletLiquid.P, OutletLiquid.z);
@@ -405 +375 @@
 
 #"Vapour viscosity"
-#       miV = PP.VapourViscosity(InletVapour.T, InletVapour.P, InletVapour.z);
-
-"Area occupied by the liquid"
-        Al = ML*vL/hs;
+ #      miV = PP.VapourViscosity(InletVapour.T, InletVapour.P, InletVapour.z);
 
 "Volume flow rate of liquid, m^3/m^2/s"
         uL * Al = OutletLiquid.F * vL;
         
-"Volume flow rate of vapour, m^3/m^2/s"
-        uV * (V*VoidFraction/hs - Al) = InletVapour.F * vV;
-        
-"Liquid holdup"
-    hl*V*VoidFraction = ML*vL;
-
-"Liquid velocity as a function of liquid holdup, Billet (4-27)"
-        hl^3 = (12/Gconst) * AreaPerPackingVol^2 * (miL/rhoL) * uL;
-               
-        switch VapourFlow
-                case "on":
-                        "Pressure drop and Vapor flow, Billet (4-58)"
-                                deltaP/hs  = ResistanceCoeff * (AreaPerPackingVol/2 + 2/ColumnDiameter) * 1/((VoidFraction-hl)^3) * (uV^2) * rhoV;
-    
-            when InletVapour.F < 1e-6 * 'kmol/h' switchto "off";
-
-
-        case "off":
-                        InletVapour.F = 0 * 'mol/s';
-                        when deltaP > 1e-4 * 'atm' switchto "on";
-    end
-        
-"Pressure profile"
-    deltaP = InletVapour.P - OutletVapour.P;
+        deltaP = InletVapour.P - OutletVapour.P;
         
 end