Globalization technique for projected Newton–Krylov methods

Jinhai Chen, Kees Vuik

Research output: Contribution to journalArticleScientificpeer-review

31 Downloads (Pure)

Abstract

Large-scale systems of nonlinear equations appear in many applications. In various applications, the solution of the nonlinear equations should also be in a certain interval. A typical application is a discretized system of reaction diffusion equations. It is well known that chemical species should be positive otherwise the solution is not physical and in general blow up occurs. Recently, a projected Newton method has been developed, which can be used to solve this type of problems. A drawback is that the projected Newton method is not globally convergent. This motivates us to develop a new feasible projected Newton–Krylov algorithm for solving a constrained system of nonlinear equations. Combined with a projected gradient direction, our feasible projected Newton–Krylov algorithm circumvents the non-descent drawback of search directions which appear in the classical projected Newton methods. Global and local superlinear convergence of our approach is established under some standard assumptions. Numerical experiments are used to illustrate that the new projected Newton method is globally convergent and is a significate complementarity for Newton–Krylov algorithms known in the literature.
Original languageEnglish
Pages (from-to)661-674
Number of pages14
JournalInternational Journal for Numerical Methods in Engineering
Volume110
Issue number7
DOIs
Publication statusPublished - 11 Oct 2017

Keywords

  • nonlinear equations with constraints
  • projected Newton–Krylov method
  • projected gradient direction
  • global convergence

Fingerprint Dive into the research topics of 'Globalization technique for projected Newton–Krylov methods'. Together they form a unique fingerprint.

  • Cite this