Two-Dimensional Simulation Model of Sediment Removal and Flow in Rectangular Sedimentation Basin

A steady, two-dimensional numerical model was created to study the hydrodynamics of a rectangular sedimentation basin under turbulent conditions. The strip integral method was used to formulate the flow equations, using a forward marching scheme for solving the governing partial differential equations of continuity, momentum, advection–diffusion, turbulent kinetic energy, and its dissipation. In this way the flow equations were converted to a set of ordinary differential equations (ODEs) in terms of the key physical parameters. These parameters, along with a set of shape functions, describe flow variables including the velocity, the concentration of suspended sediments, and both the kinetic energy and its dissipation rate. Four Gaussian distributions were investigated, one corresponding to each flow parameter. In order to calculate the turbulent shear stresses, a two-equation turbulence model (i.e., k-ε model) was used. A fourth order Runge–Kutta method numerically integrates the set of ODEs. Simulation results were compared with experimental data, and close agreement (generally within 5–10%) was observed.