Changeset 768 for branches/gui/eml

Timestamp:

Jun 17, 2009, 6:45:31 PM (14 years ago)

Author:

gerson bicca

Message:

added pressure drop equation

Location:

branches/gui/eml/stage_separators

Files:

• condenser.mso (modified) (5 diffs)
• reboiler.mso (modified) (2 diffs)

branches/gui/eml/stage_separators/condenser.mso

@@ -294 +294 @@
         outer PP                        as Plugin       (Brief = "External Physical Properties", Type="PP");
         outer NComp     as Integer (Brief="Number of Components");
+        
+        Mw(NComp)       as molweight    (Brief = "Component Mol Weight",Hidden=true);
+
+        
+        VapourFlow              as Switcher     (Brief="Vapour Flow", Valid = ["on", "off"], Default = "on",Hidden=true);
 
         V                       as volume       (Brief="Condenser total volume");
         Across  as area                 (Brief="Cross Section Area of condenser");
+        Kfactor as positive     (Brief="K factor for pressure drop", Lower = 1E-8, Default = 1E-3);
         
         Initial_Level                                                           as length                               (Brief="Initial Level of liquid phase");
@@ -313 +319 @@
         out     PCI as control_signal   (Brief="Pressure  Indicator of Condenser", Protected = true, PosX=1, PosY=0.10);
 
-        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");
-
+        M(NComp)        as mol                                  (Brief="Molar Holdup in the tray", Protected = true);
+        ML                              as mol                                  (Brief="Molar liquid holdup", Protected = true);
+        MV                              as mol                                  (Brief="Molar vapour holdup", Protected = true);
+        E                                       as energy                       (Brief="Total Energy Holdup on tray", Protected = true);
+        vL                              as volume_mol   (Brief="Liquid Molar Volume", Protected = true);
+        vV                              as volume_mol   (Brief="Vapour Molar volume", Protected = true);
+        rho                             as dens_mass            (Brief ="Inlet Vapour Mass Density",Hidden=true);
+        Level                   as length                               (Brief="Level of liquid phase", Protected = true);
+        Pdrop                   as press_delta          (Brief = "Pressure Drop", DisplayUnit = 'kPa', Symbol ="\Delta P", Protected=true);
+
+SET
+
+        Mw   = PP.MolecularWeight();
+        
 INITIAL
 
@@ -331 +343 @@
 "Initial Composition"
         OutletLiquid.z(1:NComp-1)   = Initial_Composition(1:NComp-1)/sum(Initial_Composition);
-        
+
 EQUATIONS
+
+switch VapourFlow
+
+        case "on":
+        InletVapour.F*sum(Mw*InletVapour.z) = Kfactor *sqrt(Pdrop*rho)*'m^2';
+
+        when InletVapour.F < 1E-6 * 'kmol/h' switchto "off";
+
+        case "off":
+        InletVapour.F = 0 * 'mol/s';
+
+        when InletVapour.P > OutletLiquid.P switchto "on";
+
+end
+
 "Component Molar Balance"
         diff(M) = InletVapour.F*InletVapour.z - OutletLiquid.F*OutletLiquid.z- OutletVapour.F*OutletVapour.z;
@@ -357 +384 @@
         vV = PP.VapourVolume(OutletVapour.T, OutletVapour.P, OutletVapour.z);
 
+"Inlet Vapour Density"
+        rho = PP.VapourDensity(InletVapour.T, InletVapour.P, InletVapour.z);
+        
 "Chemical Equilibrium"
         PP.LiquidFugacityCoefficient(OutletLiquid.T, OutletLiquid.P, OutletLiquid.z)*OutletLiquid.z = 
@@ -367 +397 @@
         OutletVapour.P = OutletLiquid.P;
 
+"Pressure Drop"
+        OutletLiquid.P  = InletVapour.P - Pdrop;
+
 "Geometry Constraint"
         V = ML*vL + MV*vV;

branches/gui/eml/stage_separators/reboiler.mso

@@ -385 +385 @@
         out     PCI as control_signal   (Brief="Pressure Indicator of Reboiler", Protected = true, PosX=1, PosY=0.10);
         
-        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");
-        rhoV                    as dens_mass            (Brief="Vapour Density");
-        Level                   as length                               (Brief="Level of liquid phase");
+        M(NComp)        as mol                                  (Brief="Molar Holdup in the tray", Protected = true);
+        ML                              as mol                                  (Brief="Molar liquid holdup", Protected = true);
+        MV                              as mol                                  (Brief="Molar vapour holdup", Protected = true);
+        E                                       as energy                       (Brief="Total Energy Holdup on tray", Protected = true);
+        vL                              as volume_mol   (Brief="Liquid Molar Volume", Protected = true);
+        vV                              as volume_mol   (Brief="Vapour Molar volume", Protected = true);
+        rhoV                    as dens_mass            (Brief="Vapour Density", Protected = true);
+        Level                   as length                               (Brief="Level of liquid phase", Protected = true);
+        Pdrop                   as press_delta          (Brief = "Pressure Drop", DisplayUnit = 'kPa', Symbol ="\Delta P", Protected=true);
 
 INITIAL
@@ -437 +438 @@
         OutletLiquid.T = OutletVapour.T;
 
+"Pressure Drop"
+        OutletLiquid.P  = InletLiquid.P - Pdrop;
+
 "Geometry Constraint"
         V = ML*vL + MV*vV;