Ergun Equation: A Thorough British Guide to the Classic Pressure‑Drop Correlation

PortalAdmin Misc 6. May 2025 | 0

Introduction to the Ergun Equation

The Ergun equation stands as a cornerstone of chemical engineering, tying together fluid properties, particle geometry and bed structure to predict how difficult it is for a fluid to pass through a packed bed. Often introduced as the Ergun equation, this compact correlation blends a linear (viscous) term with a quadratic (inertial) term to describe the pressure drop per unit length in packed beds of spheres or similar shapes. In Australia, the United States, and across Europe, engineers use this relation to design catalytic reactors, filtration columns, and packed-bed absorbers. In the context of literature and lecture notes, you may encounter the Ergun relation, the Ergun correlation, or the Ergun equation; all refer to the same foundational model, albeit expressed with varying emphasis.

Put simply, the Ergun equation provides a practical bridge between fundamental fluid mechanics and the real world of packed beds. It recognises that at low velocities the flow is dominated by viscous forces (like climbing through a narrow corridor), while at higher velocities inertial effects take over (akin to pushing through a crowd). The elegance of the Ergun equation lies in its compact form: a sum of two terms, each scaling differently with velocity and bed properties, that together capture a broad range of operating conditions.

Historical context and foundational concepts

The origin of the Ergun equation traces back to a 1952 publication by Mehmet Ergun and their co‑authors, who sought a unified description for pressure drop in packed beds comprising spherical or near‑spherical particles. Before this, separate Darcy‑type relationships described laminar flow in porous media, while other models handled high‑velocity regimes with varying degrees of confidence. The achievement of Ergun was to propose a single, practical formula that could be applied across Reynolds numbers from a few tenths to well above one, incorporating porosity, particle size and density alongside viscosity and flow rate. Since then, the Ergun equation has become ubiquitous in process design courses, supplier handbooks and practical design codes in chemical engineering.

In modern practice, it is common to present the formula in the form most often used for gas-phase flows, but the underlying principles are equally valid for liquids, provided the properties are adjusted appropriately. The classic Ergun equation is robust yet simple enough to be implemented in spreadsheet tools, process simulators, and quick‑check calculations during the early design stages of a project.

The mathematical formulation: the Ergun equation explained

The Ergun equation links the pressure drop ΔP across a packed bed to the superficial velocity of the gas or liquid, the bed porosity, particle size and the properties of the fluid. The conventional expression is written as:

ΔP / L = 150 μ (1 − ε)^2 / (ε^3 D_p^2) · v + 1.75 ρ (1 − ε) / (ε^3 D_p) · v^2

where the variables denote:

ΔP: pressure drop across the bed (Pa)
L: bed length (m)
μ: dynamic viscosity of the fluid (Pa·s)
ε: void fraction or bed porosity (dimensionless)
ρ: density of the fluid (kg/m^3)
D_p: characteristic particle diameter (m)
v: superficial velocity, defined as volumetric flow rate divided by the cross‑sectional area of the bed (m/s)

Two distinct terms make up the right‑hand side: the first term is linear in velocity and arises from viscous shear within the pore channels, while the second term is quadratic in velocity, reflecting inertial losses associated with flow separation, eddies and jetting around particles. The balance between these terms shifts as flow rate increases, which is why the Ergun equation remains reliable across a broad Reynolds number spectrum.

For completeness, some texts present the Ergun equation with the two components grouped as:

ΔP / L = f_L · v + f_Q · v^2

where f_L and f_Q are coefficients dependent on the bed porosity ε, particle diameter D_p and fluid properties μ and ρ. This alternative framing emphasises the additive nature of the viscous and inertial penalties to pressure drop.

How to interpret the Ergun equation: physical meaning of the terms

The viscous term: the linear component

The viscous term, 150 μ (1 − ε)^2 / (ε^3 D_p^2) · v, dominates at low superficial velocities. Here, μ represents the resistance imposed by the fluid’s internal friction as it negotiates the tiny passages through the packed bed. Smaller particles (reducing D_p) and lower porosities (smaller ε) intensify viscous losses, because the flow path becomes more tortuous and constrained. In practice, if you design a bed with fine particles, you should expect higher pressure drop for a given superficial velocity due to this linear contribution.

The inertial term: the quadratic component

The quadratic term, 1.75 ρ (1 − ε) / (ε^3 D_p) · v^2, captures the energy dissipation caused by inertial effects as the fluid accelerates and decelerates around particles. At higher flow rates, this term grows in importance. Larger particles (larger D_p) tend to lower the inertial losses because the flow lines can pass more smoothly around each particle, reducing turbulence and jetting. Likewise, higher porosities decrease the density of obstruction, moderating the inertial penalties.

Key parameters and how to obtain them

To implement the Ergun equation accurately, you need reliable values for several parameters:

Porosity, ε: the fraction of bed volume that is void. This is not simply the bulk porosity of the packing; it reflects how tightly packed the bed is in its actual operating state. Typical values for random close packings range around 0.35–0.40, but catalytic pellets or structured packings may vary widely.
Particle diameter, D_p: a characteristic length that represents the size of the packing elements. For non-uniform particles, an equivalent diameter is used to approximate the effective flow path length.
Viscosity, μ: a property of the fluid, often measured at operating temperature. In gas mixtures, a viscosity correction with temperature is essential.
Density, ρ: the fluid’s mass per unit volume. Gas densities can vary significantly with pressure and temperature, whereas liquids are more stable but still require temperature corrections.
Superficial velocity, v: calculated from the volumetric flow rate and the cross‑sectional area of the bed. In packed beds, the superficial velocity is typically used rather than the actual seepage velocity, because the Ergun equation is framed for the bulk flow that would exist if the bed were fully open.

Accurate parameter estimation is essential for meaningful predictions. In practice, ε is often measured using packed beds after packing and expansion, or estimated from packing method and particle shape using empirical correlations. For irregular or non‑spherical particles, effective diameters and porosity can diverge from those of ideal spheres, and practitioners adjust D_p or employ correction factors to reflect the true flow resistance.

Practical uses: when and where the Ergun equation shines

The Ergun equation is especially valuable in the following areas of chemical engineering practice:

Design of packed-bed reactors and catalytic beds where gas or liquid must diffuse through a bed of solid catalysts. The equation provides a first‑order estimate of pressure drop to inform pump sizing and energy requirements.
Filtration and liquid‑solid separation processes, where a packed bed of media controls flow resistance and process efficiency. Here, the viscous and inertial components help predict how filter media will perform under different operational loads.
Gas distribution in absorber columns and contactors, where evaluating pressure drop ensures uniform distribution and avoids channeling or bypassing within the bed.
Scale‑up studies, where data from a labbed bed are translated to pilot or industrial scales. The Ergun equation, with bed property adjustments, supports the translation from bench to plant.

Applying the Ergun equation in practice: a step‑by‑step guide

To use the Ergun equation effectively, follow these practical steps:

Characterise the packing: determine ε and D_p. If the bed uses non‑uniform pellets or structured packings, obtain the equivalent diameter and porosity from supplier data or experimental measurements.
Measure operating fluid properties: obtain μ and ρ for your fluid at the process temperature. For gas mixtures, a mixture viscosity and a density derived from composition and temperature may be needed.
Determine the bed length L and superficial velocity v: calculate v = Q / A, where Q is the volumetric flow rate and A is the cross‑sectional area of the bed. Ensure Q reflects the desired operating conditions.
Compute the two terms: substitute ε, D_p, μ, ρ, and v into the viscous and inertial components of the Ergun equation. Add them to obtain ΔP / L, the pressure drop per unit length.
Verify units and consistency: Pa/m is a common unit for pressure drop per length in the metric system. If other units are used, convert carefully to obtain the correct ΔP/L value.
Analyse sensitivity: explore how changes in ε, D_p, or v influence ΔP/L. This helps in design optimization and in assessing how close you are to limits for pumping power and energy efficiency.

Modified formulations and extensions of the Ergun equation

In practice, engineers sometimes adapt the Ergun equation to fit specific systems or regimes. Some common modifications include:

Gas–liquid systems: When a two‑phase flow occurs, the Ergun equation can be applied to each phase with phase‑specific viscosities and densities, or with an effective mixture property approach. You may also encounter the Ergun equation used in conjunction with holdup correlations to account for phase distribution within the bed.
Structured packings and non‑spherical particles: For uniform, non‑spherical particles or structured packing elements, coefficients may be adjusted to reflect different flow resistance patterns. Empirical fits or literature data are often used to refine the constants 150 and 1.75 for a given packing.
Low Reynolds number regimes: In very fine packings, Reynolds numbers may be small enough that the viscous term dominates. The equation still holds, but validation against experimental data remains important to ensure accuracy.
High Reynolds number regimes: At very high superficial velocities, inertial effects become more pronounced, and practice may include correction factors or a transition to turbulence models in more advanced simulations. The Ergun equation remains a starting point, with possible refinements for exact predictions.

Limitations and caveats of the Ergun equation

While the Ergun equation is widely used and highly valuable, it is not without limitations. Understanding these helps engineers avoid over‑reliance on a single equation at the design stage:

Assumes steady, incompressible flow: In highly compressible gas flows or rapidly changing conditions, more sophisticated models may be required to capture density changes and transient effects.
Homogeneous bed assumption: The equation presumes a uniform packing. Beds with large gradients in particle size, porosity or wetting can yield deviations from predicted pressure drops.
Porosity and particle size measurement challenges: Real beds often exhibit local variations in ε or D_p. Averaged values may smooth out critical local changes that influence pressure drop.
Temperature effects: Varying temperatures influence viscosity and density, changing the balance of terms. Consistent temperature control and property correlations are essential for accuracy.
Scale effects: Extrapolating from lab or pilot‑scale data to industrial scale should be done with care, as packing behaviour, channeling and bed uniformity can change with size.

Comparisons with related models: where the Ergun equation sits in the toolbox

Several families of correlations are used to describe flow through porous media. The Ergun equation sits alongside the Darcy and Kozeny–Couette relations, but offers a broader Reynolds number applicability due to its inertial term. Here are key comparisons to keep in mind:

Darcy’s law: Darcy’s law is a linear relation between pressure gradient and superficial velocity, valid primarily for laminar flow in highly porous media. The Ergun equation extends beyond this laminar regime by incorporating an inertial term that captures non‑linear losses at higher flow rates.
Kozeny–Cuppe: Kozeny–Cuppe provides a relation that links permeability to porosity and specific surface area, often used for packed beds with homogeneous particles. The Ergun equation complements this by directly relating pressure drop to flow, allowing practical design calculations in conjunction with Kozeny‑Cuppe derived parameters.
Ergun vs other modern correlations: For structured packings or specialized catalysts, designers might adopt alternative correlations or computational tools that account for complex geometries. However, the Ergun equation remains a robust first‑order predictor in many systems.

Illustrative example: applying the Ergun equation to a practical scenario

Imagine a gas‑solid catalytic bed used for a hydrocarbon oxidation process. The bed consists of spherical catalyst pellets with an average diameter D_p of 4 mm and a porosity ε of 0.38. The gas flowing through the bed has a dynamic viscosity μ of 2.0 × 10−5 Pa·s and a density ρ of 0.8 kg/m^3. The bed length is L = 2.5 m and the cross‑sectional area is 0.5 m^2. The desired superficial velocity is v = 0.25 m/s.

First, compute the viscous term:

Coefficient for viscous term = 150 × μ × (1 − ε)^2 / (ε^3 × D_p^2)

Plug in numbers: ε = 0.38, D_p = 0.004 m, μ = 2.0 × 10−5 Pa·s

(1 − ε) = 0.62; (1 − ε)^2 = 0.3844; ε^3 ≈ 0.054872; D_p^2 = 1.6 × 10−5 m^2

Coefficient ≈ 150 × 2.0×10−5 × 0.3844 / (0.054872 × 1.6×10−5) ≈ 0.00001? (illustrative) In practice, you would compute precisely; this coefficient multiplies v = 0.25 m/s to obtain the viscous contribution to ΔP/L.

Next, the inertial term:

Coefficient for inertial term = 1.75 × ρ × (1 − ε) / (ε^3 × D_p)

Plug in numbers: ρ = 0.8 kg/m^3; (1 − ε) = 0.62; ε^3 ≈ 0.054872; D_p = 0.004 m

Coefficient ≈ 1.75 × 0.8 × 0.62 / (0.054872 × 0.004) ≈ 0.868 / 0.0002195 ≈ 3953 (units: Pa·s^2/m^2 per m/s^2?)

Now, insert v = 0.25 m/s and then ΔP/L = viscous term × v + inertial term × v^2. In practice, you would perform the exact arithmetic with a calculator or a software tool to obtain ΔP/L. Finally, multiply by bed length L = 2.5 m to obtain ΔP. This simple calculation yields a preliminary estimate of pressure drop across the bed, which informs pump sizing and energy requirements for the process. This example illustrates how the two terms interact: at this operating point, both viscous and inertial effects contribute meaningfully to the total pressure drop.

Choosing the right approach for your project

When selecting the Ergun equation for a project, consider the following guidelines:

If your bed uses uniform, spherical particles with a well‑defined diameter and porosity, the standard Ergun equation is typically adequate for initial design work.
For beds with a mix of particle sizes or non‑spherical shapes, use an effective D_p and ε based on experimental data or industry references, and validate with pilot data where possible.
For gas‑liquid systems, apply the equation separately to each phase or adopt a refined approach that accounts for holdup and phase distribution within the bed.
Always perform sensitivity analyses: small changes in porosity or particle size can lead to large changes in ΔP/L, especially in the viscous term. This helps identify robust operating windows and safe design margins.

Common pitfalls and how to avoid them

To ensure robust results, be mindful of these typical pitfalls:

Misidentifying porosity: avoid using the bulk porosity of the packing material without considering actual bed packing. Measure or estimate ε for the operational bed to reduce error.
Using inconsistent units: ensure consistency across all terms. Converting μ, ρ, D_p and v to compatible units prevents dimensional analysis mistakes that can produce huge errors in ΔP/L.
Ignoring temperature effects: viscosity is highly temperature‑dependent. In processes with significant temperature gradients, use temperature‑corrected μ values for each segment of the bed.
Neglecting transition regimes: in some flows, the Reynolds number spans transitional regimes. Where possible, compare Ergun predictions with experimental data or supplementary correlations to confirm accuracy.

Why engineers still rely on the Ergun equation today

Despite advances in computational fluid dynamics and numerical modelling, the Ergun equation remains a practical favourite for several reasons. It is:

Simple and implementable: a compact formula with a small set of clearly defined parameters.]
Versatile: applicable to gas and liquid flows, for a range of bed types from catalytic reactors to filtration media.
Empirically grounded: backed by decades of experimental data, making it a reliable starting point for design and optimisation.
Easy to calibrate: coefficients can be adjusted based on target bed configurations and operating conditions, allowing alignment with specific plant data.

Subsections and further considerations: exploring deeper nuances

Porosity versus void fraction: clarifying the terms

In the Ergun equation context, ε denotes porosity or the fraction of bed volume that is void. This can be confused with the void fraction in other formulations, but in practice, ε is the parameter that controls how congested the bed channels are. The smaller the porosity, the more constricted the flow paths, and the larger the pressure drop for a given superficial velocity. In some cases, a distinction is made between bed porosity and emulsion fraction in multiphase systems; ensure you are using the correct one for your problem setup.

Particle size distribution and shape effects

Real beds seldom contain perfectly uniform spheres. If your packing features a range of particle sizes or non‑spherical shapes, you may need to adapt the effective diameter or employ empirical correction factors. For highly irregular media, manufacturers’ data or literature reviews can guide you in selecting an appropriate effective D_p that reflects the mean flow resistance.

Using the Ergun equation with liquids

When liquids are involved, similarly to gases, the Ergun equation is applicable but with liquid properties. Because liquids often exhibit higher densities and viscosities that differ from gases, the balance between viscous and inertial losses shifts. In some liquids, particularly Newtonian fluids with constant properties, the same form holds. For non‑Newtonian liquids, more complex models may be needed to capture shear‑thinning effects or viscoelastic behaviours.

Wrapping up: a practical roadmap for engineers and students

The Ergun equation remains a foundational tool in the chemical engineer’s toolkit. It provides a clear, testable, and implementable approach to predicting pressure drop through packed beds, bridging the gap between simple Darcy flow and the complexities of high‑velocity regimes. By understanding the two terms—the viscous and inertial components—you gain insight into how bed porosity, particle size and flow rate shape the pressure landscape inside a packed bed. While modern design often leans on computational simulations for detailed insights, the Ergun equation continues to offer a rapid, transparent framework for preliminary design, feasibility studies and decision‑making in process development.

Key takeaways: the Ergun equation in one page

– The Ergun equation combines a viscous term and an inertial term to predict pressure drop in packed beds.

– Porosity, particle size, fluid viscosity and density all influence ΔP/L through the two terms.

– It applies to both gas and liquid flows, with appropriate property inputs.

– For complex beds, use effective parameters and validate with experimental data.

– It is a timeless tool: simple to apply, widely validated, and a reliable starting point for design decisions.

Final reflections: the enduring value of the Ergun equation

In a field that continually introduces new models and numerical methods, the Ergun equation endures because of its blend of physical intuition, practicality, and empirical grounding. It is the hinge point between theory and practice, enabling engineers to estimate energy requirements, set performance targets, and explore design spaces with confidence. For students and professionals alike, mastering the Ergun equation — and knowing when to apply its assumptions and limitations — is a pivotal step in developing robust, efficient, and safe process systems.

Glossary: quick definitions for the Ergun equation’s terms

Ergun equation: the pressure drop correlation for packed beds, combining viscous and inertial effects.

ε (epsilon): bed porosity or void fraction.

D_p: particle diameter, an effective diameter for flow through the bed.

μ: dynamic viscosity of the fluid.

ρ: fluid density.

v: superficial velocity through the bed.