An ISPH with k-ε closure for simulating turbulence under solitary waves

In this paper an Incompressible Smoothed Particle Hydrodynamics (ISPH) method solving the 2D RANS (Reynolds Averaged Navier-Stokes) equations with the k-epsilon turbulence closure is constructed. In the present model, the concept of massless ISPH utilizing the definition of particle density (number of computational particles within unit volume) is stressed. The skills of this numerical model are tested by applying to two laboratory experiments: (1) A non-breaking solitary wave propagating over a bottom-mounted barrier and (2) a solitary wave breaking on a 1 on 50 slope. In the former case flow separation occurs behind the barrier as the wave crest passes by and a vortex is generated, which later interacts with free surface causing breaking. In the latter case wave breaking and bottom friction both generate significant turbulence. For both cases, the effects of initial seeding of turbulent kinetic energy, required in the k-epsilon model, are studied and it is concluded that initial values of O (10(-18)) to O (10(-8)) m(2)/s(2) should be used. An adaptive wall boundary condition for k-epsilon turbulence model is employed to avoid the unrealistic production of turbulence near the wall boundary. The numerical results, in terms of free surface profile, mean velocity field, vorticity field, turbulent kinetic energy and turbulent shear stress, are compared with experimental data. Very reasonable agreement is observed. This paper presents the first comprehensively validated 2D ISPH model with the k-epsilon turbulence closure, which can be applied to transient free surface wave problems.