Wind tunnel test research on the characteristics of rotor blade-vortex interaction noise
-
摘要: 在中国航空工业空气动力研究院FL–10风洞中开展了旋翼桨–涡干扰噪声传播特性试验,对BO–105 主旋翼40% 缩比模型中等前飞速度爬升、平飞、斜下降状态的气动噪声进行了测量。首先采用Heyson洞壁干扰修正方法确定风洞试验时的旋翼下滑角,通过气流内测量阵列移动获得了桨盘平面下方完整的噪声辐射场,然后对不同飞行状态下的桨–涡干扰噪声传播特性进行了分析,得到了典型状态的声压–时间历程、频谱和声压级云图。结果表明:旋翼斜下降飞行状态出现了明显的桨–涡干扰噪声,干扰较强时桨叶前行侧和后行侧都会产生桨–涡干扰噪声,且其传播具有明显的方向性,即前行侧指向桨盘上游和桨盘下方,后行侧指向桨盘下游。
-
关键词:
- 直升机 /
- 风洞试验 /
- 桨–涡干扰;气动噪声 /
- 旋翼 /
- 斜下降飞行
Abstract: The propagation characteristics of the blade–vortex interaction (BVI) noise were tested in the FL–10 wind tunnel of AVIC Aerodynamics Research Institute. The aerodynamic noise of climb, level flight and descent at medium forward speed was measured on a 40% scale model of BO−105 main rotor. Firstly, the “Heyson” wall interference correction method was used to determine the descent angle of the rotor in the wind tunnel, and the complete noise radiation field under the rotor tip–path–plane was obtained through the movement of the measurement array in the airflow. Furthermore, the BVI noise characteristics of the rotor under different flight condition were studied, and the sound pressure time history, spectrum and sound pressure level contour were given. The results indicate that the BVI phenomena occur on both the advancing and retreating side under the descent flight condition. The noise has strong directivity, and radiates toward upstream under the rotor disk toward the advancing side and the downstream on the retreating side.-
Keywords:
- helicopter /
- wind tunnel test /
- blade–vortex interaction /
- aerodynamic noise; rotor /
- descent flight
-
0 引言
超空泡流动现象常见于水中高速鱼雷、潜水艇发射战略战术导弹和利用高速射弹消灭水雷等军事技术中,因此一直是国内外的热门研究课题之一[1-2],本课题组对此也进行了研究[3-5]。超空泡带迎角地穿过自由面问题源于从潜水艇鱼雷发射管发射战斧式巡航导弹,为此美国海军进行了大量的实验[6]。我国在这方面的研究还不够充分,无论是实验还是数值计算。附着在潜射导弹上的空泡一般是自然空化空泡(Natural Cavitation),空泡内部压力约为水的饱和蒸气压(2.35×103Pa)。导弹出水时,空泡与压力为1个标准大气压(105Pa)的空气接触,因此空气必然进入空泡,使其成为充气空泡(Ventilated Cavitation)。然而必须明确的是,与用人工方法给空泡充气的过程不同(诸如高速鱼雷),在出水时空泡获得的充气仍然是在自然状态下进行的。
胡影影等人[7]使用改进的Youngs算法对航行体在自由面附近出水进行了数值模拟,着重研究了韦伯数We和弗汝德数Fr在出水运动过程中对液面的影响。黄海龙等人[8]通过划分结构网格,对圆盘空化器进行了三维数值仿真,模拟中应用了k-ε湍流模型求解N-S方程,通过改变不同的空化数大小,得到相应的计算参数,然后将计算参数的结果与前人的经验公式相对比。王一伟等人[9]通过将结构力学和计算流体力学2种方法相结合,建立了弹体在带迎角出水时的数值模拟方法。Chu等人[10]模拟了圆柱航行体带空泡的出水过程,根据测量得出的压力数据能够得出空泡溃灭的大致时间和位置,发现空泡溃灭时会产生水拍打液面的现象,通过这种方式产生了高压,并且通过实验拍摄的图片与数值模拟产生的尾部空泡相对比,发现数值模拟产生的尾空泡更小一点;据此推测,这可能是由于在实际压缩气体喷射航行体的过程中,有部分气体会沿着航行体喷出,从而使实验中的尾空泡变得更大一些。张军等人[11]和赵蛟龙等人[12]分别对低速全沾湿物体和局部超空泡倾斜出水流场进行了PIV和高速摄影研究。
本课题组进行了细长钉体(横截面尺寸约为1mm)在小迎角α=5.5°~10°范围的出水实验和数值模拟[13]。本文将报道直径为6mm的钝体在较大迎角α=57°下的出水实验和数值模拟。这里使用“钝体”一词,是为了与以前的尖头钉体有所区别,但其长径比为8~12,实际也属于细长体。
1 实验装置与方法
图 1示出了实验装置。一台倾斜放置的一级轻气炮装置(见图 1中的部件4~12)向水箱自下而上地发射试验模型,模型冲击进入水箱1产生超空泡,然后穿过自由面2。圆柱体模型直径为6mm,长度分别为48、60和72mm(质量为3.69~5.50g),材料为铝镁合金。使用Photron Fastcam SA5高速摄影机记录超空泡穿过自由面的过程,然后将照片导入AutoCAD软件;测量位移等数据,每点测量3次,取平均值,最后计算出速度、阻力系数等物理量。根据冲击工程惯例,定义超空泡运动方向与水平液面的垂直方向之间的夹角为出水迎角α。有关实验技术的更多细节,可参见文献[14-15]。
2 实验结果与理论分析
图 2示出了长径比为10的平头圆柱体带迎角地穿过自由面的高速摄影照片。超空泡在图 2(1)中出现在视场中的左下角,其速度为V0=23.02 m/s;流场的标尺也在图 2(1)中给出。超空泡向右上方运动,逐渐靠近自由面(见图 2(2)~(7))。在图 2(8)所示的时刻,超空泡与自由面接触;然后在图 2(9)~(12)的时间段里,圆柱体伸出水面。从图 2(13)所示的时刻开始,圆柱体彻底脱离水面,进入空气中自由飞行。超空泡与自由面倾斜地相互作用后(见图 2(9)~(10)),出现了2个特殊的现象:一是在水面上沿着运动方向、即沿着超空泡的轴线,出现了不对称的水花。二是负压的超空泡从大气中倒吸空气,形成了充气空泡。如图 2(11)~(15)所示,水花先在轴线的上部形成并向右上方运动,水花头部的速度与圆柱体在空中的飞行速度接近(水花头部速度大于物体速度也是可能的)。从图 2(16)开始,轴线下部的水花才明显地形成,然后逐渐发展,与上部水花合并形成喇叭口的水花(见图 2(23)~(24))。水花以如此的方式形成,就使得空泡在水面上难以产生表面闭合(Surface Closure),空气被吸入负压的空泡中,因此空泡能在自由面附近长期存留而不会破灭。从图 2(11)~(24),空泡已存活了20ms的时间,而且空泡的最大直径从32mm增加到了46mm,尽管空泡后面的空化尾迹已经十分明显(见图 2(21)~(24))。如果对比本文第4节的数值计算结果,就可以理解上述的自然空化空泡向充气空泡转化的过程。
图 3(a)和(b)分别给出了测得的长径比为8、10和12的圆柱体的速度V随时间t的变化,以及阻力系数CD随空化数σ的变化。其中空化数σ是物体在每个时刻下计算出的,因此图 3(b)的横轴不需用时间t表示。从图 3(a)中可以看出,3个圆柱体在出水前后,即在t=10.00、12.86和14.27ms时刻,其速度都有一个跃增。这是因为物体从水中突然进入空气中甩负荷的缘故[5]。如果在穿过水面前后,在水中和在空气中的物体速度分别为Vw和Va,那么实验得到的Va/Vw=1.06~1.11。图 3(b)中给出的σ的最大值即是出水前的空化数,3个圆柱体的该空化数分别为σ=1.134、0.969和0.496。
圆柱体在水中和空气中的运动方程分别为:
(1) (2) 式中:m、aw、A、CD分别为物体的质量、加速度、截面积和阻力系数,ρ、g和α分别是流体密度、重力加速度和出水迎角,下标w和a分别表示在水中和空气中。将两式相除,整理后得到
(3) 根据图 3(a)可知,物体在出水前呈单调地减速,计算出的3个工况的加速度非常接近,其绝对值为1.04×103m/s2。运动物体在空气中的雷诺数Re=(5.93~7.01)×103,所以物体在空气中的阻力系数CDa约为1.00[16-17]。参考图 3(b),得知在出水前,物体的阻力系数CDw约为0.90~1.10,因此不妨认为CDw/CDa →1。然后从式(3)中算出Va/Vw=2.03。这个理论值大于实验值,这是因为物体刚出水时,虽然上部已露在空气中,但下部仍与水沾湿,此时阻力较大,会出现CDw/CDa < 1,这样计算出的Va/Vw的理论值会接近于实验值。
根据图 3(a)还可知,物体在出水后基本上保持恒速,因此式(2)中本应该出现的加速度项maa→0。上述有关出水前后物体速度变化的机理,如果从动量守恒定律来分析,就可以理解为:既然空泡引起的附加质量m*在出水后自动消失,而物体在穿过自由面时又没有受到其他外力作用,那么显然(m+m*)Vw=mVa。所以,
(4) 3 网格划分及计算方法设置
本文需要关注水中物体产生的空泡,所以在物体周围的网格划分较为密集,对物体周围加密时采用了尺寸函数,定义开始尺寸为1mm,增长率为1.2,最大尺寸为2mm。用Gambit软件对计算区域进行划分网格,且对物体周围进行局部加密,计算区域为1500mm×400mm,本文采用三角形非结构网格,网格总数为508 499。边界条件设置为:计算区域入口设置为wall边界,出口设置为压力出口,航行体边界设置为wall边界,如图 4所示。
将模型导入到Ansys14.0中,进行计算参数设置,具体步骤为:计算采用VOF多相流模型(水、空气和水蒸气),使用压力基求解器,计算考虑了重力的影响,重力加速度设置为-9.81m/s2。由于超空泡在运动过程中对时间的依赖性很强,所以选择非定常流动,空化模型采用的是Schnerr and Sauer模型,水的饱和蒸气压为3540Pa。求解方法采用SIMPLE算法。网格运动过程中的更新方法为弹簧近似光滑模型和网格重构模型实现,将编好的profile文件导入Ansys14.0中实现航行体的变速运动。
4 计算结果
图 5给出了水下超空泡物体倾斜出水过程相图的数值模拟结果,包括超空泡从生长到消失的全过程。已计算了各种不同的工况[15],这里只给出一组典型的数据。空泡开始出生时的位置主要发生在物体的尾部和肩部(见图 5(1)~(2)),随着时间的增大,空泡慢慢生长,直到将物体完全包裹、形成超空泡(见图 5(3)~(6))。然后物体带超空泡继续前进,直到穿出水面,并且在物体尾部带出水花(见图 5(15)~(16))。数值计算揭示了当水下超空泡物体接近自由面时(见图 5(5)~(8)),水面必然抬升隆起的、由水的不可压缩性所导致的后果。数值计算的另一个有意思的结果是,从图 5(14)~(16)的不到2ms的时间里,超空泡一接触水面就消失了,大大小于在图 2的实验中观察到的、大于20ms的超空泡的寿命。这是因为计算软件中的空化模型是完全自然空化模型,没有考虑到空泡从大气中吸取空气的机理。因此在今后建立充气空泡模型并改进软件的算法是十分必要的。另外,在第2节中讨论式(3)的计算结果时,认为物体刚出水时上半部分露在空气中,而下半部分与水沾湿,这个推断也被数值计算证实了。今后要注意的是:根据空化数的定义,当空泡浸没在水面下时,空泡壁面上的压力是水的饱和蒸气压;当空泡露出水面进行通气后,空泡壁面上的压力是大气压,有限体积算法不一定能处理好这种问题。Young和Savander[18]用边界元方法计算了部分浸没的空泡以及表面划水螺旋桨叶片的流场,值得参考。
图 6给出了3个不同时刻水下物体周围的压力图(压力的单位是Pa),其相图分别是图 5(3)、(4)和(8)。从图 6(a)中可以看出,在t=2.57ms时,超空泡离自由面还较远,此时从相对应的压力图中可以看出,物体的周围压力对称分布。随着超空泡靠近自由面(见图 6(b)和(c)),物体头部下方的压力比上方大。在研究浅水区的超空泡运动特性时,也发现了空泡下部的压力大于上部的压力[19]。这说明当超空泡靠近自由面时,浮力发挥作用,这对空泡后期的发展,包括形状和溃灭特性,产生一定的影响。
5 结论与讨论
(1) 超空泡倾斜地与自由面相互作用时,会形成非轴对称的、开口喇叭形的水花,因此原来处于负压的空泡更容易地吸取大气中的空气,空泡从自然空化向充气空泡转变,空泡寿命增加(大于20ms)。
(2) 物体出水前后,速度有一跃增,速度增加率达到6%~11%。通过分析运动方程,从理论上证明了当物体从水中进入空气中时,其速度的确是增加的。运用动量守恒定律,也得出相同结论。对于长径比为10的物体,如图 3(a)中所示的、在t=21.43ms时刻出现的二次增速,是由测量误差造成的;更多的测量结果,可见文献[15]。
(3) 当水中超空泡物体接近自由面时,自由面将隆起抬高,这是后期表面波浪形成的初始原因;自由面对超空泡的影响之一,是空泡下部的压力大于上部的压力。严格地说,用本文的方法产生的超空泡里,会夹带空气的进入。但是,实验中发射管出口离水箱侧面隔膜较近,这样对于空气进入有抑制作用。而数值计算的空泡形状与实验结果比较接近,说明即使有空气进入,其量也是有限的。
-
表 1 六分量应变天平和单分量应变天平载荷
Table 1 Measurement range of six-component and single-component strain-gage balance
Fx/N Fy/N Fz/N Mx/(N·m) My/(N·m) Mz/(N·m) Mk/(N·m) 1950 9500 850 400 490 620 1800 表 2 风洞试验状态
Table 2 Matrix of test in wind tunnel
前进比μ θFT /(°) αs/(°) Cw 0.092 6 8.1 0.0044 0.138 6 5.4 0.0044 0.150 6 4.9 0.0044 0.150 3 1.9 0.0044 0.150 0 −1.1 0.0044 0.150 −3 −4.1 0.0044 0.150 −6 −7.1 0.0044 -
[1] YU Y H. Rotor blade–vortex interaction noise[J]. Progress in Aerospace Sciences, 2000, 36(2): 97–115. doi: 10.1016/S0376-0421(99)00012-3
[2] Du Vall T, Sim B, Schmitz F. Cabin versus Far-Field Blade- Vortex Interaction Noise Level Trends[C]// Proc of American Helicopter Society Aerodynamics, Acoustics and Test Evaluation Technical Specialists Meeting. 2002.
[3] GENNARETTI M, BERNARDINI G. Novel boundary integral formulation for blade-vortex interaction aerodynamics of helicopter rotors[J]. AIAA Journal, 2007, 45(6): 1169–1176. doi: 10.2514/1.18383
[4] FOGARTY D E, WILBUR W L, SEKULA M K. Prediction of BVI noise for an active twist rotor using a loosely coupled CFD/CSD method and comparison to experimental data[R]// NF1676L-14452, 2012.
[5] 史勇杰, 苏大成, 徐国华. 桨叶气动外形对直升机桨–涡干扰噪声影响研究[J]. 南京航空航天大学学报, 2015, 47(2): 235–242. DOI: 10.16356/j.1005-2615.2015.02.009 SHI Y J, SU D C, XU G H. Research on influence of shape parameters on blade-vortex interaction noise of helicopter rotor[J]. Journal of Nanjing University of Aeronautics & Astronautics, 2015, 47(2): 235–242. doi: 10.16356/j.1005-2615.2015.02.009
[6] 史勇杰, 徐国华, 王菲. 直升机旋翼桨–涡干扰脉冲噪声传播特性研究[J]. 南京航空航天大学学报, 2014, 46(2): 212–217. SHI Y J, XU G H, WANG F. Propagation characteristics of helicopter rotor blade-vortex interaction noise[J]. Journal of Nanjing University of Aeronautics & Astronautics, 2014, 46(2): 212–217.
[7] 王菲, 徐国华, 胡志远. 大气环境对直升机旋翼桨–涡干扰噪声辐射特性的影响[J]. 南京航空航天大学学报, 2020, 52(2): 304–310. WANG F, XU G H, HU Z Y. Effects of atmospheric environment on helicopter blade-vortex interaction noise radiation characteristics[J]. Journal of Nanjing University of Aeronautics & Astronautics, 2020, 52(2): 304–310.
[8] BROOKS T F, JOLLY R J, MARCOLINI M A. Determination of noise source contributions using scaled model rotor acoustic data[R]. NASATP-2825, 1998.
[9] HELLER H, SPLETTSTOESSER W, KLOEPPEL V, et al. HELINOISE—The European Community rotor acoustics research program[C]//Proc of the 15th Aeroacoustics Conference. 1993. doi: 10.2514/6.1993-4358
[10] SPLETTSTOESSER W R, NIESL G, CENEDESE F, et al. Experimental results of the European HELINOISE aeroacoustic rotor test[J]. Journal of the American Helicopter Society, 1995, 40(2): 3–14. doi: 10.4050/jahs.40.2.3
[11] HELLER H, BUCHHOLZ H, SCHULTZ K, et al. Helicopter rotor blade aeroacoustics: a comparison of model-scale wind tunnel and full-scale flight test results [C]// Proc of 20th ICAS Congress. 1996.
[12] YU Y H, GMELIN B, HELLER H, et al. HHC aeroacoustics rotor test at the DNW—the joint German/French/US HART project[C]. Proceedings of the 20th European Rotorcraft Forum. 1994.
[13] SPLETTSTOESSER W R, KUBE R, WAGNER W, et al. Key results from a higher harmonic control aeroacoustic rotor test (HART)[J]. Journal of the American Helicopter Society, 1997, 42(1): 58–78. doi: 10.4050/jahs.42.58
[14] YU Y H, TUNG C, VAN DER WALL B G. The HART-II Test: Rotor Wakes and Aeroacoustics with Higher-Harmonic Pitch Control (HHC) Inputs —The Joint German/French/Dutch/US Project[C]//Proc of 58th Annual Forum of the American Helicopter Society. 2002.
[15] VAN DER WALL B G, BURLEY C L, YU Y, et al. The HART II test - measurement of helicopter rotor wakes[J]. Aerospace Science and Technology, 2004, 8(4): 273–284. doi: 10.1016/j.ast.2004.01.001
[16] JAYARAMAN B, WISSINK A, LIM J, et al. Helios prediction of blade-vortex interaction and wake of the HART II rotor[C]//Proc of the 50th AIAA Aerospace Sciences Meeting including the New Horizons Forum and Aerospace Exposition, Nashville, Tennessee. 2012: 714. doi: 10.2514/6.2012-714
[17] ARUN KUMAR A, VISWAMURTHY S R, GANGULI R. Correlation of helicopter rotor aeroelastic response with HART-II wind tunnel test data[J]. Aircraft Engineering and Aerospace Technology, 2010, 82(4): 237–248. doi: 10.1108/00022661011082713
[18] YIN J P, WALL B V D, OERLEMANS S. Representative test results from HeliNOVI aeroacoustic main rotor/tail rotor/fuselage test in DNW[C]// Proc of 31th European Rotercraft Forum. 2005.
[19] YIN J P, A DUMMEL, D FALCHERO. Analysis of Tail Rotor noise reduction benefits using HELINOVI aeroacoustic main/tail rotor test and posttest prediction results[C]// Proc of 32th European Rotercraft Forum. 2006.
[20] 徐国华, 高正. 悬停状态下模型旋翼噪声试验的初步研究[J]. 空气动力学学报, 1996(1): 68–72. XU G H, GAO Z. A preliminary investigation of noise experiment for helicopter model rotor in hover[J]. Acta Aerodynamica Sinica, 1996(1): 68–72.
[21] 曹亚雄, 樊枫, 林永峰, 等. 带先进桨尖的模型旋翼悬停噪声计算与试验[J]. 南京航空航天大学学报, 2016, 48(2): 180–185. DOI: 10.16356/j.1005-2615.2016.02.005 CAO Y X, FAN F, LIN Y F, et al. Numerical calculations and test research on aeroacoustics characteristics of model rotors with advanced blade tip in hover[J]. Journal of Nanjing University of Aeronautics & Astronautics, 2016, 48(2): 180–185. doi: 10.16356/j.1005-2615.2016.02.005
[22] 刘正江, 黄建萍, 陈焕, 等. 旋翼桨涡干扰噪声特性试验技术研究[J]. 直升机技术, 2019(1): 43–47,42. DOI: 10.3969/j.issn.1673-1220.2019.01.010 LIU Z J, HUANG J P, CHEN H, et al. Study on characteristic of rotor-blade vortex interaction noise[J]. Helicopter Technique, 2019(1): 43–47,42. doi: 10.3969/j.issn.1673-1220.2019.01.010
[23] 唐朝, 招启军, 王博, 等. 用于BVI噪声试验的新型涡发生器设计与分析[J]. 南京航空航天大学学报, 2018, 50(2): 157–166. DOI: 10.16356/j.1005-2615.2018.02.002 TANG C, ZHAO Q J, WANG B, et al. Design and analysis of new type vortex generator for BVI noise experiment[J]. Journal of Nanjing University of Aeronautics & Astronautics, 2018, 50(2): 157–166. doi: 10.16356/j.1005-2615.2018.02.002
[24] HEYSON H H. Linearized theory of wind-tunnel jet-boundary corrections and ground effect for VTOL-STOL aircraft [R]. NASA TR R-124, 1962.
[25] HEYSON H H. FORTRAN programs for calculating wind-tunnel boundary interference [R]. NASA TM X-1740, 1969.
[26] HEYSON H H. Use of superposition in digital computers to obtain wind tunnel interference factors for arbitrary configurations, with particular reference to V/STOL models[R]. NASA TR R-302, 1969.
[27] LANGER H, PETERSON R L, MAIER T H. An experimental evaluation of wind tunnel wall correction methods for helicopter performance [C]// Proceedings of the AHS 52nd Annual Forum. 1996.
[28] 李元首, 陈宝, 张雪, 等. 传声器阵列校准技术研究[J]. 现代电子技术, 2014, 37(24): 94–97. DOI: 10.3969/j.issn.1004-373X.2014.24.026 LI Y S, CHEN B, ZHANG X, et al. Calibration technology of microphone array[J]. Modern Electronics Technique, 2014, 37(24): 94–97. doi: 10.3969/j.issn.1004-373X.2014.24.026
-
期刊类型引用(3)
1. 田建辉,胡晨明. 高速弹体出水过程数值模拟. 兵器装备工程学报. 2023(12): 73-78+176 . 百度学术
2. 李易松,袁航,张建宇,冯会方. 航行体跌落试验研究. 舰船科学技术. 2023(24): 218-221 . 百度学术
3. 侯健,施红辉,孙亚亚,高聚瑞. 高速射弹出水过程中水弹道问题研究. 弹道学报. 2017(01): 51-56 . 百度学术
其他类型引用(2)