Accelerating Surface Tension Calculation in SPH via Particle Classification and Monte Carlo Integration

Abstract

Surface tension has a strong influence on the shape of fluid interfaces. We propose a method to calculate the corresponding forces efficiently. In contrast to several previous approaches, we discriminate to this end between surface and non-surface SPH particles. Our method effectively smooths the fluid interface, minimizing its curvature. We make use of an approach inspired by Monte Carlo integration to estimate local normals as well as curvatures, based on which the force can be calculated. We compare different sampling schemes for the Monte Carlo approach, for which a Halton sequence performed best. Our overall technique is applicable, but not limited to 2D and 3D simulations, and can be coupled with any common SPH formulation. It outperforms prior approaches with regard to total computation time per time step in dynamic scenes. Additionally, it is adjustable for higher quality in small scale scenes with dominant surface tension effects.

Keywords

• SPH fluid sim.
• particle classification
• estimation of normal, curvature, surface tension