Changeset 698 for branches/gui/eml

Timestamp:

Dec 11, 2008, 7:30:42 PM (15 years ago)

Author:

gerson bicca

Message:

added heat_flow model for control purposes - starting to modify the reactive distillation column (BRANCH)

Location:

branches/gui/eml

Files:

• controllers/heat_flow.mso (added)
• controllers/icon/heat_flow.png (added)
• controllers/icon/heat_flow.svg (added)
• stage_separators/column.mso (modified) (2 diffs)
• stage_separators/condenser.mso (modified) (1 diff)
• stage_separators/icon/DistillationReac.png (added)
• stage_separators/icon/DistillationReac.svg (added)
• stage_separators/reboiler.mso (modified) (1 diff)
• stage_separators/tray.mso (modified) (4 diffs)

branches/gui/eml/stage_separators/column.mso

@@ -1317 +1317 @@
 *  Reactive Distillation Column
 *
-* ------------------------------------------------------------------*
+* ------------------------------------------------------------------*#
 Model ReactiveDistillation
         
 ATTRIBUTES
         Pallete         = true;
-        Icon            = "icon/DistillationKettleCond"; 
+        Icon            = "icon/DistillationReac"; 
         Brief           = "Model of a reactive distillation column with dynamic condenser and reboiler.";
         Info            =
@@ -1349 +1349 @@
 ";
         
-        PARAMETERS
+PARAMETERS
+        
         outer PP as Plugin(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");
+
+        FeedTrayIndex(NumberOfTrays)                    as Integer                      (Brief="Number of trays", Default=0,Hidden=true);
+        LiqSideTrayIndex(NumberOfTrays)                 as Integer                      (Brief="Number of trays", Default=0,Hidden=true);
+        VapSideTrayIndex(NumberOfTrays)         as Integer                      (Brief="Number of trays", Default=0,Hidden=true);
+        FeedTrayLocation                                                                        as Integer                      (Brief="Feed tray Location", Default=2);
+        LiquidSideStreamLocation                                                as Integer                      (Brief="Liquid Side Stream Location", Default=2);
+        VapourSideStreamLocation                                        as Integer                      (Brief="Vapour Side Stream Location", Default=2);
+
         alfacond as Real;
 
-        VapourFlow as Switcher(Valid = ["on", "off"], Default = "off");
-        
-        SET
-        top = (NumberOfTrays-1)*(1-topdown)/2+1;
-        bot = NumberOfTrays/top;
-        
-        VARIABLES
-        TRAYS(NumberOfTrays) as trayReactTeste;
-        cond as condenserReact;
-        reb as reboilerReact;
-        sp as splitter;
-        p as pump;
-        
-        EQUATIONS
-        
-        switch VapourFlow
+        VapourFlow                      as Switcher     (Valid = ["on", "off"], Default = "off",Hidden=true);
+        TrayVapourFlow  as Switcher     (Valid = ["on", "off"], Default = "off",Hidden=true);
+        TrayLiquidFlow          as Switcher     (Valid = ["on", "off"], Default = "off",Hidden=true);
+
+        V       as volume                               (Brief="Total Volume of the tray",Hidden=true);
+        Q       as power                                (Brief="Rate of heat supply",Hidden=true); 
+        Ap      as area                                         (Brief="Plate area = Atray - Adowncomer",Hidden=true);
+        Ah      as      area                            (Brief="Total holes area",Hidden=true);
+        lw      as      length                          (Brief="Weir length",Hidden=true);
+        g               as acceleration                 (Brief="Gravity Acceleration",Default=9.81,Hidden=true);
+        hw      as      length                          (Brief="Weir height",Hidden=true);
+        beta as         fraction                        (Brief="Aeration fraction");
+        alfa    as      fraction                        (Brief="Dry pressure drop coefficient");
+
+        VolumeOfTray            as volume               (Brief="Total Volume of the tray");
+        HeatSupply                      as heat_rate    (Brief="Rate of heat supply"); 
+        PlateArea                               as area                         (Brief="Plate area = Atray - Adowncomer");
+        HolesArea                               as area                         (Brief="Total holes area");
+        WeirLength                              as length               (Brief="Weir length");
+        WeirHeight                              as length               (Brief="Weir height");
+
+SET
+
+        FeedTrayIndex(FeedTrayLocation)                 =1;
+        VapSideTrayIndex(FeedTrayLocation)      =1;
+        LiqSideTrayIndex(FeedTrayLocation)      =1;
+
+        V=VolumeOfTray;
+        Q=HeatSupply;
+        Ap=PlateArea;
+        Ah=HolesArea;
+        lw=WeirLength;
+        hw=WeirHeight ;
+        
+VARIABLES
+        
+        in      FeedTray                        as stream                               (Brief="Feed stream", PosX=0, PosY=0.55);
+        
+        VapourDrawOffFlow  as flow_mol          (Brief = "Stream Molar Flow Rate");
+        LiquidDrawOffFlow   as flow_mol                         (Brief = "Stream Molar Flow Rate");
+        
+        out     VapourDrawOff   as vapour_stream                (Brief="Vapour Outlet in the section", PosX=1, PosY=0.46,Protected = true);
+        out     LiquidDrawOff           as liquid_stream                        (Brief="Liquid Outlet in the section", PosX=1, PosY=0.58,Protected = true);
+        
+        RebNoFlow       as sourceNoFlow (Brief="No Inlet Flow to Reboiler",Hidden=true);
+        
+        in              HeatToReboiler                          as power                                (Brief="Heat supplied to Reboiler", PosX=1, PosY=0.93);
+        in              HeatToCondenser                 as power                                (Brief="Heat supplied to Condenser", PosX=1, PosY=0.10);
+
+        TRAYS(NumberOfTrays)    as trayReac;
+        CONDENSER                                       as condenserReact;
+        REBOILER                                                as reboilerReact;
+        SPLITTER                                                        as splitter;
+        PUMP                                                                    as pump;
+
+out     VapourDistillate                                as vapour_stream        (Brief="Vapour outlet stream From Top Condenser", PosX=0.67, PosY=0);
+in              ConnectorCondenserVout  as stream       (Brief="Connector for Vapour outlet stream From Top Condenser", Hidden=true);
+
+out     LiquidDistillate                                as liquid_stream        (Brief="Liquid outlet stream From Top Splitter", PosX=1, PosY=0.33);
+in              ConnectorSplitterOut    as stream       (Brief="Connector for Liquid outlet stream From Top Splitter", Hidden=true);
+
+out     BottomProduct                           as liquid_stream        (Brief="Liquid outlet stream From Reboiler", PosX=1, PosY=0.99);
+in              ConnectorReboilerLout   as stream       (Brief="Connector for Liquid outlet stream From Reboiler", Hidden=true);
+
+
+CONNECTIONS
+#vapor
+        REBOILER.OutletV                                                                to TRAYS(NumberOfTrays).InletV;
+        TRAYS(1).OutletV                                                                        to CONDENSER.InletV;
+        TRAYS([2:NumberOfTrays]).OutletV        to TRAYS([1:NumberOfTrays-1]).InletV;
+
+#liquid
+        CONDENSER.OutletL                                                               to SPLITTER.Inlet;      
+        SPLITTER.Outlet2                                                                        to PUMP.Inlet;
+        PUMP.Outlet                                                                                             to TRAYS(1).InletL;
+        TRAYS([1:NumberOfTrays-1]).OutletL      to TRAYS([2:NumberOfTrays]).InletL;
+        TRAYS(NumberOfTrays).OutletL                    to REBOILER.InletL;
+
+#Connectors
+RebNoFlow.Outlet  to REBOILER.Inlet;
+
+CONDENSER.OutletV to ConnectorCondenserVout;
+SPLITTER.Outlet1                to ConnectorSplitterOut;
+REBOILER.OutletL        to ConnectorReboilerLout;
+
+EQUATIONS
+
+for i in [1:NumberOfTrays] do
+
+"Energy Holdup"
+        TRAYS(i).E = TRAYS(i).ML*TRAYS(i).OutletL.h + TRAYS(i).MV*TRAYS(i).OutletV.h - TRAYS(i).OutletL.P*V;
+
+"Energy Balance"
+        diff(TRAYS(i).E) = ( TRAYS(i).Inlet.F*TRAYS(i).Inlet.h + TRAYS(i).InletL.F*TRAYS(i).InletL.h + TRAYS(i).InletV.F*TRAYS(i).InletV.h- TRAYS(i).OutletL.F*TRAYS(i).OutletL.h - TRAYS(i).OutletV.F*TRAYS(i).OutletV.h
+        -TRAYS(i).VapourSideStream.F*TRAYS(i).VapourSideStream.h - TRAYS(i).LiquidSideStream.F*TRAYS(i).LiquidSideStream.h + Q ) + TRAYS(i).Hr * TRAYS(i).r3 *TRAYS(i). vL*TRAYS(i).ML;
+
+"Level of clear liquid over the weir"
+        TRAYS(i).Level = TRAYS(i).ML*TRAYS(i).vL/Ap;
+
+"Geometry Constraint"
+        V = TRAYS(i).ML* TRAYS(i).vL + TRAYS(i).MV*TRAYS(i).vV;
+
+end
+
+switch VapourFlow
                 case "on":
                 "Pressure Drop through the condenser"
-                cond.InletV.F*TRAYS(top).vV / 'm^2' =
-                        sqrt((TRAYS(top).OutletV.P - cond.OutletL.P + 1e-8 * 'atm')/(TRAYS(top).rhoV*alfacond));
-                when TRAYS(top).OutletV.P < cond.OutletL.P switchto "off";
+                CONDENSER.InletV.F*TRAYS(1).vV / 'm^2' =
+                        sqrt((TRAYS(1).OutletV.P - CONDENSER.OutletL.P + 1e-8 * 'atm')/(TRAYS(1).rhoV*alfacond));
+                when TRAYS(1).OutletV.P < CONDENSER.OutletL.P switchto "off";
                 
                 case "off":
                 "Prato selado"
-                cond.InletV.F = 0.0 * 'mol/s';
-                when TRAYS(top).OutletV.P > cond.OutletL.P + 1e-3 * 'atm' switchto "on";
+                CONDENSER.InletV.F = 0.0 * 'mol/s';
+                when TRAYS(1).OutletV.P > CONDENSER.OutletL.P + 1e-3 * 'atm' switchto "on";
         end
 
-        CONNECTIONS
-        #vapor
-        reb.OutletV to TRAYS(bot).InletV;
-        TRAYS([top+topdown:topdown:bot]).OutletV to TRAYS([top:topdown:bot-topdown]).InletV;
-        TRAYS(top).OutletV to cond.InletV;
-        
-        #liquid
-        cond.OutletL to sp.Inlet;       
-        sp.Outlet2 to p.Inlet;
-        p.Outlet to TRAYS(top).InletL;
-        TRAYS([top:topdown:bot-topdown]).OutletL to TRAYS([top+topdown:topdown:bot]).InletL;
-        TRAYS(bot).OutletL to reb.InletL;
-        
+for i in [1:NumberOfTrays] do
+
+        switch TrayLiquidFlow
+                case "on":
+                "Francis Equation"
+                TRAYS(i).OutletL.F*TRAYS(i).vL = 1.84*'1/s'*lw*((TRAYS(i).Level-(beta*hw)+1e-6*'m')/(beta))^2;
+                when TRAYS(i).Level < (beta * hw) switchto "off";
+                
+                case "off":
+                "Low level"
+                TRAYS(i).OutletL.F = 0 * 'mol/h';
+                when TRAYS(i).Level > (beta * hw) + 1e-6*'m' switchto "on";
+        end
+
+        switch TrayVapourFlow
+                case "on":
+                TRAYS(i).InletV.F*TRAYS(i).vV = sqrt((TRAYS(i).InletV.P - TRAYS(i).OutletV.P - TRAYS(i).Level*g*TRAYS(i).rhoL + 1e-8 * 'atm')/(TRAYS(i).rhoV*alfa))*Ah;
+                when TRAYS(i).InletV.P < TRAYS(i).OutletV.P + TRAYS(i).Level*g*TRAYS(i).rhoL switchto "off";
+                
+                case "off":
+                TRAYS(i).InletV.F = 0 * 'mol/s';
+                when TRAYS(i).InletV.P > TRAYS(i).OutletV.P + TRAYS(i).Level*g*TRAYS(i).rhoL + 3e-2 * 'atm' switchto "on";
+        end
+
+end
+
+# Connecting Trays
+        FeedTray.F*FeedTrayIndex= TRAYS.Inlet.F;
+        FeedTray.T = TRAYS.Inlet.T;
+        FeedTray.P = TRAYS.Inlet.P;
+        FeedTray.z = TRAYS.Inlet.z;
+        FeedTray.v = TRAYS.Inlet.v;
+        FeedTray.h = TRAYS.Inlet.h;
+
+        HeatToReboiler = REBOILER.InletQ;
+        HeatToCondenser         = CONDENSER.InletQ;
+
 end
 

branches/gui/eml/stage_separators/condenser.mso

@@ -187 +187 @@
 out     OutletL as liquid_stream(Brief="Liquid outlet stream", PosX=0.4513, PosY=1, Symbol="_{outL}");
 out     OutletV as vapour_stream(Brief="Vapour outlet stream", PosX=0.4723, PosY=0, Symbol="_{outV}");
-in      InletQ as power (Brief="Cold supplied", PosX=1, PosY=0.6311, Symbol="_{in}");
+        InletQ as power (Brief="Cold supplied", PosX=1, PosY=0.6311, Symbol="_{in}");
 
         M(NComp) as mol (Brief="Molar Holdup in the tray");

branches/gui/eml/stage_separators/reboiler.mso

@@ -253 +253 @@
 out     OutletL as liquid_stream(Brief="Liquid outlet stream", PosX=0.2413, PosY=1, Symbol="_{outL}");
 out     OutletV as vapour_stream(Brief="Vapour outlet stream", PosX=0.5079, PosY=0, Symbol="_{outV}");
-in      InletQ as power (Brief="Heat supplied", PosX=1, PosY=0.6123, Symbol="_{in}");
+        InletQ as power (Brief="Heat supplied", PosX=1, PosY=0.6123, Symbol="_{in}");
 
         M(NComp) as mol (Brief="Molar Holdup in the tray");

branches/gui/eml/stage_separators/tray.mso

@@ -43 +43 @@
 VARIABLES
 
-        Inlet                   as stream                               (Brief="Feed stream", Hidden=true, PosX=0, PosY=0.4932, Symbol="_{in}");
-        LiquidSideStream        as liquid_stream        (Brief="liquid Sidestream", Hidden=true, Symbol="_{outL}");
+        Inlet                                                   as stream                               (Brief="Feed stream", Hidden=true, PosX=0, PosY=0.4932, Symbol="_{in}");
+        LiquidSideStream                as liquid_stream        (Brief="liquid Sidestream", Hidden=true, Symbol="_{outL}");
         VapourSideStream        as vapour_stream        (Brief="vapour Sidestream", Hidden=true, Symbol="_{outV}");
-        
+
 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}");
@@ -61 +61 @@
         Level                                   as length                       (Brief="Height of clear liquid on plate");
         yideal(NComp)           as fraction;
-        
-        
+
+
 EQUATIONS
 "Component Molar Balance"
@@ -161 +161 @@
 * Model of a tray with reaction
 *-------------------------------------------------------------------*#
-Model trayReactTeste
+Model trayReac
         ATTRIBUTES
         Pallete         = false;
@@ -187 +187 @@
 ";
 
-        PARAMETERS
-        outer PP as Plugin(Type="PP");
-        outer NComp as Integer;
-        V as volume(Brief="Total Volume of the tray");
-        Q as power (Brief="Rate of heat supply"); 
-        Ap as area (Brief="Plate area = Atray - Adowncomer");
-        
-        Ah as area (Brief="Total holes area");
-        lw as length (Brief="Weir length");
-        g as acceleration (Default=9.81);
-        hw as length (Brief="Weir height");
-        beta as fraction (Brief="Aeration fraction");
-        alfa as fraction (Brief="Dry pressure drop coefficient");
-
-        stoic(NComp) as Real(Brief="Stoichiometric matrix");
-        Hr as energy_mol;
-        Pstartup as pressure;
-        
-        VapourFlow as Switcher(Valid = ["on", "off"], Default = "off");
-        LiquidFlow as Switcher(Valid = ["on", "off"], Default = "off");
-        
-        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}");
-
-        yideal(NComp) as fraction;
-        Emv as Real (Brief = "Murphree efficiency");
-
-        M(NComp) as mol (Brief="Molar Holdup in the tray");
-        ML as mol (Brief="Molar liquid holdup");
-        MV as mol (Brief="Molar vapour holdup");
-        E as energy (Brief="Total Energy Holdup on tray");
-        vL as volume_mol (Brief="Liquid Molar Volume");
-        vV as volume_mol (Brief="Vapour Molar volume");
-        Level as length (Brief="Height of clear liquid on plate");
-        Vol as volume;
-        
-        rhoL as dens_mass;
-        rhoV as dens_mass;
-        r3 as reaction_mol (Brief = "Reaction resulting ethyl acetate", DisplayUnit = 'mol/l/s');
-        C(NComp) as conc_mol (Brief = "Molar concentration", Lower = -1); #, Unit = "mol/l");
-        
-        EQUATIONS
-        "Molar Concentration"
+PARAMETERS
+
+        outer PP                        as Plugin(Type="PP");
+        outer NComp     as Integer;
+        
+        stoic(NComp)    as Real(Brief="Stoichiometric matrix");
+        Hr                                              as energy_mol;
+        Pstartup                                as pressure;
+        
+VARIABLES
+
+        Inlet                                           as stream                               (Brief="Feed stream", PosX=0, PosY=0.4932, Symbol="_{in}");
+        LiquidSideStream        as liquid_stream        (Brief="liquid Sidestream", Hidden=true, Symbol="_{outL}");
+        VapourSideStream as vapour_stream       (Brief="vapour Sidestream", Hidden=true, Symbol="_{outV}");
+        
+
+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}");
+
+        yideal(NComp)   as fraction;
+        Emv                                     as Real (Brief = "Murphree efficiency");
+
+        M(NComp)                as mol                          (Brief="Molar Holdup in the tray");
+        ML                                              as mol                          (Brief="Molar liquid holdup");
+        MV                                              as mol                          (Brief="Molar vapour holdup");
+        E                                               as energy                       (Brief="Total Energy Holdup on tray");
+        vL                                              as volume_mol   (Brief="Liquid Molar Volume");
+        vV                                              as volume_mol   (Brief="Vapour Molar volume");
+        Level                                   as length                       (Brief="Height of clear liquid on plate");
+        Vol                                             as volume;
+        
+        rhoL                    as dens_mass;
+        rhoV                    as dens_mass;
+        r3                              as reaction_mol         (Brief = "Reaction resulting ethyl acetate", DisplayUnit = 'mol/l/s');
+        C(NComp)        as conc_mol             (Brief = "Molar concentration", Lower = -1); 
+        
+EQUATIONS
+
+"Molar Concentration"
         OutletL.z = vL * C;
         
-        "Reaction"
+"Reaction"
         r3 = exp(-7150*'K'/OutletL.T)*(4.85e4*C(1)*C(2) - 1.23e4*C(3)*C(4))*'l/mol/s';
         
-        "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 + stoic*r3*ML*vL;
-        
-        "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 ) + Hr * r3 * vL*ML;
-        
-        "Molar Holdup"
+"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-
+        LiquidSideStream.F*LiquidSideStream.z-VapourSideStream.F*VapourSideStream.z + stoic*r3*ML*vL;
+
+"Molar Holdup"
         M = ML*OutletL.z + MV*OutletV.z;
-        
-        "Energy Holdup"
-        E = ML*OutletL.h + MV*OutletV.h - OutletL.P*V;
-        
-        "Mol fraction normalisation"
+
+"Thermal Equilibrium Vapour Side Stream"
+        OutletV.T = VapourSideStream.T;
+
+"Thermal Equilibrium Liquid Side Stream"
+        OutletL.T = LiquidSideStream.T;
+
+"Mechanical Equilibrium Vapour Side Stream"
+        OutletV.P= VapourSideStream.P;
+
+"Mechanical Equilibrium Liquid Side Stream"
+        OutletL.P = LiquidSideStream.P;
+
+"Composition Liquid Side Stream"
+        OutletL.z= LiquidSideStream.z;
+        
+"Composition Vapour Side Stream"
+        OutletV.z= VapourSideStream.z;
+
+"Mol fraction normalisation"
         sum(OutletL.z)= 1.0;
-        
-        "Liquid Volume"
+
+"Liquid Volume"
         vL = PP.LiquidVolume(OutletL.T, OutletL.P, OutletL.z);
-        "Vapour Volume"
+
+"Vapour Volume"
         vV = PP.VapourVolume(OutletV.T, OutletV.P, OutletV.z);
 
-        "Thermal Equilibrium"
+"Thermal Equilibrium"
         OutletV.T = OutletL.T;
-        
-        "Mechanical Equilibrium"
+
+"Mechanical Equilibrium"
         OutletV.P = OutletL.P;
-        
-        "Level of clear liquid over the weir"
-        Level = ML*vL/Ap;
 
         Vol = ML*vL;
         
-        "Liquid Density"
+"Liquid Density"
         rhoL = PP.LiquidDensity(OutletL.T, OutletL.P, OutletL.z);
-        "Vapour Density"
+
+"Vapour Density"
         rhoV = PP.VapourDensity(InletV.T, InletV.P, InletV.z);
 
-        switch LiquidFlow
-                case "on":
-                "Francis Equation"
-                OutletL.F*vL = 1.84*'1/s'*lw*((Level-(beta*hw)+1e-6*'m')/(beta))^2;
-                when Level < (beta * hw) switchto "off";
-                
-                case "off":
-                "Low level"
-                OutletL.F = 0 * 'mol/h';
-                when Level > (beta * hw) + 1e-6*'m' switchto "on";
-        end
-
-        switch VapourFlow
-                case "on":
-                #InletV.P = OutletV.P + Level*g*rhoL + rhoV*alfa*(InletV.F*vV/Ah)^2;
-                InletV.F*vV = sqrt((InletV.P - OutletV.P - Level*g*rhoL + 1e-8 * 'atm')/(rhoV*alfa))*Ah;
-                when InletV.P < OutletV.P + Level*g*rhoL switchto "off";
-                
-                case "off":
-                InletV.F = 0 * 'mol/s';
-                when InletV.P > OutletV.P + Level*g*rhoL + 3e-2 * 'atm' switchto "on";
-                #when InletV.P > OutletV.P + Level*beta*g*rhoL + 1e-2 * 'atm' switchto "on";
-        end
-
-        "Chemical Equilibrium"
-        PP.LiquidFugacityCoefficient(OutletL.T, OutletL.P, OutletL.z)*OutletL.z = 
-                PP.VapourFugacityCoefficient(OutletV.T, OutletV.P, yideal)*yideal;
-        
+"Chemical Equilibrium"
+        PP.LiquidFugacityCoefficient(OutletL.T, OutletL.P, OutletL.z)*OutletL.z =       PP.VapourFugacityCoefficient(OutletV.T, OutletV.P, yideal)*yideal;
+
         OutletV.z = Emv * (yideal - InletV.z) + InletV.z;
-        
+
         sum(OutletL.z)= sum(OutletV.z);
-        
-        "Geometry Constraint"
-        V = ML* vL + MV*vV;
+
 end