Changeset 802 for branches

Timestamp:

Jul 20, 2009, 12:51:45 PM (14 years ago)

Author:

gerson bicca

Message:

vapour flow equation in condenser model

File:

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

branches/gui/eml/stage_separators/condenser.mso

@@ -47 +47 @@
         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);
+        low_flow        as flow_mol     (Brief = "Low Flow",Default = 1E-6, Hidden=true);
+        zero_flow       as flow_mol     (Brief = "No Flow",Default = 0, Hidden=true);
+        KfConst         as area                 (Brief="Constant for K factor pressure drop", Default = 1, Hidden=true);
+        
+        VapourFlow      as Switcher     (Brief="Vapour Flow", Valid = ["on", "off"], Default = "on",Hidden=true);
+
+        Kfactor as positive (Brief="K factor for pressure drop", Lower = 1E-8, Default = 1E-3);
         Pdrop   as press_delta  (Brief="Pressure Drop in the condenser",Default=0, Symbol="\Delta _P");
+
+SET
+
+        Mw   = PP.MolecularWeight();
+        low_flow = 1E-6 * 'kmol/h';
+        zero_flow = 0 * 'kmol/h';
+        KfConst = 1*'m^2';
 
 VARIABLES
@@ -54 +69 @@
         out     OutletLiquid    as liquid_stream        (Brief="Liquid outlet stream", PosX=0.53, PosY=1, Symbol="_{out}^{Liquid}");
         in      InletQ                  as power                        (Brief="Heat Duty", PosX=1, PosY=0.08, Symbol="Q_{in}",Protected=true);
-
+        
+        rho     as dens_mass    (Brief ="Inlet Vapour Mass Density",Hidden=true, Symbol ="\rho");
+        
         out     TI as control_signal    (Brief="Temperature  Indicator of Condenser", Protected = true, PosX=0.50, PosY=0);
         out     PI as control_signal    (Brief="Pressure  Indicator of Condenser", Protected = true, PosX=0.32, PosY=0);
@@ -60 +77 @@
 EQUATIONS
 
+switch VapourFlow
+
+case "on":
+        InletVapour.F*sum(Mw*InletVapour.z) = Kfactor *sqrt(Pdrop*rho)*KfConst;
+
+        when InletVapour.F < low_flow switchto "off";
+
+case "off":
+        InletVapour.F = zero_flow;
+
+        when InletVapour.P > OutletLiquid.P switchto "on";
+
+end
+
 "Molar Flow Balance"
         InletVapour.F = OutletLiquid.F;
@@ -68 +99 @@
 "Energy Balance"
         InletVapour.F*InletVapour.h = OutletLiquid.F*OutletLiquid.h + InletQ;
+
+"Inlet Vapour Density"
+        rho = PP.VapourDensity(InletVapour.T, InletVapour.P, InletVapour.z);
 
 "Pressure Drop"