Second-Order Reaction with Diffusion in a Liquid Film

Initializing live version
Download to Desktop

Requires a Wolfram Notebook System

Interact on desktop, mobile and cloud with the free Wolfram Player or other Wolfram Language products.

Gas absorption is often enhanced by a chemical reaction. For instance, acid gases ( and ) are usually eliminated from natural gas by absorption using ethanolamine () as a basic solvent.


Consider the absorption of species with a solvent containing a species such as a second-order irreversible chemical reaction, , that takes place in a liquid film. Only species is present in the gas phase since has a very low vapor pressure (i.e., is a high boiling component). Species is not present in the bulk liquid since all of reacts with component in the liquid film.

he steady state material balances within the film are given by and , where the binary diffusion approximation for and in has been used. These equations simply state that the rates of diffusion of species and are equal to the rate of the chemical reaction. The concentration of in the bulk liquid is arbitrarily set to 2 while the concentration of at the gas-liquid interface is set to 1.

This Demonstration displays the liquid film concentrations of species and (blue and orange curves, respectively) as a function of position. You can change the values of the diffusivities, and , the reaction rate constant, , as well as the number of Chebyshev collocation points, . Excellent agreement is obtained between the numerical solutions given by Chebyshev orthogonal collocation (blue and orange dots) and by NDSolve (blue and orange curves).


Contributed by: Housam Binousand Brian G. Higgins (July 2013)
Open content licensed under CC BY-NC-SA



In the discrete Chebyshev–Gauss–Lobatto case, the interior points are given by . These points are the extrema of the Chebyshev polynomials of the first kind, .

The Chebyshev derivative matrix at the quadrature points is an matrix given by

, , for , and for and ,

where for and .

The matrix is then used as follows: and , where is a vector formed by evaluating at , , and and are the approximations of and at the .


[1] P. Moin, Fundamentals of Engineering Numerical Analysis, Cambridge, UK: Cambridge University Press, 2001.

[2] L. N. Trefethen, Spectral Methods in MATLAB, Philadelphia: SIAM, 2000.

[3] M. B. Cutlip and M. Shacham, Problem Solving in Chemical Engineering with Numerical Methods, Upper Saddle River, NJ: Prentice Hall, 1999.

Feedback (field required)
Email (field required) Name
Occupation Organization
Note: Your message & contact information may be shared with the author of any specific Demonstration for which you give feedback.