Ergun equation - Packed Bed Head Loss

Select calculation method

Input data

Pressure drop

Cross-sectional area

Bed length

Particle diameter

Sphericity

Bed porosity

Fluid density

Fluid viscosity

Output data

Volumetric flow

Superficial velocity

Modified particle Reynolds number

Description

Pressure drop of a single-phase fluid flowing through a fixed bed of particles, by the Ergun equation. The equation adds a viscous term, linear in velocity, to an inertial term, quadratic in velocity, and so covers the laminar, transitional and turbulent regimes with one expression. It applies to packed columns, catalytic reactors, adsorbers, granular filters and beds at the onset of fluidization.

Ergun equationErgun equation

Calculation methods

  • Pressure drop: pressure drop and pressure gradient of a bed for a given flow, with the superficial velocity, the modified Reynolds number and the share of the viscous term.
  • Flow for a given pressure drop: flow that a bed passes under a given pressure drop.
  • Minimum fluidization: superficial velocity and flow at which the bed's pressure drop equals its buoyant weight, with the bed pressure drop at that point.

Model

ΔPL=150μ(1−ε)2ϕ2dp2ε3u+1.75ρ(1−ε)ϕdpε3u2u=QA\frac{\Delta P}{L} = \frac{150 \mu (1 - \varepsilon)^2}{\phi^2 d_p^2 \varepsilon^3} u + \frac{1.75 \rho (1 - \varepsilon)}{\phi d_p \varepsilon^3} u^2 \qquad u = \frac{Q}{A}

Re′=ρuϕdpμ(1−ε)Re' = \frac{\rho u \phi d_p}{\mu (1 - \varepsilon)}

where ΔP\Delta P is the pressure drop across a bed of length LL and cross-sectional area AA, QQ the volumetric flow, uu the superficial velocity, ε\varepsilon the bed porosity, dpd_p the volume-equivalent particle diameter, ϕ\phi the particle sphericity (ϕ=1\phi = 1 for spheres), and ρ\rho and μ\mu the fluid density and viscosity. The viscous term governs for Re′<10Re' < 10 and the inertial term for Re′>1000Re' > 1000.

At minimum fluidization the pressure drop equals the buoyant weight of the bed per unit area:

ΔPmf=(1−εmf)(ρp−ρ)gLmf\Delta P_{mf} = (1 - \varepsilon_{mf}) (\rho_p - \rho) g L_{mf}

1.75εmf3ϕRemf2+150(1−εmf)εmf3ϕ2Remf=ArAr=dp3ρ(ρp−ρ)gμ2Remf=ρumfdpμ\frac{1.75}{\varepsilon_{mf}^3 \phi} Re_{mf}^2 + \frac{150 (1 - \varepsilon_{mf})}{\varepsilon_{mf}^3 \phi^2} Re_{mf} = Ar \qquad Ar = \frac{d_p^3 \rho (\rho_p - \rho) g}{\mu^2} \qquad Re_{mf} = \frac{\rho u_{mf} d_p}{\mu}

The correlation of Wen and Yu, Remf=33.72+0.0408Ar−33.7Re_{mf} = \sqrt{33.7^2 + 0.0408 Ar} - 33.7, which requires neither εmf\varepsilon_{mf} nor ϕ\phi, is shown in the calculation steps for comparison.

Range of application

  • Incompressible flow. For a gas, the result holds while the pressure drop is small relative to the absolute pressure (roughly below 10 %); otherwise the bed should be divided into segments evaluated at their mean density.
  • The porosity has the strongest influence on the result: ΔP\Delta P varies approximately with (1−ε)2/ε3(1 - \varepsilon)^2 / \varepsilon^3, so an error of 0.02 in ε\varepsilon changes the pressure drop by about 20 %. Randomly packed spheres have ε≈0.36\varepsilon \approx 0.36 to 0.400.40.
  • Wall effects are neglected; they become significant when the bed diameter is less than about ten particle diameters.
  • Default values: in the first two methods, a bed of 3 mm spheres in air at 20 °C and 1 atm; in the minimum fluidization method, round sand (ϕ=0.86\phi = 0.86, 2650 kg/m³) in air.

References

  • ERGUN, S. Fluid flow through packed columns. Chemical Engineering Progress, v. 48, n. 2, p. 89-94, 1952.
  • KUNII, D.; LEVENSPIEL, O. Fluidization Engineering. 2. ed. Boston: Butterworth-Heinemann, 1991.
  • WEN, C. Y.; YU, Y. H. A generalized method for predicting the minimum fluidization velocity. AIChE Journal, v. 12, n. 3, p. 610-612, 1966.
  • MCCABE, W. L.; SMITH, J. C.; HARRIOTT, P. Unit Operations of Chemical Engineering. 7. ed. New York: McGraw-Hill, 2005.
  • BIRD, R. B.; STEWART, W. E.; LIGHTFOOT, E. N. Transport Phenomena. 2. ed. New York: Wiley, 2002.