Hydrodynamics of coupled flow above and below a sediment–water interface with triangular bedforms

The hydrodynamics of a system where there is a coupled flow above and below a sediment–water interface (SWI) are not completely understood. We numerically simulate mean two-dimensional, unidirectional, steady, viscous flow in these systems using a sequentially coupled formulation. Simulations were conducted to determine fundamental relationships between bedform geometry, Reynolds number for the water-column flow (Re), interfacial exchange zone depth (d_z) in the sediments, and flux through the SWI (q_int); the latter two parameters play a significant role in biogeochemical and aquatic-life processes across the SWI. d_z and Re are functionally related through an asymptotic growth-curve model while q_int and Re follow a power function. These relationships are dynamically explained by the manner in which pressure gradients along the SWI develop due to current–bedform interactions at different Res and by Darcy’s Law. We found that the coupling between water column and exchange zone flow is controlled by the behavior of the water-column eddy. The eddy detaches at or near the point of minimum pressure along the interface, and reattaches near the point of maximum pressure. These two critical points determine the pressure gradient along the bed surface that controls the exchange zone flow field. Moreover, the reattachment point corresponds to flow divides within the sediments. Lastly, pore-water velocities drop with depth below the SWI, and are larger below the bedform crests than below the troughs.