ToolFusion Chemical
ToolFusion Chemical
Calculator D3

Fick’s Laws of Diffusion in Gases and Liquids

Fick’s Laws describe how particles like salt or oxygen spread out naturally from areas where they’re crowded to areas where they’re sparse — like ink spreading in water.

Industry Applications
Chemical synthesis, pharmaceutical dissolution, wastewater aeration, semiconductor doping, food packaging barrier design
Key Standards
ISO 11357-3 (polymer diffusion), ASTM E3039 (membrane permeability), API RP 14E (mass transfer in offshore systems)
Typical Scale
Microscale (µm): drug delivery; Process scale (m): distillation columns; System scale (10s m): flue-gas scrubbers

⚠️ Why It Matters

1
Inaccurate diffusion coefficient estimation
↓
2
Underpredicted solute breakthrough in solvent extraction
↓
3
Excessive column length or residence time
↓
4
Reduced separation efficiency and yield
↓
5
Higher energy and capital cost for purification units

📘 Definition

Fick’s First Law quantifies the steady-state diffusive flux as proportional to the concentration gradient, with proportionality given by the diffusion coefficient. Fick’s Second Law describes how concentration changes over time due to non-uniform diffusion, expressed as a partial differential equation governing transient mass transport. Together, they form the foundational phenomenological framework for molecular diffusion in isotropic, non-reactive, single-phase media.

🎨 Concept Diagram

Fick’s First Law:J = −D dc/dxFick’s Second Law:∂c/∂t = D ∂²c/∂x²

AI-generated illustration for visual understanding

💡 Engineering Insight

Fick’s Laws assume no bulk motion, no chemical reaction, and isotropic media — yet real equipment violates all three. Always assess the Damköhler number (Da = reaction rate / diffusion rate) and Péclet number (Pe = convective flux / diffusive flux); if Da > 0.1 or Pe > 10, Fick’s First Law alone is insufficient — you must couple it with momentum and reaction terms in your model.

📖 Detailed Explanation

Fick’s First Law states that the diffusive flux J (mol·m⁻²·s⁻¹) equals −D∇c — a linear response to concentration slope, analogous to Fourier’s Law for heat. This holds only under steady-state, dilute, non-interacting conditions and forms the basis for designing simple dialysis membranes or gas-absorption trays.

Fick’s Second Law (∂c/∂t = D∇²c) emerges from mass conservation applied to the First Law. It predicts how concentration evolves over time — essential for modeling batch leaching, drug release from polymers, or oxygen ingress into packaged food. Analytical solutions exist only for idealized geometries (slab, cylinder, sphere), while industrial cases require numerical methods (e.g., finite-volume discretization in ANSYS Fluent or COMSOL).

At engineering scales, molecular diffusion rarely acts alone. In turbulent flows, eddy diffusivity dominates; in porous catalysts, Knudsen diffusion or surface diffusion may exceed bulk diffusion; in electrolytes, migration couples with diffusion. Advanced treatments replace D with a tensor (for anisotropic media), incorporate activity gradients instead of concentration, or embed Fickian terms within porous-media formulations like the dusty-gas model for fuel-cell electrodes.

🔄 Engineering Workflow

Step 1
Step 1: Identify diffusing species and phase (gas/liquid/solid matrix)
→
Step 2
Step 2: Determine thermodynamic state (T, P, composition) and estimate D using validated method (e.g., Chapman-Enskog or Wilke-Chang)
→
Step 3
Step 3: Characterize geometry and flow regime to compute Sc, Re, and Pe numbers
→
Step 4
Step 4: Select appropriate model (Fickian, effective-medium, or coupled convection-diffusion) and discretize domain
→
Step 5
Step 5: Solve transient or steady-state mass balance with boundary conditions (e.g., constant surface concentration or flux)
→
Step 6
Step 6: Validate against experimental concentration profiles or breakthrough curves
→
Step 7
Step 7: Iterate design parameters (e.g., column height, particle size, flow rate) to meet purity/yield targets

📋 Decision Guide

Rock/Field Condition Recommended Design Action
High Schmidt Number liquid (Sc > 500, e.g., aqueous organics) Use packed beds or static mixers to enhance interfacial area; avoid plug-flow assumptions without axial dispersion correction.
Low-pressure gas mixture (P < 100 kPa) with large molecular weight difference Apply Chapman-Enskog theory for D calculation; avoid empirical correlations like Fuller-Schettler-Giddings which overpredict by up to 40%.
Diffusion-limited reaction in porous catalyst (Thiele modulus ϕ > 3) Reduce pellet diameter or increase temperature to lower ϕ; verify D_eff via pulse-response experiments, not literature values alone.

📊 Key Properties & Parameters

Diffusion Coefficient (D)

1e−9 to 1e−5 m²/s (gases: ~1e−5; liquids: ~1e−9–1e−10; supercritical CO₂: ~1e−7)

A material-specific constant quantifying the rate at which molecules diffuse through a medium under a unit concentration gradient.

⚡ Engineering Impact:

Directly governs design scale of absorbers, extractors, and membrane contactors — errors >20% cause >30% oversizing.

Concentration Gradient (∇c)

0.1–100 mol/m⁴ (e.g., 5 mol/m³ across 0.01 m yields ∇c = 500 mol/m⁴)

The spatial rate of change of species concentration, driving the diffusive flux per Fick’s First Law.

⚡ Engineering Impact:

Determines local mass transfer rate; steep gradients near interfaces demand fine meshing in CFD simulations.

Schmidt Number (Sc)

0.1–1000 (gases: 0.2–2; water: ~600; glycerol/water: ~10⁴)

Dimensionless ratio of kinematic viscosity to molecular diffusivity, characterizing relative thickness of hydrodynamic vs. concentration boundary layers.

⚡ Engineering Impact:

Dictates whether mass transfer is controlled by diffusion (Sc ≫ 1) or convection (Sc ≪ 1), guiding correlation selection (e.g., Sherwood number models).

Effective Diffusivity (D_eff)

1e−11 to 1e−8 m²/s (e.g., D_eff ≈ 0.2×D₀ in catalyst pellets with ε=0.4, τ=3)

Diffusion coefficient corrected for tortuosity and porosity in porous media, representing actual transport resistance.

⚡ Engineering Impact:

Critical for sizing fixed-bed reactors and predicting intraparticle limitations — neglect causes >50% error in conversion prediction for exothermic reactions.

📐 Key Formulas

Fick’s First Law

J = -D \frac{\partial c}{\partial x}

Steady-state diffusive molar flux in one dimension

Variables:
Symbol Name Unit Description
J Diffusive molar flux mol/(m²·s) Steady-state molar flux in one dimension
D Diffusion coefficient m²/s Proportionality constant relating flux to concentration gradient
c Concentration mol/m³ Molar concentration of the diffusing species
x Position m Spatial coordinate in the direction of diffusion
Typical Ranges:
O₂ in water at 25°C
1e−9 – 1e−8 mol·m⁻²·s⁻¹
CO₂ in air at 25°C
1e−5 – 1e−4 mol·m⁻²·s⁻¹
⚠️ Valid only when |∂c/∂x| < 0.1·c_avg to avoid nonlinearity effects
🔧 Open Calculator

Fick’s Second Law (1D Cartesian)

\frac{\partial c}{\partial t} = D \frac{\partial^2 c}{\partial x^2}

Transient concentration evolution in a homogeneous medium

Variables:
Symbol Name Unit Description
c concentration mol/m³ or kg/m³ Concentration of the diffusing species
t time s Time variable
D diffusion coefficient m²/s Material property quantifying the rate of diffusion
x position m Spatial coordinate in one dimension
Typical Ranges:
Ethanol diffusion in water (0.1 mm slab)
t_90% ≈ 1–5 s
H₂ in Ni catalyst pellet (1 mm radius)
t_90% ≈ 0.2–2 s
⚠️ Requires D constant; invalid if T or composition varies >5% across domain
🔧 Open Calculator

Effective Diffusivity (porous solid)

D_{\text{eff}} = D \frac{\varepsilon}{\tau}

Corrected diffusion coefficient accounting for porosity (ε) and tortuosity (τ)

Variables:
Symbol Name Unit Description
D_{\text{eff}} Effective Diffusivity m^2/s Corrected diffusion coefficient accounting for porosity and tortuosity
D Diffusion Coefficient m^2/s Intrinsic diffusion coefficient in the fluid phase
\varepsilon Porosity dimensionless Fraction of void volume in the porous solid
\tau Tortuosity dimensionless Measure of the sinuosity of diffusion paths in the porous medium
Typical Ranges:
Activated carbon (ε=0.6, τ=2.5)
D_eff ≈ 0.24·D
Ceramic membrane (ε=0.35, τ=4.0)
D_eff ≈ 0.088·D
⚠️ Only valid when pore diameter > 10× mean free path (continuum regime)
🔧 Open Calculator

🏭 Engineering Example

BASF Ludwigshafen Ammonia Synthesis Loop

N/A — application in gas-phase H₂/N₂ diffusion through Fe₃O₄ catalyst pores
D_eff
2.4e−6 m²/s
Porosity
0.42
D_H2_in_N2
7.2e−5 m²/s at 450°C, 15 MPa
Tortuosity
3.1
Thiele_modulus
12.7
Pellet_diameter
3.2 mm

🏗️ Applications

  • Design of absorption columns for CO₂ capture
  • Modeling drug release kinetics from polymeric implants
  • Predicting oxygen depletion in anaerobic bioreactors
  • Optimizing dopant diffusion in silicon wafers

🔧 Try It: Interactive Calculator

🔧
Fluid Flow & Transport Phenomena Calculator
Multi-tab calculator for Fluid Flow & Transport Ph...
Open →

📋 Real Project Case

Hydrocarbon Separation in Offshore Gas Processing Skid

Integrated gas processing module for North Sea platform

Challenge: Insufficient liquid carryover removal causing downstream compressor fouling
Vertical Separator Skid LayoutInletVaneSeparatorGas OutQ_g = 12,500 Sm³/hLiq Outv_t = 0.18 m/sCarryover160 mm120 mmHydrocarbon Separation Skid
Read full case study →
Slurry Transport Optimization in Copper Mine Tailings Pipeline Heat Integration in Ethylene Oxide Absorption Column CFD-Guided Mixer Redesign in Pharmaceutical Bioreactor

🎨 Technical Diagrams

High concentrationLow concentrationDiffusion direction →
c(x,0) = c₀c(0,t) = c_sc(x,t) solution curve

📚 References

[1]
[2]
[3]