# Reaction-Diffusion in a Two-Dimensional Catalyst Pellet

Reaction-Diffusion in a Two-Dimensional Catalyst Pellet

Consider the reaction-diffusion in a two-dimensional catalyst pellet with governing equations and boundary condition:

Le=++b(1-v)

u

t

u

xx

u

yy

2

ϕ

γu/(γ+u)

e

v

t

v

xx

v

yy

2

ϕ

γu/(γ+u)

e

u=v=0

x=y=±1

t=0

where is the Thiele modulus, is the adiabatic temperature rise (the Prater temperature), is the activation energy, and is the Lewis number.

2

ϕ

b

γ

Le

The steady-state temperature and reaction conversion along are plotted in magenta and blue, respectively. You can vary the parameters , , and as well as the number of Chebyshev collocation points, .

y=0

b

γ

2

ϕ

N+1