A simple shock‐capturing technique for high‐order discontinuous Galerkin methods

This article presents a novel shock‐capturing technique for the discontinuous Galerkin (DG) method. The technique is designed for compressible flow problems, which are usually characterized by the presence of strong shocks and discontinuities. The inherent structure of standard DG methods seems to suggest that they are especially adapted to capture shocks because of the numerical fluxes based on suitable approximate Riemann solvers, which, in practice, introduces some stabilization. However, the usual numerical fluxes are not sufficient to stabilize the solution in the presence of shocks for large high‐order elements. Here, a new basis of shape functions is introduced. It has the ability to change locally between a continuous or discontinuous interpolation depending on the smoothness of the approximated function. In the presence of shocks, the new discontinuities inside an element introduce the required stabilization because of numerical fluxes. Large high‐order elements can therefore be used and shocks captured within a single element, avoiding adaptive mesh refinement and preserving the locality and compactness of the DG scheme. Several numerical examples for transonic and supersonic flows are studied to demonstrate the applicability of the proposed approach. Copyright © 2011 John Wiley & Sons, Ltd.     

Higher‐order and adaptive discontinuous Galerkin methods with shock‐capturing applied to transonic turbulent delta wing flow

Discontinuous Galerkin (DG) methods allow high‐order flow solutions on unstructured or locally refined meshes by increasing the polynomial degree and using curved instead of straight‐sided elements. DG discretizations with higher polynomial degrees must, however, be stabilized in the vicinity of discontinuities of flow solutions such as shocks. In this article, we device a consistent shock‐capturing method for the Reynolds‐averaged Navier–Stokes and k‐ ω turbulence model equations based on an artificial viscosity term that depends on element residual terms. Furthermore, the DG method is combined with a residual‐based adaptation algorithm that targets at resolving all flow features. The higher‐order and adaptive DG method is applied to a fully turbulent transonic flow around the second Vortex Flow Experiment (VFE‐2) configuration with a good resolution of the vortex system.Copyright © 2013 John Wiley & Sons, Ltd.     

High‐order continuous and discontinuous Galerkin methods for wave problems

Three Galerkin methods—continuous Galerkin, Compact Discontinuous Galerkin, and hybridizable discontinuous Galerkin—are compared in terms of performance and computational efficiency in 2‐D scattering problems for low and high‐order polynomial approximations. The total number of DOFs and the total runtime are used for this correlation as well as the corresponding precision. The comparison is carried out through various numerical examples. The superior performance of high‐order elements is shown. At the same time, similar capabilities are shown for continuous Galerkin and hybridizable discontinuous Galerkin, when high‐order elements are adopted, both of them clearly outperforming compact discontinuous Galerkin. Copyright © 2013 John Wiley & Sons, Ltd.     

Assessment of shock capturing schemes for discontinuous Galerkin method

This paper carries out systematical investigations on the performance of several typical shock-capturing schemes for the discontinuous Galerkin （DG） method, including the total variation bounded （TVB） limiter and three artificial diffusivity schemes （the basis function-based （BF） scheme, the face residual-based （FR） scheme, and the element residual-based （ER） scheme）. Shock-dominated flows （the Sod problem, the Shu- Osher problem, the double Mach reflection problem, and the transonic NACA0012 flow） are considered, addressing the issues of accuracy, non-oscillatory property, dependence on user-specified constants, resolution of discontinuities, and capability for steady solutions. Numerical results indicate that the TVB limiter is more efficient and robust, while the artificial diffusivity schemes are able to preserve small-scale flow structures better. In high order cases, the artificial diffusivity schemes have demonstrated superior performance over the TVB limiter.     

Revisit of dilation-based shock capturing for discontinuous Galerkin methods

The idea of using velocity dilation for shock capturing is revisited in this paper, combined with the discontinuous Galerkin method. The value of artificial viscosity is determined using direct dilation instead of its higher order derivatives to reduce cost and degree of difficulty in computing derivatives. Alternative methods for estimating the element size of large aspect ratio and smooth artificial viscosity are proposed to further improve robustness and accuracy of the model. Several benchmark tests are conducted, ranging from subsonic to hypersonic flows involving strong shocks. Instead of adjusting empirical parameters to achieve optimum results for each case, all tests use a constant parameter for the model with reasonable success, indicating excellent robustness of the method. The model is only limited to third-order accuracy for smooth flows. This limitation may be relaxed by using a switch or a wall function. Overall, the model is a good candidate for compressible flows with potentials of further improvement.     

A high‐order discontinuous Galerkin method for extension problems

We present a novel technique for solving extension problems such as the extension velocity, by reformulating the problem into an elliptic differential equation. We introduce a novel discretization using an upwind flux without any additional stabilization. This leads to a triangular matrix structure, which can be solved using a marching algorithm and high‐order accuracy, even in the presence of singularities.     

Boundary discretization for high‐order discontinuous Galerkin computations of tidal flows around shallow water islands

In this paper some preliminary results concerning the application of the high‐order discontinuous Galerkin (DG) method for the resolution of realistic problems of tidal flows around shallow water islands are presented. In particular, tidal flows are computed around the Rattray island located in the Great Barrier Reef. This island is a standard benchmark problem well documented in the literature providing useful in situ measurements for validation of the model. Realistic elements of the simulation are a tidal flow forcing, a variable bathymetry and a non‐trivial coastline. The computation of tidal flows in shallow water around an island is very similar to the simulation of the Euler equations around bluff bodies in quasi‐steady flows. The main difference lies in the high irregularity of islands' shapes and in the fact that, in the framework of large‐scale ocean models, the number of elements to represent an island is drastically limited compared with classical engineering computations. We observe that the high‐order DG method applied to shallow water flows around bluff bodies with poor linear boundary representations produces oscillations and spurious eddies. Surprisingly those eddies may have the right size and intensity but may be generated by numerical diffusion and are not always mathematically relevant. Although not interested in solving accurately the boundary layers of an island, we show that a high‐order boundary representation is mandatory to avoid non‐physical eddies and spurious oscillations. It is then possible to parametrize accurately the subgrid‐scale processes to introduce the correct amount of diffusion in the model. The DG results around the Rattray island are eventually compared with current measurements and reveal good agreement. Copyright © 2008 John Wiley & Sons, Ltd.     

Evaluation of some high‐order shock capturing schemes for direct numerical simulation of unsteady two‐dimensional free flows

The present study addresses the capability of a large set of shock‐capturing schemes to recover the basic interactions between acoustic, vorticity and entropy in a direct numerical simulation (DNS) framework. The basic dispersive and dissipative errors are first evaluated by considering the advection of a Taylor vortex in a uniform flow. Two transonic cases are also considered. The first one consists of the interaction between a temperature spot and a weak shock. This test emphasizes the capability of the schemes to recover the production of vorticity through the baroclinic process. The second one consists of the interaction of a Taylor vortex with a weak shock, corresponding to the framework of the linear theory of Ribner. The main process in play here is the production of an acoustic wave. The results obtained by using essentially non‐oscillatory (ENO), total variation diminishing (TVD), compact‐TVD and MUSCL schemes are compared with those obtained by means of a sixth‐order accurate Hermitian scheme, considered as reference. The results are as follows; the ENO schemes agree pretty well with the reference scheme. The second‐order accurate Upwind‐TVD scheme exhibits a strong numerical diffusion, while the MUSCL scheme behavior is very sensitive to the value on the parameter β in the limiter function minmod. The compact‐TVD schemes do not yield improvement over the standard TVD schemes. Copyright © 2000 John Wiley & Sons, Ltd.     

Efficiency of high‐performance discontinuous Galerkin spectral element methods for under‐resolved turbulent incompressible flows

The present paper addresses the numerical solution of turbulent flows with high‐order discontinuous Galerkin methods for discretizing the incompressible Navier‐Stokes equations. The efficiency of high‐order methods when applied to under‐resolved problems is an open issue in the literature. This topic is carefully investigated in the present work by the example of the three‐dimensional Taylor‐Green vortex problem. Our implementation is based on a generic high‐performance framework for matrix‐free evaluation of finite element operators with one of the best realizations currently known. We present a methodology to systematically analyze the efficiency of the incompressible Navier‐Stokes solver for high polynomial degrees. Due to the absence of optimal rates of convergence in the under‐resolved regime, our results reveal that demonstrating improved efficiency of high‐order methods is a challenging task and that optimal computational complexity of solvers and preconditioners as well as matrix‐free implementations are necessary ingredients in achieving the goal of better solution quality at the same computational costs already for a geometrically simple problem such as the Taylor‐Green vortex. Although the analysis is performed for a Cartesian geometry, our approach is generic and can be applied to arbitrary geometries. We present excellent performance numbers on modern cache‐based computer architectures achieving a throughput for operator evaluation of 3·10⁸ up to 1·10⁹ DoFs/s (degrees of freedom per second) on one Intel Haswell node with 28 cores. Compared to performance results published within the last five years for high‐order discontinuous Galerkin discretizations of the compressible Navier‐Stokes equations, our approach reduces computational costs by more than one order of magnitude for the same setup.     

Numerical assessment of a shock‐detecting sensor for low dissipative high‐order simulation of shock–vortex interactions

Considering the importance of high‐order schemes implementation for the simulation of shock‐containing turbulent flows, the present work involves the assessment of a shock‐detecting sensor for filtering of high‐order compact finite‐difference schemes for simulation of this type of flows. To accomplish this, a sensor that controls the amount of numerical dissipation is applied to a sixth‐order compact scheme as well as a fourth‐order two‐register Runge–Kutta method for numerical simulation of various cases including inviscid and viscous shock–vortex and shock–mixing‐layer interactions. Detailed study is performed to investigate the performance of the sensor, that is, the effect of control parameters employed in the sensor are investigated in the long‐time integration. In addition, the effects of nonlinear weighting factors controlling the value of the second‐order and high‐order filters in fine and coarse non‐uniform grids are investigated. The results indicate the accuracy of the nonlinear filter along with the promising performance of the shock‐detecting sensor, which would pave the way for future simulations of turbulent flows containing shocks. Copyright © 2014 John Wiley & Sons, Ltd.     

Dilation‐based shock capturing for high‐order methods

In this paper, we introduce a shock‐capturing artificial viscosity technique for high‐order unstructured mesh methods. This artificial viscosity model is based on a non‐dimensional form of the divergence of the velocity. The technique is an extension and improvement of the dilation‐based artificial viscosity methods introduced in Premasuthan et al., 15 and further extended in Nguyen and Peraire 27 . The approach presented has a number attractive properties including non‐dimensional analytical form, sub‐cell resolution, and robustness for complex shock flows on anisotropic meshes. We present extensive numerical results to demonstrate the performance of the proposed approach. Copyright © 2016 John Wiley & Sons, Ltd.     

High‐order accurate p‐multigrid discontinuous Galerkin solution of the Euler equations

Discontinuous Galerkin (DG) methods have proven to be perfectly suited for the construction of very high‐order accurate numerical schemes on arbitrary unstructured and possibly nonconforming grids for a wide variety of applications, but are rather demanding in terms of computational resources. In order to improve the computational efficiency of this class of methods a p‐multigrid solution strategy has been developed, which is based on a semi‐implicit Runge–Kutta smoother for high‐order polynomial approximations and the implicit Backward Euler smoother for piecewise constant approximations. The effectiveness of the proposed approach is demonstrated by comparison with p‐multigrid schemes employing purely explicit smoothing operators for several 2D inviscid test cases. Copyright © 2008 John Wiley & Sons, Ltd.     

A high‐order discontinuous Galerkin method for all‐speed flows

In this article, we present a discontinuous Galerkin (DG) method designed to improve the accuracy and efficiency of steady solutions of the compressible fully coupled Reynolds‐averaged Navier–Stokes and k ? ω turbulence model equations for solving all‐speed flows. The system of equations is iterated to steady state by means of an implicit scheme. The DG solution is extended to the incompressible limit by implementing a low Mach number preconditioning technique. A full preconditioning approach is adopted, which modifies both the unsteady terms of the governing equations and the dissipative term of the numerical flux function by means of a new preconditioner, on the basis of a modified version of Turkel's preconditioning matrix. At sonic speed the preconditioner reduces to the identity matrix thus recovering the non‐preconditioned DG discretization. An artificial viscosity term is added to the DG discretized equations to stabilize the solution in the presence of shocks when piecewise approximations of order of accuracy higher than 1 are used. Moreover, several rescaling techniques are implemented in order to overcome ill‐conditioning problems that, in addition to the low Mach number stiffness, can limit the performance of the flow solver. These approaches, through a proper manipulation of the governing equations, reduce unbalances between residuals as a result of the dependence on the size of elements in the computational mesh and because of the inherent differences between turbulent and mean‐flow variables, influencing both the evolution of the Courant Friedrichs Lewy (CFL) number and the inexact solution of the linear systems. The performance of the method is demonstrated by solving three turbulent aerodynamic test cases: the flat plate, the L1T2 high‐lift configuration and the RAE2822 airfoil (Case 9). The computations are performed at different Mach numbers using various degrees of polynomial approximations to analyze the influence of the proposed numerical strategies on the accuracy, efficiency and robustness of a high‐order DG solver at different flow regimes. Copyright © 2014 John Wiley & Sons, Ltd.     

An improvement of classical slope limiters for high‐order discontinuous Galerkin method

In this paper, we describe some existing slope limiters (Cockburn and Shu's slope limiter and Hoteit's slope limiter) for the two‐dimensional Runge–Kutta discontinuous Galerkin (RKDG) method on arbitrary unstructured triangular grids. We describe the strategies for detecting discontinuities and for limiting spurious oscillations near such discontinuities, when solving hyperbolic systems of conservation laws by high‐order discontinuous Galerkin methods. The disadvantage of these slope limiters is that they depend on a positive constant, which is, for specific hydraulic problems, difficult to estimate in order to eliminate oscillations near discontinuities without decreasing the high‐order accuracy of the scheme in the smooth regions. We introduce the idea of a simple modification of Cockburn and Shu's slope limiter to avoid the use of this constant number. This modification consists in: slopes are limited so that the solution at the integration points is in the range spanned by the neighboring solution averages. Numerical results are presented for a nonlinear system: the shallow water equations. Four hydraulic problems of discontinuous solutions of two‐dimensional shallow water are presented. The idealized dam break problem, the oblique hydraulic jump problem, flow in a channel with concave bed and the dam break problem in a converging–diverging channel are solved by using the different slope limiters. Numerical comparisons on unstructured meshes show a superior accuracy with the modified slope limiter. Moreover, it does not require the choice of any constant number for the limiter condition. Copyright © 2008 John Wiley & Sons, Ltd.     

An hp‐adaptive discontinuous Galerkin method for shallow water flows

An adaptive spectral/hp discontinuous Galerkin method for the two‐dimensional shallow water equations is presented. The model uses an orthogonal modal basis of arbitrary polynomial order p defined on unstructured, possibly non‐conforming, triangular elements for the spatial discretization. Based on a simple error indicator constructed by the solutions of approximation order p and p?1, we allow both for the mesh size, h, and polynomial approximation order to dynamically change during the simulation. For the h‐type refinement, the parent element is subdivided into four similar sibling elements. The time‐stepping is performed using a third‐order Runge–Kutta scheme. The performance of the hp‐adaptivity is illustrated for several test cases. It is found that for the case of smooth flows, p‐adaptivity is more efficient than h‐adaptivity with respect to degrees of freedom and computational time. Copyright © 2010 John Wiley & Sons, Ltd.     

Connections between the discontinuous Galerkin method and high‐order flux reconstruction schemes

With high‐order methods becoming more widely adopted throughout the field of computational fluid dynamics, the development of new computationally efficient algorithms has increased tremendously in recent years. One of the most recent methods to be developed is the flux reconstruction approach, which allows various well‐known high‐order schemes to be cast within a single unifying framework. Whilst a connection between flux reconstruction and the more widely adopted discontinuous Galerkin method has been established elsewhere, it still remains to fully investigate the explicit connections between the many popular variants of the discontinuous Galerkin method and the flux reconstruction approach. In this work, we closely examine the connections between three nodal versions of tensor‐product discontinuous Galerkin spectral element approximations and two types of flux reconstruction schemes for solving systems of conservation laws on quadrilateral meshes. The different types of discontinuous Galerkin approximations arise from the choice of the solution nodes of the Lagrange basis representing the solution and from the quadrature approximation used to integrate the mass matrix and the other terms of the discretization. By considering both linear and nonlinear advection equations on a regular grid, we examine the mathematical properties that connect these discretizations. These arguments are further confirmed by the results of an empirical numerical study. Copyright © 2014 John Wiley & Sons, Ltd.     

Compressive advection and multi‐component methods for interface‐capturing

This paper develops methods for interface‐capturing in multiphase flows. The main novelties of these methods are as follows: (a) multi‐component modelling that embeds interface structures into the continuity equation; (b) a new family of triangle/tetrahedron finite elements, in particular, the P₁DG‐P₂(linear discontinuous between elements velocity and quadratic continuous pressure); (c) an interface‐capturing scheme based on compressive control volume advection methods and high‐order finite element interpolation methods; (d) a time stepping method that allows use of relatively large time step sizes; and (e) application of anisotropic mesh adaptivity to focus the numerical resolution around the interfaces and other areas of important dynamics. This modelling approach is applied to a series of pure advection problems with interfaces as well as to the simulation of the standard computational fluid dynamics benchmark test cases of a collapsing water column under gravitational forces (in two and three dimensions) and sloshing water in a tank. Two more test cases are undertaken in order to demonstrate the many‐material and compressibility modelling capabilities of the approach. Numerical simulations are performed on coarse unstructured meshes to demonstrate the potential of the methods described here to capture complex dynamics in multiphase flows. Copyright © 2015 John Wiley & Sons, Ltd.     

A high‐order discontinuous Galerkin solver for low Mach number flows

In this work, we present a high‐order discontinuous Galerkin method (DGM) for simulating variable density flows at low Mach numbers. The corresponding low Mach number equations are an approximation of the compressible Navier–Stokes equations in the limit of zero Mach number. To the best of the authors'y knowledge, it is the first time that the DGM is applied to the low Mach number equations. The mixed‐order formulation is applied for spatial discretization. For steady cases, we apply the semi‐implicit method for pressure‐linked equation (SIMPLE) algorithm to solve the non‐linear system in a segregated manner. For unsteady cases, the solver is implicit in time using backward differentiation formulae, and the SIMPLE algorithm is applied to solve the non‐linear system in each time step. Numerical results for the following three test cases are shown: Couette flow with a vertical temperature gradient, natural convection in a square cavity, and unsteady natural convection in a tall cavity. Considering a fixed number of degrees of freedom, the results demonstrate the benefits of using higher approximation orders. Copyright © 2015 John Wiley & Sons, Ltd.