Ignore:
Timestamp:
Jul 15, 2009, 5:28:09 PM (13 years ago)
Author:
gerson bicca
Message:

updates

File:
1 edited

Legend:

Unmodified
Added
Removed

  • branches/gui/eml/stage_separators/condenser.mso

    r793 r794  
    2020using "streams";
    2121
    22 Model condenser2
    23         ATTRIBUTES
    24         Pallete         = true;
    25         Icon            = "icon/Condenser";
    26         Brief   = "Model of a dynamic condenser.";
    27         Info            =
    28 "== Assumptions ==
    29 * perfect mixing of both phases;
    30 * thermodynamics equilibrium.
    31        
    32 == Specify ==
    33 * the inlet stream;
    34 * the outlet flows: OutletVapour.F and OutletLiquid.F;
    35 * the InletQ (the model requires an energy stream).
    36        
    37 == Initial Conditions ==
    38 * Initial_Temperature :  the condenser temperature (OutletLiquid.T);
    39 * Initial_Level : the condenser liquid level (Level);
    40 * Initial_Composition : (NoComps) OutletLiquid compositions.
    41 ";     
    42        
    43 PARAMETERS
    44         outer PP                        as Plugin       (Brief = "External Physical Properties", Type="PP");
    45         outer NComp     as Integer(Brief = "Number of Components");
    46 
    47         V                       as volume       (Brief="Condenser total volume");
    48         Across  as area                         (Brief="Cross Section Area of condenser");
    49        
    50         Initial_Level                           as length                       (Brief="Initial Level of liquid phase");
    51         Initial_Temperature                     as temperature          (Brief="Initial Temperature of Condenser");
    52         Initial_Composition(NComp)      as fraction             (Brief="Initial Liquid Composition");
    53        
    54 VARIABLES
    55 in              InletVapour             as stream                               (Brief="Vapour inlet stream", PosX=0.15, PosY=0, Symbol="_{inV}");
    56 out     OutletLiquid    as liquid_stream                (Brief="Liquid outlet stream", PosX=0.4513, PosY=1, Symbol="_{outL}");
    57 out     OutletVapour    as vapour_stream        (Brief="Vapour outlet stream", PosX=0.4723, PosY=0, Symbol="_{outV}");
    58 in              InletQ          as power                                (Brief="Cold supplied", PosX=1, PosY=0, Symbol="_{in}",Protected=true);
    59 
    60         M(NComp)        as mol                          (Brief="Condenser Total Molar Holdup",Protected=true);
    61         ML                              as mol                          (Brief="Molar liquid holdup",Protected=true);
    62         MV                              as mol                          (Brief="Molar vapour holdup",Protected=true);
    63         E                                       as energy                       (Brief="Total Energy Holdup",Protected=true);
    64         vL                              as volume_mol   (Brief="Liquid Molar Volume",Protected=true);
    65         vV                              as volume_mol   (Brief="Vapour Molar volume",Protected=true);
    66         Level                   as length                       (Brief="Level of liquid phase",Protected=true);
    67 
    68 INITIAL
    69 
    70         Level                                   = Initial_Level;
    71         OutletLiquid.T                          = Initial_Temperature;
    72         OutletLiquid.z(1:NComp-1)       = Initial_Composition(1:NComp-1)/sum(Initial_Composition);
    73        
    74 EQUATIONS
    75 "Component Molar Balance"
    76         diff(M) = InletVapour.F*InletVapour.z - OutletLiquid.F*OutletLiquid.z- OutletVapour.F*OutletVapour.z;
    77 
    78 "Energy Balance"
    79         diff(E) = InletVapour.F*InletVapour.h - OutletLiquid.F*OutletLiquid.h- OutletVapour.F*OutletVapour.h + InletQ;
    80 
    81 "Molar Holdup"
    82         M = ML*OutletLiquid.z + MV*OutletVapour.z;
    83        
    84 "Energy Holdup"
    85         E = ML*OutletLiquid.h + MV*OutletVapour.h - OutletVapour.P*V;
    86        
    87 "Liquid Mol fraction normalisation"
    88         sum(OutletLiquid.z)=1.0;
    89 
    90 "Mol fraction constraint"
    91         sum(OutletLiquid.z)=sum(OutletVapour.z);
    92 
    93 "Liquid Volume"
    94         vL = PP.LiquidVolume(OutletLiquid.T, OutletLiquid.P, OutletLiquid.z);
    95        
    96 "Vapour Volume"
    97         vV = PP.VapourVolume(OutletVapour.T, OutletVapour.P, OutletVapour.z);
    98 
    99 "Chemical Equilibrium"
    100         PP.LiquidFugacityCoefficient(OutletLiquid.T, OutletLiquid.P, OutletLiquid.z)*OutletLiquid.z =
    101                 PP.VapourFugacityCoefficient(OutletVapour.T, OutletVapour.P, OutletVapour.z)*OutletVapour.z;
    102 
    103 "Thermal Equilibrium"
    104         OutletLiquid.T = OutletVapour.T;
    105 
    106 "Mechanical Equilibrium"
    107         OutletVapour.P = OutletLiquid.P;
    108 
    109 "Geometry Constraint"
    110         V = ML*vL + MV*vV;
    111 
    112 "Level of liquid phase"
    113         Level = ML*vL/Across;
    114 
    115 end
    116 
    11722Model condenserSteady
    11823
Note: See TracChangeset for help on using the changeset viewer.