Finite Element Numerical Simulation of Local Scour of a Three-Dimensional Cylinder under Steady Flow

Pile safety has received increasing attention in marine engineering, especially in the field of local scour. In this paper, a finite element numerical model is established for local scour around a cylinder in steady currents. The flow is described by unsteady Reynolds–averaged Navier–Stokes equations with a traditional k−ε turbulent closure model. The proposed scour model takes bed load into account. The Exner equation is solved to determine the bed variation and the moving mesh approach is used to capture the evolution of the bed. When the resulting slope exceeds the angle of repose, a novel sand-slide model based on Rodrigues' rotation formula is used to prevent simulation distortion. All the equations are discretized by the two-step Taylor–Galerkin algorithm, and the resulting approach is fast to implement with second-order accuracy in space. The numerical results are found to be in good agreement with the experimental data.