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

航空科学技术Aeronautical Science&TechnologyJun.25 2023 Vol.34 No.06 11-19不同翼型厚度和来流马赫数下的桨涡干扰噪声分析相倩1，余启梁2，王钢林1，刘勇2，叶海剑21.中国航空研究院，北京 1000292.南昌航空大学，江西 南昌 330063摘 要：桨涡干扰噪声是直升机气动噪声主要组成之一，为了正确预测和降低直升机噪声，必须开展气动噪声相关物理参数研究。在对声场进行计算流体力学（CFD）直接数值模拟的基础上，分析了不同厚度和来流马赫数下二维平行桨涡干扰噪声传播特性和声源位置，分析了翼型厚度和来流马赫数对桨涡干扰噪声的影响，并得到了可压缩情况下远场声压预测公式。研究表明，低马赫数下，翼型厚度对噪声指向性影响不大，高马赫数下，翼型厚度对噪声指向性影响程度增大；噪声强弱主要随来流马赫数变化，翼型厚度对其影响较小；翼型厚度和来流马赫数变化不会改变声源点位置。开展不同翼型厚度和来流马赫数下的桨涡干扰噪声分析可以为进一步了解并控制直升机桨涡干扰噪声提供一定的参考。关键词：翼型厚度；来流马赫数；二维平行BVI；气动噪声；直接法中图分类号：V211 文献标识码：A DOI：10.19452/j.issn1007-5453.2023.06.002 在直升机低速下降、机动等飞行过程中，前行桨叶产生的桨尖涡脱落后会靠近或通过后行桨叶，形成桨涡干扰（BVI）现象。当桨尖涡轴线与桨叶平行时，可以简化为二维平行BVI。BVI噪声一旦出现，会显著增大直升机的总体噪声水平，带来严重的噪声污染1。因此，对影响BVI气动噪声的物理参数的研究，一直是预测和降低直升机噪声的前提2。先前，在理论和试验方面研究者们已开展了许多与影响BVI气动噪声的物理参数相关的研究工作3-7。孟晓伟等3通过建立悬停状态直升机倾斜式尾桨诱导速度的方法有效地对尾桨BVI进行了预测。赵俊等4基于计算流体力学（CFD）方法对多旋翼飞行器气动力和噪声特性进行了分析。H.Y.Yung等5通过理论计算、全机试验和风洞试验等对BVI噪声声压级和飞行状态影响进行了研究，但对于影响BVI的物理参数、噪声波的形成和传播过程等方面没有深入研究。F.Caradonna等6在风洞中用旋转机翼研究了BVI噪声，分析了机翼与涡流相互作用，比较了Kirchhoff方法和FW-H方法等各种预测BVI噪声的方法。远场压力的计算结果和测量结果对比表明，当涡特征（涡核位置、涡强度和涡核半径）已知时，可以较准确地预测BVI噪声。C.Kitaplioglu 等7为减少旋翼转动时自身尾流的影响，利用外部产生的涡流与表面装有压力传感器的静止旋翼相互作用，并通过固定的近场麦克风和可移动的远场麦克风阵列收集声学数据。试验研究了翼型与涡干扰垂直间距、涡核强度以及叶尖马赫数等几个影响BVI的参数，该试验方法虽准确可靠，但试验周期长且成本高。旋翼BVI现象复杂且随机性强，在旋翼各个角度和展向位置上可能同时发生BVI现象，通过试验方式测量噪声存在局限。随着计算机的发展，CFD方法逐渐成为直升机流场分析最有效的方法之一8。J.M.Chen等9通过在涡流发生器的下游放置NACA0012翼型，对二维BVI进行了试验研究，研究结果表明，当旋涡与非升力翼型相互作用时，它在前缘附近产生相当剧烈的局部表面压力变化。对于迎角参数的影响，该研究者先后试验了迎角为5和10的工况，试验结果显示，相比于5迎角，10迎角下翼型上的瞬态升力幅度更大，瞬态升力的幅值也随着涡流发生器和机翼收稿日期：2023-03-07；退修日期：2023-03-30；录用日期：2023-04-25基金项目：国家自然科学基金（11962018，12262023）引用格式：Xiang Qian，Yu Qiliang，Wang Ganglin，et al.Analysis on rotor vortex interference noise under different airfoil thickness and free-stream Mach numberJ.Aeronautical Science&Technology，2023，34（06）：11-19.相倩，余启梁，王钢林，等.不同翼型厚度和来流马赫数下的桨涡干扰噪声分析 J.航空科学技术，2023，34（06）：11-19.航空科学技术Jun.25 2023 Vol.34 No.06之间垂直间隔的减小而增大。I.Marcel等10-12使用大涡模拟（LES）方法对不同迎角、不同涡特征以及翼型与涡不同干扰垂直间距下的BVI流场进行了数值模拟，减小了数值耗散，并捕捉到了干扰时的流场信息，研究结果表明，气动系数的大小随着干扰垂直间距的增加和涡核尺寸的减小而减小，LES方法为获得非定常壁面压力场提供了一种很有前景的工具。在国内，史勇杰等13通过Euler方程分析BVI流场的特征，将自由涡在流场的诱导速度等效成网格畸变运动，并入真实网格运动中保持涡核结构不变。描述了干扰流场的升力变化特性，然后分析了涡核强度、翼型与涡干扰垂直间距和干扰角度等参数对BVI流场特性的影响。乔渭阳等14采用有限翼展桨叶模型计算叶片表面的非定常力，对涡强度、涡核半径、翼型与涡干扰的垂直间距和夹角、马赫数等影响参数进行实例计算。声场的计算使用了声比拟法。结果表明，该方法得出的BVI噪声声压数据与试验结果吻合良好。目前，对于影响BVI气动噪声的物理参数研究主要围绕流场分析，而物理参数对BVI气动噪声产生和传播方面的影响还需要进一步计算和分析。BVI气动噪声生成过程是一个强烈的非线性过程15，对于气动噪声的数值计算，大多采用声比拟方法预测远场噪声；而LBM-LES相结合的直接计算气动噪声方法，不仅能准确预测远场噪声，还能捕捉气动噪声的生成、传播和衰减实时变化过程，有利于认识不同翼型厚度和来流马赫数对BVI脉冲气动噪声的影响问题，所以本文采用LBM-LES气动噪声直接计算方法，对不同翼型厚度和来流马赫数下，二维平行BVI产生的脉冲非线性气动噪声进行直接计算。1 数值方法及边界条件1.1 直接法直接法计算气动噪声具有同时得到精确流场和声场的优点，可用于对气动噪声进行传播过程、产生机理等研究，但需要采用高精度的时空离散格式和无反射边界条件两个关键技术来实施。1.2 LBM-LES方法格子玻耳兹曼（LBM）方法可高效求解非定常低速不可压的湍流问题16，多松弛时间（MRT）-LBM方法具有四阶空间精度和二阶时间精度，时间离散采用显式格式，具有数值稳定性高、低耗散、精度高等特性17。在相同的耗散情况下，该方法计算效率高，因此逐渐发展出了基于LBM直接法计算气动噪声。本文数值模拟采用D2Q9格子离散速度模型18，离散速度为e=()0,0,=0,()cos,sinc,=()-1 2,=1,2,3,42()cos,sinc,=()-1 2+4,=5,6,7,8（1）式中，e为离散速度，为角度，为节点，c为格子声速。LES方法是通过亚格子模型求解小尺度涡，对含能较高的大尺度涡则直接求解。BVI噪声属于中低频窄带噪声，具有含能高且尺度大的特点，和高雷诺数低马赫数湍流噪声相比，对计算精度、时间尺度、空间尺度、计算机及耗散性等的要求可以降低，因而比较容易实施，这使采用LES方法直接计算BVI近场和远场噪声成为可能。1.3 声场计算方法因为流场与声场的基本方程是相同的，所以可通过求解LBM方程直接得出流场和声场的统一解。通过LBM-LES方法可直接得出流体宏观压力p和流场密度，p=c2s，cs为流体声速。对声场进行直接计算时，格子声速等于流体声速。声压p为介质中压力的扰动值19，考虑到流体的静密度0和静压0c2s，p可以表示为p=c2s-0c2s（2）1.4 网格分布及边界条件由于BVI干扰主要作用于翼型附近区域，所以采用多尺度模式来划分计算网格。将计算区域划分成壁面附近、近场、远场和声波吸收4个区域。图1所示为区域划分及网格分布示意图，图1（a）翼型附近区域放大为图1（b）。计算域大小为80C 80C，其中C为翼型弦长。上述4个边界均采用求解局部一维无黏方程20-21实现无反射边界条件，避免了边界处波向声场内反射，远场外有网格尺度较大的声波吸收区域，LES方法可过滤掉来自边界的小尺度伪反射波。2 计算结果与分析本 文 算 例 分 别 采 用 厚 度 为 12%、15%、18%的NACA0012、NACA0015、NACA0018翼型，翼型迎角=0，翼型弦长C=0.46m；来流马赫数Ma分别为0.14、0.4、0.6；旋涡模型为Lamb-Oseen22，其环量=-0.16UC，负号表示涡的旋转方向为顺时针方向，涡核半径rc=0.05C，涡的初始位置(xv,yv)=(-4C,0)，其中xv为涡相对于翼型前缘的位置，12相倩 等：不同翼型厚度和来流马赫数下的桨涡干扰噪声分析yv为涡与翼型之间的垂直距离。2.1 网格无关性验证首先进行网格无关性分析，以保证计算结果的可靠性。基于BVI引起的脉冲声压峰值建立网格数目影响分析。表1为网格数目和尺寸，N表示网格数目。图2表示三种不同网格数目在监测点(0,20C)处的声压时间历程。由图2可以发现，N=107万和N=220万之间的脉冲声压曲线更为接近，而N=70万的脉冲声压曲线与另两条脉冲声压曲线存在差异；随着网格数目的增加，N=220万和N=107万的声压峰值差异更小。考虑到节省计算资源以及计算精度，本文采用N=107万网格进行数值计算。图2(0,20C)处监测点声压时间历程Fig.2Time image of sound pressure at the monitoring point(0,20C)2.2 噪声传播通过直接法计算气动噪声可以直接获得声场信息，一般用胀量=V显示声波23：胀量值大小可表示声波强弱，胀量值最大值连线可表示辐射方向。图3所示胀量云图表示不同翼型厚度和来流马赫数下桨涡干扰脉冲噪声声波向远场传播的过程，其中相同来流马赫数下脉冲声波截取时刻相同，即图3（a）图3（c）为相同时刻胀量云图，图3（d）图3（f）为相同时刻胀量云图，图3（g）图3（i）为相同时刻胀量云图。以图3（a）为例，t1t4代表脉冲噪声声波传播的不同时刻；角度表示声波上下辐射角度值（辐射角以翼型下游对应的方向为0，逆时针为正）；正负数字表示t1时刻声波沿辐射角方向胀量值，即t1时刻胀量最大值。由图3可知，不存在因为网格分区和远场边界产生与脉冲噪声同量级强度伪反射波；从上下远场胀量的变化规律可知，上下远场噪声幅值大小相近，但相位相反，胀量图呈8字形，说明声场的辐射特性符合偶极子声源的辐射特性；上远场辐射角绝对值（简称上辐射角）和下远场辐射角绝对值（简称下辐射角）表示脉冲噪声声源指向性；上远场胀量绝对值（简称上胀量值）均略小于下远场胀量绝对值（简称下胀量值），说明来流涡为顺时针旋转时，上源场噪声比下源场噪声稍弱，是旋涡的下洗导致；相同马赫数下比较，随着厚度的增加，相同时刻相同位置处胀量绝对值减小，说明声波随厚度增加而变弱，但变化幅度不大。如图3（a）图3（b）所示，在来流马赫数Ma为0.14时，表1 网格数目和尺寸Table 1 Number and size of meshesN=70万N=107万N=220万壁面附近区域/m3.125010-31.562510-31.562510-3近场区域/m3.125010-33.125010-31.562510-3远场区域/m0.0500.0500.02540202040402000.100.050.050.100.100.10 xCxCyCyC20400NACAaC图1区域划分及网格分布示意图Fig.1Schematic diagram of zoning and grid distribution0.40.2pPats0.20.40.600.120.150.180.21N70N107N22013航空科学技术Jun.25 2023 Vol.34 No.06t1时刻