Changeset 467

Timestamp:

Feb 28, 2008, 3:49:14 PM (16 years ago)

Author:

Paula Bettio Staudt

Message:

More shots

Location:

branches/packed

Files:

• eml/stage_separators/tray.mso (modified) (5 diffs)
• sample/stage_separators/packedStage_BilletSchultes_Test_1.rlt (added)
• sample/stage_separators/sample_column.mso (modified) (7 diffs)
• sample/stage_separators/sample_tray.mso (modified) (2 diffs)

branches/packed/eml/stage_separators/tray.mso

@@ -358 +358 @@
         Mw(NComp)       as molweight    (Brief = "Component Mol Weight");
         hs as length (Brief="Height of the packing stage");
+        Qsio as Real (Brief="Resistance coefficient", Lower = 0);
         
         VARIABLES
@@ -369 +370 @@
         invK as Real (Brief="Wall factor");
         Rev as Real (Brief="Reynolds number of the vapor stream", Lower = 0, Default=100);
-        Qsio as Real (Brief="Resistance coefficient", Lower = 0);
+        #Qsio as Real (Brief="Resistance coefficient", Lower = 0);
         
         SET
@@ -385 +386 @@
 
         "Volume flow rate of liquid, m^3/m^2/s"
-        uL * Ap = OutletL.F * vL;
+        uL * Ap * e = OutletL.F * vL;
         "Volume flow rate of vapor, m^3/m^2/s"
-        uV * Ap = OutletV.F * vV;
+        uV * Ap * e = OutletV.F * vV;
         
         "Liquid holdup and Liquid flow"
@@ -401 +402 @@
         Rev*invK = dp*uV*rhoV / (miV*(1-e));
 
-        if Rev > 1e-4 then
+#*      if Rev > 1e-4 then
                 "Resistance Coefficient"        
                 Qsio = Cpo * (64/Rev + 1.8/Rev^0.08);
@@ -407 +408 @@
                 Qsio = 1;
         end
-
+*#
         "Pressure drop and Vapor flow"
-        (InletV.P - OutletV.P)*0.0091/hs  = Qsio*a*uV^2*rhoV*invK / (2*e^3);
+        (InletV.P - OutletV.P)/hs  = Qsio*a*uV^2*rhoV*invK / (2*e^3);
         
         "Efficiency"

branches/packed/sample/stage_separators/sample_column.mso

@@ -531 +531 @@
         reb to sec.stage(8).InletV;
         cond to sec.stage(1).InletL;
-
         
         SPECIFY
@@ -556 +555 @@
         reb.z = [0.16, 0.542, 0.013, 0.008, 0.277];
         
-        #sec.stage.Qsio = 0.8;
-        
-        SET
-        sec.H = 6 * 'ft';
+        SET
+        sec.H = 16 * 'm';
         sec.NStages = 8;
         sec.stage.Q = 0 * 'kW';
@@ -567 +564 @@
         sec.stage.e = 0.951;
         sec.stage.a = 100*112.6 * 'm^2/m^3';
+        sec.stage.Qsio = 1;
 
         INITIAL
         sec.stage.OutletL.T =[283:(325-283)/(sec.NStages-1):325] *'K';
-        sec.stage.ML = 0.5 * 'kmol';
+        sec.stage.ML = 0.3 * 'kmol';
         sec.stage.OutletL.z([1:4]) = [0.2, 0.2, 0.2, 0.2];
 
@@ -608 +606 @@
         Qr.OutletQ to col.reb.InletQ;
         
+        VARIABLES
+        deltaP(col.NStages) as Real (Unit = 'atm/m');
+        
         SPECIFY
         feed.Outlet.F = 0.98 * 'mol/min';
@@ -628 +629 @@
         Qc.OutletQ.Q = -700 * 'cal/min';
         col.pump1.dP = 0.1 * 'atm';
-        col.stage.Qsio = 0.05;
-        
-#       EQUATIONS
-#       col.reb.OutletV.F = col.K * sqrt(Qr.OutletQ.Q);
-        
-        SET
-        col.H = 1 * 'm';
-        col.NStages = 5;
+        
+        EQUATIONS
+        col.reb.OutletV.F = 0.2 * 'mol*min^0.5/kg^0.5/m' * sqrt(Qr.OutletQ.Q);
+        deltaP = (col.stage.InletV.P - col.stage.OutletV.P)/col.stage.hs;
+        
+        SET
+        col.H = 2 * 'm';
+        col.NStages = 2;
 
         col.cond.V = 1 * 'l';
@@ -650 +651 @@
         col.stage.a = 185.4 * 'm^2/m^3';
         col.stage.Cpo = 0.763;
+        col.stage.Qsio = 1;
         
         INITIAL
@@ -660 +662 @@
         col.reb.OutletL.T = (82+273.15) *'K';
         col.reb.Level = 2 * 'cm';
-        col.reb.OutletL.z(1) = 0.16;
+        col.reb.OutletL.z(1) = 0.7;
 
         # column stages
         col.stage.OutletL.T = [(63.5+273.15):((82+273.15)-(63.5+273.15))/(col.NStages-1):(82+273.15)] * 'K';
-        col.stage.Level = 2 * 'cm';
-        #col.stage.ML = 0.2 * 'mol';
-        col.stage.OutletL.z(1) = 0.5; #0.16;
+        #col.stage.Level = 2 * 'cm';
+        col.stage.ML = 0.2 * 'mol';
+        col.stage.OutletL.z(1) = 0.5; #[0.79, 0.65, 0.5, 0.25, 0.16];
 
         OPTIONS

branches/packed/sample/stage_separators/sample_tray.mso

@@ -96 +96 @@
 end
 
-FlowSheet packedStage_BilletSchultes_Test
+FlowSheet packedStage_BilletSchultes_Test_1
+        PARAMETERS
+        PP      as Plugin(Brief="Physical Properties",
+                Type="PP",
+                Components = [ "n-pentane", "benzene"],
+                LiquidModel = "PR",
+                VapourModel = "PR"
+        );
+        NComp   as Integer;
+
+        SET
+        NComp = PP.NumberOfComponents;
+        
+        VARIABLES
+        deltaP as Real (Unit='atm/m');
+        
+        DEVICES
+        t1 as packedStage_BilletSchultes;
+        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)/t1.hs;
+        
+        SPECIFY
+        feed.Outlet.F = 113.4 * 'kmol/h';
+        feed.Outlet.T = 291 * 'K';
+        feed.Outlet.P = 1.66 * 'atm';
+        feed.Outlet.z = [0.5, 0.5];
+        
+        inL.P = 165 * 'kPa';
+        inL.T = 310 * 'K';
+        inL.F = 61.99 * 'kmol/h';
+        inL.z = [0.1641, 0.8359];
+
+        inV.F = 201.25 * 'kmol/h';
+        inV.P = 150 * 'kPa';
+        inV.T = 321 * 'K';
+        inV.z = [0.0584, 0.9416];
+
+        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' * 1 * '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.Qsio = 1;
+        
+        INITIAL
+        t1.OutletL.T = 320 *'K';
+        t1.ML = 0.25 * 'kmol';
+        t1.OutletL.z(1) = 0.5;
+        
+        OPTIONS
+        InitialFile = "packedStage_BilletSchultes_Test_1.rlt";
+        DAESolver(File="sundials");
+        TimeStep = 100;
+        TimeEnd = 1000;
+end
+
+FlowSheet packedStage_BilletSchultes_Test_2
         PARAMETERS
         PP      as Plugin(Brief="Physical Properties",
@@ -155 +225 @@
         t1.a = 112.6 * 'm^2/m^3';
         t1.hs = 0.4 * 'm';
+        t1.Qsio = 1;
 
         INITIAL