Sign up and NOW to receive the latest news, updates and technological advancements made for the Special Core Analysis & Enhanced Oil Recovery Industry

If you don't sign up, you won't know when the next breakthrough occurs!

Fick’s Law of Binary Diffusion

Fick’s Law of Binary Diffusion 2016-10-25T11:54:45+00:00

Fundamentals of Fluid Flow in Porous Media


Chapter 3

Fick’s Law of Binary Diffusion

The description of diffusion involves a mathematical model based on a fundamental hypothesis or ‘‘law’’.  Imagine two large bulbs connected by a long thin capillary (Figure 3‑1). Both of bulbs are at the same pressure and temperature but are filled with two different gases (N2 upper bulb and CO2 lower one). The concentration of the carbon dioxide in the initially N2 filled bulb is measured with time to find how fast these two gases will mix. The concentration of CO2 varies linearly with time (Figure 3‑1). So the amount of CO2 transferred could be determined from this graph at each time step. Carbon dioxide flux is defined as follows:
Carbon Dioxide Flux

Simple Diffusion Experiment
Figure 3-1: Simple Diffusion Experiment

We can assume that the flux is proportional to the gas concentration difference and we can recognize that increasing the capillary tube length will decrease the flux, so:
Carbon Dioxide Flux Equation

The new proportionality constant D is the diffusion coefficient. Its introduction implies a model for diffusion, the model often called Fick’s law. There is a similarity between Fick’s law and Ohm’s law for flux of electrons:
Electron Flux Equation

The diffusion coefficient D is proportional to the reciprocal of the resistivity. So the general form of Fick’s law is
Fick's Law (general form)

Where, j is the diffusion flux, and the minus appears because of the opposite directions of diffusion flux and concentration gradient.

From eq. (3‑2), we see that D has units ( L2 / T ). Since Fick’s law is derived for molecules moving in Brownian motion, D is a molecular diffusion coefficient, which is called Do to be specific. The intensity (energy and freedom of motion) of these Brownian motions controls the value of D. Thus, D depends on the phase (solid, liquid or gas), temperature, and molecule size.

It should be considered that the Fick’s Law could not be applied when the capillary is very thin or two gases react.

Parallel to Fourier’s law for heat conduction Fick’s second law is developed as;
Fick's Second Law for Heat Conduction

for one dimensional unsteady state diffusion, and for the constant area, A, it becomes the more known Fick’s second law equation:
Fick's Second Law Equation

This equation can be applied only for isotropic media and when the potential for diffusion is only given by concentration gradients.  The diffusion coefficient is also independent of concentration.

Example 3-1

Find the diffusion flux and concentration profile in a steady diffusion across a thin film.


The objective is to determine how much solvent moves across the film and how the solvent concentration changes within the film. On each side of the film is a well-mixed solution of one solvent. The solvent diffuses from the fixed higher concentration, located at x < 0 on the left-hand side of the film, into the fixed, less concentrated solution, located at x > l on the right-hand side. As the first step mass balance on a thin layer located at some arbitrary position x within the film is written:
mass balance on thin layer located at some arbitrary position x within the film

Diffusion Across a Thin Film
Figure 3-2: Diffusion Across a Thin Film

Because the process is in steady state, the accumulation is zero. The diffusion rate is the diffusion flux times the film’s area A. Thus
Diffusion rate

Divide this equation by the layer volume: A x Δx
Diffusion rate

Combining eq. (3‑6) with Fick’s law equation (3‑2) yields the following equation:
Fick's Law combined with another equation

Where C = solvent concentration in the layer. There are two boundary conditions for this differential equation:
Boundary conditions for differential equation

Analytical solution of eq. (3‑7), according to these two boundary conditions, will be the concentration profile:
Concentration profile

The resulted solution for the concentration profile shows that the profile is independent of the diffusion coefficient. Base on the Fick’s law the diffusion flux can be found by differentiation of concentration profile:
Fick's Law the diffusion flux found by differentiation of concentration


[1] “Diffusion, Mass Transfer in Fluid Systems”, E.L. Custler


If you have any questions at all, please feel free to ask PERM!  We are here to help the community.