High-order Discontinuous Galerkin solver based on Jacobi polynomial expansion for compressible flows on unstructured meshes

TL;DR

This paper introduces a high-order Discontinuous Galerkin solver based on Jacobi polynomial expansion for simulating compressible flows on unstructured meshes, enhancing accuracy and robustness in shock capturing.

Contribution

The work develops a novel high-order DG solver using Jacobi polynomials for arbitrary element types, with improved shock-capturing and verified efficiency and accuracy.

Findings

Verified convergence rates for different orders.

Demonstrated solver's effectiveness on complex flow problems.

Showed high computational efficiency at high precision.

Abstract

Based on the Jacobi polynomial expansion, an arbitrary high-order Discontinuous Galerkin solver for compressible flows on unstructured meshes is proposed in the present work. First, we construct orthogonal polynomials for 2D and 3D isoparametric elements using the 1D Jacobi polynomials. We perform modal expansions of the state variables using the orthogonal polynomials, enabling arbitrary high-order spatial discretization of these variables. Subsequently, the discrete governing equations are derived by considering the orthogonality of the Euler equations' residuals and the test functions. On this basis, we develop a high-order Discontinuous Galerkin solver that supports various element types, including triangles, quadrilaterals, tetrahedra, hexahedra, etc. An improved shock-capturing scheme has been adopted to capture shock discontinuities within the flow field. The variable's gradients at the discontinuous elements are reconstructed by its adjacent elements, and the slope limiter is applied to modify the state variables, smoothing the state variables and enhancing the robustness of the solver. The convergence rates of solvers of different orders have been verified by a benchmark case, and the CPU costs are given to prove that high-precision algorithms have higher computational efficiency under the same error level. Finally, several two- and three-dimensional compressible fluid dynamics problems are studied, compared with literature and experimental results, the effectiveness and accuracy of the solver were verified.