一种高效的组合动导数分量数值模拟算法

基于传统网格速度法,发展了在流场网格不变的情况下,数值模拟沉浮运动、纯俯仰运动和俯仰振荡运动的算法,提出了一套直接求解旋转阻尼导数和洗流时差导数的方法,该方法不需要实时更新网格,减少了计算时间和所需内存,避免了负体积的出现。首先通过NACA 0006的阵风响应和NACA 0012翼型的强迫俯仰振荡,验证了程序和方法的正确性;然后数值模拟了BFM(basic finner missile)标准模型的沉浮运动、纯俯仰运动和俯仰振荡运动,并计算了其旋转阻尼导数、洗流时差导数和组合动导数,所得结果与动网格方法和风洞试验值吻合,最大误差不超过4%,进一步验证了该方法的正确性。     

飞行器动态稳定性参数计算方法研究进展

动态稳定性参数（简称动导数）是飞行器控制系统设计、飞行器动不稳定发生边界分析及相应动态稳定性判据研究的关键气动参数。在对飞行稳定性问题进行概述的基础上,介绍飞行器动态稳定性参数数值模拟的国内外研究进展。并按照理论方法、工程近似方法及计算流体力学（CFD）模拟方法的动导数发展方向对近年来主要的动导数计算方法进行了综述,评价了各种动导数预测方法的优缺点,指出了动导数数值模拟在理论基础、非定常气动力建模、预测方法精度和效率等方面存在的问题。最后对动导数数值模拟的发展趋势进行了展望。     

战斗机垂尾脉动压力数值模拟

在亚跨超计算流体力学（CFD）软件平台（TRIP）上开发了基于RANS/LES混合思路的IDDES流动模拟技术,并通过NACA0021翼型60°大迎角分离流动与串列圆柱绕流模拟对RANS/LES混合方法的精确度进行了验证,针对某战斗机外形的垂尾脉动压力开展了数值模拟研究。战斗机来流马赫数为0.1,基于全机长度的雷诺数为2×10⁶,模型迎角为20°、30°和40°。分别通过脉动压力系数、脉动压力功率谱密度、空间流动结构以及侧向力响应曲线等对战斗机的垂尾脉动压力进行了分析。脉动压力模拟结果表明：当垂尾完全沉浸在边条翼脱体涡破碎后的宽频湍流脉动气流中时,垂尾翼梢位置的脉动压力会发生明显的增大。     

弹性飞机跨声速机动载荷计算方法

基于跨声速非定常气动力求解的重叠场源法,开展了弹性飞机跨声速机动载荷计算方法研究,为现代飞机结构强度设计提供更加可靠和精确的临界载荷计算方法。首先通过采用重叠场源法求解关于62%根弦俯仰振荡的LANN机翼在马赫数为0.822的非定常气动力并与实验结果进行对比,验证了重叠场源法对跨声速激波效应预测和非定常气动力计算的能力;其次应用频域气动力有理近似技术拟合场源法计算得到的广义气动力系数矩阵,建立了机动载荷分析的状态空间模型;然后完成了某型民用飞机俯仰机动载荷分析,研究了俯仰机动飞行情况下飞机机体状态量及飞机部件载荷响应规律。计算结果表明：考虑机体弹性变形后,机翼和平尾气动载荷响应最大值分别减小了5.1%和10.6%,升降舵气动载荷响应最大值增大了16.2%,在飞机结构强度设计中必须考虑机体弹性效应对飞机部件载荷的影响。     

临界冰形确定方法及其对气动特性影响研究

临界冰形是指在适航规章结冰包线内,每个可用飞行构型下,对飞机操纵性和稳定性影响最严重的冰形,临界冰形分析是飞机适航取证中的重要工作。对临界冰形确定方法,及临界冰形对气动特性的影响进行了研究。发展了基于 CFD 方法计算临界冰形的一般方法,包括临界冰形分析状态、敏感性分析截面确定、结冰参数敏感性分析、临界结冰条件确定、临界冰形确定等。流场计算采用中航工业空气动力研究院气动力计算平台(UNSMB),基于 Jameson 中心格式的有限体积法求解 N-S 方程;水滴撞击特性计算采用 Eulerian 方法求解水滴轨迹运动方程;结冰计算采用经典的 Messinger 热力学模型。选取 CRM 飞机为研究对象,以机翼外翼50%展长处为敏感性分析截面,在典型飞行条件下,分析了结冰对环境温度、水滴直径、飞行速度、飞行迎角等参数的敏感性。利用“几何外形敏感性分析方法”,即通过对比冰形的上下冰角角度和冰角厚度等冰形几何参数来确定最严重冰形,得到了CRM 飞机的临界结冰条件和临界冰形,其中敏感性分析截面在水滴直径为30μm 时上冰角厚度和下冰角厚度最大,冰角最大厚度约41 mm。计算了结冰后的气动性能衰减规律,临界冰形对飞机气动性能影响严重,导致升力降低6.7%~23.8%,阻力增加17%~70.9%。发展了45min 待机临界冰形确定方法,基于几何外形敏感性分析方法进行环境温度、水滴直径、飞行条件等各类参数的结冰敏感性分析,得到飞机的临界结冰条件和临界冰形,对于民用飞机设计和适航取证具有一定的工程应用价值。     

槽壁几何参数对跨声速风洞流场品质的影响

采用数值模拟方法研究了再入片开角、槽腔扩开角和开槽形状对跨声速槽壁风洞流场品质的影响。在数值模拟方面,介绍了对应的出入口边界条件设置方法,通过量纲分析得出适用于内流问题的相似参数。研究结果表明:再入片关闭工况下,其前缘处形成了声速喉道,难以建立稳定的试验段超声速流场;再入片打开工况下,试验段流场的不利扰动主要来自于从驻室到试验段的槽壁出流现象,其根源是驻室涡系产生的垂直于槽壁方向的冲击气流,较大的槽腔扩开角对冲击气流有一定的阻挡作用。槽壁各部件对试验段流场品质的影响均是通过增强或抑制槽壁出流来实现的,其中开槽平面形状的影响尤为突出。采用基于自由变形(FFD)的优化设计方法得到了适用于试验段马赫数0.8到1.2状态的开槽形状,优化结果显示初始构型中存在的槽壁出流大部分得到消除,模型区马赫数均方差降低了一个量级,证明了该方法的合理性与实用性。     

基于熵产最小曲线型散热器参数优化设计及数值模拟

基于参数化设计和模拟方法,采用VB开发了放射状散热器参数化设计模拟软件,可自动生成高质量的六面体网格和进行数值计算.利用该软件快速对风扇散热器组件的结构参数进行性能分析和讨论,加深了对风扇散热器组件各主要参数在散热器整体性能中影响程度的认识.采用并行计算流体力学方法,基于组合优化策略以熵产最小为目标函数对多参数结构风扇散热器组件进行优化,曲线型散热器优化结果比初步设计方案的性能提高了7％.     

箔条云跨流域整体气动特性计算研究

研究随机动力学理论与统计散布取样技术,提出对飞行器抛撤"海量"箔条及形成箔条云有限分组随机统计模拟方法;采用所建立空气动力学理论与计算技术,提出可有效模拟稀薄流到连续流不同高度、马赫数、攻角与侧滑角及极高长细比箔条气动特性计算方法;研制飞行力学数值方法,计算确定任意时刻箔条个体团空间位置、速度;应用所发展随机动力学加权技术及概率统计分布原理,进行箔条云整体性能分析,建立箔条云跨流域飞行整体运动性能与覆盖范围统计模拟方法.对飞行器抛撒"海量"箔条形成箔条云跨流域飞行气动特性计算分析,模拟仿真结果与实验测试相吻合,证实所建立箔条及箔条云跨流域整体气动特性计算分析方法的工程实用与可靠性.     

悬停及小速度前飞状态下旋翼对机身及平尾气动特性的影响

介绍了一副研究型旋翼 Ｚ９ 机身组合模型在 Ｃ Ａ Ｒ Ｄ Ｃ８ｍ ×６ｍ 低速风洞的试验概况,给出了悬停及小速度前飞时旋翼对机身气动特性的影响,着重讨论了水平尾翼对机身纵向气动特性的贡献。机身横截面上的不对称压力分布测量结果,揭示了旋翼在机身表面诱导的附加横侧向气动载荷的机理。     

LU-AUSMPW混合格式在可压缩无粘和粘性流动数值分析中的应用

本文首先将 AUSMPW格式与三阶 MUSCL格式融合, 给出了其在任意曲线坐标下的三维形式, 并将其与 LU-SGS格式结合, 应用于可压缩 Euler和 Navier-Stokes方程的求解。其次, 构造了一种新的通量限制器。最后, 为了验证 LU -AUSMPW混合格式的性能, 将平面叶栅跨音速无粘流动以及喷管超音速粘性流动作为算例。本文计算结果与文献计算结果和实验数据相符很好, 表明采用 LU-AUSMPW混合格式数值模拟可压缩流场, 具有较高的计算精度、较快的收敛速度和良好的稳定性。