During the last years, various algorithms solving the transport equation in slab geometry have been proposed in the literature [2–4,6,7]. Most of them use a convergence acceleration based on a Picard type algorithm called the source iteration method [4,7]. New non-accelerated algorithms have been recently introduced [1,2]. We develop here a non-accelerated method based on two points: a splitting of the collision operator and an infinite dimensional adaptation of the minimal residual method. We first prove the theoretical convergence of the method in the frame of non-reflecting boundary conditions. Then we compare numerically this method with existing non accelerated schemes. It gives good results which could even be further improved by adding a DSA kind acceleration.
[1]
Edward W. Larsen,et al.
Synthetic acceleration methods for linear transport problems with highly anisotropic scattering
,
1991
.
[2]
P. Lascaux,et al.
Analyse numérique matricielle appliquée a l'art de l'ingénieur
,
1987
.
[3]
Edward W. Larsen,et al.
UNCONDITIONALLY STABLE DIFFUSION-SYNTHETIC ACCELERATION METHODS FOR THE SLAB GEOMETRY DISCRETE ORDINATES EQUATIONS. PART I: THEORY.
,
1982
.
[4]
M. Nicolet,et al.
Nouveaux algorithmes performants en théorie du transport
,
1998
.
[5]
Ivo Marek,et al.
Frobenius Theory of Positive Operators: Comparison Theorems and Applications
,
1970
.