A Hybrid Particle Level Set Method for Improved Interface Capturing

Abstract

In this paper, we propose a new numerical method for improving the mass conservation properties of the level set method when the interface is passively advected in a flow field. Our method uses Lagrangian marker particles to rebuild the level set in regions which are under-resolved. This is often the case for flows undergoing stretching and tearing. The overall method maintains a smooth geometrical description of the interface and the implementation simplicity characteristic of the level set method. Our method compares favorably with volume of fluid methods in the conservation of mass and purely Lagrangian schemes for interface resolution. The method is presented in three spatial dimensions.

Level Set Method for Simulating the Interface Kinematics: Application of a Discontinuous Galerkin Method

“…Although this method is highly robust, but there is no guaranty for the area/volume conservation if the numerical method applied to solve the level set advection equation is not preciseness enough, see [5]. The particle level set method which was introduced in [16], implements a different way for maintaining the area/volume conservation. Each side of the interface in this method is assigned to a distinguished set of the massless particles which can be advected through the domain in a Lagrangian way.…”

Particle level set implementation on the finite volume method

“…In recent years an attempt has been made to improve the LS method in a way that it can be used for free surface prediction. Particle Level Set (PLS) is the result of recent researches in this field [14,15]. The first implementation of the PLS Scheme was on the Finite Difference method.…”

“…The step by step method of correction can be summarized as [14]: 1) Those particles which have distance more than their radius outside of their own fluid phase with free surface, considered as an escaped particle. These particles are used to correct the free surface position.…”

Fast Sparse Level Sets on Graphics Hardware

“…9. Performance comparison between our GPU-STL method and the DT-Grid [24] for the "Enright test" [8], using the re-initialization PDE with HJ-WENO numerical scheme in both cases.…”

“…The so-called "Enright test" [8] (pure advection) was also performed, to further investigate the performance of the GPU-STL method. The test was performed by both the DT-Grid method [24] and ours, on a grid of size 512 3 voxels, using the HJ-WENO scheme with identical time steps and an even larger narrow band, obtained by setting γ = 5.…”