基于非结构网格的TTGC有限元格式的实现及在超声速流动中的应用

基于非结构网格系统,实现了时空三阶精度的TTGC有限元格式,并在三阶TTGC格式上发展了基于人工粘性的激波捕捉技术。在非结构网格下,采用这种方法对若干典型的超声速流动问题(SOD激波管、马赫数为3的前台阶流动以及马赫数为8的高超声速圆柱流动)进行了验证计算。结果表明,TTGC格式分辨率高,在粗糙网格下能够准确的模拟超声速流场中的激波、接触间断等复杂流动现象,并且能有效的控制间断附近的数值色散现象。与传统的有限体积方法相比,本文实现的TTGC有限元格式在模拟超声速流动问题方面具有格式精度高、数值耗散小等优点。     

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

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

一种健壮的混合Roe黎曼求解器

使用Roe格式计算多维流动问题时,在强激波附近会出现数值激波不稳定现象。带有剪切粘性的HLLEC格式不仅可以捕捉接触间断,而且表现出很好的稳定性。混合Roe格式和HLLEC格式来消除数值激波不稳定性。在强激波附近,通过激波面法向和网格界面法向的夹角来定义开关函数,使得数值通量在激波面横向切换成HLLEC格式。在其余地方,数值通量依然使用Roe格式来计算。数值试验表明,混合格式不仅消除了Roe格式的数值激波不稳定性,还最大程度地减少了HLLEC格式所带来的剪切耗散,保留了Roe格式高分辨率的优点。     

一种简单的精确捕捉接触间断的黎曼求解器

基于Godunov型数值格式的有限体积法是求解双曲型守恒律系统的主流方法,其中用来计算界面数值通量的黎曼求解器在很大程度上决定了数值格式在计算中的表现。单波的Rusanov求解器和双波的HLL求解器具有简单、高效和鲁棒性好等优点,但是在捕捉接触间断时耗散太大。全波的HLLC格式能够精确捕捉接触间断,但是在计算中出现的激波不稳定现象限制了其在高马赫数流动问题中的应用。本文利用双曲正切函数和五阶WENO格式来重构界面两侧的密度值,并且结合边界变差下降算法来减小Rusanov格式耗散项中的密度差,从而提高格式对于接触间断的分辨率。研究表明,相比于全波的HLLC求解器,本文构造的黎曼求解器不仅具有更高的接触分辨率,而且还具有更好的激波稳定性。     

网格自适应技术在复杂外形流场模拟中的应用

建立了一套适用于非结构混合网格自适应方法,针对激波和涡的不同特征采用不同加密探测器,各向异性加密棱柱单元并沿物面法向方向剖分所有棱柱层,各向异性剖分四面体单元,并保证四面体与棱柱交界面上网格协调。构造Hermit插值近似投影物面新加网格点和基于Laplacian光滑方法对空间网格进行优化。通过网格自适应加密,使用Roe格式计算高超声速球头绕流的红玉现象得到明显减轻。F16飞机含激波和脱体涡的流场自适应计算表明,网格加密集中在激波面和涡核附近区域,激波和涡计算更准确。     

超声速流动中非线性EASM湍流模式应用研究

针对超声速复杂流动区域精确模拟的需要,发展了基于k-ω可压缩修正形式的非线性显式代数雷诺应力模式（EASM）,提高了该模式对超声速复杂流动的数值模拟精度。通过对二维超声速凹槽和三维双椭球的数值计算表明,与SA和SST常规线性涡黏性湍流模式比较,非线性的EASM模式对大分离以及剪切层流动结构的刻画能力更精细,对剪切层再附区的压力及摩擦系数分布模拟更加精确;EASM模式能够准确地模拟二次激波引起的压强和热流分布情况。     

超声速中等欠膨胀冲击射流的涡声关系研究

使用大涡模拟方法对冲击面为平面的超声速中等欠膨胀冲击射流进行了数值模拟.利用三阶迎风和四阶对称紧致格式对无量纲化轴对称可压缩滤波N-S方程进行空间离散,时间上推进采用的是三阶精度的TVD型Rugge-Kutta法. 通过与经典的冲击射流实验比较,证实了程序的可靠性. 数值模拟得到了流场中不同尺度的涡结构和激波结构,观察到了上行声波和反射波以及流场中不同位置的声源,分析了冲击区剪切层附近区域的压强和涡旋转强度变化的频率、冲击平板上的压强变化频率以及射流剪切层中不同位置的涡合并出现的频率,发现冲击区剪切层附近区域的压强和涡强度变化以及射流剪切层中的涡合并现象和离散频率的冲击单音有重要关联.      

一种改进的中心-WENO混合格式

基于中心差分与WENO格式混合可以改善WENO格式耗散特性的思想,在理论推导的基础上,给出了一种用于激波捕捉计算的守恒型中心-WENO混合格式,该混合格式可视为三阶WENO格式和二阶中心差分格式的加权平均。在数值研究现有加权函数的基础上,给出了适用于该混合格式的加权函数,使其能够自适应地调整数值耗散以捕捉激波间断。数值结果表明:与三阶WENO格式相比,混合格式HY3_4能够降低数值耗散,更陡峭地捕捉间断,对复杂流场结构具有较高的分辨率;混合格式HY3_5对于包含高压比激波间断流场结构,能给出无振荡、低耗散的结果。     

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

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

超声速平面混合层小激波的形成与演变

为了揭示超声速混合层中小激波形成机理及其与涡相互作用的演变过程,本文基于大涡模拟（LES）方法,结合五阶精度混合TCD／WENO格式,对超声速平面混合层在对流马赫数为Mc=0．65条件下的流场结构进行了数值模拟,数值结果详细描述了超声速混合层中小激波的形成过程。研究了小激波形成后,随涡运动而产生的变形、脱落及发展过程。同时,对混合层双涡合并过程中,小激波与相邻涡相互作用所产生的变形与演变过程进行了讨论。     

STUDY ON THE FUEL AIR MIXING INDUCED BY A SHOCK WAVE PROPAGATING INTO A H_2-AIR INTERFACE

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｆｕｅｌａｉｒｍｉｘｉｎｇｉｎａｓｈｅａｒｌａｙｅｒｆｌｏｗｉｓａｎｉｍｐｏｒｔａｎｔｐｒｏｂｌｅｍｉｎｓｔｕｄｙｉｎｇｃｏｍｐｒｅｓｓｉｂｌｅｆｌｏｗａｎｄｓｕｐｅｒｓｏｎｉｃｃｏｍｂｕｓｔｉｏｎ ,ｓｕｃｈａｓｆｕｅｌｄｉｆｆｕｓｉｏｎａｎｄｍｉｘｉｎｇｉｎａＳｃｒａｍｊｅｔ[1].Ｒｏｓｈｋｏ[2 ]ｓｔｕｄｉｅｄｅｘｐｅｒｉｍｅｎｔａｌｌｙｔｈｅｐｈｅｎｏｍｅｎｏｎｏｆｆｕｅｌａｉｒｍｉｘｉｎｇｉｎａｓｕｂｓｏｎｉｃｓｈｅａｒｆｌｏｗａｎｄｆｏｕｎｄｌａｒｇｅ ,ｃｏｈ…     

三阶WENO格式精度分析与改进

高分辨率激波捕捉格式对含激波流场的数值模拟具有重要意义。通过理论推导分析了WENO-JS3格式和WENO-Z3格式的精度,发现两种格式在光滑流场区域（包含极值点处）具有相同的理论精度且均低于三阶设计精度,WENO-Z3格式由于增大了非光滑模板的非线性权重使其计算精度有所提高。在理论推导的基础上,提出了WENO-Z+3格式及其改进格式（WENO-Z+3P1和WENO-Z+3P32）,且改进格式在光滑流场区域能满足所设计的三阶精度要求。选用一维平面黎曼问题及双马赫反射等经典算例,验证了本文提出的WENO-Z+3格式及其改进格式相较其他格式具有耗散低和对流场结构分辨率高的特性。     

The formation of a secondary shock wave behind a shock wave diffracting at a convex corner

This paper deals with the formation of a secondary shock wave behind the shock wave diffracting at a two-dimensional convex corner for incident shock Mach numbers ranging from 1.03 to 1.74 in air. Experiments were carried out using a 60 mm 150 mm shock tube equipped with holographic interferometry. The threshold incident shock wave Mach number () at which a secondary shock wave appeared was found to be = 1.32 at an 81° corner and = 1.33 at a 120° corner. These secondary shock waves are formed due to the existence of a locally supersonic flow behind the diffracting shock wave. Behind the diffracting shock wave, the subsonic flow is accelerated and eventually becomes locally supersonic. A simple unsteady flow analysis revealed that for gases with specific heats ratio the threshold shock wave Mach number was = 1.346. When the value of is less than this, the vortex is formed at the corner without any discontinuous waves accompanying above the slip line. The viscosity was found to be less effective on the threshold of the secondary shock wave, although it attenuated the pressure jump at the secondary shock wave. This is well understood by the consideration of the effect of the wall friction in one-dimensional duct flows. In order to interpret the experimental results a numerical simulation using a shock adaptive unstructured grid Eulerian solver was also carried out. Received 1 May 1996 / Accepted 12 September 1996     

High‐order shock capturing schemes for turbulence calculations

This work investigates high‐order central compact methods for simulating turbulent supersonic flows that include shock waves. Several different types of previously proposed characteristic filters, including total variation diminishing, monotone upstream‐centered scheme for conservation laws, and weighted essentially non‐oscillatory filters, are investigated in this study. Similar to the traditional shock capturing schemes, these filters can eliminate the numerical instability caused by large gradients in flow fields, but they also improve efficiency compared with classical shock‐capturing schemes. Adding the nonlinear dissipation part of a classical shock‐capturing scheme to a central scheme makes the method suitable for incorporation into any existing central‐based high‐order subsonic code. The amount of numerical dissipation to add is sensed by means of the artificial compression method switch. In order to improve the performance of the characteristic filters, we propose a hybrid approach to minimize the dissipation added by the characteristic filter. Through several numerical experiments (including a shock/density wave interaction, a shock/vortex interaction, and a shock/mixing layer interaction) we show that our hybrid approach works better than the original method, and can be used for future turbulent flow simulations that include shocks. Copyright © 2009 John Wiley & Sons, Ltd.     

Assessment of TVD schemes for inviscid and turbulent flow computation

A systematic study has been conducted to assess the performance of the TVD schemes for practical flow computation. The viewpoint adopted here is to treat the TVD schemes as a combination of the standard central difference scheme with numerical dissipation terms. The controlled amount of numerical dissipation modifies the computed fluxes to ensure that the solution is oscillation-free. Four variants of TVD schemes, two with upwind dissipation terms and two with symmetric dissipation terms, have been studied and compared with the conventional Beam-Warming scheme for inviscid and turbulent axisymmetric flow computations. The results obtained show that all four variants can accurately resolve the shock and flow profiles with fewer grid points than the Beam-Warming scheme. The convergence rates of the TVD schemes are also substantially superior to that of the Beam-Warming scheme. The combination of high accuracy, good robustness and improved computational efficiency offered by the TVD schemes makes them attractive for computing high-speed flow with shocks. In terms of the relative performances it is found that the symmetric schemes converge slightly faster but that the upwind schemes are less sensitive to the number of grid points being employed.     

加权型紧致格式与加权本质无波动格式的比较

线性紧致格式和加权本质无波动格式是两种典型的高阶精度数值格式,它们各有优缺点.线性紧致格式在具有高阶精度的同时,格式的分辨率也比较高,耗散低,是计算多尺度流场结构的较好格式,但是不能计算具有强激波的流场.加权本质无波动格式是一种高阶精度捕捉激波格式,鲁棒性好,但耗散比较高,分辨率也不理想.近年来,在莱勒的线性紧致格式基础上,采用加权本质无波动格式捕捉激波思想,发展了一系列加权型紧致格式.本文较全面地比较了加权型紧致格式和加权本质无波动格式,包括构造方法、鲁棒性、分辨率、耗散特性、收敛特性以及并行计算效率.结果表明,现有的加权型紧致格式基本保持了加权本质无波动格式的性质,对于气动力等宏观量的计算,比加权本质无波动格式没有明显的优势.      

Role of the momentum interpolation mechanism of the Roe scheme in shock instability

The shock instability phenomenon is a well‐known problem for hypersonic flow computation by the shock‐capturing Roe scheme. The pressure checkerboard is another well‐known problem for low‐Mach‐number flow computation. The momentum interpolation method (MIM) is necessary for low‐Mach‐number flows to suppress the pressure checkerboard problem, and the pressure‐difference‐driven modification for cell face velocity can be regarded as a version of the MIM by subdividing the numerical dissipation of the Roe scheme. In this paper, MIM has been discovered through analysis and numerical tests to have the most important function in shock instability. MIM should be completely removed for nonlinear flows. However, the unexpected MIM is activated on the cell face nearly parallel to the flow for the high‐Mach‐number flows or low‐Mach‐number cells in numerical shock. Therefore, MIM should be retained for low‐Mach‐number flows and be completely removed for high‐Mach‐number flows and low‐Mach‐number cells in numerical shock. For such conditions, two coefficients are designed on the basis of the local Mach number and a shock detector. Thereafter, the improved Roe scheme is proposed. This scheme considers the requirement of MIM for incompressible and compressible flows, and is validated for good performance of numerical tests. An acceptable result can also be obtained with only the Mach number coefficient for general practical computation. Therefore, the objective of decreasing rather than increasing numerical dissipation to cure shock instability can be achieved with simple modification. Moreover, the mechanism of shock instability has been profoundly understood, in which MIM plays the most important role, although it is not the only factor. Copyright © 2016 John Wiley & Sons, Ltd.     

二维激波与剪切层相互作用的直接数值模拟研究

采用五阶weighed esseritially non-oscillatory (WENO) 格式和三阶total variation diminishing (TVD) Runge-Kutta 格式, 通过求解二维非定常Navier-Stokes 方程, 直接数值模拟了激波与剪切层相互作用, 目的在于揭示激波与剪切层相互作用过程中噪声产生的机理. 研究发现:(1) 当入射激波穿过剪切层时, 剪切层中心位置向下层区域偏移;(2) 入射激波穿过剪切层产生小激波, 在小激波与剪切层接触点处产生声波并向外辐射;(3) 反射激波穿过剪切层后形成了分段弧状激波;(4) 当反射激波穿过剪切层时, 激波在鞍点处泄漏并向外辐射声波, 这是一种激波泄漏机制.     

Mixing Characteristics in a Supersonic Combustor with Gaseous Fuel Injection Upstream of a Cavity Flameholder

Non-reacting experiments and numerical simulations have been performed to investigate the mixing characteristics in a supersonic combustor with gaseous fuel injection upstream of a flameholding cavity in a supersonic vitiated air flow with stream Mach number 1.7. Using helium as simulated fuel, the acetone vapor is adulterated into the fuel jet. The fuel distribution in spanwise and streamwise direction is imaged by the planar laser-induced fluorescence (PLIF) measurement. According to the similarity of experimental observations with different cavities, the typical L/D = 7 cavity with aft wall angle 45° is chosen and the flowfield with the injection is calculated by Large Eddy Simulation. Experimental and numerical results have shown that most of the fuel flow away upon the open cavity with the lifting counter-rotating vortex structures induced by the transverse jet. Only a small portion of the fuel is convected into the cavity shear layer by the vortex interaction of the jet with cavity shear layer, and then transported into the cavity due to the cavity shear layer motion and the interaction of the shear layer with the cavity trailing edge.