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

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

模型大迎角高速动态特性与数据精度分析

为满足新一代高机动飞机气动性能评估、控制系统精确设计与高机动作战指标实现的需求,模型高速风洞大迎角俯仰动态特性探索及其试验数据精度的确定势在必行,且具有十分重要的工程意义。选取70°三角翼模型、SDM和Su-27飞机模型,在FL-24风洞的大振幅俯仰动态试验技术平台上对动态气动特性与试验数据精度进行了研究,获取了70°三角翼模型、SDM和Su-27飞机模型动态气动特性与重复性试验结果。研究结果表明：试验条件下,3种模型的动态数据精度较高,基本达到了高速风洞大迎角常规测力试验数据的精度水平。     

一种复杂空间飞网系统参数优化设计方法

由轻质软绳索编织而成的空间飞网是为非合作目标捕获而提出的空间系统概念,在空间碎片和废弃航天器处理方面具有很大的应用潜力。从平台抛射出后,飞网在空间形成不稳定的网形,且网形变化规律受初始参数设计的影响较大。针对空间飞网系统设计与试验中系统参数匹配问题,本文提出以容错值作为飞网展开性能的定量描述,从捕获任务的层面,建立面向捕获容错的空间飞网系统参数优化数学模型;以抓捕固定距离、确定大小的目标任务为算例,联合Isight优化平台与ANSYS/LS-Dyna求解,得到飞网系统最优参数匹配,算例仿真结果表明结果的适用性;最后,利用试验设计和极差分析方法验证最优点的稳定性。研究的模型与方法为开展空间飞网系统地面及空间试验等工程应用提供理论依据。     

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

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

针刺C/C复合材料拉伸强度及渐进失效数值预测

基于金相显微镜观测的针刺C/C复合材料细观结构,考虑材料内部纤维的真实分布,建立了针刺C/C复合材料单胞模型。采用Linde失效准则,考虑纤维渐进损伤对材料进行了刚度折减,通过引入周期性位移边界条件,对针刺C/C复合材料的拉伸破坏进行了有限元法（FEM）数值模拟,分析了单胞模型的渐进失效过程,并预测了材料的拉伸强度。开展了针刺C/C复合材料拉伸试验,采用扫描电子显微镜（SEM）观察了材料的断口形貌,数值预测结果与试验吻合良好。针刺C/C复合材料受拉伸载荷后发生准脆性断裂,拉伸应力-应变曲线呈双线性,0°无纬碳布发生纤维断裂和基体开裂破坏,90°无纬碳布出现横向基体劈裂,最终断裂发生在针刺纤维与面内无纬碳布交叉区域。     

含激波的湍流流动高精度大涡数值模拟方法

针对采用亚格子模型进行含激波的湍流流动模拟时会面临激波附近的精度损失问题,考虑从通过亚格子模型以及数值模拟方法两方面的改进来实现湍流流动大涡模拟的精度提高.大涡模拟采用了Yee及Sj(o)green (2009)提出的高阶低耗散方法.该方法采用自适应的流场探测器以控制计算中所需区域的数值耗散,并考虑对动力学模型采用在激波位置使用Sagaut和Germano(2005)提出的单边亚格子过滤器和(或)直接禁用亚格子项等方法加以改进.对于标准的马赫数1.5和3条件下的激波-湍流干扰问题,上述新方法相较于全区域采用亚格子模型的方法均表现出了相似的精度提升.同时实现的数值精度改进方案采用了Harten的亚单元分辨过程来定位和锐化激波,并在精确激波位置附近的网格点处采用了单边测试滤波.     

一种鲁棒的概率核主成分分析模型

大数据时代面临的数据维数越来越高,对数据降维处理越发显得重要。经典的主成分分析模型已被证明是一种有效的数据降维方法。但它在处理非线性、存在噪声和异常点的数据时存在效果较差的问题。对此,文章提出了一种鲁棒概率核主成分分析模型。该模型将核方法与基于高斯隐变量模型的极大似然框架相结合,用多元 t分布作为先验分布,以同时解决主成分分析在这 3个方面的弊端。提出混合鲁棒概率核主成分分析模型,使其可直接用于对混合的非线性数据进行降维和聚类分析。在不同数据集上进行的实验结果表明,与标准的混合概率核主成分分析模型相比,文中模型在数据聚类方面有更高的准确率。     

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

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

第一届航空CFD可信度研讨会总结

总结了第一届航空可信度研讨会(AeCW-1)的数值模拟结果。AeCW-1会议由中国空气动力学会计算空气动力学专委会主办,并在第一届空气动力学大会上做了专题研讨。本次研讨会吸引了14家单位参与,共提交了26组计算结果。活动的主要目的是针对典型运输机巡航构型,评估国内CFD当前技术状态,探索CFD进一步发展的研究方向,促进CFD验证和确认工作的开展,为大型客机的研制提供技术支撑。AeCW-1选择的研究模型为中国空气动力研究与发展中心(CARDC)设计的单通道运输机模型(CHN-T1),CHN-T1包括机身/机翼/平尾/立尾等部件,分别在CRADC的2.4m×2.4m风洞(FL-26)、欧洲的2.4m×2.0m风洞中(European Transonic Wind Tunnel,ETW)开展了风洞试验。AeCW-1组委会提供了系列结构网格和非结构网格,研究工况主要包括网格收敛性研究、抖振特性研究、雷诺数影响研究,风洞模型支撑装置和静气动弹性变形对数值模拟结果的影响是重要的研究内容。本文对提交组委会的数值模拟结果进行统计分析,并与FL-26及部分ETW的风洞试验结果进行了对比,对未来的研究工作提出了意见和建议。     

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

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