Skip to main content

Shields Diagram Sediment Incipient Motion & Critical Shear Stress Calculator

Predicting whether riverbed sediments will remain stable or mobilize under flowing water is essential for bridge pier scour protection, channel design, and alluvial sediment transport modeling.

Median sieve diameter of bed sediments (e.g. 0.2mm fine sand, 2.0mm gravel).

Specific gravity of mineral particles (typically 2.65 for quartz sand).

Cross-sectional area divided by wetted perimeter R = A / P.

Channel bed slope or energy slope S (e.g. 0.002 = 2 m per km).

Ambient water temperature affecting kinematic viscosity ν.

Calculated Result
23.54 Pa

Bed Shear Stress (τ₀)

Critical Threshold (τ_cr)

1.29 Pa

Shields Parameter (θ)

0.727

Critical Shields (θ_cr)

0.04

Bed Motion State

Intense Bedload & Suspended Sediment Transport (Sheet flow regime)

Calculation Breakdown

  1. Boundary Bed Shear Stress (τ₀)τ₀ = γ_w · R_h · S₀ = 9810 · 1.2 · 0.002 = 23.54 Pa
  2. Shields Dimensionless Entrainment Parameter (θ)θ = τ₀ / [(ρ_s - ρ_w)·g·d] = 0.727; Critical Threshold θ_cr = 0.04
  3. Critical Threshold Shear Stress (τ_cr)τ_cr = θ_cr · (Δρ·g·d) = 1.29 Pa (Mobility Ratio τ₀/τ_cr = 18.25)

Sediment Shear Stress Comparison

Interactive visualization based on your current inputs

Value
0.05.9121824Bed Shear τ0Critical Shear τcrMobility RatioParameterValue (Pa)

What Is the Shields Diagram Sediment Incipient Motion & Critical Shear Stress Calculator?

Predicting whether riverbed sediments will remain stable or mobilize under flowing water is essential for bridge pier scour protection, channel design, and alluvial sediment transport modeling.

The Shields diagram relates the dimensionless boundary shear stress (Shields parameter θ) to the boundary Reynolds number Re* or dimensionless particle diameter D*.

This calculator utilizes the Soulsby-van Rijn continuous analytical approximation of the Shields curve to evaluate critical shear stress τ_cr and bed mobility status.

How Does the Shields Diagram Sediment Incipient Motion & Critical Shear Stress Calculator Work?

The calculation evaluates user-provided measurements using recognized domain equations, converts between measurement units, and adjusts for real-world efficiency factors.

Shields Diagram Sediment Incipient Motion & Critical Shear Stress Calculator Formula & Variables

The core mathematical equation utilized by this calculator is expressed as:

\tau_0 = \rho g R S, \quad \theta = \frac{\tau_0}{(\rho_s - \rho) g d_{50}}, \quad \theta_{cr} = \frac{0.30}{1 + 1.2 D_*} + 0.055 \left(1 - e^{-0.02 D_*}\right)

Soulsby continuous analytical formulation of the empirical Shields curve for incipient sediment motion.

How to Use the Shields Diagram Sediment Incipient Motion & Critical Shear Stress Calculator

  1. Enter your primary measurements in the input fields above.
  2. Select your preferred units (e.g. metric or imperial) if applicable.
  3. Review or adjust operational assumptions such as field efficiency.
  4. Click Calculate to instantly generate the full results breakdown and visual chart.
  5. Use the Reset button at any time to clear the form and test a new scenario.

Step-by-Step Example Calculation

Gravel Bed Stream (d50 = 2.0mm, R = 1.2m, Slope = 0.002)

Input Values:

medianGrainDiameterD50Mm:2.0
sedimentSpecificGravityGs:2.65
channelHydraulicRadiusM:1.2
energyBedSlopeMPerM:0.002
waterTemperatureCelsius:20
Worked Steps: Bed shear τ₀ = 23.54 Pa exceeds critical shear τ_cr = 1.15 Pa (Mobility ratio = 20.5), indicating intense active sediment motion.

Understanding Your Result

Your calculated result represents the realistic operational capacity or baseline output under the specified conditions. Comparing theoretical and effective outputs reveals the direct impact of turns, overlap, and practical downtime.

Factors That Affect the Result

Field terrain, operator experience, equipment maintenance, overlap margin, and weather conditions can significantly influence real-world output.

When Should You Use This Calculator?

Use this calculator whenever you need quick, verified estimates for job planning, budgeting, equipment sizing, or project timelines.

Assumptions & Limitations

  • Actual Bed Shear Stress: τ₀ = ρ_w · g · R · S_0 (Pa).
  • Dimensionless Grain Size: D* = d₅₀ · [(s - 1)·g / ν²]^(1/3).
  • Critical Shields Parameter: θ_cr = 0.30 / (1 + 1.2·D*) + 0.055 · [1 - exp(-0.02·D*)].
  • Critical Shear Stress: τ_cr = θ_cr · (ρ_s - ρ_w) · g · d₅₀.
  • Mobility Ratio: M = τ₀ / τ_cr. When M > 1.0, bed motion is active.

Frequently Asked Questions

Calculation Accuracy & Reference Note

This calculator implements verified, deterministic mathematical equations based on published standards. Results should be treated as professional engineering estimates; always verify critical operations with local equipment manuals and site inspections.

Explore more tools and calculators in Math Calculators