Tetrahedral grids in Monte Carlo radiative transfer
arXiv:2407.20216 · doi:10.1051/0004-6361/202450658
Abstract
Context. 3D numerical simulations of radiative transfer are crucial for understanding complex astrophysical objects. For Monte Carlo radiative transfer, the spatial grid design is critical yet complex. Common grids include hierarchical octree and unstructured Voronoi grids, each with its own strengths and weaknesses. Tetrahedral grids, widely used in ray-tracing graphics, are a potential alternative. Aims. We explore the possibilities, advantages, and limitations of tetrahedral grids for Monte Carlo radiative transfer, comparing their performance with other grid structures. Method. We integrated a tetrahedral grid structure, using the TetGen library, into the SKIRT Monte Carlo radiative transfer code. Tetrahedral grids can be imported or adaptively constructed and refined within SKIRT. We implemented an efficient grid traversal method using Plücker coordinates and Plücker products. Results. We validated the tetrahedral grid construction and traversal algorithm with 2D radiative transfer benchmarks. In a simple 3D model, we compared the performance of tetrahedral, octree, and Voronoi grids. The octree grid outperformed the others in traversal speed, while the tetrahedral grid had the lowest grid quality. Overall, tetrahedral grids performed worse than octree and Voronoi grids. Conclusion. While tetrahedral grids may not be ideal for most astrophysical simulations, they offer a viable unstructured alternative to Voronoi grids for specific applications, such as post-processing hydrodynamical simulations on tetrahedral or unstructured grids.
10 pages, 8 figures, abstract was shortened to fit the arxiv abstract requirements, article is accepted to Astronomy & Astrophysics (in production)
References in corpus (15)
- Monte Carlo radiative transfer in protoplanetary disks
- SKIRT: an Advanced Dust Radiative Transfer Code with a User-Friendly Architecture
- Monte Carlo Radiative Transfer
- IllustrisTNG in the HSC-SSP: image data release and the major role of mini mergers as drivers of asymmetry and star formation
- Rich: Open Source Hydrodynamic Simulation on a Moving Voronoi Mesh
- SOC program for dust continuum radiative transfer
- Morphological signatures of mergers in the TNG50 simulation and the Kilo-Degree Survey: the merger fraction from dwarfs to Milky Way-like galaxies
- SKIRT: hybrid parallelization of radiative transfer simulations
- The atomic-to-molecular hydrogen transition in the TNG50 simulation: Using realistic UV fields to create spatially resolved HI maps
- Self-consistent dust and non-LTE line radiative transfer with SKIRT
- Effects of Spatial Discretization in Lyman-alpha Line Radiation Transfer Simulations
- A Moving Mesh Hydrodynamic Solver for ChaNGa
- Smart detectors for Monte Carlo radiative transfer
- Atomic line radiative transfer with MCFOST I. Code description and benchmarking
- Monte Carlo radiative transfer with explicit absorption to simulate absorption, scattering, and stimulated emission