A Finite Difference Domain Decomposition Method Using Local Corrections for the Solution of Poisson's Equation
β Scribed by Gregory T. Balls; Phillip Colella
- Publisher
- Elsevier Science
- Year
- 2002
- Tongue
- English
- Weight
- 349 KB
- Volume
- 180
- Category
- Article
- ISSN
- 0021-9991
No coin nor oath required. For personal study only.
β¦ Synopsis
We present a domain decomposition method for computing finite difference solutions to the Poisson equation with infinite domain boundary conditions. Our method is a finite difference analogue of Anderson's Method of Local Corrections. The solution is computed in three steps. First, fine-grid solutions are computed in parallel using infinite domain boundary conditions on each subdomain. Second, information is transferred globally through a coarse-grid representation of the charge, and a global coarse-grid solution is found. Third, a fine-grid solution is computed on each subdomain using boundary conditions set with the global coarse solution, corrected locally with fine-grid information from nearby subdomains. There are three important features of our algorithm. First, our method requires only a single iteration between the local fine-grid solutions and the global coarse representation. Second, the error introduced by the domain decomposition is small relative to the solution error obtained in a single-grid calculation. Third, the computed solution is second-order accurate and only weakly dependent on the coarse-grid spacing and the number of subdomains. As a result of these features, we are able to compute accurate solutions in parallel with a much smaller ratio of communication to computation than more traditional domain decomposition methods. We present results to verify the overall accuracy, confirm the small communication costs, and demonstrate the parallel scalability of the method.
π SIMILAR VOLUMES
A domain decomposition method (DDM) is presented for the solution of the timeharmonic electromagnetic scattering problem by inhomogeneous 3-D objects. The computational domain is partitioned into concentric subdomains on the interfaces of which Robin-type transmission conditions are prescribed. On t
This paper describes a numerical method for the solution of a system of plasma fluid equations. The fluid model is similar to those employed in the simulation of high-density, low-pressure plasmas used in semiconductor processing. The governing equations consist of a drift-diffusion model of the ele
## Abstract A hybrid method for solution of Maxwell's equations of electromagnetics in the frequency domain is developed as a combination between the method of moments and the approximation in physical optics. The equations are discretized by a Galerkin method and solved by an iterative block Gauss