基于三维空域网格的低空飞行航迹战略规划方法

针对复杂低空环境下航空应急救援飞行安全问题,提出了一种基于三维空域网格的飞行航迹战略规划方法。将空域网格化分成多个飞行航迹节点,同时考虑地形、气象、飞行规则以及航空器性能等多种约束,使用改进的A*算法搜索单个航空器的最优飞行航迹。在单个航空器初始飞行航迹基础上,结合时间窗原理,提出两种方法解决了多机无冲突航迹战略规划问题。最后,通过案例仿真验证了方法的有效性。     

可压缩粘性流动笛卡尔网格虚拟单元方法研究

以二维高雷诺数可压缩粘性流动问题为背景,提出了一种全新的笛卡尔网格虚拟单元方法。基于壁面函数基本假设,构造了壁面函数-虚拟单元方法(WF-GCM),用于定义湍流壁面边界条件。引入参考点的概念计算虚拟单元上的基本变量与湍流变量值,定义了"非贴体"笛卡尔网格下的湍流壁面边界条件,并通过壁面函数模型修正近壁面单元与界面单元。基于自适应笛卡尔网格体系,采用发展的具有二阶精度的格心格式有限体积求解器,数值模拟了跨音速RAE2822翼型绕流问题与超音速圆柱绕流问题,计算结果与实验值吻合良好,显示了WF-GCM对高雷诺数可压缩粘性问题是有效的。     

某中空长航时无人机机载浮标分离安全性研究

研究机载浮标投放过程中气动干扰特性及浮标运动规律,对机载浮标安全投放方案的设计具有重要意义。采用计算流体力学方法对机载浮标投放过程进行数值模拟。通过动态嵌套网格技术模拟浮标运动,结合刚体六自由度运动方程求解三维非定常 N-S ( Navier-Stokes )方程,模拟基于察/打一体无人机平台的浮标自由投放过程,获得详细的包括浮标运动姿态、运动轨迹和受力情况等在内的数据信息;通过算例,验证本文数值模拟程序的准确性。结果表明:浮标投放过程中,轻质浮标口盖自由分离会对螺旋桨后推式无人机的安全性造成重大影响;考虑浮标口盖固连在浮标降落伞上的状态下,初速度为 0m / s 时投放浮标整体(包含浮标、降落伞、口盖)离螺旋桨最近距离为 0.151m ,安全余量不够;初速度为 1.8m / s 时投放浮标整体离螺旋桨最近距离为 0.671m ,满足投放安全性要求。     

约束阻尼结构力学性能研究与参数优化

飞行过程中的气流扰动、机械振动使得飞机典型薄壁结构承受宽频振动,产生振动疲劳从而影响飞机结构的寿命和安全性,约束阻尼结构处理是一种经济有效的抑振手段。采用黏弹性阻尼材料和单向玻璃纤维复合材料制备两种约束阻尼结构,研究它们的阻尼和力学性能;利用有限元仿真手段对约束阻尼结构修补飞机典型悬臂梁结构的参数进行设计和分析,考察阻尼层厚度、约束层厚度以及补片尺寸对结构振动响应的影响,得到约束阻尼参数规律,并以参数优化后的约束结构为基础设计振动测试实验。结果表明：约束阻尼处理后的悬臂梁结构能有效地降低结构危险点的应变幅值;三层阻尼结构抑振效果优于双层约束阻尼,能够显著降低结构危险点在承受振动加载时的应力值,提高结构振动疲劳寿命。     

基于径向基函数与混合背景网格的动态网格变形方法

发展了一种基于径向基函数(RBF)插值和混合背景网格映射的动态网格变形方法。混合背景网格由运动物面附近的各向异性网格和偏离运动物面的各向同性网格组成。物面变形或运动后,包含在各向异性背景网格中的所有物面网格点的位移由指定的物面运动规律得到,而各向同性背景网格点的位移通过径向基函数插值获得。然后变形后的计算网格利用计算网格与混合背景网格之间固定的映射关系代数插值得到。各向同性网格距离物面一定的距离,消除了物面保形对径向基函数插值的要求,从而可以减少径向基函数基点的数量,以提高插值的效率和可靠性。最后,通过二维和三维算例对该方法的变形网格质量、参数影响、网格拓扑依赖性、复杂外形变形能力及网格变形效率进行了测试。测试结果表明,该方法能够高效生成大变形下高质量的变形网格。     

一种适用于三维混合网格的GMRES加速收敛新方法

为提高流场计算收敛效率,发展了一套适用于三维混合网格Naiver-Stokes方程求解的并行广义最小残差（GMRES）隐式时间推进方法。该方法由科学计算可移植扩展工具包（PETSc）中的Krylov子空间求解器实现,线性方程系统中的系数矩阵直接以显式给出以提高算法的稳定性。为进一步提高GMRES方法的收敛速度,对非结构网格的序号进行了重排序,使得系数矩阵的非零元素尽量向主对角线靠近。利用所发展的GMRES方法,完成了对ONERA-M6机翼、AIAA阻力预测会议通用研究模型（CRM）等算例的计算,计算结果与试验结果吻合良好。通过与其他隐式推进方法进行比较,对算法的收敛特性进行了研究。结果表明,所发展的GMRES方法计算更加稳定,残差下降速度相对LU-SGS（Lower-Upper Symmetric Gauss-Seidel）方法更快,尤其是气动力系数向着收敛解逼近的速度更加明显,提高了计算效率。     

带反馈控制的基本带翼外形的耦合气动/运动数值模拟研究

介绍了一种先进的、时间精度的多学科数值模拟平台的开发和应用,该平台可用于预测飞行器机动过程中的非定常气动特性和飞行力学特性.通过耦合求解非定常Navier-Stokes方程、刚体动力学方程及控制律,对基本带翼外形进行了虚拟飞行数值模拟研究.求解器采用基于结构网格的有限体积法,时间方法采用双时间步方法.使用反馈控制器获得舵偏控制指令,重叠网格方法用来模拟物体间的相对运动.模拟结果显示控制器的增益必须正确选择才能获得正确的响应.仿真结果表明本文发展的虚拟飞行数值模拟方法能够预测包含动边界的复杂非定常流动,具有较强的工程适用性.     

基于非结构网格流场超大规模并行计算

大规模并行的计算流体力学已成为现代航空工业研发的核心手段之一。基于非结构混合网格和有限体积法,发展了适用于工业级复杂外形气动计算的并行流动数值模拟方法。文中首先介绍了紧致数值离散格式、基于Metis的分布式多核系统网格分区技术、并行边界虚拟单元技术和MPI并行实现等相关算法。采用网格量相对较小的旋成体构型绕流模型对比分析多核并行计算结果与单核计算结果以验证并行计算的正确性,比较了不同并行规模下并行效率和残差收敛情况。然后通过对上亿网格单元的运输机复杂构型绕流模拟,开展并行效率的测试,结果表明,本文方法并行加速性能高,直到多达18816核并行效率都保持在80%以上。     

基于非结构混合网格的CHN-T1标模气动特性预测

本文使用自行研制的基于非结构混合网格的亚跨超声速流场解算器MFlow,对AeCW-1提供的客机标模CHN-T1进行了数值模拟研究。介绍了非结构混合网格的生成情况,重点分析了网格收敛特性、压力分布、气动特性曲线、流动分离等。计算得到了近似线性的网格收敛特性。随着网格加密,对激波和分离气泡的模拟更精细。尾支撑对气动特性的影响非常明显,特别是对平尾气动特性有很大影响。机翼静气动弹性变形的影响主要是使升力系数和阻力系数减小。湍流模型的QCR修正对大迎角计算结果有较大的影响。计算结果表明,MFlow程序能够准确地预测客机标模的气动特性。     

基于嵌套网格的CHN-T1标模数值模拟

为响应第一届航空CFD可信度研讨会(AeCW-1),对CHN-T1标模翼身-尾翼组合体构型采用两种RANS求解器进行了基于结构化嵌套网格的网格收敛性研究和抖振特性计算,并与对接网格结果进行了对比分析。参考AeCW-1网格生成指南自主研发了一族以约3.3倍规模增长的稀-中等-密-特密嵌套网格系列。相比于分区点对接网格,融合了贴体网格和笛卡尔网格优势的嵌套网格从根本上降低了网格生成难度,单块网格拓扑形式更为合理且正交性较好,同时能够有效地平衡近场和远场的网格量。对于跨音速运输机构型,同等规模下嵌套网格的物面网格密度较对接网格更大,能够对激波和分离等复杂流动进行更理想的预测。对比计算过程中,自研求解器OFS3D表现出了较高的嵌套网格计算效率和可靠性,而NASA CFL3D求解器的嵌套网格计算能力有待进一步挖掘和验证。计算还发现,不同湍流模型得到的表面压力分布、分离区大小以及最大升力系数存在差别,而风洞模型支撑和机翼静气弹变形对飞行器力矩特性预测影响较大。