uz 2*u0*(1-(r/R)^2) Axial flow velocity xA (cA0-cA)/cA0 Conversion, species A cB cB0-cA0*xA Concentration, species B cC 2*cA0*xA Concentration, species C rA -A*exp(-E/(Rg*T))*(cA*cB-cC^2/Keq) Reaction rate Keq Keq0*exp(dHRx/Rg*(1/303[K]-1/T)) Equilibrium constant Q (-rA)*(-dHRx) Heat production