Solving Stefan-Maxwell Equations Using Orthogonal Collocation and Shooting Method

Consider acetone and methanol diffusing concurrently through a stagnant air column. The mole fractions of all three components are obtained by solving the species balance equation augmented with the Stefan–Maxwell equations for the diffusion fluxes. The resulting nonlinear boundary value problem can be solved using either a shooting method or by implementing an orthogonal collocation technique. In this Demonstration we use both methods. The plot shows the steady state concentration profiles of the three species in the diffusion column. The solid curves were obtained with a shooting technique using the terminal concentrations , , , and , using Mathematica's built-in function NDSolve. The blue, magenta, and brown colors correspond to acetone, methanol, and air, respectively. The colored dots shows the solution obtained using the orthogonal collocation technique for user-set values for the number of interior collocation points. As expected, excellent agreement is seen between both numerical methods.
To quantify the accuracy of the collocation method, we compute the Euclidean norm of the error for collocation points: where and are the mole fraction of a species with obtained with orthogonal collocation and NDSolve, respectively. This error is displayed for any value of . As the number of collocation points is increased, the Euclidean norm is reduced. For , .


  • [Snapshot]
  • [Snapshot]
  • [Snapshot]
  • [Snapshot]


The Stefan–Maxwell relationships for unidirectional diffusion are given by:
for , where is the flux of species (expressed in ). From the calculation we find that and , while is assumed equal to zero. Binary diffusivities are equal to , , and (all expressed in ).
[1] J. Benitez, Principles and Modern Applications of Mass Transfer Operations, New York: Wiley–Interscience, 2002.
[2] R. Carty and T. Schrodt, "Concentration Profiles in Ternary Gaseous Diffusion," Industrial Engineering and Chemistry Fundamentals, 14(3), 1975 pp. 276–278.
    • Share:

Embed Interactive Demonstration New!

Just copy and paste this snippet of JavaScript code into your website or blog to put the live Demonstration on your site. More details »

Files require Wolfram CDF Player or Mathematica.

Mathematica »
The #1 tool for creating Demonstrations
and anything technical.
Wolfram|Alpha »
Explore anything with the first
computational knowledge engine.
MathWorld »
The web's most extensive
mathematics resource.
Course Assistant Apps »
An app for every course—
right in the palm of your hand.
Wolfram Blog »
Read our views on math,
science, and technology.
Computable Document Format »
The format that makes Demonstrations
(and any information) easy to share and
interact with.
STEM Initiative »
Programs & resources for
educators, schools & students.
Computerbasedmath.org »
Join the initiative for modernizing
math education.
Step-by-Step Solutions »
Walk through homework problems one step at a time, with hints to help along the way.
Wolfram Problem Generator »
Unlimited random practice problems and answers with built-in step-by-step solutions. Practice online or make a printable study sheet.
Wolfram Language »
Knowledge-based programming for everyone.
Powered by Wolfram Mathematica © 2018 Wolfram Demonstrations Project & Contributors  |  Terms of Use  |  Privacy Policy  |  RSS Give us your feedback
Note: To run this Demonstration you need Mathematica 7+ or the free Mathematica Player 7EX
Download or upgrade to Mathematica Player 7EX
I already have Mathematica Player or Mathematica 7+