A high-order semi-explicit discontinuous Galerkin solver for 3D incompressible flow with application to DNS and LES of turbulent channel flow

Abstract

We present an efficient discontinuous Galerkin scheme for simulation of the incompressible Navier-Stokes equations including laminar and turbulent flow. We consider a semi-explicit high-order velocity-correction method for time integration as well as nodal equal-order discretizations for velocity and pressure. The non-linear convective term is treated explicitly while a linear system is solved for the pressure Poisson equation and the viscous term. The key feature of our solver is a consistent penalty term reducing the local divergence error in order to overcome recently reported instabilities in spatially under-resolved high-Reynolds-number flows as well as small time steps. This penalty method is similar to the grad-div stabilization widely used in continuous finite elements. We further review and compare our method to several other techniques recently proposed in literature to stabilize the method for such flow configurations. The solver is specifically designed for large-scale computations through matrix-free linear solvers including efficient preconditioning strategies and tensor-product elements, which have allowed us to scale this code up to 34.4 billion degrees of freedom and 147,456 CPU cores. We validate our code and demonstrate optimal convergence rates with laminar flows present in a vortex problem and flow past a cylinder and show applicability of our solver to direct numerical simulation as well as implicit large-eddy simulation of turbulent channel flow at $R{e}_{\tau }=180$ as well as $590$.

28 pages, in preparation for submission to Journal of Computational Physics

References in corpus (3)

• On the stability of projection methods for the incompressible Navier-Stokes equations based on high-order discontinuous Galerkin discretizations
• Stabilised dG-FEM for incompressible natural convection flows with boundary and moving interior layers on non-adapted meshes
• Wall modeling via function enrichment within a high-order DG method for RANS simulations of incompressible flow

Cited by in corpus (19)

• Fast matrix-free evaluation of discontinuous Galerkin finite element operators
• On the stability of projection methods for the incompressible Navier-Stokes equations based on high-order discontinuous Galerkin discretizations
• Robust and efficient discontinuous Galerkin methods for under-resolved turbulent incompressible flows
• Efficiency of high-performance discontinuous Galerkin spectral element methods for under-resolved turbulent incompressible flows
• Hybrid multigrid methods for high-order discontinuous Galerkin discretizations
• High-order DG solvers for under-resolved turbulent incompressible flows: A comparison of $L^2$ and $H$(div) methods
• A matrix-free high-order discontinuous Galerkin compressible Navier-Stokes solver: A performance comparison of compressible and incompressible formulations for turbulent incompressible flows
• A weakly compressible hybridizable discontinuous Galerkin formulation for fluid-structure interaction problems
• The analogue of grad-div stabilization in DG methods for incompressible flows: Limiting behavior and extension to tensor-product meshes
• Wall modeling via function enrichment within a high-order DG method for RANS simulations of incompressible flow
• An interior penalty discontinuous Galerkin approach for 3D incompressible Navier--Stokes equation for permeability estimation of porous media
• High-order arbitrary Lagrangian-Eulerian discontinuous Galerkin methods for the incompressible Navier-Stokes equations
• A spectral deferred correction method for incompressible flow with variable viscosity
• Fluid-structure interaction with $H(\text{div})$-conforming finite elements
• A multiscale approach to hybrid RANS/LES wall modeling within a high-order discontinuous Galerkin scheme using function enrichment
• Wall modeling via function enrichment: extension to detached-eddy simulation
• Multigrid methods for the Stokes problem on GPU systems
• Modern discontinuous Galerkin methods for the simulation of transitional and turbulent flows in biomedical engineering: A comprehensive LES study of the FDA benchmark nozzle model
• A high-order discontinuous Galerkin pressure robust splitting scheme for incompressible flows