气相爆轰高阶中心差分-WENO组合格式自适应网格方法

研究一种高阶中心差分-WENO组合格式,并采用自适应网格方法进行二维和三维气相爆轰波的数值模拟.采用ZND爆轰模型的控制方程为包含化学反应源项的Euler方程组.组合格式在大梯度区采用WENO格式捕捉间断,在光滑区采用高阶中心差分格式提高计算效率.采用一种基于流场结构特征的自适应网格.计算结果,表明这种方法同时具有高精度、高分辨率和高效率的特点.     

An adaptive central-upwind weighted essentially non-oscillatory scheme

In this work, an adaptive central-upwind 6th-order weighted essentially non-oscillatory (WENO) scheme is developed. The scheme adapts between central and upwind schemes smoothly by a new weighting relation based on blending the smoothness indicators of the optimal higher order stencil and the lower order upwind stencils. The scheme achieves 6th-order accuracy in smooth regions of the solution by introducing a new reference smoothness indicator. A number of numerical examples suggest that the present scheme, while preserving the good shock-capturing properties of the classical WENO schemes, achieves very small numerical dissipation.     

A Fifth-Order Low-Dissipative Conservative Upwind Compact Scheme Using Centered Stencil

In this paper, a conservative fifth-order upwind compact scheme using centered stencil is introduced. This scheme uses asymmetric coefficients to achieve the upwind property since the stencil is symmetric. Theoretical analysis shows that the proposed scheme is low-dissipative and has a relatively large stability range. To maintain the convergence rate of the whole spatial discretization, a proper non-periodic boundary scheme is also proposed. A detailed analysis shows that the spatial discretization implemented with the boundary scheme proposed by Pirozzoli [J. Comput. Phys., 178 (2001), pp. 81-117] is approximately fourth-order. Furthermore, a hybrid methodology, coupling the compact scheme with WENO scheme, is adopted for problems with discontinuities. Numerical results demonstrate the effectiveness of the proposed scheme.     

Generalized finite compact difference scheme for shock/complex flowfield interaction

A class of generalized high order finite compact difference schemes is proposed for shock/vortex, shock/boundary layer interaction problems. The finite compact difference scheme takes the region between two shocks as a compact stencil. The high order WENO fluxes on shock stencils are used as the internal boundary fluxes for the compact scheme. A lemma based on the property of smoothness estimators on a 5-points stencil is given to detect the shock position. There is no free parameter introduced to switch the compact scheme and the WENO scheme. Some numerical experiments are given and they demonstrate that the present scheme has low dissipation due to the compact central differencing scheme used in the smooth regions.     

A General-Purpose Finite-Volume Advection Scheme for Continuous and Discontinuous Fields on Unstructured Grids

We present a new general-purpose advection scheme for unstructured meshes based on the use of a variation of the interface-tracking flux formulation recently put forward by O. Ubbink and R. I. Issa (J. Comput. Phys.153, 26 (1999)), in combination with an extended version of the flux-limited advection scheme of J. Thuburn (J. Comput. Phys.123, 74 (1996)), for continuous fields. Thus, along with a high-order mode for continuous fields, the new scheme presented here includes optional integrated interface-tracking modes for discontinuous fields. In all modes, the method is conservative, monotonic, and compatible. It is also highly shape preserving. The scheme works on unstructured meshes composed of any kind of connectivity element, including triangular and quadrilateral elements in two dimensions and tetrahedral and hexahedral elements in three dimensions. The scheme is finite-volume based and is applicable to control-volume finite-element and edge-based node-centered computations. An explicit–implicit extension to the continuous-field scheme is provided only to allow for computations in which the local Courant number exceeds unity. The transition from the explicit mode to the implicit mode is performed locally and in a continuous fashion, providing a smooth hybrid explicit–implicit calculation. Results for a variety of test problems utilizing the continuous and discontinuous advection schemes are presented.     

A class of finite difference schemes with low dispersion and controllable dissipation for DNS of compressible turbulence

In this paper, a class of finite difference schemes which achieves low dispersion and controllable dissipation in smooth region and robust shock-capturing capabilities in the vicinity of discontinuities is presented. Firstly, a sufficient condition for semi-discrete finite difference schemes to have independent dispersion and dissipation is derived. This condition enables a novel approach to separately optimize the dissipation and dispersion properties of finite difference schemes and a class of schemes with minimized dispersion and controllable dissipation is thus obtained. Secondly, for the purpose of shock-capturing, one of these schemes is used as the linear part of the WENO scheme with symmetrical stencils to constructed an improved WENO scheme. At last, the improved WENO scheme is blended with its linear counterpart to form a new hybrid scheme for practical applications. The proposed scheme is accurate, flexible and robust. The accuracy and resolution of the proposed scheme are tested by the solutions of several benchmark test cases. The performance of this scheme is further demonstrated by its application in the direct numerical simulation of compressible turbulent channel flow between isothermal walls.     

An improved weighted essentially non-oscillatory scheme for hyperbolic conservation laws

In this article we develop an improved version of the classical fifth-order weighted essentially non-oscillatory finite difference scheme of [G.S. Jiang, C.W. Shu, Efficient implementation of weighted ENO schemes, J. Comput. Phys. 126 (1996) 202–228] (WENO-JS) for hyperbolic conservation laws. Through the novel use of a linear combination of the low order smoothness indicators already present in the framework of WENO-JS, a new smoothness indicator of higher order is devised and new non-oscillatory weights are built, providing a new WENO scheme (WENO-Z) with less dissipation and higher resolution than the classical WENO. This new scheme generates solutions that are sharp as the ones of the mapped WENO scheme (WENO-M) of Henrick et al. [A.K. Henrick, T.D. Aslam, J.M. Powers, Mapped weighted essentially non-oscillatory schemes: achieving optimal order near critical points, J. Comput. Phys. 207 (2005) 542–567], however with a 25% reduction in CPU costs, since no mapping is necessary. We also provide a detailed analysis of the convergence of the WENO-Z scheme at critical points of smooth solutions and show that the solution enhancements of WENO-Z and WENO-M at problems with shocks comes from their ability to assign substantially larger weights to discontinuous stencils than the WENO-JS scheme, not from their superior order of convergence at critical points. Numerical solutions of the linear advection of discontinuous functions and nonlinear hyperbolic conservation laws as the one dimensional Euler equations with Riemann initial value problems, the Mach 3 shock–density wave interaction and the blastwave problems are compared with the ones generated by the WENO-JS and WENO-M schemes. The good performance of the WENO-Z scheme is also demonstrated in the simulation of two dimensional problems as the shock–vortex interaction and a Mach 4.46 Richtmyer–Meshkov Instability (RMI) modeled via the two dimensional Euler equations.     

Balanced finite volume WENO and central WENO schemes for the shallow water and the open-channel flow equations

The goal of this work is to extend finite volume WENO and central WENO schemes to the hyperbolic balance laws with geometrical source term and spatially variable flux function. In particular, we apply proposed schemes to the shallow water and the open-channel flow equations where the source term depends on the channel geometry. For obtaining stable numerical schemes that are free of spurious oscillations, it becomes crucial to use the decomposed source term evaluation, which maintains the balancing between the flux gradient and the source term. In addition, the open-channel flow equations contain spatially variable flux function. The appropriate definitions of the terms that arise in the source term decomposition, in combination with the Roe approximate Riemann solver that includes the spatial derivative of the flux function, lead to the finite volume WENO scheme that satisfies the exact conservation property – the property of preserving the quiescent flow exactly. When the central WENO schemes are applied, additional reformulations are introduced for the transition from the staggered values to the nonstaggered ones and vice versa by using the WENO reconstruction procedure. The proposed central WENO schemes also preserve the quiescent flow, but only in prismatic channels. In various test problems the obtained balanced schemes show improvements in comparison with the standard versions of the proposed type schemes, as well as with some other first- and second-order numerical schemes.     

A direct Eulerian GRP scheme for relativistic hydrodynamics: One-dimensional case

The paper proposes a direct Eulerian generalized Riemann problem (GRP) scheme for one-dimensional relativistic hydrodynamics. It is an extension of the Eulerian GRP scheme for compressible non-relativistic hydrodynamics proposed in [M. Ben-Artzi, J.Q. Li, G. Warnecke, A direct Eulerian GRP scheme for compressible fluid flows, J. Comput. Phys. 218 (2006) 19–43]. Two main ingredients, the Riemann invariant and the Rankine–Hugoniot jump condition, are directly used to resolve the local GRP in the Eulerian formulation, and thus the crucial and delicate Lagrangian treatment in the original GRP scheme [3] can be avoided. Several numerical examples are given to demonstrate the accuracy and effectiveness of the proposed GRP scheme.     

A High-Order WENO Finite Difference Scheme for the Equations of Ideal Magnetohydrodynamics

We present a high-order accurate weighted essentially non-oscillatory (WENO) finite difference scheme for solving the equations of ideal magnetohydrodynamics (MHD). This scheme is a direct extension of a WENO scheme, which has been successfully applied to hydrodynamic problems. The WENO scheme follows the same idea of an essentially non-oscillatory (ENO) scheme with an advantage of achieving higher-order accuracy with fewer computations. Both ENO and WENO can be easily applied to two and three spatial dimensions by evaluating the fluxes dimension-by-dimension. Details of the WENO scheme as well as the construction of a suitable eigen-system, which can properly decompose various families of MHD waves and handle the degenerate situations, are presented. Numerical results are shown to perform well for the one-dimensional Brio–Wu Riemann problems, the two-dimensional Kelvin–Helmholtz instability problems, and the two-dimensional Orszag–Tang MHD vortex system. They also demonstrate the importance of maintaining the divergence free condition for the magnetic field in achieving numerical stability. The tests also show the advantages of using the higher-order scheme. The new 5th-order WENO MHD code can attain an accuracy comparable with that of the second-order schemes with many fewer grid points.     

Multi-domain Fourier-continuation/WENO hybrid solver for conservation laws

We introduce a multi-domain Fourier-continuation/WENO hybrid method (FC–WENO) that enables high-order and non-oscillatory solution of systems of nonlinear conservation laws, and which enjoys essentially dispersionless, spectral character away from discontinuities, as well as mild CFL constraints (comparable to those of finite difference methods). The hybrid scheme employs the expensive, shock-capturing WENO method in small regions containing discontinuities and the efficient FC method in the rest of the computational domain, yielding a highly effective overall scheme for applications with a mix of discontinuities and complex smooth structures. The smooth and discontinuous solution regions are distinguished using the multi-resolution procedure of Harten [J. Comput. Phys. 115 (1994) 319–338]. We consider WENO schemes of formal orders five and nine and a FC method of order five. The accuracy, stability and efficiency of the new hybrid method for conservation laws is investigated for problems with both smooth and non-smooth solutions. In the latter case, we solve the Euler equations for gas dynamics for the standard test case of a Mach three shock wave interacting with an entropy wave, as well as a shock wave (with Mach 1.25, three or six) interacting with a very small entropy wave and evaluate the efficiency of the hybrid FC–WENO method as compared to a purely WENO-based approach as well as alternative hybrid based techniques. We demonstrate considerable computational advantages of the new FC-based method, suggesting a potential of an order of magnitude acceleration over alternatives when extended to fully three-dimensional problems.     

A systematic methodology for constructing high-order energy stable WENO schemes

A third-order Energy Stable Weighted Essentially Non-Oscillatory (ESWENO) finite difference scheme developed by the authors of this paper [N.K. Yamaleev, M.H. Carpenter, Third-order energy stable WENO scheme, J. Comput. Phys. 228 (2009) 3025–3047] was proven to be stable in the energy norm for both continuous and discontinuous solutions of systems of linear hyperbolic equations. Herein, a systematic approach is presented that enables “energy stable” modifications for existing WENO schemes of any order. The technique is demonstrated by developing a one-parameter family of fifth-order upwind-biased ESWENO schemes including one sixth-order central scheme; ESWENO schemes up to eighth order are presented in the Appendix. We also develop new weight functions and derive constraints on their parameters, which provide consistency, much faster convergence of the high-order ESWENO schemes to their underlying linear schemes for smooth solutions with arbitrary number of vanishing derivatives, and better resolution near strong discontinuities than the conventional counterparts.     

An improved gas-kinetic BGK finite-volume method for three-dimensional transonic flow

During the past decade gas-kinetic methods based on the BGK simplification of the Boltzmann equation have been employed to compute fluid flow in a finite-difference or finite-volume context. Among the most successful formulations is the finite-volume scheme proposed by Xu [K. Xu, A gas-kinetic BGK scheme for the Navier–Stokes equations and its connection with artificial dissipation and Godunov method, J. Comput. Phys. 171 (48) (2001) 289–335]. In this paper we build on this theoretical framework mainly with the aim to improve the efficiency and convergence of the scheme, and extend the range of application to three-dimensional complex geometries using general unstructured meshes. To that end we propose a modified BGK finite-volume scheme, which significantly reduces the computational cost, and improves the behavior on stretched unstructured meshes. Furthermore, a modified data reconstruction procedure is presented to remove the known problem that the Chapman–Enskog expansion of the BGK equation fixes the Prandtl number at unity. The new Prandtl number correction operates at the level of the partial differential equations and is also significantly cheaper for general formulations than previously published methods. We address the issue of convergence acceleration by applying multigrid techniques to the kinetic discretization. The proposed modifications and convergence acceleration help make large-scale computations feasible at a cost competitive with conventional discretization techniques, while still exploiting the advantages of the gas-kinetic discretization, such as computing full viscous fluxes for finite volume schemes on a simple two-point stencil.     

Construction of low dissipative high-order well-balanced filter schemes for non-equilibrium flows

The goal of this paper is to generalize the well-balanced approach for non-equilibrium flow studied by Wang et al. (2009) [29] to a class of low dissipative high-order shock-capturing filter schemes and to explore more advantages of well-balanced schemes in reacting flows. More general 1D and 2D reacting flow models and new examples of shock turbulence interactions are provided to demonstrate the advantage of well-balanced schemes. The class of filter schemes developed by Yee et al. (1999) [33], Sjögreen and Yee (2004) [27] and Yee and Sjögreen (2007) [38] consist of two steps, a full time step of spatially high-order non-dissipative base scheme and an adaptive non-linear filter containing shock-capturing dissipation. A good property of the filter scheme is that the base scheme and the filter are stand-alone modules in designing. Therefore, the idea of designing a well-balanced filter scheme is straightforward, i.e. choosing a well-balanced base scheme with a well-balanced filter (both with high-order accuracy). A typical class of these schemes shown in this paper is the high-order central difference schemes/predictor–corrector (PC) schemes with a high-order well-balanced WENO filter. The new filter scheme with the well-balanced property will gather the features of both filter methods and well-balanced properties: it can preserve certain steady-state solutions exactly; it is able to capture small perturbations, e.g. turbulence fluctuations; and it adaptively controls numerical dissipation. Thus it shows high accuracy, efficiency and stability in shock/turbulence interactions. Numerical examples containing 1D and 2D smooth problems, 1D stationary contact discontinuity problem and 1D turbulence/shock interactions are included to verify the improved accuracy, in addition to the well-balanced behavior.     

基于WENO重构的高阶移动网格动理学格式

提出一种基于WENO重构的高阶(至少三阶)移动网格动理学格式.利用流体力学方程的积分形式得到移动网格上离散格式,再利用自适应移动网格方法移动网格,进而得到网格速度,利用WENO重构得到高阶插值多项式,最后使用时间方向上精确的动理学数值方法构造数值通量,得到移动网格单元上新的物理量.数值实验表明这种格式同时具有高精度、高分辨率的特点.     

Further improvement and analysis of CCD scheme: Dissipation discretization and de-aliasing properties

In this paper, we further analyze a combined compact difference (CCD) scheme proposed recently [T.K. Sengupta, V. Lakshmanan, V.V.S.N. Vijay, A new combined stable and dispersion relation preserving compact scheme for non-periodic problems, J. Comput. Phys. 228 (8) (2009) 3048–3071] for its dissipation discretization properties to show that its superiority also helps in controlling aliasing error for a benchmark internal flow. However, application of the same CCD method to study the receptivity of a boundary layer experiencing adverse pressure gradient is not successful. This is traced to the nature of the equilibrium flow where the better dissipation property is not helpful in the inviscid part of the flow, where the aliasing problems continue to persist. A further modification is proposed to the CCD method here to solve complex physical problems requiring information on higher order disturbance quantities – as in problems of flow receptivity and instability.     

迎风紧致格式求解Hamilton-Jacobi方程

基于Hamilton-Jacobi(H-J)方程和双曲型守恒律之间的关系,将三阶和五阶迎风紧致格式推广应用于求解H-J方程,建立了高精度的H-J方程求解方法.给出了一维和二维典型数值算例的计算结果,其中包括一个平面激波作用下的Richtmyer Meshkov界面不稳定性问题.数值试验表明,在解的光滑区域该方法具有高精度,而在导数不连续的不光滑区域也获得了比较好的分辨效果.相比于同阶精度的WENO格式,本方法具有更小的数值耗散,从而有利于多尺度复杂流动的模拟中H-J方程的求解.     

New High-Resolution Semi-discrete Central Schemes for Hamilton–Jacobi Equations

We introduce a new high-resolution central scheme for multidimensional Hamilton–Jacobi equations. The scheme retains the simplicity of the non-oscillatory central schemes developed by C.-T. Lin and E. Tadmor (in press, SIAM J. Sci. Comput.), yet it enjoys a smaller amount of numerical viscosity, independent of 1/Δt. By letting Δt↓0 we obtain a new second-order central scheme in the particularly simple semi-discrete form, along the lines of the new semi-discrete central schemes recently introduced by the authors in the context of hyperbolic conservation laws. Fully discrete versions are obtained with appropriate Runge–Kutta solvers. The smaller amount of dissipation enables efficient integration of convection-diffusion equations, where the accumulated error is independent of a small time step dictated by the CFL limitation. The scheme is non-oscillatory thanks to the use of nonlinear limiters. Here we advocate the use of such limiters on second discrete derivatives, which is shown to yield an improved high resolution when compared to the usual limitation of first derivatives. Numerical experiments demonstrate the remarkable resolution obtained by the proposed new central scheme.     

A direct Eulerian GRP scheme for relativistic hydrodynamics: Two-dimensional case

This paper develops a direct Eulerian generalized Riemann problem (GRP) scheme for two-dimensional (2D) relativistic hydrodynamics (RHD). It is an extension of the GRP scheme for one-dimensional (1D) RHDs [Z.C. Yang, P. He, H.Z. Tang, J. Comput. Phys. 230 (2011) 7964–7987] and the GRP scheme for the non-relativistic hydrodynamics [M. Ben-Artzi, J.Q. Li, G. Warnecke, J. Comput. Phys. 218 (2006) 19–43]. In order to derive the direct Eulerian GRP scheme, the (local) GRP of the split 2D RHD equations in the Eulerian formulation has to be directly resolved by using corresponding Riemann invariants and Rankine–Hugoniot jump conditions so that the crucial and delicate Lagrangian treatment in the original GRP scheme [M. Ben-Artzi, J. Falcovitz, J. Comput. Phys. 55 (1984) 1–32] may be avoided. An important difference of resolving the GRP of the split 2D RHD equations from the GRP of the 1D RHD equations or the non-relativistic hydrodynamical equations is coming from the fact that the flow regions across the shock or rarefaction wave in the GRP of the split 2D RHD equations are nonlinearly coupled through the Lorentz factor which is also built in terms of the tangential velocities. It is a purely multi-dimensional relativistic feature. Several numerical examples are given to demonstrate the accuracy and effectiveness of the proposed 2D GRP scheme.