Changeset 666 for branches

Timestamp:

Oct 13, 2008, 2:57:59 PM (14 years ago)

Author:

gerson bicca

Message:

added column model

Location:

branches/gui

Files:

• eml/stage_separators/column.mso (modified) (2 diffs)
• sample/stage_separators/Column_Rectifier_Subcooling_Diagram.pfd (added)

branches/gui/eml/stage_separators/column.mso

@@ -893 +893 @@
 *       - a pump in reflux stream.
 *
-* ------------------------------------------------------------------*
-Model Rectifier_subcoolingTeste
-        ATTRIBUTES
+* ------------------------------------------------------------------*#
+
+Model Rectifier_subcooling as Section_ColumnBasic
+
+ATTRIBUTES
         Pallete         = true;
         Icon            = "icon/RefluxedSubcooling"; 
@@ -920 +922 @@
 ";
         
-        PARAMETERS
-        outer PP as Plugin(Brief = "External Physical Properties", Type="PP");
-        outer NComp as Integer;
-        NumberOfTrays 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");
-        VapourFlow as Switcher(Valid = ["on", "off"], Default = "on");
-
-        SET
-        top = (NumberOfTrays-1)*(1-topdown)/2+1;
-        bot = NumberOfTrays/top;
-        
-        VARIABLES
-        trays(NumberOfTrays) as trayTeste;
-        cond as condenserSteady;
-        ttop as tank_cylindrical;
-        sptop as splitter;
-        pump1 as pump;
-        alfaTopo as Real;
-
-        EQUATIONS
-        switch VapourFlow
-                case "on":
-                cond.InletV.F*trays(top).vV = alfaTopo * trays(top).Ah * sqrt(2*(trays(top).OutletV.P -
-                cond.OutletL.P + 1e-8 * 'atm') / (trays(top).alfa*trays(top).rhoV));
-                when cond.InletV.F < 1e-6 * 'kmol/h' switchto "off";
+PARAMETERS
+        CondenserVapourFlow             as Switcher     (Valid = ["on", "off"], Default = "on",Hidden=true);
+
+VARIABLES
+        CondenserUnity as condenserSteady;
+        TopVessel               as tank_cylindrical;
+        SplitterTop             as splitter;
+        PumpUnity               as pump;
+        alfaTopo                        as Real;
+
+        in              VapourInlet             as stream                               (Brief="Vapour Inlet in the section", PosX=0.18, PosY=1);
+        out     LiquidOutlet            as liquid_stream                (Brief="Liquid Outlet in the section", PosX=0.25, PosY=1);
+
+        out             HeatToCondenser         as power                                (Brief="Heat supplied to Condenser",Hidden=true);
+        out     LiquidDistillate                                as liquid_stream                (Brief="Liquid outlet stream From Top Splitter", PosX=1, PosY=0.47);
+        in              ConnectorSplitterTop            as stream                               (Brief="Connector for Liquid outlet stream From Top Splitter", Hidden=true);
+        out             HeatToTopVessel         as power                                (Brief="Heat supplied to Top Vessel",Hidden=true);
+
+        VapourConnector as stream               (Brief="Vapour connection at the middle trays", PosX=0.55, PosY=0,Hidden=true);
+
+EQUATIONS
+
+        LiquidOutlet.F= trays(NumberOfTrays).OutletL.F;
+        LiquidOutlet.T = trays(NumberOfTrays).OutletL.T;
+        LiquidOutlet.P = trays(NumberOfTrays).OutletL.P;
+        LiquidOutlet.z = trays(NumberOfTrays).OutletL.z;
+
+        VapourConnector.F= VapourInlet.F;
+        VapourConnector.T = VapourInlet.T;
+        VapourConnector.P = VapourInlet.P;
+        VapourConnector.z = VapourInlet.z;
+        VapourConnector.v = VapourInlet.v;
+        VapourConnector.h = VapourInlet.h;
+        
+# Splitter Connector Equations
+        ConnectorSplitterTop.T = LiquidDistillate.T;
+        ConnectorSplitterTop.P = LiquidDistillate.P;
+        ConnectorSplitterTop.F = LiquidDistillate.F;
+        ConnectorSplitterTop.z = LiquidDistillate.z;
+
+switch CondenserVapourFlow
+
+        case "on":
+                CondenserUnity.InletV.F*trays(1).vV = alfaTopo *Ah * sqrt(2*(trays(1).OutletV.P -
+                CondenserUnity.OutletL.P + 1e-8 * 'atm') / (alfa*trays(1).rhoV));
+        when CondenserUnity.InletV.F < 1e-6 * 'kmol/h' switchto "off";
                 
-                case "off":
-                cond.InletV.F = 0 * 'mol/s';
-                when trays(top).OutletV.P > cond.OutletL.P + 1e-1 * 'atm' switchto "on";
-        end     
-        
-        CONNECTIONS
-        #vapor
-        trays([top+topdown:topdown:bot]).OutletV to trays([top:topdown:bot-topdown]).InletV;
-        trays(top).OutletV to cond.InletV;
-        
-        #liquid
-        cond.OutletL to ttop.Inlet;     
-        ttop.Outlet to sptop.Inlet;
-        sptop.Outlet2 to pump1.Inlet;   
-        pump1.Outlet to trays(top).InletL;
-        trays([top:topdown:bot-topdown]).OutletL to trays([top+topdown:topdown:bot]).InletL;
+        case "off":
+                CondenserUnity.InletV.F = 0 * 'mol/s';
+        when trays(1).OutletV.P > CondenserUnity.OutletL.P + 1e-1 * 'atm' switchto "on";
+
+end     
+        
+CONNECTIONS
+#vapor
+        trays(1).OutletV to CondenserUnity.InletV;
+        
+#liquid
+        CondenserUnity.OutletL to TopVessel.Inlet;      
+        TopVessel.Outlet                        to SplitterTop.Inlet;
+        SplitterTop.Outlet2             to PumpUnity.Inlet;     
+        PumpUnity.Outlet                        to trays(1).InletL;
+
+#Connectors
+        VapourConnector                 to trays(NumberOfTrays).InletV;
+        HeatToCondenser                 to CondenserUnity.InletQ;
+        SplitterTop.Outlet1     to ConnectorSplitterTop;
+        HeatToTopVessel         to TopVessel.InletQ;
+        
 end