Fick's Laws of Diffusion

Whether you're designing a membrane separator, modeling drug delivery, or optimizing a chemical reactor, understanding how molecules move from regions of high concentration to low concentration is fundamental. This spontaneous process, known as molecular diffusion, is the primary mechanism for mass transfer in stagnant fluids or solids. The quantitative framework governing this movement was established by Adolf Fick in 1855, and his laws remain the cornerstone for analyzing countless engineering processes, from gas absorption to semiconductor doping.

The Driving Force: Concentration Gradient

At its heart, diffusion is driven by a concentration gradient. Imagine opening a bottle of perfume in a still room. The scent molecules, initially at a very high concentration near the bottle, will randomly collide and migrate into the surrounding air where their concentration is essentially zero. This net movement—from high to low concentration—continues until the molecules are uniformly distributed and the gradient disappears. The steeper this gradient, the faster the net movement. This intuitive concept is precisely captured by Fick's First Law.

Fick's First Law: The Steady-State Flux

Fick's First Law provides the equation for the steady-state rate of diffusion. It states that the molar flux $J_A$ of species A (the amount of A crossing a unit area per unit time) is proportional to its concentration gradient. In one-dimensional form (direction z), it is written as:

$$J_A=-D_{AB}\frac{d{C}_A}{dz}$$

Here, $J_A$ has typical units of $\text{mol/(m²cdotps)}$. $C_A$ is the molar concentration of A ($\text{mol/m³}$), and $d{C}_A/dz$ is its concentration gradient. The negative sign indicates that diffusion occurs down the concentration gradient (from high to low $C_A$). The proportionality constant $D_{AB}$ is the binary diffusion coefficient (or diffusivity) for species A diffusing through medium B. Its units are $\text{m²/s}$. This law is directly analogous to Fourier's law for heat conduction (flux proportional to temperature gradient) and Ohm's law for electricity (current proportional to voltage gradient).

Fick's Second Law: Transient Diffusion

Real-world processes are rarely at perfect steady state. Concentration changes with time at a given location, such as when a sugar cube dissolves in a cup of coffee. Fick's Second Law describes this time-dependent, or transient, diffusion. It is derived by performing a mass balance on a differential volume element, combining the First Law with the principle of conservation of mass. For diffusion in one dimension with a constant diffusivity $D_{AB}$, it takes the form:

$$\frac{\partial C_A}{\partial t}=D_{AB}\frac{{\partial }^2{C}_A}{\partial z^2}$$

This is a partial differential equation. The term $\partial C_A/\partial t$ represents the rate of accumulation (or depletion) of species A at a point. The equation shows that the rate of concentration change is proportional to the curvature of the concentration profile. Solutions to this equation, often involving error functions for semi-infinite media, are used to model case hardening of metals, drug release from pills, and transient diffusion in catalysts.

Important Diffusion Geometries and Flux Relations

Applying Fick's First Law requires integrating it across the diffusion path for specific boundary conditions. Two classic scenarios in chemical engineering are particularly important.

Equimolar Counter-diffusion occurs in a binary mixture where the molar flux of A in one direction is equal and opposite to the molar flux of B ($J_A=-J_B$). This is a good approximation for the distillation of an ideal binary mixture. For steady-state diffusion across a film of thickness $\delta z$, integration yields: $$J_A=\frac{D_{AB}}{\delta z}(C_{A1}-C_{A2})$$ where $C_{A1}$ and $C_{A2}$ are the concentrations at the boundaries. The flux is directly proportional to the concentration difference.

Diffusion Through a Stagnant (Non-diffusing) Film is even more common. Here, species A diffuses through a stationary layer of species B. A common example is a gas like water vapor diffusing through a stagnant layer of air to a condensing surface. Because B is stagnant, a bulk flow is induced in the direction of A's diffusion. This convective flow carries additional A, increasing the flux. The integrated form of Fick's Law for this case uses a log-mean concentration driving force: $$J_A=\frac{D_{AB}}{\delta z}\frac{C_T}{(C_B{)}_{lm}}(C_{A1}-C_{A2})$$ Here, $C_T$ is the total molar concentration, and $(C_B{)}_{lm}$ is the log-mean concentration of the stagnant gas B. The ratio $C_T/(C_B{)}_{lm}$ is always greater than 1, reflecting the enhancement of flux due to the induced bulk flow.

Estimating the Diffusion Coefficient

The binary diffusion coefficient $D_{AB}$ is a crucial property that depends on temperature, pressure, and the nature of the diffusing species and medium. For gases, kinetic theory provides a good foundation. A widely used empirical correlation is the Fuller-Schettler-Giddings equation, which estimates $D_{AB}$ as increasing with temperature ($T^{1.75}$) and decreasing with pressure. For liquids, diffusion is much slower (orders of magnitude smaller than in gases) and more complex due to solvent-solute interactions. The Wilke-Chang equation is a common estimation method, relating $D_{AB}$ to the solvent viscosity and molar volume. For solids, diffusivity varies enormously with structure (e.g., crystalline vs. polymeric) and is highly temperature-sensitive, typically following an Arrhenius-type relationship.

Common Pitfalls

1. Confusing Molar Flux with Mass Flux: Fick's Laws are most commonly expressed using molar units (concentration $C_A$, flux $J_A$). In systems with varying molecular weights or when using mass balances, you might need mass flux ($N_A$ in $\text{kg/(m²cdotps)}$). Always check your units and convert carefully using molecular weights.
2. Applying the Simplified Flux Equation Incorrectly: The simple equation $J_A=(D_{AB}/\delta z)\Delta C_A$ is only valid for equimolar counter-diffusion. The most frequent error is using it for diffusion through a stagnant film, where the log-mean correction factor must be included. Always identify which physical scenario you are modeling first.
3. Assuming Constant Diffusivity: While it simplifies the math, $D_{AB}$ can be a strong function of concentration, especially in liquids and polymers. For precise design work, especially involving large concentration changes, you may need to account for this variation, which makes solving Fick's Second Law more complex.
4. Neglecting the Analogy to Other Transport Processes: The mathematical forms of Fick's Law (mass), Fourier's Law (heat), and Newton's Law of viscosity (momentum) are identical. Failing to recognize this analogy means missing out on powerful problem-solving techniques. Solutions and dimensionless groups (like the Schmidt number $Sc=\nu /D_{AB}$, analogous to the Prandtl number for heat transfer) developed for one transport mode can often be directly adapted for another.

Summary

• Fick's First Law ($J_A=-D_{AB}(d{C}_A/dz)$) quantifies the steady-state diffusive flux, which is proportional to the concentration gradient. The constant of proportionality is the binary diffusion coefficient $D_{AB}$.
• Fick's Second Law ($\partial C_A/\partial t=D_{AB}({\partial }^2{C}_A/\partial z^2)$) governs transient diffusion, describing how concentration profiles evolve over time.
• Two key integrated solutions are for equimolar counter-diffusion (simple linear profile) and diffusion through a stagnant film (flux enhanced by bulk flow, requiring a log-mean concentration correction).
• The diffusion coefficient can be estimated using correlations like Fuller-Schettler-Giddings for gases and Wilke-Chang for liquids, and it is strongly dependent on temperature, pressure, and medium.
• The mathematical structure of Fick's Laws is directly analogous to the laws of heat and momentum transfer, enabling the use of similar solution techniques and dimensionless analysis across transport phenomena.