Introduction & Context

Molecular diffusivity quantifies how fast a solute migrates through a solvent due to random thermal motion. In process engineering it governs the rate of mass‑transfer limited steps such as gas absorption, liquid‑liquid extraction, crystallization, membrane separation, and heterogeneous catalysis. Accurate estimates allow engineers to size equipment (column heights, residence times, film thicknesses) and to interpret lab‑scale kinetic data, as well as to predict the diffusion layer thickness for design and scale‑up, and to complement these analyses with a thermal diffusivity calculation.

Methodology & Formulas

Convert temperature from Celsius to absolute scale \[T\,[\mathrm{K}]=T\,[^\circ\mathrm{C}]+273.15\]
Convert dynamic viscosity from centipoise to SI units \[\mu\,[\mathrm{Pa\,s}]=\mu\,[\mathrm{cP}]\times10^{-3}\]
Convert solute radius from nanometres to metres \[R\,[\mathrm{m}]=R\,[\mathrm{nm}]\times10^{-9}\]
Apply the Einstein–Stokes relation for the diffusion coefficient of a spherical particle in a continuum fluid \[D=\frac{k_{\mathrm{B}}T}{6\pi\mu R}\] where \(k_{\mathrm{B}}\) is the Boltzmann constant, \(1.380649\times10^{-23}\ \mathrm{J\,K^{-1}}\).

Continuum & creeping-flow regime limits
Parameter	Condition	Interpretation
Particle Reynolds number	\(\displaystyle Re_{\mathrm{p}}=\frac{2R\rho v}{\mu}\ll1\)	Inertial effects negligible; Stokes drag valid
Knudsen number	\(\displaystyle Kn=\frac{\lambda}{2R}\ll1\)	Fluid behaves as continuum; no slip at surface
Schmidt number	\(\displaystyle Sc=\frac{\mu}{\rho D}\gg1\)	Momentum diffusivity dominates mass diffusivity

When the above criteria are met, the Einstein–Stokes estimate provides a reliable order-of-magnitude value for design calculations and scale‑up analyses, and it can be further refined by applying steady‑state diffusion principles from Fick's law for more precise modeling.

The relation is reliable when all of the following are true:

Solute diameter > 5× solvent diameter (large spherical or quasi-spherical molecule)
Temperature is well below the solvent normal boiling point (viscosity is Newtonian)
System is dilute (< 1 mol % solute) so coupling effects are negligible
No specific interactions such as hydrogen bonding or micelle formation occur

Outside these limits use a semi-empirical correlation (Wilke-Chang, Hayduk-Minhas, etc.) or obtain experimental data.

Use the hydrodynamic (Stokes) radius, not the van-der-Waals or covalent radius:

Estimate from the molar volume at the normal boiling point: r = (3 V_b / 4πN_A)^1/3
If V_b is unknown, group-contribution methods (e.g., Joback) give ±10 % accuracy
For long-chain or disk-like molecules, take the radius of the equivalent sphere that has the same hydrodynamic drag; aspect-ratio corrections are available in literature

Use the solution viscosity at the composition of interest, not the pure-solvent value:

Measure it with a viscometer or retrieve from DIPPR/Dechema if available
If data are missing, blend pure-component viscosities with a mixing rule (e.g., Grunberg-Nissan) and validate against a single measurement point
Remember that even trace water or cosolvents can change μ by 20–30 %, so book values for “dry” solvent may mislead

For rigid hydrophobic solutes 0.5–2 nm in diameter, expect ±30 % of experimental D_AB; deviations rise to ±50–80 % for:

Small ions (hydration shell changes effective radius)
Flexible polymers (segmental motion not captured)
Supercritical solvents (density, not viscosity, controls friction)

Always benchmark against a lab measurement if the mass-transfer design margin is < 50 %.

Worked Example: Estimating Molecular Diffusivity for a Trace Contaminant in Water

A process engineer needs to estimate how fast a trace pharmaceutical contaminant (approximated as a spherical molecule with radius 0.45 nm) will diffuse through a 25 °C water stream. The stream viscosity is 0.89 cP. Use the Einstein–Stokes equation to obtain the molecular diffusivity.

Knowns

Temperature: 25.0 °C (298.15 K)
Dynamic viscosity of water: 0.89 cP = 0.00089 Pa·s
Molecular radius: 0.45 nm = 4.5 × 10⁻¹⁰ m
Boltzmann constant: 1.381 × 10⁻²³ J·K⁻¹

Step-by-Step Calculation

Convert temperature to kelvin: \( T = 25.0 + 273.15 = 298.15\ \text{K} \)
Convert viscosity to SI units: \( \mu = 0.89\ \text{cP} = 0.00089\ \text{Pa·s} \)
Convert radius to metres: \( r = 0.45\ \text{nm} = 4.5 \times 10^{-10}\ \text{m} \)
Compute the denominator of the Einstein–Stokes equation: \[ 6\pi\mu r = 6 \times 3.142 \times 0.00089\ \text{Pa·s} \times 4.5 \times 10^{-10}\ \text{m} = 7.55 \times 10^{-12}\ \text{kg·m·s}^{-2} \]
Calculate diffusivity: \[ D = \frac{kT}{6\pi\mu r} = \frac{(1.381 \times 10^{-23}\ \text{J·K}^{-1})(298.15\ \text{K})}{7.55 \times 10^{-12}\ \text{kg·m·s}^{-2}} = 5.45 \times 10^{-10}\ \text{m}^2\text{·s}^{-1} \]

Final Answer

The estimated molecular diffusivity of the contaminant in water at 25 °C is 5.45 × 10⁻¹⁰ m²·s⁻¹.

"Un projet n'est jamais trop grand s'il est bien conçu." — André Citroën

"La difficulté attire l'homme de caractère, car c'est en l'étreignant qu'il se réalise." — Charles de Gaulle

Molecular Diffusivity Estimation (Einstein-Stokes)

Introduction & Context

🚀 Skip the Manual Math!

Methodology & Formulas

Worked Example: Estimating Molecular Diffusivity for a Trace Contaminant in Water