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

some adaption in models

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

Legend:

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

    r755 r757  
    273273        ATTRIBUTES
    274274        Pallete         = true;
    275         Icon            = "icon/Condenser_column";
    276         Brief           = "Model of a  dynamic condenser with splitter.";
     275        Icon            = "icon/Condenser_column2";
     276        Brief           = "Model of a  dynamic condenser with control.";
    277277        Info            =
    278278"== Assumptions ==
     
    286286       
    287287== Initial Conditions ==
    288 * the condenser temperature (OutletLiquid.T);
    289 * the condenser liquid level (Level);
    290 * (NoComps - 1) OutletLiquid (OR OutletVapour) compositions.
     288* Initial_Temperature :  the condenser temperature (OutletLiquid.T);
     289* Initial_Level : the condenser liquid level (Level);
     290* Initial_Composition : (NoComps) OutletLiquid compositions.
    291291";     
    292292       
    293293PARAMETERS
    294294        outer PP                        as Plugin       (Brief = "External Physical Properties", Type="PP");
    295         outer NComp     as Integer;
     295        outer NComp     as Integer (Brief="Number of Components");
    296296
    297297        V                       as volume       (Brief="Condenser total volume");
    298         Across  as area                         (Brief="Cross Section Area of reboiler");
    299        
    300         Initial_Level                           as length                       (Brief="Initial Level of liquid phase");
    301         Initial_Temperature                     as temperature          (Brief="Initial Temperature of Condenser");
    302         Initial_Composition(NComp)      as fraction             (Brief="Initial Liquid Composition");
     298        Across  as area                 (Brief="Cross Section Area of condenser");
     299       
     300        Initial_Level                                                           as length                               (Brief="Initial Level of liquid phase");
     301        Initial_Temperature                                     as temperature (Brief="Initial Temperature of Condenser");
     302        Initial_Composition(NComp)      as fraction                     (Brief="Initial Liquid Composition");
    303303       
    304304VARIABLES
    305 in              InletVapour             as stream                               (Brief="Vapour inlet stream", PosX=0, PosY=0.5, Symbol="_{inV}");
     305
     306in              InletVapour             as stream                                       (Brief="Vapour inlet stream", PosX=0, PosY=0.5, Symbol="_{inV}");
    306307out     OutletLiquid    as liquid_stream                (Brief="Liquid outlet stream", PosX=0.5, PosY=1, Symbol="_{outL}");
    307 out     OutletVapour    as vapour_stream        (Brief="Vapour outlet stream", PosX=0.5, PosY=0, Symbol="_{outV}");
    308 in              InletQ          as power                                (Brief="Cold supplied", PosX=1, PosY=0.6, Symbol="_{in}");
    309 
    310         M(NComp)        as mol                          (Brief="Molar Holdup in the tray");
    311         ML                              as mol                          (Brief="Molar liquid holdup");
    312         MV                              as mol                          (Brief="Molar vapour holdup");
     308out     OutletVapour    as vapour_stream                (Brief="Vapour outlet stream", PosX=0.5, PosY=0, Symbol="_{outV}");
     309in              InletQ                          as power                                        (Brief="Heat supplied", Protected = true, PosX=1, PosY=0.6, Symbol="_{in}");
     310
     311        out     TCI as control_signal   (Brief="Temperature  Indicator of Condenser", Protected = true, PosX=1, PosY=0.40);
     312        out     LCI as control_signal   (Brief="Level  Indicator of Condenser", Protected = true, PosX=1, PosY=0.25);
     313
     314        M(NComp)        as mol                                  (Brief="Molar Holdup in the tray");
     315        ML                              as mol                                  (Brief="Molar liquid holdup");
     316        MV                              as mol                                  (Brief="Molar vapour holdup");
    313317        E                                       as energy                       (Brief="Total Energy Holdup on tray");
    314318        vL                              as volume_mol   (Brief="Liquid Molar Volume");
    315319        vV                              as volume_mol   (Brief="Vapour Molar volume");
    316         Level                   as length                       (Brief="Level of liquid phase");
     320        Level                   as length                               (Brief="Level of liquid phase");
    317321
    318322INITIAL
    319323
    320         Level                                   = Initial_Level;
    321         OutletLiquid.T                          = Initial_Temperature;
    322         OutletLiquid.z(1:NComp-1)       = Initial_Composition(1:NComp-1)/sum(Initial_Composition);
     324"Initial Level"
     325        Level                                                                   = Initial_Level;
     326
     327"Initial Temperature"
     328        OutletLiquid.T                                          = Initial_Temperature;
     329
     330"Initial Composition"
     331        OutletLiquid.z(1:NComp-1)   = Initial_Composition(1:NComp-1)/sum(Initial_Composition);
    323332       
    324333EQUATIONS
     
    361370        Level = ML*vL/Across;
    362371
     372"Temperature indicator"
     373        TCI * 'K' = OutletLiquid.T;
     374
     375"Level indicator"
     376        LCI*V = Level*Across;
     377       
    363378end
  • 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.