超燃进气道激波/湍流边界层干扰

针对超燃进气道湍流边界层/激波干扰引起的分离问题,采用基于5阶WENO数值格式的大涡模拟(LES)方法开展流场湍流非定常预测,旨在分析进气道湍流化技术实现进气道起动的可行性。研究表明,平板激波/湍流边界层干扰(STBLI)问题,LES方法能够清晰、可靠预测反射、分离激波形成过程及激波与充分发展湍流边界层的相互干扰,定量结果与试验一致;进气道研究方面,层流状态下,激波干扰产生强分离,导致进气道堵塞,而采用湍流化控制后试验和计算均表明流场分离明显减小,流场稳定且无明显堵塞现象,进气道可以起动,总压恢复系数达到要求,该结果表明,利用强湍流化减弱分离,实现进气道起动思想是可行的。      

基于LES方法的平板非定常激波/湍流边界层干扰研究

以高超声速发动机进气道湍流分离控制为应用背景,采用大涡模拟(LES)方法进行马赫数为3.0(唇口附近马赫数约为3.0)的激波/湍流边界层干扰(SWTBLI)流场机理研究.利用扰动循环引入的方法,先得到充分发展湍流场,然后根据斜激波关系式引入激波的方法进行激波/湍流干扰模拟.研究结果显示:充分发展湍流场在激波作用下产生逆...     

TSTO并联分离激波/边界层干扰流动特性分析

针对两级入轨飞行器的缩比模型，通过试验与数值模拟相结合的方式，在马赫数6条件下开展典型级间距状态的激波/边界层干扰流场研究，详细分析干扰区壁面及空间的流动结构与特性。结果表明：试验中模型壁面边界层在激波入射之前为层流状态，在强激波干扰后迅速转捩为湍流状态，因此试验结果在第一道激波作用结束之前与层流计算结果吻合，而在第一道激波作用结束之后与湍流计算结果一致；激波/边界层干扰呈现复杂的三维流动特征和明显的开放结构，强激波在壁面形成的高压区呈弧状向下游展开，轨道级头部产生的入射激波在级间来回反射，强度依次递减；同时，干扰区内存在展向弯曲的主分离线与再附线、沿流向排列的二次分离线与再附线、流动剪切形成的旋涡结构以及包括鞍点、结点、焦点在内的临界点；层流边界层受到激波作用形成的分离区明显大于湍流边界层，同时开放特征更为显著。     

SST湍流模型的激波与可压缩效应改进

由于缺乏对某些重要流动特征的考虑,针对不可压流发展的标准SST湍流模型在描述超声速流场时存在明显的局限性。为改善SST模型在吸气式高超声速推进系统内流等复杂超声速流场中的预测精度,基于流动特征结构定向开展了激波和可压缩效应改进。通过激波/湍流边界层判别函数和可压缩效应判别函数定位标准SST模型参数或建模假设失效的区域,针对性地改进湍流模型。采用超声速平板边界层流动、超声速压缩拐角分离流动、超声速隔离段复杂激波串流动以及HIFiRE-2超声速内流等算例进行了测试,结果表明改进模型具有与标准SST模型一致的边界层预测能力,但显著提高了对激波干扰流动及逆压分离流的预测精度。     

运用GAO-YONG湍流模型对扩压器内跨声速流动的数值模拟

运用GAO-YONG可压缩湍流方程组,采用同位网格SIMPLE算法,对扩压器跨声速流动中的二维激波/湍流边界层干扰现象进行了数值模拟。将计算得到的流场的时均参数与实验值进行比较,数值模拟结果在激波强度、壁面压力分布以及分离点和再附点位置等方面,与实验值吻合较好,表明GAO-YONG可压缩湍流方程组能够比较准确的模拟较强激波/湍流边界层干扰流动,从而进一步为GAO-YONG湍流模型的正确性及其在可压缩流场模拟方面的适用性提供了佐证。     

混合动理学通量WENO方法拓展

在马赫数2.5来流条件下，开展了高频微秒脉冲放电控制压缩折角激波/边界层干扰非定常性的风洞实验，放电位于压缩折角上游沿流向布局的6对电极之间，所选取的放电频率为14 kHz,接近于来流边界层的特征频率。采用高速纹影成像技术记录流场的动态变化，并基于纹影图像灰度值的时间序列采用平均、均方根、本征正交分解、动态模态分解、傅里叶变换等方法进行处理，对比研究有/无控制情形下激波/边界层干扰的非定常特性。研究发现，对于无控情形的基准流场，流动的低频特性表现为分离激波的振荡及边界层大尺度涡经过激波的脱落行为，中、高频特性表现为边界层小尺度涡与激波的相互作用；对于受控情形，来流边界层内的大涡尺度在放电作用下增大，大尺度涡与分离激波相互作用使得激波的振荡转变为稀疏压缩波的脉动，流经激波的边界层脉动更强，分离激波的低频振荡(10～300 Hz)有所改善。此时，流动的低频特性主要表现为边界层大尺度涡经过激波的脱落行为，而中、高频特性与基准流场相似。     

高马赫数下激波湍流边界层干扰数值研究

应用GAO—YONG可压缩湍流方程组数值模拟了入射斜激波／平板湍流边界层相互干扰现象,计算了来流马赫数为5．0,激波入射角度分别为15．876°、23．287°两种不同激波干扰强度下的流场。计算程序中的对流项、扩散项分别采用二阶ROE格式和二阶中心差分格式离散,并用多步Runge—Kutta显式时间推进法求解空间离散后的控制方程。计算较好地模拟了高马赫数下的激波／湍流边界层干扰的流场结构,位移边界层厚度,动量损失厚度等,也比较准确地预测了平板壁面压力、摩阻系数等气动力参数的分布。     

边界层转捩与压缩拐角分离流动的非定常作用

为了认识激波与转捩边界层之间的相互作用,选取压缩拐角模型为研究对象,在高超声速风洞中开展了激波/转捩边界层干扰试验研究。试验在FD-20炮风洞中进行,试验马赫数为8,雷诺数0.5×107~2×107/m。试验采用了薄膜电阻温度计和常规压力传感器,分别测量压缩拐角周围的热流和压力分布。根据干扰区上下游的边界层流态,将试验分为两部分:层流/湍流干扰和转捩/湍流干扰。对比分析了边界层转捩发生在干扰区时的热流分布、压力分布以及脉动热流的非定常特性。研究结果表明,激波/转捩边界层干扰的热流和压力分布特征,不同于常规的层流干扰和湍流干扰,其介于二者之间。层流/湍流干扰的热流和压力分布特征类似于层流干扰;转捩/湍流干扰的热流和压力分布特征类似于湍流干扰。互双谱分析结果表明,当边界层转捩发生在分离区时,转捩脉动与分离泡脉动同时出现并增长。当两者幅值足够大时,转捩脉动会与分离泡脉动发生非线性耦合作用。这种耦合作用会诱导出新的频率特征的脉动结构,从而使分离区内的脉动量显著增加。     

管内正激波/湍流边界层干扰的数值模拟

本文采用时均Ｎ－Ｓ方程和Ｂａｌｄｗｉｎ／Ｌｏｍａｘ代数湍流模型计算了典型拉伐尔喷管内正激波与湍流边界层干扰流场。计算与实验结果的比较表明，方法可准确地预测激波结构，激波与边界层干扰区流动特征，激波位置、激波前马赫数和壁面压力分布等。     

动态微涡流发生器对激波/边界层干扰的影响研究

本文基于OPENFOAM数值仿真平台,采用动态网格技术和湍流离散涡(DES)模型,研究了微涡流发生器以一定速度向下游移动时,激波/边界层干扰(SWBLI)流场特性的变化,重点关注干扰区域内的流向和展向的流场特性。来流马赫数为4,微涡流发生器向下游移动速度为0m/s,20m/s和40m/s。研究表明：当MVG向下游移动时,SWBLI区域的“弓”形高压区会演化成“双弓”形；入射激波形成高压区的压力明显降低,同时,入射激波和反射激波形成高压区的峰值位置均会向下游移动；流场下游 “双圆弧”状高压区的高度逐渐降低；SWBLI区域边界层的高度逐渐降低,同时边界层底部的速度也有所降低；随着MVG移动速度的增加,对SWBLI流场的控制效果更加明显；动态MVG对流场的控制是通过尾迹涡和波系结构实现的。     

高能电弧等离子体激励控制双压缩拐角激波/边界层干扰实验研究

超声速/高超声速飞行器进气道入口处多采用多级压缩构型,其诱导的激波/边界层干扰严重影响进气道效率和飞行器的气动性能,因此对双压缩拐角激波/边界层干扰进行流动控制具有较强的应用背景。在来流速度为Ma=2.0的风洞内,针对三种典型的双压缩拐角构型,开展了高能流向脉冲电弧放电阵列调控双压缩拐角激波/边界层干扰的实验研究,并对激励流场的高速纹影图像进行了空间梯度阈值处理和均方根处理。结果显示,在激励的作用下两道分离激波的强度均减弱,验证了利用高能流向脉冲电弧放电阵列控制双压缩拐角激波/边界层干扰的可行性。在分析控制效果的基础上,获得了在不同构型拐角的流场中前驱冲击波列和控制气泡的演化规律,结合控制效果的时序特征,最终揭示了高能流向脉冲电弧放电阵列作用于双压缩拐角激波/边界层干扰的前驱控制和接续控制的接力控制机理。     

超声速湍流边界层与圆柱相互作用试验

在超声速风洞中开展了湍流边界层与圆柱相互作用流场研究,试验马赫数为3.4和3.8.圆柱安装在试验段底板上,安装位置的边界层为充分发展的湍流边界层,研究了圆柱直径和高度对流场结构和压力脉动的影响.采用基于纳米示踪的平面激光散射(NPLS)技术获取了流向和展向流场精细结构,激波系和马蹄涡结构均可清晰分辨.通过展向流场图像可...     

超声速压缩拐角激波/边界层干扰动力学模态分解

压缩拐角激波与边界层干扰问题广泛存在于高速飞行器的外部和内部流动中,其非定常复杂流场结构对飞行器气动性能影响显著。动力学模态分析将有助于进一步加深理解激波与边界层干扰流场不同特征频率对应的流动结构及动力学特性,为揭示其复杂流动机理提供参考。本文采用动态模态分解(DMD)方法对来流马赫数为2.9、24°压缩拐角内激波与超声速边界层干扰下的非定常流动进行了模态分析。评估了稀疏改进动态模态分解方法在压缩拐角流动中的适用性,研究了湍流干扰和转捩干扰下典型特征频率对应的动力学模态空间结构差异及其原因,分析了转捩边界层展向非均匀性对低频/高频模态动力学机制的影响规律。研究发现,湍流干扰与转捩干扰下拐角干扰区内均存在两类截然不同的动力学模态:低频模态和高频模态。低频模态结构集中在分离激波及分离泡剪切层的根部,表征为分离泡的大尺度膨胀和收缩运动;高频模态空间分布则以平均声速线附近正负交替结构为主,对应为边界层内不稳定波沿剪切层往下游的传播。转捩边界层的展向结构对低频模态运动特性影响明显,而对高频模态的影响则相对较小。     

高分辨率激波/边界层干扰时间演化过程分析

在超声速风洞中开展了湍流边界层与圆柱相互作用流场研究，试验马赫数为3.4和3.8。圆柱安装在试验段底板上，安装位置的边界层为充分发展的湍流边界层，研究了圆柱直径和高度对流场结构和压力脉动的影响。采用基于纳米示踪的平面激光散射(NPLS)技术获取了流向和展向流场精细结构，激波系和马蹄涡结构均可清晰分辨。通过展向流场图像可以发现干扰区内激波与湍流结构的相互作用具有明显的非定常性。采用动态压力传感器测量了圆柱前方相互作用区域的压力脉动特性，在激波足区域压力呈现11～38 kHz的宽频分布，推测主要由激波足与涡结构相互作用及滞止区涡结构的破碎引起。随着圆柱高度的增加，激波足附近测点对应的特征频率有所降低；上游测点则发现了0～3 kHz低频区能量的增强，这主要是由分离区引起的，表明在一定高度范围内高度的增加增强了流动分离。     

进气道中激波-边界层干扰特性研究

针对高超进气道中激波/边界层干扰等现象,选取了24°压缩拐角及二维高超进气道DLR-GK01等算例,比较研究了混合网格并行数值模拟软件中SA、SST、TNT三种湍流模型对激波/边界层干扰等现象的模拟能力,主要考察了不同湍流模型对流场结构、壁面压强以及壁面摩阻的影响。结果表明:三种湍流模型在模拟大分离流场时存在不同的流动分离特性,TNT湍流模型较SST和SA湍流模型模拟的分离起始位置、分离区大小及边界层厚度与实验吻合度最好。     

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

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

超声速膨胀角入射激波/湍流边界层干扰直接数值模拟

为了揭示膨胀效应对激波/湍流边界层干扰区内复杂流动现象的影响规律,采用直接数值模拟方法对来流马赫数2.9、30°激波角的入射激波与10°膨胀角湍流边界层相互作用问题进行了数值研究。系统地探讨了激波入射点分别位于膨胀角上游、膨胀角角点和膨胀角下游3种工况下膨胀角干扰区内若干基本流动现象,如分离泡、物面压力脉动及激波非定常运动、湍流边界层统计特性和相干结构动力学过程等。结果表明,激波入射点流向位置改变对分离区流向和法向尺度的影响显著,尤其是当激波入射点位于角点及其下游区域。研究发现,膨胀角干扰区内物面压力脉动强度急剧减小,分离区内压力波向下游传播速度将降低而在膨胀区内将升高,膨胀效应极大地抑制了分离激波的低频振荡运动。相较于入射激波与平板湍流边界层干扰,入射激波流向位置改变对膨胀角再附区速度剖面对数区及尾迹区影响显著,将导致其内层结构参数升高而外层降低,近壁区内将呈现远离一组元湍流状态的趋势。此外,流向速度脉动场本征正交分解分析指出,主模态空间结构集中在分离激波及剪切层根部附近而高阶模态以边界层内小尺度正负交替脉动结构为主。低阶重构流场结果表明,前者对应为分离泡低频膨胀/收缩过程而后者表征为分离泡高频脉动。     

压力恢复系统扩压器激波串现象的数值模拟

通过求解由S-A一方程湍流模式封闭的三维雷诺平均Navier-Stokes方程,采用有限体积方法对激光器压力恢复系统变截面扩压器常温状态下的三维流场进行了数值模拟.数值计算和实验符合良好,较好模拟了扩压器中由激波/边界层干扰诱导的复杂流场的流场特性并再现了流场中的"激波串"和"伪激波"现象.本文对激波串的形成机理进行了分析,发现了几种变截面扩压器特有的激波串结构,并研究了扩压器扩张角对扩压效率的影响.     

马赫数对后掠激波和湍流边界层干扰特性的影响

本文介绍了尖前缘翼诱导激波和湍流边界层干扰流场壁面特性，着重强调马赫数影响。给出２．０≤Ｍ≤８．２、ａ≤３５°分离流场中，锥型干扰区内主分离线和再附线位置与无粘激波角β０和迎角ａ的相关式，证实无粘条件是控制锥型区尺度的主要因素，面高超声速与超声速干扰流中二次分离随激波强度的不同发展，表明干扰流场的细致结构与可压缩性有关。     

应用GAO-YONG湍流模式数值模拟三维激波/湍流边界层干扰

应用GAO-YONG可压缩湍流模式数值模拟了三维激波/湍流边界层干扰算例之一——单鳍流动.攻角20°,来流马赫数2.93,雷诺数9.8×105.对流项和扩散项分别采用Roe格式和二阶中心差分格式计算.Runge-Kutta显示时间推进方法求解了半离散的控制方程.包括壁面压力分布,边界层内流动偏移角等在内的计算值与试验数据进行了比较.准确地预测出了三维激波/湍流边界层干扰流场的主要流动特性——λ波结构,主分离涡核,膨胀区,滑移线等.计算与Alvi等提出的单鳍流动的理论模型符合很好,得到了平板表面压力以及分离线、再附线等在单鳍流动中所独有的半圆锥特性.