飞机着陆过程中非定常地面效应数值模拟方法

针对飞机着陆过程的大位移运动,提出适于多块结构动网格生成的改进谢别德-无限插值混合方法,改进谢别德方法准确计算各网格块的边界网格变形,无限插值方法快速计算内部网格点变形,结合局部重构保证动态网格的合理分布.混合方法有效解决了单一无限插值法存在的动网格交错问题.进而将二维非定常地面效应N-S方程计算方法发展至三维情形,应用于带增升装置的翼身组合体着陆过程中非定常地面效应计算分析.结果表明:模型的升力系数随离地高度的减小而减小,非定常和准定常计算方法的结果存在一定差异,而考虑了非定常效应的算法更为准确.     

以风轮气动性能为目标的风力机翼型优化设计

提出了一种优化风力机翼型气动特性的设计方法。在现有风力机翼型的基础上,通过改变其最大厚度位置至后缘部分的外形进行优化计算。采用N-S方程计算翼型的流场及气动特性,以提高风力机切向力同时减小其法向力为目的,用二维线性插值的方法进行多目标优化。以DU 93W 210翼型为例,针对风力机翼型气动性能本身的矛盾性,根据实际需求提出几种方案加以优化。结果表明,在不同的性能需求下,本文所提出的设计方法能取得较好的效果。     

基于CFD的机翼突风响应计算

利用网格速度理论,计算机翼在锐边突风和1-cos突风下的响应,研究气动非线性和自由度耦合对翼尖加速度和升力系数的影响.采用中心格式有限体积法进行空间离散,并用双时间推进法求解非定常Euler方程.计算了刚性(沉浮)和弹性机翼在锐边突风下的加速度响应,以及刚性机翼(沉浮+俯仰)在1-cos突风下的升力响应过程,并分别与片...     

跨声速机翼操纵面非定常气动力Euler方程计算

采用代数法与椭圆型方程优化相结合的结构生成方法为带操纵面F5机翼生成面对接块结构网格;弹性变形技术生成动态网格;有限体积法、双时间推进法求解非定常Euler方程.给出一种适于面对接块结构网格计算的守恒量插值方法.数值计算了F5机翼操纵面作简谐运动时的非定常气动力,计算值与试验值吻合较好.     

跨音速颤振计算方法研究

采用快捷的动态弹性变形技术生成三维块结构动态贴体网格。结构运动方程采用二级精度的龙格—库塔时间推进。气动力求解则用非定常欧拉方程的双时间有限体积推进,外时间为物理时间,与结构运动方程同步;对每个真实物理时间步采用五步Runge Kutta时间推进进行虚拟定常迭代。针对可压流颤振计算中存在的质量不相似问题,采用变质量、变刚度的方法计算出质量匹配点处的颤振速压值,并根据颤振速压随质量或刚度倍数的变化趋势,得出可用的飞机跨音速颤振速压。对一飞机简化外形的跨音速颤振特性计算得到了合理的计算结果。     

基于预处理和网格自适应的风力机翼型气动力计算

用预处理后的Navier-Stokes方程数值模拟低马赫数下风力机翼型的静态和动态绕流流动。为了提高计算的精度,运用了网格自适应技术:初始网格的流场计算收敛后,采用流场速度的紊乱度捕捉湍流区域;采用直接分裂初始父单元的H型网格加密方法对湍流区域进行局部加密。基于初始网格与自适应网格的计算结果与实验结果的对比表明,网格自适应方法提高了计算的精度。     

基于速度分裂法的翼型阵风响应及减缓数值模拟

国内现阶段飞行器阵风响应CFD模拟通常采用网格速度法（FVM）,而近年提出的更为准确的速度分裂法（SVM）目前主要针对刚性飞行器阵风响应分析。本文将SVM拓展至弹性翼型的One-Minus-Cosine阵风响应计算和减缓研究。基于动态网格系统非定常Navier-Stokes方程,将阵风条件下的速度场分解为阵风速度与背景速度的叠加,理论推导出SVM阵风模拟控制方程,结果表明,FVM是SVM忽略源项后的一种近似方法。进一步建立起弹性翼型阵风响应预测的CFD/CSD时域耦合算法和基于俯仰控制的阵风响应减缓数值模拟方法。算例分析了NACA0012刚性及弹性翼型的One-Minus-Cosine阵风响应,计算结果与文献数据一致,在此基础上分析了阵风尺度、黏性和结构弹性对阵风响应的影响。开展了NACA64A010弹性翼型One-Minus-Cosine阵风减缓模拟,结果表明：单一的沉浮速度控制输入量更有利于减缓阵风载荷峰值,而沉浮速度、俯仰角相结合的控制输入量则能同时减缓翼型的沉浮和俯仰运动。本文工作可为三维飞行器阵风响应和减缓研究提供参考。     

高超声速全动翼面全时域耦合分析方法及应用

在流场-结构温度场同步计算方法的基础上,建立了多物理场全时域耦合分析方法,将方法应用于沿轨道运动的高超声速全动翼面热气动弹性稳定性分析。采用基于有限体积模型的CFD同步计算方法求解高超声速流场和结构温度场,建立映射关系实现结构有限元模型气动载荷加载和温度场赋值。采用移动坐标系和动网格相结合的方式描述变速度飞行和翼面偏转过程。通过坐标系变换将翼面偏转过程和振动过程的网格变形量叠加,考虑翼面振动和偏转的耦合非定常效应。针对沿轨道运动的高超声速飞行器,建立了同步计算方法与全时域耦合分析方法相结合的热气动弹性稳定性分析流程。研究发现,与同步计算方法相比,全时域耦合分析方法能够模拟结构振动对流场和结构温度场的影响,计算得到的监测点热流密度波动幅值占热流峰值的10%左右,而温度变化并不明显,相比于刚体模型,监测点温度只下降了0.3%左右。全时域耦合热气动弹性分析方法得到的颤振临界点在4-5号状态点之间,颤振形式为铰链扭转模态与一阶弯曲模态的耦合颤振,与"冻结"模态的热颤振方法结果一致。     

Euler及N-S方程的二维四边形网格自适应算法

描述了一种基于贴体网格的Euler和N-S方程自适应算法。采用中心有限体积法,对Euler及N-S方程进行数值求解,N-S方程的计算中加入了基于SA湍流模型的DES模型。在流场中,运用了网格自适应算法,提高了数值计算对激波、流动分离、涡等特性的捕捉和分辨能力。采用上述方法,数值分析了单段翼型的绕流问题。文末给出的2个二维问题的对比算例表明本文方法有效。