Understanding how particles settle in a fluid is fundamental to fields ranging from geology and environmental science to chemical engineering and pharmaceutics. The sedimentation velocity—the terminal speed at which a particle descends under gravity through a liquid or gas—determines the efficiency of clarifiers, the transport of sediment in rivers, the formation of sedimentary rocks, and the separation of cells in bioprocessing. An accurate prediction of this velocity requires accounting for particle size, density, shape, fluid viscosity, and the flow regime. This article provides a comprehensive guide to calculating sedimentation velocity for different particle types, from ideal spheres to irregular aggregates, and explores the governing equations, corrections, and practical measurement techniques.

The Physics of Sedimentation: Key Parameters

Sedimentation occurs when the gravitational force acting on a particle exceeds the buoyant force and the viscous drag resisting its motion. At terminal velocity, the net gravitational force equals the drag force. The primary parameters that influence this equilibrium are:

  • Particle density (ρp) and fluid density (ρf) — the density difference Δρ = ρp − ρf provides the driving force.
  • Particle size — usually expressed as an equivalent diameter. For natural sediments, the sieve diameter or Stokes diameter is commonly used.
  • Fluid viscosity (η) — the resistance to flow, which depends on temperature and fluid composition.
  • Particle shape — non-spherical particles experience different drag and often settle slower than spheres of the same mass.
  • Flow regime — characterized by the particle Reynolds number (Re = ρf v d / η). Low Re corresponds to laminar flow (Stokes regime), while high Re indicates turbulent flow.

These parameters are combined in various empirical and theoretical formulas to estimate the terminal settling velocity. The simplest and most widely used starting point is Stokes' Law, valid for small, spherical particles settling in a quiescent fluid under laminar conditions.

Stokes' Law and Its Assumptions

Stokes' Law describes the drag force on a sphere moving at low Reynolds number (Re < 0.1). The terminal settling velocity v for a spherical particle is derived by equating the net gravitational force (weight minus buoyancy) to the Stokes drag force:

v = (2/9) · (r²) · (Δρ) · g / η

where r is the particle radius, g is the gravitational acceleration, and the other symbols are as defined above. This equation applies only when the following assumptions are met:

  • The particle is a rigid sphere.
  • The fluid is incompressible and Newtonian.
  • The flow is laminar (Re < 0.1), meaning viscous forces dominate inertial forces.
  • The particles are far from walls and from each other (i.e., dilute suspension).
  • There is no particle-particle interaction, flocculation, or hindered settling.

In practice, many natural and industrial particles violate one or more of these assumptions. For example, sand grains are irregular, microbial cells may deform, and flocs have fractal structures. Nevertheless, Stokes' Law remains a valuable baseline from which corrections can be applied. A detailed derivation of Stokes' Law can be found in standard fluid mechanics textbooks or online resources such as Wikipedia’s article on Stokes’ law.

Calculating Velocity for Spherical Particles

When dealing with approximately spherical particles (e.g., glass beads, uniform latex spheres, or oil droplets), Stokes' Law gives a direct estimate. For example, consider a silica sand particle with radius 50 μm (5 × 10−5 m), density 2650 kg/m³, settling in water at 20°C (ρf = 998 kg/m³, η = 1.002 × 10−3 Pa·s). The density difference Δρ = 1652 kg/m³. Using g = 9.81 m/s²:

v = (2/9) · (5 × 10−5)² · 1652 · 9.81 / (1.002 × 10−3)

This yields v ≈ 0.009 m/s (9 mm/s). The Reynolds number Re = (998 · 0.009 · 1 × 10−4) / 0.001 ≈ 0.9, which is slightly above the Stokes regime limit (0.1), so a correction for moderate Reynolds numbers may be warranted for greater accuracy. For particles smaller than about 30 μm, Re is typically well below 0.1, and Stokes' Law applies directly.

Beyond Stokes: Transition and Turbulent Regimes

As particle size increases, inertial effects become significant, and the drag coefficient CD deviates from the Stokes value (CD = 24/Re). For intermediate Reynolds numbers (0.1 < Re < 1000), the relationship is described by the standard drag curve. In this regime, the terminal velocity can be found by solving the force balance iteratively or using empirical correlations such as the Schiller-Naumann correlation: CD = (24/Re)(1 + 0.15 Re0.687) for Re < 800.

For large, fast-settling particles (Re > 1000), the drag coefficient becomes approximately constant (CD ≈ 0.44 for smooth spheres). The settling velocity in the Newton’s law regime is independent of viscosity and given by:

v = sqrt( (4 g d (ρp − ρf)) / (3 ρf CD) )

where d is the particle diameter. For a gravel particle of diameter 2 mm sinking in water, Re will be on the order of hundreds, and the turbulent regime equation should be used. Many engineering design guides, such as those from the Engineering Toolbox, provide calculators and tables for these calculations.

Worked Example: Transition Regime

Take a quartz sphere of diameter 0.5 mm (ρp = 2650 kg/m³) in water. Start by guessing a velocity, compute Re, find CD, then solve for v. An initial guess from Stokes (ignoring Re correction) gives v ≈ 0.09 m/s, Re ≈ 45. Using Schiller-Naumann, CD ≈ 1.1. Plugging into the force balance yields v ≈ 0.074 m/s. Repeating the iteration gives convergence around v = 0.072 m/s, Re = 36, CD = 1.25. The iterative correction is essential when Stokes' Law overpredicts velocity by more than a few percent.

Adjustments for Non-Spherical Particles

Most natural and industrial particles are not perfect spheres. Their irregular shape increases the drag force, reduces the settling velocity, and makes characterization more complex. The standard approach is to define a shape factor (S or ψ) that relates the drag of the particle to that of a sphere of the same volume. Common shape factors include:

  • Sphericity (ψ): Ratio of the surface area of a sphere of the same volume as the particle to the actual surface area of the particle. Sphericity ranges from 0 (extremely elongated or flaky) to 1 (perfect sphere).
  • Corey shape factor: Used for sediment particles; defined as c / sqrt(ab) where a, b, c are the longest, intermediate, and shortest dimensions.
  • Dynamic shape factor: The ratio of the drag force on the actual particle to that on a sphere of the same volume and density at the same settling velocity. This factor is often determined empirically.

For non-spherical particles, an equivalent diameter (such as the volume-equivalent diameter dv or the Stokes diameter dst) is used in the Stokes equation, and then the calculated velocity is multiplied by a correction factor typically less than 1. For example, sand grains with a sphericity of 0.7 may settle at 60–80% of the velocity predicted for a sphere of the same mass and density. Detailed tables for shape correction factors are available in sedimentology references, such as studies on the settling of natural sediments.

Empirical Correlations for Natural Sediments

Many researchers have proposed equations for the settling velocity of natural sand and silt particles that incorporate shape. One of the most widely used is the formula by Soulsby (1997) for fine sand to coarse silt, which uses a drag-law approach but replaces the diameter with a nominal diameter and applies an empirical coefficient. Another often-cited correlation is the Ferguson-Church equation:

v = (g Δρ d²) / (C1 η + sqrt( C2 g Δρ d³ ) )

where C1 and C2 are empirical constants that depend on particle shape (e.g., C1 ≈ 18 for natural sand, C2 ≈ 0.4 to 1.0 for spherical to angular particles). This equation smoothly transitions from the Stokes regime (small d) to the turbulent regime (large d) and inherently accounts for shape via the constants.

Hindered Settling and Concentrated Slurries

In many industrial processes (e.g., sedimentation in thickeners, wastewater clarifiers), the particle concentration is too high to ignore interactions. Under hindered settling, the downward motion of each particle is impeded by the upward displacement of fluid and collisions with neighbors. The effective settling velocity vh is less than the terminal velocity of an isolated particle. The most popular model is the Richardson-Zaki equation:

vh = vt · (1 − φ)n

where φ is the volume fraction of solids, and n is an exponent that depends on the particle Reynolds number at infinite dilution (Ret). For Ret < 0.2 (laminar), n ≈ 4.65; for intermediate Re, n decreases, and for fully turbulent (Ret > 500), n ≈ 2.4. The Richardson-Zaki correlation is extensively used in design of sedimentation tanks and is discussed in detail in hindered settling theory.

Effect of Flocculation

In water treatment and many natural environments, particles are not individual but form aggregates (flocs) that have a fractal structure. Floc density decreases with increasing floc size because of internal porosity, and the floc settling velocity often follows a power-law relationship with floc diameter. For example, clay flocs may have settling velocities described by v = k dm where m is around 1–2 (much lower than the power 2 expected from Stokes' Law). The effective density of a floc ρfloc can be approximated as ρf + (ρp − ρf) (dp/dfloc)Df − 3, where Df is the fractal dimension (typically 1.7–2.3 for aggregates). This fractal nature complicates exact prediction, so experimental measurement of floc settling is often required.

Practical Methods to Measure Sedimentation Velocity

While theoretical and empirical formulas are useful, direct measurement is often necessary for heterogeneous or irregular particles. Common techniques include:

  • Andreasen pipette method: A classical technique for sediments; a sample is allowed to settle in a calibrated cylinder, and aliquots are withdrawn at different times to determine particle size distribution via cumulative mass.
  • Sedimentation balance: The mass of settled material on a submerged pan is recorded over time, and the velocity is derived from the accumulation curve.
  • Laser diffraction systems: Instruments such as the Malvern Mastersizer use the principle of laser diffraction to infer particle size distribution from which settling velocities can be calculated using the appropriate drag model.
  • Video microscopy: Direct observation of particle settling in a cell enables tracking of individual particles, particularly useful for studying irregular or flocculent particles.

When measuring, it is critical to control temperature (viscosity changes ≈2% per degree Celsius for water) and to use dilute enough suspensions to avoid hindered settling.

Applications Across Industries

Accurate sedimentation velocity calculations are vital in many fields:

  • Water and wastewater treatment: Design of primary clarifiers, secondary settling tanks, and sludge thickeners relies on settling velocity distributions to achieve required effluent quality.
  • Mining and mineral processing: Classification of ores by size in hydrocyclones and thickeners uses hindered settling models; tailings pond design requires knowledge of how fine particles consolidate.
  • Oceanography and sediment transport: Predicting the fate of sediment plumes from rivers, dredging operations, or deep-sea mining depends on settling velocities of sand, silt, and clay, including flocculation in saline water.
  • Pharmaceuticals and biotechnology: Cell harvesting in centrifugation and sedimentation; also for characterizing the size of protein aggregates or emulsions.
  • Environmental engineering: Modeling the removal of pollutants in constructed wetlands, stormwater detention basins, and oil-water separators.

In each case, the particle type dictates which model is appropriate. For example, in wastewater clarifiers where biological flocs are present, the fractal floc model combined with hindered settling is used, whereas for grit removal in the same plant, Stokes' Law with shape correction is applied.

Conclusion

Calculating sedimentation velocity begins with Stokes' Law for simple spherical particles under laminar conditions. For larger or faster-moving particles, the drag coefficient must be adjusted for transition or turbulent flow. Irregular shapes are handled by shape factors and empirical correlations such as the Ferguson-Church equation. In concentrated suspensions, hindered settling models (Richardson-Zaki) are necessary, and for flocculating systems, the fractal nature of aggregates must be considered. By systematically applying the appropriate model for the particle type and flow regime, engineers and scientists can accurately predict sedimentation behavior, optimize separation processes, and better understand natural sediment dynamics.