Role of grain dynamics in determining the onset of sediment transport

Sediment transport occurs when the nondimensional fluid shear stress Θ at the bed surface exceeds a minimum value Θc. A large collection of data, known as the Shields curve, shows that Θc is primarily a function of the shear Reynolds number Re_*. It is commonly assumed that Θ>Θc(Re_*) occurs when the Re_*-dependent fluid forces are too large to maintain static equilibrium for a typical surface grain. A complimentary approach, which remains relatively unexplored, is to identify Θc(Re_*) as the applied shear stress at which grains cannot stop moving. With respect to grain dynamics, Re_* can be viewed as the viscous time scale for a grain to equilibrate to the fluid flow divided by the typical time for the fluid force to accelerate a grain over the characteristic bed roughness. We performed simulations of granular beds sheared by a model fluid, varying only these two time scales. We find that the critical Shields number Θc( Re_*) obtained from the model mimics the Shields curve and is insensitive to the grain properties, the model fluid flow, and the form of the drag law. Quantitative discrepancies between the model results and the Shields curve are consistent with previous calculations of lift forces at varying Re_*. Grains at low Re_* find more stable configurations than those at high Re_* due to differences in the grain reorganization dynamics. Thus, instead of focusing on mechanical equilibrium of a typical grain at the bed surface, Θc(Re_*) may be better described by the stress at which mobile grains cannot find a stable configuration and stop moving.