跨音速非定常二维翼型绕流计算

引言 为了计算跨音速非定常气动载荷和了解其流场的性质,需要对翼型绕流进行数值模拟。在构造一个跨音速非定常流动的计算方法时,重要的是要考虑能够处理流场中不断运动的激波;其次要能够处理各种类型的翼型运动,并具有高的计算效率。     

前缘机动襟翼气动力设计计算研究

本文应用欧拉方程对翼身组合体加装前缘机动襟翼构型进行数值计算，通过跨音速气动特性的计算和分析，得出前缘机动襟翼对翼身组合体气动特性的影响规律，为某型机加装机动襟翼的改进设计提供数值计算依据。     

用于多级压气机非定常流动的谐波平衡方法研究

为了在保证计算精度的同时,进一步提高压气机流动非定常模拟的计算效率,在课题组自行开发的结构化有限体积解算器上实现了谐波平衡法。以1-1/2级轴流压气机为研究对象,采用谐波平衡法和相滞后法对其非定常流场进行了数值模拟,并采用谐波平衡法对该压气机的时序效应进行了初步研究。计算表明,随着谐波阶数的增加,谐波平衡法的计算结果逐渐接近于相滞后法的计算结果;在谐波阶数达到5阶以后,谐波平衡法的计算结果不再随谐波阶数的增加而变化,且与相滞后法的计算结果匹配;谐波平衡法在保证与相滞后法相同的计算精度的情况下,计算效率相比相滞后法高了近3倍;5阶谐波平衡法即可很好地模拟入口导叶-静叶时序效应。算例表明,谐波平衡法相比相滞后法具有更多的选择空间和更高的计算效率。     

导弹气动特性在亚跨音速下的风洞试验研究

在导弹的设计过程中,导弹的气动特性作为重要因素直接影响导弹飞行的动态品质。在亚跨音速段气动特性呈现剧烈非线性的情况下,工程估算以及CFD数值计算方法所能提供的气动计算精度有限,导致对舵效特性的辨识精度较低,需要进一步采用风洞试验的方法精确计算气动参数,进而确定导弹的舵效。本文应用风洞试验方法研究导弹飞行马赫数在亚跨音速段对导弹气动特性的影响。研究结果表明：亚音速时导弹的气动特性基本一致,跨音速时发生剧烈的非线性变化;导弹的俯仰舵效先增加后减小,滚转舵效先减小后增大。结论对导弹控制律的设计以及后续的工程型号研制有参考价值。     

同伦法在现代飞机尾旋预测中的应用

在利用ＢＡＣＴＭ法对高性动性能飞机尾旋预测中，首先要求解平衡解，为此需要使用静态，强迫振荡天平气动模型和旋转天平气动模型，而在两气动模型转换过程中，常常遇到难找或不断得到平衡解的问题，针对这一问题，利用同伦法对气动模型进行处理，有效地了解决了该问题，同时，结合某战斗机进行了有关尾旋的计算和预测，结果表明，所提出的方法是可靠和切实可行的。     

锥导乘波体气动代理建模方法研究

传统的气动计算方法计算繁琐、计算效率低,不适应于乘波体多学科设计优化,通过建立气动代理模型可以很好地解决气动计算精度和效率的矛盾.利用面元法进行气动估算,采集了锥导乘波体在设计点、非设计点的气动特性作为训练数据,构建了Kriging和LS-SVM代理模型,对比了两种模型对此高维问题的代理效果.结果表明,Kriging代理模型能更准确地表达锥导乘波体的气动特性,应用代理模型进行优化等工作的计算效率与传统气动计算方法相比有显著的提高.     

谐波平衡法在非定常流场中的应用

谐波平衡法（HBM）适用于模拟周期性非定常流场。求解这类问题时,只需要一个周期中几个等距时刻的流场解,即可重建整个周期流动的时间历程,计算效率较高;同时还采用谱算子精确求解流场变量对物理时间的导数项,保障了算法的计算精度。鉴于此,详细介绍了谐波平衡法的原理和实现过程,并给出了如何在现有计算程序下实现谐波平衡法的流程。取跨声速翼型和超声速钝锥两个周期性强迫俯仰振荡绕流为例,以双时间步方法为参考,对谐波平衡法的效率和精度作了详细的考核和分析。计算结果表明,谐波平衡法的内存消耗较大,谐波数N_H=2时约为双时间步方法的4倍,但谐波平衡法取得与双时间步方法相近结果的计算时间仅为后者的1/5。因此,谐波平衡法是一种有应用前景的工程实用方法。对于长周期问题,谐波平衡法的优势更加明显。     

应用遗传算法的气动优化设计

首先对简单遗传算法进行改进,并对改进前、后算法的计算效率进行了比较;然后以跨音速翼型为例,应用改进的遗传算法进行气动优化设计。经过优化的翼型气动特性垢改善,表明了遗传算法在优化设计中的有效性。在优化设计的过程中,翼型用解析函数的线性叠加表示,目标函数和个体和适应值由欧方程的流场解来提供。     

一种民机颤振模型配重设计优化方法

民机跨音速颤振模型配重计算应考虑模型肋板不平行、肋板与刚轴不垂直的情况,现提出一种基于刚轴方向的质量、惯量线性拆分方法,并将该方法拓展到模型的展向重心位置配平计算中。某型民机垂尾跨音速颤振模型运用该方法进行了配重计算,试验结果显示,该垂尾跨音速颤振模型质量、惯量配平精度高,展向重心位置配平效果好,模型的动力学特性与设计目标相符。实践表明,该方法符合民机跨音速颤振模型设计规范要求及相应适航要求,可以推广到其他跨音速颤振模型的配重计算中。     

跨音速三维多升力面气动设计

发展了一种跨音速多升力面的气动设计方法和设计程序,它基于已成功应用于亚、跨音速机翼设计和亚音速双翼面设计的"余量修正迭代"概念.当升力面上出现超音速区和激波时自动引用迎风格式对控制方程进行修正.开发了一系列接口程序,包括目标压力设计程序.由此气动设计程羊、TAU程序以及相应的接口程序建立了跨音速多升力面气动反设计软件系统.用两个鸭翼-机翼构型验证了设计方法和设计程序,结果表明在高跨音速下设计迭代有很好的收敛性.     

弹性机翼的跨音速压强分布计算

本文采用迭代求解的方法计算弹性机翼的跨音速压强分布。用跨音速小扰动势流的二级近似方法求解气动力;结构变形则采用一维简单梁理论或二维平面结构的柔度影响系数矩阵进行计算。对M6弹性机翼所作的示范性计算表明,计算弹性翼的迭代次数与刚性翼计算相当,弹性变形对冀面气动载荷分布和激波强度均有明显影响。     

用于压气机流动计算的3种模型比较

为了更好地对压气机流动进行模拟,在课题组自行开发的结构化有限体积解算器上实现了用于压气机流场计算的混合平面法、谐波平衡法及相滞后法.以NASA Stage 35为例,对3种方法的计算结果进行了比较分析.结果表明:相滞后法的计算精度最高,混合平面法的计算精度最低;相滞后法与半环的双时间推进法结果相近,计算速度提高了20倍;相比混合平面法,谐波平衡法能准确地模拟动静叶间的非定常干涉及进出口参数变化;在谐波阶数达到5阶后,谐波平衡法计算得到的结果不随阶数变化,且与相滞后法的结果基本吻合;混合平面法的计算效率远高于另外两种方法,相滞后法与谐波平衡法在谐波阶数为5阶时的计算效率相当.      

跨音速翼型气动优化设计方法

 提出了一种有效的跨音速翼型气动优化设计方法。翼型的流场解由欧拉方程的数值解提供。带约束条件的优化计算分别采用了直接法和间接法两种优化算法。算例结果表明,本方法提供了一种跨音速翼型改型设计及新翼型设计的有效工具。     

RAE2822翼型跨音速流动CFD计算的可信度分析

采用格点格式的有限体积法对RAE2822翼型的气动特性进行了粘流数值模拟,流动模型为雷诺平均NS（RANS）方程。在此基础上系统研究了不同网格类型以及不同湍流模型的计算精准度。通过与经典的试验数据进行确认研究,得出了适宜二维跨音速粘性流动计算的网格拓扑分布规律,并总结出了适宜的湍流模型。     

一种快速模拟振荡叶栅非定常流的数值方法

为了提高叶轮机内周期性非定常流数值计算效率,发展了一种计算振荡叶栅非线性非定常流的谐波平衡方法,将周期性非定常流动参数采用傅里叶级数表示,从而将传统的定常流动数值计算方法应用于非定常流问题求解中,大大提高了非定常流计算效率。计算结果表明,发展的谐波平衡方法与线性方法、非线性时域法计算结果吻合较好,且其计算效率比时域法高1个数量级。另外,系统研究了谐波阶次对计算效率与计算精度的影响,认为在气流未出现严重分离时,应采用低阶谐波平衡方法,以减少计算耗时。此外,从压气机振荡叶栅非定常流场计算结果可以看出,叶片吸力面诱导涡的生成频率与振荡频率并不一致,涡生成频率要高于振荡频率。     

高超声速热化学非平衡钝体绕流数值模拟

对高超声速热化学非平衡钝体绕流进行研究,发展了高温热化学非平衡条件下的高超声速气动热数值模拟方法。采用二阶精度迎风TVD格式的N-S方程有限体积法多块结构网格求解器,考虑了空气5组分及化学非平衡效应的影响,对组分生成源项及振动能源项采用了点隐式的处理,提高了计算效率。对球锥模型和RAM-CII飞行实验等算例进行了数值气动热模拟,分析了不同非平衡模型在热化学非平衡条件下流场的影响。计算结果表明建立的数值模拟方法具有较高的精度。     

方腔流动气动噪声边界散射数值预测方法研究北大核心CSCD

方腔流动噪声问题因其流场变化复杂且剧烈而倍受关注,本文主要开展高亚声速方腔流动气动噪声数值预测方法研究。基于任意边界条件的格林函数解和Lighthill声模拟理论,提出可以考虑空间边界影响的气动噪声积分计算方法。数值模拟包含流动和噪声计算两部分,通过二阶精度的DDES模型进行流动数值模拟,边界积分方法计算散射声场分布。数值结果显示声场分布随时间呈现周期性变化,与流场的脉动及其脉动周期一致。观察点的声压级随频率逐渐下降且在谐波频率突然增大。本文计算结果与高精度计算气动声学方法计算结果相符,表明该方法合理、可靠,并且具有较高的计算效率。     

高超声速飞行器表面温度分布与气动热耦合数值研究

针对高超声速飞行器热防护设计中的高温气体非平衡效应问题和气动热环境精确预测问题,基于流场的非平衡Navier-Stokes方程、表面的能量守恒方程和内部的热传导方程,考虑流场的非平衡效应、表面的热辐射效应、催化效应和烧蚀效应以及热防护层内部的热传导效应,建立了初步的表面温度分布与气动热的耦合计算方法,完善了高超声速飞行器气动物理流场计算软件(AEROPH_Flow)。在表面材料为碳-碳(C-C)条件下,对飞行高度为65km和飞行速度为8,10km/s的半球以及飞行高度为50km和飞行速度为8km/s的球锥模型,开展了表面温度分布与气动热的耦合计算,验证了计算方法和计算软件,分析了表面温度分布对气动热环境的影响。研究结果表明:表面温度分布对气动热的计算结果有较大影响,在气动热环境的预测中,不仅要考虑热化学非平衡效应和表面催化效应的影响,还要考虑表面温度分布的影响,最好是采用表面温度分布与气动热耦合计算的方法,以减小表面温度分布对气动热计算结果的影响。为此,需要发展完善非平衡流场/表面催化和烧蚀/热传导温度场(气/表/固)的计算模型、耦合求解技术和计算软件,实现对高超声速飞行器的真实飞行条件下高温气体非平衡效应和气动热环境的精确模拟。     

方腔流动气动噪声边界散射数值预测方法研究

方腔流动噪声问题因其流场变化复杂且剧烈而倍受关注,本文主要开展高亚声速方腔流动气动噪声数值预测方法研究。基于任意边界条件的格林函数解和Lighthill声模拟理论,提出可以考虑空间边界影响的气动噪声积分计算方法。数值模拟包含流动和噪声计算两部分,通过二阶精度的DDES模型进行流动数值模拟,边界积分方法计算散射声场分布。数值结果显示声场分布随时间呈现周期性变化,与流场的脉动及其脉动周期一致。观察点的声压级随频率逐渐下降且在谐波频率突然增大。本文计算结果与高精度计算气动声学方法计算结果相符,表明该方法合理、可靠,并且具有较高的计算效率。     

跨音速轴流压气机动叶弯曲的数值优化

采用优化软件ISIGHT作为平台,集成了三维粘性流计算软件NUMECA和轴流压气机叶片造型程序,建立了一套叶轮机械气动性能优化系统。针对跨音速轴流压气机动叶和静叶的周向弯曲进行了优化。结果在整条特性线上,该级压气机的效率和压比都得到了不同程度的提高,充分验证了该叶轮机械气动性能优化系统的可靠性。同时表明转子的尖部反弯可以有效地减小叶尖间隙的横向流动。