Radiative–conductive transfer equation in spherical geometry: arithmetic stability for decomposition using the condition number criterion

The radiative–conductive transfer equation in the S N approximation for spherical geometry is solved using a modified decomposition method. The focus of this work is to show how to distribute the source terms in the recursive equation system in order to guarantee arithmetic stability and thus numerical convergence of the obtained solution, guided by a condition number criterion. Some examples are compared with results from literature and parameter combinations are analyzed, for which the condition number analysis indicates convergence of the solutions obtained by the recursive scheme.