Tetrahedral grids in Monte Carlo radiative transfer

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.

arXiv comment10 pages, 8 figures, abstract was shortened to fit the arxiv abstract requirements, article is accepted to Astronomy & Astrophysics (in production)

References in corpus (43)

• E pur si muove: Galiliean-invariant cosmological hydrodynamical simulations on a moving mesh
• Infrared Emission from Interstellar Dust. I. Stochastic Heating of Small Grains
• The optical morphologies of galaxies in the IllustrisTNG simulation: a comparison to Pan-STARRS observations
• The Arepo public code release
• Monte Carlo radiative transfer in protoplanetary disks
• Benchmark problems for continuum radiative transfer. High optical depths, anisotropic scattering, and polarisation
• Efficient 3D NLTE dust radiative transfer with SKIRT
• SKIRT: an Advanced Dust Radiative Transfer Code with a User-Friendly Architecture
• Three-Dimensional Dust Radiative Transfer
• Radiative transfer with POLARIS: I. Analysis of magnetic fields through synthetic dust continuum polarization measurements
• Radiative transfer in disc galaxies -- III. The observed kinematics of dusty disc galaxies
• The 2D Continuum Radiative Transfer Problem: Benchmark Results for Disk Configurations
• SWIFT: A modern highly-parallel gravity and smoothed particle hydrodynamics solver for astrophysical and cosmological applications
• Monte Carlo Radiative Transfer
• The dust-continuum size of TNG50 galaxies at $z=1-5$: a comparison with the distribution of stellar light, stars, dust and H${}_2$
• Benchmarking the Calculation of Stochastic Heating and Emissivity of Dust Grains in the Context of Radiative Transfer Simulations
• Global Calculations of Density Waves and Gap Formation in Protoplanetary Disks using a Moving Mesh
• The physical origins and dominant emission mechanisms of Lyman-alpha halos: results from the TNG50 simulation in comparison to MUSE observations
• Using hierarchical octrees in Monte Carlo radiative transfer simulations
• Hierarchical octree and k-d tree grids for 3D radiative transfer simulations
• Using 3D Voronoi grids in radiative transfer simulations
• 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
• Composite biasing in Monte Carlo 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
• SOC program for dust continuum radiative transfer
• X-ray radiative transfer in full 3D with SKIRT
• A redshift-dependent IRX-$\beta$ dust attenuation relation for TNG50 galaxies
• Polarization in Monte Carlo radiative transfer and dust scattering polarization signatures of spiral galaxies
• Tree-Structured Grid Model of Line and Polarization Variability from Massive Binaries
• 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
• Arepo-MCRT: Monte Carlo Radiation Hydrodynamics on a Moving Mesh
• Effects of Spatial Discretization in Lyman-alpha Line Radiation Transfer Simulations
• The moving mesh code Shadowfax
• Observing gas and dust in simulations of star formation with Monte Carlo radiation transport on Voronoi meshes
• A Moving Mesh Hydrodynamic Solver for ChaNGa
• ART${}^2$: A 3D Parallel Multi-wavelength Radiative Transfer Code for Continuum and Atomic and Molecular Lines
• Smart detectors for Monte Carlo radiative transfer
• Insights from Synthetic Star-forming Regions: I. Reliable Mock Observations from SPH Simulations
• 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

Cited by in corpus (2)

• TODDLERS: A New UV to Millimeter Emission Library for Star-Forming Regions. II. Star Formation Rate Indicators Using Auriga Zoom Simulations
• Towards an observationally motivated AGN dusty torus model -- II. The roles of density distribution and chemical composition of the dust