Changeset 496 for branches/packed

Timestamp:

Apr 9, 2008, 4:42:26 PM (15 years ago)

Author:

Paula Bettio Staudt

Message:

Almost there!

Location:

branches/packed/eml

Files:

• stage_separators/PackedStage.mso (modified) (7 diffs)
• stage_separators/packedStage_Test_1.rlt (added)
• stage_separators/packedStage_Test_2.rlt (added)
• streams.mso (modified) (1 diff)

branches/packed/eml/stage_separators/PackedStage.mso

@@ -33 +33 @@
 out     OutletV as vapour_stream (Brief="Outlet vapour stream", PosX=0.8043, PosY=0, Symbol="_{outV}");
 
-        M(NComp) as mol (Brief="Molar Holdup in the tray");
-        ML as mol (Brief="Molar liquid holdup");
-        MV as mol (Brief="Molar vapour holdup");
-        E as energy (Brief="Total Energy Holdup on tray");
-        vL as volume_mol (Brief="Liquid Molar Volume");
-        vV as volume_mol (Brief="Vapour Molar volume");
+        M(NComp) as mol (Brief="Molar Holdup in the tray", Default=0.01, Lower=0, Upper=100);
+        ML as mol (Brief="Molar liquid holdup", Default=0.01, Lower=0, Upper=100);
+        MV as mol (Brief="Molar vapour holdup", Default=0.01, Lower=0, Upper=100);
+        E as energy (Brief="Total Energy Holdup on tray", Default=-500);
+        vL as volume_mol (Brief="Liquid Molar Volume", Lower=0);
+        vV as volume_mol (Brief="Vapour Molar volume", Lower=0);
         
         miL as viscosity (Brief="Liquid dynamic viscosity", DisplayUnit='kg/m/s');
@@ -47 +47 @@
 #       miw as viscosity (Brief="Water viscosity");
         
-        uL as velocity (Brief="volume flow rate of liquid, m^3/m^2/s", Lower = -10, Default = 0.007);
-        uV as velocity (Brief="volume flow rate of vapor, m^3/m^2/s", Lower = -10, Default = 1.14);
-        dp as length (Brief="Particle diameter", Default=1e-3);
-        invK as Real (Brief="Wall factor");
-        Rev as Real (Brief="Reynolds number of the vapor stream", Lower = 0, Default=100);
+        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");
+        Rev as Real (Brief="Reynolds number of the vapor stream", Default=100);
         Al as area;
         hl as Real (Brief="Column holdup", Unit='m^3/m^3');
@@ -151 +151 @@
         "Liquid holdup"
         hl = (12*miL*a^2*uL/rhoL/g)^1/3;
-
-        
 #*      
         "Liquid holdup at the loading point"
@@ -232 +230 @@
         inV.z = [0.0584, 0.9416];#[0.5, 0.5];#
         
-        t1.OutletV.F = 147.1 * 'kmol/h';
-        #t1.deltaP = 0.05 * 'atm';
+        #t1.OutletV.F = 180.1 * 'kmol/h';
+        t1.deltaP = 0.00189536 * 'atm';
         t1.Qsil = 1;
         
@@ -248 +246 @@
         
         INITIAL
-        t1.OutletL.T = 315 *'K';
-        t1.ML = 0.25 * 'kmol';
-        t1.OutletL.z(1) = 0.1641;
+        t1.OutletL.T = 340 *'K';#347.54 *'K';
+        t1.ML = 0.02 * 'kmol';#0.0131283 * 'kmol';
+        t1.OutletL.z(1) = 0.1;#0.112595;
         
         OPTIONS
+        #InitialFile = "packedStage_Test_1";
         DAESolver(File="sundials");
         TimeStep = 10;
@@ -309 +308 @@
         
         #t1.OutletV.P = 2.30409 * 'atm';
-        t1.OutletV.F = 165 * 'kmol/h';
-        #t1.deltaP = 0.0026 * 'atm';
+        #t1.OutletV.F = 165 * 'kmol/h';
+        t1.deltaP = 0.0214961 * 'atm';
         t1.Qsil = 10;
         
@@ -326 +325 @@
 
         INITIAL
-        t1.OutletL.T = 290 *'K';
-        t1.ML = 0.002 * 'kmol';
-        t1.OutletL.z([1:4]) = [0.226, 0.425, 0.035, 0.025];
+        t1.OutletL.T = 305.912 *'K';
+        t1.ML = 0.00327166 * 'kmol';
+        t1.OutletL.z([1:4]) = [0.169668, 0.511287, 0.0321714, 0.046279];
         
         OPTIONS
         DAESolver(File="sundials");
-        #InitialFile = "/home/paula/tray_Test_2.rlt";
+        #InitialFile = "packedStage_Test_2";
         TimeStep = 10;
         TimeEnd = 100;

branches/packed/eml/streams.mso

@@ -35 +35 @@
 
         VARIABLES
-        F as flow_mol                   (Brief = "Stream Molar Flow Rate");
+        F as flow_mol                   (Brief = "Stream Molar Flow Rate", Lower=-1e2);
         T as temperature                (Brief = "Stream Temperature");
         P as pressure                   (Brief = "Stream Pressure");