Vol. 14, No. 1, 2019

Download this article
Download this article For screen
For printing
Recent Issues
Volume 14, Issue 2
(in progress)
Volume 14, Issue 1
Volume 13, Issue 2
Volume 13, Issue 1
Volume 12, Issue 1
Volume 11, Issue 2
Volume 11, Issue 1
Volume 10, Issue 2
Volume 10, Issue 1
Volume 9, Issue 2
Volume 9, Issue 1
Volume 8, Issue 1
Volume 7, Issue 2
Volume 7, Issue 1
Volume 6, Issue 1
Volume 5, Issue 2
Volume 5, Issue 1
Volume 4, Issue 1
Volume 3, Issue 1
Volume 2, Issue 1
Volume 1, Issue 1
The Journal
About the Journal
Editorial Board
Subscriptions
Scientific Advantages
Submission Guidelines
Submission Form
Editorial Login
Ethics Statement
ISSN: 2157-5452 (e-only)
ISSN: 1559-3940 (print)
Author Index
To Appear
 
Other MSP Journals
This article is available for purchase or by subscription. See below.
Computation of volume potentials on structured grids with the method of local corrections

Chris Kavouklis and Phillip Colella

Vol. 14 (2019), No. 1, 1–32
Abstract

We present a new version of the method of local corrections (MLC) of McCorquodale, Colella, Balls, and Baden (2007), a multilevel, low-communication, noniterative domain decomposition algorithm for the numerical solution of the free space Poisson’s equation in three dimensions on locally structured grids. In this method, the field is computed as a linear superposition of local fields induced by charges on rectangular patches of size O(1) mesh points, with the global coupling represented by a coarse-grid solution using a right-hand side computed from the local solutions. In the present method, the local convolutions are further decomposed into a short-range contribution computed by convolution with the discrete Green’s function for a Q-th-order accurate finite difference approximation to the Laplacian with the full right-hand side on the patch, combined with a longer-range component that is the field induced by the terms up to order P 1 of the Legendre expansion of the charge over the patch. This leads to a method with a solution error that has an asymptotic bound of O(hP) + O(hQ) + O(ϵh2) + O(ϵ), where h is the mesh spacing and ϵ is the max norm of the charge times a rapidly decaying function of the radius of the support of the local solutions scaled by h. The bound O(ϵ) is essentially the error of the global potential computed on the coarsest grid in the hierarchy. Thus, we have eliminated the low-order accuracy of the original method (which corresponds to P = 1 in the present method) for smooth solutions, while keeping the computational cost per patch nearly the same as that of the original method. Specifically, in addition to the local solves of the original method we only have to compute and communicate the expansion coefficients of local expansions (that is, for instance, 20 scalars per patch for P = 4). Several numerical examples are presented to illustrate the new method and demonstrate its convergence properties.

PDF Access Denied

However, your active subscription may be available on Project Euclid at
https://projecteuclid.org/camcos

We have not been able to recognize your IP address 34.231.21.123 as that of a subscriber to this journal.
Online access to the content of recent issues is by subscription, or purchase of single articles.

Please contact your institution's librarian suggesting a subscription, for example by using our journal-recom­mendation form. Or, visit our subscription page for instructions on purchasing a subscription.

You may also contact us at contact@msp.org
or by using our contact form.

Or, you may purchase this single article for USD 40.00:

Keywords
Poisson solver, method of local corrections, Mehrstellen stencils, domain decomposition, parallel solvers
Mathematical Subject Classification 2010
Primary: 65N06, 65N12, 65N15, 68W10
Milestones
Received: 5 October 2016
Revised: 14 July 2018
Accepted: 16 July 2018
Published: 31 October 2018
Authors
Chris Kavouklis
Computational Engineering Division
Lawrence Livermore National Laboratory
Livermore, CA
United States
Phillip Colella
Applied Numerical Algorithms Group
Computational Research Division
Lawrence Berkeley National Laboratory
Berkeley, CA
United States