Changeset 495 for branches/packed


Ignore:
Timestamp:
Apr 7, 2008, 3:42:56 PM (14 years ago)
Author:
Paula Bettio Staudt
Message:

More implementation

Location:
branches/packed
Files:
2 edited

Legend:

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

    r493 r495  
    6363#       QsiS as Real (Brief="Resistance coefficient at the loading point", Lower = 0);
    6464#       QsiFl as Real (Brief="Resistance coefficient at the flooding point", Lower = 0);
     65        deltaP as pressure;
    6566
    6667        SET
     
    136137        Rev*invK = dp*uV*rhoV / (miV*(1-e));
    137138       
    138                 if Rev < 1e-4 then
     139#*      if Rev < 1e-4 then
    139140                Qsil = 1;
    140141        else
     
    143144                                exp(13300/'m'^1.5/a^1.5*sqrt((uL^2*a)/g));
    144145        end
    145        
     146*#     
     147        deltaP = InletV.P - OutletV.P;
    146148        "Pressure drop and Vapor flow"
    147         (InletV.P - OutletV.P)/hs  = Qsil*a*uV^2*rhoV*invK / (2*(e-hl)^3);
     149        deltaP/hs  = Qsil*a*uV^2*rhoV*invK / (2*(e-hl)^3);
    148150
    149151        "Liquid holdup"
     
    198200       
    199201        VARIABLES
    200         deltaP as Real (Unit='atm/m');
     202        dP_h as Real (Unit='atm/m');
    201203       
    202204        DEVICES
     
    212214
    213215        EQUATIONS
    214         deltaP = (t1.InletV.P - t1.OutletV.P)/t1.hs;
     216        dP_h = t1.deltaP/t1.hs;
    215217       
    216218        SPECIFY
     
    230232        inV.z = [0.0584, 0.9416];#[0.5, 0.5];#
    231233       
    232         #t1.OutletV.P = 148 * 'kPa';
    233234        t1.OutletV.F = 147.1 * 'kmol/h';
    234         #t1.Qsil = 1;
     235        #t1.deltaP = 0.05 * 'atm';
     236        t1.Qsil = 1;
    235237       
    236238        SET
     
    308310        #t1.OutletV.P = 2.30409 * 'atm';
    309311        t1.OutletV.F = 165 * 'kmol/h';
     312        #t1.deltaP = 0.0026 * 'atm';
    310313        t1.Qsil = 10;
    311 #       t1.QsiS = 20;
    312 #       t1.QsiFl = 50;
    313314       
    314315        SET
  • branches/packed/sample/stage_separators/sample_column.mso

    r493 r495  
    183183       
    184184        OPTIONS
    185         TimeStep = 10;
     185        TimeStep = 1;
    186186        TimeEnd = 100;
    187187        #GuessFile="SectionColumn_Test_with8tray.rlt";
Note: See TracChangeset for help on using the changeset viewer.