Ignore:
Timestamp:
Dec 10, 2007, 11:24:49 AM (15 years ago)
Author:
Rodolfo Rodrigues
Message:

Updated to use vol_tank.mso

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/eml/reactors/tank_basic.mso

    r422 r426  
    3030
    3131using "streams";
     32using "vol_tank";
    3233
    3334
     
    4344in      Inlet   as stream       (Brief="Inlet stream", PosX=0, PosY=0, Symbol="_{in}");
    4445        Outletm as stream       (Brief="Intermediary outlet stream", Symbol="_{outm}");
    45 
     46        Tank    as vol_tank     (Brief="Routine to volume tank calculation", Symbol="_{tank}");
     47       
    4648        M(NComp)as mol          (Brief="Component molar holdup");
    4749        Mt              as mol          (Brief="Total component molar holdup");
    48         Vr              as volume       (Brief="Volume of reactional mixture");
    4950        E               as energy       (Brief="Internal energy");
    5051        Q               as heat_rate(Brief="Reactor duty", Default=0);
    51        
    52         Across  as area         (Brief="Tank cross section area");
    53         Level   as length       (Brief="Tank level");
    5452       
    5553        EQUATIONS
     
    6563        "Energy balance"
    6664        diff(E) = Inlet.F*Inlet.h - Outletm.F*Outletm.h + Q;
    67        
    68         "Geometry"
    69         Vr = Across*Level;
    7065end
    7166
     
    8681
    8782        "Volume constraint"
    88         Vr = Mt*PP.VapourVolume(Outletm.T,Outletm.P,Outletm.z);
     83        Tank.V = Mt*PP.VapourVolume(Outletm.T,Outletm.P,Outletm.z);
    8984
    9085        "Total internal energy"
     
    108103       
    109104        "Volume constraint"
    110         Vr = Mt*PP.LiquidVolume(Outletm.T,Outletm.P,Outletm.z);
     105        Tank.V = Mt*PP.LiquidVolume(Outletm.T,Outletm.P,Outletm.z);
    111106       
    112107        "Total internal energy"
    113         E = Mt*Outletm.h - Outletm.P*Vr;
     108        E = Mt*Outletm.h - Outletm.P*Tank.V;
    114109end
    115110
     
    130125        OutletmL as liquid_stream       (Brief="Intermediary liquid outlet stream", Symbol="_{outmL}");
    131126out     OutletV  as vapour_stream       (Brief="Outlet vapour stream", Symbol="_{outV}");
    132 
     127        Tank     as vol_tank            (Brief="Routine to volume tank calculation", Symbol="_{tank}");
     128       
    133129        M(NComp)as mol                  (Brief="Component molar holdup");
    134130        ML              as mol                  (Brief="Molar liquid holdup");
    135131        MV              as mol                  (Brief="Molar vapour holdup");
    136         Vr              as volume               (Brief="Volume of reactional mixture");
    137132        E               as energy               (Brief="Internal energy");
    138133        Q               as heat_rate    (Brief="Reactor duty", Default=0);
    139134        vL              as volume_mol   (Brief="Liquid Molar Volume");
    140        
    141         Across  as area                 (Brief="Tank cross section area");
    142         Level   as length               (Brief="Tank level");
    143135       
    144136        EQUATIONS
     
    168160       
    169161        "Total internal energy"
    170         E = ML*OutletmL.h + MV*OutletV.h; #- OutletmL.P*V; P_tank*V_tank
     162        E = ML*OutletmL.h + MV*OutletV.h;
    171163       
    172164        "Geometry constraint"
    173         Vr = ML*vL + MV*PP.VapourVolume(OutletV.T,OutletV.P,OutletV.z);
     165        Tank.V = ML*vL + MV*PP.VapourVolume(OutletV.T,OutletV.P,OutletV.z);
    174166       
    175167       
     
    188180       
    189181        "Tank Level"
    190         ML*vL = Across*Level;
     182        ML*vL = Tank.V;
    191183end
Note: See TracChangeset for help on using the changeset viewer.