Krylov iterative methods and synthetic acceleration for transport in binary statistical media
- Los Alamos National Laboratory
In particle transport applications there are numerous physical constructs in which heterogeneities are randomly distributed. The quantity of interest in these problems is the ensemble average of the flux, or the average of the flux over all possible material 'realizations.' The Levermore-Pomraning closure assumes Markovian mixing statistics and allows a closed, coupled system of equations to be written for the ensemble averages of the flux in each material. Generally, binary statistical mixtures are considered in which there are two (homogeneous) materials and corresponding coupled equations. The solution process is iterative, but convergence may be slow as either or both materials approach the diffusion and/or atomic mix limits. A three-part acceleration scheme is devised to expedite convergence, particularly in the atomic mix-diffusion limit where computation is extremely slow. The iteration is first divided into a series of 'inner' material and source iterations to attenuate the diffusion and atomic mix error modes separately. Secondly, atomic mix synthetic acceleration is applied to the inner material iteration and S{sup 2} synthetic acceleration to the inner source iterations to offset the cost of doing several inner iterations per outer iteration. Finally, a Krylov iterative solver is wrapped around each iteration, inner and outer, to further expedite convergence. A spectral analysis is conducted and iteration counts and computing cost for the new two-step scheme are compared against those for a simple one-step iteration, to which a Krylov iterative method can also be applied.
- Research Organization:
- Los Alamos National Laboratory (LANL), Los Alamos, NM (United States)
- Sponsoring Organization:
- USDOE
- DOE Contract Number:
- AC52-06NA25396
- OSTI ID:
- 960633
- Report Number(s):
- LA-UR-08-05581; LA-UR-08-5581; JCTPAH; TRN: US201006%%1267
- Journal Information:
- Journal of Computational Physics, Vol. 228, Issue 22; ISSN 0021-9991
- Country of Publication:
- United States
- Language:
- English
Similar Records
Modified Closures in Monte Carlo Algorithms for Diffusive Binary Stochastic Media Transport Problems
Benchmark Solutions for Radiation Transport in Stochastic Media with Inhomogeneous Material Statistics