Changeset 470


Ignore:
Timestamp:
Mar 4, 2008, 6:20:36 PM (15 years ago)
Author:
Paula Bettio Staudt
Message:

More changes

Location:
branches/packed
Files:
2 edited

Legend:

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

    r468 r470  
    365365        rhoL as dens_mass;
    366366        rhoV as dens_mass;
    367         uL as velocity (Brief="volume flow rate of liquid, m^3/m^2/s", Lower=0, Default = 0.007);
    368         uV as velocity (Brief="volume flow rate of vapor, m^3/m^2/s", Lower=0, Default = 1.14);
     367        uL as Real (Brief="volume flow rate of liquid, m^3/m^2/s", Default = 0.007);
     368        uV as Real (Brief="volume flow rate of vapor, m^3/m^2/s", Default = 1.14);
    369369        dp as length (Brief="Particle diameter", Default=1e-3);
    370370        invK as Real (Brief="Wall factor");
     
    570570        V = ML* vL + MV*vV;
    571571end
    572 
    573572
    574573FlowSheet test
  • branches/packed/sample/stage_separators/sample_tray.mso

    r469 r470  
    112112        VARIABLES
    113113        deltaP as pressure (Lower=-10);
    114         phiL as Real;
     114        phiL as volume;
    115115       
    116116        DEVICES
     
    127127        EQUATIONS
    128128        deltaP = t1.InletV.P - t1.OutletV.P;
    129         phiL = t1.ML*t1.vL / t1.V;
     129        phiL = t1.ML*t1.vL;
    130130       
    131131        SPECIFY
     
    257257        VARIABLES
    258258        deltaP as Real (Unit='atm/m');
     259        phiL as volume;
    259260       
    260261        DEVICES
     
    271272        EQUATIONS
    272273        deltaP = (t1.InletV.P - t1.OutletV.P)/t1.hs;
     274        phiL = t1.vL * t1.ML;
    273275
    274276        SPECIFY
    275277        feed.Outlet.F = 0 * 'kmol/h';
    276278        feed.Outlet.T = 300 * 'K';
    277         feed.Outlet.P = 1 * 'atm';
    278         feed.Outlet.z = 1/NComp;
     279        feed.Outlet.P = 1.66 * 'atm';
     280        feed.Outlet.z = [0.2, 0.2, 0.2, 0.2, 0.2];
    279281
    280282        inL.F = 71.21 * 'kmol/h';       
     
    293295        #Metal Pall Ring - nominal packing size 50 mm - Billet and Schultes, 1999.
    294296        t1.Q = 0 * 'kW';
    295         t1.Ap = 0.8 * 'm^2';
    296         t1.V = 0.8 * 'm^2' * 0.4 * 'm';
    297         t1.ds = 1.009 * 'm';
     297        #t1.Ap = 0.8 * 'm^2';
     298        #t1.V = 0.8 * 'm^2' * 0.4 * 'm';
     299#       t1.ds = 1.009 * 'm';
    298300        t1.Cpo = 0.763;
    299301        t1.e = 0.951;
    300302        t1.a = 112.6 * 'm^2/m^3';
    301         t1.hs = 0.4 * 'm';
     303#       t1.hs = 0.4 * 'm';
    302304        t1.Qsio = 1;
    303305
     306        t1.V = 4 * 'ft^2' * 1 * 'ft';
     307        t1.Ap = 4 * 'ft^2';
     308        t1.hs = 1 * 'ft';
     309        t1.ds = 2.26 * 'ft';
     310
    304311        INITIAL
    305         t1.OutletL.T = 280 *'K';
    306         t1.ML = 0.05 * 'kmol';
    307         t1.OutletL.z([1:4]) = [0.286, 0.203, 0.148, 0.006];
     312        t1.OutletL.T = 304 *'K';
     313        t1.ML = 0.18 * 'kmol';
     314        t1.OutletL.z([1:4]) = [0.2, 0.2, 0.2, 0.2];
    308315       
    309316        OPTIONS
    310317        DAESolver(File="sundials");
    311318        #InitialFile = "packedStage_BilletSchultes_Test.rlt";
    312         TimeStep = 2;   
    313         TimeEnd = 50;
     319        InitialFile = "/home/paula/tray_Test_2.rlt";
     320        TimeStep = 5;   
     321        TimeEnd = 60;
    314322end
Note: See TracChangeset for help on using the changeset viewer.