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 (36)

• The cosmological simulation code GADGET-2
• HEALPix -- a Framework for High Resolution Discretization, and Fast Analysis of Data Distributed on the Sphere
• Dust Grain Size Distributions and Extinction in the Milky Way, LMC, and SMC
• E pur si muove: Galiliean-invariant cosmological hydrodynamical simulations on a moving mesh
• GADGET: A code for collisionless and gasdynamical cosmological simulations
• Structure of Stationary Photodissociation Fronts
• Gasoline: An adaptable implementation of TreeSPH
• Modeling Collapse and Accretion in Turbulent Gas Clouds: Implementation and Comparison of Sink Particles in AMR and SPH
• Tracing the Mass during Low-Mass Star Formation. II. Modelling the Submillimeter Emission from Pre-Protostellar Cores
• The SILCC (SImulating the LifeCycle of molecular Clouds) project: I. Chemical evolution of the supernova-driven ISM
• Simulating the formation of molecular clouds. I. Slow formation by gravitational collapse from static initial conditions
• A Hierarchical O(N) Force Calculation Algorithm
• The SILCC (SImulating the LifeCycle of molecular Clouds) project - II. Dynamical evolution of the supernova-driven ISM and the launching of outflows
• A Very Fast and Momentum-Conserving Tree Code
• The SILCC project: III. Regulation of star formation and outflows by stellar winds and supernovae
• On the nature of star-forming filaments: I. Filament morphologies
• TreeCol: a novel approach to estimating column densities in astrophysical simulations
• A Parallel Tree Code
• A Direct Multigrid Poisson Solver for Oct-Tree Adaptive Meshes
• How long does it take to form a molecular cloud?
• A New Parallel N-body Gravity Solver: TPM
• 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
• TreePM: A code for Cosmological N-Body Simulations
• SEREN - A new SPH code for star and planet formation simulations
• Tree-Particle-Mesh: an adaptive, efficient, and parallel code for collisionless cosmological simulation
• The Tree-Particle-Mesh N-body Gravity Solver
• VINE -- A numerical code for simulating astrophysical systems using particles I: Description of the physics and the numerical methods
• On column density thresholds and the star formation rate
• GRAPESPH with Fully Periodic Boundary Conditions: Fragmentation of Molecular Clouds
• 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
• EvoL: The new Padova T-SPH parallel code for cosmological simulations - I. Basic code: gravity and hydrodynamics
• A GPU accelerated Barnes-Hut Tree Code for FLASH4
• The Adaptive TreePM: An Adaptive Resolution Code for Cosmological N-body Simulations
• A Modified TreePM Code
• FLY: MPI-2 High Resolution code for LSS Cosmological Simulations

Cited by in corpus (64)

• Simulating cosmic structure formation with the GADGET-4 code
• Cooler and smoother -- the impact of cosmic rays on the phase structure of galactic outflows
• Impact of Cosmic Ray Transport on Galactic Winds
• SILCC VI -- Multi-phase ISM structure, stellar clustering, and outflows with supernovae, stellar winds, ionising radiation and cosmic rays
• Radiative cooling rates, ion fractions, molecule abundances and line emissivities including self-shielding and both local and metagalactic radiation fields
• The importance of magnetic fields for the initial mass function of the first stars
• From parallel to perpendicular -- On the orientation of magnetic fields in molecular clouds
• SILCC-Zoom: H${}_2$ and CO-dark gas in molecular clouds -- The impact of feedback and magnetic fields
• The SILCC project - V. The impact of magnetic fields on the chemistry and the formation of molecular clouds
• SILCC-ZOOM: The early impact of ionizing radiation on forming molecular clouds
• Role of Cosmic Ray Streaming and Turbulent Damping in Driving Galactic Winds
• SILCC VII -- Gas kinematics and multiphase outflows of the simulated ISM at high gas surface densities
• Photochemistry and Heating/Cooling of the Multiphase Interstellar Medium with UV Radiative Transfer for Magnetohydrodynamic Simulations
• A purely hyperbolic discontinuous Galerkin approach for self-gravitating gas dynamics
• The Astrochemical Impact of Cosmic Rays in Protoclusters I: Molecular Cloud Chemistry
• SILCC-Zoom: Polarisation and depolarisation in molecular clouds
• Is molecular cloud turbulence driven by external supernova explosions?
• How fast do young star clusters expel their natal gas?: Estimating the upper limit of the gas expulsion time-scale
• Synthetic observations of star formation and the interstellar medium
• Synthetic [CII] emission maps of a simulated molecular cloud in formation
• On the resolution requirements for modelling molecular gas formation in solar neighbourhood conditions
• Non-Equilibrium Chemistry and Destruction of CO by X-ray Flares
• The in situ formation of molecular and warm ionised gas triggered by hot outflows
• The Athena++ Adaptive Mesh Refinement Framework: Multigrid Solvers for Self-Gravity
• The role of initial magnetic field structure in the launching of protostellar jets
• On the accuracy of HI observations in molecular clouds -- More cold HI than thought?
• Tree-based solvers for adaptive mesh refinement code FLASH -- II: radiation transport module TreeRay
• Dust charge distribution in the interstellar medium
• Alignment of the magnetic field in star forming regions and why it might be difficult to observe
• The role of the H${}_2$ adiabatic index in the formation of the first stars
• The imprint of cosmic ray driven outflow on Lyman-alpha spectra
• On the extraction of power-law parts of the probability density functions in star-forming clouds
• The role of turbulence during the formation of circumbinary discs
• A new method for measuring the 3D turbulent velocity dispersion of molecular clouds
• Radiation transport methods in star formation simulations
• The dependence of episodic accretion on eccentricity during the formation of binary stars
• Catastrophic Cooling in Superwinds. III. Non-equilibrium Photoionization
• Unravelling the structure of magnetised molecular clouds with SILCC-Zoom: sheets, filaments and fragmentation
• Pyroclastic Blowout: Dust Survival in Isolated versus Clustered Supernovae
• Magnetic fields in star-forming systems (II): examining dust polarization, the Zeeman effect, and the Faraday rotation measure as magnetic field tracers
• TREVR: A general $Nlo{g}^{2}N$ radiative transfer algorithm
• On estimating angular momenta of infalling protostellar cores from observations
• The evolution of HCO${}^{+}$ in molecular clouds using a novel chemical post-processing algorithm
• Synthetic CO emission and the $X_{\mathrm{CO}}$ factor of young molecular clouds: a convergence study
• How do supernova remnants cool? -- I. Morphology, optical emission lines, and shocks
• Cloud properties across spatial scales in simulations of the interstellar medium
• Cosmic-ray-induced H${}_2$ line emission: Astrochemical modeling and implications for JWST observations
• Tree-based solvers for adaptive mesh refinement code FLASH -- III: a novel scheme for radiation pressure on dust and gas and radiative transfer from diffuse sources
• Dust Grain Growth & Dusty Supernovae in Low-Metallicity Molecular Clouds
• Population III star formation in the presence of turbulence, magnetic fields and ionizing radiation feedback
• The impact of cosmic ray heating on the cooling of the low-metallicity interstellar medium
• Simulating the evolution of optically dark HI clouds in the Virgo cluster : will no-one rid me of this turbulent sphere ?
• A Discontinuous Galerkin Solver in the FLASH Multi-Physics Framework
• Protostellar spin-up and fast rotator formation through binary star formation
• Gravity or turbulence? VII. The Schmidt-Kennicutt law, the star formation efficiency, and the mass density of clusters from gravitational collapse rather than turbulent support
• Rotation measure and synchrotron emission signatures in simulations of magnetized galactic discs
• The origin and evolution of the [CII] deficit in HII regions and star-forming molecular clouds
• Spectral reconstruction for radiation hydrodynamic simulations of galaxy evolution
• Tree-based solvers for adaptive mesh refinement code FLASH -- IV: An X-ray radiation scheme to couple discrete and diffuse X-ray emission sources to the thermochemistry of the interstellar medium
• A Generalized Solution Method for Parallelized Computation of the Three-dimensional Gravitational Potential on a Multi-patch grid in Spherical Geometry
• SILCC -- IX. The multi-phase interstellar medium at low metallicity
• Effect of cosmic rays and ionizing radiation on observational ultraviolet plasma diagnostics in the circumgalactic medium
• A resolution criterion based on characteristic timescales for MHD simulations of molecular clouds
• Linking simulations and observations of star-forming regions I. Masses and temperatures of ALMAGAL-like cores