Numerical Modeling of Fluid Interface Phenomena

In [A.-K. Tornberg, B. Engquist, Numerical approximations of singular source terms in differential equations, J. Comput. Phys. 200 (2004) 462-488], it was shown for simple examples that the then most common way to regularize delta functions in connection to level set methods produces inconsistent approximations with errors that are not reduced with grid refinement. Since then, several clever approximations have been derived to overcome this problem. However, the great appeal of the old method was its simplicity. In this paper it is shown that the old method - a one-dimensional delta function approximation extended to higher dimensions by a distance function - can be made accurate with a different class of one-dimensional delta function approximations. The prize to pay is a wider support of the resulting delta function approximations.