Changeset 875 for branches/gui/eml/stage_separators/column.mso

Timestamp:

Nov 3, 2009, 12:28:44 PM (13 years ago)

Author:

gerson bicca

Message:

packed column/tray revised

File:

• branches/gui/eml/stage_separators/column.mso (modified) (5 diffs)

branches/gui/eml/stage_separators/column.mso

@@ -1768 +1768 @@
 
         outer PP                as Plugin       (Brief = "External Physical Properties", Type="PP");
-        outer NComp     as Integer      (Brief="Number of components");
-        
-        low_flow                as flow_mol     (Brief ="Low Flow",Default = 1E-6, Hidden=true);
-        low_pressure    as pressure     (Brief ="Low Pressure",Default = 1E-6, Hidden=true);
-        zero_flow               as flow_mol     (Brief ="No Flow",Default = 0, Hidden=true);
-        
+        outer NComp     as Integer      (Brief="Number of components",Protected=true);
+        
+
         NumberOfFeeds                                           as Integer      (Brief="Number of Feed Stages",Default=3,Protected=true);
         NumberOfStages                                          as Integer      (Brief="Number of Stages", Default=3);
@@ -1784 +1781 @@
         ResistanceCoeff         as positive     (Brief="Resistance coefficient on the liquid load", Default=1);
         AreaPerPackingVol       as Real                 (Brief="surface area per packing volume", Unit='m^2/m^3');
-        
-        V                               as volume                       (Brief="Total Volume of the tray",Protected=true);
-        Pi                              as constant             (Brief="Pi Number",Default=3.14159265, Symbol = "\pi",Hidden=true);
-        Gconst                  as acceleration         (Brief="Gravity Acceleration",Default=9.81,Hidden=true);
-        HETP                    as length                       (Brief="the Height Equivalent to a Theoretical Plate",Protected=true);
-        ColumnArea              as area                         (Brief="Column Sectional Cross Area",Protected=true);
-        
-        VapourFlow as Switcher(Valid = ["on", "off"], Default = "on",Hidden=true);
-
+        HETP                            as length               (Brief="The Height Equivalent to a Theoretical Plate",Protected=true);
 
 VARIABLES
@@ -1802 +1791 @@
         out     PCI as control_signal   (Brief="Pressure Indicator", Protected = true, PosX=0, PosY=0.24);
         
-        in      LiquidInlet             as      stream                          (Brief="Liquid Inlet in the section", PosX=0.70, PosY=0);
-        out     VapourOutlet            as vapour_stream        (Brief="Vapour Outlet in the section", PosX=0.30, PosY=0);
-        
-        in              VapourInlet             as stream                               (Brief="Vapour Inlet in the section", PosX=0.30, PosY=1);
-        out     LiquidOutlet            as liquid_stream                (Brief="Liquid Outlet in the section", PosX=0.70, PosY=1);
+        in      LiquidInlet             as      stream                          (Brief="Liquid Inlet in the section", PosX=0.70, PosY=0,Protected=true);
+        out     VapourOutlet            as vapour_stream        (Brief="Vapour Outlet in the section", PosX=0.30, PosY=0,Protected=true);
+        
+        in              VapourInlet             as stream                               (Brief="Vapour Inlet in the section", PosX=0.30, PosY=1,Protected=true);
+        out     LiquidOutlet            as liquid_stream                (Brief="Liquid Outlet in the section", PosX=0.70, PosY=1,Protected=true);
         
         LiquidConnector as stream       (Brief="Liquid connection at the middle STAGES", PosX=0.75, PosY=1,Hidden=true);
         VapourConnector as stream       (Brief="Vapour connection at the middle STAGES", PosX=0.55, PosY=0,Hidden=true);
         
-        STAGES(NumberOfStages) as packedStage;
+        STAGES(NumberOfStages) as packedStage (Brief="Column Stages",Protected=true);
         PressureDrop               as pressure (Brief="Column Pressure Drop", Protected=true);
-        
+
 SET
-        
-        ColumnArea = 0.25*Pi*ColumnDiameter^2;
-        HETP =PackingHeight/NumberOfStages;     
-        V = HETP * ColumnArea;
-        
-        low_pressure = 1E-4 * 'atm';
-        low_flow = 1E-6 * 'kmol/h';
-        zero_flow = 0 * 'kmol/h';
-        
+
+        STAGES.Number_Stages                            = NumberOfStages;
+        STAGES.HeightOfPacking                          = PackingHeight;
+        STAGES.HeatOnStage                                      = HeatSupply;
+        STAGES.ColumnInternalDiameter           = ColumnDiameter;
+        STAGES.PackingVoidFraction                      = VoidFraction;
+        STAGES.LiquidResistanceCoeff            = ResistanceCoeff ;
+        STAGES.AreaPerPackingVolume             = AreaPerPackingVol;
+        
+        HETP = PackingHeight/NumberOfStages;    
+
 CONNECTIONS
 
@@ -1868 +1859 @@
         PressureDrop = STAGES(NumberOfStages).OutletLiquid.P - STAGES(1).OutletLiquid.P;
 
-for i in [1:NumberOfStages] do
-
-switch VapourFlow
-        
-        case "on":
-"Pressure drop and Vapor flow, Billet (4-58)"
-        STAGES(i).deltaP/HETP  = ResistanceCoeff *( 0.5*AreaPerPackingVol + 2/ColumnDiameter) * 1/((VoidFraction-STAGES(i).hl)^3) * (STAGES(i).uV^2) * STAGES(i).rhoV;
-        
-        when STAGES(i).InletVapour.F < low_flow switchto "off";
-        
-        case "off":
-"Vapour Flow"
-        STAGES(i).InletVapour.F = zero_flow;
-        
-        when STAGES(i).deltaP > low_pressure switchto "on";
-
-end
-
-"Energy Balance"
-        diff(STAGES(i).E) = (STAGES(i).Inlet.F*STAGES(i).Inlet.h + STAGES(i).InletLiquid.F*STAGES(i).InletLiquid.h + STAGES(i).InletVapour.F*STAGES(i).InletVapour.h- STAGES(i).OutletLiquid.F*STAGES(i).OutletLiquid.h
-        - STAGES(i).OutletVapour.F*STAGES(i).OutletVapour.h + HeatSupply );
-
-"Energy Holdup"
-        STAGES(i).E = STAGES(i).ML*STAGES(i).OutletLiquid.h + STAGES(i).MV*STAGES(i).OutletVapour.h - STAGES(i).OutletLiquid.P*V;
-        
-"Geometry Constraint"
-        V*VoidFraction= STAGES(i).ML*STAGES(i).vL + STAGES(i).MV*STAGES(i).vV;
-        
-"Volume flow rate of vapor, m^3/m^2/s"
-        STAGES(i).uV * (V*VoidFraction/HETP - STAGES(i).Al) = STAGES(i).InletVapour.F * STAGES(i).vV;
-
-"Liquid holdup"
-        STAGES(i).hl*V*VoidFraction = STAGES(i).ML*STAGES(i).vL;
-        
-"Liquid velocity as a function of liquid holdup, Billet (4-27)"
-        STAGES(i).hl^3 = (12/Gconst) * AreaPerPackingVol^2 * (STAGES(i).miL/STAGES(i).rhoL) * STAGES(i).uL;
-        #STAGES(i).hl = (12*STAGES(i).miL*AreaPerPackingVol^2*STAGES(i).uL/STAGES(i).rhoL/Gconst)^1/3;
-
-"Area occupied by the liquid"
-        STAGES(i).Al = STAGES(i).ML*STAGES(i).vL/HETP;
-
-end
-
 INITIAL
 
@@ -1924 +1872 @@
 
 for i in 1:NComp-1 do
-        
+
 for j in 1:NumberOfStages do