Design and implementation of DIRK integrators for stiff systems

  • Authors:
  • Roger Alexander

  • Affiliations:
  • Department of Mathematics, 400 Carver Hall, Iowa State University, Ames, IA

  • Venue:
  • Applied Numerical Mathematics
  • Year:
  • 2003

Quantified Score

Hi-index 0.01

Visualization

Abstract

This paper presents new four-stage diagonally implicit Runge-Kutta integration formulas for stiff initial value problems, designed to be L-stable and have optimal order of accuracy. The design makes estimation of local error and interpolation of the solution possible without additional stages. It also enables improved prediction of stage values in a software implementation. Correctly monitoring the convergence rate of simplified Newton iterations improves the implementation too, by providing a sound basis for deciding on Jacobian updates.