Discretization of Dirac delta functions in level set methods

Discretization of singular functions is an important component in many problems to which level set methods have been applied. We present two methods for constructing consistent approximations to Dirac delta measures concentrated on piecewise smooth curves or surfaces. Both methods are designed to be convenient for level set simulations and are introduced to replace the commonly used but inconsistent regularization technique that is solely based on a regularization parameter proportional to the mesh size. The first algorithm is based on a tensor product of regularized one-dimensional delta functions. It is independent of the irregularity relative to the grid. In the second method, the regularization is constructed from a one-dimensional regularization that is extended to multi-dimensions with a variable support depending on the orientation of the singularity relative to the computational grid. Convergence analysis and numerical results are given.