Gap-SPH

This solver programmed in C applies the Smoothed Particle Hydrodynamics method to subsonic incompressible 3d fluid flow. Viscosity is updated every half-time step to account for nonlinear behavior of one of the phases and all solid boundaries are modeled using the virtual particle method. A basic all-pair search algorithm is used to determine pair lists. The purpose of this self-contained solver is to further develop the code to optimize it to the case of river bedload dynamics. Parallel programming directives have been added where applicable. No GUI for pre- or post- processing is provided yet.