Bump进气道中鼓包诱导的激波/边界层干扰特性

为了探索Bump进气道中鼓包诱导的锥形激波和机身发展而来的湍流边界层干扰问题,分析其气动优势,首先选取了半锥和半棱锥这两种与鼓包的流场结构具有一定相似性的构型作为参照,采用数值仿真方法,分别对这三类典型的三维激波/湍流边界层干扰问题进行了流场分析。在此基础之上,设计了三个不同马赫数的鼓包,并研究了设计马赫数对鼓包流场特性的影响。结果表明:当三类构型的无黏激波强度相等时,半锥诱导产生的旋涡强度最强,鼓包次之,半棱锥最弱。尽管鼓包诱导的流场非常复杂,其干扰流场却呈现出准锥形相似的特性。虽然半锥对边界层的排移能力最强,但是综合考虑边界层排移能力及进气道出口流场畸变下,鼓包最具优势,这也是其被选为超声速进气道前缘压缩面的重要原因之一。此外,在设计状态下,适当增加设计马赫数能改善鼓包排移边界层的能力,但设计马赫数太高,边界层排移能力基本不变,反而使得进气道总压损失急剧增加。      

激波／边界层干扰脉动压力测量技术

建立了激波／边界层干扰脉动压力测试系统。对一组后掠压缩角产生的激波／边界层干扰脉动压力流场进行测量,结果表明,脉动压力时均值及其分布与静态测量得到的结果一致。测量还发现在实验的锥形干扰范围内,未出现类似于二维干扰中由非定常低频振荡激波引起的间歇现象。     

后掠激波边界层干扰中Mach数对特性区影响的研究

 本文介绍了由后掠压缩角模型引起的激波和湍流边界层干扰的实验研究。实验雷诺数Re=2.42～2.47×10~7/m,Ma_∞=1.79,2.04和2.50。模型共15个,其后掠角变化范围是0°～60°,流向压缩角变化范围为10°～30°。实验结果表明,在本实验范围内,激波边界层干扰中的上游影响区都呈现出柱形区或锥形区特性;柱形区和锥形区之间的边界随来流Mach数减小向锥形区发展。该边界主要决定于无粘激波的形式。     

Strong interactions of incident shock wave with boundary layer along compression corner

The coherent structure and instability of the interaction of incident shock wave with boundary layer developing on a compression corner are experimentally studied. The experiments are carried out in a supersonic wind tunnel of Mach number 2. Particular attention is paid to shock patterns and unsteady shock motions induced by the separation bubble. The high-speed schlieren is used to visualize the flowfield evolution and to characterize the instability. The snapshot proper orthogonal decomposition of schlieren sequences is applied to investigate the primary coherent structure in the flowfield. Fast Fourier transform and continuous wavelet transformation are applied to characterize the instability. The results show that there are large-scale low-frequency oscillations of the shock waves and small-scale high-frequency pulsations in the separation region. The peak frequency of shock oscillation is mainly concentrated in the range of 100–1000 Hz. The pulsation of the small flow structure in the separation bubble is mainly concentrated above 12.5 kHz. Based on the results of experimental analysis, the preliminary mechanism of the large-scale instability of such interaction is obtained.     

高超声速双模块内转式进气道的流动特性研究-Part I： 设计状态

本文通过特征线法设计了基于多道激波+等熵压缩波的基准流场,在此基础上通过流线追踪法设计了双模块内转式进气道。通过数值仿真和风洞试验相结合的方法,获得了内转式进气道的内外流特性。研究结果表明：在内转进气道最大半径对应的角区位置存在大量的边界层堆积,受第二道激波/边界层干扰,在激波根部卷起锥形旋涡;在内转式进气道内部,唇罩激波和管道边界层干扰显著,管道内存在自唇罩指向压缩面的强周向压力梯度,从而诱导管道内边界层均往一处汇聚,卷起大尺度流向涡。仿真和试验结果表明在来流马赫数5.74,攻角0度状态下,进气道气动性能优良,出口总压恢复系数系数达到0.58,最大抗反压为112倍。     

前后缘型线同时可控的乘波体设计

提出了一种前后缘型线同时可控的乘波体设计方法,在马赫数可控的外锥形曲面基准流场中,结合流线追踪技术和混合函数,实现了椭圆前缘转椭圆后缘的乘波体设计,并在设计点(Ma=6.0)和接力点(Ma=4.0)对其进行数值仿真研究。前后缘同时可控的乘波体在型面剧烈过渡处产生了较弱的激波,出口两侧存在高温高压区,后部对称面附近的激波形状由圆弧变为平直线且出口处流场基本均匀,非常有利于与进气道匹配设计。另外,该乘波体具有较高的容积率和预压缩效率,附面层修正后的容积率为0.24,设计点时乘波特性较好,接力点时前部完全乘波,具有较高的升阻比,有黏条件下设计点和接力点的升阻比分别为2.54和2.41。此外,与给定前缘的乘波体相比,其升力、阻力、俯仰力矩和出口增压比都有明显增加,但是升阻比和出口总压恢复系数有所降低,在设计点无黏升阻比由3.56降为3.00。以上研究表明,本文的设计方法可行且更加灵活,拓宽了乘波体的选择范围。     

在激波边界层干扰中后掠柱形区的独立性原理研究

在后掠压缩角引起的三维激波边界层干扰中,柱形的上游影响区表现出准二维的流动特性。实验结果表明,在柱形区中通常的边界层独立性原理不再适用。但是,如果考虑由后掠引起的横流效应修正,柱形区中的上游影响距离的尺度和相应的二维干扰上游影响区尺度之间可以得到很好的相关关系。     

一种基于聚类分析的二维激波模式识别算法

在激波捕捉求解器计算的可压缩无黏流场基础上,提出了一种探测并识别二维激波干扰模式的新算法,从3个层面详细介绍了该算法的实施流程。首先,采用基于当地流场参数设计的传统激波探测方法,辨识出激波附近的一系列网格单元;其次,通过经典的K-means聚类算法将这些激波单元划分成许多簇,并根据簇的相邻信息定义每个簇的类别;最后,设定相关准则对某些紧邻的簇进行合并,进而确定各个激波干扰点的位置,记录各条激波分支所对应的簇,采用Bézier曲线拟合算法分别对其聚类中心进行拟合以获取更加光滑的激波线。数值试验表明,该算法不受网格类型的限制,不仅可以保证最终拟合的激波线具有较高的位置精度,还可以清晰地识别出流场中多激波干扰的模式,同时对分析非定常流场中激波的运动与演化过程也提供了一种有效的可视化手段。     

受限空间内激波与火焰作用的三维计算

激波与火焰的相互作用常发生在超声速燃烧和燃烧转爆轰过程中,为深入了解这一现象,采用带化学反应的三维Navier-Stokes方程,对平面入射激波及其反射激波与火焰的作用进行了计算研究,其中燃烧过程使用单步反应模型描述。研究结果显示:在受限空间内模拟激波与火焰作用,能更好地符合实验结果,从而体现出受限空间的三维效应;火焰在入射激波的作用下主要经Richtmyer-Meshkov不稳定而发生变形,此时火焰变形以物理作用为主,燃烧膨胀效果相对不明显;当反射激波与变形火焰再次作用后,火焰迅速膨胀变形,放热率维持在较高水平,此时化学反应过程起主要作用;在反射激波的作用下,变形火焰复杂三维涡结构的形成能强化热量与质量的输送,提高燃烧速率。     

变截面管道中激波诱导的两相流动

本文研究扩张型与收缩型变截面管道中激波诱导的非定常两相流动。对于稀相气固悬浮体采用了双连续介质模型,对于两相流动按照准一维近似处理。控制方程则是利用算子分裂技术和二阶GRP方法数值求解。文中给出了波后气固两相的流动结构,并讨论了粒子对两相流场的影响。     

不同波系配置的鼓包压缩面引起的激波/边界层干扰特性

为了对不同波系配置的鼓包压缩面引起的锥形激波和来流边界层干扰的流场特性进行对比分析,分别以总偏转角相同的单锥、双锥和等熵锥等3种外锥流场为基准流场,在来流马赫数为2.0的条件下,采用流线追踪法设计了3种鼓包压缩面。在此基础上,采用数值仿真的方法对这3种鼓包压缩面在不同马赫数下的三维激波/边界层干扰流场进行了对比分析,结果表明:相比于单锥鼓包,双锥和等熵鼓包流场中的涡量集中区明显减小,且涡核处的涡量值也有所降低;对于双锥和等熵鼓包,由于受到后续波系的干扰,其表面流动的拓扑结构不再呈现出锥形相似性;双锥与等熵鼓包具有和单锥鼓包相同的边界层排移能力,且这一能力不随来流马赫数的变化而变化。此外,采用多波系配置能提高鼓包外压缩系统的总压恢复系数,波系数目越多,提升的程度就越大,且这一优势在高马赫数下更加明显。      

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

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

楔—后掠柱激波干扰流场和热流的工程算法

本文针对航天飞机机翼前缘的简化模型,即楔-后掠柱模型,在大后掠、高超声速条件下,提供了一种计算柱体前缘激波干扰流场和热流的工程计算方法。在中国科技大学流体力学实验室的激波风洞中取得了激波干扰的流场显示结果。计算的结果与国内外实验数据的比较表明,这是一种有效的方法。     

高马赫数来流条件下斜激波与流向涡对相互作用

针对斜激波入射流向涡对形成的复杂流动干扰现象,在来流马赫数6条件下,利用后掠斜坡式涡流发生器生成流向涡对,采用数值模拟结合激波风洞实验,研究了斜激波与流向涡对相互作用引起的激波变形和旋涡演变。结果表明,斜激波与流向涡对相互作用的剧烈程度大致分为强、中、弱三种。在强相互作用下,激波面向上游凸起的程度最剧烈,呈现类“T”字型特征,其对称面处接近于一道正激波,波后出现回流区和局部流动滞止,使得流向涡对发生破碎;在中等相互作用下,对称面处存在局部正激波,但没有出现回流区;在弱相互作用下,不再出现局部正激波。进一步地分析发现,斜激波与流向涡对相互作用产生的流向涡对破碎现象,不符合经典的激波作用下单个流向涡破碎理论的预测。采用流向涡对对称面处的气流参数,对经典的流向涡破碎理论进行修正后,能够预测本文流向涡对在斜激波作用下是否发生破碎。本文为掌握激波作用下流向涡对破碎现象的规律,提供了重要参考。     

Numerical study of compression corner flowfield using Gao-Yong turbulence model

A numerical simulation of shock wave turbulent boundary layer interaction induced by a 24° compression corner based on Gao-Yong compressible turbulence model was presented.The convection terms and the diffusion terms were calculated using the second-order AUSM (advection upstream splitting method) scheme and the second-order central difference scheme,respectively.The Runge-Kutta time marching method was employed to solve the governing equations for steady state solutions.Significant flow separation-region which indicates highly non-isotropic turbulence structure has been found in the present work due to intensity interaction under the 24° compression corner.Comparisons between the calculated results and experimental data have been carried out,including surface pressure distribution,boundary-layer static pressure profiles and mean velocity profiles.The numerical results agree well with the experimental values,which indicate Gao-Yong compressible turbulence model is suitable for the prediction of shock wave turbulent boundary layer interaction in two-dimensional compression corner flows.      

非对称来流隔离段流动特性研究

1引言隔离段是超燃冲压发动机的一个重要气动部件,在进气道与燃烧室之间构建一气动热力缓冲区域,为进气道提供一个较宽的连续工作范围。通常超燃冲压发动机位于高超声速飞行器下表面的中后部,这样自由来流需要经过一段较长的前体后进入进气道,这就造成进气道进口靠近机体一侧存     

Reflection of a moving shock wave over an oblique shock wave

The reflection of a moving shock wave over a wedge immersed in a still gas and the reflection of a wedge induced steady shock wave over symmetrical and asymmetrical reflecting surfaces have received intensive considerations since more than 70 years ago. Here we consider a different shock reflection problem—reflection of a moving shock wave over an initially steady oblique shock wave induced by a wedge immersed in supersonic flow. For the flow condition we considered, five moving triple points, with each connecting an incident shock wave, a reflected shock wave and a Mach stem, are identified. By using the reference frame co-moving with each triple point, the type of each shock wave of this triple point is clarified. The present study is significant in that it treats a new shock reflection problem leading to a new shock reflection configuration and showing potential applications in supersonic flow with unsteady shock interaction.     

隐式紧耦合SST和TNT湍流模型的高速流动数值模拟

将SST(shear stress transport)和TNT(turbulent/non-turbulent)湍流模型输运方程与平均流场控制方程进行隐式紧耦合求解,结合当地时间步长方法和湍流源项隐式处理确保求解过程的快速和稳定.采用AUSMPW+(AUSM by pressure-based weight functions)格式和LU-SGS(lower-upper symmetric Gauss-Seidel)隐式紧耦合方法对高超声速压缩拐角流动、锥柱裙流动和超声速非对称激波/边界层干扰问题进行了数值模拟.计算结果与实验值的对比表明:SST模型和TNT湍流模型可以很好地预测15°压缩拐角流动的壁面压力和热流密度;随着压缩拐角的增大,计算结果与实验值偏差增大;可压缩性修正对压缩拐角流动的压力和热流密度分布有很大影响,对超声速非对称激波/边界层干扰流动影响很小;隐式紧耦合方法比显式紧耦合方法具有更好的收敛特性.      

激波绕射方台的非定常流场研究

本文以数值方法研究了运动激波绕射矩形突台的非定常流场发展过程。计算中采用ＴＶＤ格式，得到了清晰的流场变化图谱。从流场图中可分辨出旋涡的形式过程及分离域形态。研究结果表明，本方法可以很好地描述激波绕射方台形成的复杂流动过程，对流场间断有较高的捕获力。     

气流出口型面对凹面腔内激波聚焦起爆的影响

为提高两级脉冲爆震发动机工作的可靠性,探寻强化凹面腔内激波聚焦起爆爆震波的方法,通过数值模拟,采用氢气作为燃料、空气作为氧化剂,探讨了3种凹面腔气流出口型面及出口面积对激波聚焦起爆的影响。结果表明:采用垂直出口壁面有助于提高凹面腔内的激波聚焦起爆压力和温度,使起爆时刻提前,随着出口面积的减小,激波聚焦峰值压力有较大的提高;但效果不如垂直壁面,随着倾斜壁面出口面积减小,起爆点的压力和温度只有小幅度增加,起爆时刻基本无变化;采用弧形出口壁面,不能对向凹面腔底部传播的激波产生反射强化作用,不能有效起爆爆震波。对比3种出口壁面型面,垂直壁面能更有利于凹面腔内的能量聚集,强化激波聚焦起爆爆震波。