Preventing numerical oscillations in the flux‐split based finite difference method for compressible flows with discontinuities,II

Problems in the characteristic‐wise flux‐split based finite difference method when compressible flows with contact discontinuities or material interfaces are computed were presented and analyzed. The current analysis showed the following: (i) Even with the local characteristic decomposition technique, numerical errors could be caused by point‐wise flux vector splitting (FVS) methods, such as the Steger–Warming FVS or the van Leer FVS. Therefore, the Lax–Friedrichs type FVS method is required. (ii) If the isobars of a material are vertical lines, the combination of using the local characteristic decomposition and the global Lax–Friedrichs FVS can avoid velocity and pressure oscillations of contact discontinuities in this material for weighted essentially non‐oscillatory (WENO) schemes. (iii) For problems with material interfaces, the quasi‐conservative approach can be realized using characteristic‐wise flux‐split based finite difference WENO schemes if nonlinear WENO schemes in genuinely nonlinear characteristic fields can be guaranteed to be the same and the decomposition equation representing material interfaces is discretized properly. Copyright © 2015 John Wiley & Sons, Ltd.     

基于WENO重构的真正多维黎曼求解器

具有良好守恒性与网格适应性的有限体积格式在流体力学的数值计算中占有重要地位。其中，求解数值流通量是实施有限体积法的关键步骤。一维情形下，通过求解局部黎曼问题来获得数值流通量的相关理论已经比较成熟。但是在计算多维问题时，传统的维度分裂方法仅考虑沿界面法向传播的信息，这不仅影响格式的精度，还可能会造成数值不稳定性从而诱发非物理现象。本文基于对流-压力通量分裂方法来构造真正多维的黎曼求解器，通过求解网格顶点处的多维黎曼问题来实现格式的多维特性。采用五阶WENO重构方法来获得空间的高阶精度，时间离散采用三阶TVD龙格-库塔格式。一系列数值实验的结果表明，真正多维的黎曼求解器不仅具有更高的分辨率还能有效克服多维强激波模拟中的数值不稳定性。     

Realization of contact resolving approximate Riemann solvers for strong shock and expansion flows

The Harten–Lax–van Leer contact (HLLC) and Roe schemes are good approximate Riemann solvers that have the ability to resolve shock, contact, and rarefaction waves. However, they can produce spurious solutions, called shock instabilities, in the vicinity of strong shock. In strong expansion flows, the Roe scheme can admit nonphysical solutions such as expansion shock, and it sometimes fails. We carefully examined both schemes and propose simple methods to prevent such problems. High‐order accuracy is achieved using the weighted average flux (WAF) and MUSCL‐Hancock schemes. Using the WAF scheme, the HLLC and Roe schemes can be expressed in similar form. The HLLC and Roe schemes are tested against Quirk's test problems, and shock instability appears in both schemes. To remedy shock instability, we propose a control method of flux difference across the contact and shear waves. To catch shock waves, an appropriate pressure sensing function is defined. Using the proposed method, shock instabilities are successfully controlled. For the Roe scheme, a modified Harten–Hyman entropy fix method using Harten–Lax–van Leer‐type switching is suggested. A suitable criterion for switching is established, and the modified Roe scheme works successfully with the suggested method. Copyright © 2009 John Wiley & Sons, Ltd.     

High resolution flux‐difference‐splitting scheme on adaptive grid for open‐channel flows

The effectiveness and usefulness of further enhancing the shock resolution of a second‐order accurate scheme for open‐channel flows by using an adaptive grid is investigated. The flux‐difference‐splitting (FDS) scheme based on the Lax–Wendroff numerical flux is implemented on a fixed as well as on a self‐adjusting grid for this purpose. The grid‐adjusting procedure, developed by Harten and Hyman, adjusts the grid by averaging the local characteristic velocities with respect to the signal amplitude in such a way that a shock always lies on a mesh point. This enables a scheme capable of perfectly resolving a stationary shock to capture a shock that moves from mesh point to mesh point. The Roe's approximate Jacobian is used for conservation and consistency, while theoretically sound treatment for satisfying entropy inequality conditions ensures physically realistic solutions. Details about inclusion of source terms, often left out of analyses for the homogeneous part of governing equations, are also explained. The numerical results for some exacting problems are compared with analytical as well as experimental results for examining improvements in resolution of discontinuities by the adaptive grid. Copyright © 2001 John Wiley & Sons, Ltd.     

一种激波稳定的对流-压力通量分裂格式

针对欧拉方程三种流行的对流-压力通量分裂方法(Liou-Steffen,Zha-Bilgen和Toro-Vázquez)进行特征分析,进而提出一种新的对流-压力通量分裂格式。采用Zha-Bilgen分裂方法将欧拉方程的通量分裂成对流项和压力项两部分,使用TV格式来计算这两部分的数值通量。利用压力比构造激波探测函数,并且在强激波附近的亚声速区域增加TV格式的剪切粘性来克服数值模拟中的激波不稳定性。数值算例的计算结果表明,新的对流-压力通量分裂格式不仅保留了原始TV格式精确分辨接触间断的优点,而且具有更好的鲁棒性,在数值模拟多维强激波问题时不会出现不稳定现象。因此,该格式是一种精确并且具有强鲁棒性的数值方法,可以广泛地应用于可压缩流体的数值计算中。     

Numerical solutions of Euler equations by using a new flux vector splitting scheme

A new flux vector splitting scheme has been suggested in this paper. This scheme uses the velocity component normal to the volume interface as the characteristic speed and yields the vanishing individual mass flux at the stagnation. The numerical dissipation for the mass and momentum equations also vanishes with the Mach number approaching zero. One of the diffusive terms of the energy equation does not vanish. But the low numerical diffusion for viscous flows may be ensured by using higher-order differencing. The scheme is very simple and easy to be implemented. The scheme has been applied to solve the one dimensional (1D) and multidimensional Euler equations. The solutions are monotone and the normal shock wave profiles are crisp. For a 1D shock tube problem with the shock and the contact discontinuities, the present scheme and Roe scheme give very similar results, which are the best compared with those from Van Leer scheme and Liou–Steffen's advection upstream splitting method (AUSM) scheme. For the multidimensional transonic flows, the sharp monotone normal shock wave profiles with mostly one transition zone are obtained. The results are compared with those from Van Leer scheme, AUSM and also with the experiment.     

Healing of shock instability for Roe's flux‐difference splitting scheme on triangular meshes

Healing of nonphysical flow solutions and shock instability from the use of Roe's flux‐difference splitting scheme is presented. The proposed method heals nonphysical flow solutions such as the carbuncle phenomenon, the shock instability from the odd–even decoupling problem, and the expansion shock generated from the violated entropy condition. The performance and efficiency of the proposed method are evaluated by solving several benchmark and complex high‐speed compressible flow problems. Copyright © 2008 John Wiley & Sons, Ltd.     

Flux‐difference splitting‐based upwind compact schemes for the incompressible Navier–Stokes equations

Third‐order and fifth‐order upwind compact finite difference schemes based on flux‐difference splitting are proposed for solving the incompressible Navier–Stokes equations in conjunction with the artificial compressibility (AC) method. Since the governing equations in the AC method are hyperbolic, flux‐difference splitting (FDS) originally developed for the compressible Euler equations can be used. In the present upwind compact schemes, the split derivatives for the convective terms at grid points are linked to the differences of split fluxes between neighboring grid points, and these differences are computed by using FDS. The viscous terms are approximated with a sixth‐order central compact scheme. Comparisons with 2D benchmark solutions demonstrate that the present compact schemes are simple, efficient, and high‐order accurate. Copyright © 2008 John Wiley & Sons, Ltd.     

一种新型的高分辨率通量差分裂格式

传统的Roe格式不满足熵条件并且在计算激波问题时会遭遇不同形式的不稳定现象,如慢行激波的波后振荡和红玉(carbuncle)现象.基于Zha-Bilgen对流-压力通量分裂方法,构造一种新型的通量差分裂格式.利用约旦标准型理论,通过添加广义特征向量构造通量差分裂方法来计算对流子系统.压力子系统具有一组完备的线性无关特征向量,因此可以构造传统的通量差分裂格式进行计算.为了提高接触间断的分辨率,利用界面变差下降(BVD)算法来重构对流通量耗散项中的密度差.激波稳定性分析表明,新格式可以有效地衰减数值误差,从而抑制不稳定现象的发生.一系列数值实验证明了本文构造的新型通量差分裂格式比Roe格式具有更高的分辨率和更好的鲁棒性.     

Modified Multidimensional Dissipation Scheme on Unstructured Meshes for High-speed Compressible Flow Analysis

A Roe's flux-difference splitting scheme, combining with the entropy fix method according to Van Leer et al., and the H-correction entropy fix method by Pandolfi and D'Ambrosio, is proposed. The presented scheme eliminates unphysical flow behaviors such as a spurious bump of the carbuncle phenomenon that occurs on the bow shock from flow over a blunt body, and the expansion shock generated from flow over a forward facing step. The proposed scheme is further extended to obtain high-order spatial and temporal solution accuracy. The scheme is, in addition, combined with an adaptive meshing technique that generates unstructured triangular meshes to resemble the flow phenomena for reducing computational effort. The entire procedure is evaluated by solving several benchmarks as well as complex steady-state and transient high-speed compressible flow problems.     

求解可压缩流的二维通量分裂格式

传统的一维通量分裂格式在计算界面数值通量时,只考虑网格界面法向的波系。采用传统的TV格式分别求解对流通量和压力通量。通过求解考虑了横向波系影响的角点数值通量来构造一种真正二维的TV通量分裂格式。在计算一维数值算例时,该格式与传统的TV格式具有相同的数值通量计算公式,因此其保留了传统的TV格式精确捕捉接触间断和膨胀激波的优点。在计算二维算例时,该格式比传统的TV格式具有更高的分辨率;在计算二维强激波问题时,消除了传统TV格式的非物理现象,表现出更好的鲁棒性;此外,该格式大大提高了稳定性CFL数,从而具有更高的计算效率。因此,本文方法是一种精确、高效并且具有强鲁棒性的数值方法,在可压缩流的数值模拟中具有广阔的应用前景。     

A study on the behaviour of high-order flux reconstruction method with different low-dissipation numerical fluxes for large eddy simulation

This work focuses on the numerical dissipation features of high-order flux reconstruction (FR) method combined with different numerical fluxes in turbulence flows. The famous Roe and AUSM+ numerical fluxes together with their corresponding low-dissipation enhanced versions (LMRoe, SLAU2) and higher resolution variants (HR-LMRoe, HR-SLAU2) are incorporated into FR framework, and the dissipation interplay of these combinations is investigated in implicit large eddy simulation. The numerical dissipation stemming from these convective numerical fluxes is quantified by simulating the inviscid Gresho vortex, the transitional Taylor–Green vortex and the homogenous decaying isotropic turbulence. The results suggest that low-dissipation enhanced versions are preferential both in high-order and low-order cases to their original forms, while the use of HR-SLAU2 has marginal improvements and the HR-LMRoe leads to degenerated solution with high-order. In high-order the effects of numerical fluxes are reduced, and their viscosity may not be dissipative enough to provide physically consistent turbulence when under-resolved.     

An efficient upwind/relaxation algorithm for the Euler and Navier–Stokes equations

An efficient Euler and full Navier–Stokes solver based on a flux splitting scheme is presented. The original Van Leer flux vector splitting form is extended to arbitrary body-fitted co-ordinates in the physical domain so that it can be used with a finite volume scheme. The block matrix is inverted by Gauss–Seidel iteration. It is verified that the often used reflection boundary condition will produce incorrect flux crossing the wall and cause too large numerical dissipation if flux vector splitting is used. To remove such errors, an appropriate treatment of wall boundary conditions is suggested. Inviscid and viscous steady transonic internal flows are analysed, including the case of shock-induced boundary layer separation.     

A robust low‐dissipation AUSM‐family scheme for numerical shock stability on unstructured grids

The simple low‐dissipation advection upwind splitting method (SLAU) scheme is a parameter‐free, low‐dissipation upwind scheme that has been applied in a wide range of aerodynamic numerical simulations. In spite of its successful applications, the SLAU scheme could be showing shock instabilities on unstructured grids, as many other contact resolved upwind schemes. Therefore, a hybrid upwind flux scheme is devised for improving the shock stability of SLAU scheme, without compromising on accuracy and low Mach number performance. Numerical flux function of the hybrid scheme is written in a general form, in which only the scalar dissipation term is different from that of the SLAU scheme. The hybrid dissipation term is defined by using a differentiable multidimensional‐shock‐detection pressure weight function, and the dissipation term of SLAU scheme is combined with that of the Van Leer scheme. Furthermore, the hybrid dissipation term is only applied for the solution of momentum fluxes in numerical flux function. Based on the numerical test results, the hybrid scheme is deemed to be a successful improvement on the shock stability of SLAU scheme, without compromising on the efficiency and accuracy. Copyright © 2016 John Wiley & Sons, Ltd.     

N-S方程在非结构网格下的求解

在Ｒｏｅ的矢通量差分分裂的基础上，吸收了ＮＮＤ格式的优点，提出了一种非结构网格下求解Ｅｕｌｅｒ方程和Ｎ－Ｓ方程的高分辨率高精度迎风格式．这种格式具有捕捉强激波和滑移线的良好性能．在时间方向上采用了显式和隐式两种解法．文中还给出了自适应技术．最后，成功地完成了ＧＡＭＭ超音速前台阶绕流、二维平板无粘激波反射、三维Ｈｏｂｓｏｎ叶栅流动、ＶＫＩ叶栅流动、Ｃ３Ｘ叶栅流动的数值模拟，得到了满意的结果     

Developing a hybrid flux function suitable for hypersonic flow simulation with high‐order methods

In this paper, we develop a new hybrid Euler flux function based on Roe's flux difference scheme, which is free from shock instability and still preserves the accuracy and efficiency of Roe's flux scheme. For computational cost, only 5% extra CPU time is required compared with Roe's FDS. In hypersonic flow simulation with high‐order methods, the hybrid flux function would automatically switch to the Rusanov flux function near shock waves to improve the robustness, and in smooth regions, Roe's FDS would be recovered so that the advantages of high‐order methods can be maintained. Multidimensional dissipation is introduced to eliminate the adverse effects caused by flux function switching and further enhance the robustness of shock‐capturing, especially when the shock waves are not aligned with grids. A series of tests shows that this new hybrid flux function with a high‐order weighted compact nonlinear scheme is not only robust for shock‐capturing but also accurate for hypersonic heat transfer prediction. Copyright © 2015 John Wiley & Sons, Ltd.     

Simulation of free surface flows using the flux‐difference splitting scheme on the hybrid Cartesian/immersed boundary method

In this study, a method is developed to simulate the interaction between free surface flows and moving or deforming boundaries using the flux‐difference splitting scheme on the hybrid Cartesian/immersed boundary method. At each physical time step, the boundary is defined by an unstructured triangular surface grid. Immersed boundary (IB) nodes are distributed inside an instantaneous fluid domain based on edges crossing the boundary. At an IB node, dependent variables are reconstructed along the local normal line to the boundary. Inviscid fluxes are computed using Roe's flux‐difference splitting scheme for immiscible and incompressible fluids. The free surface is considered as a contact discontinuity in the density field. The motion of free surface is captured without any additional treatment along the fluid interface. The developed code is validated by comparisons with other experimental and computational results for a piston‐type wave maker, impulsive motion of a submerged circular cylinder, flow around a submerged hydrofoil, and Rayleigh–Taylor instability. The developed code is applied to simulate wave generation due to a continuously deforming bed beneath the free surface. The violent motion of a free surface caused by sloshing in a spherical tank is simulated. In this case, the free surface undergoes breakup and reconnection. Copyright © 2011 John Wiley & Sons, Ltd.     

Calculation of Reynolds stresses in turbulent supersonic flows

The baseline numerical procedure of interest in this study combines flux vector splitting, flux difference splitting and an explicit treatment of the diffusion terms of the flow equations. The viscous terms are treated explicitly to preserve the wave propagation properties of the Euler fluxes and permit splitting. The experience with this scheme has been limited to laminar or, at best, ‘eddy viscosity’ flows. In this paper the applicability of the scheme is extended to include the calculation of turbulent Reynolds stresses in supersonic flows. The schemes and our implementation are discussed. Both laminar and turbulence subsets of the Reynolds/Favre-averaged equations are tested, with a discussion of relative performance. The test problem for turbulence consists of a zero-pressure-gradient supersonic boundary layer as well as a supersonic boundary layer experiencing the combined effects of adverse pressure gradient, bulk compression and a concave streamline curvature. Excellent agreement with experimental measurements is observed for most of the quantities compared, which suggests that the numerical procedures presented in this paper are potentially very useful.     

On the spectrum of the Steger‐Warming flux‐vector splitting scheme

The flux‐vector splitting scheme of Steger and Warming is a popular approach for the Euler equations. In this work, we consider the spectrum of the scheme and show for 1≤γ≤5/3, where γ is the ideal gas constant, that the eigenvalues are strictly real and of an appropriate sign.     

A Godunov‐type fractional semi‐implicit method based on staggered grid for dam‐break modeling

A semi‐implicit finite volume model based upon staggered grid is presented for solving shallow water equation. The model employs a time‐splitting scheme that uses a predictor–corrector method for the advection term. The fluxes are calculated based on a Riemann solver in the prediction step and a downwind scheme in the correction step. A simple TVD scheme is employed for shock capturing purposes in which the Minmond limiter is used for flux functions. As a consequence of using staggered grid, an ADI method is adopted for solving the discretized equations for 2‐D problems. Several 1‐D and 2‐D flows have been modeled with satisfactory results when compared with analytical and experimental test cases. The model is also capable of simulating supercritical as well as subcritical flow. Copyright © 2010 John Wiley & Sons, Ltd.