A Higher-order Hybridisable Discontinuous Galerkin IMEX method for the incompressible Euler equations

TL;DR

This paper introduces a high-order hybridisable discontinuous Galerkin IMEX method for the incompressible Euler equations, combining spatial and temporal discretizations to achieve efficient, accurate solutions with reduced linear system complexity.

Contribution

A novel high-order hybridisable discontinuous Galerkin IMEX method for incompressible Euler equations that reduces pressure-velocity coupling to a well-studied system with efficient multigrid solvers.

Findings

Demonstrates high accuracy in numerical tests

Achieves computational efficiency through reduced linear systems

Validates method on standard test cases

Abstract

The incompressible Euler equations are an important model system in computational fluid dynamics. Fast high-order methods for the solution of this time-dependent system of partial differential equations are of particular interest: due to their exponential convergence in the polynomial degree they can make efficient use of computational resources. To address this challenge we describe a novel timestepping method which combines a hybridised Discontinuous Galerkin method for the spatial discretisation with IMEX timestepping schemes, thus achieving high-order accuracy in both space and time. The computational bottleneck is the solution of a (block-) sparse linear system to compute updates to pressure and velocity at each stage of the IMEX integrator. Following Chorin's projection approach, this update of the velocity and pressure fields is split into two stages. As a result, the hybridised equation for the implicit pressure-velocity problem is reduced to the well-known system which arises in hybridised mixed formulations of the Poisson- or diffusion problem and for which efficient multigrid preconditioners have been developed. Splitting errors can be reduced systematically by embedding this update into a preconditioned Richardson iteration. The accuracy and efficiency of the new method is demonstrated numerically for two time-dependent testcases that have been previously studied in the literature.