Implicit high-order discontinuous Galerkin method with HWENO type limiters for steady viscous flow simulations

Two types of implicit algorithms have been im- proved for high order discontinuous Galerkin （DG） method to solve compressible Navier-Stokes （NS） equations on tri- angular grids. A block lower-upper symmetric Gauss- Seidel （BLU-SGS） approach is implemented as a nonlin- ear iterative scheme. And a modified LU-SGS （LLU-SGS） approach is suggested to reduce the memory requirements while retain the good convergence performance of the origi- nal LU-SGS approach. Both implicit schemes have the sig- nificant advantage that only the diagonal block matrix is stored. The resulting implicit high-order DG methods are applied, in combination with Hermite weighted essentially non-oscillatory （HWENO） limiters, to solve viscous flow problems. Numerical results demonstrate that the present implicit methods are able to achieve significant efficiency improvements over explicit counterparts and for viscous flows with shocks, and the HWENO limiters can be used to achieve the desired essentially non-oscillatory shock tran- sition and the designed high-order accuracy simultaneously.

Research of influence of reduced-order boundary on accuracy and solution of interior points

The flow field with a high order scheme is usually calculated so as to solve complex flow problems and describe the flow structure accurately. However, there are two problems, i.e., the reduced-order boundary is inevitable and the order of the scheme at the discontinuous shock wave contained in the flow field as the supersonic flow field is low. It is questionable whether the reduced-order boundary and the low-order scheme at the shock wave have an effect on the numerical solution and accuracy of the flow field inside. In this paper, according to the actual situation of the direct numerical simulation of the flow field, two model equations with the exact solutions are solved, which are steady and unsteady, respectively, to study the question with a high order scheme at the interior of the domain and the reduced-order method at the boundary and center of the domain. Comparing with the exact solutions, it is found that the effect of reduced-order exists and cannot be ignored. In addition, the other two model equations with the exact solutions, which are often used in fluid mechanics, are also studied with the same process for the reduced-order problem.

High-Order Schemes Combining the Modified Equation Approach and Discontinuous Galerkin Approximations for the Wave Equation

We present a new high ordermethod in space and time for solving the wave equation,based on a newinterpretation of the“Modified Equation”technique.Indeed,contrary to most of the works,we consider the time discretization before the space discretization.After the time discretization,an additional biharmonic operator appears,which can not be discretized by classical finite elements.We propose a new Discontinuous Galerkinmethod for the discretization of this operator,andwe provide numerical experiments proving that the new method is more accurate than the classicalModified Equation technique with a lower computational burden.

High-order implicit discontinuous Galerkin schemes for unsteady compressible Navier–Stokes equations

Efficient solution techniques for high-order temporal and spatial discontinuous Galerkin（DG） discretizations of the unsteady Navier–Stokes equations are developed. A fourth-order implicit Runge–Kutta（IRK） scheme is applied for the time integration and a multigrid preconditioned GMRES solver is extended to solve the nonlinear system arising from each IRK stage. Several modifications to the implicit solver have been considered to achieve the efficiency enhancement and meantime to reduce the memory requirement. A variety of time-accurate viscous flow simulations are performed to assess the resulting high-order implicit DG methods. The designed order of accuracy for temporal discretization scheme is validate and the present implicit solver shows the superior performance by allowing quite large time step to be used in solving time-implicit systems. Numerical results are in good agreement with the published data and demonstrate the potential advantages of the high-order scheme in gaining both the high accuracy and the high efficiency.

High Order Compact Schemes in Projection Methods for Incompressible Viscous Flows

Within the projection schemes for the incompressible Navier-Stokes equations(namely"pressure-correction"method),we consider the simplest method(of order one in time)which takes into account the pressure in both steps of the splitting scheme.For this scheme,we construct,analyze and implement a new high order compact spatial approximation on nonstaggered grids.This approach yields a fourth order accuracy in space with an optimal treatment of the boundary conditions(without error on the velocity)which could be extended to more general splitting.We prove the unconditional stability of the associated Cauchy problem via von Neumann analysis.Then we carry out a normal mode analysis so as to obtain more precise results about the behavior of the numerical solutions.Finally we present detailed numerical tests for the Stokes and the Navier-Stokes equations(including the driven cavity benchmark)to illustrate the theoretical results.

Recent advances of computational aeroacoustics

Computational aeroacoustics （CAA） is an interdiscipline of aeroacoustics and computational fluid dynamics （CFD） for the investigation of sound generation and propagation from various aeroacoustics problems. In this review, the foundation and research scope of CAA are introduced firstly. A review of the early advances and applications of CAA is then briefly surveyed, focusing on two key issues, namely, high order finite difference scheme and non-reflecting boundary condition. Furthermore, the advances of CAA during the past five years are highlighted. Finally, the future prospective of CAA is briefly discussed.

High order semi-implicit weighted compact nonlinear scheme for the all-Mach isentropic Euler system

The computation of compressible flows at all Mach numbers is a very challenging problem.An efficient numerical method for solving this problem needs to have shock-capturing capability in the high Mach number regime,while it can deal with stiffness and accuracy in the low Mach number regime.This paper designs a high order semi-implicit weighted compact nonlinear scheme(WCNS)for the all-Mach isentropic Euler system of compressible gas dynamics.To avoid severe Courant-Friedrichs-Levy(CFL)restrictions for low Mach flows,the nonlinear fluxes in the Euler equations are split into stiff and non-stiff components.A third-order implicit-explicit(IMEX)method is used for the time discretization of the split components and a fifth-order WCNS is used for the spatial discretization of flux derivatives.The high order IMEX method is asymptotic preserving and asymptotically accurate in the zero Mach number limit.One-and two-dimensional numerical examples in both compressible and incompressible regimes are given to demonstrate the advantages of the designed IMEX WCNS.

Revisiting the space-time gradient method:A time-clocking perspective, high order difference time discretization and comparison with the harmonic balance method

This paper revisits the Space-Time Gradient(STG) method which was developed for efficient analysis of unsteady flows due to rotor–stator interaction and presents the method from an alternative time-clocking perspective. The STG method requires reordering of blade passages according to their relative clocking positions with respect to blades of an adjacent blade row. As the space-clocking is linked to an equivalent time-clocking, the passage reordering can be performed according to the alternative time-clocking. With the time-clocking perspective, unsteady flow solutions from different passages of the same blade row are mapped to flow solutions of the same passage at different time instants or phase angles. Accordingly, the time derivative of the unsteady flow equation is discretized in time directly, which is more natural than transforming the time derivative to a spatial one as with the original STG method. To improve the solution accuracy, a ninth order difference scheme has been investigated for discretizing the time derivative. To achieve a stable solution for the high order scheme, the implicit solution method of Lower-Upper Symmetric GaussSeidel/Gauss-Seidel(LU-SGS/GS) has been employed. The NASA Stage 35 and its blade-countreduced variant are used to demonstrate the validity of the time-clocking based passage reordering and the advantages of the high order difference scheme for the STG method. Results from an existing harmonic balance flow solver are also provided to contrast the two methods in terms of solution stability and computational cost.

Large-eddy simulation of wall-bounded turbulent flow with high-order discrete unified gas-kinetic scheme

In this paper,we introduce the discrete Maxwellian equilibrium distribution function for incompressible flow and force term into the two-stage third-order Discrete Unified Gas-Kinetic Scheme(DUGKS)for simulating low-speed turbulent flows.The Wall-Adapting Local Eddy-viscosity(WALE)and Vreman sub-grid models for Large-Eddy Simulations(LES)of turbulent flows are coupled within the present framework.Meanwhile,the implicit LES are also presented to verify the effect of LES models.A parallel implementation strategy for the present framework is developed,and three canonical wall-bounded turbulent flow cases are investigated,including the fully developed turbulent channel flow at a friction Reynolds number(Re)about 180,the turbulent plane Couette flow at a friction Re number about 93 and lid-driven cubical cavity flow at a Re number of 12000.The turbulence statistics,including mean velocity,the r.m.s.fluctuations velocity,Reynolds stress,etc.are computed by the present approach.Their predictions match precisely with each other,and they are both in reasonable agreement with the benchmark data of DNS.Especially,the predicted flow physics of three-dimensional lid-driven cavity flow are consistent with the description from abundant literature.The present numerical results verify that the present two-stage third-order DUGKS-based LES method is capable for simulating inhomogeneous wall-bounded turbulent flows and getting reliable results with relatively coarse grids.

Developing shock-capturing difference methods

A new shock-capturing method is proposed which is based on upwind schemes and flux-vector splittings. Firstly, original upwind schemes are projected along characteristic directions. Secondly, the amplitudes of the characteristic decompositions are carefully controlled by limiters to prevent non-physical oscillations. Lastly, the schemes are converted into conservative forms, and the oscillation-free shock-capturing schemes are acquired. Two explicit upwind schemes （2nd-order and 3rd-order） and three compact upwind schemes （3rd-order, 5th-order and 7th-order） are modified by the method for hyperbolic systems and the modified schemes are checked on several one-dimensional and two-dimensional test cases. Some numerical solutions of the schemes are compared with those of a WENO scheme and a MP scheme as well as a compact-WENO scheme. The results show that the method with high order accuracy and high resolutions can capture shock waves smoothly.

ON NUMERICAL TECHNIQUES IN CFD

Numerical techniques play an important role in CFD. Some of them are reviewed in this paper. The necessity of using high order difference scheme is demonstrated for the study of high Reynolds number viscous how. Physical guide lines are provided for the construction of these high order schemes. To avoid unduly ad hoc treatment in the boundary region the use of compact scheme is recommended because it has a small stencil size compared with the traditional finite difference scheme. Besides preliminary Fourier analysis shows the compact scheme can also yield better space resolution which makes it more suitable to study flow with multiscales e.g. turbulence. Other approaches such as perturbation method and finite spectral method are also emphasized. Typical numerical simulations were carried out. The first deals with Euler equations to show its capabilities to capture flow discontinuity. The second deals with Navier-Stokes Equations studying the evolution of a mixing layer, the pertinent structures at different times are shown. Asymmetric break down occurs and also the appearance of small vortices.

Fourth Order Difference Approximations for Space Riemann-Liouville Derivatives Based on Weighted and Shifted Lubich Difference Operators

High order discretization schemes playmore important role in fractional operators than classical ones.This is because usually for classical derivatives the stencil for high order discretization schemes is wider than low order ones;but for fractional operators the stencils for high order schemes and low order ones are the same.Then using high order schemes to solve fractional equations leads to almost the same computational cost with first order schemes but the accuracy is greatly improved.Using the fractional linear multistep methods,Lubich obtains the n-th order(n≤6)approximations of the a-th derivative(a>0)or integral(a<0)[Lubich,SIAM J.Math.Anal.,17,704-719,1986],because of the stability issue the obtained scheme can not be directly applied to the space fractional operator with a∈(1,2)for time dependent problem.By weighting and shifting Lubich’s 2nd order discretization scheme,in[Chen&Deng,SINUM,arXiv:1304.7425]we derive a series of effective high order discretizations for space fractional derivative,called WSLD operators there.As the sequel of the previous work,we further provide new high order schemes for space fractional derivatives by weighting and shifting Lubich’s 3rd and 4th order discretizations.In particular,we prove that the obtained 4th order approximations are effective for space fractional derivatives.And the corresponding schemes are used to solve the space fractional diffusion equation with variable coefficients.

Continuous Finite Element Subgrid Basis Functions for Discontinuous Galerkin Schemes on Unstructured Polygonal Voronoi Meshes

We propose a new high order accurate nodal discontinuous Galerkin(DG)method for the solution of nonlinear hyperbolic systems of partial differential equations(PDE)on unstructured polygonal Voronoi meshes.Rather than using classical polynomials of degree N inside each element,in our new approach the discrete solution is represented by piecewise continuous polynomials of degree N within each Voronoi element,using a continuous finite element basis defined on a subgrid inside each polygon.We call the resulting subgrid basis an agglomerated finite element(AFE)basis for the DG method on general polygons,since it is obtained by the agglomeration of the finite element basis functions associated with the subgrid triangles.The basis functions on each sub-triangle are defined,as usual,on a universal reference element,hence allowing to compute universal mass,flux and stiffness matrices for the subgrid triangles once and for all in a pre-processing stage for the reference element only.Consequently,the construction of an efficient quadrature-free algorithm is possible,despite the unstructured nature of the computational grid.High order of accuracy in time is achieved thanks to the ADER approach,making use of an element-local space-time Galerkin finite element predictor.The novel schemes are carefully validated against a set of typical benchmark problems for the compressible Euler and Navier-Stokes equations.The numerical results have been checked with reference solutions available in literature and also systematically compared,in terms of computational efficiency and accuracy,with those obtained by the corresponding modal DG version of the scheme.