Changeset 451 for branches

Timestamp:

Jan 31, 2008, 3:10:36 PM (16 years ago)

Author:

Paula Bettio Staudt

Message:

Updated packedstage model

Location:

branches/packed

Files:

• eml/stage_separators/tray.mso (modified) (3 diffs)
• sample/stage_separators/sample_tray.mso (modified) (4 diffs)

branches/packed/eml/stage_separators/tray.mso

@@ -282 +282 @@
         g as acceleration;
         e as Real (Brief="Void fraction of packing, m^3/m^3");
-        V as volume;
         Across as area;
         ds as length (Brief="Column diameter");
         d as length (Brief="size of an element of packing");
-        h as length (Brief="Height of packing");
+        #h as length (Brief="Height of packing");
         C as Real (Brief="Constant for resitance factor equation");
         Cp as Real (Brief="Constant for resitance at loading point factor equation");
         Mw(NComp)       as molweight    (Brief = "Component Mol Weight");
 
-        SET
-        Mw = PP.MolecularWeight();
-        
         VARIABLES
         niL as viscosity (Brief="Liquid dynamic viscosity", DisplayUnit='kg/m/s');
@@ -309 +305 @@
         hs as length (Brief="Height of the packing stage");
         
+        SET
+        Mw = PP.MolecularWeight();
+        Ap = Across;
+        
         EQUATIONS
         "Liquid Density"
@@ -318 +318 @@
         "Vapour viscosity"
         niV = PP.VapourViscosity(InletV.T, InletV.P, InletV.z);
+
         "Liquid volume"
         VL = vL * ML;
         "Liquid holdup"
         hL = VL/V;
+
         "Volume flow rate of liquid, m^3/m^2/s"
         uL * Across = OutletL.F * vL;
-        "volume flow rate of vapor, m^3/m^2/s"
+        "Volume flow rate of vapor, m^3/m^2/s"
         uV * Across = OutletV.F * vV;
-
+        
         "Coefficient of Resistance"
         ksi * C^2 * (uL/uV * sqrt(rhoV/rhoL) * (niL/niV)^5.8)^3 = g/1*'s^2/m';
         ksil = Cp * (exp(uL*rhoL/a/niL/200)*(hL/hLs)^0.3) * (64/Rev+(1.8/Rev)) *
-                                ((e-hL)/e);#1.5
+                                ((e-hL)/e);#^1.5;
+        
         "Liquid holdup and Liquid flow"
         a^2 * niL * uL = hL^1 *(g*rhoL/3 - ksi*a*rhoV*uV^2/(4*hL*(e-hL)^2));
+
         "Pressure drop and Vapor flow"
         (InletV.P - OutletV.P)/hs = ksil *(a/2 + 2/ds)*(uV^2*rhoV/(e-hL)^3);
         
-        
+        "Vapor Reynolds number"
         Rev = uV * (d-2*hL/a) * rhoV/ niV;
-        
+
+        "Theoretical Liquid Holdup"
         hLs = (12*a^2*niL*uL/g/rhoL)^0.333;
 end

branches/packed/sample/stage_separators/sample_tray.mso

@@ -110 +110 @@
 
         VARIABLES
-        deltaP as Real (Unit='Pa/m');
+        deltaP as Real (Unit='atm/m');
         
         DEVICES
@@ -124 +124 @@
 
         EQUATIONS
-        deltaP = (t1.InletV.P - t1.OutletV.P)/ (t1.V/t1.Across);
+        deltaP = (t1.InletV.P - t1.OutletV.P)/t1.hs;
 
         SPECIFY
@@ -137 +137 @@
         inL.z = [0.1641, 0.8359];
 
-        #inV.F = 201.25 * 'kmol/h';
-        inV.F = 240 * 'kmol/h';
+        inV.F = 201.25 * 'kmol/h';
         inV.P = 1.48 * 'atm';
         inV.T = 321 * 'K';
@@ -144 +143 @@
 
         t1.Emv = 1;
+        t1.hs = 0.075 * 'm';
+#       t1.OutletV.F = 147.1 * 'kmol/h';
+#       t1.OutletL.F = 229.5 * 'kmol/h';
+#       t1.ksi = 0.805;
         
         SET
-#       t1.PressureDropModel = "Leva";
-#       t1.PackingType = "random";
-        t1.V = 4 * 'ft^3';
+        #Ceramic Pall Ring - nominal packing size 50 mm
         t1.Q = 0 * 'kW';
-        t1.Across = 1 * 'm^2';#3.94 * 'ft^2';
+        t1.Across = 0.8 * 'm^2';
+        t1.V = 0.8 * 'm^2' * 0.075 * 'm';
+        t1.ds = 1.009 * 'm';
+        t1.d = 50 * 'mm';
+        t1.C = 2.37;
+        t1.Cp = 0.662;
+        t1.e = 0.78;
+        t1.a = 120 * 'm^2/m^3';
 
         INITIAL
-        t1.OutletL.T = 320 *'K';
-        t1.Level = 1.2 * 0.125 * 'ft';
-        t1.OutletL.z(1) = 0.5;
+        t1.OutletL.T = 350 *'K';
+        t1.ML = 0.2 * 'mol';
+        t1.OutletL.z(1) = 0.1;
         
         OPTIONS
-        InitialFile = "/home/paula/tray_Test.rlt";
+        #InitialFile = "/home/paula/tray_Test.rlt";
         #DAESolver(File="dassl");
         TimeStep = 0.1;