Efficient Order-Optimal Preconditioners for Implicit Runge-Kutta and Runge-Kutta-Nyström Methods Applicable to a Large Class of Parabolic and Hyperbolic PDEs

06/17/2022
by   Michael R. Clines, et al.
0

We generalize previous work by Mardal, Nilssen, and Staff (2007, SIAM J. Sci. Comp. v. 29, pp. 361-375) and Rana, Howle, Long, Meek, and Milestone (2021, SIAM J. Sci. Comp. v. 43, p. 475-495) on order-optimal preconditioners for parabolic PDEs to a larger class of differential equations and methods. The problems considered are those of the forms u_t=-𝒦u+g and u_tt=-𝒦u+g, where the operator 𝒦 is defined by 𝒦u:=-∇·(α∇ u)+β u and the functions α and β are restricted so that α>0, and β≥0. The methods considered are A-stable implicit Runge–Kutta methods for the parabolic equation and implicit Runge–Kutta–Nyström methods for the hyperbolic equation. We prove the order optimality of a class of block preconditioners for the stage equation system arising from these problems, and furthermore we show that the LD and DU preconditioners of Rana et al. are in this class. We carry out numerical experiments on several test problems in this class – the 2D diffusion equation, Pennes bioheat equation, the wave equation, and the Klein–Gordon equation, with both constant and variable coefficients. Our experiments show that these preconditioners, particularly the LD preconditioner, are successful at reducing the condition number of the systems as well as improving the convergence rate and solve time for GMRES applied to the stage equations.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
01/19/2023

Convergence Analysis of a Krylov Subspace Spectral Method for the 1-D Wave Equation in an Inhomogeneous Medium

Krylov subspace spectral (KSS) methods are high-order accurate, explicit...
research
06/24/2020

Semi-implicit Taylor schemes for stiff rough differential equations

We study a class of semi-implicit Taylor-type numerical methods that are...
research
02/22/2022

Acceleration of Gossip Algorithms through the Euler-Poisson-Darboux Equation

Gossip algorithms and their accelerated versions have been studied exclu...
research
02/05/2021

Convergence analysis of explicit stabilized integrators for parabolic semilinear stochastic PDEs

Explicit stabilized integrators are an efficient alternative to implicit...
research
01/21/2022

Sparse grid implementation of a fixed-point fast sweeping WENO scheme for Eikonal equations

Fixed-point fast sweeping methods are a class of explicit iterative meth...
research
01/10/2023

Exponential Runge-Kutta Parareal for Non-Diffusive Equations

Parareal is a well-known parallel-in-time algorithm that combines a coar...
research
10/03/2022

Nonlinear Reconstruction for Operator Learning of PDEs with Discontinuities

A large class of hyperbolic and advection-dominated PDEs can have soluti...

Please sign up or login with your details

Forgot password? Click here to reset