Abstract
Nonlinear preconditioning is a globalization technique for Newton’s method applied to systems of equations with unbalanced nonlinearities, in which nonlinear residual norm reduction stagnates due to slowly evolving subsets of the degrees of freedom. Even though the Newton corrections may effectively be sparse, a standard Newton method still requires large ill-conditioned linear systems resulting from global linearizations of the nonlinear residual to be solved at each step. Nonlinear preconditioners may enable faster global convergence by shifting work to where it is most strategic, on subsets of the original system. They require additional computation per outer iteration while aiming for many fewer outer iterations and correspondingly fewer global synchronizations. In this work, we improve upon previous nonlinear preconditioning implementations by introducing parameters that allow turning off nonlinear preconditioning during outer Newton iterations where it is not needed. Numerical experiments show that the adaptive nonlinear preconditioning algorithm has performance similar to monolithically applied nonlinear preconditioning, preserving robustness for some challenging problems representative of several PDE-based applications while saving work on nonlinear subproblems.
Original language | English (US) |
---|---|
Pages (from-to) | A1171-A1186 |
Journal | SIAM Journal on Scientific Computing |
Volume | 40 |
Issue number | 2 |
DOIs | |
State | Published - 2018 |
Bibliographical note
Publisher Copyright:© 2018 Society for Industrial and Applied Mathematics.
Keywords
- Domain decomposition
- Field splitting
- Newton’s method
- Nonlinear equations
- Nonlinear preconditioning
ASJC Scopus subject areas
- Computational Mathematics
- Applied Mathematics