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

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

窄带信号小样本谱估计算法

针对窄带平稳随机信号的小样本谱估计问题,结合均匀线性阵列的空间谱估计算法,提出了一种高性能小样本谱估计算法。该算法首先利用采样数据构造一个Toeplitz矩阵;然后对该矩阵进行特征值分解得到信号子空间和噪声子空间;再结合MUSIC、ESPRIT等子空间算法,即可实现对小样本采样数据的谱估计。数值仿真验证了该算法的有效性。     

二维抛物化稳定性方程的特征分析

在对抛物化稳定性方程(PSE)的基本流场没有做任何近似假定的情况下,分析了PSE的特征性质。分析表明,当法向速度不为零时,PSE有一个非零主特征值,其余主特征值都为零。PSE的次特征值与扰动波的空间波数α有关,α的实部代表扰动波的波动情况,它可以直接导致复特征值出现;α的虚部表示扰动波的增长(衰减)情况,当它的绝对值超过一定范围时,也会在边界层内亚声速区的局部区域导致复特征值出现。增大求解PSE的空间推进步长,可以克服PSE的椭圆性。     

飞机纵向动稳定性的奇异摄动理论

本文应用奇异摄动理论分析飞机的纵向扰动运动。首先利用匹配的渐近展开法将扰动运动近似分解为长、短周期两个模态,并利用Van Dyke匹配原理及加法法则给出了合成解,它给出了在任意的起始扰动条件下整个扰动运动的发展过程。其次,利用坐标变形法对特征值进行了修正。算例表明,经过修正后的长周期特征值和精确解相当吻合。理论分析还表明,短周期特征值的修正项是高阶小量,因此通常零级近似的特征值已经与精确解相差无几了。     

Wilkinson定理扰动界的改进估计

Wilkinson定理是代数特征值问题中的一个经典定理,文章给出了Wilkinson定理的结果中关于扰动矩阵上界的另一种估计形式,并指出该形式对Wilkinson定理具有一定的改进。     

孤立翼尖涡模态演化规律的实验研究

采用主动流动控制方法加快翼尖涡衰减破碎是提升机场起降频率、保证飞机飞行安全的最具潜力的技术之一。由于翼尖涡不稳定性认识的不足,已有的主动控制方法常常不能获得最优的控制效果。为了揭示翼尖涡不稳定模态的演化规律,采用体视粒子图像测速技术和线性稳定性分析方法对孤立翼尖涡的不稳定模态演化特征进行研究,结果表明：孤立翼尖涡的扰动模态可以根据其在特征值谱的位置分成主扰动模态、P族次级扰动模态、A族次级扰动模态、S族次级扰动模态4种;其中主扰动模态和P族扰动模态具有两瓣式的结构特征,决定了翼尖涡摇摆的各向异性特征,A族次级扰动模态具有流向速度波动大于横向速度波动的特征,S族次级扰动模态则具有更高的切向波数和作用范围。不同族扰动模态的流向演化规律不同,翼尖涡的主扰动模态和P族扰动模态沿流向发生旋转,并且扰动幅值随着流向逐渐放大,A族次级扰动模态随着流向发展会逐渐增大扰动幅值;S族次级扰动模态随着流向会逐渐覆盖住整个涡核,这种穿透涡核的扰动会随着流向进一步放大。描述了不同翼尖涡扰动模态的扰动能量随流向的演化规律,发现S族次级扰动模态有更高的切向波数特征,也同时具有较高的扰动能量增长,意味着利用S族次级扰...     

一种基于数值计算的轴流压气机旋转失速稳定性模型——之二:简化模型及应用

提出了一个基于直接数值计算结果基础上的旋转失速线性稳定性模型,描述了如何从基本的Navier-Stokes方程出发,把轴流压气机流动稳定性问题转化为一个特征值问题的基本理论框架,对压气机不稳定边界的判断通过计算和检验特征值的方式来实现,给出了有限差分法和基于Gauss-Lobatto节点的谱配置方法两种数值离散的方式,并分析了这两种离散方式的收敛特性.     

基于特征值理论的轴流跨声速压气机失稳预测

流动失稳问题是叶轮机械所面临的严峻挑战之一,而旋转失速是轴流压气机及压缩系统中普遍存在的一种典型流动不稳定状态。基于特征值理论发展了一个轴流压气机流动失稳开始点预测模型,应用谱方法和奇异值分解(SVD)算法求解所建立的特征值方程,通过特征频率的虚部判断压缩系统的流动稳定性。并通过一个高速单转子在高低两个转速时的算例,定量分析了流动压缩性对跨声速压气机流动失稳开始点预测的重要影响。该模型可在合理的精度范围内,以工程应用可接受的计算资源消耗来预测跨声速单级压气机旋转失速的开始点流量,同时不应用大量经验数据,为压气机设计阶段的稳定性评估提供指导。     

一种基于数值计算的轴流压气机旋转失速稳定性模型——之一:理论模型及算法分析

本文提出了一个基于直接数值计算结果基础上的旋转失速线性稳定性模型,描述了如何从基本的Navier－Stokes方程出发,把轴流压气机流动稳定性问题转化为一个特征值问题的基本理论框架。对压气机不稳定边界的判断通过计算和检验特征值的方式来实现。本文给出了有限差分法和基于Gauss－Lobatto节点的谱配置方法两种数值离散的方式,并分析了这两种离散方式的收敛特性。     

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

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

求解结构动力反问题的逆特征值方法

设计一个结构,使其具有事先给定的频率和（或）振型,这类结构动力反问题有重要的应用和理论价值。本文提出了求解它的逆特征值方法,并给出了数值例子。算例一说明:逆特征值方法的计算精度要比反摄动法高得多。经过实验验证的算例二表明:计算值、目标值和实验值十分吻合。     

大型广义特征值问题的动力降阶迭代法

 本文在Guyan、Kuhar等人工作基础上,提出通过动力降阶和迭代,对大型广义特征值问题实施逼近求解的一个新方法,以能对大型、复杂动力系统实施动力分析。既解决计算机容量不足,又保证精度,亦不耗时过多。通过实际使用,证明该方法的有效性。     

低速预处理方法在悬停旋翼模拟中的应用

 在可压缩方程的求解当中引入时间导数低速预处理技术,构造了能够同时处理低速和高速流场的统一数值方法,采用该方法对Caradonna-Tung旋翼的亚声速悬停状态进行了模拟,其中方程采用在旋转坐标系下求解,网格采用重叠网格进行模拟。数值结果表明,在加入低速预处理方法以后,流场的收敛速度得到了明显提高,更为适合对下洗速度场等细节的模拟;另外,在加入低速预处理以后,方程的各个特征值都处在同一量级,避开了因特征值差异带来的数值耗散过大的问题,提高了数值计算的精度。     

基于广义Richardson外插方法的颤振模拟耦合时间精度研究

为更准确地评估颤振问题时域模拟在实际计算中表现出的时间精度,采用不同精度的耦合方法对Isogai Wing和AGARD 445.6 Wing的颤振问题进行时域模拟。参照网格收敛性分析方法,提出了基于广义Richardson外插方法的颤振问题时域模拟耦合时间精度分析方法,分析时域计算结果的时间步长收敛性,计算数值模拟结果时间精度并获得时间步长无关解。分析表明,对于Isogai Wing和AGARD 445.6 Wing颤振问题,计算结果具有良好的时间收敛性,采用广义Richardson外插方法对各耦合方法分析所得精度与理论分析结果基本一致。在合适的时间步长区间内,可忽略具体时间步长选取对广义Richardson外插方法分析结果的影响,验证了提出的气动弹性模拟耦合时间精度分析方法。     

基于滤波残差的组合导航系统故障检测研究

论文研究了两种用于组合导航系统中的故障检测方法,首先介绍了传统的残差χ^2检验法,随后给出了一种基于最大残差样本协方差阵特征值的故障检测算法。通过计算机仿真比较了两种算法在组合系统中的应用,仿真结果显示,基于特征值的故障检验算法对故障较敏感,系统组合精度略高,一定程度上优于传统残差χ^2算法。最后,文章讨论了组合导航系统中两种故障检测算法的优缺点。     

基于Kriging形函数的线性时变结构模态参数辨识

近年来,对航空航天飞行器随时间变化的动力学特性研究需求越来越迫切。仅输出参数化时域的时变时间序列模型以其结构简约、精度高且跟踪能力强而成为研究热点,尤其是泛函向量时变自回归(FS-VTAR)模型已经得到了广泛应用。然而传统的FS-VTAR模型在保证其辨识优势的同时却需要针对不同时变结构选择合适的基函数形式及较高的基函数阶数,该过程相当复杂且耗时。本文借鉴无网格法中移动最小二乘(MLS)法构造形函数的思想,提出一种基于Kriging形函数的线性时变结构模态参数辨识方法。该方法首先引入自适应于辨识信号的Kriging形函数;再把时变系数在形函数上线性展开,利用最小二乘(LS)法得到形函数的展开系数;最后把时变模型特征方程转换为广义特征值问题提取出模态参数。利用时变刚度系统非平稳振动信号验证该方法,结果表明:基于Kriging形函数的FS-VTAR模型相比于传统的FS-VTAR模型能有效地避免基函数形式的选择和较高的基函数阶数,且精度相当;相比于移动最小二乘法能有效地解决其数值条件问题且具有更高的模态参数辨识精度。     

基于Lattice Bohzmann方法的翼型绕流数值模拟

将通过物理平面的结构化C型贴体网格变换到计算平面的均匀正交网格,引入非均匀网格插值方法,采用基于Lattice-Boltzmann方程的D2Q9模型,对低雷诺数下的NACA0012翼型绕流进行了数值模拟,对五个不同位置的边界层速度型与基于有限体积法的CFL3D程序的结果进行了比较.结果表明,两种方法的计算结果一致.改进后的Lattice-Boltzmann方法适用于曲边边界,非均匀网格,同时对计算速度影响不大.Lattice-Boltzmann方法计算较简单,因而用LB方法来模拟翼型绕流是一种新的选择,并可应用于更复杂的低雷诺数流动中.     

高速粘性内流的高分辨率高精度迎风型杂交格式

本文在有限体积离散和ＬＵ分解的基础上，构造出一个新的隐式迎风型杂交格式，并用于求解定常流动的稳态解。     

非均匀网格湍流大涡模拟高精度有限体积解法

为准确预测不可压复杂湍流,提出了一种可用于大涡模拟均匀或非均匀网格上的高精度有限体积法。该方法空间离散采用有限体四阶紧致格式,时间推进采用四阶Runge-Kutta法,压力-速度耦合应用四阶紧致格式的动量插值。通过直接求解顶盖驱动方腔流动证实了该方法具有近四阶的空间精度;并在此基础上,采用动态Smagor-insky亚格子应力模式,成功地实现了充分发展槽道湍流和后台阶湍流流动的大涡模拟计算,所得结果与直接数值模拟结果吻合良好,且采用非均匀网格可在比均匀网格数少的离散系统上得到同样满意的结果。结果表明,该方法是实现高精度湍流大涡数值模拟的一个有效途径。