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

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

满足几何守恒律的WENO格式及其应用

对几何守恒律的来源进行了分析,发展了一种满足几何守恒律的WENO格式,并应用于翼型层流分离现象的数值模拟中。为消除网格质量影响,采用守恒型方法计算网格导数,并将标准的WENO格式分解为中心差分部分和数值耗散部分。算例计算结果表明,几何守恒律对高精度有限差分WENO格式计算结果具有重要影响,本文方法能够消除网格导数计算误差,保证来流保持性。将本文方法应用于SD7003翼型层流分离现象的数值模拟中,计算结果与文献中计算及试验数据吻合较好,同时能够精细捕捉小尺度流场结构,准确模拟翼型层流分离现象中的复杂流动过程。     

翼型、机翼非定常运动模式下动态数值模拟

基于RANS方程,通过刚性动网格技术实现对翼型和机翼典型运动模式的描述,采用双时间推进方法和Roe空间离散格式对流场求解,构建了一个非定常气动计算平台;以NACA0012翼型为算倒进行了动态数值模拟可信度验证。数值模拟结果与试验数据吻合较好,升力和俯仰力矩的最大计算误差分别为3％和10％,表明了该平台的可靠性。另外,还数值模拟了M6机翼的动态非定常流场,并分析了两种湍流模型对非定常流场激波的捕捉能力。结果表明非定常流动中S-A湍流模型对激波的捕捉较B-L模型更敏感。文中开发的非定常计算平台对进一步解决三维复杂流场的流动问题有很高的工程应用价值。     

近似黎曼解对高超声速气动热计算的影响研究

高超声速流场计算一般采用TVD型格式,这些格式中,大多采用了不同形式的近似黎曼解. 通过分析和数值验证,论述了激波捕捉格式中近似黎曼解的耗散性质,说明其对高超声速热流计算的影响. 数值实验证明,采用低耗散格式可大大提高热流计算精度,降低热流计算对网格的依赖程度,从而获得精确的热流数值解.      

激波与堆积粉尘相互作用的数值模拟

基于双流体模型和测定的堆积粉尘的本构方程 ,利用AUSM+ 格式 ,对激波与堆积粉尘的相互作用进行了数值模拟。计算所反映的流场结构与实验图像一致。此外还对激波强度 ,颗粒材料密度等对流场的影响进行了讨论。     

一种鲁棒的HLLC格式及其稳定性分析

精确捕捉接触波和剪切波的Godunov型数值方法，如流行的HLLC格式，在模拟高超声速流动问题时会出现激波异常现象。对HLLC格式进行稳定性分析发现，流体主流方向的扰动都能有效衰减，但是横向的密度与剪切速度的扰动不会衰减。具有特殊对称性的二维Sedov爆轰波问题证明了横向通量和不稳定现象之间的密切联系。利用压力比和马赫数来探测数值激波层亚声速区的横向网格界面，并且在该界面的数值通量上增加熵波粘性和剪切波粘性来构造一种激波稳定的HLLC格式。分析表明，在熵波粘性和剪切波粘性的作用下，横向的所有扰动都会衰减。一系列数值测试证明了新格式不仅可以成功地抑制各类激波异常现象，还保留了原HLLC格式低耗散性的优点。     

HLLE++格式在高马赫数空腔流动模拟中的应用

空腔流动存在剪切层运动、涡脱落与破裂,以及激波与激波、激波与剪切层、激波与膨胀波和激波/涡/剪切层相互干扰等现象,流动非常复杂,特别是高马赫数(M>2)时,剪切层和激波更强,激波与激波干扰更严重,对数值格式的要求更高,既需要格式耗散小,对分离涡等有很高的模拟精度,又需要格式在激波附近具有较大的耗散,可以很好地捕捉激波,防止非物理解的出现。Roe和HLLC等近似Riemann解格式在高马赫数强激波处可能会出现红玉现象,而HLLE++格式大大改善了这种缺陷,在捕捉高超声速激波时避免了红玉现象的发生,同时还保持在光滑区域的低数值耗散特性。本文在结构网格下HLLE++格式的基础上,通过改进激波探测的求解,建立了基于非结构混合网格的HLLE++计算方法,通过无粘斜坡算例,验证了HLLE++格式模拟高马赫数流动的能力,并应用于高马赫数空腔流动的数值模拟,开展了网格和湍流模型影响研究,验证了方法模拟高马赫数空腔流动的可靠性和有效性。     

尺度无关高阶WCNS格式

高速流场的数值模拟中, 既要保证对小尺度结构的高保真分辨, 又要实现对激波稳定、无振荡地捕捉.当前工程中广泛应用的高精度数值格式虽然都能一定程度地满足上述两种要求, 但仍与理想目标存在较大差距.例如, 模拟雷诺应力模型等小尺度问题时, 高精度格式在间断解附近易产生数值振荡.基于高精度格式所存在的上述问题, 本文引入去尺度函数, 探索了一种更加简单稳定的非线性权重构造方法, 并将其应用于7阶精度加权紧致非线性格式WCNS, 提出了一种尺度无关的7阶WCNS格式.该格式的性能与灵敏度参数和尺度因子的选择无关, 并且在小尺度下仍可以有效捕捉流场激波.同时, 该格式在间断处具有基本无振荡性质, 且在任意尺度函数下保持尺度无关, 并且在极值点处也能保持最优精度.本文还推导了7阶D权函数的形式.最后, 在一维线性对流方程中验证了新格式在流场光滑区能够达到设计精度, 并通过一系列数值实验证明了尺度无关的7阶WCNS格式在激波捕捉能力上具有良好表现, 为WCNS格式改进和解决可压缩湍流等非线性问题提供了一种新途径.      

数值研究激波与旋涡的相互作用

从非定常形式的Euler方程出发,数值模拟了运动激波与旋涡相互作用的非定常流动过程。为保证激波具有较高的分辨率,采用对称型TVD格式进行了数值计算。结果表明。这样可以有效地模拟流场中一些复杂的流动现象,如激波变形、激波分叉和三波点的形成,以及旋涡结构的变化过程等,并与已有的实验流动显示相符良好。同时,也是对TVD格式求解这类问题的一次初步尝试。     

环形激波绕射、反射和聚焦流场的间断有限元模拟

采用间断有限元方法对环形激波在圆柱形激波管内绕射、反射和聚焦流场进行了数值模拟。将二维守恒方程的间断有限元方法发展到轴对称Euler方程,并对环形激波绕后台阶流动进行了数值计算。计算结果表明,采用间断有限元方法能够有效地捕捉运动激波在圆柱形激波管内传播的复杂流场结构;在聚焦点附近,数值解具有较大的梯度变化,表明该方法对间断解具有较强的捕捉能力,在聚焦点附近不会产生振荡或抹平间断现象。     

Solution to Euler equations by high‐resolution upwind compact scheme based on flux splitting

A high‐resolution upwind compact method based on flux splitting is developed for solving the compressive Euler equations. The convective flux terms are discretized by using the modified advection upstream splitting method (AUSM). The developed scheme is used to compute the one‐dimensional Burgers equation and four different example problems of supersonic compressible flows, respectively. The results show that the high‐resolution upwind compact scheme based on modified AUSM+ flux splitting can capture shock wave and other discontinuities, obtain higher resolution and restrain numerical oscillation. Copyright © 2007 John Wiley & Sons, Ltd.     

Study on the numerical schemes for hypersonic flow simulation

Hypersonic flow is full of complex physical and chemical processes, hence its investigation needs careful analysis of existing schemes and choosing a suitable scheme or designing a brand new scheme. The present study deals with two numerical schemes Harten, Lax, and van Leer with Contact (HLLC) and advection upstream splitting method (AUSM) to effectively simulate hypersonic flow fields, and accurately predict shock waves with minimal diffusion. In present computations, hypersonic flows have been modeled as a system of hyperbolic equations with one additional equation for non-equilibrium energy and relaxing source terms. Real gas effects, which appear typically in hypersonic flows, have been simulated through energy relaxation method. HLLC and AUSM methods are modified to incorporate the conservation laws for non-equilibrium energy. Numerical implementation have shown that non-equilibrium energy convect with mass, and hence has no bearing on the basic numerical scheme. The numerical simulation carried out shows good comparison with experimental data available in literature. Both numerical schemes have shown identical results at equilibrium. Present study has demonstrated that real gas effects in hypersonic flows can be modeled through energy relaxation method along with either AUSM or HLLC numerical scheme.     

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.     

Evaluation of Euler fluxes by a high-order CFD scheme: shock instability

The construction of Euler fluxes is an important step in shock-capturing/upwind schemes. It is well known that unsuitable fluxes are responsible for many shock anomalies, such as the carbuncle phenomenon. Three kinds of flux vector splittings (FVSs) as well as three kinds of flux difference splittings (FDSs) are evaluated for the shock instability by a fifth-order weighted compact nonlinear scheme. The three FVSs are Steger–Warming splitting, van Leer splitting and kinetic flux vector splitting (KFVS). The three FDSs are Roe's splitting, advection upstream splitting method (AUSM) type splitting and Harten–Lax–van Leer (HLL) type splitting. Numerical results indicate that FVSs and high dissipative FDSs undergo a relative lower risk on the shock instability than that of low dissipative FDSs. However, none of the fluxes evaluated in the present study can entirely avoid the shock instability. Generally, the shock instability may be caused by any of the following factors: low dissipation, high Mach number, unsuitable grid distribution, large grid aspect ratio, and the relative shock-internal flow state (or position) between upstream and downstream shock waves. It comes out that the most important factor is the relative shock-internal state. If the shock-internal state is closer to the downstream state, the computation is at higher susceptibility to the shock instability. Wall-normal grid distribution has a greater influence on the shock instability than wall-azimuthal grid distribution because wall-normal grids directly impact on the shock-internal position. High shock intensity poses a high risk on the shock instability, but its influence is not as much as the shock-internal state. Large grid aspect ratio is also a source of the shock instability. Some results of a second-order scheme and a first-order scheme are also given. The comparison between the high-order scheme and the two low-order schemes indicates that high-order schemes are at a higher risk of the shock instability. Adding an entropy fix is very helpful in suppressing the shock instability for the two low-order schemes. When the high-order scheme is used, the entropy fix still works well for Roe's flux, but its effect on the Steger–Warming flux is trivial and not much clear.     

A simple hybrid finite volume solver for compressible turbulence

A simple, explicit, hybrid finite volume method for simulating compressible turbulence is developed by combining a fourth‐order central scheme and a shock‐capturing simple low‐dissipation advection upstream splitting method. The total flux on each of the cell faces is computed as a weighted average of central/nondissipative and upwind/dissipative fluxes. The weights are determined using an unphysical oscillation sensor in addition to a more traditional discontinuity sensor used in earlier studies. Shocks are well captured, but overshoots in density are predicted around contact discontinuities that are normal to the flow. The use of the latter sensor effectively prevents these overshoots from generating spurious oscillations that travel away from the contact lines. The efficacy of the proposed method for direct or large‐eddy simulations of supersonic turbulence is established using several canonical test problems. Copyright © 2015 John Wiley & Sons, Ltd.     

Delayed mesh relaxation for multi-material ALE formulation

Multi-Material Arbitrary Lagrangian–Eulerian (ALE) finite element methods can solve large deformations in fast dynamic problems like explosions because the mesh motion can be independent of the material motion. However materials must flow between elements and this advection involves numerical dissipations. The rezoning mesh method presented in this paper was designed to reduce these numerical errors for shock wave propagation. The mesh moves to refine the elements near the shock front. This refinement limits the advection fluxes and so the numerical diffusion. This technique is applied to the numerical simulations of airblast problems for which a parameter controlling the mesh refinement is studied.     

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

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

Shock‐capturing with natural high‐frequency oscillations

This paper explores the potential of a newly developed conjugate filter oscillation reduction (CFOR) scheme for shock‐capturing under the influence of natural high‐frequency oscillations. The conjugate low‐ and high‐pass filters are constructed based on the principle of the discrete singular convolution (DSC), a local spectral method. The accuracy and resolution of the DSC basic algorithm are accessed with a one‐dimensional advection equation. Two Euler systems, the advection of an isotropic vortex flow and the interaction of shock–entropy wave are utilized to demonstrate the utility of the CFOR scheme. Computational accuracy and order of approximation are examined and compared with the literature. Some of the best numerical results are obtained for the shock–entropy wave interaction. Numerical experiments indicate that the CFOR scheme is stable, conservative and reliable for the numerical simulation of hyperbolic conservation laws. Copyright © 2003 John Wiley & Sons, Ltd.     

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.