水下超声速气体射流线性稳定性研究

水下超音速气体射流稳定性是水下推进、水下焊接技术研究的焦点领域。通过建立线化小扰动气液混合流体控制方程,开展水下超音速气体射流线性稳定性研究。控制方程的形式表明射流基本流分布(速度与浓度分布)、雷诺数、无量纲密度、无量纲粘度能够影响水下超音速气体射流稳定性。使用Chebyshev配置点法对控制方程进行求解。计算结果表明射流基本流分布以及无量纲密度是影响射流稳定性的主要因素。射流基本流分布越平缓,射流最大扰动增长率越小;无量纲密度越小,即气液密度差越小,射流最大扰动增长率越小。随雷诺数增加,射流最大扰动增长率先减小后增加,当雷诺数大于105时,射流最大扰动增长率趋于常数。对不同无量纲粘度,射流扰动增长率-特征波数曲线高度重合。同时分析了射流不同截面的稳定性特征,计算结果表明射流最大扰动增长率随距喷口距离增加而减小。这些结论有助于理解水下超音速射流物理过程,为实际应用提供理论指导。     

计算三维振动机翼上跨-亚音速压力分布的局部线化面元法

本文以小扰动位势流理论为基础,采用局部线化假设,建立了三维振动机翼的跨音速下洗积分方程;给出了跨-亚音速流核函数的计算方法;用推广的偶极子栅格法计算了矩形谐振动机翼上的非定常压力分布。结果表明,本方法可以较准确地计算非线性的跨-亚音速非定常流动,关于跨-超音速流动的数值解法,尚需作进一步的研究。     

跨音速小扰动方程的近似因式分解法

本文对跨音速小扰动方程的近似因式分解法进行了研究。详细叙述了近似因式分解法。计算表明,α_l,α_h选择恰当,收敛速度是较快的。     

绕钝前缘翼型定常跨音速小扰动势流的一种计算方法

提高钝前缘翼型的跨音速压力分布计算方法的精度与效率对翼型设计十分重要。国内对跨音速小扰动势流方法进行了研究,但对M_∞>0.8的超临界情况尚未计算。国外文献[3]指出当M_∞>0.9时完全速势方程的一种方法失败了。 本文试图改进经典的小扰动势流方程,探索比较稳定的迭代方法和超松弛方法,以克服超临界流计算常常不易收敛的困难,使小扰动势流计算的应用范围扩大到更高的M     

射流破碎的线性不稳定性分析方法

线性不稳定性分析是目前处理柱状射流及复合射流破碎的最常用理论方法之一。该方法基于流体力学基本方程和边界条件,利用线化小扰动法和正则模法开展相关研究,能够成功预测实验中发现的不同流动现象和规律。线化小扰动法的原理是将非线性问题中涉及的物理量分解成零阶近似和小扰动,代入原非线性方程并略去高阶小量,从而使非线性问题转化为一阶近似的线性定解问题。正则模法是将小扰动量分解成一系列模态的叠加(比如Fourier级数形式),使每个模态都满足线性系统,从而将线性定解问题转化为求解广义特征值问题。正则模法适用于具有对称性的流场,根据小扰动的演化情况,可以分为时间模式、空间模式和时空模式。本文回顾了射流及复合射流线性不稳定性分析方法的主要内容和具体实施步骤,并选取流动聚焦这一代表性的微细射流生成技术,概述了线性不稳定性分析方法在实际中的应用,最后对已取得的相关成果进行了总结。     

绕翼型跨音速大扰动对称势流的插值混合差分法

 本文将跨音速定常小扰动势流混合差分法推广到跨音速大扰动定常势流,提出了在局部速度坐标系中求解跨音速精确势流方程的插值混合差分法。作为算例,计算了双圆弧翼型和NACA0015翼型对称问题压力分布,并与已知实验值和双圆弧翼型小扰动混合差分法计算值进行比较,结果接近。试算表明,本文提出的插值混合差分格式是稳定和收敛的。本文解决了M_∞趋近于1的计算难点。     

直升机抛撒地雷的稳定性和动态响应研究

建立了直升机外挂抛撒地雷装置的飞行动力学模型,外挂装置的存在改变了直升机的气动特性,对直升机的稳定性也造成较大的影响。直升机受扰动后的响应主要取决于被扰动对象的稳定性和扰动向量,而直升机抛撒地雷产生的反作用力和力矩可作为主要的扰动力向量。在小扰动假设的基础上,建立了直升机抛撒物体的扰动运动方程,通过实例计算分别得出在对称-单发、非对称-单发,非对称-连续和不同飞行速度时,直升机速度、角速度和姿态角的动态响应,并对其进行对比分析。     

多孔表面抑制第二模态失稳的最优开孔率和孔半径分析

多孔介质、微槽道和超声波吸声材料等可用于抑制高超声速边界层第二模态扰动波。通过声波在无穷长小管道中的传播模型给定多孔介质表面的扰动边界条件,采用时间模式的线性稳定性理论分析了多孔介质表面对边界层稳定性的影响。对马赫数6热壁平板边界层的考察表明,多孔介质表面不仅可以大大推迟第二模态扰动波的中性点,还可以大大抑制最不稳定第二模态扰动波的幅值增长率。为了找到最优控制参数,在较广的参数范围内考察了多孔介质表面的开孔率和孔半径对第二模态扰动波幅值增长率的影响,得到了最优开孔率和孔半径沿流向的分布。最后,还考察了基本流的当地非平行性(法向速度)对控制效果和最优控制参数的影响。     

考虑运行参数耦合作用的舰船轴系动力学特性研究

为分析清楚多动态因素作用下的轴系动力学特性，以保障轴系运行的可靠性和安全性，以某桨-轴-壳体为对象，建立其有限元模型，基于有限差分法和小扰动法求解各轴承油膜压力分布和油膜刚度，基于来流法求解艉部伴流场影响下的螺旋桨激振力。在此基础上，提出中间变量法，以转速为中间变量，在系统分析转速对油膜压力分布和螺旋桨激振力变化影响的基础上，研究得到转速变化下轴系整体动力学特性规律。研究结果表明：转速的增加使相邻轴承载荷差值分别降低了13398N和11102N，改善了轴系的载荷分布状态；轴系振动响应中的共振峰值受特定转速影响较大，在53.333Hz的轴向共振频率处以及14.667Hz,26Hz的横向共振频率处，轴系各参考点共振峰值的最大值出现在190r/min运行工况下。     

渐近展开匹配法在跨音速钝头薄翼型绕流分析中的应用

一、引言 众所周知,钝头翼型的头部区域对均匀来流的扰动不再是小量,该处流动不符合小扰动假设。因此,由跨音速小扰动方程得出的解在头部区域失效。本文将流场分为外部解区,内外解区及头部附近的内内解区（见图1）。     

壁面小折角对马赫数4.5边界层中扰动演化的影响

为了利用e N方法预测带小折角平板边界层的转捩,采用数值模拟和流动稳定性方法研究来流马赫数为4.5,扩张角和压缩角分别为1°、2°、3°的高超声速边界层中扰动的演化,得到不同角度下不同频率的扰动波沿流向的幅值及增长率分布。并对工程中常用的预测带小折角平板转捩位置的e N方法进行验证分析,指出采用e N方法预测带小折角平板转捩时需要对N值进行修正的概念。     

弹性机翼的跨音速压强分布计算

本文采用迭代求解的方法计算弹性机翼的跨音速压强分布。用跨音速小扰动势流的二级近似方法求解气动力;结构变形则采用一维简单梁理论或二维平面结构的柔度影响系数矩阵进行计算。对M6弹性机翼所作的示范性计算表明,计算弹性翼的迭代次数与刚性翼计算相当,弹性变形对冀面气动载荷分布和激波强度均有明显影响。     

计算特征向量导数的若干方法的短评和一种改进模态法

 对特征向量导数计算的若干方法作了简短地评述,同时提出了一种改进模态法。这种方法具备Nelson法、改进Nelson法和直接扰动法的优点。该方法还可以具有多种形式,但本文仅介绍了一种同于直接扰动法的主要形式,这也是精度最好的一种形式。     

旋翼非定常平行桨-涡干扰流场的数值模拟

建立了一个基于Euler方程的旋翼平行桨-涡干扰流场的数值模拟方法。该方法在时间离散上使用隐式的双时间法,以避免因物理时间步长过小而引起的数值不稳定问题;为减小由于空间离散格式精度和网格密度引起的涡数值耗散,采用了预定涡方法。以Kitaplioglu的平行桨-涡干扰试验模型为算例,计算了表面瞬时压强分布,并与可得到的试验结果相比较,验证了方法的有效性,且描述了干扰流场的气动力变化特性以及桨叶表面扰动压力波的产生和传播过程,然后分析了马赫数、涡强度、干扰距离等参数对桨-涡干扰流场特性的影响。在此基础上,得出了一些有意义的结论。     

民用飞机线性化仿真分析

线性分析是控制律设计和操纵性稳定性检查的重要手段。为了实现线性分析,需要从飞机非线性模型中提取小扰动运动方程。本文介绍了小扰动理论,提出小扰动假设,根据假设得到一种线性化方法。将小扰动方程的时域响应与飞机运动学方程的时域响应进行对比,结果表明,线性化和非线性化的飞机迎角、速度和俯仰角曲线吻合很好,该线性化方法简单有效,并具有一定的工程设计指导作用。     

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

基于对线性常微分方程组中各系数相对量级的估计,本文应用奇异摄动理论对飞机纵向小扰动运动的短周期和长周期模态的经典简化理论重新进行分析.结果表明,除长周期模态特征根的实部(阻尼部分)具育约13%的相对误差外,渐近展开式的首项所给出的短周期和长周期特征根均具有很好的精确度。然后利用坐标变形法进一步修正了特征值误差中的绝大部分。最后还讨论了在不同初始扰动下的动态反应的解.     

状态变量模型建立方法

为了克服小扰动方法精度不高、稳态终值响应法动态过程不一致以及传统拟合法的随着维数增加精度下降和拟合时间长等缺点,综合了3种算法（小扰动法、稳态终值响应法和传统拟合法）的优点,设计了基于变尺度法的混合求解方法建立了航空发动机状态变量线性模型。仿真结果表明：建立的航空发动机状态变量线性模型与非线性部件级模型在动态过程响应中吻合良好,而且具有较高的稳态精度,能够保证航空发动机最终稳定值的一致性。     

基于变尺度寻优混合求解算法的航空发动机状态变量模型建立方法

为了克服小扰动方法精度不高、稳态终值响应法动态过程不一致以及传统拟合法的随着维数增加精度下降和拟合时间长等缺点,综合了3种算法(小扰动法、稳态终值响应法和传统拟合法)的优点,设计了基于变尺度法的混合求解方法建立了航空发动机状态变量线性模型。仿真结果表明:建立的航空发动机状态变量线性模型与非线性部件级模型在动态过程响应中吻合良好,而且具有较高的稳态精度,能够保证航空发动机最终稳定值的一致性。     

亚音速复杂组合体表面压力分布及气动力特性数值计算方法

本文从亚音速小扰动位势方程出发,给出了计算任意复杂组合体表面压力分布和纵向气动力特性的基本解的数值计算方法。把全机表面划分成许多网格,在每个网格上布置基本解以模拟物体,用物面上没有法向气流穿透的边界条件决定基本解的未知强度,从而计算出物面压力系数,积分出总气动力特性。该方法适用于任意形状的机身及任意平面形状、任意厚度、弯度和扭转分布的机翼(尾翼),机翼(尾翼)可以带有两段后缘操纵面。为降低对电子计算机存贮量、速度和稳定时间的要求,本方法采用了机翼、机身、尾翼部件之间的循环迭代求解。     

准静弹性飞机动态特性的计算研究

本文在小扰动和准静弹性假设下对弹性飞机纵向动态特性进行了计算和分析。对弹性变形对动态特性的影响进行了初步讨论。本文采用格林函数法求得飞机的气动力,用正则模态法计入弹性变形影响并引入准静弹性假设,得到了研究纵向动态特性所需的全部刚体飞机气动导数和弹性飞机气动导数。利用这些导数,对准静弹性飞机动态特性进行计算研究。并将计算结果与刚化模型的同一飞机进行比较和分析。