Ignore:
Timestamp:
Feb 25, 2010, 1:45:10 PM (14 years ago)
Author:
Rafael de Pelegrini Soares
Message:

Fixed some flash and cost models

Location:
branches/gui/eml/stage_separators
Files:
2 edited

Legend:

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

    r879 r918  
    1919using "tank";
    2020
    21 Model flash
     21Model flash as VesselVolume
    2222
    2323ATTRIBUTES
     
    6161
    6262        Mw(NComp)               as molweight    (Brief="Mol Weight", Hidden=true);
    63         Gconst          as acceleration (Brief="Gravity Acceleration",Default=9.81,Hidden=true);
    6463       
    6564        Levelpercent_Initial                    as positive     (Brief="Initial liquid height in Percent", Default = 0.70);
     
    7170        Mw=PP.MolecularWeight();
    7271
    73         Gconst  = 9.81 * 'm/(s^2)';
    74 
    7572VARIABLES
    76 
    77         Geometry                as VesselVolume (Brief="Vessel Geometry", Symbol=" ");
    7873
    7974in      Inlet                   as stream                       (Brief="Feed Stream", PosX=0, PosY=0.48, Symbol="_{in}");
     
    123118       
    124119"Energy Holdup"
    125         E = LiquidHoldup*OutletLiquid.h + VapourHoldup*OutletVapour.h - OutletLiquid.P*Geometry.Vtotal;
     120        E = LiquidHoldup*OutletLiquid.h + VapourHoldup*OutletVapour.h - OutletLiquid.P*Vtotal;
    126121       
    127122"Mol fraction normalisation"
     
    151146       
    152147"Static Head"   
    153         Pstatic = PP.LiquidDensity(OutletLiquid.T, Peq, OutletLiquid.z) * Gconst * Geometry.Level;
     148        Pstatic = PP.LiquidDensity(OutletLiquid.T, Peq, OutletLiquid.z) * Gconst * Level;
    154149
    155150"Mechanical Equilibrium for the Liquid Phase"
     
    163158
    164159"Geometry Constraint"
    165         Geometry.Vtotal = LiquidHoldup * vL + VapourHoldup * vV;
     160        Vtotal = LiquidHoldup * vL + VapourHoldup * vV;
    166161
    167162"Temperature indicator"
     
    172167
    173168"Level indicator"
    174         LI*Geometry.Vtotal= Geometry.Vfilled;
     169        LI*Vtotal= Vfilled;
    175170
    176171"Liquid Level"
    177         LiquidHoldup * vL = Geometry.Vfilled;
     172        LiquidHoldup * vL = Vfilled;
    178173
    179174end
  • branches/gui/eml/stage_separators/tank.mso

    r913 r918  
    3737       
    3838PARAMETERS
    39         PI                      as positive             (Brief="Pi value", Default=3.141593,Hidden=true, Symbol="\pi");
     39        Pi                      as positive             (Brief="Pi value", Default=3.141593,Hidden=true, Symbol="\pi");
    4040        Gconst                  as acceleration (Brief="Gravity Acceleration",Default=9.81,Hidden=true);
    4141
     
    5353
    5454        Gconst = 9.81 * 'm/(s^2)';
    55         Vhead_elliptical        = (PI*Diameter^3)/12;
    56         Vhead_hemispherical = (PI*Diameter^3)/6;
    57         Vcylinder = 0.25*(PI*Diameter^2)*Lenght;
     55        Vhead_elliptical        = (Pi*Diameter^3)/12;
     56        Vhead_hemispherical = (Pi*Diameter^3)/6;
     57        Vcylinder = 0.25*(Pi*Diameter^2)*Lenght;
    5858        radius = 0.5*Diameter;
    5959       
     
    7272
    7373"Vessel Cross Section Area"
    74         Across = 0.25*(PI*Diameter^2);
     74        Across = 0.25*(Pi*Diameter^2);
    7575
    7676switch Heads
     
    8484
    8585"Vessel Filled Volume"
    86         Vfilled = 0.25*PI*(((Diameter*Level)/(0.25*Diameter))^2)*(0.25*Diameter-Level/3);
     86        Vfilled = 0.25*Pi*(((Diameter*Level)/(0.25*Diameter))^2)*(0.25*Diameter-Level/3);
    8787
    8888else
    8989
    9090"Vessel Filled Volume"
    91         Vfilled = 0.25*PI*(Diameter^2)*(Level - 0.25*Diameter/3);
     91        Vfilled = 0.25*Pi*(Diameter^2)*(Level - 0.25*Diameter/3);
    9292
    9393end
     
    101101
    102102"Vessel Filled Volume"
    103         Vfilled = 0.25*PI*(Level^2)*(2*Diameter-4*Level/3);
     103        Vfilled = 0.25*Pi*(Level^2)*(2*Diameter-4*Level/3);
    104104
    105105else
    106106
    107107"Vessel Filled Volume"
    108         Vfilled = 0.25*PI*((2/3)*((0.5*Diameter)^3) - (0.25*(Diameter)^3) + Level*Diameter^2);
     108        Vfilled = 0.25*Pi*((2/3)*((0.5*Diameter)^3) - (0.25*(Diameter)^3) + Level*Diameter^2);
    109109
    110110end
Note: See TracChangeset for help on using the changeset viewer.