Smoothed Particle Hydrodynamics (SPH) is a meshless Lagrangian numerical method ideal for simulating potentially violent freesurface phenomena such as a wave breaking, or a dam break where many Eulerian methods can be difficult to apply.Dealing with wall boundary conditions is one of the most challenging parts of the SPH method and many different approaches have been developed among (i) repulsive forces such as LennardJones one, which is efficient to give impermeable boundaries but leads to nonphysical behaviours, (ii) fictitious (or ghost) particles which provide a better physical behaviour in the vicinity of a wall but are hard to define for complex geometries and (iii) semianalytical approach such as Kulasegaram et al. (2004) which consists of renormalising the density field near a solid wall with respect to the missing kernel support area. The present work extends this semianalytical methodology, where intrinsic gradient and divergence operators are employed that ensures conservation properties. The accuracy of the physical field such as the pressure next to walls is considerably improved, and the consistent manner developed to wallcorrect operators allows us to perform simulations with turbulence models. This work will present three key advances:• The time integration scheme used for the continuity equation requires particular attention, and as already mentioned by Vila (1999), we prove there is no point in using a dependence in time of the particles' density if no kernel gradient corrections are added. Thus, by using a nearboundary kernelcorrected version of the time integration scheme of the form proposed by Vila, we are able to simulate longtime simulations ideally suited for turbulent flow in a channel in the context of accurate boundary conditions.• In order to compute the kernel correction, Feldman and Bonet (2007) use an analytical value which is computationally expensive whereas Kulasegaram et al. (2004) and De Leffe et al. (2009) use polynomial approximation which can be difficult to define for complex geometries. We propose here to compute the renormalisation term of the kernel support near a solid with a novel time integration scheme, allowing us any shape for the boundary.• All boundary terms issued from the continuous approximation are given by surface summations which only require information from a mesh file of the boundary. The technique developed here allows us to correct the pressure gradient and viscous terms and hence provide a physically correct wallshear stress so that even the diffusion equation of a scalar quantity can be solved accurately using SPH such as the turbulent kinetic energy or its dissipation in a k  ε model of turbulence.The new model is demonstrated for cases including hydrostatic conditions for a channel flow, still water in a tank of complex geometry and a dam break over triangular bed profile with sharp angle where significant improved behaviour is obtained in comparison with the conventional boundary techniques. Simulation of the benchmark test case of a square object moving in an enclosed tank is shows good agreement with the reference solution and no voids are formed within the fluid domain. The performance of the model for a 2D turbulent flow in a channel is demonstrated where the profiles of velocity, turbulent kinetic energy and its dissipation are in agreement with the theoretical ones. Finally, the performance of the model is demonstrated for flow in a fishpass where velocity field and turbulent viscosity field are satisfactory reproduced compared to meshbased codes.
Date of Award  1 Aug 2011 

Original language  English 

Awarding Institution   The University of Manchester


Supervisor  Dominique Laurence (Supervisor) 

 Turbulence modelling in SPH
 SPH
 Wall boundary conditions
Unified semianalytical wall boundary conditions for inviscid, laminar and turbulent slightly compressible flows in SPARTACUS2D combined with an improved time integration scheme on the continuity equation.
Ferrand, M. (Author). 1 Aug 2011
Student thesis: Master of Philosophy