Abstract:The Immersed Boundary (IB) method of Peskin (J. Comput. Phys., 1977) is useful for problems that involve fluid-structure interactions or complex geometries. By making use of a regular Cartesian grid that is independent of the geometry, the IB framework yields a robust numerical scheme that can efficiently handle immersed deformable structures. Additionally, the IB method has been adapted to problems with prescribed motion and other PDEs with given boundary data. IB methods for these problems traditionally involve penalty forces which 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 step in time, this method can be prohibitively inefficient without preconditioning. This work introduces a new, well-conditioned IB formulation for boundary value problems, called the Immersed Boundary Double Layer (IBDL) method. In order to lay the groundwork for similar formulations of Stokes and Navier-Stokes equations, this paper focuses on the Poisson and Helmholtz equations to introduce the methodology and to demonstrate its efficiency over the original constraint method. In this double layer 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 without preconditioning. Furthermore, the iteration count is independent of both the mesh size and spacing of the immersed boundary points. 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.

A Projection Preconditioner for Solving the Implicit Immersed Boundary Equations

An Efficient Semi-Implicit Immersed Boundary Method for the Navier-Stokes Equations

Iterative projection method for unsteady Navier-Stokes equations with high Reynolds numbers

On Preconditioned Iterative Methods for Unsteady Incompressible Navier-Stokes Equations

Some Investigations on Applications of Point-implicit Runge-Kutta Method

An implicit lattice Boltzmann flux solver with a projection-based interpolation scheme for the convection-diffusion equation

An efficient preconditioning method and its application to numerical simulation of steady and unsteady viscous flows

Analysis of some projection method based preconditioners for models of incompressible flow

Efficient Variable-Coefficient Finite-Volume Stokes Solvers

Evaluation of preconditioned Krylov subspace methods for Navier-Stokes equations

An Immersed Boundary Projection Method for Complex Fluid-Structure-Interaction Simulation

A Second-Order Decoupled Implicit/Explicit Method of the 3d Primitive Equations of Ocean Ii: Finite Element Spatial Discretization

An Immersed Boundary Solver For Inviscid Compressible Flows

A Nonlinearly Preconditioned Inexact Newton Algorithm for Steady State Lattice Boltzmann Equations.

Semi-implicit method of high-index saddle dynamics and application to construct solution landscape

Immersed Boundary Double Layer method: An introduction of methodology on the Helmholtz equation

SIMPLE-like Preconditioners for Saddle Point Problems from the Steady Navier-Stokes Equations

An Efficient Immersed Boundary Projection Method for Flow over Complex/moving Boundaries

Preconditioning Technique for A Simultaneous Solution to Wind-Membrane Interaction

Implicit interpolation method for immersed boundary methods

Preconditioning immersed isogeometric finite element methods with application to flow problems