A Direct Multigrid Poisson Solver for Oct-Tree Adaptive Meshes
Abstract
We describe a finite-volume method for solving the Poisson equation on oct-tree adaptive meshes using direct solvers for individual mesh blocks. The method is a modified version of the method presented by Huang and Greengard, which works with finite-difference meshes and does not allow for shared boundaries between refined patches. Our algorithm is implemented within the FLASH code framework and makes use of the PARAMESH library, permitting efficient use of parallel computers. We describe the algorithm and present test results that demonstrate its accuracy.
- Publication:
-
The Astrophysical Journal Supplement Series
- Pub Date:
- May 2008
- DOI:
- 10.1086/526425
- arXiv:
- arXiv:0710.4397
- Bibcode:
- 2008ApJS..176..293R
- Keywords:
-
- gravitation;
- methods: numerical;
- Astrophysics
- E-Print:
- 10 pages, 6 figures, accepted by the Astrophysical Journal