边界层稳定性的非平行性问题研究

采用抛物化稳定性方程的新方法研究稳定性的非平行问题，其控制方程从涡量形式的N－S方程导出，为提高计算精度和收敛速度，在法向采用高精度的高阶谱方法，结合边界层特征安排边界条件的方程，以及对无限区域数值问题作高效的代数变换处理，通过对预估校正迭代过程中推进步长有效的控制，以满足正规化条件，研究了不同类型压力梯度对稳定性的作用，由系列算例得到的增长率变化曲线和扰动中性曲线等，精确地给出了非平行性对流动稳定性的影响，计算的结果与相关数据符合甚好。     

二维抛物化稳定性方程计算

针对二维抛物化稳定性方程(PSE)在法向上采用正交函数展开法进行了计算.在求解PSE的初始条件以及L0的中性曲线中应用了Chebyshev多项式展开.计算结果表明可以获得非常精确的解,证实了在PSE中应用Chebyshev多项式的有效性.     

高速流非平行边界层稳定性研究

导出了可压缩流的抛物化稳定性方程（PSE）。针对高速流动,特别是超声速和高超声速流动的非平行边界层稳定性问题进行了研究。引入高效的边界层变换、全流场高精度的差分格式及预估校正迭代和推进求解法来求解PSE方程,使得PSE方法中至关重要的正规化条件得到了满足,确保了数值计算的稳定。采用高马赫数下对稳定性起支配作用的第二模式,研究了高速流边界层稳定性的演变和特征,分析了流动的非平行性、压缩性,以及壁面冷却等因素对流动稳定性的影响．所得结果与相关实验数据吻合较好。     

边界层中TS波及其高阶谐波的演化

基于抛物化稳定性方程，研究了边界层中TS波及其高阶谐波的线性和非线性演化问题。由局部法和Landau展开式导出初始条件，并计算了扰动幅值和速度型等的演化过程和特征，特别是非线性的重要作用。探讨了初始幅值、压力梯度、扰动频率对扰动演化的影响及其规律，这与边界层的稳定性和转捩研究紧密相关。算例结果与全Navier—Stokes方程的直接数值模拟结果一致。     

边界层壁面振动诱导T-S波的直接数值模拟

采用空间模式的平板边界层Blasius解为基本流,利用直接数值模拟方法求解二维不可压扰动方程．研究了边界层对二维壁面局部持续微振动的感受性。计算结果表明,当地扰动速度随时间变化显示弦振动特性,振动周期与壁面加载周期完全一致;扰动速度在空间上具有波动解的特性,流向波数与O一S方程解接近;数值计算获得的扰动幅值增长率稍大于预测值,二者运行趋势基本吻合,且对所计算的雷诺数条件,周期为30的扰动幅值逐渐增长,而周期为20的扰动帽值先增加后衰减;周期为30的扰动的二次谐波明显大于周期为20的扰动;流向扰动速度及法向扰动速度的剖面模值函数在壁面局部振动下游附近与O-S方程解不完全重合,而在壁面局部振动下游较远处则与其重合良好;边界层壁面局部微振动获得的二维扰动解具有T—S波的形式。     

具有参数激励约简的扰动KdV方程的混沌行为

研究了具有参数激励约简的扰动KdV方程的混沌动力学行为.利用改进的Melnikov方法分析了由于同宿轨道的横截相交而产生的混沌行为.对周期外激励、周期线性参数激励和周期非线性参数激励下的扰动KdV方程的混沌行为进行了比较,发现划分混沌区与非混沌区的临界曲线是互不相同的.尤其是对非线性参数激励系统,存在"死频率".当这类系统受到该频率激励时,不论激励的振幅多大,混沌也不会发生.用时间积分法对上述系统进行了数值计算,结果与理论分析一致.     

带压力梯度的非平行流边界层稳定性研究

发展高精度的高阶谱方法，通过抛物化稳定性方程及其局部法，研究在压力梯度作用下的非平行流边界层稳定性问题。与经典的平行流边界层稳定性结果相比，显示了在某些条件下非平行性对稳定性的关键作用。提出了压力梯度对临界雷诺数的影响研究，探讨了不同压力场对流动稳定性作用的规律性。     

软管-锥套式空中加油动态建模与性能分析

提出了一种基于凯恩方程的自主空中加油软管-锥套动态模型及性能分析方法。该方法将软管视为由有限段以铰链形式连接的刚性 杆组成，锥套为软管末端的一个质点。本文定义了描述系统状态的广义坐标及广义速率，导出了软管段位置多级递推公式和系统动态方程，估算了软管在加油机尾流、定常流和大气扰动下的气动载荷。通过数值仿真分析了平稳大气中锥套在加油机不同飞行条件下的稳态阻力和软管拖拽轨迹，验证了模型的正确性及系统稳定性。最后研究了大气扰动对锥套运动的影响及不同软管段的受扰运动。     

均匀受压粘弹性圆环稳定性分析

采用Boltzmann线粘弹性元件本构理论,运用假设模态法、Galerkin方法及数值计算相结合的方法,建立了均匀受压粘弹性圆环的屈曲模型,研究了粘弹性系统在外在载荷的作用下的极限载荷及稳定性位移问题。发现了粘弹性结构蠕变屈曲的典型特征：延迟失稳。通过极限分析法和数值方法研究了环失稳后的后屈曲模态及相关参数对系统后屈曲解的影响。     

高阶(2+1)维色散长波方程的孤波解

为了求高阶(2 1)维色散长波方程的精确解，先通过变换将其变为一个简单的方程．然后利用各种方法求该方程的精确解，最后再通过变换得到了高阶(2 1)维色散长波方程的若干新孤波解，本文简化了求非线性模型精确解的双曲正切函数法、双曲函数法、平衡齐次法及扩展的双曲正切函数法等方法。     

三维扰动波在局部粗糙边界层中的稳定性研究

采用局部喷、吸、喷和吸组合来模拟局部粗糙壁面,数值计算获得稳定的三维边界层的基本流.在此基础上研究了三维扰动波在该基本流中的空间演化问题,讨论了局部粗糙的形式、分布结构对三维扰动波的幅值增长率及流动稳定性影响.计算结果表明:三维局部粗糙对三维扰动波的增长、涡的形成都起着激励的作用.扰动波演化产生的平均流修正及局部粗糙诱导展向速度的存在,影响着流体运动稳定性.与光滑壁面相比,三维局部粗糙作用下扰动波的传播角度与相位角发生明显变化,而不同形式的二维局部粗糙壁面边界层显示出不同的稳定特性.     

局部粗糙边界层流中二维T-S波的非线性演化

采用局部喷、吸、喷与吸组合的结构来模拟局部粗糙壁面,通过数值计算获得了稳定的三维边界层基本流的数值解.在此基础上,研究了最不稳定的和稳定的二维扰动T-S波在时、空间演化过程中非线性问题;探讨了局部粗糙的形式、强度大小及分布结构对二维T-S波的非线性演化的影响及其对增长速率贡献的大小.数值结果表明,三维局部粗糙壁面诱导产生的基本流是扰动波得以增长的关键因素;最不稳定的二维扰动T-S波在非线性演化过程中,由于非线性相互作用的不断增强,逐渐诱导激发出三维扰动波及高次谐波,其三维扰动波的流向波数和频率与二维扰动波的流向波数和频率大致相同.此外,展向速度的大小激励着二维扰动波的快速增长,随后发生初次失稳、流向涡形成、正负相间剪切层逐渐增强、诱导产生三维不稳定扰动波及高阶失稳等机制.上述结论与实验结果相吻合.     

前轮非线性摆振稳定性分析

飞机前起落架支柱上的几个非线性项,例如扭转间隙、库仑摩擦以及速度平方阻尼等,对前轮摆振稳定区域有着重要的影响.文中假设前起落架支柱上端固支,并且给出了用以描述五自由度的前轮摆振运动的非线性微分方程组,使用描述函数法研究了上述三个非线性项对前轮摆振稳定区域的影响,将准线性化的方程组代入状态空间内,利用优化特征值实部的方法确定摆振系统的极限环幅值、摆振频率以及临界参数曲线,这种方法可以免去以往摆振稳定性分析时冗长的代数推导.作为比较,文中还计算了前轮线性摆振系统的临界频率以及临界参数曲线,两种结果比较分析表明,支柱内的几个非线性项对前轮摆振稳定区域有着重要的影响,由于扭转间隙的存在,使得飞机临界滑跑速度大大降低.