In this paper, a three-dimensional two-phase flow solver, SedFoam-2.0, is presented for sediment transport applications. The solver is extended from twoPhaseEulerFoam available in the 2.1.0 release of the open-source CFD (computational fluid dynamics) toolbox OpenFOAM. In this approach the sediment phase is modeled as a continuum, and constitutive laws have to be prescribed for the sediment stresses. In the proposed solver, two different intergranular stress models are implemented: the kinetic theory of granular flows and the dense granular flow rheology μ(I). For the fluid stress, laminar or turbulent flow regimes can be simulated and three different turbulence models are available for sediment transport: a simple mixing length model (one-dimensional configuration only), a k − ε, and a k − ω model. The numerical implementation is demonstrated on four test cases: sedimentation of suspended particles, laminar bed load, sheet flow, and scour at an apron. These test cases illustrate the capabilities of SedFoam-2.0 to deal with complex turbulent sediment transport problems with different combinations of intergranular stress and turbulence models.

The mobile layer of a granular bed composed of spherical particles is experimentally investigated in a laminar rectangular channel flow. Both particle and fluid velocity profiles are obtained using particle image velocimetry for different index-matched combinations of particles and fluid and for a wide range of fluid flow rates above incipient motion. A full three-dimensional investigation of the flow field inside the mobile layer is also provided. These experimental observations are compared to the predictions of a two-phase continuum model having a frictional rheology to describe particle–particle interactions. Different rheological constitutive laws having increasing degrees of sophistication are tested and discussed.

The local granular rheology is investigated numerically in turbulent bedload transport. Considering spherical particles, steady uniform configurations are simulated using a coupled fluid–discrete-element model. The stress tensor is computed as a function of the depth for a series of simulations varying the Shields number, the specific density and the particle diameter. The results are analysed in the framework of the $\unicode[STIX]{x1D707}(I)$ rheology and exhibit a collapse of both the shear to normal stress ratio and the solid volume fraction over a wide range of inertial numbers. Contrary to expectations, the effect of the interstitial fluid on the granular rheology is shown to be negligible, supporting recent work suggesting the absence of a clear transition between the free-fall and turbulent regimes. In addition, data collapse is observed up to unexpectedly high inertial numbers $I\sim 2$ , challenging the existing conceptions and parametrisation of the $\unicode[STIX]{x1D707}(I)$ rheology. Focusing upon bedload transport modelling, the results are pragmatically analysed in the $\unicode[STIX]{x1D707}(I)$ framework in order to propose a granular rheology for bedload transport. The proposed rheology is tested using a 1D volume-averaged two-phase continuous model, and is shown to accurately reproduce the dense granular flow profiles and the sediment transport rate over a wide range of Shields numbers. The present contribution represents a step in the upscaling process from particle-scale simulations towards large-scale applications involving complex flow geometry.

A new dataset of uniform and steady sheet-flow experiments is presented in this paper. An acoustic concentration and velocity profiler (ACVP) is used to measure time-resolved profiles of collocated 2C velocity ( $u,w$ ) and sediment concentration and to measure the time evolution of the bed interface position. Ensemble averaging over 11 similar experiment realisations is done to evaluate the mean profiles of streamwise velocity, concentration, sediment flux and Reynolds shear stress. The repeatability, stationarity and uniformity of the flow are carefully checked for a Shields number ${\it\theta}\approx 0.5$ and a suspension number of $S=1.1$ . The mean profile analysis allows to separate the flow into two distinct layers: a suspension layer dominated by turbulence and a bed layer dominated by granular interactions. The bed layer can be further subdivided into a frictional layer capped by a collisional layer. In the suspension layer, the mixing length profile is linear with a strongly reduced von Karman parameter equal to 0.225. The Schmidt number is found to be constant in this region with a mean value of ${\it\sigma}_{s}=0.44$ . The present results are then interpreted in terms of existing modelling approaches and the underlying assumptions are discussed. In particular, the well-known Rouse profile is shown to predict the concentration profile adequately in the suspension layer provided that all the required parameters can be evaluated separately. However, the strong intermittency of the flow observed in the bed layer under the impact of turbulent large-scale coherent flow structures suggests the limitations of averaged steady two-phase flow models.