Changeset 453 for branches/packed


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

More tests

Location:
branches/packed
Files:
4 edited

Legend:

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

    r450 r453  
    13061306        stage([top+topdown:topdown:bot]).OutletV to stage([top:topdown:bot-topdown]).InletV;
    13071307        stage([top:topdown:bot-topdown]).OutletL to stage([top+topdown:topdown:bot]).InletL;
    1308        
    1309        
    1310 end
    1311 
     1308end
     1309
     1310Model PackedDistillation_kettle_cond
     1311        PARAMETERS
     1312        outer PP as Plugin(Brief = "External Physical Properties", Type="PP");
     1313        outer NComp as Integer;
     1314        NStages as Integer(Brief="Number of stages", Default=2);
     1315        topdown as Integer(Brief="Trays counting (1=top-down, -1=bottom-up)", Default=1);
     1316        top as Integer(Brief="Number of top tray");
     1317        bot as Integer(Brief="Number of bottom tray");
     1318       
     1319        SET
     1320        top = (NStages-1)*(1-topdown)/2+1;
     1321        bot = NStages/top;
     1322       
     1323        VARIABLES
     1324        stage(NStages) as packedStage;
     1325        cond as condenser;
     1326        reb as reboiler;
     1327        sptop as splitter;
     1328        pump1 as pump;
     1329       
     1330        CONNECTIONS
     1331        #vapor
     1332        reb.OutletV to stage(bot).InletV;
     1333        stage([top+topdown:topdown:bot]).OutletV to stage([top:topdown:bot-topdown]).InletV;
     1334        stage(top).OutletV to cond.InletV;
     1335       
     1336        #liquid
     1337        cond.OutletL to sptop.Inlet;   
     1338        sptop.Outlet2 to pump1.Inlet;
     1339        pump1.Outlet to stage(top).InletL;
     1340        stage([top:topdown:bot-topdown]).OutletL to stage([top+topdown:topdown:bot]).InletL;
     1341        stage(bot).OutletL to reb.InletL;
     1342       
     1343        EQUATIONS
     1344        reb.OutletV.F = stage(bot).OutletV.F;
     1345end
  • branches/packed/eml/stage_separators/tray.mso

    r451 r453  
    286286        d as length (Brief="size of an element of packing");
    287287        #h as length (Brief="Height of packing");
    288         C as Real (Brief="Constant for resitance factor equation");
    289         Cp as Real (Brief="Constant for resitance at loading point factor equation");
     288        C as Real (Brief="Constant for resitance factor equation"); #Billet pp 80 (Cf+Cfl)/2
     289        Cp as Real (Brief="Constant for resitance at loading point factor equation"); # Billet pp94
    290290        Mw(NComp)       as molweight    (Brief = "Component Mol Weight");
    291291
     
    329329        uV * Across = OutletV.F * vV;
    330330       
    331         "Coefficient of Resistance"
    332         ksi * C^2 * (uL/uV * sqrt(rhoV/rhoL) * (niL/niV)^5.8)^3 = g/1*'s^2/m';
     331#       "Coefficient of Resistance"
     332#       ksi * C^2 * (uL/uV * sqrt(rhoV/rhoL) * (niL/niV)^5.8)^3 = g/1*'s^2/m';
     333        #ksil = Cp * (exp(uL*rhoL/a/niL/200)*(hL/hLs)^0.3) * (64/Rev+(1.8/Rev)) *
    333334        ksil = Cp * (exp(uL*rhoL/a/niL/200)*(hL/hLs)^0.3) * (64/Rev+(1.8/Rev)) *
    334                                 ((e-hL)/e);#^1.5;
    335        
    336         "Liquid holdup and Liquid flow"
    337         a^2 * niL * uL = hL^1 *(g*rhoL/3 - ksi*a*rhoV*uV^2/(4*hL*(e-hL)^2));
    338 
    339         "Pressure drop and Vapor flow"
    340         (InletV.P - OutletV.P)/hs = ksil *(a/2 + 2/ds)*(uV^2*rhoV/(e-hL)^3);
     335                                        ((e-hL)/e);
     336       
     337#       "Liquid holdup and Liquid flow"
     338#       a^2 * niL * uL = hL^1 *(g*rhoL/3 - ksi*a*rhoV*uV^2/(4*hL*(e-hL)^2));
     339
     340#       "Pressure drop and Vapor flow"
     341#       (InletV.P - OutletV.P)/hs = ksil *(a/2 + 2/ds)*(uV^2*rhoV/(e-hL)^3);
    341342       
    342343        "Vapor Reynolds number"
     
    344345
    345346        "Theoretical Liquid Holdup"
    346         hLs = (12*a^2*niL*uL/g/rhoL)^0.333;
     347        #hLs = (12*a^2*niL*uL/g/rhoL)^0.333;
     348        hLs = (12*a^2*niL*uL/g/rhoL);
    347349end
    348 
    349350
    350351#*-------------------------------------------------------------------
  • 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
  • branches/packed/sample/stage_separators/sample_tray.mso

    r451 r453  
    143143
    144144        t1.Emv = 1;
    145         t1.hs = 0.075 * 'm';
    146 #       t1.OutletV.F = 147.1 * 'kmol/h';
    147 #       t1.OutletL.F = 229.5 * 'kmol/h';
    148 #       t1.ksi = 0.805;
     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;
    149149       
    150150        SET
     
    152152        t1.Q = 0 * 'kW';
    153153        t1.Across = 0.8 * 'm^2';
    154         t1.V = 0.8 * 'm^2' * 0.075 * 'm';
     154        t1.V = 0.8 * 'm^2' * 0.4 * 'm';
    155155        t1.ds = 1.009 * 'm';
    156156        t1.d = 50 * 'mm';
     
    167167        OPTIONS
    168168        #InitialFile = "/home/paula/tray_Test.rlt";
     169        InitialFile = "/home/paula/packedStage_Test.rlt";
    169170        #DAESolver(File="dassl");
    170         TimeStep = 0.1;
    171         TimeEnd = 10;
     171        TimeStep = 2;   
     172        TimeEnd = 50;
    172173end
Note: See TracChangeset for help on using the changeset viewer.