Changeset 461

Timestamp:

Feb 15, 2008, 1:00:49 PM (15 years ago)

Author:

Paula Bettio Staudt

Message:

More changes

Location:

branches/packed

Files:

• eml/stage_separators/tray.mso (modified) (2 diffs)
• sample/stage_separators/sample_tray.mso (modified) (1 diff)

branches/packed/eml/stage_separators/tray.mso

@@ -364 +364 @@
         rhoL as dens_mass;
         rhoV as dens_mass;
-        uL as velocity (Brief="volume flow rate of liquid, m^3/m^2/s", Default = 0.007);
-        uV as velocity (Brief="volume flow rate of vapor, m^3/m^2/s", Default = 1.14);
+        uL as velocity (Brief="volume flow rate of liquid, m^3/m^2/s", Lower=0, Default = 0.007);
+        uV as velocity (Brief="volume flow rate of vapor, m^3/m^2/s", Lower=0, Default = 1.14);
         dp as length (Brief="Particle diameter", Default=1e-3);
         invK as Real (Brief="Wall factor");
@@ -400 +400 @@
         "Reynolds number of the vapor stream"
         Rev*invK = dp*uV*rhoV / (miV*(1-e));
-        
-        "Resistance Coefficient"
-        Qsio = Cpo * (64/Rev + 1.8/Rev^0.08);
+
+        if Rev > 1e-4 then
+                "Resistance Coefficient"        
+                Qsio = 1;#Cpo * (64/Rev + 1.8/Rev^0.08);
+        else
+                Qsio = 1;
+        end
 
         "Pressure drop and Vapor flow"

branches/packed/sample/stage_separators/sample_tray.mso

@@ -144 +144 @@
         
         #t1.OutletV.F = 177.9 * 'kmol/h';
-        #t1.Qsio = 0.61;
 
         SET
         #Metal Pall Ring - nominal packing size 50 mm - Billet and Schultes, 1999.
         t1.Q = 0 * 'kW';
-        t1.Ap = 0.8 * 'm^2';
-        t1.V = 0.8 * 'm^2' * 0.4 * 'm';
-        t1.ds = 1.009 * 'm';
+        t1.Ap = 0.2 * 'm^2'; #0.8 * 'm^2';
+        t1.V = 0.2 * 'm^2' * 0.1 * 'm';#0.8 * 'm^2' * 0.4 * 'm';
+        t1.ds = 0.5 * 'm'; #1.009 * 'm';
         t1.Cpo = 0.763;
         t1.e = 0.951;
         t1.a = 112.6 * 'm^2/m^3';
-        t1.hs = 0.4 * 'm';#0.075 * 'm';
+        t1.hs = 0.1 * 'm'; #0.4 * 'm';
 
         INITIAL
-        t1.OutletL.T = 290 *'K';
-        t1.ML = 0.2 * 'kmol';
+        t1.OutletL.T = 280 *'K';
+        t1.ML = 0.01 * 'kmol';
         t1.OutletL.z([1:4]) = [0.286, 0.203, 0.148, 0.006];