Ignore:
Timestamp:
Feb 7, 2008, 5:26:42 PM (15 years ago)
Author:
Paula Bettio Staudt
Message:

More tests

File:
1 edited

Legend:

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

    r448 r453  
    582582end
    583583
     584FlowSheet Packed_kettle_cond_Test
     585        PARAMETERS
     586        PP      as Plugin(Brief="Physical Properties",
     587                Type="PP",
     588                Components = [ "methanol", "water"],
     589                LiquidModel = "IdealLiquid",
     590                VapourModel = "Ideal"
     591        );
     592        NComp   as Integer;
     593       
     594        VARIABLES
     595        Qc as energy_source (Brief="Heat rate removed from condenser");
     596        Qr as energy_source (Brief="Heat rate supplied to reboiler");
     597       
     598        SET
     599        NComp = PP.NumberOfComponents;
     600
     601        DEVICES
     602        col as PackedDistillation_kettle_cond;
     603        feed as source;
     604        zero as stream;
     605       
     606        CONNECTIONS
     607        feed.Outlet to col.reb.Inlet;
     608        zero to col.stage([1:col.NStages]).Inlet;
     609        Qc.OutletQ to col.cond.InletQ;
     610        Qr.OutletQ to col.reb.InletQ;
     611       
     612        SPECIFY
     613        feed.Outlet.F = 0.98 * 'mol/min';
     614        feed.Outlet.T = (50+273.15) * 'K';
     615        feed.Outlet.P = 0.91 * 'kPa';
     616        feed.Outlet.z = [0.16, 0.84];
     617       
     618        zero.F = 0 * 'kmol/h';
     619        zero.T = 300 * 'K';
     620        zero.P = 1 * 'atm';
     621        zero.z = 1/NComp;
     622        zero.v = 0;
     623        zero.h = 0 * 'J/mol';
     624       
     625        col.sptop.frac = 0.4;
     626        col.sptop.Outlet1.F = 0.10 * 'mol/min';
     627        col.reb.OutletL.F = 0.88 * 'mol/min';
     628        col.cond.OutletV.F = 0 * 'kmol/h';
     629        Qr.OutletQ.Q = 700 * 'cal/min';
     630        Qc.OutletQ.Q = -700 * 'cal/min';
     631        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.cond.V = 1 * 'l';
     638        col.cond.Across = 100 * 'cm^2';
     639        col.reb.V = 2 * 'l';
     640        col.reb.Across = 200 * 'cm^2';
     641        col.stage.V = 576.4 * 'cm^3';
     642        col.stage.Q = 0 * 'kW';
     643        col.stage.Ap = 56.74 * 'cm^2';
     644        col.stage.ds = 85 * 'mm';
     645        col.stage.d = 20 * 'mm';
     646        col.stage.C = 2.18;
     647        col.stage.Cp = 1.329;
     648        col.stage.e = 0.662;
     649        col.stage.a = 185.4 * 'm^2/m^3';
     650
     651        INITIAL
     652        # condenser
     653        col.cond.OutletL.T = (63.5+273.15) *'K';
     654        col.cond.Level = 2 * 'cm';
     655        col.cond.OutletL.z(1) = 0.16;
     656
     657        # reboiler
     658        col.reb.OutletL.T = (82+273.15) *'K';
     659        col.reb.Level = 2 * 'cm';
     660        col.reb.OutletL.z(1) = 0.16;
     661
     662        # column stages
     663        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;
     666
     667        OPTIONS
     668        TimeStep = 0.1;
     669        TimeEnd = 50;
     670        #GuessFile = "/home/paula/SectionColumn_Test_with8tray.rlt";
     671        #GuessFile="Distillation_kettle_cond_Test.rlt";
     672        #Dynamic = false;       
     673end
Note: See TracChangeset for help on using the changeset viewer.