共振腔吸声机理的数值模拟研究.pdfVIP

CSAA06．P．109 共振腔吸声机理的数值模拟研究 孙晓涛，李晓东 (北京航空航天大学能源与动力工程学院，北京100083) 摘要：由共振腔按一定规则排列的声衬被广泛应用于航守发动机的降噪中，所以研究，E振腔的吸声过程 在工程上有很重要意义。由于共振腔尺寸非常小．通过实验无法观测到其内部及附近的流动情况。理论 上的研究广泛采用假设，难于达到准确。计算气动声学为解决这一问题提供了一种手段。本文对不同入 射频率下的共振腔吸声过程进行了数值模拟。在此基础1．对共振腔的吸声机理进行分析，认为腔门处 粘性耗散和涡脱落现蒙是其吸声的形式。 关键词：声衬：共振腔；吸声：计算气动声学 1引言 噪声的危害不言而喻，而且噪声问题几乎涉及到我们生活的各个方面。气动噪声占了 相当大的比重，航空噪声就是其中之一。随着民用飞机数量的快速增加和普及，航空噪声 问题日益引起世界范围的关注。随着国际民航组织(ICAO)和各国民航组织(如：FAA 等)的噪声适航条例对民用航空飞机的要求日益苛刻，民用b机制造厂商目前已经把噪声 放在与性能、安全、可靠、工艺等同等重要的地位。从民用市场来看，为了提高自己的竞 争力，必须降低飞机的近场与远场噪声．以满足旅客的舒适性要求同时减少对环境的危 害。 如今，声衬被广泛用于各种商业飞机的喷气发动机中来降低风扇噪声。安装在航空发 动机管道中的声衬。经过适当调整，会对噪声抑制起到十分明显的作用。声衬多采用按一 定形式排列的亥姆霍兹共振腔，是由共振腔口形成的震荡喷流来形成阻尼作用的，这种喷 流是在声衬表而伴随着声扰动形成的压力波动，从而将声能转化成内能消耗。声衬中的共 振腔附近流场和声场非常复杂。目前对声衬的流场和吸声机理的研究丰要是通过半经验的 理论研究和基于大尺度模型的实验研究。理论．^的研究广泛采用假设，难于达到准确。由 于共振腔开口处尺寸非常小，但是大部分声能耗散发生在共振腔的开u处，所以需要基于 大尺度模型的实验研究。但是重要的雷诺数很难和真实共振腔相一致，这是目前实验研究 最大的缺点。 自上个世纪90年代以来，计算气动声学(Computational 发展，为解决这一问题提供了一种新的方法。CAA格式本身就具有低频散和低耗散特点， 所以适合解决声音的产生和传播问题。而且如果要获得相同的精度，CAA方法要比cFD方法 需要更少的网格点，这样就减少了计算量，缩短了计算时间，刊以更快更准的获得计算结 果。所以本研究采用计算气动声学的方法来研究共振腔的吸声机理。获得了包括腔口处的 共振腔整体的流动情况，从而可以很好的分析共振腔的吸声机理，为实际工程中更好的设 计声衬提供理论依据。 2主控方程 兆振腔尺寸非常小，在腔口粘性作用对流动有重要影响。所以主控方程采用可压 Navier-Stokes方程。二维可压雷诺平均Navier—stokes方程形式如下： 百Op+昙(瓦)=。 昙(赢)+毒(瓦巧)+瓦ap=一毒(瓦) 昙@)+毒(石砒)=一瓦詈+言l_(芑+Vj毒j-压 ，i：j于 E2南+翔2膏卜 勺=詈后岛叱w)【嚣+等上丝毛1 二三篡k2川菪氧3cqxk] q=0．09；el=1．44；t2=1．92；o-k=1．o；t=1．217 3数值方法 本文采用的数值模拟方法和程序是在已有研究工作[8儿9]基础上的进一步发展，下面 分别简要介绍所采用的空间格式、时问格式、人工粘性和边界条件。 3．1空间格式 在实际计算中，为求解偏微分方程，需要将对空间偏导项在空间上离散。研究采用频 散相关保持(DRP)格式剥空间进行离散[3]。这种格式是一种高阶有限差分格式，将传统 CFD方法中泰勒展丌与傅立叶一拉普拉斯变换相结合米对主控方程空间偏导项进行展开，通 过选定一个波长和频率范围使得差分方程和微分方程之间的频散联系差别尽量小，所以采 用较少网格点就可以达到较高精度。 3．2时间格式 在实际计算rrl，为求解偏微分方程，需要埘时问进行推进。研究采用低耗散低频散56 步龙格库塔(LDDRK56)方法[4]对时间进行推进。与传统龙格库塔法相比，这种格式不但 具备低耗散低频散特性，还允许采用较大的时削步长而不影响结果的精度，因此在计算气 动声学中得到了广泛应用。 3．3人工粘性 有限差分格式在高波数(短波)范围内与偏微分方程并不能满足很好的近似。当波长 只有几个网格宽