Uniform convergence of diffusion synthetic acceleration for heterogeneous slab transport
Abstract
The diffusion synthetic accelerated (DSA) source iteration is a standard solver for the radiative transfer equation. Using Fourier analysis, a convergence rate ρ∞(c)≤0.2247 c with c the maximum ratio of scattering to total cross section has been established for an infinite homogeneous medium. For slab geometry with inflow boundary conditions and arbitrary bounded cross sections we prove that the spectral radius of the DSA iteration is at most ρ∞(c). We show that this convergence rate carries over to a variational discretization of the DSA iteration on every conforming tensor-product Galerkin space whose angular factor contains the constants. The analysis rests on an exact min--max characterization of the spectral radius. From it we derive a checkable sufficient condition for a given rate. We verify this condition using a suitable energy-stable projection, a weighted angular average whose weight is chosen so that the condition holds with the rate ρ∞(c). Since the projection maps discrete spaces into discrete spaces, the argument applies verbatim to the discrete iteration. As a consequence, the condition number of the preconditioned system is uniformly bounded by 1+0.29 c, which implies rapid convergence of the preconditioned conjugate gradients method.
Related
- Transport-Matched Penalties for Diffusion Synthetic Acceleration of Polytopic Discontinuous Galerkin Discretisations
- Convergence and acceleration of a nonlinear fixed-point iteration for computing the Fitness Centrality of general graphs
- Asymptotic behavior of Eckhoff's method for convergence acceleration of Dirac eigenfunction expansions
- Convergence rates for generative drifting flows: fixed-scale obstructions and multihead acceleration
- Polynomial Acceleration of Alternating Projections: Chebyshev Acceleration, Convergence Rates, and Retractions
BibTeX
@misc{schlottbom2026uniform,
title = {{Uniform convergence of diffusion synthetic acceleration for heterogeneous slab transport}},
author = {Matthias Schlottbom},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2609.31346},
eprint = {2609.31346},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2609.31346},
}
RIS
TY - UNPB
TI - Uniform convergence of diffusion synthetic acceleration for heterogeneous slab transport
AU - Matthias Schlottbom
PY - 2026
DA - 2026/09/25
PB - arXiv
DO - 10.48550/arxiv.2609.31346
UR - https://arxiv.org/abs/2609.31346
AB - The diffusion synthetic accelerated (DSA) source iteration is a standard solver for the radiative transfer equation. Using Fourier analysis, a convergence rate ρ_∞(c)≤0.2247 c with c the maximum ratio of scattering to total cross section has been established for an infinite homogeneous medium. For slab geometry with inflow boundary conditions and arbitrary bounded cross sections we prove that the spectral radius of the DSA iteration is at most ρ_∞(c). We show that this convergence rate carries over to a variational discretization of the DSA iteration on every conforming tensor-product Galerkin space whose angular factor contains the constants. The analysis rests on an exact min--max characterization of the spectral radius. From it we derive a checkable sufficient condition for a given rate. We verify this condition using a suitable energy-stable projection, a weighted angular average whose weight is chosen so that the condition holds with the rate ρ_∞(c). Since the projection maps discrete spaces into discrete spaces, the argument applies verbatim to the discrete iteration. As a consequence, the condition number of the preconditioned system is uniformly bounded by 1+0.29 c, which implies rapid convergence of the preconditioned conjugate gradients method.
ER -