Journal
SIAM JOURNAL ON SCIENTIFIC COMPUTING
Volume 44, Issue 2, Pages A636-A663Publisher
SIAM PUBLICATIONS
DOI: 10.1137/21M1390438
Keywords
Runge-Kutta; numerical PDEs; preconditioning; Navier--Stokes; time integration
Categories
Funding
- U.S. Department of Energy by Lawrence Livermore National Laboratory [DE-AC52-07NA27344 (LLNL-JRNL817953), LA-UR-20-30412]
- Lawrence Livermore National Laboratory [B639443]
- Australian Government Research Training Program (RTP) scholarship
Ask authors/readers for more resources
This paper presents a theoretical and algorithmic framework for solving the nonlinear equations arising from IRK methods, offering new linearizations and demonstrating their effectiveness and stability. By applying these methods to challenging fluid flow problems, the high accuracy and efficiency of high-order IRK methods are showcased.
Fully implicit Runge--Kutta (IRK) methods have many desirable accuracy and stability properties as time integration schemes, but high-order IRK methods are not commonly used in practice with large-scale numerical PDEs because of the difficulty of solving the stage equations. This paper introduces a theoretical and algorithmic framework for solving the nonlinear equations that arise from IRK methods (and discontinuous Galerkin discretizations in time) applied to nonlinear numerical PDEs, including PDEs with algebraic constraints. Several new linearizations of the nonlinear IRK equations are developed, offering faster and more robust convergence than the often-considered simplified Newton, as well as an effective preconditioner for the true Jacobian if exact Newton iterations are desired. Inverting these linearizations requires solving a set of block 2 \times 2 systems. Under quite general assumptions, it is proven that the preconditioned 2 \times 2 operator's condition number is bounded by a small constant close to one, independent of the spatial discretization, spatial mesh, and time step, and with only weak dependence on the number of stages or integration accuracy. Moreover, the new method is built using the same preconditioners needed for backward Euler-type time stepping schemes, so can be readily added to existing codes. The new methods are applied to several challenging fluid flow problems, including the compressible Euler and Navier-Stokes equations, and the vorticity-streamfunction formulation of the incompressible Euler and Navier-Stokes equations. Up to 10th-order accuracy is demonstrated using Gauss IRK, while in all cases fourth-order Gauss IRK requires roughly half the number of preconditioner applications as required by standard Singly diagonally implicit Runge--Kutta methods.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available