Changeset 470 for branches/packed/sample/stage_separators/sample_tray.mso

Timestamp:

Mar 4, 2008, 6:20:36 PM (15 years ago)

Author:

Paula Bettio Staudt

Message:

More changes

File:

• branches/packed/sample/stage_separators/sample_tray.mso (modified) (5 diffs)

branches/packed/sample/stage_separators/sample_tray.mso

@@ -112 +112 @@
         VARIABLES
         deltaP as pressure (Lower=-10);
-        phiL as Real;
+        phiL as volume;
         
         DEVICES
@@ -127 +127 @@
         EQUATIONS
         deltaP = t1.InletV.P - t1.OutletV.P;
-        phiL = t1.ML*t1.vL / t1.V;
+        phiL = t1.ML*t1.vL;
         
         SPECIFY
@@ -257 +257 @@
         VARIABLES
         deltaP as Real (Unit='atm/m');
+        phiL as volume;
         
         DEVICES
@@ -271 +272 @@
         EQUATIONS
         deltaP = (t1.InletV.P - t1.OutletV.P)/t1.hs;
+        phiL = t1.vL * t1.ML;
 
         SPECIFY
         feed.Outlet.F = 0 * 'kmol/h';
         feed.Outlet.T = 300 * 'K';
-        feed.Outlet.P = 1 * 'atm';
-        feed.Outlet.z = 1/NComp;
+        feed.Outlet.P = 1.66 * 'atm';
+        feed.Outlet.z = [0.2, 0.2, 0.2, 0.2, 0.2];
 
         inL.F = 71.21 * 'kmol/h';       
@@ -293 +295 @@
         #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.8 * 'm^2';
+        #t1.V = 0.8 * 'm^2' * 0.4 * 'm';
+#       t1.ds = 1.009 * 'm';
         t1.Cpo = 0.763;
         t1.e = 0.951;
         t1.a = 112.6 * 'm^2/m^3';
-        t1.hs = 0.4 * 'm';
+#       t1.hs = 0.4 * 'm';
         t1.Qsio = 1;
 
+        t1.V = 4 * 'ft^2' * 1 * 'ft';
+        t1.Ap = 4 * 'ft^2';
+        t1.hs = 1 * 'ft';
+        t1.ds = 2.26 * 'ft';
+
         INITIAL
-        t1.OutletL.T = 280 *'K';
-        t1.ML = 0.05 * 'kmol';
-        t1.OutletL.z([1:4]) = [0.286, 0.203, 0.148, 0.006];
+        t1.OutletL.T = 304 *'K';
+        t1.ML = 0.18 * 'kmol';
+        t1.OutletL.z([1:4]) = [0.2, 0.2, 0.2, 0.2];
         
         OPTIONS
         DAESolver(File="sundials");
         #InitialFile = "packedStage_BilletSchultes_Test.rlt";
-        TimeStep = 2;   
-        TimeEnd = 50;
+        InitialFile = "/home/paula/tray_Test_2.rlt";
+        TimeStep = 5;   
+        TimeEnd = 60;
 end