数值模拟方法对NASA CRM模型阻力预测的影响

在复杂工程外形的数值模拟中，网格类型、规模和分布、湍流模型、数值格式等都会不同程度影响模拟结果。如何评估这些模拟方法对模拟结果的影响，并识别对模拟结果影响较显著的因素，对于CFD的发展方向有积极的借鉴指导意义。为了综合研究不同因素对商用运输机外形阻力预测的影响，以AIAA第六届阻力预测会议外形NASA CRM为研究对象，同时考虑了网格、湍流模型、无黏通量格式和体心梯度求解方法等因素对阻力预测的影响。分析中采用枚举法和正交试验设计两种方法，并使用了基于聚类分析的定性敏感性分析方法和基于Mckay主影响分析的定量方法，识别出对阻力预测影响较大的因素，这为数值模拟方法的发展指明了方向。     

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

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

小展弦比飞翼标模雷诺数影响数值模拟研究

采用数值模拟方法开展小展弦比飞翼标模的雷诺数影响研究。使用自行研制的多块对接结构网格亚跨超声速流场解算器程序Mbflow,计算了试验雷诺数、二倍试验雷诺数和飞行雷诺数等三种雷诺数情况下小展弦比飞翼标模的流场。通过对计算结果的分析研究,得到了不同马赫数(Ma=0.2、0.8和1.5)和攻角情况下雷诺数对小展弦比飞翼标模的气动特性和流场特征的影响规律。小攻角情况下,雷诺数主要影响摩阻的大小;而大攻角时,雷诺数对压阻也有明显影响。在亚声速时,雷诺数主要影响分离涡的起始位置和强度;跨声速和超声速时,雷诺数还会影响到激波的位置和强度。进一步研究了小展弦比飞翼标模的自准区雷诺数问题,发现试验雷诺数接近于小展弦比飞翼标模的自准区雷诺数。     

面向非结构混合网格高精度阻力预测的梯度求解方法

针对非结构混合网格的特点,通过改进传统的Green-Gauss梯度求解方法,提出了一种可提高非结构混合网格黏性计算精度的节点型Green-Gauss梯度求解方法。利用改进后的方法,完成了DLR-F4翼身组合体算例的计算和对比分析。改进后的梯度求解方法残差收敛更好,下降量级更多,阻力系数和试验吻合更好,激波区域压力分布和分离区域流场细节的模拟更精确,说明改进后的梯度求解方法有效提高了程序的鲁棒性和阻力预测精度,验证了方法的有效性。采用改进后的方法对第5届AIAA阻力预测研讨会的通用研究模型（CRM）进行了详细的模拟分析,结果表明：改进的梯度求解方法更加适用于非结构混合网格的黏性计算,计算精准度达到国际同类CFD软件水平,进一步验证了改进方法的可靠性。     

非结构网格体心梯度求解方法的精度分析

对于空间二阶精度的非结构网格求解器,无黏通量和黏性通量的计算都需要流场变量在单元体心处的梯度值。体心梯度求解的精度很大程度上决定了求解器整体的精度。本文推导了几种常用的梯度求解方法(包括格林高斯法,节点格林高斯法和最小二乘法)求解体心梯度的理论精度,并从数值方面对上述理论精度进行了验证。在数值测试过程中,通过以当前单元体心为基准进行坐标局部缩放的做法,保证了在非结构网格上做精度测试时,网格拓扑能够严格保持不变。理论推导和数值测试一致发现:使用保线性权插值到节点的节点格林高斯法、使用共面单元的最小二乘法和使用共点单元的最小二乘法,不管网格拓扑关系如何,都能保证梯度求解为一阶精度。而格林-高斯法和使用距离倒数权插值到节点的节点格林-高斯法在一般的网格中只有零阶精度,只有在特定的网格和权函数关系下,才能有一阶精度。最后用NACA 0012翼型亚声速绕流模拟进一步验证了梯度求解方法的收敛精度。     

飞机结冰的不确定性量化研究进展

为深入认识影响飞机结冰的不确定性及其研究方法,从自然条件结冰、冰风洞试验、数值模拟等方面,介绍了飞机结冰不确定性来源,以蒙特卡洛法、多项式混沌法和随机配置法等为例,系统分析了各种不确定性量化方法在计算能力、求解精度等方面的优劣。考虑飞机结冰不确定性量化处于起步时期,重点综述了数值模拟中结冰条件不确定性对冰形和气动特性的影响。从结合单步法和多步法确定最佳结冰时间步长来提升结冰计算精度和效率、对其他关键部件进行结冰不确定性量化从而为更精细化的防/除冰系统的设计提供支撑,以及通过建立高精度的代理模型来代替原本复杂的数值模拟系统以应对考虑多个不确定性因素共同作用所带来的计算挑战等多个方面,全面展望不确定性量化方法及其在飞机结冰应用中的发展方向。     

混合型三阶格式及关于两种DES算法的比较计算研究

为开展RANS/LES混合模拟,在传统三阶迎风偏置格式的基础上,提出一类混合型的三阶计算格式,格式在中高波数范围具有可调的耗散水平.在此基础上,文中通过计算比较研究两种RANS/LES混合算法:一种是基于Spalart-Allmaras一方程模型的DES模型(DES-SA),另一种是基于混合长度模型的DES模型(DES-ML),使用 发展的三阶格式对圆柱绕流进行模拟,并将得到的结果与文献比较,对计算格式、DES-ML算法进行了初步探讨.     

SA一方程湍流模型参数影响分析与辨识

SA(Spalart-Allmaras)一方程湍流模型是目前工程湍流计算中主要采用的湍流模型之一。首先,针对某典型战斗机的小攻角、中等攻角、大攻角流动工况,利用均匀试验设计方法分析SA一方程模型中8个参数取值对上述工况下飞机升力和阻力系数计算结果的影响规律;然后,建立工程湍流模型参数的辨识方法,并将其用于该战斗机大攻角工况下湍流模型参数cb1值的辨识调整。结果表明:不同工况下,湍流模型参数对计算结果的影响规律不同;在附着流状态下,对升力和阻力影响较大的参数是cv1和cb1;在中等攻角和较大攻角下,对升力和阻力影响较大的参数是cb1;适当减小参数cb1的取值后,升力和阻力系数的计算结果有较明显的改善,这可能与飞机大攻角分离流场中涡粘系数和剪切应力的发展与自由剪切流存在一定差异有关。     

基于数据挖掘的飞行器气动布局设计知识提取

为了更深入地理解飞行器气动布局设计优化中多目标/多设计变量间的影响关系,提高优化模型的科学性及优化效率,对基于数据挖掘技术的飞行器气动布局隐含设计知识提取问题开展了探索研究。以高升阻比滑翔飞行器布局设计优化问题为例,基于当前比较有代表性的方差分析、等度量映射、决策树、自组织映射4类机器学习算法对气动布局优化设计中产生的中间数据进行了挖掘分析。对不同方法得到的升阻比、横/侧向稳定性及容积率4种目标性能间的权衡关系,目标性能与设计变量间的敏感性关系及产生较优布局外形的设计变量取值规则进行了综合对比分析,凝练形成了适用于该类飞行器的设计知识,同时对4种方法的特点及适用性进行了总结分析,给出了相关结论。