Numerical and experimental analysis of shallow turbulent flow over complex roughness beds

A set of shallow-water equations (SWEs) based on a k−ε Reynold stress model is established to simulate the turbulent flows over a complex roughness bed. The fundamental equations are discretized by the second-order finite-difference method (FDM), in which spatial and temporal discretization are conducted by staggered-grid and leap-frog schemes, respectively. The turbulent model in this study stems from the standard k−ε model, but is enhanced by replacing the conventional vertical production with a more rigorous and precise generation derived from the energy spectrum and turbulence scales. To verify its effectiveness, the model is applied to compute the turbulence in complex flow surroundings (including a rough bed) in an abrupt bend and in a natural waterway. The comparison of the model results against experimental data and other numerical results shows the robustness and accuracy of the present model in describing hydrodynamic characteristics, especially turbulence features on the complex roughness bottom.