Application of the Jacobian-Free Newton-Krylov Method to Nonlinear Acceleration of Transport Source Iteration in Slab Geometry
The use of the Jacobian-free Newton-Krylov (JFNK) method within the context of nonlinear diffusion acceleration (NDA) of source iteration is explored. The JFNK method is a synergistic combination of Newton's method as the nonlinear solver and Krylov methods as the linear solver. JFNK methods do not form or store the Jacobian matrix, and Newton's method is executed via probing the nonlinear discrete function to approximate the required matrix-vector products. Current application of NDA relies upon a fixed-point, or Picard, iteration to resolve the nonlinearity. We show that the JFNK method can be used to replace this Picard iteration with a Newton iteration. The Picard linearization is retained as a preconditioner. We show that the resulting JFNK-NDA capability provides benefit in some regimes. Furthermore, we study the effects of a two-grid approach, and the required intergrid transfers when the higher-order transport method is solved on a fine mesh compared to the low-order acceleration problem.
- Research Organization:
- Idaho National Laboratory (INL)
- Sponsoring Organization:
- USDOE
- DOE Contract Number:
- AC07-05ID14517
- OSTI ID:
- 1016204
- Report Number(s):
- INL/JOU-11-22416
- Journal Information:
- Nuclear Science and Engineering, Journal Name: Nuclear Science and Engineering Journal Issue: 2 Vol. 167; ISSN NSENAO; ISSN 0029-5639
- Country of Publication:
- United States
- Language:
- English
Similar Records
Implementation of the Jacobian-free Newton–Krylov method for solving the first-order ice sheet momentum balance
A Parallel, Fully Coupled, Fully Implicit Solution to Reactive Transport in Porous Media Using the Preconditioned Jacobian-Free Newton-Krylov Method