Investigation of acoustic liner vibroacoustic response and its influence on impedance
-
摘要: 声衬在高声压级声波激发下产生声振响应,刚性假设不再成立,其结构振动会对吸声产生一定影响。本文针对振动对声衬吸声的影响和声振响应开展实验研究,通过参数化研究,获得了不同工况和不同穿孔板几何参数下声振响应对声阻抗的影响规律。实验结果表明:穿孔板振动会导致声阻在结构共振频率处出现波峰或波谷,吸声系数出现“额外”的吸声波峰或吸声波谷;穿孔率和声压级的增大会减弱振动的影响,且存在一个临界穿孔率;穿孔板参数会影响高声压级下结构振动导致声阻变化的特征;在结构共振频率附近,小孔–面板速度相位差会发生突变,导致相对速度增大,吸声效果改变。Abstract: The acoustic liners produce vibroacoustic response under the excitation of sound waves at high sound pressure level, and the rigid structure assumption is no longer applied. Their structural vibrations have a certain impact on sound absorption performance. The work presented here is an experimental study on the influence of panel vibrations on sound absorption and vibroacoustic response, and the influence law of vibroacoustic response on acoustic impedance under different perforated plate geometric parameters is obtained through parametric research. The experimental results show that the vibration of the perforated plate causes resistance to the generation of peaks or dips at the structural resonance frequency, and the sound absorption coefficient generates extra absorption peaks or dips that cannot be understood assuming rigid acoustic liners. The increase of the perforation rate and sound pressure level suppresses the influence of vibration, and there is a critical perforation rate. Perforated plate parameters affect the characteristics of resistance changes caused by structural vibration at high sound pressure levels. The phase difference between the small holes and the panel near the structure resonance frequency changes abruptly, resulting in an increase in the relative velocity and a change in the sound absorption performance.
-
Keywords:
- acoustic liner /
- acoustic impedance /
- vibroacoustic
-
0 引 言
高速飞行器设计日趋精细化,对风洞试验的精准度提出了更高的要求,为此,需研制高性能连续式跨声速风洞,以提高飞行器跨声速试验的流场品质[1-3]。其中,湍流度(速度脉动)是流场先进性的重要指标之一,较高的气流脉动会影响模型的绕流,对非定常试验、边界层转捩试验和洞体结构疲劳都有着非常不利的影响[3-4]。为此,针对压缩机、阀门、换热器、扩散段和试验段产生的气流脉动,需采取相应的技术措施以降低风洞气流脉动[5-6]。研究风洞降湍特性对流场湍流度进行准确评估是当前一项重要工作[7],然而,如何测量低湍流度尤其是跨声速可压流场低湍流度仍有待研究[8-9]。
目前,风洞流场湍流度主要使用热线风速仪(简称热线)测量。其中,恒温式热线风速仪(Constant Temperature Anemometer,CTA)商业化程度高、可及性较好,是国内外广泛使用的主流热线风速仪。CTA基于King公式的低速不可压流的湍流度测量方法已广泛用于评估低速风洞流场湍流度水平及其对模型试验的影响[10-11]。在常规工况下,连续式风洞试验段的上游回路风速较低,是不可压流动,也可采用King方法进行热线测量。
连续式跨声速风洞湍流度测量的主要难点在于试验段高速可压流场的测量。首先,需要解决热丝易断的问题,连续式风洞缺乏理想的洞体清洁和吹扫手段,因此,其气体洁净度远低于使用洁净气罐的暂冲式风洞,风洞回路中的微尘、金属屑等颗粒物随着气流循环经过试验段,在高速状态下极易打断热丝。其次,需要解决试验段可压流场数据分析的难点,相关测试方法一直是国内外学者不断探讨的研究内容。Kovasznay最早建立了基于连续变热线过热比的超声速流场测量方法[12],这一方法在超声速流场脉动量测量领域得到了不断优化和发展 [13-14]。在跨声速流场测量中的研究也在不断完善[9]。近年,King等[15]开展了热线低温校准,完成最低温度116 K的跨声速流场测试,根据湍流度测量结果评估了美国NTF风洞的层流试验可行性。Bauinger等[16]在涡轮机试验中采用三维热线进行了跨声速流场测量,根据傅里叶滤波方法获得湍流度分布。Andrews[17]开展了热线测量高超声速流场附面层转捩研究。国内早期关于跨声速流场湍流度测量方面的研究较少,近年来,流场测量精细化要求促进了相关技术的发展。杜钰锋等[9]优化了热线变过热比测量方法,完成了马赫数小于0.7的风洞可压流场湍流度测量。马护生等[18]建立了可压流热线校准数学模型,并利用自建的热线校准风洞开展了热线校准试验和湍流度测量试验。袁湘江等[19]分析了超声速流场中脉动压力与脉动速度的本质联系,从理论上得到了声压级与声致湍流度的关系。
以此可见,热线技术是测量风洞流场湍流度的主要方法,国外较早将其用于超声速流场湍流度的测量,近年关注于跨声速和高超声速流场精细化测量。在国内,随着高性能风洞研制需求越来越迫切,近年对可压流场热线测试的理论研究较多。但是,对马赫数大于0.7的流场湍流度测量理论和试验研究较少,尤其对连续式跨声速风洞的湍流度。
为了解连续式跨声速风洞流场湍流特征,评估和优化关键部段的降湍效果,本文在0.6 m连续式跨声速风洞换热器入口至试验段开展了湍流度测试试验研究,采用二维热线测量低速流场三维湍流度分布,采用一维热线连续变过热比测试方法测量马赫数为0.2~1.5可压流场,获得了试验段跨声速流场湍流度及风洞关键部段的整流特性。
1 风洞设备及测试仪器
1.1 风洞设备
本文试验在中国空气动力研究与发展中心的0.6 m连续式跨声速风洞中进行。热线测点分布如图1所示,包括换热器入口(前测点A)、换热器出口(后测点B)、第四拐角段入口(测点C)、稳定段入口(测点D)、稳定段出口(测点E)和试验段核心流(测点F)。
换热器前测点A,后测点B:换热器采用椭圆翅片管型式[7];第四拐角段入口测点C:拐角段设计采用棉式阻性消声器方案[3];稳定段前测点D,后测点E:稳定段安装1层蜂窝器、3层阻尼网,可增加至5层阻尼网 [3];收缩段收缩比为12;试验段测点F:试验段上下壁板采用低噪声槽壁[3],开闭比为6%。试验段试验风速范围覆盖马赫数为0.2~1.5,试验过程为常压运行工况,试验段总压为100 kPa。
低速回流段热线测量采用L型支架,该种支架安装简便可靠,支架安装于洞壁,测点距离洞壁约60 cm。由于回流段流速较低,试验过程中无支架振动影响。典型实验现场稳定段出口测点如图2所示。
试验段试验现场如图3所示,热线测量支架安装于弯刀机构上,热线探头位于支架前端,是弯刀机构的旋转中心位置,处于试验段核心流区域内,测量支架稳定,在试验过程中未观察到支架振动。
1.2 热线风速仪
实验使用丹麦丹迪公司的Streamline热线风速仪系统,换热器入口至稳定段采用55P61二维探头进行测量,试验段采用55P11一维探头进行测量,系统配置的采集卡分辨率为16位。数据采样频率为20 kHz,单点采样时间为5 s。Streamline系统为CTA热线,在低过热比状态下频率响应较低,在高过热比状态下频率响应较高。因此,在测量前对热线变过热比的增益和滤波参数进行调节,将不同过热比下的通道频响调节在10 kHz的一致水平。
连续变热线过热比的测量流程是:先通过程序控制系统测量热丝在当前流场的冷丝电阻;再根据冷丝电阻设置热丝过热比的桥电阻,依次设置0.1、0.2、0.3、0.4、0.5、0.6、0.7、0.8等8个过热比;最后进行数据采集。以单点采样5 s计,Streamline系统完成8个过热比设置和数据采集大约需要3.5~4 min。
探头热丝材料采用直径5 μm的镀铂钨丝,长度1.25 mm,电阻温度系数0.0036 K−1。钨丝强度高但是易氧化,镀铂可以降低钨丝的高温氧化效应,保证热丝在较长工作时间内的稳定性。
重点对探头焊点进行了加固,以提高探头在高速流场中的适应性。图4为高速流场中吹断了热丝的55P11的2个针头,探头左侧针头上的热丝被完全吹掉,右侧针头还留有一段残留热丝。一般情况下,若针头焊点牢固,在高速流场中的热丝会在气动载荷较大的中间某点断裂,并且热丝会在针头留下残丝,如图4中的右侧针头;若针头焊点不牢,则热丝不会在针头留下残丝,如图4中的左侧针头。热丝极细,工艺上是用瞬间产生的电弧将热丝熔焊于针头,焊接时电弧过大或者过小、环境湿度的不同,都可能导致电弧焊接不牢,使得热丝在较高风速中首先在焊点处发生脱落,这是许多原厂的丹迪热线探头容易在高速流场中吹断的主要原因。为此,对原厂的55P11探头进行了重新焊接加固,加固后的探头能够适用于本文高速流场。
2 热线测量和数据处理方法
连续式跨声速风洞存在低速不可压流场和高速可压流场。针对低速不可压流场,采用单过热比测试方法和二维热线探头进行测量;针对高速可压流场,采用变过热比测试方法和一维热线探头进行测量。
2.1 低速不可压流体单过热比测试方法
单过热比测试方法通常用于低速不可压流体测量,在测试过程中保持热丝工作温度不变,通过热线输出电压与流场速度的线性关系来计算流场参数。
换热器入口至稳定段出口的流场是常温低速运行,流体具有不可压缩性特征,热丝对流体温度、密度的灵敏度可以忽略不计,其主要对速度敏感,因此,热丝固定一个高过热比的工作状态就可以完成测试过程。此时,热线输出电压与流场速度满足King公式:
$$ {U^2} = A + B \cdot {u^n} $$ (1) 式中:$ U $为热线风速仪输出电压,V;$ u $为流向速度,m/s;A、B为热线的校准系数;。
流场流向湍流度计算为:
$$ \varepsilon = \frac{{{u_\sigma }}}{{\bar u}} \times 100\% $$ (2) 式中:$ {u_\sigma } $为沿流向速度标准差;$ \bar u $为流向平均速度。
采用多维探头测量时,流场三维综合湍流度值计算为:
$$ \varepsilon = \frac{1}{{\bar u}}\sqrt {\frac{{u_\sigma ^2 + v_\sigma ^2 + w_\sigma ^2}}{3}} $$ (3) 式中:$ {v_\sigma } $ 是垂直径向速度标准差,$ w_\sigma ^{} $是水平径向速度标准差。
2.2 试验段可压流体连续变热线过热比测试方法
连续变热线过热比测试方法通常用于可压流体测量,在测试过程中连续改变热丝工作温度,每改变一次热丝工作温度采集一次数据。通过热丝在不同工作温度的一组测试数据来分析流场参数。
本文试验段流场属于跨声速范围,流体具有可压缩性,存在速度、密度和温度脉动,因此,热丝的相应灵敏度系数值不可忽略,热线输出电压与流体速度u、密度$\; \rho , $、总温T0有关:
$$ U=f(\rho \text{,}u\text{,}{T}_{\text{0}}) $$ (4) 由于$ m = \rho u $,将热丝对流体速度和密度的灵敏度统一考虑为质量流灵敏度,即
$$ U=f(m\text{,}{T}_{\text{0}}) $$ (5) 式中:m为质量流量。
热丝的质量流灵敏度系数和总温灵敏度系数可以通过热丝的热平衡关系[12]求解:
$$ \frac{{{U^{\text{2}}}}}{{{R_{\text{w}}}}} = \pi l{\lambda _0}({T_{\text{w}}} - {T_{\text{e}}})(A + B\sqrt {Re} )(1 - k{a_{\text{w}}}) $$ (6) 式中:$ {R_{\text{w}}} $为热丝的工作电阻,$ l $为热丝长度,$ {\lambda _0} $为气体导热率,$ {T_{\text{w}}} $为热丝的工作温度,$ {T_{\text{e}}} $为热丝在流场中的非工作温度,Re为雷诺数,$ k $为热丝变过热比响应系数,$ {a_{\text{w}}} $为热丝过热比。
由雷诺数的定义,有:
$$ {{\rm{Re}}} = \frac{{md}}{\mu } $$ (7) 式中: d为特征长度,$ {\mu _{}} $为流体动力黏性系数。
将式(7)代入式(6),并对公式两边取自然对数后求偏导,整理可得:
$$ \begin{gathered} \frac{{\Delta U}}{U} = \frac{{\Delta m}}{m}\left[ {\frac{{B\sqrt {{{\rm{Re}}} } }}{{4(A + B\sqrt {{{\rm{Re}}} } )}} - \frac{{k{a_w}}}{{2(1 - k{a_w})}}\frac{{\partial \ln k}}{{\partial \ln {{\rm{Re}}} }}} \right] - \\ \frac{{\Delta {T_0}}}{{{T_0}}}\left\{ \begin{gathered} \frac{{{\alpha _ * }{R_ * }}}{{{R_e}}}\frac{{\eta {T_0}\left[ {1 - k{a_w}({a_w} + 2)} \right]}}{{2{a_w}(1 - k{a_w})}} - \\ 0.38\left[ {{\text{1}} - \frac{{B\sqrt {{{\rm{Re}}} } }}{{2(A + B\sqrt {{{\rm{Re}}} } )}} + \frac{{k{a_w}}}{{1 - k{a_w}}}\frac{{\partial \ln k}}{{\partial \ln {{\rm{Re}}} }}} \right] \\ \end{gathered} \right\} \\ \end{gathered} $$ 即:
$$ \frac{{\Delta U}}{U} = {F_{CTA}}\frac{{\Delta m}}{m} - {G_{CTA}}\frac{{\Delta {T_0}}}{{{T_0}}} $$ (8) $$ {F_{CTA}} = \frac{{B\sqrt {{{\rm{Re}}} } }}{{4(A + B\sqrt {{{\rm{Re}}} } )}} - \frac{{k{a_w}}}{{2(1 - k{a_w})}}\frac{{\partial \ln k}}{{\partial \ln {{\rm{Re}}} }} $$ $$ {G_{CTA}}{\text{ = }}\frac{{{\alpha _ * }{R_ * }}}{{{R_e}}}\frac{{\eta {T_0}\left[ {1 - k{a_w}({a_w} + 2)} \right]}}{{2{a_w}(1 - k{a_w})}} - \frac{\omega }{2}{\text{(1}} - {\text{2}}{F_{CTA}}) $$ 式中:$ {F_{CTA}} $为恒温式热线的流量灵敏度系数,$ {G_{CTA}} $为恒温式热线的总温灵敏度系数,$ {\alpha _ * } $为热丝的电阻温度系数,$ {R_ * } $为热丝在参考温度时的电阻值,$ \eta $为可压流温度恢复系数。
对(8)式两边除以$ {G_{CTA}} $,再取平方,得到
$$ {\theta ^{{2}}} = {r^{{2}}}{\left\langle m \right\rangle ^{{2}}}{{ - 2}}r{R_{m{T_0}}}\left\langle m \right\rangle \left\langle {{T_0}} \right\rangle + {\left\langle {{T_0}} \right\rangle ^2} $$ (9) 式中:
$$ \theta = {{{\left\langle U \right\rangle } \mathord{\left/ {\vphantom {{\left\langle U \right\rangle } G}} \right. } G}_{{\text{CTA}}}} $$ $$ r = {{{{F_{{\text{CTA}}}}} \mathord{\left/ {\vphantom {{{F_{{\text{CTA}}}}} G}} \right. } G}_{{\text{CTA}}}} $$ $$ {R_{m{T_0}}} = {{\left\langle {m{T_0}} \right\rangle } \mathord{\left/ {\vphantom {{\left\langle {m{T_0}} \right\rangle } {\left\langle m \right\rangle \left\langle {{T_0}} \right\rangle }}} \right. } {\left\langle m \right\rangle \left\langle {{T_0}} \right\rangle }} $$ $$ \left\langle m \right\rangle = {( {{{\partial \theta } \mathord{\left/ {\vphantom {{\partial \theta } {\partial r}}} \right. } {\partial r}}} )_{r \to \infty }} $$ $$ \left\langle {{T_0}} \right\rangle = \theta ( 0 ) $$ $\left\langle x \right\rangle = \dfrac{{\sqrt {{{( {x'} )}^{\text{2}}}} }}{{\bar x}}$,$ x' $为物理量x的瞬态值与平均值之差,$ \bar x $为物理量x的平均值,$ \left\langle x \right\rangle $为物理量$ x $的脉动量。
基于公式(9),作$ \theta $值与$ r $值的拟合曲线,可以求解质量流脉动值$ \left\langle m \right\rangle $、总温脉动值$ \left\langle {{T_0}} \right\rangle $及相关系数$ {R_{m{T_0}}} $。
对$ m = \rho u $方程两边取自然对数后求偏导,可得:
$$ \frac{{\rho '}}{{\bar \rho }} + \frac{{u'}}{{\bar u}} = \frac{{(\rho u)'}}{{\overline {\rho u} }} $$ (10) $$ \frac{{\overline {\mathop {\rho '}\nolimits^2 } }}{{\mathop {\bar \rho }\nolimits^2 }}{\text{ + 2}}\frac{{\overline {\rho 'u'} }}{{\overline {\rho u} }}{\text{ + }}\frac{{\overline {\mathop {u'}\nolimits^2 } }}{{\mathop {\bar u}\nolimits^2 }} = \frac{{\overline {\mathop {(\rho u)'}\nolimits^2 } }}{{\mathop {\overline {\rho u} }\nolimits^2 }} $$ (11) 由一维等熵关系式,有:
$$ {T_0} = T\left( {{{1 + }}\frac{{\gamma {{ - 1}}}}{{{2}}}M{a^{{2}}}} \right) $$ (12) 式中:$ T $为流体静温,$ \gamma $为比热比,$ Ma $为马赫数。
由马赫数定义,有:
$$ Ma = \frac{u}{{\sqrt {\gamma RT} }} $$ (13) 式中:$ R $为气体常数。
由压力与密度的关系,有:
$$ P = C{\rho ^\gamma } $$ (14) 式中:$ P $为流体静压,$ C $为常数。
由理想气体状态方程,有:
$$ P = \rho RT $$ (15) 对公式(12)~(15)取对数求偏导后,代入公式(11),可得:
$$ \frac{{\overline {\mathop {u'}\nolimits^2 } }}{{\mathop {\bar u}\nolimits^2 }} = \frac{1}{{\mathop {(\alpha + \beta )}\nolimits^2 }}\left[ {{\alpha ^2}\frac{{\overline {\mathop {(\rho u)'}\nolimits^2 } }}{{\mathop {\overline {\rho u} }\nolimits^2 }} + 2\alpha {R_{m{T_0}}}\frac{{\overline {(\rho u)'{{T'_0}}} }}{{\overline {\rho u \cdot } {{\bar T}_0}}}{\text{ + }}\frac{{\overline {\mathop {{{{T}'^2_0}}} } }}{{\mathop {{{\bar T}_0}}\nolimits^2 }}} \right] $$ (16) $$ \alpha = {\left[ {1 + ( {\gamma - 1} )\overline {M{a^2}} /2} \right]^{ - 1}} $$ $$ \beta = \alpha ( {\gamma - 1} )\overline {M{a^2}} $$ 由式(16)可计算跨声速流场湍流度。当采用一维探头测量流场时,由于一维探头的热丝对垂直于热丝的流动敏感(即不仅对流向速度敏感,还同时对一个径向速度敏感),此时,由式(16)得到的流场湍流度是流向和径向的合速度的湍流度。
3 实验及测量结果分析
3.1 换热器入口至稳定段流场湍流度测试
采用55P61二维热线探头对回流段流场湍流度测量,探头迎来流安装,先测量流向(X向)和垂直径向(Y向),之后再旋转90°,测量流向(X向)和水平径向(Z向),从而获得X、Y、Z各向湍流度[20]。
换热器入口湍流度分布见图5,由图可见,流场湍流度量值较大,根据公式(3)计算该测点各流速状态下三维综合湍流度的平均值为22%。湍流度波动范围也较大,约50%(湍流度值最大60%,最小10%),湍流度值明显随着当地流速的增加而降低。流向湍流度值略大于径向湍流度值。
换热器出口(三拐入口)湍流度分布见图6,三维综合湍流度平均值为1.6%,比换热器入口的湍流度(22%)降低约93%,湍流度值明显随着当地流速的增加而降低。经过换热器整流后,流向湍流度值转变为略小于径向湍流度值。
四拐入口湍流度分布见图7,由于Z向穿插四拐声学实验,所以未测试。由图可见,经过流道流动发展,流向湍流度值转变为略大于Y向湍流度值,湍流度值略微随着当地流速的增加而降低,综合湍流度平均值为3.6%,比前一测点(换热器出口)的量值增加约1倍。
稳定段入口(四拐出口)湍流度分布见图8,流向湍流度值保持略大于径向湍流度值,湍流度值随着当地流速的增加而略微降低,三维综合湍流度平均值为16.7%,比前一测点(四拐入口)的量值增加约4倍。稳定段入口湍流度值虽然上升了,但是波动范围仅为4%(湍流度值最大约19%,最小约15%),明显小于换热器入口的湍流度波动范围(50%),说明流场稳定性提高了。
三层阻尼网条件下的稳定段出口湍流度分布见图9,经过稳定段整流,流向湍流度值转变为略小于径向湍流度值。湍流度值基本不跟随当地流速而变化,三维综合湍流度平均值为1.4%,比稳定段入口的湍流度(16.7%)降低约92%。
比较三层阻尼网与五层阻尼网条件下,稳定段出口湍流度测量结果见图10。在三层阻尼网条件下,稳定段出口湍流度均值为1.238%,增加至五层阻尼网后,稳定段出口湍流度均值为0.623%,湍流度下降约50%,而且湍流度值随当地流速变化而产生的波动更小,流场更为稳定。
风洞回路流向湍流度分布见图11。由图可见,连续式风洞的换热器和稳定段是关键的降湍部段,气流经过换热器和稳定段后,流场湍流度分别产生两次明显的下降。气流经过换热器后,湍流度第一次显著下降,但是,在四拐后,即稳定段入口前又恢复较高湍流度,从湍流度平均值的量级上看,换热器后与三层阻尼网稳定段出口的湍流度基本在同一量级,但是,由前文分析可见,稳定段出口流场的湍流度波动范围明显更小,比换热器后的流场更为稳定。因此,经过稳定段第二次降湍后,流场湍流度值不仅下降了,而且流场波动随当地速度的变化更小。
3.2 试验段流场湍流度测试
试验段可压流场采用连续变热线过热比方法进行测量。
在试验中首先要降低热丝断裂的概率,重点考虑的因素主要包括:气流中固体颗粒浓度、热丝强度、热丝焊点强度和热丝工作参数等。相对来说,连续式风洞的气源不如暂冲式风洞的气源干净,因此,试验前特别进行了洞体清洁,采用面团粘和酒精洗等处理方法,尽量减少风洞中的粉尘等固体颗粒。为了避免热丝焊点在高速气流中断裂,对购入的标准探头进行焊点加固。在满足频率响应要求的情况下,优化热丝变温范围和信号频率,信号通过频率范围0~20 kHz,低通滤波10 kHz。变热丝过热比0.1~0.8获得测量数据,根据式(9)作流场扰动图,根据式(16)计算流场湍流度。
在三层阻尼网条件下,风洞试验段典型扰动图见图12和图13。由图12可见,马赫数0.4流场扰动图呈现一阶线性特征,意味着式(9)的二次项变量对流场的扰动贡献较小,流场扰动由一种扰动变量为主导,其他扰动量与主导扰动量的相关性系数可以忽略不计,因此,扰动图中的θ值与r值为一阶线性关系。由图13可见,在马赫数0.7,扰动图呈现双曲线特征,这时,式(9)的二次项变量不可忽视,流场中至少有两种扰动量的量级接近,扰动量的相关性较强,无单一主导扰动量。
图12流场对应的热线脉动电压幅值谱,见图14,可见频谱总体上是从低频向高频呈现指数衰减,引起流场扰动的主要能量在70 Hz以下的低频区域。一般流场中的主要扰动量有速度脉动、温度脉动、密度脉动和静压脉动量,但是在马赫数0.4,流动特性为低可压缩,温度脉动和密度脉动极低可忽略,考虑频谱贡献量以低频为主的特征,可判断流场扰动主要为试验段下游向前传的低频压力脉动。从文献[3]对本文风洞噪声测量的结果表明:在低马赫数范围(0.2≤Ma<0.5),二喉道节流调节范围有限,下游噪声(压缩机、再导入段等)前传导致试验段压力脉动较大,从而验证了本文测量结果的合理性。
风洞在五层阻尼网条件下,试验段典型扰动图见图15和图16。由图15可见,在马赫数0.4扰动图也呈现一阶线性特征,流场速度脉动(湍流度)较低,流场扰动主要为低频静压脉动。由图16可见,在马赫数0.7扰动图也呈现双曲线特征。
图17是根据式(16)获得的试验段在三层阻尼网和五层阻尼网条件下,马赫数0.2~1.5流场湍流度测量结果。根据图17中的湍流度数据分别统计两种工况的湍流度平均值,可得:在三层阻尼网条件下,试验段马赫数0.2~1.5湍流度平均值为0.081%,增加至五层阻尼网后,试验段马赫数0.2~1.5湍流度平均值为0.067%,湍流度下降约17%。由图可见,流场湍流度在马赫数0.95附近出现一个峰值,这是由于流场速度接近马赫数1时,风洞流场较难稳定控制。湍流度在马赫数1.3以后逐渐递增,这是由于附面层噪声向试验段自由射流区域的发散可导致流场扰动量增加[17],超声速流场噪声随着马赫数增加而增加,从而导致流场湍流度增加。
由图13与图16可见,扰动图拟合曲线具有非线性特征,因此,为评估实测值与公式(9)预测值的拟合程度,按照非线性方程计算拟合优度G:
$$ G = 1 - {\left( {\frac{{\displaystyle\sum {{{\left( {z - {z^ * }} \right)}^2}} }}{{\displaystyle\sum {{z^2}} }}} \right)^{1/2}} $$ 式中:$ z $表示实测值;$ {z^ * } $表示预测值。
图12、图13、图15和图16数据的非线性拟合优度计算结果见表1,数据非线性拟合优度为0.969~0.987,说明实测数据与预测值吻合度较高。按照蒙特卡洛模拟方法[9],以实测值与预测值的差值为标准偏差,按正态分布随机产生1000组数据,进而计算湍流度的测量不确定度,见表1,蒙特卡洛模拟不确定度为0.003%~0.0014%,说明测量结果可靠性较高。
表 1 不同阻尼网层数试验段湍流度测量结果Table 1 Turbulence level result for different screen layers in test section阻尼网 马赫数 湍流度
<u>, %拟合优度 不确定度/
%三层 0.4 0.070 0.987 0.0012 三层 0.7 0.054 0.981 0.0003 五层 0.4 0.056 0.982 0.0014 五层 0.7 0.047 0.969 0.0013 4 结 论
本文采用恒温式热线风速仪完成了0.6m连续式跨声速风洞换热器入口至试验段流场湍流度测量,实验研究表明:
1)换热器段和稳定段是重要的降湍部段,均可降低湍流度90%以上。换热器段出口的湍流度平均值较低,换热器起到了初步整流效果,但是波动范围随当地流速变化较大。稳定段出口的湍流度不仅平均值较低,而且随流速波动范围小,流场稳定性较好。
2)稳定段阻尼网从三层增加至五层,可降低稳定段湍流度50%,可降低试验段湍流度17%。
3)采用CTA连续变热线过热比方法,可以测量获得试验段可压流场的扰动图和湍流度值。扰动图特征反映了试验段扰动模态,马赫数0.4的扰动图呈现一阶线性特征,这是由于流场湍流度较低,流场扰动以低频静压脉动为主;马赫数0.7的扰动图呈现双曲线特征;流场湍流度在马赫数0.95附近和马赫数1.3以后逐渐递增出现峰值。
-
表 1 声衬结构参数
Table 1 Summary of liner structure parameters
声衬编号 孔径d/mm 穿孔率σ/% 板厚t/mm 腔深D/mm L1 1 8.4 0.5 50 L2 1 8.4 0.8 50 L3 1 8.4 1.0 50 L4 1 19.6 1.0 50 L5 2 9.4 0.5 50 L6 2 9.4 0.8 50 L7 2 9.4 1.0 50 L8 2 19.6 0.5 50 L9 2 19.6 0.8 50 L10 2 19.6 1.0 50 -
[1] NAYFEH A H, KAISER J E, TELIONIS D P. Acoustics of aircraft engine-duct systems[J]. AIAA Journal, 1975, 13(2): 130–153. doi: 10.2514/3.49654
[2] CUMPSTY N A. Review—A critical review of turbomachinery noise[J]. Journal of Fluids Engineering, 1977, 99(2): 278–293. doi: 10.1115/1.3448745
[3] MAA D Y. Potential of microperforated panel absorber[J]. The Journal of the Acoustical Society of America, 1998, 104(5): 2861–2866. doi: 10.1121/1.423870
[4] SIVIAN L J. Acoustic impedance of small orifices[J]. The Journal of the Acoustical Society of America, 1935, 7(2): 94–101. doi: 10.1121/1.1915795
[5] INGARD U, LABATE S. Acoustic circulation effects and the nonlinear impedance of orifices[J]. The Journal of the Acoustical Society of America, 1950, 22(2): 211–218. doi: 10.1121/1.1906591
[6] INGARD U, ISING H. Acoustic nonlinearity of an orifice[J]. The Journal of the Acoustical Society of America, 1967, 42(1): 6–17. doi: 10.1121/1.1910576
[7] MELLING T H. The acoustic impendance of perforates at medium and high sound pressure levels[J]. Journal of Sound and Vibration, 1973, 29(1): 1–65. doi: 10.1016/S0022-460X(73)80125-7
[8] GUESS A W. Calculation of perforated plate liner parameters from specified acoustic resistance and reactance[J]. Journal of Sound and Vibration, 1975, 40(1): 119–137. doi: 10.1016/S0022-460X(75)80234-3
[9] YU J, RUIZ M, KWAN H W. Validation of Goodrich perforate liner impedance model using NASA langley test data[C]//Proc of the 14th AIAA/CEAS Aeroacoustics Conference (29th AIAA Aeroacoustics Conference). 2008. doi: 10.2514/6.2008-2930
[10] KOOI J W, SARIN S L, FOKKER B V. An experimental study of the acoustic impedance of Helmholtz resonator arrays under a turbulent boundary layer[C]//Proc of the 7th AIAA Aeroacoustics Conference. 1981. doi: 10.2514/6.1981-1998
[11] CUMMINGS A. The effects of grazing turbulent pipe-flow on the impedance of an orifice[J]. Acta Acustica United With Acustica, 1986, 61(4): 233–242.
[12] NARK D M, JONES M G, SUTLIFF D L. Modeling of broadband liners applied to the advanced noise control fan[C]//Proc of the 21st AIAA/CEAS Aeroacoustics Conference. 2015. doi: 10.2514/6.2015-2693
[13] NARK D M, JONES M G, SUTLIFF D L. Further development and assessment of a broadband liner optimization process[C]//Proc of the 22nd AIAA/CEAS Aeroacoustics Conference. 2016. doi: 10.2514/6.2016-2784
[14] SYED A A, ICHIHASHI F. The modeling and experimental validation of the acoustic impedance of multi-degrees-of-freedom liners[C]//Proc of the 14th AIAA/CEAS Aeroacoustics Conference (29th AIAA Aeroacoustics Conference). 2008. doi: 10.2514/6.2008-2927
[15] JONES M G, HOWERTON B M, AYLE E. Evaluation of parallel-element, variable-impedance, broadband acoustic liner concepts[C]//Proc of the 18th AIAA/CEAS Aeroacoustics Conference (33rd AIAA Aeroacoustics Conference). 2012. doi: 10.2514/6.2012-2194
[16] JONES M G, WATSON W R, NARK D M, et al. Evaluation of a variable-impedance ceramic matrix composite acoustic liner[C]//Proc of the 20th AIAA/CEAS Aeroacoustics Conference. 2014. doi: 10.2514/6.2014-3352
[17] JONES M G, WATSON W R, NARK D M, et al. Optimization of variable-depth liner configurations for increased broadband noise reduction[C]//Proc of the 22nd AIAA/CEAS Aeroacoustics Conference. 2016. doi: 10.2514/6.2016-2783
[18] BECK B S, SCHILLER N H, JONES M G. Impedance assessment of a dual-resonance acoustic liner[J]. Applied Acoustics, 2015, 93(4): 15–22. doi: 10.1016/j.apacoust.2015.01.011
[19] MEI J, MA G C, YANG M, et al. Dark acoustic metamaterials as super absorbers for low-frequency sound[J]. Nature Communications, 2012, 3: 756. doi: 10.1038/ncomms1758
[20] LEE J, SWENSON G W. Compact sound absorbers for low frequencies[J]. Noise Control Engineering Journal, 1992, 38(3): 109–117. doi: 10.3397/1.2827811
[21] TOYODA M, MU R L, TAKAHASHI D. Relationship between Helmholtz-resonance absorption and panel-type absorption in finite flexible microperforated-panel absorbers[J]. Applied Acoustics, 2010, 71(4): 315–320. doi: 10.1016/j.apacoust.2009.10.007
[22] SAKAGAMI K, MORIMOTO M, YAIRI M. A note on the effect of vibration of a microperforated panel on its sound absorption characteristics[J]. Acoustical Science and Technology, 2005, 26(2): 204–207. doi: 10.1250/ast.26.204
[23] LEE Y Y, LEE E W M, NG C F. Sound absorption of a finite flexible micro-perforated panel backed by an air cavity[J]. Journal of Sound and Vibration, 2005, 287: 227–243. doi: 10.1016/j.jsv.2004.11.024
[24] TAKAHASHI D, TANAKA M. Flexural vibration of perforated plates and porous elastic materials under acoustic loading[J]. The Journal of the Acoustical Society of America, 2002, 112(4): 1456–1464. doi: 10.1121/1.1497624
[25] TOYODA M, TAKAHASHI D. Sound transmission through a microperforated-panel structure with subdivided air cavities[J]. The Journal of the Acoustical Society of America, 2008, 124(6): 3594–3603. doi: 10.1121/1.3001711
[26] BRAVO T, MAURY C, PINHÈDE C. Enhancing sound absorption and transmission through flexible multi-layer micro-perforated structures[J]. The Journal of the Acoustical Society of America, 2013, 134(5): 3663–3673. doi: 10.1121/1.4821215
[27] BRAVO T, MAURY C, PINHÈDE C. Sound absorption and transmission through flexible micro-perforated panels backed by an air layer and a thin plate[J]. The Journal of the Acoustical Society of America, 2012, 131(5): 3853–3863. doi: 10.1121/1.3701987
[28] CHUNG J Y, BLASER D A. Transfer function method of measuring in-duct acoustic properties. I. Theory[J]. The Journal of the Acoustical Society of America, 1980, 68(3): 907–913. doi: 10.1121/1.384778
[29] CHUNG J Y, BLASER D A. Transfer function method of measuring in-duct acoustic properties. II. Experiment[J]. The Journal of the Acoustical Society of America, 1980, 68(3): 914–921. doi: 10.1121/1.384779