A 3D radiative transfer framework: I. non-local operator splitting and continuum scattering problems

TL;DR

This paper introduces a flexible 3D radiative transfer framework using operator splitting and iterative solvers, enabling accurate modeling of scattering-dominated environments with potential applications in static media.

Contribution

The paper presents a novel 3D radiative transfer method combining long characteristics, operator splitting, and iterative solvers, with validation against 1D solutions.

Findings

Efficient solution of linear systems with Gauss-Seidel and Jordan methods

Accurate 3D results comparable to 1D benchmarks

Framework applicable to static media with low velocity fields

Abstract

We describe a highly flexible framework to solve 3D radiation transfer problems in scattering dominated environments based on a long characteristics piece-wise parabolic formal solution and an operator splitting method. We find that the linear systems are efficiently solved with iterative solvers such as Gauss-Seidel and Jordan techniques. We use a sphere-in-a-box test model to compare the 3D results to 1D solutions in order to assess the accuracy of the method. We have implemented the method for static media, however, it can be used to solve problems in the Eulerian-frame for media with low velocity fields.