A High-Order Kernel-Free Boundary Integral Method for the Biharmonic Equation on Irregular Domains

• Published in: Journal of scientific computing Vol. 80; no. 3; pp. 1681 - 1699
• Main Authors: Xie, Yaning, Ying, Wenjun, Wang, Wei-Cheng
• Format: Journal Article
• Language: English
• Published: New York Springer US 01.09.2019; Springer Nature B.V
• Subjects: Accuracy; Algorithms; Biharmonic equations; Boundary conditions; Boundary integral method; Boundary value problems; Cartesian coordinates; Computational Mathematics and Numerical Analysis; Decomposition; Fast Fourier transformations; Finite difference method; Fluid mechanics; Fourier transforms; Free boundaries; Green's functions; Integral equations; Mathematical analysis; Mathematical and Computational Engineering; Mathematical and Computational Physics; Mathematics; Mathematics and Statistics; Methods; Numerical analysis; Operators (mathematics); Poisson equation; Theoretical; Biharmonic equation; Boundary integral method; Interface problem; Cartesian grid method; Nine-point compact finite difference scheme
• ISSN: 0885-7474; 1573-7691
• DOI: 10.1007/s10915-019-01000-6

This work proposes second-order and fourth-order versions of a Cartesian grid based kernel-free boundary integral (KFBI) method for the biharmonic equation on both bounded irregular domains and singly periodic irregular domains. It is further development of the previous KFBI method for second-order elliptic PDEs. It reformulates boundary value problems of the fourth-order PDE as boundary integral equations of the first kind but the solution never needs to know the fundamental solution or Green’s function of the elliptic operator. Evaluation of boundary or volume integrals in the solution of boundary integral equations is made by solving equivalent interface problems on Cartesian grids with standard finite difference methods and fast Fourier transform based solvers. The work decomposes the biharmonic equation into two Poisson equations. It assumes the solution to one Poisson equation, which has no boundary conditions, as the sum of a volume integral with a double layer boundary integral, and applies Green’s third identity to derive a scalar boundary integral equation from the other Poisson equation that are subject to two boundary conditions. In the solution of the scalar boundary integral equation, each volume or boundary integral is evaluated with the KFBI method. Numerical examples are presented to demonstrate the solution accuracy and algorithm efficiency. A remarkable point of the work is that the nine-point compact difference scheme in dealing with each split second-order elliptic interface problem on irregular domains yields fourth-order accurate solution for the biharmonic equation.