A Parallel, Fully Coupled, Fully Implicit Solution to Reactive Transport in Porous Media Using the Preconditioned Jacobian-Free Newton-Krylov Method
Journal Article
·
· Advances in Water Resources
OSTI ID:1062744
Modeling large multicomponent reactive transport systems in porous media is particularly challenging when the governing partial differential algebraic equations (PDAEs) are highly nonlinear and tightly coupled due to complex nonlinear reactions and strong solution-media interactions. Here we present a preconditioned Jacobian-Free Newton-Krylov (JFNK) solution approach to solve the governing PDAEs in a fully coupled and fully implicit manner. A well-known advantage of the JFNK method is that it does not require explicitly computing and storing the Jacobian matrix during Newton nonlinear iterations. Our approach further enhances the JFNK method by utilizing physics-based, block preconditioning and a multigrid algorithm for efficient inversion of the preconditioner. This preconditioning strategy accounts for self- and optionally, cross-coupling between primary variables using diagonal and off-diagonal blocks of an approximate Jacobian, respectively. Numerical results are presented demonstrating the efficiency and massive scalability of the solution strategy for reactive transport problems involving strong solution-mineral interactions and fast kinetics. We found that the physics-based, block preconditioner significantly decreases the number of linear iterations, directly reducing computational cost; and the strongly scalable algebraic multigrid algorithm for approximate inversion of the preconditioner leads to excellent parallel scaling performance.
- Research Organization:
- Idaho National Laboratory (INL)
- Sponsoring Organization:
- DOE - SC
- DOE Contract Number:
- AC07-05ID14517
- OSTI ID:
- 1062744
- Report Number(s):
- INL/JOU-11-22645
- Journal Information:
- Advances in Water Resources, Journal Name: Advances in Water Resources Vol. 53
- Country of Publication:
- United States
- Language:
- English
Similar Records
REACTIVE TRANSPORT MODELING USING A PARALLEL FULLY-COUPLED SIMULATOR BASED ON PRECONDITIONED JACOBIAN-FREE NEWTON-KRYLOV
Acceleration of k-Eigenvalue / Criticality Calculations using the Jacobian-Free Newton-Krylov Method
Physics-based preconditioning and the Newton-Krylov method for non-equilibrium radiation diffusion
Conference
·
Fri Jun 01 00:00:00 EDT 2012
·
OSTI ID:1047191
Acceleration of k-Eigenvalue / Criticality Calculations using the Jacobian-Free Newton-Krylov Method
Journal Article
·
Mon Jan 31 23:00:00 EST 2011
· Nuclear Science and Engineering
·
OSTI ID:1010691
Physics-based preconditioning and the Newton-Krylov method for non-equilibrium radiation diffusion
Journal Article
·
Sat May 20 00:00:00 EDT 2000
· Journal of Computational Physics
·
OSTI ID:20067701