高升阻比自然层流翼型多点/多目标优化设计

研究并发展了一套进行高升阻比自然层流翼型多设计点/多设计目标的优化设计方法。为减少设计中的计算量,使用了XFOIL程序进行流场计算。使用基于响应面方法的优化方法进行最优化计算,并使用了不含二阶交叉项的二阶多项式模型作为响应面模型,大大减少了构造模型所需的试验次数,使进行更多设计变量和多点设计成为可能。试验点的选择满足D-优化准则。研究了设计点及目标函数的选择,进行了单点/单目标及多点/多目标的设计,结果表明:多点/多目标设计可以很好的改善单点设计中非设计点性能变差的缺点,设计结果有工程实用价值。     

Kriging模型在翼型反设计中的应用研究

目前,基于Kriging代理模型的优化方法在气动优化设计中得到了广泛的应用,然而却并未能应用于气动反设计中。本文将翼型气动反设计问题转化为优化问题,利用基于Kriging模型的优化方法进行了翼型的单目标、多目标反设计,验证了Kriging模型在翼型单目标、多目标反设计中的适应性,并研究了模型加点准则对设计结果的影响;此外,将本文方法与多项式响应面模型优化方法、基于Adjoint的优化方法进行了比较,结果表明,该方法明显优于基于多项式响应面的优化方法,且与Adjoint方法相比,设计效率相当且具有更好的通用性。最后将本文方法应用于接近工程实际的在已有翼型基础上修改压力分布进行反设计,验证了本文方法的可行性。     

无尾布局翼型的DISC设计研究

采用DISC方法耦合N-S方程求解程序研究了无尾布局翼型的设计。通过NLF1015,RAE2822翼型的反设计表明,该方法可用于亚、跨声速翼型设计,具有简单、计算量小的优点。完成了翼型改进和新翼型设计,对某BWB飞机内翼段翼型进行了改进设计,采用人工干预对现有翼型压力分布进行局部修正,达到快速削弱激波、改善翼型气动性能的目的;针对小展弦比翼身融合布局内翼段新翼型设计,采用基于多目标多约束压力分布优化的设计思想,得到满足气动和几何约束的目标压力分布。设计结果表明,数值优化目标压力分布方法较易满足总体与气动的设计目标和约束。     

跨声速层流翼型的混合反设计/优化设计方法

跨声速层流翼型设计须兼顾优良的超临界特性和自然层流特性,因而对设计方法提出了更高的要求。针对现有反设计方法和直接优化设计方法的不足,发展了一种适用于跨声速层流翼型的混合反设计/优化设计方法。该方法引入了基于经验的局部流场特征作为反设计目标,翼型性能指标作为直接优化设计目标,然后加权形成了混合反设计/优化设计总目标,并同时考虑了气动和几何约束。优化算法采用基于自适应并行加点技术的代理优化,流动数值模拟采用耦合基于线性稳定性理论的e^N转捩自动判定的雷诺平均Navier-Stokes（RANS）方程求解器。针对现代中短程民用客机需求,以NPU-LSC-72613翼型为基准,开展了层流翼型减阻的混合反设计/优化设计。分别将局部目标压力分布、总阻力作为反设计和直接优化设计目标,得到了较好的优化结果,验证了方法的有效性。经过2轮优化结果显示混合反设计/优化设计总目标显著下降。所设计翼型吸力面局部压力分布与目标压力分布基本一致,总阻力下降15.5%;吸力面和压力面层流范围均大于55%倍弦长,激波强度显著减弱,说明所设计翼型同时具有优良的超临界和层流特性。将所设计翼型配置到机翼上,通过三维数值模拟进行校验,结果显示所设计跨声速层流机翼升阻比提高了6.64%,在一定升力系数范围内,气动性能均有显著提高,验证了所设计跨声速层流翼型在机翼设计中的适用性。     

基于响应面法的低速翼型气动优化设计

响应面方法相较于其它直接优化方法有其高效、实用的优势,此前的研究更多地将响应面方法用于超音速和跨音速翼型的减阻优化设计中。本文将此方法应用于低速翼型优化设计中,进行了基于RANS(Reynolds-Aver-aged Navier-Stokes)方程和自由转捩预测耦合求解的低速翼型气动优化设计。通过计算附面层方程得到附面层参数并用en方法计算转捩位置,并考虑了T-S波和层流分离造成的转捩。RANS方程计算中,使用了转捩过渡区模型,以保证附面层外边界压力分布的精度。RANS方程和转捩预测迭代进行至转捩位置收敛。在响应面模型计算中,使用不含二阶交叉项的二阶多项式模型,减少了构造模型所需的计算量;合理的选择设计空间保证了构造的响应面模型具有较高的精度。使用三个设计点的多目标优化设计,保证了设计的合理性。通过对NACA64(1)-112翼型优化计算结果表明,本文的方法可以有效地进行低速翼型的气动优化,各设计点上转捩位置也得到了改善,有较好的工程实用前景。     

低速翼型多点／多目标优化设计方法研究

进行了基于响应面优化方法的低速翼型多点／多目标多约束设计方法研究。分析了完全二阶多项式模型（完全模型）和不含二阶交叉项的多项式模型（简化模型）对优化的影响,并分别进行了单目标和多目标的优化设计。结果表明多目标优化可以改善单目标设计的翼型性能在设计点附近取得极值的局限,而且用简化模型进行的更多设计变量的设计结果要优于完整模型的设计结果。     

现代自然层流翼型的设计方法

本文介绍一种设计跨声速自然层流翼型的计算流体力学(CFD)方法。本方法采用“正反迭代、余量修正”设计原理,通过将跨声速翼型设计软件NPU-TD2D中的反设计程序进行改进,并与含有层、湍流混合边界层修正的跨声速层流翼型计算程序DLRBGKWALZ耦合,实现了在跨声速粘性流动条件下直接设计层流翼型。亚、超临界的设计实例和风洞验证表明,本方法可以在几个设计迭代内设计出压力分布、转捩位置及气动参数均准确收敛于设计目标的新翼型,是一种设计现代自然层流翼型的有效而实用的CFD方法。     

基于响应面方法的风力机叶片多目标优化设计研究

运用基于响应面方法的优化设计技术,于径向使用NPU-WA-风力机专用翼型族的某1.5MW水平轴风力机叶片的多目标、多约束优化设计研究中。风力机气动性能使用基于叶素-动量理论的风力机性能分析和设计软件PROPID51。设计变量为叶片径向外形参数,包括弦长和扭转角分布,但是相对厚度保持不变;设计目标为年发电量和功率系数的最大化;在多目标优化中,使用"统一目标函数"法将多个设计目标函数通过加权求和统一到一个目标函数中。为减小计算量,响应面模型使用不含二阶交叉项的二阶多项式模型;构建模型中试验点的选择满足D-优化准则。以某1.5MW变速变矩型风力机叶片为例,进行了优化设计研究。叶片径向使用西北工业大学翼型研究中心设计的NPU-WA-风力机专用翼型族,使用CFD计算的气动性能数据作为输入进行了设计,分析了目标函数的权值分配对设计结果的影响;使用风洞测量自由转捩的气动性能数据进行了设计并分析了表面光滑条件对气动性能的影响。     

考虑三维流动效应的自然层流短舱压力分布反设计

针对三维流动效应下非轴对称自然层流（NLF）短舱设计困难的问题，发展了一种基于机器学习技术的压力分布反设计方法。短舱参数化建模使用自由曲面变形技术（FFD），通过求解RANS方程和基于SST(Shear Stress Trans‐port）湍流模型的■转捩模型实现自然转捩预测，利用生成拓扑映射（GTM）模型建立短舱外形及其压力分布数据集与低维隐空间变量的映射关系，全局优化算法在隐空间高效寻优，获得与目标压力分布匹配的短舱气动外形，实现自然层流短舱的反设计。GTM模型建立了数据集在高低维度间的高精度映射关系，优化设计过程中无需反复调用CFD求解器，设计效率大幅提升。利用该方法对通气短舱进行三维优化设计，设计所得非轴对称短舱外表面的自然层流区最大长度达当地弦长的40.5%，比优化前延长了12.2%，验证了该方法在考虑三维流动效应下NLF短舱的优化设计能力。     

高雷诺数下跨声速自然层流短舱优化设计

为解决高雷诺数下大涵道比发动机自然层流短舱高维优化设计问题,提取短舱3个基准面实现轴对称自然层流优化设计。通过获得较大范围层流区域,从而降低短舱表面摩擦阻力。采用类别形状函数（CST）参数化、γ-Re_θt转捩模型和遗传算法建立自然层流（NLF）短舱自动优化流程。表明通过优化短舱基准面建立三维NLF短舱设计方法的可行性。进而通过CATIA二次开发构建三维非轴对称NLF短舱,解决了基准面优化后大量数据点的高效导入和曲面生成问题。针对设计的非轴对称NLF短舱,进行了设计点附近迎角、侧滑角、来流马赫数以及湍流度的转捩敏感性分析。结果表明:跨声速状态下,迎角增大层流范围减小;来流马赫数增大层流范围扩大;侧滑角和湍流度对层流范围影响很小。      

行星际小推力轨道Lambert解及应用

 为提高行星际小推力转移轨道初始设计精度,提出了基于N次逆多项式逼近的半解析Lambert算法,并基于该算法发展了一种转移轨道初始设计方法。首先,采用N次逆多项式近似小推力轨道形状,应用推力方向假设和位置速度边界条件推导出部分系数及推力大小解析式。接着,分析了飞行时间约束和轨道动力学约束下解的存在性,并给出了关键系数的可行域。然后,利用探测器质量消耗方程建立了Lambert问题求解模型并加以解决。最后,基于所提Lambert算法,通过对连续推力约束进行降维,提出一种求解多圈非固定时间的行星际小推力转移轨道初始设计方法。分别以固定和非固定时间转移任务为例对所提Lambert算法和初始轨道设计方法进行了数学仿真,数值结果表明：相比传统6阶方法,所提Lambert算法在目标轨道半长轴为5 AU时可减少速度增量需求36.63%;所提初始设计方法与最优化方法设计结果接近,可为转移轨道的精确设计提供可行的设计初值。     

改进型Gappy POD翼型反设计方法

 为了提高基本Gappy本征正交分解(POD)翼型反设计方法的精度,在原始方法快照采样过程中,调整参数化方法,并用已产生翼型中压力分布最接近目标压力分布的翼型替换基础扰动翼型,形成最优快照替换采样法。在迭代求解阶段,根据迭代产生的压力分布与迭代产生的翼型实际压力分布之间的误差,引入校正快照,并据此调整目标压力分布,形成校正迭代法。实际算例表明,最优快照替换采样法所采集快照张成的空间较原始采样法更接近设计目标。而校正迭代法较原始迭代法能明显提高反设计精度。但最佳模态数量的选择对于Gappy POD翼型反设计方法仍然是一个难点。     

现代超临界翼型设计及其风洞试验

开展了现代超临界翼型的设计研究,对现役飞机的压力分布形态进行了分析,针对现役飞机在巡航状态和阻力发散点的压力分布进行对比,提取了现役飞机超临界剖面设计的要点。采用类函数/型函数变换(CST)参数化方法、基于二阶震荡及自然选择的随机权重混合粒子群算法(RwSecSelPSO)、雷诺平均Navier-Stokes(RANS)方程、Kriging代理模型结合定期望值型的目标函数建立了优化设计系统。针对提高阻力发散马赫数和降低巡航低头力矩的设计指标,利用优化设计系统通过调整目标期望值设计了一系列满足设计指标但阻力发散马赫数不同的超临界翼型,并选择了其中具有典型特性的翼型进行了对比分析,验证了提高阻力发散马赫数和低速失速特性的设计方法,指出了在阻力发散点形成平顶形压力分布的超临界翼型具有较好的综合性能。对设计的超临界翼型进行了高、低速风洞试验验证,试验结果表明:设计结果达到了设计指标要求,提出的低速改进方案有效,层流对超临界翼型失速特性影响较大。     

一种基于动模态分解的翼型流动转捩预测新方法

考虑自由转捩的定常/非定常流动Navier-Stokes方程数值求解,对于翼型流动细节的精确模拟和气动力的精确预测均具有十分重要的意义。采用动模态分解(DMD)方法进行流动稳定性分析,再结合e N方法,提出了一套适用于翼型绕流的转捩预测新方法,称为DMD/e N方法。相比于传统的线性稳定性分析方法,DMD方法不需要求解附面层方程和线性稳定性方程,也没有引入平行流假设,具有更好的理论适用性和算法鲁棒性。开展了NLF0416、S809和SD7003等翼型的转捩预测数值验证研究,通过与实验结果以及与传统的基于线性稳定性分析的e N方法的比较,验证了本文所发展的转捩预测新方法在预测翼型的定常流动和非定常流动转捩方面的正确性,也表明了该方法具有解决含层流分离泡的翼型绕流转捩预测的能力。     

Double-decoupled inverse design of natural laminar flow nacelle under transonic conditions

The inverse design based on the pressure distribution is an essential approach to realize the improvement of Natural Laminar Flow(NLF) performance for nacelles. However, the direct definition of target pressure distribution at design point is challenging for the dilemma to consider the constraints of shock wave and laminar flow at the same time. In addition, the universality of method will be limited when the inverse design is strongly coupled with the solver. Thus, a double-decoupled methodolog...     

多段翼型前缘缝翼位置的优化设计研究

本文讨论优化前缘缝翼位置的位流设计方法。优化变量为缝翼相对于主翼的缝隙(Gap)、覆盖量(Ouerlap)和偏角δ_s,目标函数为主翼上的压力峰值。应用高阶面元法计算多段翼型压强分布。用Powell优化法使主翼上压力峰值减至最小,以延迟多段翼型的失速,增大最大升力系数。本方法已用于计算NACA64A010两段和四段翼型以及Foster三段翼型,所得结果与实验数据和位流/边界层耦合设计法的结果有很好的一致性。     

基于欧拉方程的二维振荡机翼非定常气动设计反命题方法（Ⅱ）

建立了二维振荡机翼的非定常反问题的有限差分解法,并通过若干算例对求解方法进行了验证.首先根据基于欧拉方程的非定常气动反问题数学模型确定了具体求解方案,接着重点讨论了壁面非定常反问题边界条件和远场非定常无反射边界条件的处理方法,编写了相应的计算程序,并以绕1/4弦长点做小幅俯仰周期振动的NACA-65-0012翼型为初始机翼,在不同目标压强的条件下进行了气动反问题计算得到新翼型,它们的压强分布与预设目标相吻合,这表明本文的反问题求解方案与计算方法都是可行的.