Changeset 506 for branches/rate


Ignore:
Timestamp:
Apr 28, 2008, 3:46:20 PM (14 years ago)
Author:
Gabriela Sporleder Straatmann
Message:

Adding a example for trayRate.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/rate/sample/stage_separators/sample_tray.mso

    r498 r506  
    215215        TimeEnd = 100;
    216216end
     217
     218FlowSheet trayRate_Test
     219        PARAMETERS
     220        PP      as Plugin(Brief="Physical Properties",
     221                Type="PP",
     222                Components = [ "n-pentane", "benzene"],
     223                LiquidModel = "PR",
     224                VapourModel = "PR"
     225        );
     226        NComp   as Integer;
     227
     228        SET
     229        NComp = PP.NumberOfComponents;
     230       
     231        DEVICES
     232        t1 as trayRate;
     233        feed as source;
     234        feedV as source;
     235        inL as liquid_stream;
     236        inV as vapour_stream;
     237       
     238        CONNECTIONS
     239        feed.Outlet to t1.Inlet;
     240        feedV.Outlet to t1.InletFV;
     241        inL to t1.InletL;
     242        inV to t1.InletV;
     243       
     244        SPECIFY
     245        feed.Outlet.F = 113.4 * 'kmol/h';
     246        feed.Outlet.T = 291 * 'K';
     247        feed.Outlet.P = 1.66 * 'atm';
     248        feed.Outlet.z = [0.5, 0.5];
     249       
     250        feedV.Outlet.F = 0 * 'kmol/h';
     251        feedV.Outlet.T = 291 * 'K';
     252        feedV.Outlet.P = 1.66 * 'atm';
     253        feedV.Outlet.z = [0.5, 0.5];
     254       
     255        inL.P = 165 * 'kPa';
     256        inL.T = 310 * 'K';
     257        inL.F = 61.99 * 'kmol/h';
     258        inL.z = [0.1641, 0.8359];
     259
     260        inV.P = 150 * 'kPa';
     261        inV.T = 321 * 'K';
     262        inV.z = [0.0584, 0.9416];
     263        inV.F = 61.99 * 'kmol/h';
     264       
     265        #t1.Emv = 1;
     266        t1.OutletV.F = 147.1 * 'kmol/h';
     267        t1.OutletL.F = 147.1 * 'kmol/h';
     268        t1.interf.a=0.01*'m^2';
     269        t1.interf.htL=0.4*'kW/m^2/K';
     270        t1.interf.htV=0.4*'kW/m^2/K';
     271        t1.interf.kL=0.001*'m/s';
     272        t1.interf.kV=0.001*'m/s';
     273        t1.ML=3*'mol';
     274        t1.MV=3*'mol';
     275
     276        SET
     277        t1.V = 4 * 'ft^3';
     278        #*t1.Ah = 0.394 * 'ft^2';
     279        t1.lw = 20.94 * 'in';
     280        t1.hw = 0.125 * 'ft';
     281        t1.beta = 0.6;
     282        t1.alfa = 4;*#
     283       
     284        t1.Ap = 3.94 * 'ft^2';
     285        t1.Q = 0 * 'kW';
     286       
     287       
     288        INITIAL
     289        t1.OutletL.T = 290 *'K';
     290        t1.OutletV.T = 290 *'K';
     291        #t1.Level = 0.9*'m' ;
     292        t1.OutletL.z(1) = 0.5;
     293        t1.OutletV.z(1) = 0.5;
     294        #t1.MV = 0.01*'mol' ;
     295       
     296        OPTIONS
     297        TimeEnd = 100;
     298end
Note: See TracChangeset for help on using the changeset viewer.