Acceleration of implicit schemes for large systems of nonlinear differential-algebraic equations

Abstract

When solving large systems of nonlinear differential-algebraic equations by implicit schemes, each integration step requires the solution of a system of large nonlinear algebraic equations. The latter is solved by an inexact Newton method which, in its turn, leads to a set of large linear systems commonly solved by a Krylov subspace iterative method. The efficiency of the whole process depends on the initial guesses for the inexact Newton and the Krylov subspace methods. An inexpensive approach is proposed and justified that computes good initial guesses for these methods. It requires a subspace of small dimension and the use of line search and trust region for the inexact Newton method and Petrov-Galerkin for the Krylov subspace method. Numerical examples are included to illustrate the effectiveness of the proposed approach.