An Efficient Boundary Handling with a Modified Density Calculation for SPH

Abstract

We propose a new boundary handling method for smoothed particle hydrodynamics (SPH). Previous approaches required the use of boundary particles to prevent particles from sticking to the boundary. We address this issue by correcting the fundamental equations of SPH with the integration of a kernel function. Our approach is able to directly handle triangle mesh boundaries without the need for boundary particles. We also show how our approach can be integrated into a position-based fluid framework.

Semi‐analytical Solid Boundary Conditions for Free Surface Flows

“…= 967 / 3 = 3079 / 3 = 7008 / 3 = 988 / 3 = 952 / 3 = 954 / 3 ℎ = 1.5 ( )ℎ = 2 ℎ = 2.5 Figure 10: A comparison between our method (bottom) and Fujisawa and Miura [FM15](top). By assuming the rest density to be 1000kg/m 3 , we calculate the average density for the one layer of particles that are closest to the boundary.…”

“…Compared to the one using a signed distance field, our method only introduces 26% additional computational cost while the simulation with ghost particles introduce 52%. In addition, Figure 10 shows a comparison between our method and the one proposed by Fujisawa and Miura [FM15] in calculating the density for particles near the boundary. Note the unsatisfactory results Fujisawa and Miura's method gets when a large kernel size is set for the density calculation.…”

“…To reduce the computational cost, semi‐analytical methods aim to convert the volume integral over the boundary into a surface integral. However, existing methods either cannot handle complex boundary geometries [KBLP04,FM15] or are only suited to take the boundary integral for the gradient term [FLR*13, MFK*15].…”

Angiography Simulation and Planning Using a Multi-Fluid Approach

“…The simulation of a SPH fluid always shows artifacts at the boundary. Therefore, many boundary handling techniques [26,27] for SPH-based fluid simulations have been proposed. The blood vessel is an elastic object and its deformation is influenced by the heart beat and movements of nearby tissues.…”

Density maps for improved SPH boundary handling

“…Similar gradient correction terms for weakly compressible fluids were introduced by Ferrand et al [2013] for two-dimensional and by Mayrhofer et al [2011] for three-dimensional domains. The renormalization approach was later adopted by Fujisawa and Miura [2015] in the field of computer graphics. They represent boundaries using triangular meshes and compute the renormalization factors during runtime.…”