Hi. It seems that everyone these days make the sand solvers. This is the demonstration of a custom sand solver I have developed for Melnitsa Animation. It is erosion based solver. Mostly implemented with vex. Part of collision algorithm is writen in c++. There are also vex and c++ implementation of erosion nodes. Vex is faster but c++ is more accurate.
It can't produce small nice details as sideFx pbd solver in houdini 14, but it is faster and works in certan situations.
Pros: - Fast - Stable - Produce predictable results Cons: - not physically correct - not all types of concave collision objects are supported - it is not all-purpose solver