Krylov Iterative Methods and the Degraded Effectiveness of Diffusion Synthetic Acceleration for Multidimensional S{sub N} Calculations in Problems with Material Discontinuities
Journal Article
·
· Nuclear Science and Engineering
- Los Alamos National Laboratory (United States)
A loss in the effectiveness of diffusion synthetic acceleration (DSA) schemes has been observed with certain S{sub N} discretizations on two-dimensional Cartesian grids in the presence of material discontinuities. We will present more evidence supporting the conjecture that DSA effectiveness will degrade for multidimensional problems with discontinuous total cross sections, regardless of the particular physical configuration or spatial discretization. Fourier analysis and numerical experiments help us identify a set of representative problems for which established DSA schemes are ineffective, focusing on diffusive problems for which DSA is most needed. We consider a lumped, linear discontinuous spatial discretization of the S{sub N} transport equation on three-dimensional, unstructured tetrahedral meshes and look at a fully consistent and a 'partially consistent' DSA method for this discretization. The effectiveness of both methods is shown to degrade significantly. A Fourier analysis of the fully consistent DSA scheme in the limit of decreasing cell optical thickness supports the view that the DSA itself is failing when material discontinuities are present in a problem. We show that a Krylov iterative method, preconditioned with DSA, is an effective remedy that can be used to efficiently compute solutions for this class of problems. We show that as a preconditioner to the Krylov method, a partially consistent DSA method is more than adequate. In fact, it is preferable to a fully consistent method because the partially consistent method is based on a continuous finite element discretization of the diffusion equation that can be solved relatively easily. The Krylov method can be implemented in terms of the original S{sub N} source iteration coding with only slight modification. Results from numerical experiments show that replacing source iteration with a preconditioned Krylov method can efficiently solve problems that are virtually intractable with accelerated source iteration.
- OSTI ID:
- 20804942
- Journal Information:
- Nuclear Science and Engineering, Journal Name: Nuclear Science and Engineering Journal Issue: 3 Vol. 147; ISSN NSENAO; ISSN 0029-5639
- Country of Publication:
- United States
- Language:
- English
Similar Records
Krylov iterative methods applied to multidimensional S[sub n] calculations in the presence of material discontinuities
Fully Consistent Diffusion Synthetic Acceleration of Linear Discontinuous S{sub N} Transport Discretizations on Unstructured Tetrahedral Meshes
On the degraded effectiveness of diffusion synthetic acceleration for multidimensional sn calculations in the presence of material discontinuities
Conference
·
Mon Dec 31 23:00:00 EST 2001
·
OSTI ID:976408
Fully Consistent Diffusion Synthetic Acceleration of Linear Discontinuous S{sub N} Transport Discretizations on Unstructured Tetrahedral Meshes
Journal Article
·
Mon Jul 15 00:00:00 EDT 2002
· Nuclear Science and Engineering
·
OSTI ID:20804791
On the degraded effectiveness of diffusion synthetic acceleration for multidimensional sn calculations in the presence of material discontinuities
Conference
·
Mon Dec 31 23:00:00 EST 2001
·
OSTI ID:976407