Changeset 495 for branches/packed

Timestamp:

Apr 7, 2008, 3:42:56 PM (15 years ago)

Author:

Paula Bettio Staudt

Message:

More implementation

Location:

branches/packed

Files:

• eml/stage_separators/PackedStage.mso (modified) (7 diffs)
• sample/stage_separators/sample_column.mso (modified) (1 diff)

branches/packed/eml/stage_separators/PackedStage.mso

@@ -63 +63 @@
 #       QsiS as Real (Brief="Resistance coefficient at the loading point", Lower = 0);
 #       QsiFl as Real (Brief="Resistance coefficient at the flooding point", Lower = 0);
+        deltaP as pressure;
 
         SET
@@ -136 +137 @@
         Rev*invK = dp*uV*rhoV / (miV*(1-e));
         
-                if Rev < 1e-4 then
+#*      if Rev < 1e-4 then
                 Qsil = 1;
         else
@@ -143 +144 @@
                                 exp(13300/'m'^1.5/a^1.5*sqrt((uL^2*a)/g));
         end
-        
+*#      
+        deltaP = InletV.P - OutletV.P;
         "Pressure drop and Vapor flow"
-        (InletV.P - OutletV.P)/hs  = Qsil*a*uV^2*rhoV*invK / (2*(e-hl)^3);
+        deltaP/hs  = Qsil*a*uV^2*rhoV*invK / (2*(e-hl)^3);
 
         "Liquid holdup"
@@ -198 +200 @@
         
         VARIABLES
-        deltaP as Real (Unit='atm/m');
+        dP_h as Real (Unit='atm/m');
         
         DEVICES
@@ -212 +214 @@
 
         EQUATIONS
-        deltaP = (t1.InletV.P - t1.OutletV.P)/t1.hs;
+        dP_h = t1.deltaP/t1.hs;
         
         SPECIFY
@@ -230 +232 @@
         inV.z = [0.0584, 0.9416];#[0.5, 0.5];#
         
-        #t1.OutletV.P = 148 * 'kPa';
         t1.OutletV.F = 147.1 * 'kmol/h';
-        #t1.Qsil = 1;
+        #t1.deltaP = 0.05 * 'atm';
+        t1.Qsil = 1;
         
         SET
@@ -308 +310 @@
         #t1.OutletV.P = 2.30409 * 'atm';
         t1.OutletV.F = 165 * 'kmol/h';
+        #t1.deltaP = 0.0026 * 'atm';
         t1.Qsil = 10;
-#       t1.QsiS = 20;
-#       t1.QsiFl = 50;
         
         SET

branches/packed/sample/stage_separators/sample_column.mso

@@ -183 +183 @@
         
         OPTIONS
-        TimeStep = 10;
+        TimeStep = 1;
         TimeEnd = 100;
         #GuessFile="SectionColumn_Test_with8tray.rlt";