A New Block Preconditioner for Implicit Runge-Kutta Methods for Parabolic PDE

10/22/2020
by   Md Masud Rana, et al.
0

A new preconditioner based on a block LDU factorization with algebraic multigrid subsolves for scalability is introduced for the large, structured systems appearing in implicit Runge-Kutta time integration of parabolic partial differential equations. This preconditioner is compared in condition number and eigenvalue distribution, and in numerical experiments with others in the literature: block Jacobi, block Gauss-Seidel, and the optimized block Gauss-Seidel method of 10.4173/mic.2006.2.3. Experiments are run with implicit Runge-Kutta stages up to s=7, and it is found that the new preconditioner outperforms the others, with the improvement becoming more pronounced as spatial discretization is refined and as temporal order is increased.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
06/27/2020

Tusas: A fully implicit parallel approach for coupled nonlinear equations

We develop a fully-coupled, fully-implicit approach for phase-field mode...
research
11/02/2020

Constructing Polynomial Block Methods

The recently introduced polynomial time integration framework proposes a...
research
09/03/2019

Numerical valuation of Bermudan basket options via partial differential equations

We study the principal component analysis (PCA) based approach introduce...
research
02/04/2020

Compress-and-restart block Krylov subspace methods for Sylvester matrix equations

Block Krylov subspace methods (KSMs) comprise building blocks in many st...
research
12/23/2020

Optimal and Low-Memory Near-Optimal Preconditioning of Fully Implicit Runge-Kutta Schemes for Parabolic PDEs

Runge-Kutta (RK) schemes, especially Gauss-Legendre and some other fully...
research
05/19/2023

Implicit low-rank Riemannian schemes for the time integration of stiff partial differential equations

We propose two implicit numerical schemes for the low-rank time integrat...

Please sign up or login with your details

Forgot password? Click here to reset