Numerical modeling investigation on turbulent oscillatory flow over a plane rough bed composed by randomly arrayed particles

A three-dimensional numerical model is established to simulate the turbulent oscillatory boundary layer over a fixed and rough bed composed by randomly arrayed solid spheres based on the lattice Boltzmann method and the large eddy simulation model. The equivalent roughness height, the location of the theoretical bed and the time variation of the friction velocity are investigated using the log-fit method. The time series of turbulent intensity and Reynolds stress are also investigated. The equivalent roughness height of cases with Reynolds numbers of 1x10(4)-6x10(4) is approximately 2.81 d (grain size). The time variation of the friction velocity in an oscillatory cycle exhibits sinusoidal-like behavior. The friction factor depends on the relative roughness in the rough turbulent regime, and the pattern of solid particles arrayed as the rough bed in the numerical simulations has no obvious effect on the friction factor.