The Immersed Boundary Double Layer (IBDL) Method

TL;DR

The paper introduces the Immersed Boundary Double Layer (IBDL) method, a new well-conditioned formulation for boundary value problems that improves efficiency and convergence over traditional IB methods, applicable to various PDEs and boundary conditions.

Contribution

It proposes the IBDL method, a novel formulation that results in a well-conditioned second-kind integral equation, enhancing computational efficiency for boundary value problems.

Findings

Efficient solution with fewer Krylov iterations.

Independence of iteration count from mesh size.

Applicable to Dirichlet and Neumann boundary conditions.

Abstract

The Immersed Boundary (IB) method of Peskin (J. Comput. Phys., 1977) is useful for problems involving fluid-structure interactions or complex geometries. By making use of a regular grid that is independent of the geometry, the IB framework yields a robust numerical scheme that can efficiently handle immersed deformable structures. The IB method has also been adapted to problems with prescribed motion and other PDEs with given boundary data. IB methods for these problems traditionally involve penalty forces that only approximately satisfy boundary conditions, or they are formulated as constraint problems. In the latter approach, one must find the unknown forces by solving an equation that corresponds to a poorly conditioned first-kind integral equation. This operation can therefore require a large number of iterations of a Krylov method, and since a time-dependent problem requires this solve at each time step, this method can be prohibitively inefficient without preconditioning. This dissertation introduces a new, well-conditioned IB formulation for boundary value problems, called the Immersed Boundary Double Layer (IBDL) method. We formulate it for Poisson, Helmholtz, Brinkman, Stokes, and Navier-Stokes equations and demonstrate its efficiency over the other constraint method. In this new formulation, the equation for the unknown boundary distribution corresponds to a well-conditioned second-kind integral equation that can be solved efficiently with a small number of iterations of a Krylov method. Furthermore, the iteration count is independent of both the mesh size and the boundary point spacing. The method converges away from the boundary, and when combined with a local interpolation, it converges in the entire PDE domain. Additionally, while the original constraint method applies only to Dirichlet problems, the IBDL formulation can also be used for Neumann boundary conditions.