Discretization of Dirac delta functions in level set methods

  • Authors:

  • Björn Engquist;Anna-Karin Tornberg;Richard Tsai

  • Affiliations:

  • Department of Mathematics and PACM, Princeton University, Fine Hall, Princeton, NJ 08544-1000, USA;Courant Institute of Mathematical Sciences, New York University, 251 Mercer Street, New York, 10012-1185, USA;Institute for Advanced Study and Department of Mathematics, Princeton University. Fine Hall, Princeton, 08544-1000, USA

  • Venue:
  • Journal of Computational Physics
  • Year:
  • 2005

Quantified Score

Hi-index 31.59

Visualization

Abstract

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.