Modified Bessel functions provide the exact radial concentration profile for steady diffusion with a first-order reaction in a solid cylinder. For a cylindrical catalyst pellet or a radially modeled tubular reactor, the dimensionless reactant concentration is typically

[ \Theta(\rho)=\frac{C_A(r)}{C_{As}} =\frac{I_0(\phi \rho)}{I_0(\phi)}, \qquad \rho=\frac{r}{R}, ]

where (I_0) is the modified Bessel function of the first kind, (R) is the cylinder radius, and (\phi=R\sqrt{k/D_e}) is the cylindrical Thiele modulus. The corresponding reaction-diffusion solution reveals internal concentration gradients, diffusion limitations, and the effectiveness of the catalyst.

Core takeaway: In a solid cylinder, the first-order reaction-diffusion balance produces a modified Bessel equation. The physically acceptable solution uses (I_0) because it remains finite at the center; the resulting profile and effectiveness factor show how strongly intraparticle diffusion suppresses the observable reaction rate.

How the Cylindrical Reaction-Diffusion Model Is Formed

The steady-state shell balance

Consider a cylindrical pellet with radial coordinate (r), radius (R), effective diffusivity (D_e), and a first-order reaction

[ -r_A=kC_A. ]

At steady state, radial diffusion into the pellet must balance reactant consumption. Assuming constant (D_e), no axial variation, and no reaction-rate variation other than through (C_A), the governing equation is

[ D_e\left[\frac{1}{r}\frac{d}{dr} \left(r\frac{dC_A}{dr}\right)\right]-kC_A=0. ]

Equivalently,

[ \frac{d^2C_A}{dr^2} +\frac{1}{r}\frac{dC_A}{dr} -\frac{k}{D_e}C_A=0. ]

Why the equation is a modified Bessel equation

Introduce the dimensionless radius

[ \rho=\frac{r}{R} ]

and define the cylindrical Thiele modulus

[ \phi=R\sqrt{\frac{k}{D_e}}. ]

The equation becomes

[ \frac{d^2\Theta}{d\rho^2} +\frac{1}{\rho}\frac{d\Theta}{d\rho} -\phi^2\Theta=0, ]

where

[ \Theta=\frac{C_A}{C_{As}}. ]

This is the modified Bessel equation of order zero.

Applying the Boundary Conditions

Centerline symmetry

At the cylinder center, there is no preferred radial direction. Therefore,

[ \left.\frac{dC_A}{dr}\right|_{r=0}=0. ]

The concentration must also remain finite at (r=0).

Surface concentration

For a pellet whose external mass-transfer resistance is neglected, the surface concentration is prescribed:

[ C_A(R)=C_{As}, ]

or, in dimensionless form,

[ \Theta(1)=1. ]

If an external fluid-film resistance is important, the surface concentration is not known in advance and must instead be determined from a flux-matching condition.

Selecting the physically valid Bessel function

The general solution is

[ \Theta(\rho)=A I_0(\phi\rho)+B K_0(\phi\rho), ]

where (I_0) and (K_0) are modified Bessel functions of the first and second kinds.

The function (K_0) diverges as its argument approaches zero. It is therefore inadmissible for a solid cylinder containing the centerline, which requires a finite concentration. Setting (B=0) and applying the surface condition gives

[ \boxed{ \Theta(\rho)= \frac{I_0(\phi\rho)}{I_0(\phi)} } ]

and hence

[ \boxed{ C_A(r)=C_{As} \frac{I_0\left(r\sqrt{k/D_e}\right)} {I_0\left(R\sqrt{k/D_e}\right)} }. ]

The concentration is lowest at the center and highest at the surface because the reaction consumes reactant throughout the pellet.

What the Concentration Profile Reveals

Weak diffusion limitation

When (\phi\ll 1), diffusion is fast relative to reaction. The concentration remains nearly uniform:

[ C_A(r)\approx C_{As}. ]

The pellet operates close to its intrinsic kinetic rate, so increasing diffusivity has little additional benefit.

Strong diffusion limitation

When (\phi\gg 1), reaction is fast relative to radial diffusion. Reactant penetrates only a limited distance from the surface, producing a steep concentration gradient toward the center.

In this regime, much of the pellet interior contributes little to the overall reaction rate because its local reactant concentration is low.

Interpreting the Damköhler number

The relevant dimensionless group is often written as

[ \mathrm{Da}=\frac{kR^2}{D_e}=\phi^2. ]

It compares the characteristic reaction rate with the characteristic diffusion rate. Large (\mathrm{Da}) indicates stronger internal diffusion limitation, while small (\mathrm{Da}) indicates a nearly uniform concentration field.

The precise numerical value depends on the selected characteristic length and the definition used by the reactor model, so the geometry and convention should always be stated.

Calculating the Effectiveness Factor

Definition

The internal effectiveness factor compares the actual reaction rate integrated over the pellet with the rate that would occur if the entire pellet were at the surface concentration:

[ \eta= \frac{\text{actual total reaction rate}} {\text{rate at } C_A=C_{As}\text{ everywhere}}. ]

For a first-order reaction in a solid cylinder, integration of the Bessel-function profile gives

[ \boxed{ \eta= \frac{2I_1(\phi)} {\phi I_0(\phi)} } ]

where (I_1) is the modified Bessel function of the first kind of order one.

Limiting behavior

For (\phi\ll 1),

[ \eta\approx 1, ]

so nearly the entire pellet is kinetically active.

For (\phi\gg 1),

[ \eta\approx \frac{2}{\phi}, ]

which shows that the effective reaction rate decreases as the Thiele modulus increases.

Relating the factor to the measured rate

The total reaction rate per pellet volume can be written as

[ \overline{r}_A

\eta k C_{As}. ]

Thus, a pilot plant model can use (I_0) to calculate the internal profile and (\eta) to convert the intrinsic kinetic constant into an experimentally observable pellet-scale rate.

Using the Model in Pilot Plant Analysis

Estimating internal concentration gradients

Given (k), (D_e), (R), and (C_{As}), engineers can calculate (\phi), evaluate (I_0), and determine (C_A(r)) throughout the pellet.

This identifies whether the pellet center remains active or whether the reaction is concentrated near the outer region.

Comparing catalyst geometries

The radius enters the Thiele modulus through (R\sqrt{k/D_e}). Reducing the characteristic radial dimension lowers (\phi), improves reactant penetration, and generally increases the effectiveness factor.

This is why catalyst geometry, pellet size, and structured supports are important design variables in pilot-scale reactors.

Separating internal and external limitations

The Bessel-function solution describes internal diffusion and reaction. It does not, by itself, account for resistance in the surrounding fluid film.

With an external film, the model combines the internal effectiveness factor with a flux condition such as

[ N_A=k_g(C_{Ab}-C_{As})

-D_e\left.\frac{dC_A}{dr}\right|_{r=R}, ]

where (C_{Ab}) is the bulk-fluid concentration and (k_g) is the external mass-transfer coefficient.

This distinction prevents external film resistance from being incorrectly interpreted as poor intrapellet diffusivity.

Related Geometries and Bessel Functions

Hollow cylindrical reactors and annular pellets

If the domain is a hollow cylinder or annulus, the centerline is excluded. Both independent solutions may then be retained:

[ C_A(r)=A I_0\left(r\sqrt{k/D_e}\right) +B K_0\left(r\sqrt{k/D_e}\right). ]

The coefficients are determined from the conditions at the inner and outer radii.

Comparison with spherical pellets

For a sphere, the governing equation contains the spherical radial term

[ \frac{1}{r^2}\frac{d}{dr} \left(r^2\frac{dC_A}{dr}\right) -\frac{k}{D_e}C_A=0. ]

Its solution can be expressed using hyperbolic functions, giving

[ \Theta(\rho)= \frac{\sinh(\phi\rho)} {\rho\sinh(\phi)} ]

under a fixed surface-concentration condition. This is not the same profile as the cylindrical (I_0) solution because the geometry changes the diffusion area available at each radius.

Understanding the Trade-offs

A larger pellet is not always more productive

Increasing pellet radius can increase the amount of catalyst per particle, but it also increases (\phi) and makes internal diffusion more limiting.

The resulting rate per catalyst volume may be lower than expected from intrinsic kinetics alone.

Effective diffusivity is model-dependent

The value of (D_e) incorporates effects such as pore structure, tortuosity, porosity, and sometimes adsorption or multicomponent transport assumptions.

Using a molecular diffusivity in place of an effective porous-medium diffusivity can substantially understate the predicted concentration gradient.

Isothermal, constant-property assumptions may fail

The standard profile assumes constant (D_e), constant (k), and a uniform temperature. Strong heat effects, variable porosity, deactivation, nonlinear kinetics, or concentration-dependent diffusivity require a modified model.

A fixed surface concentration may be unrealistic

The normalized (I_0) profile assumes that (C_{As}) is known. In an operating reactor, the surface concentration may be reduced by an external film or by axial depletion in the bulk fluid.

The internal solution should therefore be coupled to reactor-scale balances when those effects are significant.

Making the Right Choice for Your Goal

Use the cylindrical modified-Bessel solution as the baseline model, then expand it only when the operating conditions require additional transport or reaction effects.

If your primary focus is estimating internal concentration profiles: Calculate (\phi=R\sqrt{k/D_e}) and use (\Theta(\rho)=I_0(\phi\rho)/I_0(\phi)) with finite-center and fixed-surface boundary conditions.
If your primary focus is predicting catalyst performance: Use (\eta=2I_1(\phi)/[\phi I_0(\phi)]) to convert intrinsic first-order kinetics into the pellet-scale reaction rate.
If your primary focus is diagnosing mass-transfer limitations: Compare the internal Thiele modulus with external film resistance rather than attributing every rate loss to intraparticle diffusion.
If your primary focus is analyzing hollow cylindrical equipment: Retain both (I_0) and (K_0), because the excluded centerline no longer eliminates the (K_0) solution.

The modified Bessel formulation turns the competition between radial diffusion and first-order reaction into a quantitative basis for interpreting and designing cylindrical pilot-plant systems.

Summary Table:

Key element	Practical insight
Governing model	Steady radial diffusion with a first-order reaction
Concentration profile	`Theta(rho) = I_0(phi rho) / I_0(phi)`
Thiele modulus	`phi = R sqrt(k / D_e)`; larger values indicate stronger diffusion limitation
Physical solution	`I_0` is retained because it remains finite at the cylinder center; `K_0` diverges there
Effectiveness factor	`eta = 2 I_1(phi) / [phi I_0(phi)]`
Design implications	Smaller pellets and higher effective diffusivity improve reactant penetration
Extended models	External film resistance, annular geometries, variable properties, and heat effects may require additional equations

Improve Pilot-Scale Reactor and Catalyst Performance

Need to evaluate internal diffusion, reaction rates, or scale-up risks in cylindrical systems? Contact LABPARK to discuss the right Educational and Vocational Unit Operations Pilot Plant for your application.

LABPARK supports universities, research institutes, and enterprises with pilot plants for chemical engineering, bioprocess and biotechnology, and environmental and water treatment. Our systems help you investigate mass-transfer limitations, validate kinetic models, compare catalyst geometries, and train engineers using practical operating data. Tell us your reactor type, catalyst or process, operating conditions, and testing goals, and our team can help identify a suitable pilot-scale solution.