中航工业高性能计算和网格应用系统

简要介绍了网格计算在航空工业的应用情况,提出了中航工业高性能计算和网格应用系统AVICGrid的总体框架,及其在飞机高精度大规模气动计算、气动弹性优化设计和数字化装配等三个典型应用方向的网格实现技术,实现了AVICGrid原型系统,为通过网格计算技术提供全行业高性能计算服务打下基础。     

基于网格聚类的雷达信号预分选

针对现有网格聚类算法需要人为确定网格划分、边界处理精度低的问题,提出了基于固定网格划分和动态网格划分的雷达信号预分选算法。基于固定的网格划分算法对输入信号脉冲顺序不敏感,根据网格数据压缩率自适应确定网格划分和密度阈值。基于移动技术的动态网格聚类算法能够识别任意形状和大小的聚类。仿真实验表明,这两种方法能有效适用于雷达信号分选,且能很好地识别出孤立点和噪声。     

动能BGK格式模拟高雷诺数定常流的收敛性研究

为克服动能BGK格式计算耗时多、收敛慢的缺点,在格式中引入当地时间步长、隐式LU-SGS方法和多重网格技术,基于RAE2822翼型粘性跨声速绕流的数值模拟,对其收敛性的改善进行了研究。结果表明:上述3种加速手段对求解定常问题时的收敛性有明显的改善,达到定常解所需的迭代步数随着CFL数的提高而减少;因而求解某类问题时,隐式格式比显式格式更具优势;多重网格既适用于显式格式,也适用于隐式格式,加速性更加突出。     

同向旋转螺旋桨对着陆构型气动影响分析

基于分区对接网格技术分别生成翼身区域和螺旋桨区域网格,并通过求解定常N-S方程模拟螺旋桨滑流影响下全机流动。通过对某螺旋桨飞机的着陆构型(机翼+短舱+螺旋桨+机身+平垂尾)进行数值模拟,结果与实验值吻合较好。模拟结果表明:螺旋桨滑流会改变机翼当地绕流迎角、影响压力分布;当螺旋桨同向旋转时,机身两侧襟翼的分离会出现在不同位置并改变升力在机翼展向的分布,进而对全机产生一个滚转力矩。     

一种基于三阶WENO重构的无网格算法研究

提出了基于三阶WENO重构的无网格算法,以期替代传统的无网格线性逼近重构,提高无网格算法的求解精度。基于无网格点云卫星点分布特征,通过沿点云中心点与卫星点连线方向引入局部一维坐标系,并结合虚拟点的设置,构成了在点云中实施三阶WENO重构所要求的模板。算法涉及的卫星点连线中点处的左右状态值,则利用模板上各点的流场值,采用WENO重构方法计算给出。对于模板中设置的虚拟点上的流场值,本文则利用已有的无网格点云信息,基于最近节点流场值插值得到。文中给出了采用上述WENO重构方法,结合Roe的近似Riemann求解器确定通量,并采用三阶TVD Runge-Kutta方法进行时间推进求解Euler方程的具体实施过程。对模型问题的典型流动进行了计算分析,通过比较数值解和精确解,验证了算法获得的数值解能逼近三阶精度。给出的激波管流动和超声速半圆柱绕流算例展示出所发展的算法捕捉激波等间断问题具有更高分辨率,能体现出WENO重构"基本无振荡"的特性。文中最后给出了非定常激波过弯道绕双圆柱流场的数值模拟算例,展示了所发展的算法处理含非定常激波演化的复杂流动问题的效果。     

基于动网格高空长航时机翼优化

分析了机翼设计优化的一般策略,针对高空长航时飞行器的气动性能需求和机翼构型特点,提出以机翼剖面翼型沿展向分布的设计优化;针对Hicks-Henne翼型参数化方法后缘不光顺的缺陷,提出了改进Hicks-Henne翼型参数化方法;针对通常基于线性插值的三维动网格生成方式的缺陷,提出采用基于B样条函数的三维动网格生成;通过对商业软件Pointwise和Fluent的二次开发实现数据流和工作流的定制及自动化,并将所有流程在iSIGHT平台下集成;优化算例表明,基于动网格的机翼优化实现提高功率因子的目标,是高空长航时机翼优化中的一种有效方法。     

风力机标模非定常数值模拟中的影响因素研究

基于自主研发的“亚跨超 CFD 软件平台”(TRIP3.0),采用刚性运动网格技术和动态拼接网格技术,开发了针对旋转类机械的非定常求解模块。本文开展了 NREL Phase VI 风力机标模非定常数值模拟中的影响因素研究,影响因素主要包括不同时间步长的影响、不同湍流模型的影响、刚性动网格技术和动态拼接网格技术的影响三个方面。本文数值模拟采用的控制方程为雷诺平均 N-S 方程,采用有限体积法离散控制方程,离散方程的时间方向采用“双时间步”方法求解,空间方向无粘项离散采用 MUSCL-Roe 格式,湍流模型采用 SA 和 SST 湍流模型,并引入多重网格和并行计算技术加速求解。数值模拟结果表明：时间步长、湍流模型、网格方案等因素主要影响风力机叶片吸力面的流动结构,进而影响吸力面的压力分布,而对压力面的流动结构和压力分布基本没有影响;采用刚性运动网格结合 SA 湍流模型所得到的压力分布更接近实验值。     

适用于涡激振荡问题研究的并行高精度方法

提出了一种用于研究大振幅及物体间小间距等极端情况下涡激振荡（VIV）问题的高精度方法。该方法针对三角形/四边形非结构混合网格,采用高精度谱差分（SD）格式对Navier-Stokes方程进行空间离散。通过引入非均匀滑移网格方法,将计算域分割成多个互相不重合的子区域,从而实现了子区域网格的独立变形,子区域之间通过粘接元进行信息传递。采用全耦合方法精确求解VIV问题中流体和固体间的相互作用。通过研究并行策略以及使用消息传递接口,实现了求解器的并行计算能力。数值模拟表明：对于无黏和有黏流动问题,该方法均能保持SD方法的高精度特性;对于动网格下的定常均匀来流,求解器能够做到自由流保持;单个圆柱VIV问题仿真与现有文献结果符合较好,验证了方法的可靠性;对于流场变形情况比较复杂的涡激振荡问题,求解器可以有效实现网格变形,并具有理想的并行效率。     

HiLiftPW-3高升力构型数值模拟

基于雷诺平均Navier-Stokes方程和拼接网格技术,采用MUSCL-Roe格式和Spalart-Almaras一方程湍流模型,对第3届高升力构型性能预测会议提供的两组高升力标模进行了数值模拟,主要目的是确认本文计算方法模拟复杂高升力构型的能力。研究内容主要包括高升力构型网格生成技术、网格收敛性研究及气动特性数值模拟。通过与JAXA （Japan Aerospace eXploration Agency）提供的测压、测力试验结果的对比分析,表明,在失速迎角之前,数值模拟得到的气动力系数和压力分布均与试验结果吻合;较好地模拟了局部外形变化引起的气动特性差量。本文建立的数值模拟方法对典型运输机三段翼布局的低速问题具有良好的适用性,可以为大飞机低速构型的气动设计及评估提供技术支撑。     

喷流流动控制增升效应数值模拟

飞机的增升系统是影响其起降性的重要因素,为了研究飞机下表面喷流流动控制增升效应,采用数值模拟方法与结构网格求解不可压Navier-Stokes方程,对NACA0012翼型和机翼进行研究。通过改变喷气口的参数来研究其对翼型和机翼气动特性的影响。对于翼型,主要研究喷气口速度和位置对翼型气动特性的影响;对于机翼,主要研究喷气口区域沿展向分布和喷气速度对机翼气动特性和能量利用率的影响。计算结果表明:采用喷气控制可以获得较好的气动特性,增升效果也比较明显,控制效果和上述几个参数均有一定关系。本次研究为机翼喷气口设计以及机翼局部喷气控制展向分布的选取提供了依据。