Chemical Reaction and Diffusion in a Spherical Catalyst Pellet (the Nonisothermal Effectiveness Factor)

Consider a nonisothermal first-order reaction taking place in a spherical catalyst pellet. The governing equations, obtained from mole and enthalpy balances, and the boundary conditions, obtained from the symmetry condition and the known values at the surface, are the following:
and ,
where is the radius of the spherical pellet, is the effective thermal conductivity of the pellet, the reaction rate constant, the effective binary diffusivity of within the pellet, the heat of reaction, the concentration, and the temperature.
The rate constant is a function of temperature: .
In terms of dimensionless quantities, we have the following boundary value problem:
with and ,
where and .
Three dimensionless parameters are introduced in the new governing equation:
, and , where , the Thiele modulus, is a measure of internal mass transfer resistance. When , mass transfer resistance is negligible and diffusion is very fast. Large values of , the dimensionless activation energy, mean that the reaction rate is very sensitive to temperature. In this Demonstration is taken equal to 1.
The Demonstration displays the effectiveness factor, , a measure of the importance of the total reaction rate inside the pellet compared to its value at the pellet surface, versus the Thiele modulus , for various values of parameter .
The effectiveness factor is given by the formula:
where .
The parameter , which you can set, is a measure of the relative importance of the heat of reaction relative to conduction.
Large values of mean that there is significant internal heating; thus for the effectiveness factor will be greater than 1. This means that the rate of reaction is higher inside the pellet compared to its value at the surface of the pellet in the absence of internal diffusion resistance.
Negative values of the parameter correspond to endothermic reactions () and smaller values of the effectiveness factor, .
When , the internal transport resistance becomes large (a depleted zone inside the pellet forms and the concentration of species is low). Thus, the reaction rate inside the pellet is smaller than its value at the pellet surface and consequently the effectiveness factor becomes very small compared to 1.


  • [Snapshot]
  • [Snapshot]
  • [Snapshot]
  • [Snapshot]


K. J. Beers, Numerical Methods for Chemical Engineering, Cambridge: Cambridge University Press, 2007.
    • Share:

Embed Interactive Demonstration New!

Just copy and paste this snippet of JavaScript code into your website or blog to put the live Demonstration on your site. More details »

Files require Wolfram CDF Player or Mathematica.

Mathematica »
The #1 tool for creating Demonstrations
and anything technical.
Wolfram|Alpha »
Explore anything with the first
computational knowledge engine.
MathWorld »
The web's most extensive
mathematics resource.
Course Assistant Apps »
An app for every course—
right in the palm of your hand.
Wolfram Blog »
Read our views on math,
science, and technology.
Computable Document Format »
The format that makes Demonstrations
(and any information) easy to share and
interact with.
STEM Initiative »
Programs & resources for
educators, schools & students.
Computerbasedmath.org »
Join the initiative for modernizing
math education.
Step-by-Step Solutions »
Walk through homework problems one step at a time, with hints to help along the way.
Wolfram Problem Generator »
Unlimited random practice problems and answers with built-in step-by-step solutions. Practice online or make a printable study sheet.
Wolfram Language »
Knowledge-based programming for everyone.
Powered by Wolfram Mathematica © 2018 Wolfram Demonstrations Project & Contributors  |  Terms of Use  |  Privacy Policy  |  RSS Give us your feedback
Note: To run this Demonstration you need Mathematica 7+ or the free Mathematica Player 7EX
Download or upgrade to Mathematica Player 7EX
I already have Mathematica Player or Mathematica 7+