超声速平头圆柱绕流DSMC/EPSM混合算法研究

本文将EPSM算法与DSMC方法结合，构造了可模拟含近连续流区及过渡流区的DSMC／EPSM混合算法。运用混合算法模拟马赫数等于5．37时超声速平头圆柱绕流，并与DSMC结果进行比较，证明了DSMC／EPSM混合算法的有效性，并对两种算法的计算效率进行了比较。     

二维不可压流程序在PVM环境下的并行计算

本文在PVM环境下进行了二维不可压流的并行数值模拟,应用并行虚拟机模型对机器环境和并行程序的特性进行了分析,研究了并行程序和并行机器的匹配问题.在MPP环境下的测试表明,改进后的并行程序并行效率高.     

燃气源复杂通道三维冷流数值模拟

采用两方程RNGκ-ε湍流模型，通过非结构网格上的SIMPLE算法，对某远程火箭弹控制系统燃气源复杂通道内流场进行了冷流数值模拟，获得了通道内各气流参量的三维分布，了解了燃气源内各处特别是控制气流流经的迷宫形多层导热套环结构内的气流流动特性，计算结果与实验结果基本吻合。     

不可压扰动方程高精度对称紧致差分数值解法及应用

从粘性不可压扰动方程一阶改型形式出发,对其实现了高精度对称紧致差分离散,就导出的扰动线性特征值问题给出了一个高铲双重迭代局部解法,以相同精度将特征值和特征函数同时得到。通过不可平面Ｐｏｉｓｅｕｉｌｌｅ流时间稳定性算例详细对比显示了算法良好的谱分辨能力和较弱的网格依赖性,并结合复矩阵广义特征隐式移ＱＺ算法获取了一个 扰动特征值谱计算结果。     

喷管潜入段二维两相冷流实验

建立了喷管潜入段二维两相冷流实验设备,进行了相应的实验。实验清晰地显示了在潜入喷管进口背壁凹区旋涡的产生、稳定、消失,以及粒子的冲刷、粘贴、脱落、沉积等反复作用的复杂现象,它们与流动的雷诺数、粒子的含量和尺寸、潜入喷管进口处的形状和潜入深度、背壁型面等因素有关。     

迟滞非线性二元机翼颤振特性分析

 采用多项式迟滞非线性模型建立二元机翼气动弹性运动方程，并用数值积分法进行求解。通过系统响应振幅随来流速度变化的分叉图和频谱分析发现，俯仰方向由于含有非线性因素，振动中的高阶分量随速度提高不断增加，并引起高次分叉。重点研究“机翼/空气”质量比以及“沉浮/俯仰”两个自由度的自然振动频率比对非线性颤振速度边界的影响，并提出可以通过提高自然振动频率比来减小迟滞非线性因素的不利影响。     

考虑旋翼调制影响的直升机RCS特性分析及评估

考虑旋翼调制的影响,建立了适合于直升机全机雷达散射截面(RCS)特性计算的“面元-边缘法”.兼顾计算效率和精度,在外形变化剧烈地区域,网格进行加密处理,并保证电磁网格密度和尺度满足雷达波波长大小的比例关系.以桨盘倾倒方式计入桨叶的挥舞和变距运动,采用准静态法模拟旋翼对全机RCS特性的影响.分析了某直升机RCS极化、姿态、频率响应特性,并根据直升机RCS和雷达探测距离关系,提出了4级预警机制和角域范围.研究表明:旋翼转动时全机RCS动态响应具有连续性和对称性,振荡区散射水平强,RCS幅值为-5~12dB ·m2;奇数片桨叶比偶数片桨叶的RCS减缩2~5dB ·m2,且有利于控制直升机RCS包络线和散射峰值的时域响应,增强雷达隐身性能.      

航空发动机加力燃烧室跨流气冷稳定器流场研究

针对1种航空发动机加力燃烧室用引气冷却的跨流稳定器结构形式,采用标准k-ε湍流模型数值模拟方法,对其流场结构进行了分析和研究,将数值计算结果与封闭风洞中PIV流场测试结果进行了对比。结果表明:跨流气冷稳定器流场结构复杂,呈3维分布,环形稳定器后方的流场结构与普通钝体的类似,带后掠角的径向稳定器后方回流区与环形稳定器回流区互相耦合为整体,与沿环形稳定器展向的回流区形成高度一体的复杂回流区结构,形成稳定的低速区,起到良好的火焰稳定作用。     

解耦N-S/DSMC方法计算推力器真空羽流的边界条件研究

准确模拟羽流流场是合理评估空间推力器真空羽流效应的基础。因羽流流场同时包含连续和稀薄两种流动机理,通常采用解耦方式的Navier-Stokes(N-S)方程和直接模拟蒙特卡罗(DSMC)混合方法对其进行模拟。为保证解耦N-S/DSMC方法应用于羽流计算时的准确性并尽量提高其计算效率,对计算中的DSMC入口和喷管壁面等边界条件设置问题开展了研究。通过与文献中低总压工况的羽流试验数据比较,确定了合理设置DSMC入口边界位置、壁面反射类型、壁温等边界条件的方法。针对将此方法应用于实际推力器工况计算时效率过低的问题,通过多种数值试验,表明从喷管出口处开始沿KnGL为0.05的等值线作为DSMC入口界面可同时保证仿真精度和较高的计算效率;并证明实际工况下壁温设置对羽流场仿真结果影响不大。     

Ma5斜激波串动态特性实验研究

为加深对激波串动态特性的认识,进行了简单管道模型在Ma风洞中的实验研究,探讨了两种不同背压变化条件下斜激波串的流动结构和动态特性。背压的生成和控制通过模型尾部两块斜板的闭合运动实现。为了更好观察实验现象,实验过程中动态压力传感器和高速纹影实行同步采集。研究发现:随着背压升高,激波串逐渐前移,由对称形态发展为非对称形态,并持续整个移动过程。背压增加的速度对激波串前缘的前移速度没有影响。激波串在管道内的前移过程并不是一个恒定运动过程,而是存在有稳定前移和急剧前移两种状态,其特性与管道内流场结构、壁面压强分布、背压大小以及距离管道出口的距离等因素相关。在不同的背压条件下,管道内的斜激波串表现出不同的频谱特征,但在同一工况中,激波串区域内不同位置处壁面压强的频谱特性相似。在Case 1状态中,压强振荡主频为f512Hz在Case 2状态中,振荡主频为f578Hz次级频率约为f260Hz。两次实验的主频均大于Pipon.niau模型理论计算值,而Case 2中频率f和声学振荡频率相近。