Changeset 498

Timestamp:

Apr 11, 2008, 3:21:24 PM (15 years ago)

Author:

Paula Bettio Staudt

Message:

First official version of packed stage and packed column with samples

Files:

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

branches/packed/eml/stage_separators/column.mso

@@ -1263 +1263 @@
 
 
-Model Section_Column_Packed
+Model Packed_Section_Column
         
         PARAMETERS
@@ -1275 +1275 @@
 
         VARIABLES
-        stage(NStages) as packedStage;#_BilletSchultes;
+        stage(NStages) as packedStage;
         
         SET
@@ -1288 +1288 @@
 end
 
-Model PackedDistillation_kettle_cond as Section_Column_Packed
+Model PackedDistillation_kettle_cond as Packed_Section_Column
         PARAMETERS
         VapourFlow as Switcher(Valid = ["on", "off"], Default = "on");

branches/packed/eml/stage_separators/tray.mso

@@ -348 +348 @@
         Mw(NComp)       as molweight    (Brief = "Component Mol Weight");
         hs as length (Brief="Height of the packing stage");
+        Qsil as positive (Brief="Resistance coefficient on the liquid load", Default=1);
 
         VARIABLES
@@ -377 +378 @@
         Al as area (Brief="Area occupied by the liquid", Default=0.001, Upper=1);
         hl as positive (Brief="Column holdup", Unit='m^3/m^3', Default=0.01,Upper=10);
-        Qsil as positive (Brief="Resistance coefficient on the liquid load");
 
         SET

branches/packed/sample/stage_separators/sample_column.mso

@@ -519 +519 @@
         
         DEVICES
-        sec as Section_Column_Packed;
+        sec as Packed_Section_Column;
         feed as liquid_stream;
         reb as vapour_stream;
@@ -555 +555 @@
         reb.z = [0.16, 0.542, 0.013, 0.008, 0.277];
 
-        sec.stage.Qsil = 0.1;
-        #sec.stage(1).OutletV.F = 150 * 'kmol/h';
         sec.stage.deltaP = 0.001 * 'atm';
 
@@ -567 +565 @@
         sec.stage.e = 0.951;
         sec.stage.a = 112.6 * 'm^2/m^3';
+        sec.stage.Qsil = 0.1;
 
         INITIAL

branches/packed/sample/stage_separators/sample_tray.mso

@@ -371 +371 @@
         t1.OutletV.P = 145 * 'kPa';
         #t1.OutletV.F = 190 * 'kmol/h';
-        t1.Qsil = 10;
         
         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' * 0.4 * 'm';
-        t1.ds = 1.009 * 'm';
+        t1.d = 1.009 * 'm';
         t1.Cpo = 0.763;
         t1.e = 0.951;
         t1.a = 112.6 * 'm^2/m^3';
         t1.hs = 0.4 * 'm';
+        t1.Qsil = 10;
         
         INITIAL
@@ -447 +446 @@
         t1.deltaP = 0.01 * 'atm';
         #t1.OutletV.F = 165 * 'kmol/h';
-        t1.Qsil = 10;
         
         SET
@@ -461 +459 @@
 
         t1.V = 4 * 'ft^2' * 1 * 'ft';
-        t1.Ap = 4 * 'ft^2';
         t1.hs = 1 * 'ft';
-        t1.ds = 2.26 * 'ft';
+        t1.d = 2.26 * 'ft';
+        t1.Qsil = 10;
 
         INITIAL

trunk/eml/stage_separators/column.mso

@@ -1261 +1261 @@
         
 end
+
+#*----------------------------------------------------------------------
+* Model of a  packed column section with:
+*       - NStages = theoretical number of equilibrium stages.
+* 
+*---------------------------------------------------------------------*# 
+Model Packed_Section_Column
+        ATTRIBUTES
+        Pallete         = true;
+        Icon            = "icon/SectionColumn"; 
+        Brief           = "Model of a packed column section.";
+        Info            =
+"== Model of a packed column section containing ==
+* NStages theoretical stages.
+        
+== Specify ==
+* the feed stream of each tray (Inlet);
+* the InletL stream of the top tray;
+* the InletV stream of the bottom tray;
+* the pressure drop (deltaP) of each tray.
+        
+== Initial Conditions ==
+* the stages temperature (OutletL.T);
+* the stages liquid holdup;
+* (NoComps - 1) OutletL (OR OutletV) compositions for each tray.
+";
+        
+        PARAMETERS
+        outer PP as Plugin(Brief = "External Physical Properties", Type="PP");
+        outer NComp as Integer;
+        NStages as Integer(Brief="Number of trays", 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");
+        H as length (Brief="Height of packing");
+
+        VARIABLES
+        stage(NStages) as packedStage;
+        
+        SET
+        top = (NStages-1)*(1-topdown)/2+1;
+        bot = NStages/top;
+        stage.hs = H/NStages;   
+        stage.V = stage.hs * stage.d^2*3.14159/4;
+        
+        CONNECTIONS
+        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
+
+#*----------------------------------------------------------------------
+* Model of a  packed distillation column containing:
+*       - a section with NStages theoretical stages;
+*       - a kettle reboiler;
+*       - dymamic condenser;
+*       - a splitter which separate reflux and distillate;
+*       - a pump in reflux stream;
+*---------------------------------------------------------------------*# 
+Model PackedDistillation_kettle_cond as Packed_Section_Column
+        PARAMETERS
+        VapourFlow as Switcher(Valid = ["on", "off"], Default = "on");
+        
+        VARIABLES
+        cond as condenser;
+        reb as reboiler;
+        sptop as splitter;
+        pump1 as pump;
+        
+        CONNECTIONS
+        #vapor
+        reb.OutletV to stage(bot).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(bot).OutletL to reb.InletL;
+        
+        EQUATIONS
+        switch VapourFlow
+                case "on":
+                stage(bot).InletV.F*stage(bot).vV = sqrt((reb.OutletV.P - stage(bot).OutletV.P)/
+                                        (stage(bot).rhoV*stage(bot).Qsil*20))*(stage(bot).d^2*3.14159/4*stage(bot).e - stage(bot).Al);
+                when stage(bot).InletV.F < 1e-6 * 'kmol/h' switchto "off";
+                
+                case "off":
+                stage(bot).InletV.F = 0 * 'mol/s';
+                when reb.OutletV.P > stage(bot).OutletV.P + 1e-1 * 'atm' switchto "on";
+        end
+
+end

trunk/eml/stage_separators/tray.mso

@@ -323 +323 @@
         V = ML* vL + MV*vV;
 end
+
+#*-------------------------------------
+* Model of a packed column stage
+-------------------------------------*#
+Model packedStage
+        ATTRIBUTES
+        Pallete         = false;
+        Icon            = "icon/Tray"; 
+        Brief           = "Complete model of a packed column stage.";
+        Info            =
+"== Specify ==
+* the Feed stream
+* the Liquid inlet stream
+* the Vapour inlet stream
+* the stage pressure drop (deltaP)
+        
+== Initial ==
+* the plate temperature (OutletL.T)
+* the liquid molar holdup ML
+* (NoComps - 1) OutletL compositions
+";      
+        
+        PARAMETERS
+outer PP as Plugin(Brief = "External Physical Properties", Type="PP");
+outer NComp as Integer;
+        PPwater as Plugin(Brief="Physical Properties",
+                Type="PP",
+                Components = [ "water" ],
+                LiquidModel = "PR",
+                VapourModel = "PR"
+        );
+
+        V as volume(Brief="Total Volume of the tray");
+        Q as heat_rate (Brief="Rate of heat supply"); 
+        d as length (Brief="Column diameter");  
+
+        a as Real (Brief="surface area per packing volume", Unit='m^2/m^3');
+        g as acceleration;
+        e as Real (Brief="Void fraction of packing, m^3/m^3");
+        Cpo as Real (Brief="Constant for resitance equation"); # Billet and Schultes, 1999.
+        Mw(NComp)       as molweight    (Brief = "Component Mol Weight");
+        hs as length (Brief="Height of the packing stage");
+        Qsil as positive (Brief="Resistance coefficient on the liquid load", Default=1);
+
+        VARIABLES
+in      Inlet as stream (Brief="Feed stream", PosX=0, PosY=0.4932, Symbol="_{in}");
+in      InletL as stream (Brief="Inlet liquid stream", PosX=0.5195, PosY=0, Symbol="_{inL}");
+in      InletV as stream (Brief="Inlet vapour stream", PosX=0.4994, PosY=1, Symbol="_{inV}");
+out     OutletL as liquid_stream (Brief="Outlet liquid stream", PosX=0.8277, PosY=1, Symbol="_{outL}");
+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", 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");
+        vV as volume_mol (Brief="Vapour Molar volume");
+        
+        miL as viscosity (Brief="Liquid dynamic viscosity", DisplayUnit='kg/m/s');
+        miV as viscosity (Brief="Vapor dynamic viscosity", DisplayUnit='kg/m/s');
+        rhoL as dens_mass;
+        rhoV as dens_mass;
+        
+        deltaP as pressure;
+        
+        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", Default=1, Upper=10);
+        Rev as Real (Brief="Reynolds number of the vapor stream", Default=4000);
+        Al as area (Brief="Area occupied by the liquid", Default=0.001, Upper=1);
+        hl as positive (Brief="Column holdup", Unit='m^3/m^3', Default=0.01,Upper=10);
+
+        SET
+        Mw = PP.MolecularWeight();
+
+        EQUATIONS
+        "Component Molar Balance"
+        diff(M)=Inlet.F*Inlet.z + InletL.F*InletL.z + InletV.F*InletV.z
+                - OutletL.F*OutletL.z - OutletV.F*OutletV.z;
+
+        "Energy Balance"
+        diff(E) = ( Inlet.F*Inlet.h + InletL.F*InletL.h + InletV.F*InletV.h
+                - OutletL.F*OutletL.h - OutletV.F*OutletV.h + Q );
+        
+        "Molar Holdup"
+        M = ML*OutletL.z + MV*OutletV.z;
+        
+        "Energy Holdup"
+        E = ML*OutletL.h + MV*OutletV.h - OutletL.P*V;
+        
+        "Mol fraction normalisation"
+        sum(OutletL.z)= 1.0;
+        
+        "Liquid Volume"
+        vL = PP.LiquidVolume(OutletL.T, OutletL.P, OutletL.z);
+        "Vapour Volume"
+        vV = PP.VapourVolume(OutletV.T, OutletV.P, OutletV.z);
+        
+        "Chemical Equilibrium"
+        PP.LiquidFugacityCoefficient(OutletL.T, OutletL.P, OutletL.z)*OutletL.z = 
+                PP.VapourFugacityCoefficient(OutletV.T, OutletV.P, OutletV.z)*OutletV.z;
+        
+        "Thermal Equilibrium"
+        OutletV.T = OutletL.T;
+        
+        "Mechanical Equilibrium"
+        OutletL.P = OutletV.P;
+        
+        "Geometry Constraint"
+        V*e = ML*vL + MV*vV;
+        
+        "Liquid Density"
+        rhoL = PP.LiquidDensity(OutletL.T, OutletL.P, OutletL.z);
+        "Vapour Density"
+        rhoV = PP.VapourDensity(InletV.T, InletV.P, InletV.z);
+        "Liquid viscosity"
+        miL = PP.LiquidViscosity(OutletL.T, OutletL.P, OutletL.z);
+        "Vapour viscosity"
+        miV = PP.VapourViscosity(InletV.T, InletV.P, InletV.z);
+
+        "Area occupied by the liquid"
+        Al = ML*vL/hs;
+
+        "Volume flow rate of liquid, m^3/m^2/s"
+        uL * Al = OutletL.F * vL;
+        "Volume flow rate of vapor, m^3/m^2/s"
+        uV * ((d^2*3.14159/4)*e - Al) = OutletV.F * vV;
+        
+        "Liquid holdup"
+        hl = ML*vL/V/e;
+        
+        "Particle diameter"
+        dp = 6 * (1-e)/a;
+        
+        "Wall Factor"
+        invK = (1 + (2*dp/(3*d*(1-e))));
+        
+        "Reynolds number of the vapor stream"
+        Rev*invK = dp*uV*rhoV / (miV*(1-e));
+        
+        deltaP = InletV.P - OutletV.P;
+        
+        "Pressure drop and Vapor flow"
+        deltaP/hs  = Qsil*a*uV^2*rhoV*invK / (2*(e-hl)^3);
+
+        "Liquid holdup"
+        hl = (12*miL*a^2*uL/rhoL/g)^1/3;
+end

trunk/sample/stage_separators/sample_column.mso

@@ -503 +503 @@
 end
 
-
+# packed column section with 8 trays
+FlowSheet PackedSectionColumn
+        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;
+        
+        DEVICES
+        sec as Packed_Section_Column;
+        feed as liquid_stream;
+        reb as vapour_stream;
+        cond as liquid_stream;
+        zero as stream;
+        
+        CONNECTIONS
+        feed to sec.stage(5).Inlet;
+        zero to sec.stage([1:4]).Inlet;
+        zero to sec.stage([6:sec.NStages]).Inlet;
+        reb to sec.stage(8).InletV;
+        cond to sec.stage(1).InletL;
+        
+        SPECIFY
+        feed.F = 113.4 * 'kmol/h';
+        feed.T = 291 * 'K';
+        feed.P = 1.66 * 'atm';
+        feed.z = [0.2, 0.2, 0.2, 0.2, 0.2];
+
+        zero.F = 0 * 'kmol/h';
+        zero.T = 300 * 'K';
+        zero.P = 1 * 'atm';
+        zero.z = [0.2, 0.2, 0.2, 0.2, 0.2];
+        zero.v = 0;
+        zero.h = 0 * 'J/mol';
+        
+        cond.F = 85 * 'kmol/h';
+        cond.P = 2.33 * 'atm';
+        cond.T = 283.5 * 'K';
+        cond.z = [0.599, 0.044, 0.035, 0.007, 0.315];
+
+        reb.F = 137.57 * 'kmol/h';
+        reb.P = 2.46 * 'atm';
+        reb.T = 325 * 'K';
+        reb.z = [0.16, 0.542, 0.013, 0.008, 0.277];
+
+        sec.stage.deltaP = 0.001 * 'atm';
+
+        SET
+        sec.H = 4 * 'm';
+        sec.NStages = 8;
+        sec.stage.Q = 0 * 'kW';
+        sec.stage.d = 1.009 * 'm';
+        sec.stage.Cpo = 0.763;
+        sec.stage.e = 0.951;
+        sec.stage.a = 112.6 * 'm^2/m^3';
+        sec.stage.Qsil = 0.1;
+
+        INITIAL
+        sec.stage.OutletL.T =[283:(325-283)/(sec.NStages-1):325] *'K';
+        sec.stage.ML = 0.01 * 'kmol';
+        sec.stage.OutletL.z([1:4]) = [0.3, 0.2, 0.002, 0.1];
+
+        OPTIONS
+        DAESolver(File="dassl");
+        NLASolver(File="nlasolver");
+        TimeStep = 10;
+        TimeEnd = 100;
+end
+
+FlowSheet Packed_kettle_cond_Test
+        PARAMETERS
+        PP      as Plugin(Brief="Physical Properties",
+                Type="PP",
+                Components = [ "isobutane", "n-pentane", "propylene",
+                "benzene", "isobutene" ],
+                LiquidModel = "PR",
+                VapourModel = "PR"
+        );
+        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.stage(5).Inlet;
+        zero to col.reb.Inlet;
+        zero to col.stage([1:4]).Inlet;
+        zero to col.stage([6:col.NStages]).Inlet;
+        Qc.OutletQ to col.cond.InletQ;
+        Qr.OutletQ to col.reb.InletQ;
+        
+        SPECIFY
+        feed.Outlet.F = 113.4 * 'kmol/h';
+        feed.Outlet.T = 291 * 'K';
+        feed.Outlet.P = 168.3 * 'kPa';
+        feed.Outlet.z = 1/NComp;
+        
+        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.Outlet2.F = 85 * 'kmol/h';
+        col.reb.OutletL.F = 28.4 * 'kmol/h';
+        col.sptop.frac = 0.444445;
+        col.cond.OutletV.F = 0 * 'kmol/h';
+        Qr.OutletQ.Q = 3.7743e6 * 'kJ/h';
+        Qc.OutletQ.Q = -3.71e6 * 'kJ/h';
+        col.pump1.dP = 16 * 'kPa';
+
+        col.stage.deltaP = 0.003 * 'atm';
+
+        SET
+        col.H = 3.5 * 'm';
+        col.NStages = 8;
+        col.stage.Q = 0 * 'kW';
+        col.stage.d = 2.24 * 'ft';
+        col.stage.Cpo = 0.763;
+        col.stage.e = 0.951;
+        col.stage.a = 112.6 * 'm^2/m^3';
+
+        col.cond.V = 2 * 'm^3';
+        col.cond.Across = 1 * 'm^2';
+        col.reb.V = 2 * 'm^3';
+        col.reb.Across = 1 * 'm^2';
+        col.stage.Qsil = 2;
+        
+        INITIAL
+        # condenser
+        col.cond.OutletL.T = 260 *'K';
+        col.cond.Level = 1 * 'm';
+        col.cond.OutletL.z([1:4]) = [0.65, 0.05, 0.01, 0.01];
+
+        # reboiler
+        col.reb.OutletL.T = 330 *'K';
+        col.reb.Level = 1 * 'm';
+        col.reb.OutletL.z([1:4]) = [0.1, 0.7, 0.01, 0.01];
+
+        # column trays
+        col.stage.OutletL.T = [290:(330-290)/(col.NStages-1):330] * 'K';
+        col.stage.ML = 0.1 * 'kmol';
+        col.stage.OutletL.z([1:4]) = [0.15, 0.5, 0.001, 0.1];
+
+        OPTIONS
+        TimeStep = 10;
+        TimeEnd = 500;
+end
+
+FlowSheet PackedColumn_ctrl
+        PARAMETERS
+        PP      as Plugin(Brief="Physical Properties",
+                Type="PP",
+                Components = [ "isobutane", "n-pentane", "propylene",
+                        "benzene", "isobutene" ],
+                LiquidModel = "PR",
+                VapourModel = "PR"
+        );
+        NComp   as Integer;
+        
+        Qcmin as heat_rate (Brief="Minimum Condenser Heat supplied");
+        Qcmax as heat_rate (Brief="Maximum Condenser Heat supplied");
+        Qrmin as heat_rate (Brief="Minimum Reboiler Heat supplied");
+        Qrmax as heat_rate (Brief="Maximum Reboiler Heat supplied");
+        Frmin as flow_mol (Brief="Minimum bottom flow rate");
+        Frmax as flow_mol (Brief="Maximum bottom flow rate");
+        Fcmin as flow_mol (Brief="Minimum reflux flow rate");
+        Fcmax as flow_mol (Brief="Maximum reflux flow rate");
+        Hmint as length (Brief="Minimum liquid level in top tank");
+        Hmaxt as length (Brief="Maximum liquid level in top tank");
+    Hminb as length (Brief="Minimum liquid level in reboiler");
+        Hmaxb as length (Brief="Maximum liquid level in reboiler");
+        Pmax as pressure (Brief="Maximum column pressure");
+        Pmin as pressure (Brief="Minimum column pressure");
+        Tmax as temperature (Brief="Maximum column temperature");
+        Tmin as temperature (Brief="Minimum column temperature");
+
+        VARIABLES
+        Qc as energy_source (Brief="Heat rate removed from condenser");
+        Qr as energy_source (Brief="Heat rate supplied to reboiler");
+        Had_top as Real (Brief="Dimensionless condenser level");
+        Had_bot as Real (Brief="Dimensionless reboiler level");
+        Pad as Real (Brief="Dimensionless pressure");
+        Tad as Real (Brief="Dimensionless temperature");
+        RR      as positive (Brief="Reflux ratio");
+
+        SET
+        NComp = PP.NumberOfComponents;
+
+        DEVICES
+        col as PackedDistillation_kettle_cond;
+        feed as source;
+        zero as stream;
+        TCcond as PIDIncr;
+        LCtop as PIDIncr;
+        LCbot as PIDIncr;
+        PC as PIDIncr;
+
+        CONNECTIONS
+        feed.Outlet to col.stage(5).Inlet;
+        zero to col.reb.Inlet;
+        zero to col.stage([1:4]).Inlet;
+        zero to col.stage([6:col.NStages]).Inlet;
+        Qc.OutletQ to col.cond.InletQ;
+        Qr.OutletQ to col.reb.InletQ;
+
+        EQUATIONS
+   "Temperature Controller"
+        TCcond.Parameters.tau = 0*'s';  
+        TCcond.Parameters.tauSet = 0*'s';       
+        TCcond.Parameters.alpha = 0.3;
+        TCcond.Parameters.bias = 0.5;   
+        TCcond.Parameters.gamma = 1;
+        TCcond.Parameters.beta = 1;
+        TCcond.Options.action = 1;
+        TCcond.Options.clip = 1;
+        TCcond.Options.autoMan = 0;
+        TCcond.Parameters.intTime = 60*'s';
+        TCcond.Parameters.gain = 0.6;
+        TCcond.Parameters.derivTime = 1*'s';
+        TCcond.Ports.setPoint = ((15+273.15) * 'K' - Tmin)/(Tmax-Tmin);
+        TCcond.Ports.input = Tad;
+        Tad = (col.cond.OutletL.T-Tmin)/(Tmax-Tmin);
+        Qc.OutletQ.Q = Qcmin+(Qcmax-Qcmin)*TCcond.Ports.output; 
+
+        "Pressure Controller"
+        PC.Parameters.tau = 0*'s';      
+        PC.Parameters.tauSet = 0*'s';   
+        PC.Parameters.alpha = 0.3;
+        PC.Parameters.bias = 0; 
+        PC.Parameters.gamma = 1;
+        PC.Parameters.beta = 1;
+        PC.Options.action = -1;
+        PC.Options.clip = 1;
+        PC.Options.autoMan = 0;
+        PC.Parameters.intTime = 50*'s';
+        PC.Parameters.gain = 0.5;
+        PC.Parameters.derivTime = 1*'s';
+        PC.Ports.setPoint = (4.0*'bar'-Pmin)/(Pmax-Pmin);
+        PC.Ports.input = Pad;
+        Pad = (col.cond.OutletV.P-Pmin)/(Pmax-Pmin);
+        col.cond.OutletV.F = (Fcmin+(Fcmax-Fcmin)*PC.Ports.output);     
+        
+        "Ttop Level Controller"
+        LCtop.Parameters.tau = 0*'s';   
+        LCtop.Parameters.tauSet = 0*'s';        
+        LCtop.Parameters.alpha = 0.3;
+        LCtop.Parameters.bias = 0.5;    
+        LCtop.Parameters.gamma = 1;
+        LCtop.Parameters.beta = 1;
+        LCtop.Options.action = -1;
+        LCtop.Options.clip = 1;
+        LCtop.Options.autoMan = 0;
+        LCtop.Parameters.intTime = 10*'s';
+        LCtop.Parameters.gain = 1;
+        LCtop.Parameters.derivTime = 0*'s';
+        LCtop.Ports.setPoint = (1.0 * 'm' - Hmint)/(Hmaxt-Hmint);
+        LCtop.Ports.input = Had_top;
+        Had_top = (col.cond.Level-Hmint)/(Hmaxt-Hmint);
+        col.sptop.Outlet1.F = Fcmin + (Fcmax-Fcmin) * LCtop.Ports.output;
+
+        "Tbottom Level Controller"
+        LCbot.Parameters.tau = 0*'s';   
+        LCbot.Parameters.tauSet = 0*'s';        
+        LCbot.Parameters.alpha = 0.3;
+        LCbot.Parameters.bias = 0.5;    
+        LCbot.Parameters.gamma = 1;
+        LCbot.Parameters.beta = 1;
+        LCbot.Options.action = -1;
+        LCbot.Options.clip = 1;
+        LCbot.Options.autoMan = 0;
+        LCbot.Parameters.intTime = 100*'s';
+        LCbot.Parameters.gain = 1;
+        LCbot.Parameters.derivTime = 0*'s';
+        LCbot.Ports.setPoint = (1.0 * 'm' - Hminb)/(Hmaxb-Hminb);
+        LCbot.Ports.input = Had_bot;
+        Had_bot = (col.reb.Level-Hminb)/(Hmaxb-Hminb);
+        col.reb.OutletL.F = Frmin + (Frmax-Frmin) * LCbot.Ports.output;
+
+        RR * (col.cond.OutletV.F + col.sptop.Outlet1.F) = col.sptop.Outlet2.F;
+        
+        if time < 1 * 'h' then
+                col.sptop.Outlet2.F = 75 * 'kmol/h'; # reflux
+        else
+                col.sptop.Outlet2.F = 85 * 'kmol/h'; # reflux
+        end
+
+        SPECIFY
+        feed.Outlet.F = 113.4 * 'kmol/h';
+        feed.Outlet.T = 291 * 'K';
+        feed.Outlet.P = 5 * 'bar';
+        feed.Outlet.z = 1/NComp;
+        
+        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';
+        
+        Qr.OutletQ.Q = 4e6 * 'kJ/h';
+        col.pump1.dP = 16 * 'kPa';
+        
+        col.stage.deltaP = 0.003 * 'atm';
+
+        SET
+        col.H = 3.5 * 'm';
+        col.NStages = 8;
+        col.stage.Q = 0 * 'kW';
+        col.stage.d = 2.24 * 'ft';
+        col.stage.Cpo = 0.763;
+        col.stage.e = 0.951;
+        col.stage.a = 112.6 * 'm^2/m^3';
+
+        col.cond.V = 2 * 'm^3';
+        col.cond.Across = 1 * 'm^2';
+        col.reb.V = 2 * 'm^3';
+        col.reb.Across = 1 * 'm^2';
+
+        col.stage.Qsil = 2;
+
+        # Controllers type
+        TCcond.PID_Select = "Ideal_AW";
+        PC.PID_Select = "Ideal_AW";
+        LCtop.PID_Select = "Ideal_AW";
+        LCbot.PID_Select = "Ideal_AW";
+
+        Qrmax = 5e6 * 'kJ/h';
+        Qrmin = 1e6 * 'kJ/h';
+        Frmin = 0 * 'kmol/h';
+        Frmax = 60 * 'kmol/h';
+        Fcmin = 0 * 'kmol/h';
+        Fcmax = 120 * 'kmol/h';
+        Hmint = 0 * 'm';
+        Hmaxt = 2 * 'm';
+        Hminb = 0 * 'm';
+        Hmaxb = 2 * 'm';
+        Pmin = 0.5 * 'bar';
+        Pmax = 6 * 'bar';
+        Qcmax = -5e5 * 'kJ/h';
+        Qcmin = -5e6 * 'kJ/h';
+        Tmax = (30+273.15) * 'K';
+        Tmin = (-20+273.15) * 'K';
+        
+        INITIAL
+        # condenser
+        col.cond.OutletL.T = 260 *'K';
+        col.cond.Level = 1 * 'm';
+        col.cond.OutletL.z([1:4]) = [0.2, 0.2, 0.4, 0.05];
+
+        # reboiler
+        col.reb.OutletL.T = 350 *'K';
+        col.reb.Level = 1 * 'm';
+        col.reb.OutletL.z([1:4]) = [0.1, 0.4, 0.1, 0.3];
+
+        # column trays
+        col.stage.OutletL.T = [290:(330-290)/(col.NStages-1):330] * 'K';
+        col.stage.ML = 0.1 * 'kmol';
+        col.stage.OutletL.z([1:4]) = [0.15, 0.3, 0.25, 0.2];
+
+        OPTIONS
+        TimeStep = 0.1;
+        TimeEnd = 5;
+        TimeUnit = 'h';
+end

trunk/sample/stage_separators/sample_tray.mso

@@ -85 +85 @@
         TimeEnd = 100;
 end
+
+FlowSheet packedStage_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;
+        
+        DEVICES
+        t1 as packedStage;
+        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;
+
+        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 = 315 * '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 = 315 * 'K';
+        inV.z = [0.0584, 0.9416];
+        
+        t1.OutletV.P = 145 * 'kPa';
+        
+        SET
+        #Metal Pall Ring - nominal packing size 50 mm - Billet and Schultes, 1999.
+        t1.Q = 0 * 'kW';
+        t1.V = 0.8 * 'm^2' * 0.4 * 'm';
+        t1.d = 1.009 * 'm';
+        t1.Cpo = 0.763;
+        t1.e = 0.951;
+        t1.a = 112.6 * 'm^2/m^3';
+        t1.hs = 0.4 * 'm';
+        t1.Qsil = 10;
+        
+        INITIAL
+        t1.OutletL.T = 315 *'K';
+        t1.ML = 0.25 * 'kmol';
+        t1.OutletL.z(1) = 0.1641;
+        
+        OPTIONS
+        DAESolver(File="sundials");
+        TimeStep = 0.001;
+        TimeEnd = 0.1;
+end
+
+FlowSheet packedStage_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;
+
+        DEVICES
+        t1 as packedStage;
+        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;
+
+        SPECIFY
+        feed.Outlet.F = 0 * 'kmol/h';
+        feed.Outlet.T = 300 * 'K';
+        feed.Outlet.P = 1.66 * 'atm';
+        feed.Outlet.z = [0.226, 0.425, 0.035, 0.025, 0.289];
+
+        inL.F = 71.21 * 'kmol/h';       
+        inL.P = 2.22 * 'atm';
+        inL.T = 297.6 * 'K';
+        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.deltaP = 0.01 * 'atm';
+        
+        SET
+        #Metal Pall Ring - nominal packing size 50 mm - Billet and Schultes, 1999.
+        t1.Q = 0 * 'kW';
+        t1.Cpo = 0.763;
+        t1.e = 0.951;
+        t1.a = 112.6 * 'm^2/m^3';
+
+        t1.V = 4 * 'ft^2' * 1 * 'ft';
+        t1.hs = 1 * 'ft';
+        t1.d = 2.26 * 'ft';
+        t1.Qsil = 10;
+
+        INITIAL
+        t1.OutletL.T = 290 *'K';
+        t1.ML = 0.02 * 'kmol';
+        t1.OutletL.z([1:4]) = [0.226, 0.425, 0.035, 0.025];
+        
+        OPTIONS
+        DAESolver(File="sundials");
+        TimeStep = 10;
+        TimeEnd = 100;
+end