Repository logo
Log in(current)
Repository logoMIT Open ScholarshipDSpace@MIT
  1. Home
  2. MIT Open Access Articles
  3. MIT Open Access Articles
  4. Lattice Green’s Functions for High-Order Finite Difference Stencils

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

Thumbnail Image
Download
Name

2309.13503v1.pdf

Description
Accepted version
Size

522.42 KB

Format

Adobe PDF

Checksum (MD5)

fccb1348d87956561b6520823c42f3e0

sword-2024-05-09T21:02:06.original.xml (130 B)
Original SWORD entry document
Author(s)
Gabbard, James
•
van Rees, Wim M.
Date Issued
January 4, 2024
Journal
SIAM Journal on Numerical Analysis
Publisher
Society for Industrial & Applied Mathematics (SIAM)
Citation
Gabbard, James and van Rees, Wim M. 2024. "Lattice Green’s Functions for High-Order Finite Difference Stencils." SIAM Journal on Numerical Analysis, 62 (1).
Version
Author's final manuscript
Abstract
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 2D or 3D LGFs in linear time, avoiding the need for brute-force multi-dimensional 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.
MIT Department
Massachusetts Institute of Technology. Department of Mechanical Engineering
Terms of Use
Creative Commons Attribution-Noncommercial-ShareAlike
http://creativecommons.org/licenses/by-nc-sa/4.0/
Persistent DSpace Link
https://hdl.handle.net/1721.1/154893
DOI of Published Version
https://doi.org/10.1137/23m1573872
Repository logo
PrivacyPermissionsAccessibilityContact us
Repository logo
Notify us about copyright concerns.