wumpf / yasph2d

some 2d sph fluid experimentation using rust

YaSPH2D

"Yet another SPH" .. thing.. in 2D!
Playing around with Smoothed Particle Hydrodynamics simulation in Rust. All 2d for now.

Implements solvers using

• Weakly Compressible SPH (WCSPH)
• DFSPH
  • Bender & Koschier 2015, Divergence-Free Smoothed Particle Hydrodynamicss
  • Bender & Koschier 2017, Divergence-Free SPH for Incompressible and Viscous Fluids

Nearest neighbor search using ideas from Compressed Neighbour Lists for SPH, Stefan Band et al.. Actual compression is WIP (see #3)

Some more links to resources in the code.

To find even more resources about fluid simulation in general check out my gist on CFD where I continously gather links and short descriptions on various concepts.

yasph2d's Issues

Polish WCSPH

• Fix hack around equation of state
• ?

Implement DFSPH

• pressure solver
• divergence solver solver
• warmup mechanism
• adaptive time steps
• optimizations..?

implement Interlinked SPH Pressure Solvers for Strong Fluid-Rigid Coupling

https://animation.rwth-aachen.de/publication/0563/

towards floating ducks! (without needing to implement a separate rigid body solver :P)

record and compute in parallel

writing images will always be slow, gotta pipeline this a bit

[perf] try kernel lookup tables

as described in https://animation.rwth-aachen.de/media/papers/2015-SCA-DFSPH.pdf chapter 3.4

Investigate more on how to do threading right in this case

There's reason to believe rayon is not the right tool for the parallism in yasph. Work loads are not dynamic but very predictable - i.e. static scheduling without any workstealing etc. should be perfectly fine, see also:
https://www.reddit.com/r/rust/comments/bto10h/update_a_scaling_comparison_between_rustrayon_and/
another fun topic to experiment on. lot to be learned here :)

profile using Microprofile

It's out for rust now! :)
https://github.com/jonasmr/microprofile-rust

Try compressed neighbour lists

Compressed Neighbour Lists for SPH, Stefan Band et al.