Ignore:
Timestamp:
Jan 31, 2008, 3:10:36 PM (15 years ago)
Author:
Paula Bettio Staudt
Message:

Updated packedstage model

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/packed/eml/stage_separators/tray.mso

    r450 r451  
    282282        g as acceleration;
    283283        e as Real (Brief="Void fraction of packing, m^3/m^3");
    284         V as volume;
    285284        Across as area;
    286285        ds as length (Brief="Column diameter");
    287286        d as length (Brief="size of an element of packing");
    288         h as length (Brief="Height of packing");
     287        #h as length (Brief="Height of packing");
    289288        C as Real (Brief="Constant for resitance factor equation");
    290289        Cp as Real (Brief="Constant for resitance at loading point factor equation");
    291290        Mw(NComp)       as molweight    (Brief = "Component Mol Weight");
    292291
    293         SET
    294         Mw = PP.MolecularWeight();
    295        
    296292        VARIABLES
    297293        niL as viscosity (Brief="Liquid dynamic viscosity", DisplayUnit='kg/m/s');
     
    309305        hs as length (Brief="Height of the packing stage");
    310306       
     307        SET
     308        Mw = PP.MolecularWeight();
     309        Ap = Across;
     310       
    311311        EQUATIONS
    312312        "Liquid Density"
     
    318318        "Vapour viscosity"
    319319        niV = PP.VapourViscosity(InletV.T, InletV.P, InletV.z);
     320
    320321        "Liquid volume"
    321322        VL = vL * ML;
    322323        "Liquid holdup"
    323324        hL = VL/V;
     325
    324326        "Volume flow rate of liquid, m^3/m^2/s"
    325327        uL * Across = OutletL.F * vL;
    326         "volume flow rate of vapor, m^3/m^2/s"
     328        "Volume flow rate of vapor, m^3/m^2/s"
    327329        uV * Across = OutletV.F * vV;
    328 
     330       
    329331        "Coefficient of Resistance"
    330332        ksi * C^2 * (uL/uV * sqrt(rhoV/rhoL) * (niL/niV)^5.8)^3 = g/1*'s^2/m';
    331333        ksil = Cp * (exp(uL*rhoL/a/niL/200)*(hL/hLs)^0.3) * (64/Rev+(1.8/Rev)) *
    332                                 ((e-hL)/e);#1.5
     334                                ((e-hL)/e);#^1.5;
     335       
    333336        "Liquid holdup and Liquid flow"
    334337        a^2 * niL * uL = hL^1 *(g*rhoL/3 - ksi*a*rhoV*uV^2/(4*hL*(e-hL)^2));
     338
    335339        "Pressure drop and Vapor flow"
    336340        (InletV.P - OutletV.P)/hs = ksil *(a/2 + 2/ds)*(uV^2*rhoV/(e-hL)^3);
    337341       
    338        
     342        "Vapor Reynolds number"
    339343        Rev = uV * (d-2*hL/a) * rhoV/ niV;
    340        
     344
     345        "Theoretical Liquid Holdup"
    341346        hLs = (12*a^2*niL*uL/g/rhoL)^0.333;
    342347end
Note: See TracChangeset for help on using the changeset viewer.