WOLFRAM|DEMONSTRATIONS PROJECT

Diffusion and Reaction in a Falling Liquid Film

​
% depth into liquid
20.
rate constant
0.5
c
A
at interface
1.
axial distance
50
diffusivity
0.0025
axial velocity
5
This Demonstration considers the adsorption of a gas that undergoes a first-order irreversible chemical reaction in a falling film of liquid in laminar flow.
The equations describing this system for the reaction
A
k
→
B
are:
v
max
1-
2
x
δ

∂
c
A
∂z
=
2
∂
c
A
∂
2
x
-k
c
A
,
v
max
1-
2
x
δ

∂
c
B
∂z
=
2
∂
c
B
∂
2
x
+k
c
A
,
with the following boundary conditions:
c
A
(x,0)=
c
B
(x,0)=0
,
c
A
(0,z)=
c
0
,
c
A
(0,z)=0
,
∂
c
A
(δ,z)
∂x
=
∂
c
B
(δ,z)
∂x
=0
,
where
c
A
and
c
B
are the concentrations of the adsorbed and product gases, respectively,
v
max
is the maximum fluid film velocity in the
z
direction,
x
represents the coordinate along the thickness of the fluid film,
δ
is the thickness of the film,

is the diffusion coefficient,
k
is the reaction constant and
c
0
is the solubility of
A
as determined by the partial pressure of
A
in the contiguous gas phase.
The first boundary condition corresponds to the fact that the film has no gases at the top
(z=0)
, the second indicates that at the liquid-gas interface only
A
is present at a concentration determined by its solubility and the third boundary condition states that the gases cannot diffuse through the solid wall.
The system of equations is solved using Mathematica's built-in NDSolve. The solution is shown for various values of the rate constant
k
, film depth
x
, axial distance
z
, diffusivity

, maximum axial velocity
v
max
and the solubility of gas
A
,
c
0
.