miluphcuda
SPH code for the simulation of fluids and solid bodies using CUDA(TM).
For a description, please read:
Schäfer, C.; Riecker, S.; Maindl, T.; Speith, R.; Scherrer, S.; Kley, W.: A Smooth Particle Hydrodynamics Code to Model Collisions Between Solid, Self-Gravitating Objects, doi link.
github site [miluphcuda is open source]

some examples

Little boat on water column, colour code shows absolute value of velocity

Colliding rubber rings

Rotating rubber stick, two densities

Pile of sand simulation. Factor a is the ratio between initial height and initial radius of the pile.

Pile of sand under gravity, Drucker-Prager yield criterion, a = 0.55

Pile of sand under gravity, Mohr-Coulomb yield criterion, a = 0.55

Pile of sand under gravity, Drucker-Prager yield criterion, a = 2.75

Pile of sand under gravity, Mohr-Coulomb yield criterion, a = 2.75

Collision between two Ceres-sized (mass) objects, blue are ice/water particles, yellow are basalt particles
The radii are 437 km and 509 km, the collision speed is 630 m/s. The model includes: von Mises yield strength, tensile
damage model (following Benz & Asphaug 1995), self-gravity and Tillotson equation of state. Total simulation time is
24 hours. The final outcome is merging.