A Fast Algorithm for Solving the Poisson Equation on a Nested Grid

Abstract

We present a numerical method for solving the Poisson equation on a nested grid. The nested grid consists of uniform grids having different grid spacing and is designed to cover the space closer to the center with a finer grid. Thus our numerical method is suitable for computing the gravity of a centrally condensed object. It consists of two parts: the difference scheme for the Poisson equation on the nested grid and the multi-grid iteration algorithm. It has three advantages: accuracy, fast convergence, and scalability. First it computes the gravitational potential of a close binary accurately up to the quadraple moment, even when the binary is resolved only in the fine grids. Second residual decreases by a factor of 300 or more by each iteration. We confirmed experimentally that the iteration converges always to the exact solution of the difference equation. Third the computation load of the iteration is proportional to the total number of the cells in the nested grid. Thus our method gives a good solution at the minimum expense when the nested grid is large. The difference scheme is applicable also to the adaptive mesh refinement in which cells of different sizes are used to cover a domain of computation.

22 pages 21 figures. To appear in ApJ

Cited by in corpus (27)

• Radiation Magnetohydrodynamic Simulations of Protostellar Collapse: Protostellar Core Formation
• Radiation Magnetohydrodynamic Simulations of Protostellar Collapse: Non-Ideal Magnetohydrodynamic Effects and Early Formation of Circumstellar Disks
• Fragmentation of a Molecular Cloud Core versus Fragmentation of the Massive Protoplanetary Disk in the Main Accretion Phase
• Directions of Outflows, Disks, Magnetic Fields, and Rotation of YSOs in Collapsing Molecular Cloud Cores
• Conditions for Circumstellar Disk Formation: Effects of Initial Conditions and Sink Treatment
• A Direct Multigrid Poisson Solver for Oct-Tree Adaptive Meshes
• Collapse and Fragmentation of Rotating Magnetized Clouds. I. Magnetic Flux - Spin Relation
• Formation Scenario for Wide and Close Binary Systems
• First MHD Simulation of Collapse and Fragmentation of Magnetized Molecular Cloud Cores
• Self-gravitational Magnetohydrodynamics with Adaptive Mesh Refinement for Protostellar Collapse
• The First Two Thousand Years of Star Formation
• The Formation of Population III Binaries
• Origin of Misalignments: Protostellar Jet, Outflow, Circumstellar Disc, and Magnetic Field
• The Effect of Poloidal Magnetic Field on Type I Planetary Migration: Significance of Magnetic Resonance
• A Fast Poisson Solver of Second-Order Accuracy for Isolated Systems in Three-Dimensional Cartesian and Cylindrical Coordinates
• The Origin of the Rotation Profiles in Star Forming Clouds
• An Extension of the Athena++ Framework for Fully Conservative Self-Gravitating Hydrodynamics
• Hyperbolic Self-Gravity Solver for Large Scale Hydrodynamical Simulations
• Radiation Magnetohydrodynamic Simulations of Protostellar Collapse: Low-Metallicity Environments
• Self-gravity in curved mesh elements
• Self-gravitational force calculation of infinitesimally thin gaseous disks
• Self-Gravitational Force Calculation of Infinitesimally Thin Gaseous Disks on Nested Grids
• Differentiable Cosmological Hydrodynamics for Field-Level Inference and High Dimensional Parameter Constraints
• Computing the gravitational potential on nested meshes using the convolution method
• Exact, singularity-free recasting of the Newtonian potential in continuous media
• Self-gravity at the scale of the polar cell
• Synthetic Polarization Maps of an Outflow Zone from Magnetohydrodynamic Simulations