+高级检索
网刊加载中。。。

使用Chrome浏览器效果最佳,继续浏览,你可能不会看到最佳的展示效果,

确定继续浏览么?

复制成功,请在其他浏览器进行阅读

小尺寸切开式尾缘锯齿对翼型气动噪声的影响  PDF

  • 张喆 1
  • 陈涛 1
  • 张英朝 1
  • 王中检 1
  • 张成春 2
  • 沈淳 1
  • 任露泉 2
1. 吉林大学 汽车底盘集成与仿生全国重点实验室,吉林 长春 130022; 2. 吉林大学 工程仿生教育部重点实验室,吉林 长春 130022

中图分类号: U461.1

最近更新:2024-10-28

DOI: 10.16339/j.cnki.hdxbzkb.2024203

  • 全文
  • 图表
  • 参考文献
  • 作者
  • 出版信息
EN
目录contents

摘要

研究小尺寸切开式尾缘锯齿对翼型气动噪声的影响,并分析齿高与尾缘边界层厚度之间的关系对降噪效果的影响. 以NACA0012翼型为研究对象,采用混合计算气动声学方法进行仿真.在20 m/s风速下对NACA0012翼型进行噪声实验,验证仿真方法的准确性. 通过控制齿高和齿宽研究不同形状尾缘锯齿分别在0°和10°攻角下对翼型气动噪声的影响. 研究结果表明,尾缘锯齿的降噪效果与齿高和齿宽均成正比. 0°攻角下最大降噪量为3.3 dB,10°攻角下最大降噪量为2.6 dB. 切开式尾缘锯齿结构会增加流向噪声的高频分量,尤其在0°攻角下显著. 相比嵌入式尾缘锯齿,切开式尾缘锯齿具有更小的极限降噪齿高.

翼型广泛运用于航空、风电等领域,对翼型进行降噪研究具有十分重要的现实意义. 现有关于翼型的降噪研究主要以被动控制技术为主,其中翼型前缘处

1-2、尾缘处3、表面微结4以及翼型穿5等研究占大多数.

自Howe

6提出尾缘锯齿降噪理论后,对翼型尾缘锯齿的降噪研究众多.大体分为嵌入式尾缘锯齿研究和切开式尾缘锯齿研7. 嵌入式尾缘锯齿的研究以Gruber8的实验研究最为著名. 以上研究提出了一些重要结论,其中包括当h/δ<0.25时(h为锯齿长度的一半,δ为尾缘边界层厚度,在原文中由XFoil估计的边界层位移厚度乘以固定比值获得),尾缘锯齿基本没有降噪效果,这一结论对后续研究产生了很大的影响,以至于在之后研究中,尾缘锯齿长度多使用较大尺寸. 在当前的研究7,尾缘锯齿长度均要高于0.05CC为翼型弦长. Gruber8研究中,尾缘锯齿最小相对齿高为0.006 7C. Gelot9、汪瑞欣10分别对切开式尾缘锯齿茹科夫斯基翼型以及风力机专用翼型DU91-W2-250气动噪声的影响进行研究,相对锯齿高度均为0.025C. Gelot9的研究表明,在0°攻角下,0.025C的尾缘锯齿可以有效降低翼型气动噪声,运行工况雷诺数为2.5×105,马赫数为0.4;汪瑞欣10的研究表明,在5°攻角下,0.025C的尾缘锯齿会增大翼型周向噪声,运行工况雷诺数为2×106,马赫数约为0.14. Gelot9、汪瑞欣10的研究均未对翼型尾缘边界层厚度进行说明,研究结果有较大的差异.

本文将锯齿高度小于0.05C的锯齿定义为小尺寸锯齿,利用声类比方法探究小尺寸切开式尾缘锯齿对翼型气动噪声的影响,通过实验数据验证仿真方法的准确性,并研究尾缘锯齿参数对翼型降噪效果的影响,为确定切开式尾缘锯齿降噪极限齿长给予一定参考.

1 噪声实验

实验翼型为NACA0012,弦长C=0.2 m,展长为0.4 m. 如图1(a)所示,噪声实验在中国空气动力研究与发展中心低速空气动力研究所进行,喷口面积为0.55 m×0.4 m,实验风速为20 m/s,翼型攻角为0°. 图1(b)为麦克风监测点布置,监测点距离翼型中心点为R=1 m. 将监测到的声压信息进行傅里叶变换,获得麦克风监测点频谱信息.

fig

(a)  翼型噪声实验布置

fig

(b)  麦克风监测点布置

图1  实验示意图

Fig.1  Experimental diagram

2 数值仿真

2.1 几何模型

锯齿参数结构示意图如图2所示.本文对NACA0012翼型的尾缘进行了改进,通过调整齿高2h和齿宽λ来控制尾缘形状.设置4组算例,如表1所示,其中包含3种不同齿高(2h=0.01C、0.02C、0.03C)、2种不同齿宽(λ=0.04C、0.08C)的锯齿.

fig

图2  锯齿参数结构示意图

Fig.2  Serration parameter structure schematic

表1  各算例参数设置表
Tab.1  Parameter settings list of each case
算例2h/Cλ/h
基础翼型 0
2h=2 mm,λ=4 mm 0.01 4
2h=4 mm,λ=4 mm 0.02 2
2h=6 mm,λ=4 mm 0.03 1.33
2h=6 mm,λ=8 mm 0.03 2.67

2.2 计算域与网格

计算域示意图如图3所示.计算域由前半圆柱和方形柱组成.圆柱半径0.9 m,方形柱长度为 1.8 m,翼型中点处于半圆柱圆心处.计算域左右边界条件设置为对称面,翼型表面设置为无滑移壁面,来流进出口设置为自由流.采用非结构化网格对计算域进行离散,将壁面第一层网格厚度设置为3.3×10-6 m,壁面延伸率为1.2,总边界层数为35层,壁面xz方向网格尺度为5×10-4 m. 此外,对网格进行逐层加密,以过渡到内部网格基本与边界层网格尺寸相当,网格总数为1 468万.最终生成的网格示意图如图4所示.

fig

图3  计算域示意图

Fig.3  Computing domain diagram

fig

图4  网格示意图

Fig.4  Grid diagram

2.3 物理模型

本文利用大涡模拟(Large Eddy Simulation, LES)对流场时域信息进行捕捉,通过壁面自适应局部涡黏(Wall-Adapting Local Eddy-Viscosity,WALE)模型计算亚格子黏度.仿真在STARCCM+中进行. 式(1)式(2)式(3)分别为LES连续性方程、动量方程和能量方程的张量形式,是由N-S方程滤波后的结果,需要注意的是,其中压力以及速度分量是滤波后的量. 黏性应力可以通过涡黏假设用速度表示,若要封闭方程需要对亚格子黏度进行求解.从式(3)式(4)中可以看出,亚格子应力的求解可以转化为对亚格子黏度的求解.

ρt+xiρU˜i=0 (1)
ρU˜it+xjρU˜iU˜j=-P˜xi+xjτ˜ij+τsgs (2)
ρE˜t+UjρE˜xj=-U˜iP˜xi+U˜iτ˜ij+τsgsxj-q˜jxj (3)
τsgs=2μsgsSij-23μsgsU˜kkδij (4)
Sij=12U˜ixj+U˜jxi (5)

式中:ρ为密度;U˜i为滤波速度分量;P˜为流体微元所受滤波表面力;τ˜ij为滤波黏性应力张量;τsgs为亚格子应力张量;E˜为单位质量的滤波总能量;q˜j为滤波热通量;μsgs为亚格子黏度;δij为克罗内克符号;Sij为应变率张量.

在WALE亚格子模型中,对亚格子黏度进行如下建模:

μsgs=ρΔ2Sw (6)
Δ=CsV1/3通常处理minkd,CsV1/3近壁处理 (7)
Sw=sd:sd3/2sd:sd5/4+sij:sij3/2 (8)
Sd=12UkxiUjxk+UkxjUixk-13δijUkxiUjxk (9)

式中:Δ表示网格尺度,对于近壁网格通常会特殊处理;Sw为变形参数;Sd为与速度梯度的二阶导数相关的张量;Cs为模型系数,在STARCCM+中取0.544;k为冯卡门常数,取0.41.

为加快瞬态收敛,在此之前对算例进行稳态仿真,选用K-Omega SST湍流模

11.将稳态收敛结果用于LES初始化,物理总时长为0.2 s,求解器时间步长为1×10-5 s,时间离散采用二阶格式,在0.1 s后开始对FW-H (Fowcs Williams-Hawkings)方程进行求解,内部最大迭代步数为5,声学采样步长为2×10-5,对应频率分辨率为10 Hz. 式(10)式(11)为对应FW-H方程的解,分别表示单极子声源和偶极子声源. 积分下限符f=0,表明单极子声源和偶极子声源仅出现在固体表面. LES计算完成后,可以获得固壁表面流动时域信息,即噪声源. FW-H监测点布置如图5所示,由此可获得监测点时域声压信息.

fig

图5  FW-H监测点布置

Fig.5  FW-H monitoring site layout

4πpT'x,t=f=0ρ0U˙n+Un˙r(1-Mr)2retdS+          f=0ρ0UnrM˙r+c(Mr-M2)r3(1-Mr)3retdS (10)
4πpL'x,t=1cf=0L˙rr1-Mr2retdS+        f=0Lr-LMcr(1-Mr)2retdS+        1cLrrM˙r+cMr-M2r2(1-Mr)3retdS (11)
Un=Uini,Lr=Liri,LM=LiMi (12)
Ui=vi+ρρ0ui-vi (13)
Li=Pijn^j+ρuiun-vn (14)
Pij=pδij-μuixj+ujxi-23ukxkδij (15)

式中:下标ret表示迟至时间;ρ0为无穷远处空气密度;r为监测点距离声源位置;c为声速;M为声源位置马赫数,M=u/cMr为声源位置在辐射方向上的马赫数;下标rn分别表示该向量与辐射方向r上的单位向量或表面法线方向上的单位向量的点积.

2.4 数值验证

对如图5所示的监测点进行FW-H监测,监测点实验与仿真噪声频谱对比结果如图6所示. 由图6可知,在低频阶段(低于200 Hz)内,实验与仿真结果吻合较差,这是由于风洞中背景噪声所致;在高频阶段(高于6 000 Hz)难以捕捉到信息,这是由于受到数值方法限制;在中间频段吻合较好,在一定程度上说明仿真的准确性.

fig

图6  实验与仿真噪声频谱对比结果

Fig.6  Noise spectrum comparison results between experimental and simulated noise

3 结果分析

3.1 0°攻角工况

0°攻角下小尺寸尾缘锯齿参数对监测点总声压级的影响如图7所示.从图7中可以看出,在50 m/s工况下,翼型产生的噪声以二极子声源为主.仿真结果显示,在0°攻角下,单极子噪声占比几乎为0.图7(a)对比了锯齿相对长度对各个监测点总声压级(Overall Sound Pressure Level, OASPL)的影响. 小尺寸锯齿可以明显降低垂直于流动方向的声压级,且在0.01C~0.03C内,锯齿长度越大,降噪效果越明显,在135°监测点位置相对尺高为0.03的锯齿降噪3.1 dB. 在Gruber的理论中,δ值由XFoil估计的边界层位移厚度δ*≈1获取,δ/δ*≈8. 以当地速度47.3 m/s为参考速度,在仿真计算结果中测得边界层厚度约为8.6 mm,与估计值相差7%. 这里以估计值为δ大小,各算例对应的h/δ值分别为0.12、0.23、0.35. 一方面,切开式尾缘锯齿可以在h/δ<0.25时仍然降低翼型垂向噪声;另一方面,流向声压级在小尺寸锯齿的影响下开始增加,且锯齿长度越大,增量越大,在0°监测点位置,相对尺高为0.03的锯齿使得总声压级增加3.4 dB. 由图7(b)可知,齿宽的增加使得流向声压级增幅降低,相比基础翼型,0°监测点位置仅增加0.8 dB. 在此基础上,仍然能够降低垂向噪声,并在135°监测点达到3.3 dB的降噪量.

fig

(a)  锯齿相对长度

fig

(b)  锯齿宽度

图7  0°攻角下小尺寸尾缘锯齿参数对监测点

Fig.7  Effect of small size trailing edge serration parameters on overall sound pressure level at 0° angle of attack

总声压级的影响

图8为0°攻角下锯齿参数对135°监测点频谱的影响(1/3倍频程).由图8可知,小尺寸锯齿在该监测点并没有影响频谱图的趋势,添加小尺寸锯齿的翼型各频段声压级均有所降低,集中体现在中频段区域.

fig

(a)  锯齿长度频谱对比图

fig

(b)  锯齿宽度频谱对比图

图8  0°攻角下锯齿参数对135°监测点频谱的影响(1/3倍频程)

Fig.8  Effect of serration parameters on the frequency spectrum of 135° monitoring point at 0° angle of attack (1/3 octave)

图9为0°攻角下锯齿参数对0°监测点频谱的影响(1/3倍频程).由图9可知,添加锯齿结构后,监测点高频成分增加,高频增量随锯齿长度增加而增加,随锯齿宽度增加而减小.

fig

(a)  锯齿长度频谱对比图

fig

(b)  锯齿宽度频谱对比图

图9  0°攻角下锯齿参数对0°监测点频谱的影响(1/3倍频程)

Fig.9  Effect of serration parameters on the frequency spectrum of 0° monitoring point at 0° angle of attack (1/3 octave)

二极子噪声源是由于翼型受力的波动产生,主要以翼型升力波动为主.而翼型升力的波动与其尾部涡旋的周期性脱落相关.0°攻角下,锯齿对翼型尾缘涡量的影响如图10所示,在小尺度锯齿的作用下,翼型尾部涡量强度明显减小. 锯齿加快了翼型尾缘的三维流动,破坏了翼型横向流动的一致性,降低了流体至翼型脱落对其产生的作用力,使得声压级降低. 噪声流向分量是翼型阻力波动的体现. 由于该工况下翼型阻力波动极小,声压级量级也很小,锯齿形状导致流动在翼型尾缘复杂度增加,高频成分相应增加.

fig

图10  0°攻角下锯齿对翼型尾缘涡量的影响

Fig.10  Effect of serration on trailing edge vorticity of airfoil at 0° angle of attack

(a)基础翼型 (b)2h=2 mm (c)2h=4 mm (d)2h=6 mm (e)2h=6 mm

λ=4 mm λ=4 mm λ=4 mm λ=8 mm

3.2 10°攻角工况

图11为10°攻角下尾缘锯齿参数对监测点总声压级的影响.相比0°攻角,噪声流向成分明显增加,这是由于翼型阻力波动幅值增大的结果.此外,锯齿垂向降噪作用降低.可以看到2h=2 mm,λ=4 mm锯齿已经没有降噪效果,随相对齿高增加,锯齿降噪效果逐渐显著,流向噪声也有所降低.在2h=6 mm条件下,增加齿宽齿高比可以进一步降低周向声压级.10°攻角下,随着相对齿高的增加,流向声压级逐渐降低,并开始显现出降低流向噪声的效果,降噪效果随着齿宽齿高比的增大而增大.在135°监测点位置,最大降噪量达2.6 dB.

fig

(a)  锯齿长度

fig

(b)  锯齿宽度

图11  10°攻角下尾缘锯齿参数对监测点总声压级的影响

Fig.11  Effect of small size trailing edge serration parameters on overall sound pressure level at 10° angle of attack

图12为尾缘锯齿参数对翼型中心截面尾缘处的速度分布.yu分别表示纵向坐标位置和流向速度,利用翼型特征长度与流速U=50 m/s对纵向位置与流向速度进行归一化.10°攻角下翼型吸力面与压力面边界层位移厚度分别为3.41 mm和2.83 mm.翼型在攻角较大的情况下,尾缘边界层位移厚度更大,这一特性使得小锯齿(2h=2 mm,λ=4 mm)在10°攻角下降噪效果降低.

fig

(a)  0°攻角翼型中心截面尾缘处的速度分布

fig

(b)  10°攻角翼型中心截面尾缘处的速度分布

图12  尾缘锯齿参数对翼型中心截面尾缘处的速度分布

Fig.12  Trailing edge serration parameters on velocity distribution at trailing edge of center section of airfoil

图13为10°攻角下锯齿参数对135°监测点频谱的影响(1/3倍频程),由图13可知,尾缘锯齿对整个频域范围内的降噪作用没有0°工况下显著.

fig

(a)  锯齿长度频谱对比图

fig

(b)  锯齿宽度频谱对比图

图13  10°攻角下锯齿参数对135°监测点频谱的影响(1/3倍频程)

Fig.13  Effect of serration parameters on the frequency spectrum of 135° monitoring point at 10°

angle of attack(1/3 octave)

4 结 论

本文基于声类比方法研究了小尺寸切开式尾缘锯齿对翼型气动噪声的影响. 结论如下:

1) 在本文参数范围内,锯齿对垂直于来流方向的噪声的降噪作用与锯齿长度成正比,与锯齿宽度成正比. 0°攻角最大降噪量为3.3 dB,10°攻角下最大降噪量为2.6 dB.

2) 锯齿尾缘可以降低翼型尾缘垂向噪声的产生,当锯齿尺寸较小时,尾部的复杂流动产生高频的声学分量会使得流向声压级增大.

3) 相比嵌入式尾缘锯齿,切开式尾缘锯齿具有更小的降噪极限齿高,当h/δ<0.12时,尾缘锯齿仍然能较好地降低翼型垂向噪声.

参考文献

1

任露泉孙少明徐成宇. 鸮翼前缘非光滑形态消声降噪机 理[J]. 吉林大学学报(工学版)200838增刊1): 126-131 [百度学术] 

REN L QSUN S MXU C YNoise reduction mechanism of non-smooth leading edge of owl wing[J].Journal of Jilin University (Engineering and Technology Edition)200838Sup.1):126-131(in Chinese) [百度学术] 

2

杨文琪赛庆毅郭丹阁轴流通风机仿生降噪方法[J].动力工程学报2022426):530-536 [百度学术] 

YANG W QSAI Q YGUO D Get alA bionic approach for axial flow fans to reduce the noise[J].Journal of Chinese Society of Power Engineering2022426):530-536(in Chinese) [百度学术] 

3

刘宴利仝帆王勇带仿鳍式尾缘结构有限长翼型的噪声特性及降噪实验研究[J].西安交通大学学报2022569):20-27 [百度学术] 

LIU Y LTONG FWANG Yet alExperimental study on noise characteristic and noise reduction of finite airfoil with fin-like trailing edge structure[J].Journal of Xi’an Jiaotong University2022569):20-27(in Chinese) [百度学术] 

4

党志高毛昭勇田文龙基于微沟槽表面海流发电叶片翼型降噪特性研究[J].工程热物理学报2020416):1367-1374 [百度学术] 

DANG Z GMAO Z YTIAN W LNoise reduction characteristics of hydrofoils of marine current turbines with microgrooves on the surface[J].Journal of Engineering Thermophysics2020416):1367-1374(in Chinese) [百度学术] 

5

马扬杨爱玲陈二云尾缘双穿孔对翼型性能影响研究[J].热能动力工程20193412):122-130 [百度学术] 

MA YYANG A LCHEN E YResearch of effect on the performance of airfoil with double-punched trailing edge[J].Journal of Engineering for Thermal Energy and Power20193412):122-130(in Chinese) [百度学术] 

6

HOWE M SAerodynamic noise of a serrated trailing edge[J].Journal of Fluids and Structures199151):33-45 [百度学术] 

7

乔渭阳仝帆陈伟杰. 仿生学气动噪声控制研究的历史、现状和进展[J]. 空气动力学学报2018361): 98-121 [百度学术] 

QIAO W YTONG FCHEN W Jet alReview on aerodynamic noise reduction with bionic configuration[J]. Acta Aerodynamica Sinica2018361): 98-121(in Chinese) [百度学术] 

8

GRUBER M. Airfoil noise reduction by edge treatments[D].Southampton,South East England,UKUniversity of Southampton2012. [百度学术] 

9

GELOT M B RKIM J WEffect of serrated trailing edges on aerofoil tonal noise[J].Journal of Fluid Mechanics2020904A30 [百度学术] 

10

汪瑞欣赵振宙王同光锯齿尾缘DU91-W2-250风力机翼型气动及噪声特性分析[J].太阳能学报20204112):221-228. [百度学术] 

WANG R XZHAO Z ZWANG T Get alResearch on aerodynamics and noise characteristics of DU91-W2-250 airfoil of wind blade with serrated trailing edge[J].Acta Energiae Solaris Sinica20204112):221-228(in Chinese) [百度学术] 

11

WANG H PJIANG XCHAO Yet alEffects of leading edge slat on flow separation and aerodynamic performance of wind turbine[J].Energy2019182988-998 [百度学术] 

作者稿件一经被我刊录用,如无特别声明,即视作同意授予我刊论文整体的全部复制传播的权利,包括但不限于复制权、发行权、信息网络传播权、广播权、表演权、翻译权、汇编权、改编权等著作使用权转让给我刊,我刊有权根据工作需要,允许合作的数据库、新媒体平台及其他数字平台进行数字传播和国际传播等。特此声明。
关闭