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.     

Divergence‐free discontinuous Galerkin schemes for the Stokes equations and the MAC scheme

We show that recently studied discontinuous Galerkin discretizations in their lowest order version are very similar to the marker and cell (MAC) finite difference scheme. Indeed, applying a slight modification, the exact MAC scheme can be recovered. Therefore, the analysis applied to the DG methods applies to the MAC scheme as well and the DG methods provide a natural generalization of the MAC scheme to higher order and irregular meshes. Copyright © 2007 John Wiley & Sons, Ltd.     

Analysis of a discontinuous Galerkin approximation of the Stokes problem based on an artificial compressibility flux

In this work, we propose and analyse a discontinuous Galerkin (DG) method for the Stokes problem based on an artificial compressibility numerical flux. A crucial step in the definition of a DG method is the choice of the numerical fluxes, which affect both the accuracy and the order of convergence of the method. We propose here to treat the viscous and the inviscid terms separately. The former is discretized using the well‐known BRMPS method. For the latter, the problem is locally modified by adding an artificial compressibility term of the form (1/c²)(?p/?t) for the sole purpose of interface flux computation. The flux is obtained as the exact solution of a local Riemann problem. The analysis of the method extends the well‐established strategies for the DG discretization of the Laplacian to the resulting partially coercive problem. Copyright © 2007 John Wiley & Sons, Ltd.     

INTERPOLATION PERTURBATION METHOD FOR SOLVING THE BOUNDARY LAYER TYPE PROBLEMS

ＹｕａｎＹｉｗｕ（袁镒吾）（ＲｅｃｅｉｖｅｄＯｃｔ．２，１９９４；ＣｏｍｍｕｎｉｃａｔｅｄｂｙＣｈｉｅｎＷｅｉｚａｎｇ）ＩＮＴＥＲＰＯＬＡＴＩＯＮＰＥＲＴＵＲＢＡＴＩＯＮＭＥＴＨＯＤＦＯＲＳＯＬＶＩＮＧＴＨＥＢＯＵＮＤＡＲＹＬＡＹＥＲＴＹＰＥＰＲＯＢ...     

Numerical simulations of compressible mixing layers with a discontinuous Galerkin method

Discontinuous Galerkin(DG) method is known to have several advantages for flow simulations,in particular,in fiexible accuracy management and adaptability to mesh refinement. In the present work,the DG method is developed for numerical simulations of both temporally and spatially developing mixing layers. For the temporally developing mixing layer,both the instantaneous fiow field and time evolution of momentum thickness agree very well with the previous results. Shocklets are observed at higher convective Mach numbers and the vortex paring manner is changed for high compressibility. For the spatially developing mixing layer,large-scale coherent structures and self-similar behavior for mean profiles are investigated. The instantaneous fiow field for a three-dimensional compressible mixing layer is also reported,which shows the development of largescale coherent structures in the streamwise direction. All numerical results suggest that the DG method is effective in performing accurate numerical simulations for compressible shear fiows.     

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.     

An immersed discontinuous Galerkin method for compressible Navier-Stokes equations on unstructured meshes

We introduce an immersed high-order discontinuous Galerkin method for solving the compressible Navier-Stokes equations on non–boundary-fitted meshes. The flow equations are discretised with a mixed discontinuous Galerkin formulation and are advanced in time with an explicit time marching scheme. The discretisation meshes may contain simplicial (triangular or tetrahedral) elements of different sizes and need not be structured. On the discretisation mesh, the fluid domain boundary is represented with an implicit signed distance function. The cut-elements partially covered by the solid domain are integrated after tessellation with the marching triangle or tetrahedra algorithms. Two alternative techniques are introduced to overcome the excessive stable time step restrictions imposed by cut-elements. In the first approach, the cut-basis functions are replaced with the extrapolated basis functions from the nearest largest element. In the second approach, the cut-basis functions are simply scaled proportionally to the fraction of the cut-element covered by the solid. To achieve high-order accuracy, additional nodes are introduced on the element faces abutting the solid boundary. Subsequently, the faces are curved by projecting the introduced nodes to the boundary. The proposed approach is verified and validated with several two- and three-dimensional subsonic and hypersonic low Reynolds number flow applications, including the flow over a cylinder, a space capsule, and an aerospace vehicle.     

A discontinuous Galerkin method for the Navier–Stokes equations

The foundations of a new discontinuous Galerkin method for simulating compressible viscous flows with shocks on standard unstructured grids are presented in this paper. The new method is based on a discontinuous Galerkin formulation both for the advective and the diffusive contributions. High‐order accuracy is achieved by using a recently developed hierarchical spectral basis. This basis is formed by combining Jacobi polynomials of high‐order weights written in a new co‐ordinate system. It retains a tensor‐product property, and provides accurate numerical quadrature. The formulation is conservative, and monotonicity is enforced by appropriately lowering the basis order and performing h‐refinement around discontinuities. Convergence results are shown for analytical two‐ and three‐dimensional solutions of diffusion and Navier–Stokes equations that demonstrate exponential convergence of the new method, even for highly distorted elements. Flow simulations for subsonic, transonic and supersonic flows are also presented that demonstrate discretization flexibility using hp‐type refinement. Unlike other high‐order methods, the new method uses standard finite volume grids consisting of arbitrary triangulizations and tetrahedrizations. Copyright © 1999 John Wiley & Sons, Ltd.     

Evaluation of discontinuous Galerkin and spectral volume methods for scalar and system conservation laws on unstructured grids

The discontinuous Galerkin (DG) and spectral volume (SV) methods are two recently developed high‐order methods for hyperbolic conservation laws capable of handling unstructured grids. In this paper, their overall performance in terms of efficiency, accuracy and memory requirement is evaluated using a 2D scalar conservation laws and the 2D Euler equations. To measure their accuracy, problems with analytical solutions are used. Both methods are also used to solve problems with strong discontinuities to test their ability in discontinuity capturing. Both the DG and SV methods are capable of achieving their formal order of accuracy while the DG method has a lower error magnitude and takes more memory. They are also similar in efficiency. The SV method appears to have a higher resolution for discontinuities because the data limiting can be done at the sub‐element level. Copyright © 2004 John Wiley & Sons, Ltd.     

On the discontinuous Galerkin computation of N‐waves focusing

Sonic boom focusing phenomenon can be predicted using the solution to the nonlinear Tricomi equation which is a hybrid (hyperbolic‐elliptic) second‐order partial differential equation. In this paper, the hyperbolic conservation law form is derived, which is valid in the entire domain. In this manner, the presence of two regions where the equation behaves differently (hyperbolic in the upper and elliptic in the lower half‐plane) is avoided. On the upper boundary, a new mixed boundary condition for the acoustic pressure is employed. The discretization is carried out using a discontinuous Galerkin (DG) method combined with a Runge–Kutta total‐variation diminishing scheme. The results show the accuracy of DG methods to solve problems involving sharp gradients and discontinuities. Comparisons with analytical results for the linear case, and other numerical results using classical explicit and compact finite difference schemes and weighted essentially non‐oscillatory schemes are included. Copyright © 2010 John Wiley & Sons, Ltd.     

A moving discontinuous Galerkin finite element method for flows with interfaces

A moving discontinuous Galerkin finite element method with interface condition enforcement is formulated for flows with discontinuous interfaces. The underlying weak formulation enforces the interface condition separately from the conservation law, so that the residual only vanishes upon satisfaction of both. In this formulation, the discrete grid geometry is treated as a variable, so that, in contrast to the standard discontinuous Galerkin method, this method has both the means to detect interfaces, via interface condition enforcement, and to satisfy, via grid movement, the conservation law and its associated interface condition. The method therefore directly fits interfaces, including shocks, preserving a high-order representation up to the interface without requiring shock capturing or an upwind numerical flux to achieve stability. It can be generalized to flows with a priori unknown interfaces with nontrivial topology and curved interface geometry as well as to an arbitrary number of spatial dimensions. Unsteady flows are represented in a manner similar to steady flows using a space-time formulation. In addition to computing flows with interfaces, the method can represent point singularities in a flow field by degenerating cuboid elements. In general, the method works in conjunction with standard local grid operations, including edge collapse, to ensure that degenerate cells are removed. Test cases are presented for up to three-dimensional flows that provide an initial assessment of the stability and accuracy of the method.     

An efficient discontinuous Galerkin method for aeroacoustic propagation

An efficient discontinuous Galerkin formulation is applied to the solution of the linearized Euler equations and the acoustic perturbation equations for the simulation of aeroacoustic propagation in two‐dimensional and axisymmetric problems, with triangular and quadrilateral elements. To improve computational efficiency, a new strategy of variable interpolation order is proposed in addition to a quadrature‐free approach and parallel implementation. Moreover, an accurate wall boundary condition is formulated on the basis of the solution of the Riemann problem for a reflective wall. Time discretization is based on a low dissipation formulation of a fourth‐order, low storage Runge–Kutta scheme. Along the far‐field boundaries a perfectly matched layer boundary condition is used. For the far‐field computations, the integral formulation of Ffowcs Williams and Hawkings is coupled with the near‐field solver. The efficiency and accuracy of the proposed variable order formulation is assessed for realistic geometries, namely sound propagation around a high‐lift airfoil and the Munt problem. Copyright © 2011 John Wiley & Sons, Ltd.     

A discontinuous Galerkin scheme for the numerical solution of flow problems with discontinuities

In this paper, a new discontinuous Galerkin finite element method for the numerical solution of flow problems with discontinuities is presented. The method is based on the limitation in every cell of the difference between the extrema values and the mean value of the numerical solution. The algorithm and technical details for the implementation of the method are presented in one‐and two‐dimensional problems. Numerical experiments for classical test problems are solved on unstructured triangulations to demonstrate the performance of the proposed method. Copyright © 2011 John Wiley & Sons, Ltd.     

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.     

Assessment of a discontinuous Galerkin method for the simulation of vortical flows at high Reynolds number

This paper focuses on the assessment of a discontinuous Galerkin method for the simulation of vortical flows at high Reynolds number. The Taylor–Green vortex at Re = 1600 is considered. The results are compared with those obtained using a pseudo‐spectral solver, converged on a 512³ grid and taken as the reference. The temporal evolution of the dissipation rate, visualisations of the vortical structures and the kinetic energy spectrum at the instant of maximal dissipation are compared to assess the results. At an effective resolution of 288³, the fourth‐order accurate discontinuous Galerkin method (DGM) solution (p = 3) is already very close to the pseudo‐spectral reference; the error on the dissipation rate is then essentially less than a percent, and the vorticity contours at times around the dissipation peak overlap everywhere. At a resolution of 384³, the solutions are indistinguishable. Then, an order convergence study is performed on the slightly under‐resolved grid (resolution of 192³). From the fourth order, the decrease of the error is no longer significant when going to a higher order. The fourth‐order DGM is also compared with an energy conserving fourth‐order finite difference method (FD4). The results show that, for the same number of DOF and the same order of accuracy, the errors of the DGM computation are significantly smaller. In particular, it takes 768³ DOF to converge the FD4 solution. Finally, the method is also successfully applied on unstructured high quality meshes. It is found that the dissipation rate captured is not significantly impacted by the element type. However, the element type impacts the energy spectrum in the large wavenumber range and thus the small vortical structures. In particular, at the same resolution, the results obtained using a tetrahedral mesh are much noisier than those obtained using a hexahedral mesh. Those obtained using a prismatic mesh are already much better, yet still slightly noisier. Copyright © 2013 John Wiley & Sons, Ltd.     

适用于间断Galerkin方法的限制器研究

开发了一种适用于高精度间断Galerkin方法的斜率（多项式系数）限制器。与现有的斜率限制器不同,该限制器实施过程不考虑网格单元类型（三角形或四边形）,通过全微分方法构造新的多项式系数,因此,该限制器能够适用于各种类型网格——结构化网格、具有单一单元的非结构化网格和具有混合单元的非结构化网格。由于该限制器能够方便地应用于具有混合单元的非结构化网格,因此,本文使用的程序能够方便地求解具有复杂几何结构的流动问题。本文利用一些典型算例对其性能进行了验证,表明该限制器适用于不同类型的网格单元,能够在光滑解区保证高的精度,并能够在阊断区抑帛3非物理振荡。     

Stability/dispersion analysis of the discontinuous Galerkin linearized shallow‐water system

The frequency or dispersion relation for the discontinuous Galerkin mixed formulation of the 1‐D linearized shallow‐water equations is analysed, using several basic DG mixed schemes. The dispersion properties are compared analytically and graphically with those of the mixed continuous Galerkin formulation for piecewise‐linear bases on co‐located grids. Unlike the Galerkin case, the DG scheme does not exhibit spurious stationary pressure modes. However, spurious propagating modes have been identified in all the present discontinuous Galerkin formulations. Numerical solutions of a test problem to simulate fast gravity modes illustrate the theoretical results and confirm the presence of spurious propagating modes in the DG schemes. Copyright © 2005 John Wiley & Sons, Ltd.     

An adaptive discretization of shallow‐water equations based on discontinuous Galerkin methods

In this paper, we present a discontinuous Galerkin formulation of the shallow‐water equations. An orthogonal basis is used for the spatial discretization and an explicit Runge–Kutta scheme is used for time discretization. Some results of second‐order anisotropic adaptive calculations are presented for dam breaking problems. The adaptive procedure uses an error indicator that concentrates the computational effort near discontinuities like hydraulic jumps. Copyright © 2006 John Wiley & Sons, Ltd.     

Singular perturbation of second-order nonlinear system with boundary perturbation

Ｗｅｃｏｎｓｉｄｅｒｉｎｔｈｉｓｐａｐｅｒｔｈｅｓｉｎｇｕｌａｒｐｅｒｔｕｒｂａｔｉｏｎｏｆｓｅｃｏｎｄ＿ｏｒｄｅｒｎｏｎｌｉｎｅａｒｓｙｓｔｅｍｉｎｖｏｌｖｉｎｇｉｎｔｅｒｇｒａｌｏｐｅｒａｔｏｒεｙ″=ｆ(ｔ,ｙ,Ｔｙ,ε)ｙ′ ｇ(ｔ,ｙ,Ｔｙ,ε),(1)ｗｉｔｈｂｏｕｎｄａｒｙｐｅｒｔｕｒｂａｔｉｏｎｙ(ｔ,ε)｜ｔ=φ(ε)=α(ε),ｙ(ｔ,ε)｜ｔ=1 ψ(ε)=β(ε),(2)ｗｈｅｒｅε>0ｉｓａｓｍａｌｌｐａｒａｍｅｔｅｒ,ａｎｄφ(ε),ψ(ε)ａｒｅｂｏｔｈ,ｗｉｔｈｒｅｓｐｅｃｔｔｏε,ｓｕｆｆｉｃｉｅｎｔｌｙｓｍｏ…