Ignore:
Timestamp:
Jul 25, 2008, 11:18:50 AM (14 years ago)
Author:
Rafael de Pelegrini Soares
Message:

Updated the models to work with some language constraints

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/gui/eml/reactors/equil.mso

    r530 r574  
    9797#               + PP.IdealGasEnthalpyOfFormationAt25C()*(1 - Outlet.T/To);
    9898
    99         for j in [1:NReac]
     99        for j in [1:NReac] do
    100100        "Gibbs free energy of reaction"
    101101                sumt(G*stoic(:,j)) = -Rg*Outlet.T*ln(K(j));
     
    106106        end
    107107
    108         for i in [1:NComp]
     108        for i in [1:NComp] do
    109109        "Outlet molar fraction"
    110110                Outlet.F*Outlet.z(i) = (Inlet.F*Inlet.z(i) + sumt(stoic(i,:)*extent));
    111111        end     
    112112
    113         for i in [1:NComp]
     113        for i in [1:NComp] do
    114114          if (Outletm.z(i) > 1e-16) then
    115115            "Molar conversion"
     
    192192#       K = exp(-sumt(G*stoic)/(Rg*Outlet.T));
    193193       
    194         for j in [1:NReac]
     194        for j in [1:NReac] do
    195195        "Equilibrium constant"
    196196                K(j) = prod(activ^stoic(:,j));
    197197        end
    198198       
    199         for i in [1:NComp]
     199        for i in [1:NComp] do
    200200        "Outlet molar fraction"
    201201                Outlet.F*Outlet.z(i) = (Inlet.F*Inlet.z(i) + sumt(stoic(i,:)*extent));
    202202        end     
    203203
    204         for i in [1:NComp]
     204        for i in [1:NComp] do
    205205          if (Outletm.z(i) > 1e-16) then
    206206            "Molar conversion"
Note: See TracChangeset for help on using the changeset viewer.