Ignore:
Timestamp:
Jun 17, 2009, 6:45:31 PM (13 years ago)
Author:
gerson bicca
Message:

added pressure drop equation

File:
1 edited

Legend:

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

    r767 r768  
    294294        outer PP                        as Plugin       (Brief = "External Physical Properties", Type="PP");
    295295        outer NComp     as Integer (Brief="Number of Components");
     296       
     297        Mw(NComp)       as molweight    (Brief = "Component Mol Weight",Hidden=true);
     298
     299       
     300        VapourFlow              as Switcher     (Brief="Vapour Flow", Valid = ["on", "off"], Default = "on",Hidden=true);
    296301
    297302        V                       as volume       (Brief="Condenser total volume");
    298303        Across  as area                 (Brief="Cross Section Area of condenser");
     304        Kfactor as positive     (Brief="K factor for pressure drop", Lower = 1E-8, Default = 1E-3);
    299305       
    300306        Initial_Level                                                           as length                               (Brief="Initial Level of liquid phase");
     
    313319        out     PCI as control_signal   (Brief="Pressure  Indicator of Condenser", Protected = true, PosX=1, PosY=0.10);
    314320
    315         M(NComp)        as mol                                  (Brief="Molar Holdup in the tray");
    316         ML                              as mol                                  (Brief="Molar liquid holdup");
    317         MV                              as mol                                  (Brief="Molar vapour holdup");
    318         E                                       as energy                       (Brief="Total Energy Holdup on tray");
    319         vL                              as volume_mol   (Brief="Liquid Molar Volume");
    320         vV                              as volume_mol   (Brief="Vapour Molar volume");
    321         Level                   as length                               (Brief="Level of liquid phase");
    322 
     321        M(NComp)        as mol                                  (Brief="Molar Holdup in the tray", Protected = true);
     322        ML                              as mol                                  (Brief="Molar liquid holdup", Protected = true);
     323        MV                              as mol                                  (Brief="Molar vapour holdup", Protected = true);
     324        E                                       as energy                       (Brief="Total Energy Holdup on tray", Protected = true);
     325        vL                              as volume_mol   (Brief="Liquid Molar Volume", Protected = true);
     326        vV                              as volume_mol   (Brief="Vapour Molar volume", Protected = true);
     327        rho                             as dens_mass            (Brief ="Inlet Vapour Mass Density",Hidden=true);
     328        Level                   as length                               (Brief="Level of liquid phase", Protected = true);
     329        Pdrop                   as press_delta          (Brief = "Pressure Drop", DisplayUnit = 'kPa', Symbol ="\Delta P", Protected=true);
     330
     331SET
     332
     333        Mw   = PP.MolecularWeight();
     334       
    323335INITIAL
    324336
     
    331343"Initial Composition"
    332344        OutletLiquid.z(1:NComp-1)   = Initial_Composition(1:NComp-1)/sum(Initial_Composition);
    333        
     345
    334346EQUATIONS
     347
     348switch VapourFlow
     349
     350        case "on":
     351        InletVapour.F*sum(Mw*InletVapour.z) = Kfactor *sqrt(Pdrop*rho)*'m^2';
     352
     353        when InletVapour.F < 1E-6 * 'kmol/h' switchto "off";
     354
     355        case "off":
     356        InletVapour.F = 0 * 'mol/s';
     357
     358        when InletVapour.P > OutletLiquid.P switchto "on";
     359
     360end
     361
    335362"Component Molar Balance"
    336363        diff(M) = InletVapour.F*InletVapour.z - OutletLiquid.F*OutletLiquid.z- OutletVapour.F*OutletVapour.z;
     
    357384        vV = PP.VapourVolume(OutletVapour.T, OutletVapour.P, OutletVapour.z);
    358385
     386"Inlet Vapour Density"
     387        rho = PP.VapourDensity(InletVapour.T, InletVapour.P, InletVapour.z);
     388       
    359389"Chemical Equilibrium"
    360390        PP.LiquidFugacityCoefficient(OutletLiquid.T, OutletLiquid.P, OutletLiquid.z)*OutletLiquid.z =
     
    367397        OutletVapour.P = OutletLiquid.P;
    368398
     399"Pressure Drop"
     400        OutletLiquid.P  = InletVapour.P - Pdrop;
     401
    369402"Geometry Constraint"
    370403        V = ML*vL + MV*vV;
Note: See TracChangeset for help on using the changeset viewer.