Integral Flow Modelling Approach for Surface Water-Groundwater Interactions along a Rippled Streambed

Exchange processes of surface and groundwater are important for the management of water quantity and quality as well as for the ecological functioning. In contrast to most numerical simulations using coupled models to investigate these processes, we present a novel integral formulation for the sediment-water-interface. The computational fluid dynamics (CFD) model OpenFOAM was used to solve an extended version of the three-dimensional Navier–Stokes equations which is also applicable in non-Darcy-flow layers. Simulations were conducted to determine the influence of ripple morphologies and surface hydraulics on the flow processes within the hyporheic zone for a sandy and for a gravel sediment. In- and outflowing exchange fluxes along a ripple were determined for each case. The results indicate that larger grain size diameters, as well as ripple distances, increased hyporheic exchange fluxes significantly. For higher ripple dimensions, no clear relationship to hyporheic exchange was found. Larger ripple lengths decreased the hyporheic exchange fluxes due to less turbulence between the ripples. For all cases with sand, non-Darcy-flow was observed at an upper layer of the ripple, whereas for gravel non-Darcy-flow was recognized nearly down to the bottom boundary. Moreover, the sediment grain sizes influenced also the surface water flow significantly.