Changeset 574 for branches/gui/eml/reactors/equil.mso
- Timestamp:
- Jul 25, 2008, 11:18:50 AM (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/gui/eml/reactors/equil.mso
r530 r574 97 97 # + PP.IdealGasEnthalpyOfFormationAt25C()*(1 - Outlet.T/To); 98 98 99 for j in [1:NReac] 99 for j in [1:NReac] do 100 100 "Gibbs free energy of reaction" 101 101 sumt(G*stoic(:,j)) = -Rg*Outlet.T*ln(K(j)); … … 106 106 end 107 107 108 for i in [1:NComp] 108 for i in [1:NComp] do 109 109 "Outlet molar fraction" 110 110 Outlet.F*Outlet.z(i) = (Inlet.F*Inlet.z(i) + sumt(stoic(i,:)*extent)); 111 111 end 112 112 113 for i in [1:NComp] 113 for i in [1:NComp] do 114 114 if (Outletm.z(i) > 1e-16) then 115 115 "Molar conversion" … … 192 192 # K = exp(-sumt(G*stoic)/(Rg*Outlet.T)); 193 193 194 for j in [1:NReac] 194 for j in [1:NReac] do 195 195 "Equilibrium constant" 196 196 K(j) = prod(activ^stoic(:,j)); 197 197 end 198 198 199 for i in [1:NComp] 199 for i in [1:NComp] do 200 200 "Outlet molar fraction" 201 201 Outlet.F*Outlet.z(i) = (Inlet.F*Inlet.z(i) + sumt(stoic(i,:)*extent)); 202 202 end 203 203 204 for i in [1:NComp] 204 for i in [1:NComp] do 205 205 if (Outletm.z(i) > 1e-16) then 206 206 "Molar conversion"
Note: See TracChangeset
for help on using the changeset viewer.