固定转捩在改善振荡来流下低雷诺数翼型气动性能中的应用

低雷诺数下,翼型表面易于出现的分离转捩现象会降低翼型的气动性能,采用数值计算方法探讨了固定转捩在改善低雷诺数下翼型气动性能中的应用。定常来流下,雷诺数为4×104时,在翼型E387上表面自由转捩位置之前的一定位置处使翼型固定转捩,计算结果表明分离区减小了,升力系数和升阻比明显提高。振荡来流下,雷诺数在8×104的50%范围内变化时,分别计算了NACA4412在0°迎角和E387在3°迎角时自由转捩及在60%弦长处固定转捩的升阻特性,结果显示自由转捩的翼型当雷诺数减小时,性能急剧恶化,而采用固定转捩的翼型受其影响要小得多,具有更加稳定的气动性能。     

不可压缩势流绕多段翼型流动的快速迭代解法

利用保角变换计算不可压缩势流绕翼型的流动,具有计算速度快和要求计算机内存小的优点,但主要用于单段翼型问题,而对于多段翼型,尤其是三段和三段以上的翼型则很难应用。本文在单段翼型的保角变换的基础上,提出了一种利用迭代算法计算任意多段翼型的不可压缩势流的绕流问题的方法。本文方法和目前常用的边界元法相比,具有计算速度快的优点。本方法还可进一步推广到考虑无粘流/边界层的迭代计算中。     

不可压缩实际粘性流体绕流孤立翼型的数值研究

本文提出了基于微分-积分方程组求解N-S方程的有限差分法求解不可压缩实际粘性流体绕孤立翼型流动。计算中采用了K-ε湍流模型结合考虑逆压梯度的壁面函数法,并导出任意非正交曲线坐标系下的压力修正方法和适用于大流通量的混合格式。分别计算了层流和湍流情况下带分离翼型绕流,与流动显示及实验数据比较表明,在预测分离点位置,翼型表面压力、平均速度分布、雷诺切应力uv方面,两者吻合,结果是满意的。     

用粘流－无粘流相互作用方法计算翼型的分离和失速

本文用势流－边界层相互作用方法计算低速翼型的分离和失速。势流用对称面元法。边界层用改进的滞后掺混法，考虑了高阶项影响，适用于计算分离。文中对粘流－无粘流耦合方法作了改进。改进的半反－局部联立耦合方法，考虑了相邻点之间的作用，收敛性较好。计算了ＮＡＣＡ４４１２翼型在不同迎角下的压力分布和气动力。计算结果与实验符合良好。算例表明，翼型高升力状态计算必须包括尾流的作用，也应当对势流计算压力与实际压力之间的差别进行修正。     

M_∞=1及其附近翼型绕流的计算结果

本文全速势方程的有限差分数值计算结果,给出了0°和2°攻角下NACA0012翼型当地马赫数分布冻结时的自由流马赫数范围,以及M_∞稍大于1变至1.30时前方脱体激波的变化位置。     

NACA4412翼型低速绕流数值计算中湍流模型对比

使用Spalart-Allmaras(S-A)、SSTk-ω、Gao-Yong 3种湍流模型对NACA4412翼型低速绕流进行了数值计算,研究了尾迹流动松弛效应。对流项采用Roe格式离散,扩散项采用二阶中心格式离散,离散后的控制方程用多步Runge-Kutta显示时间推进法求解。计算中对翼型尾缘流松弛效应进行了分析,比较了翼型表面压力系数、速度剖面、雷诺应力等的分布,3种湍流模型总体上能够较好地模拟NACA4412翼型低速绕流。SSTk-ω模型对流动细节及升力系数计算最好,Gao-Yong模型对翼型平均速度剖面及雷诺剪切应力分布计算最准确。     

跨音速翼型非守恒全位势流AF3高效算法的应用

本文利用Baker的AF3高效迭代算法,改进了Jameson提出的计算跨音速翼型非守恒全位势流的有限差分方法,研制了任意翼型跨音速分析的高效算法程序。许多试验算例表明,本算法,收敛速度快,且计算结果与Jameson的原方法非常一致。     

等速上仰翼型绕流结构的数值模拟

通过数值求解涡量-流函数形式的N-S方程研究等速上仰翼型的分离流动结构。用ADI方法解涡量方程;用Poisson方程直接法解流函数方程。在计算得到的流场中除可看到实验中观察到的多涡结构外,还发现一些实验中未观察到的现象。将分离流动结构的计算结果与翼型气动力变化相联系,并结合涡动力学理论,定性地解释了上仰翼型产生高升力的原因。     

不同翼型厚度和来流马赫数下的桨涡干扰噪声分析

桨涡干扰噪声是直升机气动噪声主要组成之一，为了正确预测和降低直升机噪声，必须开展气动噪声相关物理参数研究。在对声场进行计算流体力学（CFD）直接数值模拟的基础上，分析了不同厚度和来流马赫数下二维平行桨涡干扰噪声传播特性和声源位置，分析了翼型厚度和来流马赫数对桨涡干扰噪声的影响，并得到了可压缩情况下远场声压预测公式。研究表明，低马赫数下，翼型厚度对噪声指向性影响不大，高马赫数下，翼型厚度对噪声指向性影响程度增大；噪声强弱主要随来流马赫数变化，翼型厚度对其影响较小；翼型厚度和来流马赫数变化不会改变声源点位置。开展不同翼型厚度和来流马赫数下的桨涡干扰噪声分析可以为进一步了解并控制直升机桨涡干扰噪声提供一定的参考。     

旋涡对翼型的气动干扰

用数值模拟方法研究了旋涡对翼型的气动干扰。出发方程是N-S方程。数值格式为Beam-Warming格式。结果表明,旋涡位置改变,其对翼型的气动(激波和升力)干扰随之有显著变化;涡的有利干扰可使翼面上流动分离受到一定程度的控制,导致翼型升力增大;涡的位置和强度的合理组合能有效地提高升力。计算还表明,单一旋涡扫过翼型时平均升力将提高50％左右。     

仿生学翼型尾缘锯齿降噪机理

采用大涡模拟与声类比的方法研究了尾缘锯齿对翼型自噪声的影响。以SD2030翼型为研究对象,设计的尾缘锯齿幅值为10%弦长,周期为4%弦长。模拟了来流速度为31 m/s、0° 攻角下直尾缘翼型与锯齿尾缘翼型的流场,对应的基于弦长的雷诺数约为310 000。详细分析了尾缘锯齿对翼型尾缘湍流流场的影响,并通过FW-H方程计算大涡模拟提取的声源项,得到直尾缘翼型与锯齿尾缘翼型的声场。研究发现,锯齿尾缘可以明显降低翼型中低频范围内的噪声,在4 000 Hz以下,窄带噪声最多可降低约16 dB。但尾缘锯齿对翼型气动性能有着不利影响。进一步研究表明,该状态下翼型噪声主要由层流边界层引起的涡脱落噪声主导,尾缘锯齿可以抑制层流边界层引起的涡脱落现象,降低翼型升力脉动与尾缘附近的表面压力脉动,减弱尾缘处的低频湍流脉动与涡量,并有效降低尾缘附近涡的展向相关性,这些因素的综合作用使得翼型自噪声降低。     

降雨条件下翼型的气动特性

采用计算流体力学方法对降雨条件下翼型的气动特性进行了研究.通过求解定常可压的Navier-Stokes方程来计算流场,采用拉格朗日法对流场中的雨滴轨迹进行了模拟跟踪,得到了翼型表面的雨滴收集率;建立了水膜层数学模型,并假定水膜层沿着翼型表面的法向方向增长,预测了降雨引起的翼型外形的变化,并得到了水膜层表面的粗糙度高度;采用k-ωSST(shear-stress transport)两方程湍流模型结合增强的壁面函数,研究了降雨对翼型气动特性的影响.结果表明,在低迎角范围内,降雨对翼型的升力系数和阻力系数影响很小;当达到失速迎角后,降雨会引起上翼面边界层气流的提前分离,造成翼型气动性能的严重损失.      

锯齿型翼型尾缘噪声控制实验研究

翼型湍流边界层与尾缘相互作用产生的尾缘噪声是翼型自噪声的最主要分量,多年来研究者们已经在理论、数值和实验方面开展了多方面研究,但有关翼型尾缘噪声的产生机理和抑制方法仍有待深化与发展.本文利用具有全消声环境的低速开口风洞研究了采用锯齿型翼型尾缘来控制翼型噪声的方法,重点研究了不同攻角情况下不同锯齿形对翼型远声场气动噪声的影响以及翼型表面压力的影响.实验结果表明,翼型尾缘附加锯齿是一种可行的降噪方案,尤其对中低频段的远场气动噪声有比较明显的降低效果;而且,降噪效果与锯齿的齿数和齿间倒角有关.附加锯齿对翼型壁面动态载荷的影响较小,基本不影响翼型的气动性能.     

涡襟翼在不同雷诺数下的控制分离特性研究

受鸟类抬起羽毛控制分离流的启发,涡襟翼成为翼型大迎角分离流的控制措施之一。采用数值模拟方法研究不同雷诺数下涡襟翼在控制翼型大迎角分离流动时的气动特性及其物理机制。结果表明：涡襟翼在低雷诺数下能够极大地改善翼型的大迎角升力特性,其物理机理是涡襟翼将翼型主分离涡的涡心位置控制在离翼型更近的区域,且涡心位置的涡量得到大幅提升,使得涡心附近的低压特性影响到翼型上表面,而且涡襟翼能够将翼型上方前区的低压与下游的高压隔开;但是在高雷诺数（对应常规飞机雷诺数）下涡襟翼改善翼型大迎角气动特性的效果远不如低雷诺数情况,由此解释了为什么鸟类能够通过羽毛抬起提高升力特性,而常规飞机的涡襟翼只能作为阻力板使用的原因。     

跨声速翼型上激波／边界层干扰的自适应控制计算

在激波区使用自适应壁对跨声速翼型上的激波/边界层干扰进行控制,可改变机翼的气动性能。这种被动控制可通过在翼型的激波区开一凹腔,其上覆盖一弹性橡胶膜柔壁来实现。本文给出用N-S方程数值模拟这一自适应控制翼型的跨声速粘性绕流,提出了一个适用于本特殊情况(物面边界局部地区在求解过程中有变化)的处理办法。并探讨了自适应柔壁对当代跨声速翼型绕流的影响。     

A boundary surrogate model for micro/nano grooved surface structure applied in turbulence flow control over airfoil

The application of grooved surface structure is an emerging and effective means in turbulence flow control. However, for a realistic configuration, the global flow field described directly by simple application of massive grids makes it unfeasible to simulate. In this paper, a boundary surrogate model reproducing the effect of microscopic near-wall region is proposed to improve computational efficiency. The surrogate model trained with Lattice Boltzmann Method(LBM) considering the rarefied effec...     

任意翼型有小分离气泡跨音速绕流的迭代算法

 本文给出一种带小分离气泡的任意翼型粘性跨音速绕流的计算方法,采用有粘-无粘干扰迭代的概念。无粘流的全速势方程用AF差分格式在保角变换法生成的计算网格中求解,粘流附面层方程用C-S盒式法求解,用逆算法消除分离点处的奇性。本文对Ma_∞=0.8,Re_∞=2×10~6,迎角α=3.5°和4°的NACA64A010翼型粘性绕流进行了计算,结果与实验相比较,吻合良好。     

一种新的边界处理方法在笛卡尔网格中的应用

采用自适应笛卡尔网格方法求解Euler方程,固壁边界条件通过一种离散强迫浸入边界方法引入,边界切割网格与流场网格保持一致,从而有效解决了小切割网格单元的时间步长限制问题.针对自适应笛卡尔网格的特点提出了新的虚拟点反射方法和插值方法,解决了数值吸力峰值问题.对RAE2822翼型、双NACA0012翼型和Suddhoo三元翼型的流动状况进行了无黏数值模拟,并与现有的理论解、单域以及分区结构网格解进行了对比.结果表明:该离散强迫浸入边界方法有效地引进了固壁边界条件,结合自适应笛卡尔网格技术能够准确模拟复杂几何外形下的流动状况.      

基于特征正交分解的跨声速流场重构和翼型反设计方法研究

在二维流场的重构问题中应用特征正交分解(Proper Orthogonal Decomposition,POD)数据处理方法。利用CFD技术计算得到的流场快照对气动力模型进行降阶,然后利用基于POD的降阶模型(Reduced Order Model,ROM)对所需的流场参数进行重构,在快照范围内可以得到高精度的结果,且具有一定的外插能力。在翼型反设计问题中该方法仍然是成功的,通过修正快照向量,利用基于POD降阶模型的数据重构方法,由已知的翼型表面压力分布通过反设计就能够高效精确地得到对应的最优翼型形状,这极大地简化了翼型反设计问题。本文分别在跨声速范围对RAE 2822翼型的流场重构和Korn翼型及NACA 63212翼型的反设计进行了验证,证明了基于POD的流场重构和翼型反设计方法的有效性和高效性。     

基于Lattice Bohzmann方法的翼型绕流数值模拟

将通过物理平面的结构化C型贴体网格变换到计算平面的均匀正交网格,引入非均匀网格插值方法,采用基于Lattice-Boltzmann方程的D2Q9模型,对低雷诺数下的NACA0012翼型绕流进行了数值模拟,对五个不同位置的边界层速度型与基于有限体积法的CFL3D程序的结果进行了比较.结果表明,两种方法的计算结果一致.改进后的Lattice-Boltzmann方法适用于曲边边界,非均匀网格,同时对计算速度影响不大.Lattice-Boltzmann方法计算较简单,因而用LB方法来模拟翼型绕流是一种新的选择,并可应用于更复杂的低雷诺数流动中.