Tree-based solvers for adaptive mesh refinement code FLASH - I: gravity and optical depths

Abstract

We describe an OctTree algorithm for the MPI-parallel, adaptive mesh-refinement code {\sc FLASH}, which can be used to calculate the gas self-gravity, and also the angle-averaged local optical depth, for treating ambient diffuse radiation. The algorithm communicates to the different processors only those parts of the tree that are needed to perform the tree walk locally. The advantage of this approach is a relatively low memory requirement, important in particular for the optical depth calculation, which needs to process information from many different directions. This feature also enables a general tree-based radiation transport algorithm that will be described in a subsequent paper, and delivers excellent scaling up to at least 1500 cores. Boundary conditions for gravity can be either isolated or periodic, and they can be specified in each direction independently, using a newly developed generalisation of the Ewald method. The gravity calculation can be accelerated with the {\em adaptive block update} technique by partially re-using the solution from the previous time-step. Comparison with the {\sc Flash} internal multi-grid gravity solver shows that tree based methods provide a competitive alternative, particularly for problems with isolated or mixed boundary conditions. We evaluate several multipole acceptance criteria (MACs) and identify a relatively simple APE MAC which provides high accuracy at low computational cost. The optical depth estimates are found to agree very well with those of the {\sc RADMC-3D} radiation transport code, with the tree solver being much faster. Our algorithm is available in the standard release of the {\sc FLASH} code in version 4.0 and later.

26 pages, 19 figures, 6 tables. Accepted by MNRAS

References in corpus (8)

• Modeling Collapse and Accretion in Turbulent Gas Clouds: Implementation and Comparison of Sink Particles in AMR and SPH
• Simulating the formation of molecular clouds. I. Slow formation by gravitational collapse from static initial conditions
• On the nature of star-forming filaments: I. Filament morphologies
• A Direct Multigrid Poisson Solver for Oct-Tree Adaptive Meshes
• The SILCC project --- IV. Impact of dissociating and ionising radiation on the interstellar medium and Halpha emission as a tracer of the star formation rate
• VINE -- A numerical code for simulating astrophysical systems using particles I: Description of the physics and the numerical methods
• GANDALF - Graphical Astrophysics code for N-body Dynamics And Lagrangian Fluids
• VINE -- A numerical code for simulating astrophysical systems using particles II: Implementation and performance characteristics

Cited by in corpus (7)

• Radiative cooling rates, ion fractions, molecule abundances and line emissivities including self-shielding and both local and metagalactic radiation fields
• The Astrochemical Impact of Cosmic Rays in Protoclusters I: Molecular Cloud Chemistry
• How fast do young star clusters expel their natal gas?: Estimating the upper limit of the gas expulsion time-scale
• The role of initial magnetic field structure in the launching of protostellar jets
• Dust charge distribution in the interstellar medium
• TREVR: A general $Nlo{g}^{2}N$ radiative transfer algorithm
• A Generalized Solution Method for Parallelized Computation of the Three-dimensional Gravitational Potential on a Multi-patch grid in Spherical Geometry