Changeset 918 for branches/gui/eml


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

Fixed some flash and cost models

Location:
branches/gui/eml
Files:
4 edited

Legend:

Unmodified
Added
Removed
  • branches/gui/eml/costs/flash_cost.mso

    r732 r918  
    6060        Ts                              as length (Brief="Thickness");
    6161        dens_mass_material as dens_mass (Brief = "Mass Density of the Material");
    62         Pi                              as constant (Brief="Pi Number",Default=3.14159265);
    6362
    6463        VARIABLES
  • branches/gui/eml/costs/reboiler_cost.mso

    r331 r918  
    6767        EQUATIONS
    6868        "Average pressure"
    69         P = 0.5*(InletL.P + OutletL.P);
     69        P = 0.5*(InletLiquid.P + OutletLiquid.P);
    7070       
    7171        "Capital Cost"
  • 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.