Changeset 468 for branches/packed

Timestamp:

Mar 3, 2008, 3:49:53 PM (16 years ago)

Author:

Paula Bettio Staudt

Message:

More modifications

Location:

branches/packed

Files:

• eml/stage_separators/tray.mso (modified) (2 diffs)
• sample/stage_separators/sample_column.mso (modified) (5 diffs)
• sample/stage_separators/sample_tray.mso (modified) (4 diffs)

branches/packed/eml/stage_separators/tray.mso

@@ -369 +369 @@
         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);
+        #Rev as Real (Brief="Reynolds number of the vapor stream", Lower = 0, Default=100);
         #Qsio as Real (Brief="Resistance coefficient", Lower = 0);
         
@@ -399 +399 @@
         invK = (1 + (2*dp/(3*ds*(1-e))));
         
-        "Reynolds number of the vapor stream"
-        Rev*invK = dp*uV*rhoV / (miV*(1-e));
+        #"Reynolds number of the vapor stream"
+        #Rev*invK = dp*uV*rhoV / (miV*(1-e));
 
 #*      if Rev > 1e-4 then

branches/packed/sample/stage_separators/sample_column.mso

@@ -503 +503 @@
 end
 
-
+# section based on the Distillation_kettle_cond_Test
 FlowSheet SectionColumnPacked
         PARAMETERS
@@ -551 +551 @@
 
         reb.F = 137.57 * 'kmol/h';
-        reb.P = 7.46 * 'atm';
+        reb.P = 2.46 * 'atm';
         reb.T = 325 * 'K';
         reb.z = [0.16, 0.542, 0.013, 0.008, 0.277];
@@ -563 +563 @@
         sec.stage.Cpo = 0.763;
         sec.stage.e = 0.951;
-        sec.stage.a = 100*112.6 * 'm^2/m^3';
+        sec.stage.a = 112.6 * 'm^2/m^3';
         sec.stage.Qsio = 1;
 
@@ -574 +574 @@
         #InitialFile = "/home/paula/SectionColumn_Test_with8tray.rlt";
         DAESolver(File="dassl");
-        TimeStep = 0.01;
-        TimeEnd = 1;
+        TimeStep = 1;
+        TimeEnd = 10;
 end
 
@@ -635 +635 @@
         
         SET
-        col.H = 2 * 'm';
+        col.H = 1 * 'm';
         col.NStages = 2;
 

branches/packed/sample/stage_separators/sample_tray.mso

@@ -25 +25 @@
 using "stage_separators/tray";
 
-FlowSheet tray_Test
+FlowSheet tray_Test_1
         PARAMETERS
         PP      as Plugin(Brief="Physical Properties",
@@ -96 +96 @@
 end
 
+FlowSheet tray_Test_2
+        PARAMETERS
+        PP      as Plugin(Brief="Physical Properties",
+                Type="PP",
+                Components = [ "isobutane", "n-pentane", "propylene",
+                "benzene", "isobutene" ],
+                LiquidModel = "PR",
+                VapourModel = "PR"
+        );
+        NComp   as Integer;
+
+        SET
+        NComp = PP.NumberOfComponents;
+        
+        VARIABLES
+        deltaP as pressure;
+        phiL as Real;
+        
+        DEVICES
+        t1 as tray;
+        feed as source;
+        inL as liquid_stream;
+        inV as vapour_stream;
+        
+        CONNECTIONS
+        feed.Outlet to t1.Inlet;
+        inL to t1.InletL;
+        inV to t1.InletV;
+        
+        EQUATIONS
+        deltaP = t1.InletV.P - t1.OutletV.P;
+        phiL = t1.ML*t1.vL / t1.V;
+        
+        SPECIFY
+        feed.Outlet.F = 0 * 'kmol/h';
+        feed.Outlet.T = 300 * 'K';
+        feed.Outlet.P = 1.66 * 'atm';
+        feed.Outlet.z = [0.5, 0.5];
+        
+        inL.P = 2.22 * 'atm';
+        inL.T = 297.6 * 'K';
+        inL.F = 71.21 * 'kmol/h';
+        inL.z = [0.226, 0.425, 0.035, 0.025, 0.289];
+
+        inV.F = 175.3 * 'kmol/h';
+        inV.P = 2.3062 * 'atm';
+        inV.T = 308.3 * 'K';
+        inV.z = [0.265, 0.233, 0.150, 0.014, 0.338];
+
+        t1.Emv = 1;
+
+        SET
+        t1.V = 4 * 'ft^3';
+        t1.Ah = 0.394 * 'ft^2';
+        t1.lw = 20.94 * 'in';
+        t1.hw = 0.125 * 'ft';
+        t1.Q = 0 * 'kW';
+        t1.beta = 0.6;
+        t1.alfa = 4;
+        t1.Ap = 3.94 * 'ft^2';
+        
+        INITIAL
+        t1.OutletL.T = 304 *'K';
+        t1.Level = 1.2 * t1.hw;
+        t1.OutletL.z([1:4]) = [0.2, 0.2, 0.2, 0.2];
+        
+        OPTIONS
+        #InitialFile = "/home/paula/tray_Test.rlt";
+        TimeStep = 10;
+        TimeEnd = 100;
+end
+
+#to compare with tray_Test_1
 FlowSheet packedStage_BilletSchultes_Test_1
         PARAMETERS
@@ -146 +219 @@
         t1.Q = 0 * 'kW';
         t1.Ap = 0.8 * 'm^2';
-        t1.V = 0.8 * 'm^2' * 1 * 'm';
+        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 = 1 * 'm';
+        t1.hs = 0.4 * 'm';
         t1.Qsio = 1;
         
@@ -166 +239 @@
 end
 
+#the 3rd stage of the Distillation_kettle_cond_Test
 FlowSheet packedStage_BilletSchultes_Test_2
         PARAMETERS