Distributed and shared memory parallelism with a smoothed particle hydrodynamics code

Abstract

The Smoothed Particle Hydrodynamics (sph) method, being Lagrangian in nature, provides advantages in modelling flows containing interfaces. The sph method is completely mesh-free, modelling the fluid as a collection of N particles which move with the fluid velocity. The continuum fluid properties at a particular location are interpolated as weighted sums of the properties of surrounding particles in a process known as kernel interpolation. In its simplest form the sph method requires that every particle is used in the updating of every other particle. This leads to a solution