Ignore:
Timestamp:
Jun 3, 2009, 5:07:22 PM (13 years ago)
Author:
gerson bicca
Message:

some adaption in models

File:
1 edited

Legend:

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

    r735 r757  
    337337       
    338338end
     339
     340Model reboiler_column
     341       
     342        ATTRIBUTES
     343        Pallete         = true;
     344        Icon            = "icon/reboiler_column";
     345        Brief           = "Model of a dynamic reboiler - kettle with control.";
     346        Info            =
     347"== Assumptions ==
     348
     349* perfect mixing of both phases;
     350* thermodynamics equilibrium;
     351* no liquid entrainment in the vapour stream.
     352       
     353== Specify ==
     354
     355* the inlet stream;
     356* the liquid inlet stream;
     357* the outlet flows: OutletVapour.F and OutletLiquid.F;
     358* the heat supply.
     359       
     360== Initial Conditions ==
     361
     362* Initial_Temperature :  the reboiler temperature (OutletLiquid.T);
     363* Initial_Level : the reboiler liquid level (Level);
     364* Initial_Composition : (NoComps) OutletLiquid compositions.
     365";     
     366       
     367PARAMETERS
     368        outer PP                        as Plugin               (Brief = "External Physical Properties", Type="PP");
     369        outer NComp     as Integer      (Brief="Number of Components");
     370        Across                          as area                 (Brief="Cross Section Area of reboiler");
     371        V                                               as volume       (Brief="Total volume of reboiler");
     372       
     373        Initial_Level                                                           as length                               (Brief="Initial Level of liquid phase");
     374        Initial_Temperature                                     as temperature  (Brief="Initial Temperature of Reboiler");
     375        Initial_Composition(NComp)      as fraction                     (Brief="Initial Liquid Composition");
     376
     377VARIABLES
     378
     379in      InletLiquid                     as stream                               (Brief="Liquid inlet stream", PosX=0, PosY=0.80, Symbol="_{inL}");
     380out     OutletLiquid    as liquid_stream                (Brief="Liquid outlet stream", PosX=0.50, PosY=1, Symbol="_{outL}");
     381out     OutletVapour    as vapour_stream        (Brief="Vapour outlet stream", PosX=0.50, PosY=0, Symbol="_{outV}");
     382in      InletQ                  as power                                (Brief="Heat supplied", Protected = true, PosX=1, PosY=0.55, Symbol="_{in}");
     383
     384        out     TCI as control_signal   (Brief="Temperature  Indicator of Reboiler", Protected = true, PosX=1, PosY=0.40);
     385        out     LCI as control_signal   (Brief="Level Indicator of Reboiler", Protected = true, PosX=1, PosY=0.25);
     386       
     387        M(NComp)        as mol                                  (Brief="Molar Holdup in the tray");
     388        ML                              as mol                                  (Brief="Molar liquid holdup");
     389        MV                              as mol                                  (Brief="Molar vapour holdup");
     390        E                                       as energy                       (Brief="Total Energy Holdup on tray");
     391        vL                              as volume_mol   (Brief="Liquid Molar Volume");
     392        vV                              as volume_mol   (Brief="Vapour Molar volume");
     393        Level                   as length                               (Brief="Level of liquid phase");
     394        rhoV                    as dens_mass            (Brief="Vapour Density");
     395
     396INITIAL
     397
     398        Level                                                                           = Initial_Level;
     399        OutletLiquid.T                                                  = Initial_Temperature;
     400        OutletLiquid.z(1:NComp-1)       = Initial_Composition(1:NComp-1)/sum(Initial_Composition);
     401
     402EQUATIONS
     403"Component Molar Balance"
     404        diff(M)= InletLiquid.F*InletLiquid.z    - OutletLiquid.F*OutletLiquid.z - OutletVapour.F*OutletVapour.z;
     405
     406"Energy Balance"
     407        diff(E) = InletLiquid.F*InletLiquid.h   - OutletLiquid.F*OutletLiquid.h - OutletVapour.F*OutletVapour.h + InletQ;
     408
     409"Molar Holdup"
     410        M = ML*OutletLiquid.z + MV*OutletVapour.z;
     411
     412"Energy Holdup"
     413        E = ML*OutletLiquid.h + MV*OutletVapour.h - OutletLiquid.P*V;
     414
     415"Mol fraction normalisation"
     416        sum(OutletLiquid.z)=1.0;
     417        sum(OutletLiquid.z)=sum(OutletVapour.z);
     418
     419"Vapour Density"
     420        rhoV = PP.VapourDensity(OutletVapour.T, OutletVapour.P, OutletVapour.z);
     421
     422"Liquid Volume"
     423        vL = PP.LiquidVolume(OutletLiquid.T, OutletLiquid.P, OutletLiquid.z);
     424
     425"Vapour Volume"
     426        vV = PP.VapourVolume(OutletVapour.T, OutletVapour.P, OutletVapour.z);
     427
     428"Chemical Equilibrium"
     429        PP.LiquidFugacityCoefficient(OutletLiquid.T, OutletLiquid.P, OutletLiquid.z)*OutletLiquid.z = PP.VapourFugacityCoefficient(OutletVapour.T, OutletVapour.P, OutletVapour.z)*OutletVapour.z;
     430
     431"Mechanical Equilibrium"
     432        OutletLiquid.P = OutletVapour.P;
     433
     434"Thermal Equilibrium"
     435        OutletLiquid.T = OutletVapour.T;
     436
     437"Geometry Constraint"
     438        V = ML*vL + MV*vV;
     439
     440"Level of liquid phase"
     441        Level = ML*vL/Across;
     442
     443end
Note: See TracChangeset for help on using the changeset viewer.