Diffusion through a Membrane

This Demonstration shows the effect of time delay on the dynamics of a simple diffusion model.
Consider a system of two compartments with lengths and containing the same species separated by a permeable membrane of unit area. Each compartment is assumed to be well stirred and homogeneous, so each can be represented by a single concentration variable, and .
Fick's first law of diffusion gives
, with ,
where represents time, is the diffusion coefficient of , and is the thickness of the membrane. Implicit in this equation is the assumption that the time required for diffusion through the membrane is zero; in reality, the time is not zero, in fact there is a distribution of times for molecules to cross the membrane. If we consider a model [1] in which all molecules take the same time to transverse the membrane, then Fick's law becomes
with initial history functions
This system of delayed differential equations is solved with . The solution shows remarkable features of the model. If the compartments have equal volumes and different initial concentrations, the system approaches equilibrium with damped oscillations. Further, if the volumes of the compartments are different, we can start with equal concentrations and the system will approach equilibrium in a damped oscillatory fashion. Even more surprising, we can start with equal concentrations and equal volumes and equilibrium will be approached with oscillations.


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


What is the meaning of the highly counterintuitive and apparently nonphysical behavior? The key lies in the fact that the final equilibrium value of lies significantly below the starting concentrations in the two compartments. The oscillations occur while the initially empty membrane is filling up. Consider the case of equal initial concentrations with cell much longer than cell . In the first interval , the same number of molecules leaves each cell, but, because cell is smaller, its concentration drops further. At time , when molecules begin to arrive, the effect is reversed, and the concentration in cell rises more rapidly than that in cell . Since molecules are now flowing in both directions, the effect will soon damp out. Readers may find the above result more convincing if they try to picture how the initially equal densities of people in two rooms of very different sizes connected by identical doors through an anteroom will change in time if, when the doors are opened, people move from room to room with equal probability and equal velocity without collisions. It is also reassuring to observe that when the results are averaged over a realistic distribution of (which results, among other things, from molecular collisions within the membrane), the oscillations disappear and the concentrations behave monotonically as expected.
[1] I. R. Epstein, "Delay Effects and Differential Delay Equations in Chemical Kinetics," International Reviews in Physical Chemistry, 11(1), 1992 pp. 135–160.
    • 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+