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

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

直升机新型旋翼翼型气动特性与布局分析

利用结构网格计算流体力学(Computational fluid dynamics ,CFD)的翼型气动特性分析方法开展旋翼翼型气动特性计算。通过RAE282，NACA0012，OA212，OA207等翼型压力分布、升力和阻力等特性计算结果与试验结果的对比分析，验证了计算方法的准确性，并进一步完成了HF系列旋翼翼型的气动特性计算。基于翼型的气动特性，采用时间步进自由尾迹的旋翼气动性能分析方法开展旋翼桨叶翼型的气动布局优化设计，对悬停和前飞条件下的旋翼开展计算分析，得到两种条件下的旋翼气动特性。而后通过本文建立的优化方法开展旋翼翼型布局优化设计。     

连续变弯度翼型动态气动特性数值模拟

针对翼型后缘连续变弯度运动中后缘边界精确数值模拟的问题，提出了基于二维多项式的时空曲面拟合方法，实现了对后缘边界时空位置的精确模拟。在此基础上，基于OpenFOAM发展了翼型后缘连续变弯度与大幅度俯仰运动耦合的运动边界数值模拟，并计算了翼型耦合运动的气动力，讨论了后缘线性/非线性变形对翼型大迎角动态气动特性的影响规律。结果表明:后缘运动对翼型俯仰运动的升阻特性影响显著，特别在翼型大幅度俯仰时后缘非线性变形对升阻特性改善效果比线性变形大6%~10%。同时还研究了翼型俯仰与后缘变形运动相位差对气动特性的影响。特别地，当相位差为180°时，后缘运动使动态失速时的最大升力提高50.3%，平均升力提高34.6%；当2种运动相位差为0°时，后缘运动使动态失速时的最大阻力降低39.7%，平均阻力降低30.2%，最大升阻比提高22.3%，平均升阻比提高16.8%；同时，翼型在动态俯仰过程中出现负阻力现象，对产生负阻力的原因进行了分析。这些结果可用于指导连续变弯度后缘控制律的设计。      

基于温敏漆的边界层转捩测量技术研究

以自然层流翼型RAE5243模型为研究对象，在0.6m跨超声速风洞进行温敏漆（Temperature Sensitive Paint，TSP）转捩测量技术研究，在Ma0.73和Ma0.75条件下开展了模型基本外形和鼓包外形的转捩测量试验。针对缺乏定量分析手段的问题，提出基于温度梯度的转捩位置自动判定算法，包括图像预处理、转捩点定位与筛选和转捩位置计算3个步骤。模型温度分布及转捩测量结果表明:重复性试验结果偏差较小，验证了转捩测量结果的可靠性；相同马赫数条件下，鼓包外形转捩位置相对基本外形向后缘移动；相同外形条件下，Ma0.75的转捩位置相对Ma0.73向后缘移动。TSP试验结果与CFD计算结果吻合较好，变化趋势一致，检验了数值模拟方法的有效性。     

热/力耦合作用下基于应力分析的冰破坏准则

根据飞机热除冰的物理过程,考虑外部空气动力和蒙皮表面加热的作用,建立了NACA 0012翼型前缘冰层应力计算模型。采用有限元方法和平面三角形单元对控制方程组进行了求解,获得了外部空气动力和蒙皮表面加热对冰层黏附界面应力的影响规律。研究表明:蒙皮表面不加热时,来流速度影响了黏附界面应力的强度,来流攻角影响了黏附界面应力的分布,冰-蒙皮间黏附界面切应力最大值随来流速度呈近似线性增大趋势,但外部空气动力很难造成冰层破坏。蒙皮表面加热时,冰-蒙皮间黏附界面的耦合应力和冰层内部的主应力随着热流密度的增大而增大,很容易超过剪切强度,这是造成冰破坏的关键因素。耦合冰-蒙皮剪切强度随界面温度的变化关系,初步建立了基于应力分析和热/力耦合作用的冰破坏判断准则。外部空气动力产生的界面应力和蒙皮表面加热产生的界面热应力之和，必须大于与蒙皮表面温度相关的剪切强度，则冰层发生破坏，破坏位置是耦合应力超过剪切强度的区域。      

Generalized Kutta–Joukowski theorem for multi-vortex and multi-airfoil flow with vortex production——A general model

By using a special momentum approach and with the help of interchange between singularity velocity and induced flow velocity, we derive in a physical way explicit force formulas for twodimensional inviscid flow involving multiple bound and free vortices, multiple airfoils, and vortex production. These force formulas hold individually for each airfoil thus allowing for force decomposition, and the contributions to forces from singularities(such as bound and image vortices,sources, and doublets) and bodies out of an airfoil are related to their induced velocities at the locations of singularities inside this airfoil. The force contribution due to vortex production is related to the vortex production rate and the distance between each pair of vortices in production, thus frameindependent. The formulas are validated against a number of standard problems. These force formulas, which generalize the classic Kutta–Joukowski theorem(for a single bound vortex) and the recent generalized Lagally theorem(for problems without a bound vortex and vortex production) to more general cases, can be used to identify or understand the roles of outside vortices and bodies on the forces of the actual body, optimize arrangement of outside vortices and bodies for force enhancement or reduction, and derive analytical force formulas once the flow field is given or known.     

基于合成射流的旋翼翼型动态失速控制研究

针对直升机旋翼工作环境下来流速度和迎角(Angle of attack,AoA)耦合引起的动态失速问题,建立了基于合成射流的旋翼动态失速控制的数值模拟方法。采用运动嵌套网格方法,通过对翼型的平移和旋转实现变来流速度-变迎角的耦合。以积分形式的雷诺平均N-S方程为主控方程,空间离散使用Roe格式,时间离散为隐式LU-SGS方法,以OA209翼型为研究对象,在翼型上表面放置合成射流激振器,开展了射流位置、动量系数、无量纲频率以及偏角等参数对轻度失速、深度失速下翼型动态失速控制的研究。研究发现,轻度失速下,射流位置靠近气流分离点时(20%c附近,c为翼型弦长),对逆压梯度引起的轻度失速控制效果最佳。深度失速下气流分离点虽在5%c之前,但射流位于前缘分离泡后端(10%c附近)时控制效果较好。大迎角需要较大的动量系数才能有效控制。射流频率对涡结构的尺寸和数量会产生一定影响,能改变气动特性波动幅度。较小的射流偏角对轻度失速的控制更有效,而深度失速则需要较大的偏角。     

射频介质阻挡放电改善NACA 0015翼型气动性能的实验

介质阻挡放电(DBD)均匀稳定、易于敷设,是机翼/翼型等离子体流动控制(PFC)中最常用的激励方式。射频介质阻挡放电激励频率高、放电功率大,且能在流场中产生明显的加热,应用潜力大。采用射频电源驱动DBD激励器产生等离子体,分析放电的体积力、热特性和诱导流场特性,开展了射频介质阻挡放电改善NACA 0015翼型气动性能的实验,研究了占空比、调制频率、载波频率和电源功率等参数对流动控制效果的影响规律。结果表明:射频等离子体激励的体积力效应随激励电压的增大而增加;射频等离子体激励产生的热量在诱导的流场中进行传导,加速流场;当来流速度为20m/s,Re=3.36×10~5时,在翼型前缘施加激励,使翼型临界失速迎角推迟1°,最大升力系数增大6.43%,且在过失速迎角下仍具有流动控制效果,使升力下降变缓;调制频率越大,控制效果越好;存在最佳占空比、载波频率和功率,占空比对流场控制效果的影响最显著,最佳占空比、载波频率和功率分别为20%,460kHz和50W。射频等离子体激励以体积力效应、热效应和诱导壁面射流改善失速流场,使得NACA0015翼型气动性能极大改善,流动分离得到有效控制。     

前缘参数对翼型气动-隐身特性敏感性分析

针对翼面部件对隐身飞机气动和隐身性能影响更为突出的问题,提出了一种翼型前缘参数化修形方法,采用高精度气动和电磁数值计算方法,进行了前缘形状对翼型气动-隐身特性的敏感性分析。研究结果表明:翼型前缘半径增加,其最大升力系数、失速迎角、最大升阻比显著增加,但其前向RCS均值也增加明显,设计时需进行综合权衡;在不影响结构和装载的情况下,下表面偏角尽量小,修形长度取0.15c,对翼型的低、高速气动和隐身特性都有利。     

翼型低雷诺数层流分离现象随雷诺数的演化特征

层流分离现象是翼型低雷诺数条件下出现的典型流场特征。层流分离流动中包含流动分离、转捩、再附等非定常流动结构，层流分离流动的形成与演化会对翼型气动特性产生恶化作用。采用大涡模拟（LES）方法对低雷诺数范围内不同雷诺数下的翼型层流分离流动开展精细数值模拟，研究了雷诺数对翼型气动特性的影响规律及作用机理。LES方法采用隐式亚格子模型，基于结构化拼接网格，对流项离散和时间推进方法分别采用AUSM⁺格式以及双时间步方法。验证算例计算结果表明数值模拟方法的正确性及可靠性，雷诺数对翼型气动特性具有显著影响。随雷诺数降低，时均分离泡外形增大、位置后移，平均阻力系数增大，特别是在较低雷诺数下，翼型升阻力系数随时间出现振荡现象。进一步研究表明，造成不同时均分离泡形态和气动特性的原因在于翼型上表面分离剪切层的失稳与转捩特征。随雷诺数降低，流动黏性增大，导致分离剪切层速度梯度减小，流动发生转捩及再附位置后移，直至翼型表面不再发生转捩和再附。