可压缩剪切层与边界层流动相互作用的稳定性研究

为了探究可压缩剪切层与边界层组合流动稳定性规律,采用线性稳定性理论对这种组合流动的线性稳定性问题进行研究。平均流动模型取双曲正切函数描述的可压缩剪切层与可压缩相似边界层速度剖面的直接叠加。采用整体数值方法 (Global)求解组合流动线性稳定性方程,获得全部不稳定模态。结果表明,在这样的组合流动中,存在两个不稳定模态,一个与剪切层相关的不稳定模态,另一个是与边界层相关的不稳定模态。在此基础上,研究了剪切层与壁面之间的距离、扰动频率、雷诺数、马赫数、速度梯度、温度梯度对这种组合流动稳定性的影响。马赫数增大对不稳定模态起抑制作用,剪切层对边界层的不稳定模态的促进、抑制作用与马赫数、扰动频率有关,边界层对剪切层不稳定模态基本上起抑制作用。     

自由剪切层中的三维不稳定性

本文是在文[1]的基础上研究自由剪切层中由Kelvin-Helmholtz不稳定波发展而形成的展向大涡结构的三维不稳定性。以大涡结构为基本流动,将稳定性分析归结为二维特征值问题,用pseudo spectral(伪谱)方法数值求解。研究发现:在没有亚谐波存在的情况下,大涡结构的最不稳定的扰动波是流向波长和其相同。展向波数较高,有对流特性的三维扰动波。它在剪切层中的发展与展向涡量的分布有关,大涡结构的涡核不稳定性和辫子不稳定性是流向涡形成的主要力学机制。本文还给出了不同雷诺数下三维扰动波增长率与展向波数的关系,这些结论与实验及数值模拟结果基本一致。     

超声速自由剪切层流动的数值模拟和理论分析

本文从二维完全Ｎ－Ｓ方程出发，采用ＮＮＤ２Ｍ差分格式，计算了超声速自由剪切层流场，并且研究了上游扰动在剪切层内的传播规律。在数值模拟的基础上对流动涡结构进行了分析。     

二维超声速剪切层的非线性失稳过程分析

利用直接数值模拟方法，模拟了二维超声速剪切流动中的小扰动波随时间发展在空间上的传播过程，由此分析其稳定性，通过相图分析，李雅普诺夫指数分析和分维分析，本文得到了二维超声速自由剪切层流动由线性到非线性失稳过程的数学特征，结果分析表明，在本文的计算条件下，二维超声速混合层流动在失稳过程中出现了间歇现象，并进入了混沌运动的初期阶段。     

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

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

空气热力学参数对可压缩边界层线性流动稳定性的影响

本文详细讨论了空气热力学参数的不同近似方法对可压缩平板边界层线性稳定性的影响。数值结果表明Ｍ＝３以上的线性稳定性问题必须采用适当的热力学参数近似关系。     

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

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

高速三维边界层稳定性问题的数值方法

本文把四阶精度紧致差分格式用于求解声速、高超声速三维边界层稳定性问题。提出了适用地求解高速边界层特征值问题的Ｃｈｅｂｙｓｈｅｖ多项式配置点方法,并将各种数值方法进行了比较。数值结果表明上述两种方法是求解高速三维边界层稳定性问题的有效数值方法。     

干扰剪切流动稳定性理论及其对高雷诺数流动数值模拟方法的改进

在干扰剪切流(Interacting Shear Flow,ISF)理论的基础上,提出ISF稳定性理论并把它用于改进高雷诺(Re)数流动计算方法。(1)高Re数内外绕流的RANS计算及工业标准PNS计算中,流动转捩的预测均基于经典边界层理论;然而转捩并非总是最早发生在边界层中,例如发生在壁面小突起、小凹坑、小窄缝等局部粘性/无粘强干扰区,这些强干扰区可能位于边界层内,但边界层理论并不适用于它们,又如转捩发生在分离点邻域强干扰区等。(2)ISF理论表明:高Re数内外绕流为一复杂ISF,转捩总是最早发生在该ISF的层流区中。(3)ISF稳定性理论表明:作者提出的干扰剪切扰动流(Interacting Shear Perturbed Flow,ISPF)方程组可以计算ISF层流中非湍流扰动运动演化并预测转捩;ISF方程组和ISPF方程组分别与PNS和抛物化稳定性方程(PSE)为同类方程组,PSE分析计算边界层稳定性的众多成功实践,说明用ISPF(即PSE)方程组计算ISF层流扰动流并预测转捩完全可行。(4)RANS和PNS方法经ISF稳定性理论改进后,在转捩前用ISF方程组(即PNS)计算ISF层流基本流,用ISPF方程组(即PSE)计算ISF层流扰动流并预测转捩位置;转捩后RANS方法计算RANS或RANS/LES,PNS方法计算干扰剪切湍流(ISTF)方程组即抛物化RANS(PRANS)方程组。改进后的两方法,理论合理正确,方程体系完备、自洽,ISF方程组只能用ISPF方程组相配对,因此是高Re数内外绕流计算的理想且可持续发展的两种方法。     

考虑流向曲率和压力梯度的可压缩边界层稳定性分析

以圆弧,平板和尖锥为研究对象,基于线性稳定性理论,重点讨论了曲率和压力梯度对可压缩边界层稳定性的影响。     

Linear stability analysis of interactions between mixing layer and boundary layer flows

The linear instabilities of incompressible confluent mixing layer and boundary layer were analyzed.The mixing layers include wake,shear layer and their combination.The mean velocity profile of confluent flow is taken as a superposition of a hyperbolic and exponential function to model a mixing layer and the Blasius similarity solution for a flat plate boundary layer.The stability equation of confluent flow was solved by using the global numerical method.The unstable modes associated with both the mixing and boundary layers were identified.They are the boundary layer mode,mixing layer mode 1 (nearly symmetrical mode) and mode 2 (nearly anti-symmetrical mode).The interactions between the mixing layer stability and the boundary layer stability were examined.As the mixing layer approaches the boundary layer,the neutral curves of the boundary layer mode move to the upper left,the resulting critical Reynolds number decreases,and the growth rate of the most unstable mode increases.The wall tends to stabilize the mixing layer modes at low frequency.In addition,the mode switching behavior of the relative level of the spatial growth rate between the mixing layer mode 1 and mode 2 with the velocity ratio is found to occur at low frequency.     

基于边界层相似性解的放大因子输运模型

为了使二维线性稳定性理论能够适应现代CFD求解技术,通过求解Falkner-Skan边界层相似性方程获得各个形状因子下的相似速度型。基于线性稳定性理论,对各个速度型进行稳定性分析从而获得对应不同速度型的扰动放大因子包络线。最后使用标量输运方程的形式实现包络近似方法中放大因子的当地化求解,并结合原始γ-Re_θt转捩模型中的间歇因子输运方程,实现了自然转捩和分离泡转捩的建模。使用该输运模型对S&K平板、S809翼型、NLR7301翼型和DLR-F5机翼进行转捩预测,结果均与试验结果吻合较好,验证了该模型构建的合理性和可行性。     

尾缘锯齿结构对叶片边界层不稳定噪声的影响

实验研究了不同雷诺数（2×10⁵~8×10⁵）、不同攻角状态下,3种相同波长（4%弦长）不同振幅（分别为5%、10%、15%弦长）尾缘锯齿结构对叶片层流边界层不稳定噪声的影响。研究表明,在0°攻角状态下,尾缘锯齿会增强甚至诱导产生新的不稳定噪声,显著增大叶片自噪声;在大攻角状态下,尾缘锯齿会减弱甚至完全抑制不稳定噪声,降噪量高达40 dB,降噪机制在于尾缘锯齿结构破坏了不稳定噪声产生所需的声学反馈回路。尾缘锯齿会降低不稳定噪声频率,且锯齿振幅越大,不稳定噪声频率越低。     

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

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

横向矩形微槽对高超边界层失稳的控制作用

模拟了马赫数为6的空间发展平板边界层,通过在平板表面添加二维横向微槽研究了其对基本流及第二模态扰动波的影响。结果表明:所研究的二维微槽构形与常规多孔涂层相比具有较大的尺寸(100 μm以上),微槽对基本流的流向速度影响很小,空间采样点上的频谱并未出现新的不稳定模态。微槽有利于减小摩擦阻力,槽内流动导致的压差阻力比摩擦阻力低1~2个数量级,总阻力随着开槽率的增加而减小,当开槽率为06时,减阻效果超过40%。在入口以两种方式添加了第二模态扰动波,包括单个扰动和多个叠加扰动,下游的演化结果说明横向微槽能够在一个宽带频率范围内对第二模态的增长起到明显的抑制作用,且控制效果随着开槽率的增大而增强。     

受限混合层的线性稳定性

针对超声速边界层/混合层组合流动,利用可压缩线性稳定性理论研究了流动的线性失稳特性。基本流场选取了具有不同速度特征的两股来流,采用双曲正切的混合层剖面叠加可压缩边界层自相似性解剖面构造。重点考察了混合层中心与壁面距离、对流马赫数等参数对组合流动稳定性特征的影响,其中壁面采用绝热壁面。混合层中心与壁面的距离为5~15倍的边界层厚度,混合层的对流马赫数为0.6~1.2。结果表明:该组合流动中存在独特的多重不稳定模态,并相互影响;且其不稳定模态随着壁面距离及对流马赫数的变化呈现出不同的主导行为。      

小攻角高超声速钝锥边界层失稳特性

利用高阶紧致格式,采用直接数值模拟(DNS)和线性稳定性(LST)分析方法,对高超声速边界层的失稳机制和转捩特点进行了研究.通过对马赫数为6的2?攻角高超声速钝锥边界层的稳定性分析发现:小攻角高超声速钝锥边界层存在多枝不稳定模态;周向速度使钝锥的稳定性特征与不考虑周向速度时有本质的差别;转捩线在接近背风面处出现拐折现象是由失稳模态发生转换引起的.      

低雷诺数下跨声速转子流动失稳及附面层抽吸扩稳研究

采用数值方法模拟了低雷诺数条件下NASA Rotor 37跨声速压气机转子内部流场.结果表明,附面层径向涡是该压气机转子流动失稳的一个很重要的原因.通过在该压气机转子叶片吸力面上不同的叶高和轴向弦长处进行抽吸,发现都可以很好的抑制附面层径向涡的发展,压气机转子的稳定工作范围明显扩大.此外还比较了不同的抽吸量对压气机的性能的影响.