Ignore:
Timestamp:
Aug 11, 2009, 11:12:13 AM (13 years ago)
Author:
mamuller
Message:

Added the static head contribution to the OutletLiquid? Pressure for the Dynamic Flash Vessel model.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/gui/eml/stage_separators/flash.mso

    r797 r823  
    6262        Mw(NComp)               as molweight    (Brief="Mol Weight", Hidden=true);
    6363        pi                      as positive             (Brief="Pi value", Default=3.141593,Hidden=true, Symbol="\pi");
     64        g                               as acceleration         (Brief="Gravity Acceleration",Default=9.81,Hidden=true);
    6465       
    6566        Orientation     as Switcher     (Valid=["vertical","horizontal"],Default="vertical");
     
    8586        Vcylinder = 0.25*(pi*Diameter^2)*Lenght;
    8687        radius = 0.5*Diameter;
     88        g = 9.81 * 'm/(s^2)';
    8789
    8890VARIABLES
     
    109111        Pdrop           as press_delta  (Brief = "Pressure Drop", DisplayUnit = 'kPa', Symbol ="\Delta P", Protected=true);
    110112
     113        Peq                     as pressure             (Brief="Equilibrium pressure on the liquid surface", Protected=true, DisplayUnit = 'kPa', Symbol="\Delta P_ª{eq}");
     114        Pstatic         as pressure             (Brief="Static head at the bottom of the tank", Protected = true, DisplayUnit = 'kPa', Symbol="P_{static}^{Liquid}");
     115
    111116out     TI as control_signal    (Brief="Temperature Indicator", PosX=1, PosY=0.39, Protected=true);
    112117out     PI as control_signal    (Brief="Pressure Indicator", PosX=1, PosY=0.21, Protected=true);
     
    236241
    237242"Liquid Volume"
    238         vL = PP.LiquidVolume(OutletLiquid.T, OutletLiquid.P, OutletLiquid.z);
     243        vL = PP.LiquidVolume(OutletLiquid.T, Peq, OutletLiquid.z);
    239244
    240245"Vapour Volume"
    241         vV = PP.VapourVolume(OutletVapour.T, OutletVapour.P, OutletVapour.z);
     246        vV = PP.VapourVolume(OutletVapour.T, Peq, OutletVapour.z);
    242247       
    243248"Chemical Equilibrium"
    244         PP.LiquidFugacityCoefficient(OutletLiquid.T, OutletLiquid.P, OutletLiquid.z)*OutletLiquid.z =
    245                 PP.VapourFugacityCoefficient(OutletVapour.T, OutletVapour.P, OutletVapour.z)*OutletVapour.z;
     249        PP.LiquidFugacityCoefficient(OutletLiquid.T, Peq, OutletLiquid.z)*OutletLiquid.z =
     250                PP.VapourFugacityCoefficient(OutletVapour.T, Peq, OutletVapour.z)*OutletVapour.z;
    246251       
    247252"Thermal Equilibrium"
    248253        OutletVapour.T = OutletLiquid.T;
    249        
    250 "Mechanical Equilibrium"
    251         OutletVapour.P = OutletLiquid.P;
     254
     255"Mechanical Equilibrium for the Vapour Phase"
     256        OutletVapour.P = Peq;
     257       
     258"Static Head"   
     259        Pstatic = PP.LiquidDensity(OutletLiquid.T, Peq, OutletLiquid.z) * g * Level;
     260
     261"Mechanical Equilibrium for the Liquid Phase"
     262        OutletLiquid.P = Peq + Pstatic;
    252263
    253264"Pressure Drop"
Note: See TracChangeset for help on using the changeset viewer.