Semi-implicit discontinuous Galerkin methods for the incompressible Navier-Stokes equations on adaptive staggered Cartesian grids

TL;DR

This paper introduces a high-order semi-implicit discontinuous Galerkin method on adaptive staggered grids for solving incompressible Navier-Stokes equations efficiently, achieving accuracy and stability in complex geometries.

Contribution

The paper presents the first semi-implicit DG scheme on space-time adaptive staggered grids for incompressible flows, with a symmetric positive-definite pressure system.

Findings

Achieves arbitrary high order accuracy in space.

Produces a symmetric positive-definite linear pressure system.

Successfully verified on complex 2D and 3D test problems.

Abstract

In this paper a new high order semi-implicit discontinuous Galerkin method (SI-DG) is presented for the solution of the incompressible Navier-Stokes equations on staggered space-time adaptive Cartesian grids (AMR) in two and three space-dimensions. The pressure is written in the form of piecewise polynomials on the main grid, which is dynamically adapted within a cell-by-cell AMR framework. According to the time dependent main grid, different face-based spatially staggered dual grids are defined for the piece-wise polynomials of the respective velocity components. Arbitrary high order of accuracy is achieved in space, while a very simple semi-implicit time discretization is obtained via an explicit discretization of the nonlinear convective terms, and an implicit discretization of the pressure gradient in the momentum equation and of the divergence of the velocity field in the continuity equation. The real advantages of the staggered grid arise in the solution of the Schur complement associated with the saddle point problem of the discretized incompressible Navier-Stokes equations, i.e. after substituting the discrete momentum equations into the discrete continuity equation. This leads to a linear system for only one unknown, the scalar pressure. Indeed, the resulting linear pressure system is shown to be symmetric and positive-definite. The new space-time adaptive staggered DG scheme has been thoroughly verified for a large set of non-trivial test problems in two and three space dimensions, for which analytical, numerical or experimental reference solutions exist. To the knowledge of the authors, this is the first staggered semi-implicit DG scheme for the incompressible Navier-Stokes equations on space-time adaptive meshes in two and three space dimensions.