Changeset 468 for branches/packed


Ignore:
Timestamp:
Mar 3, 2008, 3:49:53 PM (15 years ago)
Author:
Paula Bettio Staudt
Message:

More modifications

Location:
branches/packed
Files:
3 edited

Legend:

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

    r467 r468  
    369369        dp as length (Brief="Particle diameter", Default=1e-3);
    370370        invK as Real (Brief="Wall factor");
    371         Rev as Real (Brief="Reynolds number of the vapor stream", Lower = 0, Default=100);
     371        #Rev as Real (Brief="Reynolds number of the vapor stream", Lower = 0, Default=100);
    372372        #Qsio as Real (Brief="Resistance coefficient", Lower = 0);
    373373       
     
    399399        invK = (1 + (2*dp/(3*ds*(1-e))));
    400400       
    401         "Reynolds number of the vapor stream"
    402         Rev*invK = dp*uV*rhoV / (miV*(1-e));
     401        #"Reynolds number of the vapor stream"
     402        #Rev*invK = dp*uV*rhoV / (miV*(1-e));
    403403
    404404#*      if Rev > 1e-4 then
  • branches/packed/sample/stage_separators/sample_column.mso

    r467 r468  
    503503end
    504504
    505 
     505# section based on the Distillation_kettle_cond_Test
    506506FlowSheet SectionColumnPacked
    507507        PARAMETERS
     
    551551
    552552        reb.F = 137.57 * 'kmol/h';
    553         reb.P = 7.46 * 'atm';
     553        reb.P = 2.46 * 'atm';
    554554        reb.T = 325 * 'K';
    555555        reb.z = [0.16, 0.542, 0.013, 0.008, 0.277];
     
    563563        sec.stage.Cpo = 0.763;
    564564        sec.stage.e = 0.951;
    565         sec.stage.a = 100*112.6 * 'm^2/m^3';
     565        sec.stage.a = 112.6 * 'm^2/m^3';
    566566        sec.stage.Qsio = 1;
    567567
     
    574574        #InitialFile = "/home/paula/SectionColumn_Test_with8tray.rlt";
    575575        DAESolver(File="dassl");
    576         TimeStep = 0.01;
    577         TimeEnd = 1;
     576        TimeStep = 1;
     577        TimeEnd = 10;
    578578end
    579579
     
    635635       
    636636        SET
    637         col.H = 2 * 'm';
     637        col.H = 1 * 'm';
    638638        col.NStages = 2;
    639639
  • branches/packed/sample/stage_separators/sample_tray.mso

    r467 r468  
    2525using "stage_separators/tray";
    2626
    27 FlowSheet tray_Test
     27FlowSheet tray_Test_1
    2828        PARAMETERS
    2929        PP      as Plugin(Brief="Physical Properties",
     
    9696end
    9797
     98FlowSheet tray_Test_2
     99        PARAMETERS
     100        PP      as Plugin(Brief="Physical Properties",
     101                Type="PP",
     102                Components = [ "isobutane", "n-pentane", "propylene",
     103                "benzene", "isobutene" ],
     104                LiquidModel = "PR",
     105                VapourModel = "PR"
     106        );
     107        NComp   as Integer;
     108
     109        SET
     110        NComp = PP.NumberOfComponents;
     111       
     112        VARIABLES
     113        deltaP as pressure;
     114        phiL as Real;
     115       
     116        DEVICES
     117        t1 as tray;
     118        feed as source;
     119        inL as liquid_stream;
     120        inV as vapour_stream;
     121       
     122        CONNECTIONS
     123        feed.Outlet to t1.Inlet;
     124        inL to t1.InletL;
     125        inV to t1.InletV;
     126       
     127        EQUATIONS
     128        deltaP = t1.InletV.P - t1.OutletV.P;
     129        phiL = t1.ML*t1.vL / t1.V;
     130       
     131        SPECIFY
     132        feed.Outlet.F = 0 * 'kmol/h';
     133        feed.Outlet.T = 300 * 'K';
     134        feed.Outlet.P = 1.66 * 'atm';
     135        feed.Outlet.z = [0.5, 0.5];
     136       
     137        inL.P = 2.22 * 'atm';
     138        inL.T = 297.6 * 'K';
     139        inL.F = 71.21 * 'kmol/h';
     140        inL.z = [0.226, 0.425, 0.035, 0.025, 0.289];
     141
     142        inV.F = 175.3 * 'kmol/h';
     143        inV.P = 2.3062 * 'atm';
     144        inV.T = 308.3 * 'K';
     145        inV.z = [0.265, 0.233, 0.150, 0.014, 0.338];
     146
     147        t1.Emv = 1;
     148
     149        SET
     150        t1.V = 4 * 'ft^3';
     151        t1.Ah = 0.394 * 'ft^2';
     152        t1.lw = 20.94 * 'in';
     153        t1.hw = 0.125 * 'ft';
     154        t1.Q = 0 * 'kW';
     155        t1.beta = 0.6;
     156        t1.alfa = 4;
     157        t1.Ap = 3.94 * 'ft^2';
     158       
     159        INITIAL
     160        t1.OutletL.T = 304 *'K';
     161        t1.Level = 1.2 * t1.hw;
     162        t1.OutletL.z([1:4]) = [0.2, 0.2, 0.2, 0.2];
     163       
     164        OPTIONS
     165        #InitialFile = "/home/paula/tray_Test.rlt";
     166        TimeStep = 10;
     167        TimeEnd = 100;
     168end
     169
     170#to compare with tray_Test_1
    98171FlowSheet packedStage_BilletSchultes_Test_1
    99172        PARAMETERS
     
    146219        t1.Q = 0 * 'kW';
    147220        t1.Ap = 0.8 * 'm^2';
    148         t1.V = 0.8 * 'm^2' * 1 * 'm';
     221        t1.V = 0.8 * 'm^2' * 0.4 * 'm';
    149222        t1.ds = 1.009 * 'm';
    150223        t1.Cpo = 0.763;
    151224        t1.e = 0.951;
    152225        t1.a = 112.6 * 'm^2/m^3';
    153         t1.hs = 1 * 'm';
     226        t1.hs = 0.4 * 'm';
    154227        t1.Qsio = 1;
    155228       
     
    166239end
    167240
     241#the 3rd stage of the Distillation_kettle_cond_Test
    168242FlowSheet packedStage_BilletSchultes_Test_2
    169243        PARAMETERS
Note: See TracChangeset for help on using the changeset viewer.