Changeset 467 for branches/packed


Ignore:
Timestamp:
Feb 28, 2008, 3:49:14 PM (15 years ago)
Author:
Paula Bettio Staudt
Message:

More shots

Location:
branches/packed
Files:
1 added
3 edited

Legend:

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

    r462 r467  
    358358        Mw(NComp)       as molweight    (Brief = "Component Mol Weight");
    359359        hs as length (Brief="Height of the packing stage");
     360        Qsio as Real (Brief="Resistance coefficient", Lower = 0);
    360361       
    361362        VARIABLES
     
    369370        invK as Real (Brief="Wall factor");
    370371        Rev as Real (Brief="Reynolds number of the vapor stream", Lower = 0, Default=100);
    371         Qsio as Real (Brief="Resistance coefficient", Lower = 0);
     372        #Qsio as Real (Brief="Resistance coefficient", Lower = 0);
    372373       
    373374        SET
     
    385386
    386387        "Volume flow rate of liquid, m^3/m^2/s"
    387         uL * Ap = OutletL.F * vL;
     388        uL * Ap * e = OutletL.F * vL;
    388389        "Volume flow rate of vapor, m^3/m^2/s"
    389         uV * Ap = OutletV.F * vV;
     390        uV * Ap * e = OutletV.F * vV;
    390391       
    391392        "Liquid holdup and Liquid flow"
     
    401402        Rev*invK = dp*uV*rhoV / (miV*(1-e));
    402403
    403         if Rev > 1e-4 then
     404#*      if Rev > 1e-4 then
    404405                "Resistance Coefficient"       
    405406                Qsio = Cpo * (64/Rev + 1.8/Rev^0.08);
     
    407408                Qsio = 1;
    408409        end
    409 
     410*#
    410411        "Pressure drop and Vapor flow"
    411         (InletV.P - OutletV.P)*0.0091/hs  = Qsio*a*uV^2*rhoV*invK / (2*e^3);
     412        (InletV.P - OutletV.P)/hs  = Qsio*a*uV^2*rhoV*invK / (2*e^3);
    412413       
    413414        "Efficiency"
  • branches/packed/sample/stage_separators/sample_column.mso

    r464 r467  
    531531        reb to sec.stage(8).InletV;
    532532        cond to sec.stage(1).InletL;
    533 
    534533       
    535534        SPECIFY
     
    556555        reb.z = [0.16, 0.542, 0.013, 0.008, 0.277];
    557556       
    558         #sec.stage.Qsio = 0.8;
    559        
    560         SET
    561         sec.H = 6 * 'ft';
     557        SET
     558        sec.H = 16 * 'm';
    562559        sec.NStages = 8;
    563560        sec.stage.Q = 0 * 'kW';
     
    567564        sec.stage.e = 0.951;
    568565        sec.stage.a = 100*112.6 * 'm^2/m^3';
     566        sec.stage.Qsio = 1;
    569567
    570568        INITIAL
    571569        sec.stage.OutletL.T =[283:(325-283)/(sec.NStages-1):325] *'K';
    572         sec.stage.ML = 0.5 * 'kmol';
     570        sec.stage.ML = 0.3 * 'kmol';
    573571        sec.stage.OutletL.z([1:4]) = [0.2, 0.2, 0.2, 0.2];
    574572
     
    608606        Qr.OutletQ to col.reb.InletQ;
    609607       
     608        VARIABLES
     609        deltaP(col.NStages) as Real (Unit = 'atm/m');
     610       
    610611        SPECIFY
    611612        feed.Outlet.F = 0.98 * 'mol/min';
     
    628629        Qc.OutletQ.Q = -700 * 'cal/min';
    629630        col.pump1.dP = 0.1 * 'atm';
    630         col.stage.Qsio = 0.05;
    631        
    632 #       EQUATIONS
    633 #       col.reb.OutletV.F = col.K * sqrt(Qr.OutletQ.Q);
    634        
    635         SET
    636         col.H = 1 * 'm';
    637         col.NStages = 5;
     631       
     632        EQUATIONS
     633        col.reb.OutletV.F = 0.2 * 'mol*min^0.5/kg^0.5/m' * sqrt(Qr.OutletQ.Q);
     634        deltaP = (col.stage.InletV.P - col.stage.OutletV.P)/col.stage.hs;
     635       
     636        SET
     637        col.H = 2 * 'm';
     638        col.NStages = 2;
    638639
    639640        col.cond.V = 1 * 'l';
     
    650651        col.stage.a = 185.4 * 'm^2/m^3';
    651652        col.stage.Cpo = 0.763;
     653        col.stage.Qsio = 1;
    652654       
    653655        INITIAL
     
    660662        col.reb.OutletL.T = (82+273.15) *'K';
    661663        col.reb.Level = 2 * 'cm';
    662         col.reb.OutletL.z(1) = 0.16;
     664        col.reb.OutletL.z(1) = 0.7;
    663665
    664666        # column stages
    665667        col.stage.OutletL.T = [(63.5+273.15):((82+273.15)-(63.5+273.15))/(col.NStages-1):(82+273.15)] * 'K';
    666         col.stage.Level = 2 * 'cm';
    667         #col.stage.ML = 0.2 * 'mol';
    668         col.stage.OutletL.z(1) = 0.5; #0.16;
     668        #col.stage.Level = 2 * 'cm';
     669        col.stage.ML = 0.2 * 'mol';
     670        col.stage.OutletL.z(1) = 0.5; #[0.79, 0.65, 0.5, 0.25, 0.16];
    669671
    670672        OPTIONS
  • branches/packed/sample/stage_separators/sample_tray.mso

    r462 r467  
    9696end
    9797
    98 FlowSheet packedStage_BilletSchultes_Test
     98FlowSheet packedStage_BilletSchultes_Test_1
     99        PARAMETERS
     100        PP      as Plugin(Brief="Physical Properties",
     101                Type="PP",
     102                Components = [ "n-pentane", "benzene"],
     103                LiquidModel = "PR",
     104                VapourModel = "PR"
     105        );
     106        NComp   as Integer;
     107
     108        SET
     109        NComp = PP.NumberOfComponents;
     110       
     111        VARIABLES
     112        deltaP as Real (Unit='atm/m');
     113       
     114        DEVICES
     115        t1 as packedStage_BilletSchultes;
     116        feed as source;
     117        inL as liquid_stream;
     118        inV as vapour_stream;
     119       
     120        CONNECTIONS
     121        feed.Outlet to t1.Inlet;
     122        inL to t1.InletL;
     123        inV to t1.InletV;
     124
     125        EQUATIONS
     126        deltaP = (t1.InletV.P - t1.OutletV.P)/t1.hs;
     127       
     128        SPECIFY
     129        feed.Outlet.F = 113.4 * 'kmol/h';
     130        feed.Outlet.T = 291 * 'K';
     131        feed.Outlet.P = 1.66 * 'atm';
     132        feed.Outlet.z = [0.5, 0.5];
     133       
     134        inL.P = 165 * 'kPa';
     135        inL.T = 310 * 'K';
     136        inL.F = 61.99 * 'kmol/h';
     137        inL.z = [0.1641, 0.8359];
     138
     139        inV.F = 201.25 * 'kmol/h';
     140        inV.P = 150 * 'kPa';
     141        inV.T = 321 * 'K';
     142        inV.z = [0.0584, 0.9416];
     143
     144        SET
     145        #Metal Pall Ring - nominal packing size 50 mm - Billet and Schultes, 1999.
     146        t1.Q = 0 * 'kW';
     147        t1.Ap = 0.8 * 'm^2';
     148        t1.V = 0.8 * 'm^2' * 1 * 'm';
     149        t1.ds = 1.009 * 'm';
     150        t1.Cpo = 0.763;
     151        t1.e = 0.951;
     152        t1.a = 112.6 * 'm^2/m^3';
     153        t1.hs = 1 * 'm';
     154        t1.Qsio = 1;
     155       
     156        INITIAL
     157        t1.OutletL.T = 320 *'K';
     158        t1.ML = 0.25 * 'kmol';
     159        t1.OutletL.z(1) = 0.5;
     160       
     161        OPTIONS
     162        InitialFile = "packedStage_BilletSchultes_Test_1.rlt";
     163        DAESolver(File="sundials");
     164        TimeStep = 100;
     165        TimeEnd = 1000;
     166end
     167
     168FlowSheet packedStage_BilletSchultes_Test_2
    99169        PARAMETERS
    100170        PP      as Plugin(Brief="Physical Properties",
     
    155225        t1.a = 112.6 * 'm^2/m^3';
    156226        t1.hs = 0.4 * 'm';
     227        t1.Qsio = 1;
    157228
    158229        INITIAL
Note: See TracChangeset for help on using the changeset viewer.