Settling Velocity Stokes Interactive Calculator

The Settling Velocity Stokes Calculator determines the terminal velocity at which spherical particles settle through a fluid under laminar flow conditions. This fundamental calculation is essential for designing sedimentation tanks, water treatment systems, air pollution control equipment, and mineral processing facilities. Environmental engineers, chemical process designers, and hydrologists use Stokes' law to predict particle behavior in both natural and industrial settings.

📐 Browse all free engineering calculators

Diagram

Settling Velocity Stokes Calculator

Equations

Theory & Engineering Applications

Fundamental Principles of Stokes Settling

Stokes' law, derived by George Gabriel Stokes in 1851, describes the force balance on a spherical particle settling through a viscous fluid under laminar flow conditions. When a particle is released in a quiescent fluid, three forces act upon it: gravitational force pulling it downward, buoyant force pushing it upward, and drag force resisting its motion. At terminal (settling) velocity, these forces reach equilibrium, and the particle descends at constant speed.

The gravitational force equals the particle weight minus the buoyant force, giving a net downward force of (π/6)d³g(ρ_p - ρ_f). The viscous drag force, according to Stokes, equals 3πμdv_s for creeping flow around a sphere. Setting these equal and solving for velocity yields the classical Stokes equation. This derivation assumes: spherical particles, infinite fluid extent, no-slip boundary condition at the particle surface, negligible inertial effects (Re ≪ 1), no particle-particle interaction, and rigid particles that don't deform.

Critical Non-Obvious Limitation: Wall Effects

A frequently overlooked constraint in applying Stokes' law is the wall effect in confined geometries. The standard equation assumes an infinite fluid medium, but in laboratory settling columns or industrial clarifiers with finite dimensions, the presence of walls significantly reduces settling velocity. For a sphere settling along the centerline of a cylindrical tube, the velocity reduction follows Ladenburg's correction: v_actual = v_Stokes[1 - 2.104(d/D)], where D is the tube diameter. When the particle diameter exceeds 1/10 of the container diameter, the error becomes substantial. In a tube where d/D = 0.1, the actual settling velocity is approximately 79% of the predicted Stokes velocity—a 21% reduction that can invalidate settling basin designs or particle size analyses if not accounted for.

Temperature Dependence and Viscosity Variations

Fluid viscosity exhibits strong temperature sensitivity, particularly for liquids. Water viscosity decreases from 1.307 mPa·s at 10°C to 0.653 mPa·s at 40°C—a reduction of 50%. Since settling velocity is inversely proportional to viscosity, a particle settling in 40°C water moves approximately twice as fast as in 10°C water, assuming all other parameters remain constant. This temperature sensitivity has critical implications for wastewater treatment plants experiencing seasonal temperature fluctuations, mineral processing circuits in different climates, and any process where temperature control affects separation efficiency.

For precise work, engineers must use temperature-corrected viscosity values. The Andrade equation provides reasonable approximation for water: μ = A·exp(B/T), where T is absolute temperature and A, B are empirical constants. Many computational tools include built-in viscosity databases, but field engineers should maintain awareness that a 10°C temperature change can alter settling rates by 25-35% for water-based systems.

Particle Shape Factors and Non-Spherical Corrections

Real particles rarely exhibit perfect sphericity. Clays, mineral fragments, microplastics, and biological materials possess irregular geometries that increase drag beyond the Stokes prediction. Engineers employ shape factors to account for this deviation. The sphericity (ψ) is defined as the ratio of surface area of a sphere with equivalent volume to the actual particle surface area. For irregular particles, ψ ranges from 0.6 to 0.9, while perfect spheres have ψ = 1.0.

Wadell proposed using an equivalent diameter based on volume: d_eq = (6V/π)^1/3. The settling velocity of non-spherical particles can be approximated as v_actual = v_Stokes·ψⁿ, where the exponent n typically ranges from 1.5 to 2.0 depending on particle orientation and Reynolds number. Elongated particles like fibers settle more slowly than equivalent spheres due to increased drag, while disk-shaped particles exhibit orientation-dependent settling, complicating predictions further.

Transition to Higher Reynolds Numbers

As particles grow larger, velocities increase, or viscosity decreases, the flow regime transitions from Stokes (Re < 0.1) through transitional (0.1 < Re < 1000) to Newton's regime (Re > 1000). In the intermediate range, drag coefficients deviate significantly from the Stokes prediction of C_d = 24/Re. The Schiller-Naumann correlation provides improved accuracy for 0.1 < Re < 1000: C_d = (24/Re)(1 + 0.15Re^0.687). At Re > 1000, the drag coefficient approaches a constant value near 0.44, characteristic of pressure drag dominating over viscous drag.

For engineering calculations outside the Stokes regime, iterative solutions become necessary. The general drag equation F_d = C_d·(ρ_fv_s²/2)·A requires solving implicitly since C_d depends on velocity through Reynolds number. Engineers typically employ iterative numerical methods or use charts correlating C_d versus Re for spheres.

Applications in Water and Wastewater Treatment

Sedimentation tanks in municipal water treatment plants rely fundamentally on Stokes settling to remove suspended solids. Primary clarifiers remove particles typically ranging from 10 to 100 μm in diameter. Design engineers calculate required surface area based on overflow rate (surface loading), which equals settling velocity for the target removal size. A typical design criterion might specify removal of 90% of particles with settling velocity exceeding 0.0004 m/s (1.44 m/h), translating to a maximum overflow rate of 1.44 m/h.

For a plant treating 10,000 m³/day, this yields a required surface area of approximately 289 m². The depth of the clarifier, typically 3-5 meters, must provide sufficient detention time for particles to reach the sludge zone before exiting with the effluent. The hydraulic residence time equals volume divided by flow rate, generally 1.5 to 2.5 hours for primary clarifiers.

Mineral Processing and Classification

Hydraulic classifiers, spiral classifiers, and thickeners in mining operations exploit differential settling velocities to separate particles by size and density. In a rising current classifier, upward water velocity is adjusted so particles with settling velocity less than the upflow rate are carried overhead (fines), while particles settling faster than the upflow collect as underflow (coarse fraction). This principle enables separation of silica (ρ ≈ 2650 kg/m³) from magnetite (ρ ≈ 5150 kg/m³) or separation of different size fractions of the same mineral.

For example, separating 50 μm quartz (ρ = 2650 kg/m³) from 20 μm magnetite (ρ = 5150 kg/m³) in water at 20°C (μ = 0.001 Pa·s): the quartz settling velocity is 1.13×10⁻³ m/s while magnetite settles at 5.59×10⁻⁴ m/s. Despite being smaller, the magnetite particles settle at roughly half the velocity of quartz due to lower density differential, but careful control of upflow velocity can achieve separation.

Atmospheric Particulate Matter and Air Quality

Air pollution control focuses heavily on particles below 10 μm diameter (PM10) and 2.5 μm (PM2.5) because these remain suspended for extended periods and penetrate deep into respiratory systems. A 10 μm dust particle with density 2000 kg/m³ settling through air (ρ_f = 1.2 kg/m³, μ = 1.8×10⁻⁵ Pa·s at 20°C) exhibits settling velocity of only 0.00306 m/s or 11 m/h. Such a particle released at 100 meters altitude requires approximately 9 hours to reach ground level in still air, explaining why fine dust can travel hundreds of kilometers from its source.

PM2.5 particles (2.5 μm) settle even more slowly, at approximately 0.00019 m/s or 0.69 m/h, taking over 144 hours (6 days) to settle from 100 meters. This extreme persistence in the atmosphere necessitates active filtration or electrostatic precipitation for removal rather than reliance on gravitational settling.

Worked Example: Sedimentation Tank Design

A municipal water treatment plant must design a primary sedimentation basin to remove 85% of suspended sand particles from 15,000 m³/day of river water. Laboratory analysis shows the sand particles have a median diameter of 78 μm and density of 2650 kg/m³. Water temperature averages 15°C (viscosity 1.139×10⁻³ Pa·s, density 999.1 kg/m³). Determine the required basin surface area, recommended depth, and verify the flow regime validity.

Step 1: Calculate settling velocity

Using Stokes' law with d = 78×10⁻⁶ m, ρ_p = 2650 kg/m³, ρ_f = 999.1 kg/m³, μ = 1.139×10⁻³ Pa·s, and g = 9.81 m/s²:

v_s = [9.81 × (78×10⁻⁶)² × (2650 - 999.1)] / [18 × 1.139×10⁻³]

v_s = [9.81 × 6.084×10⁻⁹ × 1650.9] / [0.02050]

v_s = [9.847×10⁻⁵] / [0.02050]

v_s = 4.803×10⁻³ m/s = 0.288 m/min = 17.29 m/h

Step 2: Verify Stokes regime validity

Calculate Reynolds number: Re = (ρ_f·v_s·d) / μ

Re = (999.1 × 4.803×10⁻³ × 78×10⁻⁶) / (1.139×10⁻³)

Re = (3.745×10⁻⁴) / (1.139×10⁻³) = 0.329

Since Re = 0.329, which is greater than 0.1 but less than 1, we are in the transitional regime where Stokes law provides approximate results. The actual settling velocity will be slightly lower than calculated, providing a conservative design margin.

Step 3: Determine required surface area

For 85% removal efficiency, the overflow rate (surface loading rate) should not exceed the settling velocity of the target particle. Using v_s = 17.29 m/h as the maximum overflow rate:

Surface Area = Flow Rate / Overflow Rate = (15,000 m³/day) / (17.29 m/h)

Converting flow rate to hourly: 15,000 / 24 = 625 m³/h

Surface Area = 625 / 17.29 = 36.15 m²

Applying a safety factor of 1.5 to account for flow distribution irregularities, short-circuiting, and the transitional Reynolds regime: Required Area = 36.15 × 1.5 = 54.2 m²

Step 4: Select basin dimensions and calculate detention time

Selecting rectangular basin dimensions of 9 m length × 6 m width gives 54 m² surface area. Choosing depth h = 3.5 m provides volume V = 54 × 3.5 = 189 m³.

Hydraulic detention time = Volume / Flow Rate = 189 m³ / 625 m³/h = 0.302 hours = 18.1 minutes

This detention time is somewhat low for a primary clarifier; typical designs range from 1.5 to 2.5 hours. To achieve 2 hours detention time: Required Volume = 625 × 2 = 1250 m³

Maintaining 54 m² surface area: Required Depth = 1250 / 54 = 23.1 m

This depth is impractical. Instead, increase surface area to achieve both appropriate overflow rate and detention time. For 2-hour detention with 3.5 m depth:

Required Surface Area = 1250 / 3.5 = 357 m²

This increases overflow rate to: 625 / 357 = 1.75 m/h

This overflow rate will capture particles with v_s ≥ 1.75 m/h. Recalculating minimum particle size captured:

1.75 m/h = 4.861×10⁻⁴ m/s

Solving Stokes equation for diameter: d = √[(18 × μ × v_s) / (g × (ρ_p - ρ_f))]

d = √[(18 × 1.139×10⁻³ × 4.861×10⁻⁴) / (9.81 × 1650.9)]

d = √[(9.967×10⁻⁶) / (16,193.3)] = √(6.155×10⁻¹⁰) = 2.481×10⁻⁵ m = 24.8 μm

With 357 m² surface area at 3.5 m depth, the basin will remove particles ≥ 24.8 μm with high efficiency. Since the target median diameter is 78 μm, this design provides substantial safety margin and should achieve well over 85% removal efficiency. Final recommended basin: 17 m × 21 m × 3.5 m deep = 357 m² × 3.5 m = 1250 m³ volume.

Hindered Settling and Particle Concentration Effects

The classical Stokes equation applies to isolated particles in dilute suspensions where particles do not interact. As particle concentration increases above approximately 0.1% by volume, particles begin to interfere with each other's settling paths, fluid displaced by settling particles creates upward flow that retards settling, and particle wakes influence nearby particles. This phenomenon, called hindered settling, can reduce settling velocity by 50% or more at high concentrations.

The Richardson-Zaki correlation provides empirical adjustment: v_hindered = v_Stokes(1 - φ)ⁿ, where φ is the volumetric particle concentration and n is an empirical exponent (typically 4.65 for Re < 0.2). For a suspension with 5% solids by volume (φ = 0.05), the hindered settling velocity becomes v_hindered = v_Stokes(0.95)^4.65 = 0.78v_Stokes, representing a 22% velocity reduction. Thickener and clarifier designs must account for this effect when dealing with concentrated slurries.

For additional engineering calculator resources covering fluid dynamics, material properties, and process design, visit the complete engineering calculators library.

Practical Applications

Frequently Asked Questions