A stable interface element scheme for the p-adaptive lifting collocation penalty formulation

This paper presents a procedure for adaptive polynomial refinement in the context of the lifting collocation penalty (LCP) formulation. The LCP scheme is a high-order unstructured discretization method unifying the discontinuous Galerkin, spectral volume, and spectral difference schemes in single differential formulation. Due to the differential nature of the scheme, the treatment of inter-cell fluxes for spatially varying polynomial approximations is not straightforward. Specially designed elements are proposed to tackle non-conforming polynomial approximations. These elements are constructed such that a conforming interface between polynomial approximations of different degrees is recovered. The stability and conservation properties of the scheme are analyzed and various inviscid compressible flow calculations are performed to demonstrate the potential of the proposed approach.     

求解多维欧拉方程的二阶非结构网格混合旋转Riemann求解器

将基于旋转近似Riemann求解器的二阶精度迎风型有限体积方法推广到非结构网格,采用基于网格中心的有限体积法,梯度的计算采用基于节点的方法引入更多的控制体模板,限制器的构造采用与非结构化网格相适应的形式.在求解Riemann问题时,沿具有一定物理意义的两个迎风方向,即控制体界面两侧速度差矢量方向及与之正交的方向.能够完全消除基于Riemann求解器的通量差分裂格式存在的激波不稳定或"红斑"现象.为减小计算量,采用HLL和Roe FDS混合旋转格式.     

Wavenumber-extended high-order oscillation control finite volume schemes for multi-dimensional aeroacoustic computations

A new numerical method toward accurate and efficient aeroacoustic computations of multi-dimensional compressible flows has been developed. The core idea of the developed scheme is to unite the advantages of the wavenumber-extended optimized scheme and M-AUSMPW+/MLP schemes by predicting a physical distribution of flow variables more accurately in multi-space dimensions. The wavenumber-extended optimization procedure for the finite volume approach based on the conservative requirement is newly proposed for accuracy enhancement, which is required to capture the acoustic portion of the solution in the smooth region. Furthermore, the new distinguishing mechanism which is based on the Gibbs phenomenon in discontinuity, between continuous and discontinuous regions is introduced to eliminate the excessive numerical dissipation in the continuous region by the restricted application of MLP according to the decision of the distinguishing function. To investigate the effectiveness of the developed method, a sequence of benchmark simulations such as spherical wave propagation, nonlinear wave propagation, shock tube problem and vortex preservation test problem are executed. Also, throughout more realistic shock–vortex interaction and muzzle blast flow problems, the utility of the new method for aeroacoustic applications is verified by comparing with the previous numerical or experimental results.     

A Switch Function-Based Gas-Kinetic Scheme for Simulation of Inviscid and Viscous Compressible Flows

In this paper, a switch function-based gas-kinetic scheme (SF-GKS) is presented for the simulation of inviscid and viscous compressible flows. With the finite volume discretization, Euler and Navier-Stokes equations are solved and the SF-GKS is applied to evaluate the inviscid flux at cell interface. The viscous flux is obtained by the conventional smooth function approximation. Unlike the traditional gas-kinetic scheme in the calculation of inviscid flux such as Kinetic Flux Vector Splitting (KFVS), the numerical dissipation is controlled with a switch function in the present scheme. That is, the numerical dissipation is only introduced in the region around strong shock waves. As a consequence, the present SF-GKS can well capture strong shock waves and thin boundary layers simultaneously. The present SF-GKS is firstly validated by its application to the inviscid flow problems, including 1-D Euler shock tube, regular shock reflection and double Mach reflection. Then, SF-GKS is extended to solve viscous transonic and hypersonic flow problems. Good agreement between the present results and those in the literature verifies the accuracy and robustness of SF-GKS.     

Three-dimensional non-free-parameter lattice-Boltzmann model and its application to inviscid compressible flows

In this Letter, a three-dimensional (3D) lattice-Boltzmann model is presented following the non-free-parameter lattice-Boltzmann method of Qu et al. [K. Qu, C. Shu, Y.T. Chew, Phys. Rev. E 75 (2007) 036706]. A simple function, which satisfies the zeroth- through third-order moments of the Maxwellian distribution function, is introduced to replace the Maxwellian distribution function as the continuous equilibrium distribution function in 3D space. The function is then discretized to discrete-velocity directions via a 25-point Lagrangian interpolation polynomial. To simulate compressible flows with shock waves, an implicit-explicit finite-difference scheme based on the total variation diminishing flux limitation is adopted to solve the discrete Boltzmann-BGK equation in order to capture the shock waves in compressible flows with a finite number of grid points. The model is validated by its application to some typical inviscid compressible flows ranging from 1D to 3D, and the numerical results are found to be in excellent agreement with the analytical solutions and/or other numerical results.     

On large time step TVD scheme for hyperbolic conservation laws and its efficiency evaluation

A large time step (LTS) TVD scheme originally proposed by Harten is modified and further developed in the present paper and applied to Euler equations in multidimensional problems. By firstly revealing the drawbacks of Harten’s original LTS TVD scheme, and reasoning the occurrence of the spurious oscillations, a modified formulation of its characteristic transformation is proposed and a high resolution, strongly robust LTS TVD scheme is formulated. The modified scheme is proven to be capable of taking larger number of time steps than the original one. Following the modified strategy, the LTS TVD schemes for Yee’s upwind TVD scheme and Yee–Roe–Davis’s symmetric TVD scheme are constructed. The family of the LTS schemes is then extended to multidimensional by time splitting procedure, and the associated boundary condition treatment suitable for the LTS scheme is also imposed. The numerical experiments on Sod’s shock tube problem, inviscid flows over NACA0012 airfoil and ONERA M6 wing are performed to validate the developed schemes. Computational efficiencies for the respective schemes under different CFL numbers are also evaluated and compared. The results reveal that the improvement is sizable as compared to the respective single time step schemes, especially for the CFL number ranging from 1.0 to 4.0.     

High order conservative differencing for viscous terms and the application to vortex-induced vibration flows

A new set of conservative 4th-order central finite differencing schemes for all the viscous terms of compressible Navier–Stokes equations are proposed and proved in this paper. These schemes are used with a 5th-order WENO scheme for inviscid flux and the stencil width of the central differencing scheme is designed to be within that of the WENO scheme. The central differencing schemes achieve the maximum order of accuracy in the stencil. This feature is important to keep the compactness of the overall discretization schemes and facilitate the boundary condition treatment. The algorithm is used to simulate the vortex-induced oscillations of an elastically mounted circular cylinder. The numerical results agree favorably with the experiment.     

Efficient kinetic schemes for steady and unsteady flow simulations on unstructured meshes

This paper presents efficient second-order kinetic schemes on unstructured meshes for both compressible unsteady and incompressible steady flows. For compressible unsteady flows, a time-dependent gas distribution function with a discontinuous particle velocity space at a cell interface is constructed and used for the evaluations of both numerical fluxes and conservative flow variables. As a result, a compact scheme on the unstructured meshes is developed. For incompressible steady flows, a continuous second-order gas-kinetic BGK type scheme is presented, for which the time-dependent gas distribution function with a continuous particle velocity is used on unstructured meshes. The efficiency of the schemes lies in the fact that the slopes of the flow variables inside each cell can be constructed using values of the flow variables within that cell only without involving neighboring cells. Therefore, even with the stencil of a first-order scheme, a high resolution method is constructed. Numerical examples are presented which are compared with the benchmark solutions and the experimental measurements.     

A discontinuous Galerkin method for inviscid low Mach number flows

In this work we extend the high-order discontinuous Galerkin (DG) finite element method to inviscid low Mach number flows. The method here presented is designed to improve the accuracy and efficiency of the solution at low Mach numbers using both explicit and implicit schemes for the temporal discretization of the compressible Euler equations. The algorithm is based on a classical preconditioning technique that in general entails modifying both the instationary term of the governing equations and the dissipative term of the numerical flux function (full preconditioning approach). In the paper we show that full preconditioning is beneficial for explicit time integration while the implicit scheme turns out to be efficient and accurate using just the modified numerical flux function. Thus the implicit scheme could also be used for time accurate computations. The performance of the method is demonstrated by solving an inviscid flow past a NACA0012 airfoil at different low Mach numbers using various degrees of polynomial approximations. Computations with and without preconditioning are performed on different grid topologies to analyze the influence of the spatial discretization on the accuracy of the DG solutions at low Mach numbers.     

A stable high-order Spectral Difference method for hyperbolic conservation laws on triangular elements

Numerical schemes using piecewise polynomial approximation are very popular for high order discretization of conservation laws. While the most widely used numerical scheme under this paradigm appears to be the Discontinuous Galerkin method, the Spectral Difference scheme has often been found attractive as well, because of its simplicity of formulation and implementation. However, recently it has been shown that the scheme is not linearly stable on triangles. In this paper we present an alternate formulation of the scheme, featuring a new flux interpolation technique using Raviart–Thomas spaces, which proves stable under a similar linear analysis in which the standard scheme failed. We demonstrate viability of the concept by showing linear stability both in the semi-discrete sense and for time stepping schemes of the SSP Runge–Kutta type. Furthermore, we present convergence studies, as well as case studies in compressible flow simulation using the Euler equations.     

On a new defect of shock-capturing methods

This paper proposes an explanation and a cure (or avoidance) to the new defect found of Eulerian shock-capturing methods in “A note on the conservative schemes for the Euler equations” by Tang and Liu [H. Tang, Tiegang Liu, A note on the conservative schemes for the Euler equations, J. Comput. Phys. 218 (2006) 451–459]. The latter gives a numerical investigation using several popular high resolution conservative schemes applied to Riemann problems of inviscid, compressible, perfect gas flows in Eulerian and Lagrangian coordinates with an initial high density ratio as well as a high pressure ratio. The results show that these methods work very inefficiently when applied to such problems and may give inaccurate numerical results, especially in shock location (or speed), even with a very fine grid.We have found that in problems of this type a strong rarefaction wave (SRW) is present adjacent to a contact line. Godunov averaging over the wave then produces large errors which, when the wave is strong, also persist for a long time. The cumulative error is thus very large which violates the strength of the contact line adjacent to it which, in turn, affects the speed and hence the location of the shock on the other side of the contact. We confirm this numerically using a method based on the unified coordinates with the shock-adaptive Godunov scheme plus contact strength preserving. The method, when applied to the Examples 2.1 and 2.2 of Tang and Liu [H. Tang, Tiegang Liu, A note on the conservative schemes for the Euler equations, J. Comput. Phys. 218 (2006) 451–459], produces high quality results even for comparatively coarse grids.     

三维高超声速无粘定常绕流的数值模拟

本文采用一种简单有效的通量分裂结合一种二阶TVD格式的数值通量的方法,提出一种隐式的迎风有限体积格式,并利用这种格式,从气体动力学非定常Euler方程组出发,数值模拟了三维不对称物体的高超声速无粘定常绕流。数值结果表明此格式具有分辨率较高和收敛速度较快的优点。     

新型上风格式(AUSM+)的性能分析及其应用

研究了新型上风格式AUSM+的分辨率、效率等性能,并用它与Roe、vanLeer上风格式数值模拟了前向台阶激波反射流动,通过对激波、膨胀波、接触间断及其间相互干扰的复杂波系的模拟对比,分析探讨了AUSM+格式的低数值耗散、间断高分辨率等特性。     

Shock-stable flux scheme for predicting aerodynamic heating load of hypersonic airliners

In the design of a hypersonic airliner that can considerably shorten the flight time, how to accurately predict the vehicle's aerodynamic heating loads is of great significance. In this study, a new shock-stable flux scheme called the simple low dissipation advection upwind splitting method(SLAU)-M1 is proposed for the prediction of hypersonic aerodynamic heating load. Based on the construction of the SLAU scheme for low-speed simulations, SLAU-M1 improves the robustness of the mass flux against shock instability. After validating the code employed, several numerical test cases are conducted. The onedimensional(1D) sod shock tube case and the two-dimensional(2D) inviscid NACA0012 airfoil case show that SLAU-M1 features a high level of accuracy at both low and high speeds. To simulate the hypersonic viscous flow over a blunt cone, we adopt different aspect ratios(ARs) of cells near the shock. The results suggest that SLAU-M1 is much less sensitive to the AR of cells near the shock in predicting hypersonic aerodynamic heating loads. Moreover, the findings show that the theoretical value is considerably better than that of the other schemes. The hypersonic viscous flow over a 2D double ellipsoid case and that over the Hypersonic Flight Experiment vehicle case also indicate that SLAU-M1 exhibits a considerably high level of accuracy in hypersonic heating predictions. These properties suggest that SLAU-M1 promises to be widely used in the accurate prediction of the aerodynamic heating loads of hypersonic airliners.     

Osher Flux with Entropy Fix for Two-Dimensional Euler Equations

We compare in this paper the properties of Osher flux with O-variant and P-variant (Osher-O flux and Osher-P flux) in finite volume methods for the two-dimensional Euler equations and propose an entropy fix technique to improve their robustness. We consider both first-order and second-order reconstructions. For inviscid hypersonic flow past a circular cylinder, we observe different problems for different schemes: a first-order Osher-O scheme on quadrangular grids yields a carbuncle shock, while a first-order Osher-P scheme results in a dislocation shock for high Mach number cases. In addition, a second-order Osher scheme can also yield a carbuncle shock or be unstable. To improve the robustness of these schemes we propose an entropy fix technique, and then present numerical results to show the effectiveness of the proposed method. In addition, the influence of grid aspects ratio, relative shock position to the grid and Mach number on shock stability are tested. Viscous heating problem and double Mach reflection problem are simulated to test the influence of the entropy fix on contact resolution and boundary layer resolution.     

VSIA:一种求解流体力学方程的多积分矩公式

使用变量间的积分矩公式给出一种构造流体力学方程数值格式的方法,即体积面积分平均方法(VSIA).该方法使用两种积分矩VIA(体积积分平均)和SIA(面积积分平均)给出一种完全的守恒型体积积分公式.基于一类守恒型半Lagrange输运算法CIP-CSL,能够清晰地构造出VSIA.对一般演化方程进行讨论,并应用于Burgers湍流、无粘可压流和不可压粘性流.     

Comparison of the generalized Riemann solver and the gas-kinetic scheme for inviscid compressible flow simulations

The generalized Riemann problem (GRP) scheme for the Euler equations and gas-kinetic scheme (GKS) for the Boltzmann equation are two high resolution shock capturing schemes for fluid simulations. The difference is that one is based on the characteristics of the inviscid Euler equations and their wave interactions, and the other is based on the particle transport and collisions. The similarity between them is that both methods can use identical MUSCL-type initial reconstructions around a cell interface, and the spatial slopes on both sides of a cell interface involve in the gas evolution process and the construction of a time-dependent flux function. Although both methods have been applied successfully to the inviscid compressible flow computations, their performances have never been compared. Since both methods use the same initial reconstruction, any difference is solely coming from different underlying mechanism in their flux evaluation. Therefore, such a comparison is important to help us to understand the correspondence between physical modeling and numerical performances. Since GRP is so faithfully solving the inviscid Euler equations, the comparison can be also used to show the validity of solving the Euler equations itself. The numerical comparison shows that the GRP exhibits a slightly better computational efficiency, and has comparable accuracy with GKS for the Euler solutions in 1D case, but the GKS is more robust than GRP. For the 2D high Mach number flow simulations, the GKS is absent from the shock instability and converges to the steady state solutions faster than the GRP. The GRP has carbuncle phenomena, likes a cloud hanging over exact Riemann solvers. The GRP and GKS use different physical processes to describe the flow motion starting from a discontinuity. One is based on the assumption of equilibrium state with infinite number of particle collisions, and the other starts from the non-equilibrium free transport process to evolve into an equilibrium one through particle collisions. The different mechanism in the flux evaluation deviates their numerical performance. Through this study, we may conclude scientifically that it may NOT be valid to use the Euler equations as governing equations to construct numerical fluxes in a discretized space with limited cell resolution. To adapt the Navier–Stokes (NS) equations is NOT valid either because the NS equations describe the flow behavior on the hydrodynamic scale and have no any corresponding physics starting from a discontinuity. This fact alludes to the consistency of the Euler and Navier–Stokes equations with the continuum assumption and the necessity of a direct modeling of the physical process in the discretized space in the construction of numerical scheme when modeling very high Mach number flows. The development of numerical algorithm is similar to the modeling process in deriving the governing equations, but the control volume here cannot be shrunk to zero.     

A Residual-Based Compact Scheme for the Compressible Navier–Stokes Equations

A simple and efficient time-dependent method is presented for solving the steady compressible Euler and Navier–Stokes equations with third-order accuracy. Owing to its residual-based structure, the numerical scheme is compact without requiring any linear algebra, and it uses a simple numerical dissipation built on the residual. The method contains no tuning parameter. Accuracy and efficiency are demonstrated for 2-D inviscid and viscous model problems. Navier–Stokes calculations are presented for a shock/boundary layer interaction, a separated laminar flow, and a transonic turbulent flow over an airfoil.     

The theoretical research of basic function method in incompressible viscous flow and its simulations in three-dimensional aneurysms

Basic function method is developed to treat the incompressible viscous flow. Artificial compressibility coefficient, the technique of flux splitting method and the combination of central and upwind schemes are applied to construct the basic function scheme of trigonometric function type for solving three-dimensional incompressible Navier-Stokes equations numerically. To prove the method, flows in finite-length-pipe are calculated, the velocity and pressure distribution of which solved by our method quite coincide with the exact solutions of Poiseuille flow except in the areas of entrance and exit. After the method is proved elementary, the hemodynamics in two- and three-dimensional aneurysms is researched numerically by using the basic function method of trigonometric function type and unstructured grids generation technique. The distributions of velocity, pressure and shear force in steady flow of aneurysms are calculated, and the influence of the shape of the aneurysms on the hemodynamics is studied. Supported by the National Natural Foundation of China (Grant Nos. 40874077, 40504020, and 40536029) and the National Basic Research Program of China (Grant No. 2006CB806304)     

Numerical Stability Analysis of FDLBM

We analyze the numerical stability of Finite Difference Lattice Boltzmann Method (FDLBM) by means of von Neumann stability analysis. The stability boundary of the FDLBM depends on the BGK relaxation time, the CFL number, the mean flow velocity, and the wavenumber. As the BGK relaxation time is increased at constant CFL number, the stability of the central difference LB scheme may not be ensured. The limits of maximum stable velocity are obtained around 0.39, 0.43, and 0.43 for the central difference, for the explicit upwind difference, and for the semi-implicit upwind difference schemes, respectively. We derive artificial viscosities for every difference scheme and investigate their influence on numerical stability. The requirements for artificial viscosity is consistent with the conditions derived from von Neumann stability analysis. This analysis elucidates that the upwind difference schemes are suitable for simulation of high Reynolds number flows.