Changeset 721 for branches/gui/eml/stage_separators/condenser.mso

Timestamp:

Feb 25, 2009, 7:51:53 PM (14 years ago)

Author:

gerson bicca

Message:

updates (tank/reboiler/condenser/column)

File:

• branches/gui/eml/stage_separators/condenser.mso (modified) (2 diffs)

branches/gui/eml/stage_separators/condenser.mso

@@ -162 +162 @@
 Model condenserReact
         ATTRIBUTES
-        Pallete         = true;
+        Pallete         = false;
         Icon            = "icon/Condenser"; 
         Brief           = "Model of a Condenser with reaction in liquid phase.";
@@ -183 +183 @@
 ";
         
-        PARAMETERS
-        outer PP as Plugin(Type="PP");
+PARAMETERS
+        outer PP        as Plugin(Type="PP");
         outer NComp as Integer;
-        V as volume (Brief="Condenser total volume");
-        Across as area (Brief="Cross Section Area of reboiler");
-
-        stoic(NComp) as Real(Brief="Stoichiometric matrix");
-        Hr as energy_mol;
-
-        VARIABLES
-in      InletV as stream(Brief="Vapour inlet stream", PosX=0.1164, PosY=0, Symbol="_{inV}");
-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}");
-        InletQ as power (Brief="Cold supplied", PosX=1, PosY=0.6311, Symbol="_{in}");
-
-        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="Level of liquid phase");
-        Vol as volume;
-        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"
+        
+        V               as volume (Brief="Condenser total volume");
+        Across  as area         (Brief="Cross Section Area of reboiler");
+
+        stoic(NComp)    as Real                 (Brief="Stoichiometric matrix");
+        Hr                              as energy_mol;
+        Initial_Level                           as length                       (Brief="Initial Level of liquid phase");
+        Initial_Temperature                     as temperature          (Brief="Initial Temperature of Condenser");
+        Initial_Composition(NComp)      as fraction             (Brief="Initial Liquid Composition");
+        
+VARIABLES
+
+in      InletV          as stream                       (Brief="Vapour inlet stream", PosX=0.1164, PosY=0, Symbol="_{inV}");
+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}");
+        InletQ          as power                        (Brief="Cold supplied", PosX=1, PosY=0.6311, Symbol="_{in}");
+
+        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="Level of liquid phase");
+        Vol             as volume;
+        r3                      as reaction_mol (Brief="Reaction Rates", DisplayUnit = 'mol/l/s');
+        C(NComp)        as conc_mol     (Brief="Molar concentration", Lower = -1);
+
+INITIAL
+
+        Level                                   = Initial_Level;
+        OutletL.T                               = Initial_Temperature;
+        OutletL.z(1:NComp-1)    = Initial_Composition(1:NComp-1)/sum(Initial_Composition);
+
+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) = InletV.F*InletV.z - OutletL.F*OutletL.z
-                                - OutletV.F*OutletV.z + stoic*r3*ML*vL;
-
-        "Energy Balance"
-        diff(E) = InletV.F*InletV.h - OutletL.F*OutletL.h
-                                - OutletV.F*OutletV.h + InletQ + Hr * r3 * ML*vL;
-
-        "Molar Holdup"
+"Component Molar Balance"
+        diff(M) = InletV.F*InletV.z - OutletL.F*OutletL.z - OutletV.F*OutletV.z + stoic*r3*ML*vL;
+
+"Energy Balance"
+        diff(E) = InletV.F*InletV.h - OutletL.F*OutletL.h- OutletV.F*OutletV.h + InletQ + Hr * r3 * ML*vL;
+
+"Molar Holdup"
         M = ML*OutletL.z + MV*OutletV.z; 
         
-        "Energy Holdup"
+"Energy Holdup"
         E = ML*OutletL.h + MV*OutletV.h - OutletV.P*V;
         
-        "Mol fraction normalisation"
+"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"
         OutletL.T = OutletV.T;
 
-        "Mechanical Equilibrium"
+"Mechanical Equilibrium"
         OutletV.P = OutletL.P;
 
-        "Geometry Constraint"
+"Geometry Constraint"
         V = ML*vL + MV*vV;
 
         Vol = ML*vL;
         
-        "Level of liquid phase"
+"Level of liquid phase"
         Level = ML*vL/Across;
         
-        "Chemical Equilibrium"
+"Chemical Equilibrium"
         PP.LiquidFugacityCoefficient(OutletL.T, OutletL.P, OutletL.z)*OutletL.z = 
         PP.VapourFugacityCoefficient(OutletV.T, OutletV.P, OutletV.z)*OutletV.z;