In this paper, we present a fourth-order algorithm to solve Poisson’s equation in two
and three dimensions. We use a Cartesian grid, embedded boundary method to
resolve complex boundaries. We use a weighted least squares algorithm to solve for
our stencils. We use convergence tests to demonstrate accuracy and we show the
eigenvalues of the operator to demonstrate stability. We compare accuracy and
performance with an established second-order algorithm. We also discuss in depth
strategies for retaining higher-order accuracy in the presence of nonsmooth
geometries.
Keywords
Poisson equation, finite volume methods, high order,
embedded boundary
Applied Numerical Algorithms Group
(ANAG)
Computational Research Division
Lawrence Berkeley National Laboratory
MS 50A1148
One Cyclotron Road
Berkeley, CA 94720
United States