Radiance Cascades: A Novel High-Resolution Formal Solution for Multidimensional Non-LTE Radiative Transfer
arXiv:2408.14425 · doi:10.1093/rasti/rzae062
Abstract
Non-LTE radiative transfer is a key tool for modern astrophysics: it is the means by which many key synthetic observables are produced, thus connecting simulations and observations. Radiative transfer models also inform our understanding of the primary formation layers and parameters of different spectral lines, and serve as the basis of inversion tools used to infer the structure of the solar atmosphere from observations. The default approach for computing the radiation field in multidimensional solar radiative transfer models has long remained the same: a short characteristics, discrete ordinates method, formal solver. In situations with complex atmospheric structure and multiple transitions between optically-thick and -thin regimes these solvers require prohibitively high angular resolution to correctly resolve the radiation field. Here, we present the theory of radiance cascades, a technique designed to exploit structure inherent to the radiation field, allowing for efficient reuse of calculated samples, thus providing a very high-resolution result at a fraction of the computational cost of existing methods. We additionally describe our implementation of this method in the DexRT code, and present initial results of the synthesis of a snapshot of a magnetohydrodynamic model of a solar prominence formed via levitation-condensation. The approach presented here provides a credible route for routinely performing multidimensional radiative transfer calculations free from so-called ray effects, and scaling high-quality non-LTE models to next-generation high-performance computing systems with GPU accelerators.
21 pages, 14 figures, ancillary videos. To be submitted to RASTI
References in corpus (25)
- HEALPix -- a Framework for High Resolution Discretization, and Fast Analysis of Data Distributed on the Sphere
- The stellar atmosphere simulation code Bifrost
- Physics of Solar Prominences: I - Spectral Diagnostics and Non-LTE Modelling
- MPI-AMRVAC 2.0 for Solar and Astrophysical Applications
- MPI-AMRVAC 3.0: updates to an open-source simulation framework
- DELO-Bezier formal solutions of the polarized radiative transfer equation
- PORTA: A three-dimensional multilevel radiative transfer code for modeling the intensity and polarization of spectral lines with massively parallel computers
- Simulations of solar filament fine structures and their counterstreaming flows
- Partial redistribution in 3D non-LTE radiative transfer in solar atmosphere models
- Three-dimensional modeling of chromospheric spectral lines in a simulated active region
- Prominence formation by levitation-condensation at extreme resolutions
- Thermal instabilities: Fragmentation and field misalignment of filament fine structure
- The Failure of Monte Carlo Radiative Transfer at Medium to High Optical Depths
- The Lightweaver Framework for NLTE Radiative Transfer in Python
- Magritte, a modern software library for 3D radiative transfer: I. Non-LTE atomic and molecular line modelling
- 1.5D NLTE spectral synthesis of a 3D filament/prominence simulation
- Magritte, a modern software library for 3D radiative transfer: II. Adaptive ray-tracing, mesh construction and reduction
- Flare Kernels May be Smaller than You Think: Modelling the Radiative Response of Chromospheric Plasma Adjacent to a Solar Flare
- Numerical non-LTE 3D radiative transfer using a multigrid method
- Near optimal angular quadratures for polarised radiative transfer
- COCOPLOT: COlor COllapsed PLOTting software : Using color to view 3D data as a 2D image
- Full non-LTE spectral line formation I. Setting the stage
- On the Importance of Ca II Photoionisation by the Hydrogen Lyman Transitions in Solar Flare Models
- Transverse kink oscillations of inhomogeneous prominence threads: numerical analysis and H forward modelling
- Long characteristics vs. short characteristics in 3D radiative transfer simulations of polarized radiation