三维分区对接网格Euler方程的并行数值模拟

本文根据CFD中分区对接网格的特点，应用基于工作站集群系统的MPI网络并行环境，实现了典型双三角翼模型亚声速绕流Euler方程的并行计算。文中对计算结果的收敛情况、耗时、加速比及效率等并行性能指标进行了讨论。结果表明，多域多机与单域单机的计算结果基本吻合；2节点并行效率达到88．21％，4节点并行效率达到78．99％。     

在非结构动网格上追踪多介质界面的ALE方法

在处理带有介质分界面的多介质流动问题时,把介质边界当作一种特殊的边界,并且在该边界上的点都有着双重的定义,分别对应边界两侧的不同介质状态,通过利用HLLC方法求解任意拉格朗日-欧拉方程(ALE),以及运用界面跟踪和动网格方法分别对一维和二维多介质流进行了数值模拟,并得到了满意的结果.     

非结构动网格技术及其在超声速飞行器头罩分离模拟中的应用

针对激波诱导物体运动的二维动力学模型,比较了采用流动定常、结合非定常弹道方程的所谓"准定常"计算方法和非结构动网格计算方法,通过流场和运动特性分析,指出发展有相对运动的非定常流动模拟方法的重要性;另外采用修正弹簧近似和局部重构的非结构动网格技术,耦合求解Euler方程及6DOF弹道方程,数值模拟高度H=25km、飞行马赫数Ma=6.0、迎角α=0°状态下火箭抛罩的动力学过程.研究表明,稠密大气层内头罩受到的气动力比惯性力大几个数量级,导致头罩的运动特性完全由气动力控制.分离过程中,头罩和火箭头部之间的气动干扰在时间和空间上都表现出强烈的非线性特征,头罩姿态变化剧烈.本文还对分离过程中部分位置上的流场进行定常模拟,比较了定常和非定常的气动力,所得结论对于分析地面风洞可控轨迹系统(CTS)试验具有指导意义.     

三维舵面绕流的N-S方程计算研究

基于粘性非结构混合网格,采用格心格式的有限体积法,应用Spalart-Allmaras湍流模型进行绕三维舵面流动的N-S方程数值模拟。操纵面铰链力矩的计算结果与风洞试验结果的对比分析表明,计算方法对粘性流动具有良好的适应性。对舵面铰链力矩随迎角、马赫数、雷诺数变化的规律也进行了分析研究。     

基于Euler方程结合附面层修正方法的应用研究

采用欧拉加附面层修正的方法在结构网格上对民用飞机翼身组合体流场进行了数值模拟研究.采用一种适合翼身组合体模型的的网格生成技术来生成计算用的结构网格,该技术是先采用一种基于Higenstock源项修正的椭圆型偏微分方程生成二维网格,然后沿机身轴向合理的光滑性调整成三维网格.在求解Euler方程时采用中心离散格式和显式时间推进,并利用附面层修正的方法来考虑粘性的影响.为了验证此方法的正确性,本文对DLR-F4翼身组合体进行了数值模拟,结果表明此方法和试验结果能够很好的吻合.     

无人机类脑吸引子神经网络导航技术

当前无人机在非结构化或未知环境下飞行主要采用SLAM进行导航与定位,存在如下突出问题：依赖高精度昂贵激光雷达等环境感知传感器;需要建立准确世界和无人机物理模型;受环境影响较大;自主智能水平较低,无法较好地满足无人机对导航系统的要求,需要发展自主智能的导航方式。基于吸引子神经网络的类脑导航技术,无需训练模型参数,不依赖高精度传感器,无需精确建模,且复杂环境下鲁棒性较强,具有解决上述问题的潜力。简要阐述了动物大脑导航机理,分析了吸引子神经网络和基于吸引子神经网络的类脑导航关键技术,最后讨论了吸引子类脑导航技术在无人机应用中的挑战。     

化学非平衡流动解耦算法的回顾与新进展

超声速化学非平衡流动的数值模拟一直是计算流体力学领域的难点,主要体现在如下几个方面：物理过程非常复杂,存在着激波、燃烧波等各种复杂波系的相互作用;超声速化学非平衡流动属于典型的时空多尺度物理问题,其控制方程存在严重的刚性,给数值求解带来了很大困难。对国内外的解耦算法研究现状简单综述后,主要介绍1993年刘君提出的解耦算法的理论基础,流动方程采用冻结流模型,源项方程模拟流体微团在当地绝热、定容的热力学系统内发生的化学反应过程。通过引入两个中间变量,即等效内能和等效比热比,将与温度无关的生成焓从流动方程组能量项中分离出去,源项方程组中包含等效内能,使用不同算子对流动方程和源项方程解耦求解。与传统解耦算法相比,源项方程的求解过程中包含状态参数和组元同时变化。结合刘君解耦算法机理和有限体积法空间平均特性,介绍近期在提高算法计算效率方面的研究进展,包括流动方程优化算法和耦合过程优化。采用优化算法对经典的激波诱导燃烧算例进行数值模拟,与不同文献结果进行对比,验证了优化算法的时空精度。通过总结经验发现化学非平衡反应仅发生在流场局部区域,提出质量生成率判据,结合相应模拟结果验证了方法的可行性,可以进一步提高化学反应算子的计算效率。     

基于惯性补偿的三分量测力天平瞬态特性的研究

高超声速一体化试验模型在脉冲燃烧风洞中进行测力试验时,模型振动导致天平无法准确测量模型所受到的气动力载荷。为研究脉冲风洞瞬态测力问题,采用了以下方法:根据测力天平的结构特点建立了动力学方程;对其进行了虚拟标定和模态分析;对测力天平进行了瞬态分析并对输出结果进行了惯性补偿。结果表明,单分量阶跃载荷加载时,补偿后输出结果超调量大大下降,振动衰减时间缩短,俯仰力矩会对轴向力输出产生干扰,且三分量阶跃载荷加载时,干扰降低,各分量的超调量分别为162.6N,574.4N和38N·m;单分量正弦加载时,加载分量上输入输出基本一致,俯仰力矩对轴向力的干扰作用仍然存在,3分量正弦载荷加载时,输入输出结果具有相同的周期特性,俯仰力矩结果与输入结果一致,轴向和法向输入输出间存在一定偏差,3分量的超调量分别为24.1N,375.7N和70.8N·m。     

叶轮行波振动下的气动阻尼计算

为展开运动边界下离心叶轮流场的数值分析,独立开发了网格变形程序和非定常流动模拟程序,实现了流场中振动离心叶轮的气动阻尼计算。采用紧支径向基函数法进行结构到气动表面变形的数据传递,应用二叉树技术进行壁面距离的计算,大幅提高了网格变形和流场分析中距离搜索的计算效率。通过振动叶栅和径向叶轮的算例,验证了程序应用于运动边界流场计算和叶轮流场模拟的正确性。以半开式叶轮为对象,展开其在行波振动条件下气动阻尼的计算分析。结果表明,叶轮在前、后行波振动下的气动阻尼均为正,后行波振动引起的气动推力和气动阻尼较大,叶轮表面的平均气动功率密度呈现循环对称的分布特征。     

航空发动机两级反向旋流燃烧室燃烧流场大涡模拟研究

为深入了解真实航空发动机内燃烧流场,采用动态亚网格模型结合单步快速化学反应(FC)、火焰面(FLM)及反应进度变量(FPV)等三种燃烧模型对径向两级反向旋流燃烧室单头部构型进行了大涡模拟。为避免模型简化误差,对燃烧室包括全部气膜冷却孔在内的精细结构进行了完全仿真。在已达到统计定常状态的冷态流场基础上首先模拟了燃料喷注和掺混过程,约2.6ms后燃料到达真实的点火位置,随后采用FPV模型在半径3mm的球形区域数值模拟了点火,展示了在主燃孔横向射流作用下初始火焰沿化学恰当比混合分数等值线传播并充满整个火焰筒的发展过程,结果显示火焰到达燃烧室出口的耗时约为6～7ms。不同模型算法预测的平均温度场与CARS测量结果作了对比,LES-FPV,RANS-FPV,LES-FLM以及参考文献中RANS-FLM计算平均误差分别为3.47%,4.17%,7.76%和11.22%,表明LES改进了模拟精度,且FPV模型显著优于FLM模型。RANS-FPV预测的出口存在严重热斑,导致其给出的出口温度分布因子(OTDF)及最大径向温度分布因子(RTDF)值分别达到0.593和0.313;LES-FPV结果均匀性最好,其预测值分别为0.284和0.193。