Changeset 453 for branches/packed

Timestamp:

Feb 7, 2008, 5:26:42 PM (15 years ago)

Author:

Paula Bettio Staudt

Message:

More tests

Location:

branches/packed

Files:

• eml/stage_separators/column.mso (modified) (1 diff)
• eml/stage_separators/tray.mso (modified) (3 diffs)
• sample/stage_separators/sample_column.mso (modified) (1 diff)
• sample/stage_separators/sample_tray.mso (modified) (3 diffs)

branches/packed/eml/stage_separators/column.mso

@@ -1306 +1306 @@
         stage([top+topdown:topdown:bot]).OutletV to stage([top:topdown:bot-topdown]).InletV;
         stage([top:topdown:bot-topdown]).OutletL to stage([top+topdown:topdown:bot]).InletL;
-        
-        
-end
-
+end
+
+Model PackedDistillation_kettle_cond
+        PARAMETERS
+        outer PP as Plugin(Brief = "External Physical Properties", Type="PP");
+        outer NComp as Integer;
+        NStages as Integer(Brief="Number of stages", Default=2);
+        topdown as Integer(Brief="Trays counting (1=top-down, -1=bottom-up)", Default=1);
+        top as Integer(Brief="Number of top tray");
+        bot as Integer(Brief="Number of bottom tray");
+        
+        SET
+        top = (NStages-1)*(1-topdown)/2+1;
+        bot = NStages/top;
+        
+        VARIABLES
+        stage(NStages) as packedStage;
+        cond as condenser;
+        reb as reboiler;
+        sptop as splitter;
+        pump1 as pump;
+        
+        CONNECTIONS
+        #vapor
+        reb.OutletV to stage(bot).InletV;
+        stage([top+topdown:topdown:bot]).OutletV to stage([top:topdown:bot-topdown]).InletV;
+        stage(top).OutletV to cond.InletV;
+        
+        #liquid
+        cond.OutletL to sptop.Inlet;    
+        sptop.Outlet2 to pump1.Inlet;
+        pump1.Outlet to stage(top).InletL;
+        stage([top:topdown:bot-topdown]).OutletL to stage([top+topdown:topdown:bot]).InletL;
+        stage(bot).OutletL to reb.InletL;
+        
+        EQUATIONS
+        reb.OutletV.F = stage(bot).OutletV.F;
+end

branches/packed/eml/stage_separators/tray.mso

@@ -286 +286 @@
         d as length (Brief="size of an element of packing");
         #h as length (Brief="Height of packing");
-        C as Real (Brief="Constant for resitance factor equation");
-        Cp as Real (Brief="Constant for resitance at loading point factor equation");
+        C as Real (Brief="Constant for resitance factor equation"); #Billet pp 80 (Cf+Cfl)/2
+        Cp as Real (Brief="Constant for resitance at loading point factor equation"); # Billet pp94
         Mw(NComp)       as molweight    (Brief = "Component Mol Weight");
 
@@ -329 +329 @@
         uV * Across = OutletV.F * vV;
         
-        "Coefficient of Resistance"
-        ksi * C^2 * (uL/uV * sqrt(rhoV/rhoL) * (niL/niV)^5.8)^3 = g/1*'s^2/m';
+#       "Coefficient of Resistance"
+#       ksi * C^2 * (uL/uV * sqrt(rhoV/rhoL) * (niL/niV)^5.8)^3 = g/1*'s^2/m';
+        #ksil = Cp * (exp(uL*rhoL/a/niL/200)*(hL/hLs)^0.3) * (64/Rev+(1.8/Rev)) *
         ksil = Cp * (exp(uL*rhoL/a/niL/200)*(hL/hLs)^0.3) * (64/Rev+(1.8/Rev)) *
-                                ((e-hL)/e);#^1.5;
-        
-        "Liquid holdup and Liquid flow"
-        a^2 * niL * uL = hL^1 *(g*rhoL/3 - ksi*a*rhoV*uV^2/(4*hL*(e-hL)^2));
-
-        "Pressure drop and Vapor flow"
-        (InletV.P - OutletV.P)/hs = ksil *(a/2 + 2/ds)*(uV^2*rhoV/(e-hL)^3);
+                                        ((e-hL)/e);
+        
+#       "Liquid holdup and Liquid flow"
+#       a^2 * niL * uL = hL^1 *(g*rhoL/3 - ksi*a*rhoV*uV^2/(4*hL*(e-hL)^2));
+
+#       "Pressure drop and Vapor flow"
+#       (InletV.P - OutletV.P)/hs = ksil *(a/2 + 2/ds)*(uV^2*rhoV/(e-hL)^3);
         
         "Vapor Reynolds number"
@@ -344 +345 @@
 
         "Theoretical Liquid Holdup"
-        hLs = (12*a^2*niL*uL/g/rhoL)^0.333;
+        #hLs = (12*a^2*niL*uL/g/rhoL)^0.333;
+        hLs = (12*a^2*niL*uL/g/rhoL);
 end
-
 
 #*-------------------------------------------------------------------

branches/packed/sample/stage_separators/sample_column.mso

@@ -582 +582 @@
 end
 
+FlowSheet Packed_kettle_cond_Test
+        PARAMETERS
+        PP      as Plugin(Brief="Physical Properties",
+                Type="PP",
+                Components = [ "methanol", "water"],
+                LiquidModel = "IdealLiquid",
+                VapourModel = "Ideal"
+        );
+        NComp   as Integer;
+        
+        VARIABLES
+        Qc as energy_source (Brief="Heat rate removed from condenser");
+        Qr as energy_source (Brief="Heat rate supplied to reboiler");
+        
+        SET
+        NComp = PP.NumberOfComponents;
+
+        DEVICES
+        col as PackedDistillation_kettle_cond;
+        feed as source;
+        zero as stream;
+        
+        CONNECTIONS
+        feed.Outlet to col.reb.Inlet;
+        zero to col.stage([1:col.NStages]).Inlet;
+        Qc.OutletQ to col.cond.InletQ;
+        Qr.OutletQ to col.reb.InletQ;
+        
+        SPECIFY
+        feed.Outlet.F = 0.98 * 'mol/min';
+        feed.Outlet.T = (50+273.15) * 'K';
+        feed.Outlet.P = 0.91 * 'kPa';
+        feed.Outlet.z = [0.16, 0.84];
+        
+        zero.F = 0 * 'kmol/h';
+        zero.T = 300 * 'K';
+        zero.P = 1 * 'atm';
+        zero.z = 1/NComp;
+        zero.v = 0;
+        zero.h = 0 * 'J/mol';
+        
+        col.sptop.frac = 0.4;
+        col.sptop.Outlet1.F = 0.10 * 'mol/min';
+        col.reb.OutletL.F = 0.88 * 'mol/min';
+        col.cond.OutletV.F = 0 * 'kmol/h';
+        Qr.OutletQ.Q = 700 * 'cal/min';
+        Qc.OutletQ.Q = -700 * 'cal/min';
+        col.pump1.dP = 0.1 * 'atm';
+        col.stage.Emv = 1;
+        col.stage.hs = 100 * 'mm';
+        
+        SET
+        col.NStages = 10;
+        col.cond.V = 1 * 'l';
+        col.cond.Across = 100 * 'cm^2';
+        col.reb.V = 2 * 'l';
+        col.reb.Across = 200 * 'cm^2';
+        col.stage.V = 576.4 * 'cm^3';
+        col.stage.Q = 0 * 'kW';
+        col.stage.Ap = 56.74 * 'cm^2';
+        col.stage.ds = 85 * 'mm';
+        col.stage.d = 20 * 'mm';
+        col.stage.C = 2.18;
+        col.stage.Cp = 1.329;
+        col.stage.e = 0.662;
+        col.stage.a = 185.4 * 'm^2/m^3';
+
+        INITIAL
+        # condenser
+        col.cond.OutletL.T = (63.5+273.15) *'K';
+        col.cond.Level = 2 * 'cm';
+        col.cond.OutletL.z(1) = 0.16;
+
+        # reboiler
+        col.reb.OutletL.T = (82+273.15) *'K';
+        col.reb.Level = 2 * 'cm';
+        col.reb.OutletL.z(1) = 0.16;
+
+        # 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 = 1 * 'cm';
+        col.stage.OutletL.z(1) = 0.16;
+
+        OPTIONS
+        TimeStep = 0.1;
+        TimeEnd = 50;
+        #GuessFile = "/home/paula/SectionColumn_Test_with8tray.rlt";
+        #GuessFile="Distillation_kettle_cond_Test.rlt";
+        #Dynamic = false;       
+end

branches/packed/sample/stage_separators/sample_tray.mso

@@ -143 +143 @@
 
         t1.Emv = 1;
-        t1.hs = 0.075 * 'm';
-#       t1.OutletV.F = 147.1 * 'kmol/h';
-#       t1.OutletL.F = 229.5 * 'kmol/h';
-#       t1.ksi = 0.805;
+        t1.hs = 0.4 * 'm';#0.075 * 'm';
+        t1.OutletV.F = 147.1 * 'kmol/h';
+        t1.OutletL.F = 229.5 * 'kmol/h';
+        t1.ksi = 0.805;
         
         SET
@@ -152 +152 @@
         t1.Q = 0 * 'kW';
         t1.Across = 0.8 * 'm^2';
-        t1.V = 0.8 * 'm^2' * 0.075 * 'm';
+        t1.V = 0.8 * 'm^2' * 0.4 * 'm';
         t1.ds = 1.009 * 'm';
         t1.d = 50 * 'mm';
@@ -167 +167 @@
         OPTIONS
         #InitialFile = "/home/paula/tray_Test.rlt";
+        InitialFile = "/home/paula/packedStage_Test.rlt";
         #DAESolver(File="dassl");
-        TimeStep = 0.1; 
-        TimeEnd = 10;
+        TimeStep = 2;   
+        TimeEnd = 50;
 end