Magritte, a modern software library for 3D radiative transfer: II. Adaptive ray-tracing, mesh construction and reduction
Authors:
Frederik De Ceuster,
Jan Bolte,
Ward Homan,
Silke Maes,
Jolien Malfait,
Leen Decin,
Jeremy Yates,
Peter Boyle,
James Hetherington
Abstract:
Radiative transfer is a notoriously difficult and computationally demanding problem. Yet, it is an indispensable ingredient in nearly all astrophysical and cosmological simulations. Choosing an appropriate discretization scheme is a crucial part of the simulation, since it not only determines the direct memory cost of the model but also largely determines the computational cost and the achievable…
▽ More
Radiative transfer is a notoriously difficult and computationally demanding problem. Yet, it is an indispensable ingredient in nearly all astrophysical and cosmological simulations. Choosing an appropriate discretization scheme is a crucial part of the simulation, since it not only determines the direct memory cost of the model but also largely determines the computational cost and the achievable accuracy. In this paper, we show how an appropriate choice of directional discretization scheme as well as spatial model mesh can help alleviate the computational cost, while largely retaining the accuracy. First, we discuss the adaptive ray-tracing scheme implemented in our 3D radiative transfer library Magritte, that adapts the rays to the spatial mesh and uses a hierarchical directional discretization based on HEALPix. Second, we demonstrate how the free and open-source software library Gmsh can be used to generate high quality meshes that can be easily tailored for Magritte. In particular, we show how the local element size distribution of the mesh can be used to optimise the sampling of both analytically and numerically defined models. Furthermore, we show that when using the output of hydrodynamics simulations as input for a radiative transfer simulation, the number of elements in the input model can often be reduced by an order of magnitude, without significant loss of accuracy in the radiation field. We demonstrate this for two models based on a hierarchical octree mesh resulting from adaptive mesh refinement (AMR), as well as two models based on smoothed-particle hydrodynamics (SPH) data.
△ Less
Submitted 30 November, 2020;
originally announced November 2020.
Magritte, a modern software library for 3D radiative transfer: I. Non-LTE atomic and molecular line modelling
Authors:
Frederik De Ceuster,
Ward Homan,
Jeremy Yates,
Leen Decin,
Peter Boyle,
James Hetherington
Abstract:
Radiative transfer is a key component in almost all astrophysical and cosmological simulations. We present Magritte: a modern open-source software library for 3D radiative transfer. It uses a deterministic ray-tracer and formal solver, i.e. it computes the radiation field by tracing rays through the model and solving the radiative transfer equation in its second-order form along a fixed set of ray…
▽ More
Radiative transfer is a key component in almost all astrophysical and cosmological simulations. We present Magritte: a modern open-source software library for 3D radiative transfer. It uses a deterministic ray-tracer and formal solver, i.e. it computes the radiation field by tracing rays through the model and solving the radiative transfer equation in its second-order form along a fixed set of rays originating from each point. Magritte can handle structured and unstructured input meshes, as well as smoothed-particle hydrodynamics (SPH) particle data. In this first paper, we describe the numerical implementation, semi-analytic tests and cross-code benchmarks for the non-LTE line radiative transfer module of Magritte. This module uses the radiative transfer solver to self-consistently determine the populations of the quantised energy levels of atoms and molecules using an accelerated Lambda iteration (ALI) scheme. We compare Magritte with the established radiative transfer solvers Ratran (1D) and Lime (3D) on the van Zadelhoff benchmark and present a first application to a simple Keplerian disc model. Comparing with Lime, we conclude that Magritte produces more accurate and more precise results, especially at high optical depth, and that it is faster.
△ Less
Submitted 18 December, 2019;
originally announced December 2019.