SPH solver

The SPH solver has been in development since about 2014. It the following features:

1. Fully incompressible flow – pressure Poisson equation solved using iterative solvers
2. Improved free surface capture
3. Surface tension due to pairwise force and Young-Laplace like pressure jump
4. Melting and solidifation phase change
5. Of course, heat transfer
6. Ability to read stl geometry to create domains

The code is hosted on bitbucket