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

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

高精度数值方法与非线性扰动研究

通过数值模拟扰动波在平板边界层中的传播,考察了格式对T-S波传播过程计算的精确性。本文还研究了波-波与涡-涡两种典型的非线性扰动模型,结果表明:在自由剪切层中,波-波模型能加快扰动的增长速度,并激发出了λ结构,但仍然经过了线性失稳阶段;在钝体头部的超声速边界层中,涡-涡模型通过涡-涡相互作用,产生“发卡”结构,进而出现湍流斑。在此过程中,没有线性失稳的迹象。     

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

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

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

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

超声速边界层内小激波的数值模拟

本文通过有限差分计算求解了NS方程,对超声速二维平板边界层的非线性稳定性问题进行了模拟.流场的来流马赫数为4.5.平板从尖前缘开始,流场在前缘处产生一道前缘激波.数值计算模拟了边界层对成倍频关系的双频扰动波的吸收问题.通过计算发现,在一定频率下边界层内可出现小激波.计算给出了小激波的发展过程.     

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

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

钝体头部边界层“逾越”型转捩机理研究

利用数值模拟的方法,研究了再入条件下,钝头体边界层中一种新的失稳机制.这种失稳机制无法用线性稳定性理论解释.最近的研究工作表明,流动中存在着一种与局部有限幅值扰动相关的失稳机制.它与常规的T-S波或二次失稳的过程不同.数值模拟了受空间模式的局部扰动影响的钝体头部超声速边界层的失稳过程.在模拟中,明确了不同扰动模型非线性相互作用的物理机制.      

非线性扰动方程研究平面混合层发声机制

在流动过程中由脉动引起的噪声问题在自然界和工业界中广泛存在,研究流动过程中不稳定波的发声机制对于理解、预测并最终控制气动噪声有着重要意义。本文以二维 Blasius 相似性解作为基本流场,以超声速混合层作为研究对象,采用非线性扰动方程(NLDE)研究超声速混合层不稳定波的近场动力学特性与远场声辐射之间的内在联系。针对不同类型的对流马赫数,分别研究对流马赫数 Mc =0.5(<1)和 Mc =1.2(≥1)两种情况下的扰动发声机制,结果表明：当对流马赫数 Mc <1时,流场中的发声机制主要由大尺度结构的涡条发声,且基本波与亚谐波之间的非线性作用能增强辐射强度;当对流马赫数 Mc ≥1时,根据扰动相速度是相对于上层还是下层自由流速度为超声速,可以进一步分为快慢两种模态,分别对快慢两种模态以及其相互作用模态进行了数值模拟研究,计算得到流场中的发声机制是以马赫波形式辐射向远场,即马赫波辐射,其辐角的计算结果与理论值相符。     

用弱非线性理论研究二次稳定性问题

本文把研究三波共振的弱非线性理论应用于平板边界层和平面Poiseuille流三维扰动的二次稳定性问题。其数值结果与文献[1]、[2]用Floquet形式的二次稳定性理论得出的结论是一致的。     

等熵侧压诱导的超声速平板边界层二次流研究

为了揭示等熵侧压诱导的超声速平板边界层二次流的形成机理,基于建立的简化流场模型,结合小扰动分析和数值计算,研究了来流边界层马赫数不均匀分布因素对二次流形成的贡献。研究表明,侧压诱导的二次流是流场的平均行为,和瞬时的扰动与涡的相互作用无关。压缩波在主流和边界层中的传播特征不同,导致在侧壁压缩作用下主流和边界层流体压缩不均匀,从而产生了边界层流体相对于主流的横向运动;密度梯度和压力梯度矢量方向不一致,即流动的斜压性是流向涡形成和发展的一个重要因素。     

On Nonlinear Evolution of C-type Instability in Nonparallel Boundary Layers

The process of evolution, especially that of nonlinear evolution, of C-type instability of laminar-turbulent flow transition in nonparallel boundary layers are studied by means of a newly developed method called parabolic stability equations （PSE）. Initial conditions, which are very important for the nonlinear problem, are investigated by computing initial solution of the harmonic waves, modifying the mean-flow-distortion, and giving initial value of TS wave and its subharmonic waves at initial station by solving linear PSE. A numerical method with high-order accuracy are developed in the text, the key normalization conditions in the PSE are satisfied, and nonlinear PSE are solved efficiently and implemented stably by the spatial marching. It has been shown that the computed process of nonlinear evolution of C-type instability in Blasius boundary layer is in good agreement with the experimental results.     

不可压平板边界层转捩机理

基于扰动形式N-S方程,从空间模式的角度,采用Fourier伪谱及MPI(massage passing interface)并行方法,模拟了不可压平板边界层从层流到湍流的转捩过程.通过对计算统计数据的分析,比较了不同幅值入口扰动引起的转捩过程的同异.研究结果表明:当层流中扰动幅值逐渐增大后,非线性作用将修正平均流剖面,表现为不稳定区域逐渐扩大,很多高次谐波被激发.当平均流剖面被修正到一定程度时,不稳定区域变得很大,这使得更多的高次谐波被快速激发并迅速增长,即转捩过程开始.由于大量谐波快速增长导致扰动能量快速增长,同时,扰动能量的快速增长又进一步加速了平均流剖面的快速修正,即平均流剖面的不稳定特性发生了改变,这样的相互作用使得层流快速变为湍流,因此,平均流剖面不稳定特性的改变在转捩过程中起到了关键作用.      

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

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

导电圆柱壳体磁弹性二阶谐波共振响应分析

 研究了磁场环境中受机械载荷作用圆柱壳体的二阶谐波共振响应问题。首先给出了电磁场环境中导电圆柱壳体的非线性磁弹性振动方程和电磁场方程,并应用伽辽金积分法导出了相应的无量纲化非线性振动微分方程。然后采用多尺度法对系统的二阶超谐波和亚谐波共振问题进行求解,得到了稳态运动下关于定常解的幅值响应方程及解的稳定性判别式。最后通过数值算例,给出了共振幅值随磁感应强度和外激励力等参量变化的特征曲线图,分析了解的稳定性、奇点性态及分岔特性。结果表明,通过电磁或机械参数的适当选取,能够达到激发或抑制系统共振现象的目的。     

高超声速边界层转捩机理及应用的若干进展回顾

高超声速边界层转捩对飞行器的热传递、表面摩阻和流动分离等有重要影响,尤其是再入飞行器和吸气式巡航飞行器。然而,人们对边界层转捩机理中的很多问题认识还不清楚,或存在争议。本文从扰动波演化的角度回顾了高超声速边界层感受性、线性稳定性和非线性作用的国内若干研究进展,并以基于谐波共振的人工转捩技术为例示范了这些机理认识在转捩控制上的应用。扰动的产生和发展是认识边界层转捩机理的核心。通过研究扰动波来认识边界层转捩机理,开展应用创新研究对提升飞行器性能具有重要意义。     

共轴异向旋转筒间不可压流的临界主不稳定性

基于长隙比无穷大和Tollmien-Schilchting 行谐波假设，通过导出坐标系下不可压一阶改型边疆特征值问题，实现了异向旋转筒流动临界主不稳定性的线性分析，在稳定性方程离散中引入了一类高精度对称紧致差分格式，在特征值搜索中发展了一个高效率复矩阵双重反迭代算法。双临界点分析揭示了几何参数、、流动参数、流动参数、扰动参数相互间的变化规律。理论结果同实验测量数据基本一致。     

后掠翼可压缩三维边界层稳定性计算

采用可压缩线化稳定性理论对后掠机翼的层流边界层稳定性进行了计算研究。对三维边界层中三维扰动波的可压缩稳定性方程组特征问题采用相关两步格式求解。与常用的初值法相比,具有不需预估特征值和快速收敛的特点。为得到渐近边界条件的准确表达式,导出了外边界上方程的解析解,这对精确的特征值计算和分析至关重要。为避免出现伪不稳定模态,选择了合适的稳定性方程组形式。最后,通过层流控制机翼和自然层流机翼的算例,研究了后掠翼上典型Tollmien-Schlichting(T-S)型扰动和横流(C-F)扰动稳定性问题,其结果与已有计算结果符合甚好。