On the Integro-Differential Radiative Conductive Transfer Equation: A Modified Decomposition Method

In this work we present a solution for the radiative conductive transfer equation in spherical geometry. We discuss a semi-analytical approach to the non-linear SN problem, where the solution is constructed by Laplace transform and a decomposition method. In the present discussion, we report on arithmetic stability issues for the recursive scheme. In order to stabilise convergence by virtue of limited arithmetic precision, we propose a modification of the Adomian decomposition method by splitting terms responsible for instability by the specific choice of the source terms in the recursive scheme. This procedure results then in a finite number of modified recursion steps controlled by a splitting parameter α, whereas for all subsequent recursion steps the usual recursion scheme is used. Finally, we report on some case studies with numerical results for the solutions and convergence behaviour.