DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: An adaptive fast multipole accelerated Poisson solver for complex geometries

Abstract

Here, we present a fast, direct and adaptive Poisson solver for complex two-dimensional geometries based on potential theory and fast multipole acceleration. More precisely, the solver relies on the standard decomposition of the solution as the sum of a volume integral to account for the source distribution and a layer potential to enforce the desired boundary condition. The volume integral is computed by applying the FMM on a square box that encloses the domain of interest. For the sake of efficiency and convergence acceleration, we first extend the source distribution (the right-hand side in the Poisson equation) to the enclosing box as a $$C$$0 function using a fast, boundary integral-based method. We demonstrate on multiply connected domains with irregular boundaries that this continuous extension leads to high accuracy without excessive adaptive refinement near the boundary and, as a result, to an extremely efficient “black box” fast solver.

Authors:
ORCiD logo [1]; ORCiD logo [2]
  1. Univ. of Washington, Seattle, WA (United States)
  2. New York Univ. (NYU), NY (United States)
Publication Date:
Research Org.:
New York Univ. (NYU), NY (United States)
Sponsoring Org.:
USDOE Office of Science (SC), Fusion Energy Sciences (FES); US Air Force Office of Scientific Research (AFOSR)
OSTI Identifier:
1533962
Alternate Identifier(s):
OSTI ID: 1398114
Grant/Contract Number:  
FG02-86ER53223; FG02-88ER25053; SC0012398; FA9550-10-1-0180; FA9550-15-1-0385; DEFG0288ER25053
Resource Type:
Accepted Manuscript
Journal Name:
Journal of Computational Physics
Additional Journal Information:
Journal Volume: 344; Journal ID: ISSN 0021-9991
Publisher:
Elsevier
Country of Publication:
United States
Language:
English
Subject:
97 MATHEMATICS AND COMPUTING; Poisson equation; Fast multipole method; Quadrature by expansion; Integral equations

Citation Formats

Askham, T., and Cerfon, A. J. An adaptive fast multipole accelerated Poisson solver for complex geometries. United States: N. p., 2017. Web. doi:10.1016/j.jcp.2017.04.063.
Askham, T., & Cerfon, A. J. An adaptive fast multipole accelerated Poisson solver for complex geometries. United States. https://doi.org/10.1016/j.jcp.2017.04.063
Askham, T., and Cerfon, A. J. Tue . "An adaptive fast multipole accelerated Poisson solver for complex geometries". United States. https://doi.org/10.1016/j.jcp.2017.04.063. https://www.osti.gov/servlets/purl/1533962.
@article{osti_1533962,
title = {An adaptive fast multipole accelerated Poisson solver for complex geometries},
author = {Askham, T. and Cerfon, A. J.},
abstractNote = {Here, we present a fast, direct and adaptive Poisson solver for complex two-dimensional geometries based on potential theory and fast multipole acceleration. More precisely, the solver relies on the standard decomposition of the solution as the sum of a volume integral to account for the source distribution and a layer potential to enforce the desired boundary condition. The volume integral is computed by applying the FMM on a square box that encloses the domain of interest. For the sake of efficiency and convergence acceleration, we first extend the source distribution (the right-hand side in the Poisson equation) to the enclosing box as a $C$0 function using a fast, boundary integral-based method. We demonstrate on multiply connected domains with irregular boundaries that this continuous extension leads to high accuracy without excessive adaptive refinement near the boundary and, as a result, to an extremely efficient “black box” fast solver.},
doi = {10.1016/j.jcp.2017.04.063},
journal = {Journal of Computational Physics},
number = ,
volume = 344,
place = {United States},
year = {Tue May 02 00:00:00 EDT 2017},
month = {Tue May 02 00:00:00 EDT 2017}
}

Journal Article:

Citation Metrics:
Cited by: 22 works
Citation information provided by
Web of Science

Save / Share:

Works referenced in this record:

Nonlinear, three-dimensional magnetohydrodynamics of noncircular tokamaks
journal, January 1976


A Fast Poisson Solver for Complex Geometries
journal, May 1995

  • McKenney, A.; Greengard, L.; Mayo, A.
  • Journal of Computational Physics, Vol. 118, Issue 2
  • DOI: 10.1006/jcph.1995.1104

ECOM: A fast and accurate solver for toroidal axisymmetric MHD equilibria
journal, May 2015


Fast Parallel Iterative Solution of Poisson’s and the Biharmonic Equations on Irregular Regions
journal, January 1992

  • Mayo, A.; Greenbaum, A.
  • SIAM Journal on Scientific and Statistical Computing, Vol. 13, Issue 1
  • DOI: 10.1137/0913006

Laplace's Equation and the Dirichlet-Neumann Map in Multiply Connected Domains
journal, April 1993

  • Greenbaum, A.; Greengard, L.; McFadden, G. B.
  • Journal of Computational Physics, Vol. 105, Issue 2
  • DOI: 10.1006/jcph.1993.1073

A Fast Adaptive Multipole Algorithm for Particle Simulations
journal, July 1988

  • Carrier, J.; Greengard, L.; Rokhlin, V.
  • SIAM Journal on Scientific and Statistical Computing, Vol. 9, Issue 4
  • DOI: 10.1137/0909044

A Fast Direct Solver for Structured Linear Systems by Recursive Skeletonization
journal, January 2012

  • Ho, Kenneth L.; Greengard, Leslie
  • SIAM Journal on Scientific Computing, Vol. 34, Issue 5
  • DOI: 10.1137/120866683

High-Order Corrected Trapezoidal Quadrature Rules for Singular Functions
journal, August 1997


A Direct Adaptive Poisson Solver of Arbitrary Order Accuracy
journal, May 1996

  • Greengard, Leslie; Lee, June-Yub
  • Journal of Computational Physics, Vol. 125, Issue 2
  • DOI: 10.1006/jcph.1996.0103

On the evaluation of layer potentials close to their sources
journal, February 2008


On the Convergence of Local Expansions of Layer Potentials
journal, January 2013

  • Epstein, Charles L.; Greengard, Leslie; Klöckner, Andreas
  • SIAM Journal on Numerical Analysis, Vol. 51, Issue 5
  • DOI: 10.1137/120902859

An Implementation of the Fast Multipole Method without Multipoles
journal, July 1992

  • Anderson, Christopher R.
  • SIAM Journal on Scientific and Statistical Computing, Vol. 13, Issue 4
  • DOI: 10.1137/0913055

Hybrid Gauss-Trapezoidal Quadrature Rules
journal, January 1999


Remarks on algorithm 006: An adaptive algorithm for numerical integration over an N-dimensional rectangular region
journal, December 1980


Quadrature by expansion: A new method for the evaluation of layer potentials
journal, November 2013

  • Klöckner, Andreas; Barnett, Alexander; Greengard, Leslie
  • Journal of Computational Physics, Vol. 252
  • DOI: 10.1016/j.jcp.2013.06.027

Laplace’s equation and the Dirichlet–Neumann map: a new mode for Mikhlin’s method
journal, January 2005


A fast, high-order solver for the Grad–Shafranov equation
journal, June 2013

  • Pataki, Andras; Cerfon, Antoine J.; Freidberg, Jeffrey P.
  • Journal of Computational Physics, Vol. 243
  • DOI: 10.1016/j.jcp.2013.02.045

Universal quadratures for boundary integral equations on two-dimensional domains with corners
journal, November 2010

  • Bremer, James; Rokhlin, Vladimir; Sammis, Ian
  • Journal of Computational Physics, Vol. 229, Issue 22
  • DOI: 10.1016/j.jcp.2010.06.040

An Improved Fast Multipole Algorithm for Potential Fields
journal, November 1998


On the solution of elliptic partial differential equations on regions with corners
journal, January 2016


Immersed boundary smooth extension: A high-order method for solving PDE on arbitrary smooth domains using Fourier spectral methods
journal, January 2016


Corner singularities for elliptic problems: Integral equations, graded meshes, quadrature, and compressed inverse preconditioning
journal, October 2008


A Robust and Accurate Solver of Laplace's Equation with General Boundary Conditions on General Domains in the Plane
journal, June 2012


A Nonlinear Optimization Procedure for Generalized Gaussian Quadratures
journal, January 2010

  • Bremer, James; Gimbutas, Zydrunas; Rokhlin, Vladimir
  • SIAM Journal on Scientific Computing, Vol. 32, Issue 4
  • DOI: 10.1137/080737046

Spectrally Accurate Quadratures for Evaluation of Layer Potentials Close to the Boundary for the 2D Stokes and Laplace Equations
journal, January 2015

  • Barnett, Alex; Wu, Bowei; Veerapaneni, Shravan
  • SIAM Journal on Scientific Computing, Vol. 37, Issue 4
  • DOI: 10.1137/140990826

GMRES: A Generalized Minimal Residual Algorithm for Solving Nonsymmetric Linear Systems
journal, July 1986

  • Saad, Youcef; Schultz, Martin H.
  • SIAM Journal on Scientific and Statistical Computing, Vol. 7, Issue 3
  • DOI: 10.1137/0907058

Gerris: a tree-based adaptive solver for the incompressible Euler equations in complex geometries
journal, September 2003


On Direct Methods for Solving Poisson’s Equations
journal, December 1970

  • Buzbee, B. L.; Golub, G. H.; Nielson, C. W.
  • SIAM Journal on Numerical Analysis, Vol. 7, Issue 4
  • DOI: 10.1137/0707049

Approximating real functions which possess nth derivatives of bounded variation and applications
journal, November 2008


A free-space adaptive FMM-Based PDE solver in three dimensions
journal, January 2011

  • Langston, Harper; Greengard, Leslie; Zorin, Denis
  • Communications in Applied Mathematics and Computational Science, Vol. 6, Issue 1
  • DOI: 10.2140/camcos.2011.6.79

Evaluation of Layer Potentials Close to the Boundary for Laplace and Helmholtz Problems on Analytic Planar Domains
journal, January 2014

  • Barnett, Alex H.
  • SIAM Journal on Scientific Computing, Vol. 36, Issue 2
  • DOI: 10.1137/120900253

The Fast Solution of Poisson’s and the Biharmonic Equations on Irregular Regions
journal, April 1984

  • Mayo, Anita
  • SIAM Journal on Numerical Analysis, Vol. 21, Issue 2
  • DOI: 10.1137/0721021

Accurate derivative evaluation for any Grad–Shafranov solver
journal, January 2016


The use of the ghost fluid method for Poisson's equation to simulate streamer propagation in point-to-plane and point-to-point geometries
journal, February 2009

  • Celestin, Sebastien; Bonaventura, Zdenek; Zeghondy, Barbar
  • Journal of Physics D: Applied Physics, Vol. 42, Issue 6
  • DOI: 10.1088/0022-3727/42/6/065203

Fast convolution with free-space Green's functions
journal, October 2016

  • Vico, Felipe; Greengard, Leslie; Ferrando, Miguel
  • Journal of Computational Physics, Vol. 323
  • DOI: 10.1016/j.jcp.2016.07.028

Approximating real functions which possess nth derivatives of bounded variation and applications
journal, November 2008


Laplace’s equation and the Dirichlet–Neumann map: a new mode for Mikhlin’s method
journal, January 2005


The use of the ghost fluid method for Poisson's equation to simulate streamer propagation in point-to-plane and point-to-point geometries
journal, February 2009

  • Celestin, Sebastien; Bonaventura, Zdenek; Zeghondy, Barbar
  • Journal of Physics D: Applied Physics, Vol. 42, Issue 6
  • DOI: 10.1088/0022-3727/42/6/065203

A New Fast-Multipole Accelerated Poisson Solver in Two Dimensions
journal, January 2001


Works referencing / citing this record:

Integral Equation Formulation of the Biharmonic Dirichlet Problem
journal, September 2017


Adaptive quadrature by expansion for layer potential evaluation in two dimensions
text, January 2017


Distributed and Adaptive Fast Multipole Method In Three Dimensions
preprint, January 2020