钝体边界层失稳特性分析

推导了适合可压缩二维/轴对称边界层的线性抛物化稳定性方程。在法向采用4阶中心格式,流向采用1阶向后差分,对线性抛物化稳定性方程进行离散。通过来流马赫数为4.5的平板与曲面边界层稳定性问题验证了程序的正确性。讨论了来流马赫数6的球钝锥边界层内,扰动波沿端头向下游演化的失稳特性。发现了一种新的失稳机制。     

非平行边界层C型失稳的非线性演化(英文)

提出用抛物化稳定性方程(PSE)新方法,研究非平行边界层从层流向湍流转捩的C型失稳的演化过程,特别是非线性演化问题。研究了至关重要的初始条件,包括各扰动谐波的初始解、初始站均流变形的模拟、以及通过对线性 PSE 的求解获得 TS波和C型亚谐扰动的初始站值等。发展了高精度的数值方法,满足PSE的关键性正规化条件,实现非线性抛物化稳定性方程的有效求解和稳定推进。算例展示了Blasius边界层的C型失稳的非线性演化过程,与典型实验结果相吻合。     

超声速非平行边界层的非线性稳定性研究(英文)

用非线性抛物化稳定性方程(NPSE),研究超声速非平行边界层的扰动波的非线性演化过程。发展了高精度的数值方法,对NPSE进行有效求解和稳定推进。计算得到了各谐波的振幅增长曲线和不同的扰动剖面函数,特别是展示了谐波之间的非线性作用的物理特征,包括流向涡、展向涡和Λ结构的形成和演化等,以进行转捩过程的机理研究和分析。NPSE的计算结果显示了这个新方法是研究超声速边界层非线性稳定性的有效方法。     

圆环管道中曲率对流动稳定性的影响

为了研究圆环管道中曲率对流动稳定性的影响,选取了展向带有曲率的轴向流动和轴向带有曲率的周向流动两种圆环管道模型进行分析.采用线性稳定性理论分别计算了两种模型下不同曲率的稳定性和扰动演化情况,并和槽道流进行比较,从而得到了曲率对基本流动、特征值和特征函数的影响.结果表明：曲率越大,基本流向内侧偏转越大,扰动波流向波数越大,增长率越小,中性曲线的不稳定区域越小,特征函数的对称性越差.因此曲率对边界层内的扰动起着稳定的作用,并且在一定范围内曲率越大,稳定性越好.     

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

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

边界层稳定性的三维线性抛物化算法研究

通过把扰动波在流向上的变化分解成缓慢变化的形状函数部分和快速变化的波动部分,同时把扰动波在展向上的分布也分解成形状函数部分和波动部分,得到三维抛物化稳定性方程(3D PSE).然后在流向采用一阶向后差分,法向和展向采用4阶中心差分,发展出了迭代求解三维线性PSE的方法.通过求解二维扰动波和三维斜波在三维空间中的发展,发现二维波可以自动调制出展向扰动,且展向扰动的最大值点在边界层与外流的交接面附近.三维斜波的计算结果表明,网格密度应该与波矢方向上的波长相匹配,而不是与流向波长相匹配.     

小扰动在二维超声速边界层中的弱非线性传播

采用非线性抛物化稳定性方程(PSE)研究了非平行边界层的弱非线性稳定性。通过在流向采用一阶向后差分,法向采用四阶中心差分,并用预估 校正迭代来耦合非线性项,发展出了求解非线性PSE的数值方法。研究了不同幅值扰动及其高倍谐波的演化情况,发现随着基频扰动幅值的增大,高倍谐波可能失稳,而且失稳的位置会随着基频扰动的增大而向上游移动。通过观察涡量的发展可以发现涡破碎现象。算例与Navier-Stokes方程的直接数值模拟（DNS）结果作了比较,初步检验了其正确性。     

机翼边界层的横流稳定性分析和转捩预测

通过求解经典O-S方程(LST)、扩展的O-S方程(EOS)和线性抛物化稳定性方程(LPSE),对展向无限长、后掠角25°、迎角0°、来流Mach数0.8、单位Reynolds数6.79×106/m的机翼边界层进行了稳定性分析,结合eN方法进行了转捩预测。研究发现无限长后掠机翼在(x,ω)平面上的中性曲线没有下支,在(x,β)平面上的中性曲线呈反拇指的形状,横流不稳定性在机翼前缘占主导作用。当外界扰动进入边界层后,幅值将被直接放大,对于频率相同的扰动,首先是展向波数β大的增长起来,演化到一定位置开始衰减,然后是展向波数β小的逐渐增长起来,并且增长的指数N逐步超过波数大的扰动。转捩在机翼前缘完成,引起转捩的扰动波的展向波长约为2mm。     

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

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

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

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

边界层流动的非平行稳定性研究

采用极为有效的抛物化稳定性方程（PSE）方法研究边界民支的非平行稳定性。并利用抛物线坐标变换、有限远延拓边界条件、及流向和法向全差分的数值方法，对不同频率的二维Tollmien-Schichting波的非平行稳定尾进行了计算和分析，计算结果与Orr-Sommerfeld方程（OES）的解以及谱配置方法的PSE结果作了详细的比较，得到了满意的结果。     

基于APSE确定高超声速边界层第二模态中性曲线下支

在高超声速边界层流动中,由于存在非平行效应以及激发Mack第二模态时的模态同步机制,因此在进行转捩预测研究时,中性曲线下支的确定需要考虑到这两方面的影响.全文针对高超声速平板边界层,首先采用线性抛物化稳定性方程(LPSE)对考虑感受性的中性曲线下支进行了求解,研究发现该方法存在两大难点:一是计算入口需要给定的离散模态在...     

局部扰动对平板边界层流动稳定性影响的研究

在壁面上构建局部扰动,数值研究不同特征的局部扰动对平板边界层流动稳定性影响的作用机理,进一步探讨壁面扰动的分布、类型、强度大小对三维扰动波在平板边界层流动中的非线性演化规律以及失稳机制的问题.确定什么样的局部扰动能使三维扰动波快速增长或者是什么样的局部扰动能对流动起抑制的作用,同时寻求在三维扰动波向下游发展过程中的传播方向会发生怎样的变化等.对该问题的深入研究,将使人们能深入理解认识平板边界层流动过程中的转捩发生、流动失稳及湍流形成的理论机制.