动态间断装配法模拟斜激波壁面反射

基于非结构动网格技术和边界装配思想提出了动态间断装配法,该方法能够应用于求解含有间断的流动问题。无论入射激波还是反射激波都是作为边界进行处理,激波运动速度由兰金-许贡纽(Rankine-Hugoniot)关系确定。激波作为动网格的一部分,其运动由动网格技术实现。采用该方法模拟了超声速二维流场中激波与壁面相交问题,并且与捕捉法进行比较,二者的流场结构符合良好,但是在细节上还是存在明显差异。通过对流动结构的分析,得出采用装配方法得到的流场要优于捕捉方法的结论。激波壁面反射的问题模拟,也说明了边界激波装配方法对于复杂的激波相交问题是具有处理能力的。     

计算流体动力学在航天飞行器射流逆流中的应用

超声速射流逆流通常用于导弹、航天飞机、卫星和飞船等飞行器运动状态的控制。欠膨胀超声速射流逆流的流场包含有多激波（如弓形激波和马赫盘）、接触间断和剪切层，其结构非常复杂。本文采用激波高分辨率有限差分（ＴＶＤ）格式，对恒定自由流条件，各种不同射流出出口压比的超声速轴对称逆射流进行了数值模拟，且对各种条件下的物理现象给予了分析。计算的马赫盘和弓形激波位置与实验值相吻合，为此类流动问题提供了一种有效的预测     

含激波、旋涡和声波的复杂多尺度流动数值模拟研究

近年来针对激波噪声计算,我们开展了较为系统性的数值方法研究,解决了计算激波噪声的一些方法上存在的问题。通过对典型问题的直接数值模拟,包括旋涡之间相互作用、激波与旋涡相互作用、激波与剪切层相互作用等模型问题和超声速喷流激波噪声问题,揭示了一系列激波噪声产生机理。本文对这些研究进行简要回顾。     

基于CFD的高效气动弹性计算方法

本文在高效气动弹性分析方法方面的研究工作包括:对于超声速流动范围,发展了一种基于CFD数值计算流场修正的高效、高精度工程算法,既含有CFD技术模拟复杂激波系的流场特征信息,又具备工程算法计算速度快的优点;对于失速迎角前的飞行器从亚声速、跨声速到超声速范围,发展了一种基于CFD技术的非定常气动力建模方法,由此获得气动弹性计算所需的非定常气动力,此方法可以包含物面在跨声速存在激波运动的流场信息,计算精度与直接应用CFD模拟气动弹性过程的方法相当,但运用非定常气动力建模技术途径的效率可以比直接基于CFD模拟方法提高2个数量级,对CFD技术应用于气动弹性精确预测方法具有工程实际应用价值;还研究了基于非定常气动力建模的方法用于伺服气动弹性分析、阵风响应.     

关于跨音速叶栅物面边界条件的数值试验

本文讨论了跨音速叶栅物面边界的数值处理方法,目的是寻求与我们所用内点格式相匹配的边界上流动参数的算法,以提高计算精度。本文首先推导了适合于物面边界计算的特征关系式,其次总结了五种边界处理方法,进行了数值试验。计算结果表明,配合我们的内点计算格式,以y向抛物外推的Denton修正法计算效果较好,能够较精确地反映跨音速叶栅带激波流场的物理特点。     

基于各向异性非结构网格的超声速流动自适应计算

基于各向异性非结构网格,实现了自适应求解技术,将其应用于超声速楔形体绕流及超声速横向喷流问题的计算研究,对各向异性网格自适应计算过程的收敛性和求解精度进行了对比分析,展示了各向异性网格自适应算法在降低问题求解规模。提高各向异性流场分辨率等方面的优势及存在的问题。超声速楔形体无粘绕流场的自适应计算,清晰地捕捉到了激波附近的流场信息,激波前后的马赫数、压力、密度。温度等物理量与理论分析解吻合很好。超声速横向喷流流场存在激波、分离涡、边界层等流场结构的相互干扰。计算研究表明,单纯基于Hes-sian度量张量的各向异性网格生成及自适应算法不能有效模拟边界层内的流动情况,是将来需要进一步开展的研究挑战。     

三维欧拉流动模拟的基于伴随方程的网格自适应

<正>航空工程设计中,通常需要数值模拟给出较高精度的气动力计算结果,而目标函数的误差往往是未知的。在一套极密的网格上求解控制方程当然可以保证计算精度,然而这会造成极大的计算量。自适应网格技术通过局部加密,细化关键区域的网格,保证计算精度的同时提高了计算效率。网格自适应处理的准则通常是基于流场梯度构建的,即在流场变化较快的地方进行局部加密。这将导致包含诸如激波、驻点、旋涡等物理现象的区域的网格被加密。     

MVG控制斜激波/湍流边界层干涉的大涡模拟

采用浸没边界法(IBM)对带有微型涡发生器(MVG)控制器的激波/湍流边界层干涉流动进行了大涡模拟(LES)。以来流马赫数为2.3的斜激波(由平板上方8°楔产生)为基本流动入射平板湍流边界层,通过在干涉区前布置MVG阵列来控制激波诱导的边界层分离。采用浸没边界法处理MVG的复杂几何,分析了MVG尾迹区平均流速度剖面,雷诺应力,瞬态旋涡结构。结果表明:时均流场显示MVG尾迹区存在一对对转的主流向涡,流向涡加剧了边界内的动量交换从而增加了边界层抗分离能力,而瞬态流场则反映出MVG尾迹区的剪切层由于Kelvin-Helmholtz(K-H)不稳定性会卷起为一列展向旋涡。     

弹头外形不同尺度模型对二维射流流场的影响

超声速射流装置能很好地模拟榴弹引信发电机飞行时大气环境的工况，并成为地面考核和验收必不可少的设备。从二维轴对称N-S方程出发，利用二阶精度有限体积法和S-A湍流模型，对超声速射流装置的自由射流流场以及不同模型尺寸对超声速流场结构的影响进行了模拟计算。给出了激波与射流边界之间复杂干扰流场结构图，结果证明不同模型的尺寸对超声速射流流场有较大的影响。     

斜切反喷管性能计算

本文提出了斜切反喷管流场计算方法,内部点用MacCormack二步显格式数值解控制方程,入口边界和固壁边界的参数用物理边界条件和可应用的沿双特征线的相容性方程计算.固壁边界计算中采用了激波装配技术和开关格式,有效地捕获了流场内的激波系列,计算中得到的激波情况和物理分析一致.本文在流场计算的基础上,进行了反喷管的性能计算.本文提出的方法可应用于反喷管的性能预估.     

超薄双层复合结构界面特性超声兰姆波评价方法

提出超声兰姆波方法来评价超薄双层复合结构界面联接特性 ,这里“超薄”的概念是指层厚度 h远小于层中的声波波长λ,导致采用脉冲回波方法时 ,在时域上层前后界面的各次回波信号相互混迭。分析了超薄双层复合结构界面处于刚性联接、滑移联接和完全脱层 3种情况下超声兰姆波的色散特性 ,提出了以色散曲线为依据的超声兰姆波界面特性评价方法 ,进一步分析了影响方法估计准确性的各种因素。     

嵌边法出流条件在可压缩流直接数值模拟中的应用

以嵌边法为出口边界条件试算了亚、超声速平板边界层扰动波的演化、超声速流中涡的传播、亚声速平板边界层尾缘扰动的演化等,以验证嵌边法作为出流条件在可压缩流中的有效性。然后,探讨了嵌边法有关参数的选取问题。     

激波冲击下Air/SF6斜界面不稳定性实验研究

激波在不同密度介质上的交互作用在可压缩湍流上具有重要的基础价值。激波在界面上的作用会引起Richtmyer-Meshkov不稳定性。激波不正规折射时,流场存在更多复杂的涡。研究马赫数为1.23、1.41的激波在初始倾角β=60°的Air/SF6界面上非正规折射的情况。入射激波的切向冲击和法向冲击的相互作用,在界面处产生涡,折射波在壁面发生马赫反射。利用阴影显示技术,给出了界面演化和混合的过程。     

进气道动态畸变被动控制技术的初步研究

对S型进气道扩压器在临界及超临界工况下的流动特性进行了研究,通过研究了解超临界工况下动态畸变的形成和变化特征.国内外的研究均表明:进气道内激波与附面层的相互干扰是引起动态畸变的主要原因之上,因而本文中采用壁面处理被动控制技术进行激波与附面层干扰的控制,用以抑制动态畸变.研究结果表明:在较大的超临界工况下,该控制技术可明显地抑制激波的失稳和振荡,从而使动态畸变的气流脉动峰值下降约60%.     

具有无源控制空腔时正激波/湍流附面层干扰的数值模拟

采用雷诺平均N-S方程和B/L代数湍流模型计算了具有无源控制空腔时正激波/湍流附面层干扰流场。计算与实验结果的比较表明,本文方法可较准确地预测激波结构、激波与附面层干扰区流动基本特征及波后流动分离状态、激波位置、波前马赫数等参数。      

高精度数值方法在DNS中的应用

利用不同格式模拟扰动波在平板边界层中的传播,考察了格式对T-S波传播过程计算结果的影响。研究了不同计算格式对自由剪切层扰动波增长的影响,结果表明:对于二维扰动波模型,采用二阶精度格式,扰动波在平板边界层中只能出现衰减的结果;在自由剪切层中,虽然能模拟出扰动波在传播过程中幅值的增长,但却在相当长的计算范围内看不见三维波的出现。但对于同样的扰动模型,如果采用高精度格式,在边界层中能出现增长的T-S波,并与线性理论非常吻合;在自由剪切层中,扰动波不仅迅速增长,并能在不很长的距离内,明显的看到三维波的出现,并激发出了三维结构。      

跨声速翼型绕流计算边界层反方法的研究

本文用有限差分方法,通过有粘/无粘迭代计算了二元翼型的跨声速绕流问题。在边界层粘性区域内考虑了层流、转捩及湍流流动,当边界层内出现分离时,使用边界层反方法,采用代数湍流模型。算例表明,对激波/边界层弱干扰和强干扰情况,该方法的结果与风洞实验结果吻合良好,对于求解边界层小分离流场是一种好的近似方法。     

Numerical and experimental investigation into hypersonic boundary layer transition induced by roughness elements

In this work, the Direct Numerical Simulation(DNS) and Oil-Film Interferometry(OFI)technique are used to investigate the hypersonic boundary layer transition induced by single and double roughness elements at Mach number 5. For single roughness, the DNS results showed that both horseshoe vortices and hairpin vortices caused by shear layer instability can affect the boundary layer instability. The generation of the near-wall unstable structure is the key point of boundary layer transition behind the roughness element. At the downstream of the roughness element, the interaction between horseshoe vortices and hairpin vortices will spread in the spanwise direction.For double roughness elements, the effect of the spacing between roughness elements on the transition is studied. It is found that the case of higher spacing between roughness elements is more effective for inducing transition than the lower one. The interaction between two adjacent roughness elements can suppress the evolution of horseshoe vortices in the downstream and trigger the instability of shear layer. Thus, the transition will be suppressed accordingly.     

应用S-A模型的自由射流和冲击射流数值模拟

为探索轴对称射流准确的数值模拟方法，根据可压缩的轴对称N-S方程，利用二阶精度的有限体积法对轴对称的自由射流和垂直冲击射流进行数值模拟。结果表明采用Spalart-Allmaras湍流模型，可很好描述射流流动。收缩喷嘴射流计算应从出口上游的内流开始，而不能在出口截面上给定声速条件。计算结果与实验吻合较好，所获得的喷嘴内外部流场的流动状况以及冲击射流参数分布，可用于优化喷嘴内部结构设计。     

考虑附面层影响的二元混压式进气道设计

采用等激波强度的方法,考虑附面层修正,设计了一种飞行马赫数Ma=3.0的二元混压式进气道.通过数值仿真,模拟了激波-边界层的相互影响,研究了附面层抽吸对内流场的影响,获得了进气道内部复杂的流场分布,以及不同背压下进气道的起动特性.计算表明所设计的进气道性能较好,附面层抽吸对稳定正激波有明显的作用,提高了进气道抗反压能力.给出的方法可用于二元混压式进气道的初步设计和验证.