Ignore:
Timestamp:
Feb 12, 2008, 4:55:35 PM (15 years ago)
Author:
Paula Bettio Staudt
Message:

More flowsheet tests

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/packed/sample/stage_separators/sample_column.mso

    r453 r455  
    508508        PP      as Plugin(Brief="Physical Properties",
    509509                Type="PP",
    510                 Components = [ "isobutane", "benzene"],
     510                Components = [ "isobutene", "n-pentane", "propylene", "benzene", "isobutane"],
    511511                LiquidModel = "PR",
    512512                VapourModel = "PR"
     
    525525       
    526526        CONNECTIONS
    527         feed.Outlet to sec.stage(15).Inlet;
    528        
    529         zero to sec.stage([1:14]).Inlet;
    530         zero to sec.stage([16:24]).Inlet;
    531        
    532         reb to sec.stage(24).InletV;
     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;
    533533        cond to sec.stage(1).InletL;
     534
     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;
    534542       
    535543        SPECIFY
     
    537545        feed.Outlet.T = 291 * 'K';
    538546        feed.Outlet.P = 168.3 * 'kPa';
    539         feed.Outlet.z = [0.5, 0.5];
     547        feed.Outlet.z = 1/NComp;
    540548
    541549        zero.F = 0 * 'kmol/h';
    542550        zero.T = 300 * 'K';
    543551        zero.P = 1 * 'atm';
    544         zero.z = [0.5, 0.5];
     552        zero.z = 1/NComp;
    545553        zero.v = 0;
    546554        zero.h = 0 * 'J/mol';
    547555       
    548         cond.F = 68 * 'kmol/h';
    549         cond.P = 150 * 'kPa';
    550         cond.T = 281.75 * 'K';
    551         cond.z = [0.6664, 0.3336];
    552 
    553         reb.F = 153 * 'kmol/h';
    554         reb.P = 185 * 'kPa';
    555         reb.T = 328.12 * 'K';
    556         reb.z = [0.001848, 0.9982];
    557        
    558         sec.stage.Emv = 1;
    559 #       sec.stage(1).OutletV.F = 150 * 'kmol/h';
    560 
    561         SET
    562         sec.NTrays = 24;
    563        
    564 #       sec.stage.PressureDropModel = "Leva";
    565 #       sec.stage.PackingType = "random";
    566 
    567         sec.stage.V = 4 * 'ft^3';
     556        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;
    568571        sec.stage.Q = 0 * 'kW';
    569         sec.stage.Ap = 3.94 * 'ft^2';
     572        sec.stage.Ap = 1.5 * 'm^2';
     573        sec.stage.ds = 1.009 * 'm';
     574        sec.stage.Cpo = 0.763;
     575        sec.stage.e = 0.951;
     576        sec.stage.a = 112.6 * 'm^2/m^3';
    570577
    571578        INITIAL
    572         sec.stage.OutletL.T = [290:(330-290)/(sec.NTrays-1):330] *'K';
    573         sec.stage.Level = 0.5 * 'ft';
    574         sec.stage.OutletL.z(1) = 0.5;
     579        sec.stage.OutletL.T = [273:(360-273)/(sec.NStages-1):360] *'K';
     580        sec.stage.ML = 0.5 * 'kmol';
     581        sec.stage.OutletL.z([1:4]) = 1/NComp;
    575582       
    576583        OPTIONS
    577         InitialFile = "/home/paula/SectionColumn_Test_with8tray.rlt";
    578         TimeStep = 1;
    579         TimeEnd = 100;
    580         #GuessFile="SectionColumn_Test_with8tray.rlt";
    581         #Dynamic = false;
     584        DAESolver(File="dassl");
     585        TimeStep = 0.01;
     586        TimeEnd = 1;
    582587end
    583588
    584 FlowSheet Packed_kettle_cond_Test
     589FlowSheet Packed_kettle_cond_Test_1
    585590        PARAMETERS
    586591        PP      as Plugin(Brief="Physical Properties",
     
    630635        Qc.OutletQ.Q = -700 * 'cal/min';
    631636        col.pump1.dP = 0.1 * 'atm';
    632         col.stage.Emv = 1;
    633         col.stage.hs = 100 * 'mm';
    634        
    635         SET
    636         col.NStages = 10;
     637        col.stage.Qsio = 0.05;
     638       
     639        EQUATIONS
     640        col.reb.OutletV.F = col.K * sqrt(Qr.OutletQ.Q);
     641       
     642        SET
     643        col.NStages = 5;
    637644        col.cond.V = 1 * 'l';
    638645        col.cond.Across = 100 * 'cm^2';
     646        col.K = 0.2 * 'mol*min^0.5/kg^0.5/m';
    639647        col.reb.V = 2 * 'l';
    640648        col.reb.Across = 200 * 'cm^2';
    641         col.stage.V = 576.4 * 'cm^3';
    642649        col.stage.Q = 0 * 'kW';
    643650        col.stage.Ap = 56.74 * 'cm^2';
    644651        col.stage.ds = 85 * 'mm';
    645         col.stage.d = 20 * 'mm';
    646         col.stage.C = 2.18;
    647         col.stage.Cp = 1.329;
    648652        col.stage.e = 0.662;
    649653        col.stage.a = 185.4 * 'm^2/m^3';
    650 
     654        col.H = 1 * 'm';
     655        col.stage.Cpo = 0.763;
     656       
    651657        INITIAL
    652658        # condenser
    653659        col.cond.OutletL.T = (63.5+273.15) *'K';
    654660        col.cond.Level = 2 * 'cm';
    655         col.cond.OutletL.z(1) = 0.16;
     661        col.cond.OutletL.z(1) = 0.8; #0.16;
    656662
    657663        # reboiler
     
    662668        # column stages
    663669        col.stage.OutletL.T = [(63.5+273.15):((82+273.15)-(63.5+273.15))/(col.NStages-1):(82+273.15)] * 'K';
    664         col.stage.Level = 1 * 'cm';
    665         col.stage.OutletL.z(1) = 0.16;
     670        col.stage.Level = 2 * 'cm';
     671        #col.stage.ML = 0.2 * 'mol';
     672        col.stage.OutletL.z(1) = 0.5; #0.16;
    666673
    667674        OPTIONS
    668675        TimeStep = 0.1;
    669676        TimeEnd = 50;
    670         #GuessFile = "/home/paula/SectionColumn_Test_with8tray.rlt";
    671         #GuessFile="Distillation_kettle_cond_Test.rlt";
    672         #Dynamic = false;       
    673677end
Note: See TracChangeset for help on using the changeset viewer.