Abstract:To understand the structures of complex astrophysical objects, 3D numerical simulations of radiative transfer processes are invaluable. For Monte Carlo radiative transfer, the most common radiative transfer method in 3D, the design of a spatial grid is important and non-trivial. Common choices include hierarchical octree and unstructured Voronoi grids, each of which has advantages and limitations. Tetrahedral grids, commonly used in ray-tracing computer graphics, can be an interesting alternative option. We aim to investigate the possibilities, advantages, and limitations of tetrahedral grids in the context of Monte Carlo radiative transfer. In particular, we want to compare the performance of tetrahedral grids to other commonly used grid structures. We implemented a tetrahedral grid structure, based on the open-source library TetGen, in the generic Monte Carlo radiative transfer code SKIRT. Tetrahedral grids can be imported from external applications or they can be constructed and adaptively refined within SKIRT. We implemented an efficient grid traversal method based on Pl\"ucker coordinates and Pl\"ucker products. The correct implementation of the tetrahedral grid construction and the grid traversal algorithm in SKIRT were validated using 2D radiative transfer benchmark problems. Using a simple 3D model, we compared the performance of tetrahedral, octree, and Voronoi grids. With a constant cell count, the octree grid outperforms the tetrahedral and Voronoi grids in terms of traversal speed, whereas the tetrahedral grid is poorer than the other grids in terms of grid quality. All told, we find that the performance of tetrahedral grids is relatively poor compared to octree and Voronoi grids. Although the adaptively constructed tetrahedral grids might not be favourable in most media representative of astrophysical simulation models, they still form an interesting unstructured alternative to Voronoi grids for specific applications. In particular, they might prove useful for radiative transfer post-processing of hydrodynamical simulations run on tetrahedral or unstructured grids.

Tetrahedral grids in Monte Carlo radiative transfer

Tetrahedral grids in Monte Carlo radiative transfer

A Three-Dimensional Hybrid Grid: Dragon Grid

A Novel Approach of Three-Dimensional Hybrid Grid Methodology: Part 1. Grid Generation

Development of Three-Dimensional DRAGON Grid Technology

Irregular grids for 3D NLTE radiative transfer in stellar atmospheres

HyGrid: A CPU-GPU Hybrid Convolution-Based Gridding Algorithm in Radio Astronomy.

Comparison of Ray-Tracing Methods in Self-Adaptive Grid for Inhomogeneous Medium

Topology Abstraction of Surface Models for Three-Dimensional Grid Generation

Quasi-Monte Carlo Ray Tracing Algorithm for Radiative Flux Distribution Simulation

Transformation from Octree Grid to Unstructured Hybrid Grid

A combined analytical and Monte Carlo method for detailed simulations of antiscatter grids in x-ray medical imaging: implementing scatter within the grid

Efficient 3D NLTE dust radiative transfer with SKIRT

Hybrid Radiation Hydrodynamics scheme with gravity tree-based adaptive optimization algorithm

Interfacing unstructured tetrahedron grids to structured-grid FDTD

Efficient Kd-Tree Construction for Ray Tracing Using Ray Distribution Sampling

A Monte Carlo Method for 3D Radiative Transfer Equations with Multifractional Singular Kernels

Observing gas and dust in simulations of star formation with Monte Carlo radiation transport on Voronoi meshes

Efficient radiative transfer modelling with SKIRT

A New Discrete Implicit Monte Carlo Scheme for Simulating Radiative Transfer Problems

Development of a Clear‐Sky 3D Sub‐Grid Terrain Solar Radiative Effect Parameterization Scheme Based on the Mountain Radiation Theory