Lattice Green's Functions for High-Order Finite Difference Stencils

James Gabbard,Wim M. van Rees
DOI: https://doi.org/10.1137/23m1573872
IF: 3.039
2024-01-06
SIAM Journal on Numerical Analysis
Abstract:SIAM Journal on Numerical Analysis, Volume 62, Issue 1, Page 25-47, February 2024. Lattice Green's functions (LGFs) are fundamental solutions to discretized linear operators, and as such they are a useful tool for solving discretized elliptic PDEs on domains that are unbounded in one or more directions. The majority of existing numerical solvers that make use of LGFs rely on a second-order discretization and operate on domains with free-space boundary conditions in all directions. Under these conditions, fast expansion methods are available that enable precomputation of 2-dimensional or 3-dimensional (3D) LGFs in linear time, avoiding the need for brute-force multidimensional quadrature of numerically unstable integrals. Here we focus on higher-order discretizations of the Laplace operator on domains with more general boundary conditions, by (1) providing an algorithm for fast and accurate evaluation of the LGFs associated with high-order dimension-split centered finite differences on unbounded domains, and (2) deriving closed-form expressions for the LGFs associated with both dimension-split and Mehrstellen discretizations on domains with one unbounded dimension. Through numerical experiments we demonstrate that these techniques provide LGF evaluations with near machine-precision accuracy, and that the resulting LGFs allow for numerically consistent solutions to high-order discretizations of the Poisson's equation on fully or partially unbounded 3D domains.
mathematics, applied
What problem does this paper attempt to address?