English

TREVR: A general $N log^2 N$ radiative transfer algorithm

Instrumentation and Methods for Astrophysics 2019-03-06 v1 Cosmology and Nongalactic Astrophysics Astrophysics of Galaxies

Abstract

We present TREVR (Tree-based REVerse Ray Tracing), a general algorithm for computing the radiation field, including absorption, in astrophysical simulations. TREVR is designed to handle large numbers of sources and absorbers; it is based on a tree data structure and is thus suited to codes that use trees for their gravity or hydrodynamics solvers (e.g. Adaptive Mesh Refinement). It achieves computational speed while maintaining a specified accuracy via controlled lowering of the resolution of both sources and rays from each source. TREVR computes the radiation field in order NlogNsource{N\log N_{source}} time without absorption and order NlogNsourcelogN{N \log N_{source} \log{N}} time with absorption. These scalings arise from merging sources of radiation according to an opening angle criterion and walking the tree structure to trace a ray to a depth that gives the chosen accuracy for absorption. The absorption-depth refinement criterion is unique to TREVR. We provide a suite of tests demonstrating the algorithm's ability to accurately compute fluxes, ionization fronts and shadows.

Cite

@article{arxiv.1902.08083,
  title  = {TREVR: A general $N log^2 N$ radiative transfer algorithm},
  author = {J. J. Grond and R. M. Woods and J. W. Wadsley and H. M. P. Couchman},
  journal= {arXiv preprint arXiv:1902.08083},
  year   = {2019}
}

Comments

15 pages, 15 figures and 1 table. Accepted by MNRAS

R2 v1 2026-06-23T07:47:14.420Z