Packed Bed Reactor with Three Chemical Species: Alternative Computational Methods
Packed Bed Reactor with Three Chemical Species: Alternative Computational Methods
Consider the reversible reaction . This reaction involves three chemical species, , , and , and takes place in a packed bed reactor (PBR) with axial mixing. Our objective is to compare alternative solutions of this problem using Mathematica's NDSolve and using Chebyshev's orthogonal collocation method.
A+BC
k
f
⇌
k
r
A
B
C
Introduce the dimensionless parameters , = (Péclet number), and =kτ (Damköhler number). Then the governing equations are:
τ=
L
u
P
e
uL
D
D
a
∂
C
A
∂t
1
P
e
2
∂
C
A
∂
2
x
∂
C
A
∂x
D
a
∂
C
B
∂t
1
P
e
2
∂
C
B
∂
2
x
∂
C
B
∂x
D
a
and =+ r,
∂
C
C
∂t
1
P
e
2
∂
C
C
∂
2
x
∂
C
C
∂x
D
a
where and =/.
r=
C
A
C
B
C
C
K
e
K
e
k
f
k
r
Initially, the reactor contains only species and such that (x,t=0)=0.6mol/l and (x,t=0)=0.4mol/l.
A
B
C
A
C
B
The feed to the reactor is also composed of only and such that (x=,t)=0.6mol/l and (x=,t)=0.4mol/l.
A
B
C
A

0
C
B

0
The Demonstration solves this system of partial differential equations for userset values of the Péclet and Damköhler numbers and the equilibrium constant . The Mathematica NDSolve results are shown as colored solid curves, while the Chebyshev orthogonal collocation values are shown as colored dots. As you can see, the two approaches are in excellent agreement.
K
e