基于欧拉方程的跨声速翼型设计

给出了一种基于欧拉方程的跨声速翼型设计方法。方法以Ｔａｋａｎａｓｈｉ提出的“正反迭代余量修正”设计原理为基础，在气动力分析模块中，以欧拉方程为基本控制方程，并采用Ｗａｌｚ边界层方法对其进行粘性修正；反设计模块采用经过改进的二维翼型设计方法。方法的程序经过几个设计实例主题证明是有效而实用的。     

基于多项式混沌法的翼型不确定性分析及梯度优化设计

耦合伴随方法和非嵌入式多项式混沌法，发展了高效、可靠的不确定性梯度优化设计方法。利用伴随方程法求解目标函数对不确定性变量的导数，发展了一种梯度增强型多项式混沌法。通过亚声速和跨声速下等多种算例可以证明该方法可以提高不确定性分析的效率和精度。同时，利用基于方差分解的全局敏感性分析方法对不确定性变量的敏感性进行了量化。建立了多项式混沌耦合伴随方程的统计矩梯度求解方法，并结合梯度增强型多项式混沌法搭建不确定性梯度优化设计系统。基于该优化设计系统对二维低亚声速和跨声速翼型开展确定性及不确定性优化设计研究。优化结果显示，相比于确定性优化设计，不确定性优化设计通过合理权衡确定性性能和不确定性性能，可提高抵抗马赫数和迎角不确定性扰动的能力，同时优化性能均值和标准差。其中阻力系数均值最大可降低17%，阻力系数标准差最大可降低80%。而确定性优化设计可能导致性能鲁棒性的降低。     

基于离散伴随的高超声速密切锥乘波体气动优化设计

基于定平面形状的密切锥乘波体设计方法能够显著提高传统乘波体的设计灵活性和整体升阻特性。但是该类乘波体在设计时忽略了三维效应、黏性效应以及头部/前缘的钝化效应，在设计工况下仍会出现溢流，升阻比难以达到最优；另外，这类乘波体仍具有传统乘波体在偏离设计条件下气动特性会出现恶化的不足。因此，有必要在考虑黏性的情况下，针对定平面形状的密切锥乘波体开展全机气动优化设计。结合基于全速域通量求解方法和RANS湍流模型的高精度CFD求解器、鲁棒的结构网格变形方法、自由变形参数化方法、离散伴随方法以及序列二次规划算法，实现了基于离散伴随的高超声速飞行器气动优化设计方法。基于上述方法，针对定平面形状的密切锥乘波体开展了单点和多点的三维整机气动优化设计。在400万多块结构网格、600个设计变量以及303个设计约束条件下，所采用的离散伴随优化方法仅花费2 240CPU小时和3 360CPU小时即完成了三维整机单点和多点的优化设计。结果表明，相较于初始构型，单点优化得到的构型在设计状态下的升阻比提升了近5%；多点优化得到的构型可保证在设计点状态升阻特性没有损失的同时，将非设计点的升阻比提升10%以上，进而在一定程...     

基于放大因子与Spalart-Allmaras湍流模型的转捩预测

为了验证放大因子输运方程与Spalart-Allmaras(S-A)湍流模型耦合对转捩现象的模拟精度,选取Schubauer and Klebanoff(S-K)平板、S809低速翼型、30p30n多段翼型以及复杂的三维HiLiftPW-1构型进行自由转捩计算,并将计算结果与实验进行比较分析,其中针对S809算例,还与Langtry-Menter(L-M)转捩模型进行了比较.算例结果表明:放大因子输运方程与S-A湍流模型的耦合能够较好的捕捉转捩位置以及转捩发展过程,对分离泡诱导的转捩的模拟相比L-M转捩模型更精确,转捩位置的捕捉精度提升了10%;对比实验,多段翼转捩位置的捕捉误差最大为6.5%;针对三维高升力增升构型,以实验作为参考,全湍流计算与考虑边界层转捩的对比显示考虑边界层转捩能够更加精确的模拟气动力系数,升力和表面摩擦阻力系数的模拟精度精度提升1%.      

基于参数敏感性分析的反设计方法

在基于类别形状函数变换(CST)的反设计方法中,通过考察CST参数化与反设计过程,完成了基于参数敏感性分析的反设计,实现了翼型的快速、准确反设计;研究了Bernstein多项式阶数进化对反设计过程敏感性的影响。研究结果表明,在亚音速状态下,采用反设计方法能使设计压力分布较好逼近目标压力分布;在跨音速状态出现激波时或在计算过程出现迭代发散时,加入阶数进化手段后的敏感性分析方法能改善激波区域的敏感度,使迭代较好收敛。方法能够反映各种光滑翼型压力残差和翼型几何外形改变的关系,在给定状态下用初始翼型快速且比较准确地拟合目标翼型。     

基于当地变量的横流转捩预测模型的研究与改进

Langtry和Menter提出的转捩预测模型需要改进以具备预测横流转捩的能力。当地变量Helicity参数可以指示边界层内的横流信息,因而可用来构造适用于复杂构型以及当代计算流体力学(CFD)并行计算的横流转捩预测模型。实现了基于Helicity参数的横流转捩预测模型,对于后掠角为45°的NLF(2)-0415无限展长后掠翼,模型能够预测不同雷诺数对横流转捩的影响,但是对6:1椭球的横流转捩预测结果与试验数据相差较多。针对实现的横流转捩预测模型的缺点,考虑横流速度因素进行改进。横流速度的求解经过简化近似可以当地求解,因而保证了改进的模型完全基于当地变量的优势。采用改进后的横流转捩预测模型分别对NLF(2)-0415机翼、6:1椭球以及DLR-F5机翼进行数值模拟,并与试验数据进行对比分析,结果显示改进后的横流转捩预测模型可以较为准确地捕捉横流转捩现象。     

某飞翼布局隐身飞行器的翼型优化

对某飞翼布局无人机的翼型进行考虑隐身迎角的气动、隐身综合优化研究。采用基于非均匀有理B样条曲线的自由变形方法实现翼型的参数化表达;无限插值动网格方法生成气动求解的计算网格;分别采用基于雷诺平均N—S方程求解的计算流体力学方法及物理光学法计算翼型的气动、隐身性能;选择多目标遗传算法进行综合寻优设计研究。优化翼型兼顾了气动与隐身的要求,其综合性能优于原翼型。结果表明,方法具有较好的实用性。     

考虑吸气影响的层流翼型梯度优化设计

针对层流翼优化设计问题,推导考虑吸气影响的层流-湍流转捩耦合伴随方程,并结合链式求导法则、自动微分技术及CK(Coupled Krylov)算法实现耦合伴随方程的高效求解,构建基于离散伴随理论的混合层流(HLF)翼梯度优化设计系统。转捩预测采用基于扰动放大因子模型(AFM)的e^N方法。层流飞行试验结果表明,基于AFM的转捩预测方法能够有效捕捉吸气控制影响下T-S(Tollmien-Schlichting)扰动波失稳诱导的转捩现象。利用构建的梯度优化系统开展层流翼型多点设计,并与非梯度优化方法的优化结果进行对比。自然层流设计结果对比显示,不同优化方法得到的设计结果呈现出较强的一致性。混合层流翼型多点设计结果显示,混合层流优化具有比自然层流优化更强的减阻能力。相比于自然层流优化翼型,混合层流优化翼型分别减阻6.1%、5.9%、33.3%、9.5%。本文构建的基于离散伴随的层流翼梯度优化方法能使混合层流机翼的气动性能得到大幅提升,可为未来混合层流机翼减阻设计提供方法支撑。     

亚跨声速边界层增长因子输运模式研究进展

基于线性稳定性理论的e^N方法是航空航天飞行器绕流转捩预测技术中可信度较高的一种。传统的线性稳定性理论的求解过程较为复杂,在复杂气动外形表面的推广阻力较大。随着近年来当地化转捩模式的快速发展,将传统线性稳定性理论的分析流程模式化,即将线性稳定性理论分析转化为一个CFD问题,成为研究热点。其核心思想是对层流相似性解进行大量的线性稳定性分析,然后获得扰动增长因子包络线与层流基本流特征变量之间的函数映射关系,再利用输运方程的积分特性构造对应扰动增长因子的输运模式。随后对模式中出现的非当地变量进行当地化求解,最终与湍流模式耦合形成新的转捩-湍流输运模式。经过验证,当前发展的流向N_TS输运模式与横流N_CF输运模式在广泛的算例验证中取得了较高的预测精度,但还存在需要改进和完善的地方。     

基于边界层相似性解的放大因子输运模型

为了使二维线性稳定性理论能够适应现代CFD求解技术,通过求解Falkner-Skan边界层相似性方程获得各个形状因子下的相似速度型。基于线性稳定性理论,对各个速度型进行稳定性分析从而获得对应不同速度型的扰动放大因子包络线。最后使用标量输运方程的形式实现包络近似方法中放大因子的当地化求解,并结合原始γ-Re_θt转捩模型中的间歇因子输运方程,实现了自然转捩和分离泡转捩的建模。使用该输运模型对S&K平板、S809翼型、NLR7301翼型和DLR-F5机翼进行转捩预测,结果均与试验结果吻合较好,验证了该模型构建的合理性和可行性。