Changeset 471


Ignore:
Timestamp:
Mar 5, 2008, 6:19:19 PM (15 years ago)
Author:
Rodolfo Rodrigues
Message:

Updated equil and gibbs models

Location:
trunk/eml/reactors
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • trunk/eml/reactors/equil.mso

    r426 r471  
    2222*
    2323*   Assumptions:
    24 *               * single- and two-phases involved
     24*               * single-phases involved
    2525*       * thermodynamic equilibrium
    2626*               * steady-state
    … …  
    9797#               + PP.IdealGasEnthalpyOfFormationAt25C()*(1 - Outlet.T/To);
    9898
     99        for j in [1:NReac]
    99100        "Gibbs free energy of reaction"
    100         sumt(G*stoic) = -Rg*Outlet.T*ln(K);
    101 #       K = exp(-sumt(G*stoic)/(Rg*Outlet.T));
    102        
    103         for j in [1:NReac]
     101                sumt(G*stoic(:,j)) = -Rg*Outlet.T*ln(K(j));
     102#               K(j) = exp(-sumt(G*stoic(:,j))/(Rg*Outlet.T));
     103
    104104        "Equilibrium constant"
    105105                K(j) = prod(activ^stoic(:,j));

  • trunk/eml/reactors/gibbs.mso

    r426 r471  
    2222*
    2323*   Assumptions:
    24 *               * single- and two-phases involved
     24*               * single-phases involved
    2525*       * thermodynamic equilibrium
    2626*               * steady-state
Note: See TracChangeset for help on using the changeset viewer.