High-Order Methods for Computational Fluid Dynamics: A Brief Review of Compact Differential Formulations on Unstructured Grids

Popular high-order schemes with compact stencils for Computational Fluid Dynamics (CFD) include Discontinuous Galerkin (DG), Spectral Difference (SD), and Spectral Volume (SV) methods. The recently proposed Flux Reconstruction (FR) approach or Correction Procedure using Reconstruction (CPR) is based on a differential formulation and provides a unifying framework for these high-order schemes. Here we present a brief review of recent developments for the FR/CPR schemes as well as some pacing items

De-Aliasing Through Over-Integration Applied to the Flux Reconstruction and Discontinuous Galerkin Methods

High-order methods are quickly becoming popular for turbulent flows as the amount of computer processing power increases. The flux reconstruction (FR) method presents a unifying framework for a wide class of high-order methods including discontinuous Galerkin (DG), Spectral Difference (SD), and Spectral Volume (SV). It offers a simple, efficient, and easy way to implement nodal-based methods that are derived via the differential form of the governing equations. Whereas high-order methods have enjoyed recent success, they have been known to introduce numerical instabilities due to polynomial aliasing when applied to under-resolved nonlinear problems. Aliasing errors have been extensively studied in reference to DG methods; however, their study regarding FR methods has mostly been limited to the selection of the nodal points used within each cell. Here, we extend some of the de-aliasing techniques used for DG methods, primarily over-integration, to the FR framework. Our results show that over-integration does remove aliasing errors but may not remove all instabilities caused by insufficient resolution (for FR as well as DG)

A free-energy stable nodal discontinuous Galerkin approximation with summation-by-parts property for the Cahn-Hilliard equation

We present a nodal Discontinuous Galerkin (DG) scheme for the Cahn-Hilliard equation that satisfies the summation-by-parts simultaneous-approximation-term (SBP-SAT) property. The latter permits us to show that the discrete free-energy is bounded, and as a result, the scheme is provably stable. The scheme and the stability proof are presented for general curvilinear three-dimensional hexahedral meshes. We use the Bassi-Rebay 1 (BR1) scheme to compute interface fluxes, and an IMplicit-EXplicit (IMEX) scheme to integrate in time. Lastly, we test the theoretical findings numerically and present examples for two and three-dimensional problems

A comparative study on polynomial dealiasing and split form discontinuous Galerkin schemes for under-resolved turbulence computations

This work focuses on the accuracy and stability of high-order nodal discontinuous Galerkin (DG) methods for under-resolved turbulence computations. In particular we consider the inviscid Taylor-Green vortex (TGV) flow to analyse the implicit large eddy simulation (iLES) capabilities of DG methods at very high Reynolds numbers. The governing equations are discretised in two ways in order to suppress aliasing errors introduced into the discrete variational forms due to the under-integration of non-linear terms. The first, more straightforward way relies on consistent/over-integration, where quadrature accuracy is improved by using a larger number of integration points, consistent with the degree of the non-linearities. The second strategy, originally applied in the high-order finite difference community, relies on a split (or skew-symmetric) form of the governing equations. Different split forms are available depending on how the variables in the non-linear terms are grouped. The desired split form is then built by averaging conservative and non-conservative forms of the governing equations, although conservativity of the DG scheme is fully preserved. A preliminary analysis based on Burgers’ turbulence in one spatial dimension is conducted and shows the potential of split forms in keeping the energy of higher-order polynomial modes close to the expected levels. This indicates that the favourable dealiasing properties observed from split-form approaches in more classical schemes seem to hold for DG. The remainder of the study considers a comprehensive set of (under-resolved) computations of the inviscid TGV flow and compares the accuracy and robustness of consistent/over-integration and split form discretisations based on the local Lax-Friedrichs and Roe-type Riemann solvers. Recent works showed that relevant split forms can stabilize higher-order inviscid TGV test cases otherwise unstable even with consistent integration. Here we show that stable high-order cases achievable with both strategies have comparable accuracy, further supporting the good dealiasing properties of split form DG. The higher-order cases achieved only with split form schemes also displayed all the main features expected from consistent/over-integration. Among test cases with the same number of degrees of freedom, best solution quality is obtained with Roe-type fluxes at moderately high orders (around sixth order). Solutions obtained with very high polynomial orders displayed spurious features attributed to a sharper dissipation in wavenumber space. Accuracy differences between the two dealiasing strategies considered were, however, observed for the low-order cases, which also yielded reduced solution quality compared to high-order results

Comparison of different entropy stabilization techniques for discontinuous Galerkin spectral element methods

We review and compare two techniques to get entropy stability for nodal Discontinuous Galerkin Spectral Element Methods (DG) in compressible flows. One technique is based on entropy split-forms, e.g., [8, 15, 5] and one is based on a direct algebraic correction [1]. We have implemented the Flux Differencing methodology for both, Legendre-Gauss-Lobatto (Lobatto) and Legendre-Gauss (Gauss) based spectral element basis functions. While the Lobatto operators belong to the class of diagonal norm summation-by-parts (SBP) operators, the Gauss operators belong to the generalized class of SBP operators, where it is not necessary that the boundary nodes are included. To reach entropy stability, respectively guaranteed entropy dissipation, a key ingredient is an entropy conserving numerical flux function. With this ingredient, only the volume integral term of the DG method has to be modified accordingly. We have also implemented an alternative technique, which is in general applicable for a wide range of discretizations. Abgrall [R. Abgrall, A general framework to construct schemes satisfying additional conservation relations. Application to entropy conservative and entropy dissipative schemes. Journal of Computational Physics, vol 372, 2018] introduced an algebraic correction term that retains conservation of the primary quantities and is furthermore constructed such, that an entropy (in-) equality can be shown. The second technique is at first sight a simpler alternative to the split-form based approach. Hence, questions regarding its advantages and disadvantages naturally come up