Changeset 451 for branches/packed


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

Updated packedstage model

Location:
branches/packed
Files:
2 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
  • branches/packed/sample/stage_separators/sample_tray.mso

    r449 r451  
    110110
    111111        VARIABLES
    112         deltaP as Real (Unit='Pa/m');
     112        deltaP as Real (Unit='atm/m');
    113113       
    114114        DEVICES
     
    124124
    125125        EQUATIONS
    126         deltaP = (t1.InletV.P - t1.OutletV.P)/ (t1.V/t1.Across);
     126        deltaP = (t1.InletV.P - t1.OutletV.P)/t1.hs;
    127127
    128128        SPECIFY
     
    137137        inL.z = [0.1641, 0.8359];
    138138
    139         #inV.F = 201.25 * 'kmol/h';
    140         inV.F = 240 * 'kmol/h';
     139        inV.F = 201.25 * 'kmol/h';
    141140        inV.P = 1.48 * 'atm';
    142141        inV.T = 321 * 'K';
     
    144143
    145144        t1.Emv = 1;
     145        t1.hs = 0.075 * 'm';
     146#       t1.OutletV.F = 147.1 * 'kmol/h';
     147#       t1.OutletL.F = 229.5 * 'kmol/h';
     148#       t1.ksi = 0.805;
    146149       
    147150        SET
    148 #       t1.PressureDropModel = "Leva";
    149 #       t1.PackingType = "random";
    150         t1.V = 4 * 'ft^3';
     151        #Ceramic Pall Ring - nominal packing size 50 mm
    151152        t1.Q = 0 * 'kW';
    152         t1.Across = 1 * 'm^2';#3.94 * 'ft^2';
     153        t1.Across = 0.8 * 'm^2';
     154        t1.V = 0.8 * 'm^2' * 0.075 * 'm';
     155        t1.ds = 1.009 * 'm';
     156        t1.d = 50 * 'mm';
     157        t1.C = 2.37;
     158        t1.Cp = 0.662;
     159        t1.e = 0.78;
     160        t1.a = 120 * 'm^2/m^3';
    153161
    154162        INITIAL
    155         t1.OutletL.T = 320 *'K';
    156         t1.Level = 1.2 * 0.125 * 'ft';
    157         t1.OutletL.z(1) = 0.5;
     163        t1.OutletL.T = 350 *'K';
     164        t1.ML = 0.2 * 'mol';
     165        t1.OutletL.z(1) = 0.1;
    158166       
    159167        OPTIONS
    160         InitialFile = "/home/paula/tray_Test.rlt";
     168        #InitialFile = "/home/paula/tray_Test.rlt";
    161169        #DAESolver(File="dassl");
    162170        TimeStep = 0.1;
Note: See TracChangeset for help on using the changeset viewer.