Changeset 649 for branches/gui/eml

Timestamp:

Sep 30, 2008, 6:42:41 PM (14 years ago)

Author:

gerson bicca

Message:

testing column model in the GUI (Distillation_thermosyphon_subcooling)

Location:

branches/gui/eml/stage_separators

Files:

• column.mso (modified) (2 diffs)
• condenser.mso (modified) (1 diff)
• reboiler.mso (modified) (1 diff)
• tank.mso (modified) (2 diffs)

branches/gui/eml/stage_separators/column.mso

@@ -447 +447 @@
 *       - a pump in reflux stream.
 *
-* ------------------------------------------------------------------
-Model Distillation_thermosyphon_subcoolingTeste
+* ------------------------------------------------------------------*#
+Model Distillation_thermosyphon_subcooling as Section_ColumnBasic
         ATTRIBUTES
         Pallete         = true;
@@ -478 +478 @@
 ";
         
-        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;
-        reb as reboilerSteady;
-        tbottom as tank;
-        ttop as tank_cylindrical;
-        spbottom as splitter;
-        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";
-                
-                case "off":
-                cond.InletV.F = 0 * 'mol/s';
-                when trays(top).OutletV.P > cond.OutletL.P + 1e-1 * '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 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;
-        trays(bot).OutletL to tbottom.Inlet;
-        tbottom.Outlet to spbottom.Inlet;
-        spbottom.Outlet2 to reb.InletL;
+PARAMETERS
+
+        CondenserVapourFlow as Switcher(Valid = ["on", "off"], Hidden=true, Default = "on");
+
+VARIABLES
+        CondenserUnity as condenserSteady;
+        TopVessel               as tank_cylindrical;
+        TopSplitter             as splitter;
+        PumpUnity               as pump;
+        ReboilerUnity   as reboilerSteady;
+        BottomVessel    as tank;
+        BottomSplitter  as splitter;
+
+
+        alfaTopo                        as Real;
+
+out             HeatToCondenser         as power                                (Brief="Heat supplied to Condenser",Hidden=true);
+out             HeatToReboiler          as power                                (Brief="Heat supplied to Reboiler",Hidden=true);
+out             HeatToBottomVessel      as power                                (Brief="Heat supplied to Bottom Vessel",Hidden=true);
+out             HeatToTopVessel         as power                                (Brief="Heat supplied to Top Vessel",Hidden=true);
+
+EQUATIONS
+
+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":
+                CondenserUnity.InletV.F = 0 * 'mol/s';
+                when trays(1).OutletV.P > CondenserUnity.OutletL.P + 1e-1 * 'atm' switchto "on";
+
+end     
+        
+CONNECTIONS
+#vapor
+        ReboilerUnity.OutletV   to trays(NumberOfTrays).InletV;
+        trays(1).OutletV                                to CondenserUnity.InletV;
+
+#liquid
+        CondenserUnity.OutletL          to TopVessel.Inlet;     
+        TopVessel.Outlet                                        to TopSplitter.Inlet;
+        TopSplitter.Outlet2                             to PumpUnity.Inlet;     
+        PumpUnity.Outlet                                        to trays(1).InletL;
+        trays(NumberOfTrays).OutletL    to BottomVessel.Inlet;
+        BottomVessel.Outlet                             to BottomSplitter.Inlet;
+        BottomSplitter.Outlet2                  to ReboilerUnity.InletL;
+
+#Connectors
+HeatToCondenser                         to CondenserUnity.InletQ;
+HeatToReboiler                  to ReboilerUnity.InletQ;
+HeatToBottomVessel      to BottomVessel.InletQ;
+HeatToTopVessel                 to TopVessel.InletQ;
+
 end
 

branches/gui/eml/stage_separators/condenser.mso

@@ -122 +122 @@
 * the heat supply.
 ";
-        
-        PARAMETERS
-        outer PP as Plugin(Brief = "External Physical Properties", Type="PP");
+
+PARAMETERS
+        outer PP                as Plugin       (Brief = "External Physical Properties", Type="PP");
         outer NComp as Integer;
 
-        VARIABLES
-in      InletV as stream(Brief="Vapour inlet stream", PosX=0.3431, PosY=0, Symbol="_{inV}");
-out     OutletL as liquid_stream(Brief="Liquid outlet stream", PosX=0.34375, PosY=1, Symbol="_{outL}");
-in      InletQ as power (Brief="Cold supplied", PosX=1, PosY=0.5974, Symbol="_{in}");
-        DP as press_delta (Brief="Pressure Drop in the condenser");
-
-        EQUATIONS
-        "Molar Balance"
+VARIABLES
+in      InletV          as stream                               (Brief="Vapour inlet stream", PosX=0.3431, PosY=0, Symbol="_{inV}");
+out     OutletL as liquid_stream                (Brief="Liquid outlet stream", PosX=0.34375, PosY=1, Symbol="_{outL}");
+in      InletQ          as power                                (Brief="Cold supplied", PosX=1, PosY=0.5974, Symbol="_{in}");
+        DP                      as press_delta          (Brief="Pressure Drop in the condenser");
+
+EQUATIONS
+
+"Molar Balance"
         InletV.F = OutletL.F;
         InletV.z = OutletL.z;
-                
-        "Energy Balance"
+
+"Energy Balance"
         InletV.F*InletV.h = OutletL.F*OutletL.h + InletQ;
-        
-        "Pressure"
+
+"Pressure"
         DP = InletV.P - OutletL.P;
+
 end
 

branches/gui/eml/stage_separators/reboiler.mso

@@ -130 +130 @@
 ";      
 
-        PARAMETERS
-        outer PP as Plugin(Brief = "External Physical Properties", Type="PP");
-        outer NComp as Integer;
-        DP as press_delta (Brief="Pressure Drop in the reboiler");
-
-        VARIABLES
-in      InletL as stream(Brief="Liquid inlet stream", PosX=0.3345, PosY=1, Symbol="_{inL}");
-out     OutletV as vapour_stream(Brief="Vapour outlet stream", PosX=0.3369, PosY=0, Symbol="_{outV}");
-in      InletQ as power (Brief="Heat supplied", PosX=1, PosY=0.6111, Symbol="_{in}");
-        vV as volume_mol (Brief="Vapour Molar volume");
-        rhoV as dens_mass (Brief="Vapour Density");
-
-        EQUATIONS
-        "Molar Balance"
+PARAMETERS
+        outer PP                        as Plugin                       (Brief = "External Physical Properties", Type="PP");
+        outer NComp     as Integer;
+        DP                                      as press_delta  (Brief="Pressure Drop in the reboiler");
+
+VARIABLES
+in      InletL                  as stream                                       (Brief="Liquid inlet stream", PosX=0.3345, PosY=1, Symbol="_{inL}");
+out     OutletV         as vapour_stream                (Brief="Vapour outlet stream", PosX=0.3369, PosY=0, Symbol="_{outV}");
+in      InletQ                  as power                                        (Brief="Heat supplied", PosX=1, PosY=0.6111, Symbol="_{in}");
+        vV                              as volume_mol                   (Brief="Vapour Molar volume");
+        rhoV                    as dens_mass                    (Brief="Vapour Density");
+
+EQUATIONS
+
+"Molar Balance"
         InletL.F = OutletV.F;
         InletL.z = OutletV.z;
         
-        "Vapour Volume"
+"Vapour Volume"
         vV = PP.VapourVolume(OutletV.T, OutletV.P, OutletV.z);
         
-        "Vapour Density"
+"Vapour Density"
         rhoV = PP.VapourDensity(OutletV.T, OutletV.P, OutletV.z);
 
-        "Energy Balance"
+"Energy Balance"
         InletL.F*InletL.h + InletQ = OutletV.F*OutletV.h;
         
-        "Pressure"
+"Pressure"
         DP = InletL.P - OutletV.P;
+
 end
 

branches/gui/eml/stage_separators/tank.mso

@@ -53 +53 @@
 ";
 
-        PARAMETERS
-        outer PP as Plugin(Brief = "External Physical Properties", Type="PP");
-        outer NComp as Integer;
-        Across as area (Brief="Tank cross section area", Default=2);
-        
-        VARIABLES
-in      Inlet  as stream (Brief = "Inlet stream", PosX=0.3037, PosY=0, Symbol="_{in}");
-out     Outlet as liquid_stream (Brief = "Outlet liquid stream", PosX=1, PosY=1, Symbol="_{out}");
-in      InletQ as power (Brief="Rate of heat supply", PosX=1, PosY=0.7859, Symbol="_{in}"); 
-        Level    as length(Brief="Tank level");
-        M(NComp) as mol (Brief="Molar Holdup in the tank");
-        E as energy (Brief="Total Energy Holdup on tank");
-        vL as volume_mol (Brief="Liquid Molar Volume");
-
-        EQUATIONS
-        "Mass balance"
+PARAMETERS
+        outer PP                        as Plugin               (Brief = "External Physical Properties", Type="PP");
+        outer NComp     as Integer;
+        Across                          as area                         (Brief="Tank cross section area", Default=2);
+        
+VARIABLES
+in      Inlet                   as stream                               (Brief = "Inlet stream", PosX=0.3037, PosY=0, Symbol="_{in}");
+out     Outlet          as liquid_stream        (Brief = "Outlet liquid stream", PosX=1, PosY=1, Symbol="_{out}");
+in      InletQ                  as power                                (Brief="Rate of heat supply", PosX=1, PosY=0.7859, Symbol="_{in}"); 
+        Level                   as length                               (Brief="Tank level");
+        M(NComp)        as mol                                  (Brief="Molar Holdup in the tank");
+        E                                       as energy                               (Brief="Total Energy Holdup on tank");
+        vL                              as volume_mol           (Brief="Liquid Molar Volume");
+
+EQUATIONS
+
+"Mass balance"
         diff(M) = Inlet.F*Inlet.z - Outlet.F*Outlet.z;
-        
-        "Energy balance"
+
+"Energy balance"
         diff(E) = Inlet.F*Inlet.h - Outlet.F*Outlet.h + InletQ;
 
-        "Energy Holdup"
+"Energy Holdup"
         E = sum(M)*Outlet.h;
 
-        "Mechanical Equilibrium"
+"Mechanical Equilibrium"
         Inlet.P = Outlet.P;
-        
-        "Liquid Volume"
+
+"Liquid Volume"
         vL = PP.LiquidVolume(Outlet.T, Outlet.P, Outlet.z);
-        
-        "Composition"
+
+"Composition"
         M = Outlet.z*sum(M);
-        
-        "Level of liquid phase"
+
+"Level of liquid phase"
         Level = sum(M)*vL/Across;
+
 end
 
@@ -112 +114 @@
 ";
         
-        PARAMETERS
-        outer PP as Plugin(Brief = "External Physical Properties", Type="PP");
-        outer NComp as Integer;
-        radius as length(Brief="Tank radius");
-        L as length(Brief="Tank length");
-        
-        VARIABLES
-in      Inlet  as stream (Brief = "Inlet stream", PosX=0.1825, PosY=0, Symbol="_{in}");
-out     Outlet as liquid_stream (Brief = "Outlet liquid stream", PosX=1, PosY=1, Symbol="_{out}");
-in      InletQ as power (Brief="Rate of heat supply", PosX=1, PosY=0.6160, Symbol="_{in}"); 
-        Level    as length(Brief="Tank level");
-        Across as area (Brief="Tank cross section area", Default=2);
-        M(NComp) as mol (Brief="Molar Holdup in the tank");
-        E as energy (Brief="Total Energy Holdup on tank");
-        vL as volume_mol (Brief="Liquid Molar Volume");
-
-        EQUATIONS
-        "Mass balance"
+PARAMETERS
+        outer PP                        as Plugin               (Brief = "External Physical Properties", Type="PP");
+        outer NComp     as Integer;
+        radius                          as length               (Brief="Tank radius");
+        L                                               as length               (Brief="Tank length");
+        
+VARIABLES
+in      Inlet           as stream                               (Brief = "Inlet stream", PosX=0.1825, PosY=0, Symbol="_{in}");
+out     Outlet  as liquid_stream        (Brief = "Outlet liquid stream", PosX=1, PosY=1, Symbol="_{out}");
+in      InletQ          as power                                (Brief="Rate of heat supply", PosX=1, PosY=0.6160, Symbol="_{in}"); 
+        Level           as length                               (Brief="Tank level");
+        Across          as area                                         (Brief="Tank cross section area", Default=2);
+        M(NComp) as mol                                         (Brief="Molar Holdup in the tank");
+        E                               as energy                               (Brief="Total Energy Holdup on tank");
+        vL                      as volume_mol           (Brief="Liquid Molar Volume");
+
+EQUATIONS
+
+"Mass balance"
         diff(M) = Inlet.F*Inlet.z - Outlet.F*Outlet.z;
         
-        "Energy balance"
+"Energy balance"
         diff(E) = Inlet.F*Inlet.h - Outlet.F*Outlet.h + InletQ;
 
-        "Energy Holdup"
+"Energy Holdup"
         E = sum(M)*Outlet.h;
         
-        "Mechanical Equilibrium"
+"Mechanical Equilibrium"
         Inlet.P = Outlet.P;
         
-        "Liquid Volume"
+"Liquid Volume"
         vL = PP.LiquidVolume(Outlet.T, Outlet.P, Outlet.z);
         
-        "Composition"
+"Composition"
         M = Outlet.z*sum(M);
         
-        "Cylindrical Area"
-        Across = radius^2 * (asin(1) - asin((radius-Level)/radius) ) + 
-                                (Level-radius)*sqrt(Level*(2*radius - Level));
-
-        "Level of liquid phase"
+"Cylindrical Area"
+        Across = radius^2 * (asin(1) - asin((radius-Level)/radius) ) + (Level-radius)*sqrt(Level*(2*radius - Level));
+
+"Level of liquid phase"
         L*Across = sum(M)*vL;
+
 end