Analysis of a one-dimensional model for the immersed boundary method

TLDR: In this paper, the authors studied numerical methods for the one-dimensional heat equation with a singular forcing term, where the delta function was replaced by a discrete approximation, and the resulting equation was solved by a Crank-Nicolson method on a uniform grid.

Abstract: Numerical methods are studied for the one-dimensional heat equation with a singular forcing term, $u_t = u_{xx} + c(t)\delta (x - \alpha (t)).$ The delta function $\delta (x)$ is replaced by a discrete approximation $d_h (x)$ and the resulting equation is solved by a Crank–Nicolson method on a uniform grid. The accuracy of this method is analyzed for various choices of $d_h $. The case where $c(t)$ is specified and also the case where c is determined implicitly by a constraint on the solution at the point a are studied. These problems serve as a model for the immersed boundary method of Peskin for incompressible flow problems in irregular regions. Some insight is gained into the accuracy that can be achieved and the importance of choosing appropriate discrete delta functions.