Solving Elliptic Interface Problems with Jump Conditions on Cartesian Grids

Abstract

We present a simple numerical algorithm for solving elliptic equations where the diffusion coefficient, the source term, the solution and its flux are discontinuous across an irregular interface. The algorithm produces second-order accurate solutions and first-order accurate gradients in the $L^{\infty }$-norm on Cartesian grids. The condition number is bounded, regardless of the ratio of the diffusion constant and scales like that of the standard 5-point stencil approximation on a rectangular grid with no interface. Numerical examples are given in two and three spatial dimensions.

16 pages, 9 figures, submitted to Journal of Computational Physics

Cited by in corpus (7)

• A cusp-capturing PINN for elliptic interface problems
• A multigrid ghost-point level-set method for incompressible Navier-Stokes equations on moving domains with curved boundaries
• A Shallow Ritz Method for Elliptic Problems with Singular Sources
• A High Order Cartesian Grid, Finite Volume Method for Elliptic Interface Problems
• A Numerical Method for Sharp-Interface Simulations of Multicomponent Alloy Solidification
• A Compact Coupling Interface Method with Accurate Gradient Approximation for Elliptic Interface Problems
• Immersed Boundary Method for the Complete Electrode Model in Electrical Impedance Tomography