局部凸起对可压缩平板边界层稳定性的影响

在采用直接数值模拟(DNS)方法计算得到来流马赫数为4.5的带局部凸起的可压缩平板边界层基本流的基础上,基于线性稳定性理论(LST) 进行稳定性分析,通过求解扰动方程研究了不同频率的扰动波与不同高度的局部凸起的相互作用,定义了穿透系数来研究局部凸起对流动稳定性的影响.研究结果表明:当局部凸起高度小于0.2个边界层厚度时局部凸起对来流的影响是局部的;局部凸起使频率小于最不稳定波频率的扰动波更不稳定,对频率大于最不稳定波频率的扰动波有稳定作用,其中最不稳定波频率是指局部凸起中心位置处平板边界层中第2模态最不稳定波的频率;穿透系数定量刻画了局部凸起对稳定性的影响,高度为0.2个边界层厚度的局部凸起引起e-N方法中N值的修正量达到0.8,接近转捩位置处N值的10%,在转捩预测中必须考虑.      

超声速平板边界层中流向涡的生成条件

针对马赫数为4.5的超声速平板边界层,基于线性稳定性理论（LST）选取初始扰动组合,通过直接数值模拟（DNS）,计算了第一模态不稳定波扰动组合沿流向演化生成流向涡的过程。采用改进Omega-Liutex旋涡识别方法进行涡识别,结合流向不同位置截面的流线图,分析了流向涡的生成特性。根据流向涡在zy截面内的流线特征,提出了流向涡的生成条件,研究发现:流向涡可以直接通过一对展向对称的第一模态不稳定斜波扰动与基本流叠加得到,不是必须经过非线性作用。      

抛物化稳定性方程在曲面边界层中的应用

为了研究抛物化稳定性方程在曲面边界层中的应用,选取了等曲率圆环管道、等曲率板和NACA0012翼型3种典型的曲面边界层.通过计算小幅值扰动波的演化,抛物化稳定性方程的计算结果和线性稳定性理论、数值模拟扰动方程结果相符,说明可以使用抛物化稳定性方程研究曲面边界层的小幅值扰动波的演化和稳定性分析;通过计算有限幅值扰动波的演化,线性阶段抛物化稳定性方程计算结果和扰动方程相符,在非线性很强的阶段,抛物化稳定性方程计算发散,发散的位置可作为转捩位置.      

激波边界层的相互作用对扰动波传播的影响

利用线性稳定性理论和直接数值模拟研究了带有入射斜激波的、来流马赫数Ma=4．5条件下的平板边界层的失稳特性。重点考察了在由于激波边界层相互干扰,平板边界层上形成分离区,又进而产生激波、膨胀波和旋涡等复杂流动现象的流场上游,引入小扰动的TS波后,扰动波传播通过带有这些复杂流动现象的流场时,扰动波的发展变化特点。通过对流场中扰动波（包括基本波和衍生波）演化特征的分析,研究分离和激波等复杂流动现象对平板边界层稳定性的影响特点。数值模拟发现,激波的出现不同于一般的压缩波,在亚声速区与超声速区对扰动波演化的影响是不同的,此外,分离对扰动有稳定作用。     

平板边界层中展向波包型扰动引起的转捩

过去在平板边界层转捩及湍流的研究中,主要考虑的扰动在展向是均匀分布的,这样有利于研究,但在实际问题中扰动形式是多样的,边界层可能是三维的,扰动在展向是不均匀的。对于以往研究的扰动来说,三维平板边界层中的展向非均匀扰动是比较复杂的扰动形式,更接近自然转捩,因此研究这种扰动引起的转捩和湍流具有重要的实际意义。基于此,针对平板边界层,控制方程为无量纲化的Navier-Stokes扰动方程,时间上采用三阶精度的差分格式,空间上展向采用伪谱方法,流向和法向采用高阶精度差分格式,应用数值模拟的方法研究了小幅值和有限幅值展向波包两种情况。通过数值模拟和线性稳定性理论分析小幅值波包的演化,得到了小幅值波包的演化符合线性稳定性理论(LST);分析了有限幅值的展向波包型扰动引起的转捩和湍流,描述了物理空间和谱空间中波包型扰动的演化特征;同时针对不同展向位置进行分析,结果表明不同展向位置的转捩位置不同,但转捩过程和特征是类似的。     

超声速边界层三维扰动引起小激波的数值研究

通过直接数值模拟的方法,研究了M=4.5的超声速边界层中三维扰动的演化。以某一剖面为入口,加入一个及一对三维T-S波,发现随展向波数的增加,扰动幅值的增长率逐渐减小,证实了M=4.5的超声速边界层中,当三维扰动达到一定幅值时会有小激波出现,为建立可压缩流动稳定性理论提供了依据。     

超声速绕钝板熵层不稳定性的研究

为了研究熵层的不稳定性,首先选择马赫数为4.5、前缘为圆头的0°钝板模型,利用直接数值模拟得到基本流场,然后在此基础上进行线性稳定性分析,同时通过经验方法确定熵层和边界层区域。研究发现,对于某一流向站位的基本流剖面,存在两个广义拐点,分别位于边界层和熵层。针对熵层的广义拐点,进行稳定性分析研究了熵层中不稳定模态特性,得出该模态的不稳定区域集中在前缘附近,且扰动增长率很小。通过直接数值模拟熵层扰动的演化,发现在前缘附近扰动幅值也是缓慢上升,之后缓慢下降。另外,距离前缘较远处的数值模拟得到的扰动幅值与LST结果比较,二者吻合得很好,说明熵层在下游具有很好的平行性。     

分布式粗糙度对马赫数为4.5的平板边界层稳定性的影响

采用数值模拟求解Navier-Stokes(N-S)方程的方法研究了小幅值分布式粗糙度对马赫数为4.5的平板边界层中扰动演化的影响.用泰勒展开线性化的办法将边界条件取在光滑平板壁面上,来模拟小幅值分布式粗糙度.研究了粗糙度及扰动的不同因素对边界层扰动演化的影响,其中包括粗糙度高度、波长、粗糙度区域的长度及扰动波频率等.研究结果表明:分布式粗糙度通过改变边界层基本流场影响扰动波的幅值演化,一定幅值的分布式粗糙度会使粗糙度区域前的不稳定区向低频移动,使粗糙度区域后的不稳定区向高频移动.从总体效果上看,对于通过粗糙度时,靠近中性曲线下支界相对低频的扰动,粗糙度会抑制其幅值的增长;对于靠近中性曲线上支界相对高频的扰动,粗糙度会促进其幅值的增长.      

壁面小折角对马赫数4.5边界层中扰动演化的影响

为了利用e N方法预测带小折角平板边界层的转捩,采用数值模拟和流动稳定性方法研究来流马赫数为4.5,扩张角和压缩角分别为1°、2°、3°的高超声速边界层中扰动的演化,得到不同角度下不同频率的扰动波沿流向的幅值及增长率分布。并对工程中常用的预测带小折角平板转捩位置的e N方法进行验证分析,指出采用e N方法预测带小折角平板转捩时需要对N值进行修正的概念。     

高超声速流场对自由流扰动波响应及边界层扰动波演化

采用高阶精度有限差分方法直接数值模拟了来流慢声波作用下的钝锥高超声速绕流瞬态流场,分析了自由流扰动波与流场的相互干扰,并应用Fourier频谱分析了边界层扰动波时域和空域演化.结果显示:弓形激波出现连续的“∽”变形,初始扰动被显著放大,边界层内外扰动模态存在明显差别.从空域(沿流线)演化来看,在球头附近低频扰动为主导模态,出了球头区,总扰动模态中的低频和高频成分比例迅速转变,高频模态成分显著地增大.至流场下游,大部分低频分量衰减,或者不再增长,边界层内仅存在特殊频率的不稳定波快速增长.从时域演化来看,比起其他模态,主导模态的发展对上游激励的依赖更大.无论时域还是空域演化,都存在模态竞争现象.      

基于PSE的单股剪切混合流稳定性分析

将抛物化稳定性方程(PSE)方法应用到可压缩单股剪切混合流的稳定性研究中.采用并发展了适用于自由剪切流的高精度数值方法,包括六阶紧致格式、坐标变换以及渐近边界条件等,对PSE进行有效求解.通过求解相似边界层方程得到更准确的剪切层内基本流;求解线性稳定性理论(LST)控制方程得到扰动的初始条件,并通过流向空间推进方法对扰...     

超声速平面剪切层多模态混合扰动数值分析

对Ma_c=1.2二维空间发展超声速平面自由剪切层进行扰动模态及流动结构数值分析.使用线性稳定理论和扰动直接数值模拟方法对超声速剪切层中不同流向站位的流动剖面进行扰动参数分析和局部扰动直接数值模拟,获得了剪切层不同站位处的扰动特性.通过特定单频扰动下的直接数值模拟,在剪切层的不同区域内,无相差地促发了慢模态、混合模态和快模态等多种声涡模态,得到了特定单频扰动促发的多模态无相差不稳定波.      

N-S方程数值模拟方法研究及边界层中相干结构的非线性演化

建立了一套三维不可压流动的直接数值模拟方法。该算法在x及y向构造了基于非等间距网格的紧致有限差分格式和非线性项的迎风紧致型格式,在z向采用Fourier谱方法。并研究、提出了三维、非定常流体运动下游边界的无反射出流条件。另外,本文还构建了多涡结构的初始扰动理论模型,利用直接数值模拟的方法研究了边界层中多涡结构的非线性演化特点。结果表明多涡结构非线性演化中的许多特点,如雷诺应力分布、高剪切层的形成、喷射和扫掠等与壁湍流中相干结构的发展规律和现象十分相似。同时,通过对边界层中多涡结构的非线性演化问题的直接数值模拟,验证了该算法具有计算精度高,稳定性好,收敛速度快,出流边界影响小等优点。     

剪切层与边界层组合流动的线性不稳定性分析

对不可压缩剪切层与边界层的组合流动完成了线性稳定性研究.组合流动的数学模型为Blasius边界层相似解与双曲正切函数的叠加,采用整体数值方法求解组合流动的稳定性方程,并验证了程序的准确性及网格无关性.研究给出组合流动的不稳定模态的辨识,即边界层模态和剪切层模态.在此基础上研究了剪切层对边界层模态不稳定性的影响以及壁面对剪切层模态的影响.由于剪切层的存在,使边界层模态中性曲线向左上方平移,临界雷诺数减小.此外,边界层模态不稳定性得到增强或抑制的影响,取决于扰动频率以及剪切层速度比的变化.组合流动中壁面边界层促使剪切层不稳定性减弱,主要表现在低频区域;而在高频区域,剪切层不稳定性几乎不受壁面边界层的影响.      

高超声速钝楔边界层转捩大涡模拟

 以五阶迎风和八阶对称格式混合差分格式求解三维可压缩滤波Navier-Stokes方程,对来流Mach数为6.0、半锥角5°的高超声速空间发展钝楔边界层转捩至完全湍流进行了大涡模拟。时间推进采用紧致存储三阶Runge-Kutta方法,亚格子尺度模型为Driest因子修正的Smagorinsky涡黏性模型。通过定常流场入口边界附近吹／吸引入不稳定扰动斜波的方法数值模拟得到了层流失稳转捩直至完全湍流的空间发展全过程。对扰动的线性、非线性增长以及湍流斑的形成和发展进行了分析,给出了转捩及完全湍流下的速度相关量统计并与实验、DNS结果进行了对比分析,计算结果与理论及实验吻合。     

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

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

局部凹槽对高超声速边界层中第二模态扰动演化的影响

采用直接数值模拟方法研究了局部矩形凹槽对来流马赫数为6.0的高超声速平板边界层中第二模态扰动演化的影响。引入了透射系数(定义为凹槽下游与上游扰动渐近幅值的比率)来量化凹槽的局部散射效应。数值结果表明:对于较浅的凹槽,频率较低的第二模态扰动被促进,而高频扰动的规律相反。对于大多数情况,凹槽深度的增加使得透射系数减小,这表明凹槽对扰动的促进作用减弱或抑制作用增强。当凹槽的深度超过某一临界值时,透射系数与凹槽深度的关系变为正相关,这表明了另一机制的出现。临界深度受扰动的频率影响:频率越低,临界深度越大。