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

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

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

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

基于偏相关性分析的多设计要求快速优化设计方法

针对飞行器气动外形优化设计中的多设计点、多目标、多工程约束等多设计要求问题,提出了一种新优化设计方法———多设计要求快速优化设计方法。该方法立足于统计学中的偏相关性分析和线性回归理论,通过对设计变量和设计要求进行偏相关性分析,采用线性回归的方法构建设计变量与设计要求之间的转换模型,将多个设计要求转化为对设计变量的约束,从而简化了优化设计模型,减少了优化过程中调用求解器次数,提高了优化效率。利用该方法对RAE 2822和HSNLF(1)-0213翼型进行了多设计要求优化设计算例验证,并将优化结果与Pareto多目标优化方法结果进行对比,结果证明了本文方法的有效性和可靠性。     

基于放大因子与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%.      

一种用于飞机结冰的海绵冰模型及应用

基于两相流理论,采用欧拉法对飞机结冰进行数值模拟,在空气流场速度分布的基础上,求解水滴流场,得到水滴收集率和撞击极限.结果表明:欧拉法比拉格朗日法在计算水滴收集效率和撞击极限方面有优势,尤其是在较大的水滴平均等效直径情况下.对翼型的冰形和圆柱表面的温度预测数据和实验数据进行对比,圆柱冰形厚度误差在15.0%之内,计算结果与实验结果吻合较好,同时表明结冰模型是有效可行的.      

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

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

飞翼布局无人攻击机气动设计研究

飞翼式布局在气动效率、隐身特性和结构强度等方面具有的突出优势使其成为未来无人攻击机的一种良好的布局形式.以0.7马赫为巡航设计点,设计了一架起飞总重在18吨左右,最大载荷2吨,最大航程在5000km以上的执行对地攻击任务的无人攻击机.以欧拉方程数值解法加粘性修正的方法对其气动性能进行了计算和评估,并在初始外形的基础上进行了改进,进一步提高了升阻比,并改善了力矩特性.结果表明该设计思路具有一定的可行性.     

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

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

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

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

大型水陆两栖飞机侧风起降的滑流影响分析

大型水陆两栖飞机AG600的动力装置为安装在机翼上的4台同向旋转涡轮螺旋桨发动机，针对1：15缩比模型带动力风洞试验显示的螺旋桨滑流对侧风起降状态的偏航力矩不稳定影响，对全机带动力风洞试验模型进行了大规模并行非定常数值计算，再现了风洞试验现象，通过流动机理分析明确其产生原因主要是左侧滑时右外翼分离和垂尾背鳍涡破裂，这些原因和数值模拟的准确性也为后期的风洞试验所证实。考虑到模型风洞试验中尺度限制造成的低雷诺数和高螺旋桨转速，为保证飞行安全，继续采用该非定常方法对全尺寸飞机真实侧风起降状态进行了详细数值分析和偏航稳定性评估。研究结果显示，在飞行雷诺数和螺旋桨转速下，相同侧风范围内风洞试验显示的流动不稳定因素基本消失，偏航稳定性允许的侧风范围明显增加。本研究实现了四发螺旋桨飞机起降状态横向气动特性的滑流影响非定常数值分析，建立了基于计算流体力学的风洞与飞行雷诺数效应的相互关系，进行了偏航稳定性的虚拟试飞评估，研究成果也为AG600飞机的首飞和飞行试验所验证。