Fast and high-order approximation of parabolic equations using hierarchical direct solvers and implicit Runge-Kutta methods

by   Ke Chen, et al.

A stable and high-order accurate solver for linear and nonlinear parabolic equations is presented. An additive Runge-Kutta method is used for the time stepping, which integrates the linear stiff terms by an implicit singly diagonally implicit Runge-Kutta (ESDIRK) method and the nonlinear terms by an explicit Runge-Kutta (ERK) method. In each time step, the implicit solve is performed by the recently developed Hierarchical Poincaré-Steklov (HPS) method. This is a fast direct solver for elliptic equations that decomposes the space domain into a hierarchical tree of subdomains and builds spectral collocation solvers locally on the subdomains. These ideas are naturally combined in the presented method since the singly diagonal coefficient in ESDIRK and a fixed time-step ensures that the coefficient matrix in the implicit solve of HPS remains the same for all time stages. This means that the precomputed inverse can be efficiently reused, leading to a scheme with complexity (in two dimensions) 𝒪(N^1.5) for the precomputation where the solution operator to the elliptic problems is built, and then 𝒪(N) for each time step. The stability of the method is proved for first order in time and any order in space, and numerical evidence substantiates a claim of stability for a much broader class of time discretization methods.Numerical experiments supporting the accuracy of efficiency of the method in one and two dimensions are presented.


page 1

page 2

page 3

page 4


High order pressure-based semi-implicit IMEX schemes for the 3D Navier-Stokes equations at all Mach numbers

This article aims at developing a high order pressure-based solver for t...

An implicit symplectic solver for high-precision long term integrations of the Solar System

Compared to other symplectic integrators (the Wisdom and Holman map and ...

Irksome: Automating Runge–Kutta time-stepping for finite element methods

While implicit Runge–Kutta methods possess high order accuracy and impor...

Implicit Adaptive Mesh Refinement for Dispersive Tsunami Propagation

We present an algorithm to solve the dispersive depth-averaged Serre-Gre...

Efficient implicit solvers for models of neuronal networks

We introduce economical versions of standard implicit ODE solvers that a...

Parallel Optimizations for the Hierarchical Poincaré-Steklov Scheme (HPS)

Parallel optimizations for the 2D Hierarchical Poincaré-Steklov (HPS) di...

Nonlinear p-multigrid preconditioner for implicit time integration of compressible Navier–Stokes equations

Within the framework of p-adaptive flux reconstruction, we aim to constr...

Please sign up or login with your details

Forgot password? Click here to reset