Conditioning of implicit Runge-Kutta integration for finite element approximation of linear diffusion equations on anisotropic meshes

Weizhang Huang, Lennard Kamenski, Jens Lang

Published 2017-03-19Version 1

The conditioning of implicit Runge-Kutta (RK) integration for linear finite element approximation of diffusion equations on general anisotropic meshes is investigated. Bounds are established for the condition number of the resulting linear system with and without diagonal preconditioning for the implicit Euler (the simplest implicit RK method) and general implicit RK methods. It is shown that the conditioning of an implicit RK method behaves like that of the implicit Euler method. Two solution strategies are considered for the linear system resulting from general implicit RK integration: the simultaneous solution where the system is solved as a whole and a successive solution which follows the commonly used implementation of implicit RK methods to first transform the system into a number of smaller systems using the Jordan normal form of the RK matrix and then solve them successively. The obtained bounds for the condition number have explicit geometric interpretations and take the interplay between the diffusion matrix and the mesh geometry into full consideration. They show that there are three mesh-dependent factors that can affect the conditioning: the number of elements, the mesh nonuniformity measured in the Euclidean metric, and the mesh nonuniformity with respect to the inverse of the diffusion matrix. They also reveal that the preconditioning using the diagonal of the system matrix, the mass matrix, or the lumped mass matrix can effectively eliminate the effects of the mesh nonuniformity measured in the Euclidean metric. Illustrative numerical examples are given.