Prescribed Velocity Gradients for Highly Viscous SPH Fluids with Vorticity Diffusion

Abstract

Working with prescribed velocity gradients is a promising approach to efficiently and robustly simulate highly viscous SPH fluids. Such approaches allow to explicitly and independently process shear rate, spin, and expansion rate. This can be used to, e.g., avoid interferences between pressure and viscosity solvers. Another interesting aspect is the possibility to explicitly process the vorticity, e.g., to preserve the vorticity. In this context, this paper proposes a novel variant of the prescribed-gradient idea that handles vorticity in a physically motivated way. In contrast to a less appropriate vorticity preservation that has been used in a previous approach, vorticity is diffused. The paper illustrates the utility of the vorticity diffusion. Therefore, comparisons of the proposed vorticity diffusion with vorticity preservation and additionally with vorticity damping are presented. The paper further discusses the relation between prescribed velocity gradients and prescribed velocity Laplacians which improves the intuition behind the prescribed-gradient method for highly viscous SPH fluids. Finally, the paper discusses the relation of the proposed method to a physically correct implicit viscosity formulation.

An SPH model to simulate the dynamic behavior of shear thickening fluids

“…Over the years, SPH has been established as one major approach for fluid animation in computer graphics, being applied to different applications such as for simulating foam or bubbles. 29 Several extensions have been proposed, for example, exploring vorticity diffusion with vorticity preservation for viscous fluids, 30 achieving efficient divergence-free velocity fields, 31 and addressing viscous SPH fluids that can generate coiling and buckling phenomena. 32 In addition to SPH models, position-based methods have also become popular in the computer graphics community because of the improved speed, stability, and control offered by the approach, even if the results may not be as accurate as force-based methods.…”

Incompressible Projective Smooth Particle Hydrodynamic

“…The open source library is a fluid simulation engine based on OpenGL and C++. The library contains five popular SPH fluid simulation methods (DFSPH [2], IISPH [6], PBF [8], PCISPH [5], PF [1]) and three surface tension processing methods (Becker 2007 [9], Akinci 2013 [10], He2014 [11]) and six viscosity term processing methods (standard [12], XSPH [13] [15]) and two methods of curl processing (MicropolarModel_Bender2017 [16], VorticityConfinement [19]) and two types of DragForce (Macklin 2014 [17], Gissler 2017 [18]). We use a parallel hash algorithm to speed up the domain search of particles [20], and the density calculation uses the cubic spline kernel function [21].…”

Conformation constraints for efficient viscoelastic fluid simulation

“…Even though their method computes the velocity eld very e ciently, we propose an even more e cient and versatile solution, which handles larger time steps thanks to a PBD-style constrained dynamics solver, and supports more diverse viscoelastic behaviors. Peer et al [2017] have later extended their method to improve vorticity handling.…”