Ignore:
Timestamp:
Nov 29, 2007, 10:05:51 AM (15 years ago)
Author:
Rodolfo Rodrigues
Message:

Added symbols

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/eml/reactors/equil.mso

    r415 r421  
    2222*
    2323*   Assumptions:
     24*               * single- and two-phases involved
    2425*       * thermodynamic equilibrium
    2526*               * steady-state
     
    4748        Brief           = "Model of a generic vapour-phase equilibrium CSTR";
    4849        Info            = "
    49 Requires the information of:
    50 * number of reactions
    51 * matrix of stoichiometric coefficients (components by reactions)
     50== Assumptions ==
     51* only vapour-phase
     52* thermodynamic equilibrium
     53* steady-state
     54
     55== Specify ==
     56* inlet stream
     57* stoichiometric matrix
     58* equilibrium temperature
    5259";
    5360
    5461        PARAMETERS
    5562        NReac           as Integer              (Brief="Number of reactions", Default=1);
    56     stoic(NComp,NReac) as Real  (Brief="Stoichiometric matrix");
     63    stoic(NComp,NReac) as Real  (Brief="Stoichiometric matrix", Symbol="\nu");
    5764        Rg                      as Real                 (Brief="Universal gas constant", Unit='J/mol/K', Default=8.314);
    5865        fs(NComp)       as pressure     (Brief="Fugacity in standard state", Default=1, DisplayUnit='atm');
     
    6067
    6168        VARIABLES
    62 out Outlet              as vapour_stream; # Outlet stream
     69out Outlet              as vapour_stream(Brief="Outlet stream", PosX=1, PosY=1, Symbol="_{out}");
    6370
    6471        G(NComp)        as energy_mol   (Brief="Gibbs free-energy of formation");
    6572        K(NReac)        as Real                 (Brief="Equillibrium constant",Default=1.5);
    66         activ(NComp)as Real             (Brief="Activity",Default=0.2);
     73        activ(NComp)as Real             (Brief="Activity", Symbol="\hat{a}", Default=0.2);
    6774       
    6875        rate(NComp) as reaction_mol (Brief="Overall component rate of reaction");
    69         extent(NReac) as flow_mol       (Brief="Extent of reaction");
    70         conv(NComp) as Real             (Brief="Fractional conversion of component", Default=0); # Lower=-1e3, Upper=1e3);
     76        extent(NReac) as flow_mol       (Brief="Extent of reaction", Symbol="\xi");
     77        conv(NComp) as Real             (Brief="Fractional conversion of component", Symbol="X", Default=0); # Lower=-1e3, Upper=1e3);
    7178       
    7279        EQUATIONS
    7380        "Outlet stream"
    74         Outlet.F*Outlet.z = Outletm.F*Outletm.z + rate*V;
     81        Outlet.F*Outlet.z = Outletm.F*Outletm.z + rate*Vr;
    7582       
    7683        "Mechanical equilibrium"
     
    8794
    8895#       "Gibbs free-energy of formation without Cp correction"
    89 #       G = PP.IdealGasGibbsOfFormationAt25C()*Outlet.T/To+PP.IdealGasEnthalpyOfFormationAt25C()*(1-Outlet.T/To);
     96#       G = PP.IdealGasGibbsOfFormationAt25C()*Outlet.T/To
     97#               + PP.IdealGasEnthalpyOfFormationAt25C()*(1 - Outlet.T/To);
    9098
    9199        "Gibbs free energy of reaction"
    92 #       sumt(G*stoic) = -Rg*Outlet.T*ln(K);
    93         K = exp(-sumt(G*stoic)/(Rg*Outlet.T));
     100        sumt(G*stoic) = -Rg*Outlet.T*ln(K);
     101#       K = exp(-sumt(G*stoic)/(Rg*Outlet.T));
    94102       
    95103        for j in [1:NReac]
     
    131139        Brief           = "Model of a generic liquid-phase equilibrium CSTR";
    132140        Info            = "
    133 Requires the information of:
    134 * number of reactions
    135 * matrix of stoichiometric coefficients (components by reactions)
     141== Assumptions ==
     142* only liquid-phase
     143* thermodynamic equilibrium
     144* steady-state
     145
     146== Specify ==
     147* inlet stream
     148* stoichiometric matrix
     149* equilibrium temperature
    136150";
    137151
    138152        PARAMETERS
    139153        NReac           as Integer              (Brief="Number of reactions", Default=1);
    140     stoic(NComp,NReac) as Real  (Brief="Stoichiometric matrix");
     154    stoic(NComp,NReac) as Real  (Brief="Stoichiometric matrix", Symbol="\nu");
    141155        Rg                      as Real                 (Brief="Universal gas constant", Unit='J/mol/K', Default=8.314);
    142156        Ps                      as pressure     (Brief="Standard pressure", Default=1, DisplayUnit='bar');
     
    144158
    145159        VARIABLES
    146 out Outlet              as liquid_stream; # Outlet stream
     160out Outlet              as liquid_stream(Brief="Outlet stream", PosX=1, PosY=1, Symbol="_{out}");
    147161
    148162        G(NReac)        as enth_mol     (Brief="Gibbs free-energy of formation");
    149163        K(NReac)        as fraction             (Brief="Equillibrium constant");
    150         activ(NComp)as Real             (Brief="Activity");
     164        activ(NComp)as Real             (Brief="Activity", Symbol="\hat{a}");
    151165       
    152166        rate(NComp) as reaction_mol (Brief="Overall component rate of reaction");
    153         extent(NReac)as flow_mol        (Brief="Extent of reaction");
    154         conv(NComp) as Real             (Brief="Fractional conversion of component", Default=0);
     167        extent(NReac)as flow_mol        (Brief="Extent of reaction", Symbol="\xi");
     168        conv(NComp) as Real             (Brief="Fractional conversion of component", Symbol="X", Default=0);
    155169       
    156170        EQUATIONS
    157171        "Outlet stream"
    158         Outlet.F*Outlet.z = Outletm.F*Outletm.z + rate*V;
     172        Outlet.F*Outlet.z = Outletm.F*Outletm.z + rate*Vr;
    159173       
    160174        "Mechanical equilibrium"
     
    171185
    172186#       "Gibbs free-energy of formation without Cp correction"
    173 #       G = PP.IdealGasGibbsOfFormationAt25C()*Outlet.T/To+PP.IdealGasEnthalpyOfFormationAt25C()*(1-Outlet.T/To);
     187#       G = PP.IdealGasGibbsOfFormationAt25C()*Outlet.T/To
     188#               + PP.IdealGasEnthalpyOfFormationAt25C()*(1 - Outlet.T/To);
    174189
    175190        "Gibbs free energy of reaction"
    176 #       sumt(G*stoic) = -Rg*Outlet.T*ln(K);
    177         K = exp(-sumt(G*stoic)/(Rg*Outlet.T));
     191        sumt(G*stoic) = -Rg*Outlet.T*ln(K);
     192#       K = exp(-sumt(G*stoic)/(Rg*Outlet.T));
    178193       
    179194        for j in [1:NReac]
Note: See TracChangeset for help on using the changeset viewer.