Changeset 458


Ignore:
Timestamp:
Feb 14, 2008, 5:09:53 PM (15 years ago)
Author:
Paula Bettio Staudt
Message:

More shots

Location:
branches/packed
Files:
3 edited

Legend:

Unmodified
Added
Removed

  • branches/packed/eml/stage_separators/tray.mso

    r455 r458  
    366366        uL as velocity (Brief="volume flow rate of liquid, m^3/m^2/s", Default = 0.007);
    367367        uV as velocity (Brief="volume flow rate of vapor, m^3/m^2/s", Default = 1.14);
    368         dp as length (Brief="Particle diameter");
    369         K as Real (Brief="Wall factor");
    370         Rev as Real (Brief="Reynolds number of the vapor stream");
    371         Qsio as Real (Brief="Resistance coefficient");
     368        dp as length (Brief="Particle diameter", Default=1e-3);
     369        invK as Real (Brief="Wall factor");
     370        Rev as Real (Brief="Reynolds number of the vapor stream", Lower = 0, Default=100);
     371        Qsio as Real (Brief="Resistance coefficient", Lower = 0);
    372372       
    373373        SET
    … …  
    393393       
    394394        "Particle diameter"
    395         dp = 6 * ((1-e)/a);
     395        dp = 6 * (1-e)/a;
    396396       
    397397        "Wall Factor"
    398         1 = K * (1 + (2*dp/(3*ds*(1-e))));
     398        invK = (1 + (2*dp/(3*ds*(1-e))));
    399399       
    400400        "Reynolds number of the vapor stream"
    401         Rev = K*dp*uV*rhoV / (miV*(1-e));
    402        
    403         #"Resistance Coefficient"
    404         #Qsio = Cpo * (64/Rev + 1.8/Rev^0.08);
    405         #Qsio = Cpo * (64/Rev);# + 1.8/Rev^0.08);
     401        Rev*invK = dp*uV*rhoV / (miV*(1-e));
     402       
     403        "Resistance Coefficient"
     404        Qsio = Cpo * (64/Rev + 1.8/Rev^0.08);
    406405
    407406        "Pressure drop and Vapor flow"
    408         (InletV.P - OutletV.P)/hs * 'm'/'Pa' = Qsio*a*uV^2*rhoV / (2*e^3*K) * 's^2'*'m^2'/'kg';
     407        (InletV.P - OutletV.P)/hs  = Qsio*a*uV^2*rhoV*invK / (2*e^3);
    409408       
    410409        "Efficiency"
    … …  
    414413#*-------------------------------------------------------------------
    415414* Model of a tray with reaction
    416 *-------------------------------------------------------------------*#
     415*------------------------------------------------------------------*#
    417416Model trayReact
    418417        ATTRIBUTES

  • branches/packed/sample/stage_separators/sample_column.mso

    r455 r458  
    9999       
    100100        OPTIONS
    101         TimeStep = 10;
     101        TimeStep = 100;
    102102        TimeEnd = 1000;
    103103       
    … …  
    277277
    278278        OPTIONS
    279         TimeStep = 0.1;
     279        TimeStep = 5;
    280280        TimeEnd = 50;
    281281        #GuessFile="Distillation_kettle_cond_Test.rlt";
    … …  
    508508        PP      as Plugin(Brief="Physical Properties",
    509509                Type="PP",
    510                 Components = [ "isobutene", "n-pentane", "propylene", "benzene", "isobutane"],
     510                Components = [ "isobutane", "n-pentane", "propylene",
     511                "benzene", "isobutene" ],
    511512                LiquidModel = "PR",
    512513                VapourModel = "PR"
    … …  
    519520        DEVICES
    520521        sec as Section_Column_Packed;
    521         feed as source;
     522        feed as liquid_stream;
    522523        reb as vapour_stream;
    523524        cond as liquid_stream;
    … …  
    525526       
    526527        CONNECTIONS
    527         feed.Outlet to sec.stage(10).Inlet;
    528        
    529         zero to sec.stage(1:9).Inlet;
    530         zero to sec.stage([11:sec.NStages]).Inlet;
    531        
    532         reb to sec.stage(sec.NStages).InletV;
     528        feed to sec.stage(5).Inlet;
     529        zero to sec.stage([1:4]).Inlet;
     530        zero to sec.stage([6:sec.NStages]).Inlet;
     531        reb to sec.stage(8).InletV;
    533532        cond to sec.stage(1).InletL;
    534533
    535         VARIABLES
    536         deltaP_TOP as Real (Unit='atm/m');
    537         deltaP_BOT       as Real (Unit='atm/m');
    538        
    539         EQUATIONS
    540         deltaP_TOP = (sec.stage(1).InletV.P - sec.stage(1).OutletV.P)/sec.stage(1).hs;
    541         deltaP_BOT = (sec.stage(sec.NStages).InletV.P - sec.stage(sec.NStages).OutletV.P)/sec.stage(sec.NStages).hs;
    542534       
    543535        SPECIFY
    544         feed.Outlet.F = 113.4 * 'kmol/h';
    545         feed.Outlet.T = 291 * 'K';
    546         feed.Outlet.P = 168.3 * 'kPa';
    547         feed.Outlet.z = 1/NComp;
     536        feed.F = 113.4 * 'kmol/h';
     537        feed.T = 291 * 'K';
     538        feed.P = 1.66 * 'atm';
     539        feed.z = [0.2, 0.2, 0.2, 0.2, 0.2];
    548540
    549541        zero.F = 0 * 'kmol/h';
    550542        zero.T = 300 * 'K';
    551543        zero.P = 1 * 'atm';
    552         zero.z = 1/NComp;
     544        zero.z = [0.2, 0.2, 0.2, 0.2, 0.2];
    553545        zero.v = 0;
    554546        zero.h = 0 * 'J/mol';
    555547       
    556548        cond.F = 85 * 'kmol/h';
    557         cond.P = 2.2 * 'atm';
    558         cond.T = 273 * 'K';
    559         cond.z = [0.26682, 0.19950, 0.26682, 0.000051, 0.26680];
    560 
    561         reb.F = 143.4 * 'kmol/h';
    562         reb.P = 5 * 'atm';
    563         reb.T = 361 * 'K';
    564         reb.z = [0.0001611, 0.51834, 0.00000076, 0.48115, 0.00035];
    565        
    566         #sec.stage.Qsio = 1;
    567        
    568         SET
    569         sec.H = 4 * 'm';
    570         sec.NStages = 16;
     549        cond.P = 2.33 * 'atm';
     550        cond.T = 283.5 * 'K';
     551        cond.z = [0.599, 0.044, 0.035, 0.007, 0.315];
     552
     553        reb.F = 137.57 * 'kmol/h';
     554        reb.P = 2.46 * 'atm';
     555        reb.T = 325 * 'K';
     556        reb.z = [0.16, 0.542, 0.013, 0.008, 0.277];
     557       
     558        #sec.stage.Qsio = 0.8;
     559       
     560        SET
     561        sec.H = 6 * 'ft';
     562        sec.NStages = 8;
    571563        sec.stage.Q = 0 * 'kW';
    572         sec.stage.Ap = 1.5 * 'm^2';
     564        sec.stage.Ap = 3.94 * 'ft^2';
    573565        sec.stage.ds = 1.009 * 'm';
    574566        sec.stage.Cpo = 0.763;
    … …  
    577569
    578570        INITIAL
    579         sec.stage.OutletL.T = [273:(360-273)/(sec.NStages-1):360] *'K';
     571        sec.stage.OutletL.T =[283:(325-283)/(sec.NStages-1):325] *'K';
    580572        sec.stage.ML = 0.5 * 'kmol';
    581         sec.stage.OutletL.z([1:4]) = 1/NComp;
    582        
     573        sec.stage.OutletL.z([1:4]) = [0.2, 0.2, 0.2, 0.2];
     574
    583575        OPTIONS
     576        #InitialFile = "/home/paula/SectionColumn_Test_with8tray.rlt";
    584577        DAESolver(File="dassl");
    585578        TimeStep = 0.01;
    … …  
    637630        col.stage.Qsio = 0.05;
    638631       
    639         EQUATIONS
    640         col.reb.OutletV.F = col.K * sqrt(Qr.OutletQ.Q);
    641        
    642         SET
     632#       EQUATIONS
     633#       col.reb.OutletV.F = col.K * sqrt(Qr.OutletQ.Q);
     634       
     635        SET
     636        col.H = 1 * 'm';
    643637        col.NStages = 5;
     638
    644639        col.cond.V = 1 * 'l';
    645640        col.cond.Across = 100 * 'cm^2';
    646         col.K = 0.2 * 'mol*min^0.5/kg^0.5/m';
     641#       col.K = 0.2 * 'mol*min^0.5/kg^0.5/m';
     642       
    647643        col.reb.V = 2 * 'l';
    648644        col.reb.Across = 200 * 'cm^2';
     645
    649646        col.stage.Q = 0 * 'kW';
    650647        col.stage.Ap = 56.74 * 'cm^2';
    … …  
    652649        col.stage.e = 0.662;
    653650        col.stage.a = 185.4 * 'm^2/m^3';
    654         col.H = 1 * 'm';
    655651        col.stage.Cpo = 0.763;
    656652       

  • branches/packed/sample/stage_separators/sample_tray.mso

    r455 r458  
    9696end
    9797
    98 FlowSheet packedStage_Test
    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_Billet;
    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.F = 61.99 * 'kmol/h';       
    135         inL.P = 1.63 * 'atm';
    136         inL.T = 310 * 'K';
    137         inL.z = [0.1641, 0.8359];
    138 
    139         inV.F = 201.25 * 'kmol/h';
    140         inV.P = 1.48 * 'atm';
    141         inV.T = 321 * 'K';
    142         inV.z = [0.0584, 0.9416];
    143 
    144         t1.Emv = 1;
    145         t1.hs = 0.4 * 'm';#0.075 * 'm';
    146         t1.OutletV.F = 147.1 * 'kmol/h';
    147 #       t1.OutletL.F = 229.5 * 'kmol/h';
    148         #t1.ksi = 0.805;
    149        
    150         SET
    151         #Ceramic Pall Ring - nominal packing size 50 mm
    152         t1.Q = 0 * 'kW';
    153         t1.Across = 0.8 * 'm^2';
    154         t1.V = 0.8 * 'm^2' * 0.4 * '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';
    161 
    162         INITIAL
    163         t1.OutletL.T = 350 *'K';
    164         t1.ML = 0.5 * 'mol';
    165         t1.OutletL.z(1) = 0.2;
    166        
    167         OPTIONS
    168         #InitialFile = "/home/paula/tray_Test.rlt";
    169         #InitialFile = "/home/paula/packedStage_Test.rlt";
    170         TimeStep = 2;   
    171         TimeEnd = 50;
    172 end
    173 
    17498FlowSheet packedStage_BilletSchultes_Test
    17599        PARAMETERS
    176100        PP      as Plugin(Brief="Physical Properties",
    177101                Type="PP",
    178                 Components = [ "n-pentane", "benzene"],
     102                Components = [ "isobutane", "n-pentane", "propylene",
     103                "benzene", "isobutene" ],
    179104                LiquidModel = "PR",
    180105                VapourModel = "PR"
    … …  
    203128
    204129        SPECIFY
    205         feed.Outlet.F = 113.4 * 'kmol/h';
    206         feed.Outlet.T = 291 * 'K';
    207         feed.Outlet.P = 1.66 * 'atm';
    208         feed.Outlet.z = [0.5, 0.5];
     130        feed.Outlet.F = 0 * 'kmol/h';
     131        feed.Outlet.T = 300 * 'K';
     132        feed.Outlet.P = 1 * 'atm';
     133        feed.Outlet.z = 1/NComp;
    209134
    210         inL.F = 61.99 * 'kmol/h';       
    211         inL.P = 1.63 * 'atm';
    212         inL.T = 310 * 'K';
    213         inL.z = [0.1641, 0.8359];
     135        inL.F = 7121 * 'kmol/h';       
     136        inL.P = 2.22 * 'atm';
     137        inL.T = 297.6 * 'K';
     138        inL.z = [0.226, 0.425, 0.035, 0.025, 0.289];
    214139
    215         inV.F = 201.25 * 'kmol/h';
    216         inV.P = 1.48 * 'atm';
    217         inV.T = 321 * 'K';
    218         inV.z = [0.0584, 0.9416];
     140        inV.F = 175.3 * 'kmol/h';
     141        inV.P = 2.3062 * 'atm';
     142        inV.T = 308.3 * 'K';
     143        inV.z = [0.265, 0.233, 0.150, 0.014, 0.338];
     144       
     145        #t1.OutletV.F = 177.9 * 'kmol/h';
     146        #t1.Qsio = 0.61;
    219147
    220         t1.hs = 0.4 * 'm';#0.075 * 'm';
    221 #       t1.OutletV.F = 147.1 * 'kmol/h';
    222 #       t1.OutletL.F = 229.5 * 'kmol/h';
    223 #       t1.Qsio = 2.2;
    224        
    225148        SET
    226149        #Metal Pall Ring - nominal packing size 50 mm - Billet and Schultes, 1999.
    … …  
    232155        t1.e = 0.951;
    233156        t1.a = 112.6 * 'm^2/m^3';
     157        t1.hs = 0.4 * 'm';#0.075 * 'm';
    234158
    235159        INITIAL
    236         t1.OutletL.T = 310 *'K';
    237         t1.ML = 0.5 * 'mol';
    238         t1.OutletL.z(1) = 0.05;
     160        t1.OutletL.T = 290 *'K';
     161        t1.ML = 0.2 * 'kmol';
     162        t1.OutletL.z([1:4]) = [0.286, 0.203, 0.148, 0.006];
    239163       
    240164        OPTIONS
    241         InitialFile = "packedStage_BilletSchultes_Test.rlt";
     165        #NLASolver(File="nlasolver");
     166        #InitialFile = "packedStage_BilletSchultes_Test.rlt";
    242167        TimeStep = 2;   
    243168        TimeEnd = 50;
Note: See TracChangeset for help on using the changeset viewer.