Spectral Homogenization of the Radiative Transfer Equation via Low-Rank Tensor Train Decomposition

TL;DR

This paper demonstrates that the spectral radiative transfer equation can be efficiently approximated using low-rank tensor train decompositions, significantly reducing computational cost while maintaining high spectral fidelity across various media.

Contribution

The authors introduce a tensor train-based homogenization framework that reveals the spectral complexity of radiative transfer equations has a bounded low-rank structure, enabling efficient computation.

Findings

Tensor train ranks saturate at low values regardless of spectral resolution and physical parameters.

Quantized tensor train representations achieve sub-linear storage scaling.

The homogenized approach outperforms correlated-k distribution in accuracy at equal computational cost.

Abstract

Radiative transfer in absorbing-scattering media requires solving a transport equation across a spectral domain with 10^5 - 10^6 molecular absorption lines. Line-by-line (LBL) computation is prohibitively expensive, while existing approximations sacrifice spectral fidelity. We show that the Young-measure homogenization framework produces solution tensors I that admit low-rank tensor-train (TT) decompositions whose bond dimensions remain bounded as the spectral resolution Ns increases. Using molecular line parameters from the HITRAN database for H2O and CO2, we demonstrate that: (i) the TT rank saturates at r = 8 (at tolerance e = 10^-6) from Ns = 16 to 4096, independent of single-scattering albedo, Henyey-Greenstein asymmetry, temperature, and pressure; (ii) quantized tensor-train (QTT) representations achieve sub-linear storage scaling; (iii) in a controlled comparison using identical opacity data and transport solver, the homogenized approach achieves over an order of magnitude lower L2 error than the correlated-k distribution at equal cost; and (iv) for atomic plasma opacity (aluminum at 60 eV, TOPS database), the TT rank saturates at r = 15 with fundamentally different spectral structure (bound-bound and bound-free transitions spanning 12 decades of dynamic range), confirming that rank boundedness is a property of the transport equation rather than any particular opacity source. These results establish that the spectral complexity of radiative transfer has a finite effective rank exploitable by tensor decomposition, complementing the spatial-angular compression achieved by existing TT and dynamical low-rank approaches.