Ignore:
Timestamp:
Jun 25, 2009, 6:05:45 PM (13 years ago)
Author:
gerson bicca
Message:

updated column/tank model

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/gui/eml/stage_separators/tank.mso

    r721 r781  
    5757        outer NComp     as Integer;
    5858        Across                          as area                         (Brief="Tank cross section area", Default=2);
     59        L                                       as length                       (Brief="Tank length");
     60        Vtotal                          as volume       (Brief="Tank total volume");
    5961       
    6062        Initial_Level                           as length               (Brief="Initial Level of the Tank");
    6163        Initial_Temperature                     as temperature  (Brief="Initial Temperature of Liquid");
    6264        Initial_Composition(NComp)      as fraction     (Brief="Initial Liquid Composition");
     65
     66SET
     67        Vtotal =Across*L;
    6368       
    6469VARIABLES
     
    6671out     Outlet          as liquid_stream        (Brief = "Outlet liquid stream", PosX=1, PosY=1, Symbol="_{out}");
    6772in      InletQ                  as power                                (Brief="Rate of heat supply", PosX=1, PosY=0.7859, Symbol="_{in}");
     73        Vfilled         as volume                               (Brief="Tank volume content");
    6874        Level                   as length                               (Brief="Tank level");
    6975        M(NComp)        as mol                                  (Brief="Molar Holdup in the tank");
     
    100106        Level = sum(M)*vL/Across;
    101107
     108"Volume Filled of liquid phase"
     109        Vfilled = Level*Across;
     110       
    102111end
    103112
     
    127136        outer PP                        as Plugin               (Brief = "External Physical Properties", Type="PP");
    128137        outer NComp     as Integer;
     138        pi                                      as positive     (Brief="Pi value", Default=3.141593,Hidden=true);
    129139        radius                          as length               (Brief="Tank radius");
    130140        L                                               as length               (Brief="Tank length");
     141        Vtotal                          as volume       (Brief="Tank total volume");
    131142       
    132143        Initial_Level                   as length                       (Brief="Initial Level of the Tank");
    133144        Initial_Temperature             as temperature                  (Brief="Initial Temperature of Liquid");
    134145        Initial_Composition(NComp)      as fraction                             (Brief="Initial Liquid Composition");
    135        
     146
     147SET
     148        Vtotal = (0.25*pi*(2*radius)^2)*L;
     149
    136150VARIABLES
    137151in      Inlet           as stream                               (Brief = "Inlet stream", PosX=0.1825, PosY=0, Symbol="_{in}");
     
    139153in      InletQ          as power                                (Brief="Rate of heat supply", PosX=1, PosY=0.6160, Symbol="_{in}");
    140154        Level           as length                               (Brief="Tank level");
     155        Vfilled         as volume                               (Brief="Tank volume content");
    141156        Across          as area                                         (Brief="Tank cross section area", Default=2);
    142157        M(NComp) as mol                                         (Brief="Molar Holdup in the tank");
     
    171186       
    172187"Cylindrical Area"
    173         Across = radius^2 * (asin(1) - asin((radius-Level)/radius) ) + (Level-radius)*sqrt(Level*(2*radius - Level));
     188        Across = radius^2 * (asin(1) - asin((radius-Level)/radius) ) + (Level-radius)*sqrt(Level*(2*radius - Level)+1E-8);
    174189
    175190"Level of liquid phase"
    176191        L*Across = sum(M)*vL;
     192
     193"Volume Filled of liquid phase"
     194        Vfilled = L*Across;
    177195
    178196end
Note: See TracChangeset for help on using the changeset viewer.