An implicit compressible SPH solver for snow simulation

Abstract

Snow is a complex material. It resists elastic normal and shear deformations, while some deformations are plastic. Snow can deform and break. It can be significantly compressed and gets harder under compression. Existing snow solvers produce impressive results. E.g., hybrid Lagrangian/Eulerian techniques have been used to capture all material properties of snow. The auxiliary grid, however, makes it challenging to handle small volumes. In particular, snow fall and accumulation on surfaces have not been demonstrated with these solvers yet. Existing particle-based snow solvers, on the other hand, can naturally handle small snow volumes. However, existing solutions consider simplified material properties. In particular, shear deformation and the hardening effect are typically omitted. We present a novel Lagrangian snow approach based on Smoothed Particle Hydrodynamics (SPH). Snow is modeled as an elastoplastic continuous material that captures all above-mentioned effects. The compression of snow is handled by a novel compressible pressure solver, where the typically employed state equation is replaced by an implicit formulation. Acceleration due to shear stress is computed using a second implicit formulation. The linear solvers of the two implicit formulations for accelerations due to shear and normal stress are realized with matrix-free implementations. Using implicit formulations and solving them with matrix-free solvers allows to couple the snow to other phases and is beneficial to the stability and the time step size, i.e., performance of the approach. Solid boundaries are represented with particles and a novel implicit formulation is used to handle friction at solid boundaries. We show that our approach can simulate accumulation, deformation, breaking, compression and hardening of snow. Furthermore, we demonstrate two-way coupling with rigid bodies, interaction with incompressible and highly viscous fluids and phase change from fluid to snow.

A Momentum-Conserving Implicit Material Point Method for Surface Energies with Spatial Gradients

“…Lastly, although our approach was designed for MPM, SPH is more commonly used for simulation of liquids. However, SPH and MPM have many similarities as recently shown by the work of Gissler et al [Gissler et al 2020b] and it would be interesting future work to generalize our approach to SPH.…”

A Survey on SPH Methods in Computer Graphics

“…The simulation of viscoelastic materials has been investigated by Takahashi et al [TDFN14] and Barreiro et al [BGFAO17]. Gissler et al [GHB*20] introduced an SPH formulation for elastoplastic material behavior in order to simulate snow. Yan et al [YJL*16] present a multi‐phase method for fluids and solids.…”

An Optimization‐based SPH Solver for Simulation of Hyperelastic Solids

“…On the other hand, to resolve such coupling, Dagenais et al [DGP12] used a predictor‐corrector approach with shape matching, and Huber et al [HEW15] presented a simulation method for cloth with SPH, which supports wetting and no‐slip boundary conditions. Gissler et al [GHB*20] proposed a compressible SPH pressure solver, which is coupled with a linear implicit elasticity solver, such that it can simulate plasticity.…”