螺旋桨气动噪声传播的高精度数值模拟和降噪研究

1、报告正文参照以下提纲撰写，要求内容翔实、清晰，层次分明，标题突出。请勿删除或改动下述提纲标题及括号中的文字。（一）立项依据与研究内容（4000-8000字）：1项目的立项依据（研究意义、国内外研究现状及发展动态分析，需结合科学研究发展趋势来论述科学意义；或结合国民经济和社会发展中迫切需要解决的关键科技问题来论述其应用前景。附主要参考文献目录）；1.1项目的研究意义涡桨飞机是飞机中的重要类别之一，具有燃油利用率高，经济性好等优点1如ATR72-500的燃油消耗率为2.47L/100Km,相较于Boeing737-800和A320的燃油消耗率，可分别降低8.5%和5.7%，在短途支线客机、公务机和

2、大型运输机中有重要应用。与涡扇发动机相比,涡桨飞机的桨叶没有短舱的包裹,直接暴露在空气中,因此噪声更大,是主要缺点之一。近年来,随着航空运输的持续增长,民众对噪声污染的要求越来越高2,国外政府对适航噪声的要求也更加严格3,4。整个欧洲航空工业界决定,相对于2000年生产出的典型飞机,截止2020和2050年,新型飞机的外部噪声减少50%和65%。同时,世界上超过60%的机场,美国超过50个机场将噪声视为最重要的问题1,5。这些对我国自主研制的型号MA60和MA700等涡桨飞机（如图1和图2所示）,和其它在研型号提出了更高挑战,涡桨飞机的气动噪声问题是目前航空工业界亟需研究和解决的关键技术之一。

3、鉴于此,螺旋桨气动噪声的预测和降噪控制技术的研究迫在眉睫。图1MA60涡桨客机图2MA700涡桨客机涡桨飞机气动噪声的主要噪声声源包括：涡桨桨叶、机体（如起落架和增生装置）和机体表面的湍流产生的气动噪声,其中螺旋桨桨叶旋转产生的谐波噪声6是涡桨飞机气动噪声的最大噪声源。涡桨飞机的主要飞行时间是巡航阶段，螺旋桨产生的谐波噪声传播到机体表面，透过机体传到舱内对机舱内部的舒适性产生重要影响，更重要的是声载荷造成机体的振动，影响飞机其它部件的声疲劳和电子元件的寿命，对安全性有不可忽视的影响。鉴于此，螺旋桨气动噪声在机身表面的预测和控制至关重要，直接关系到涡桨飞机的舒适性、安全性和环保性，因此，螺旋桨气

4、动噪声的数值模拟意义重大。本项目针对我国涡桨客机研制中，预测和降低螺旋桨气动噪声在机身表面分布的关键问题，发展基于高精度数值格式和声源模型的螺旋桨气动噪声传播数值模拟方法，开展大规模、精细化数值模拟研究，旨在揭示机身表面螺旋桨气动噪声的传播机理和分布特征，探索螺旋桨飞机的降噪方法，为涡桨客机机身表面的降噪设计提供技术和理论支持。1.2.国内外研究现状单个螺旋桨的气动噪声频谱如图3所示，包含了主要的谐波噪声（峰值）和宽频噪声。螺旋桨的气动噪声主要来源于低频的谐波噪声，谐波噪声由桨叶和空气的周期性运动产生，出现在桨叶旋转频率（BladePassingFrequency,BPF）的整数倍。相对而言，

5、宽频噪声的影响较小，可以忽略。因此，预测和降低螺旋桨的气动噪声首先要对谐波噪声的产生和传播开展准确的预测。关于单独螺旋桨气动噪声的产生和在均匀流动中传播机理存在较多的研究，已经形成了成熟的理论和数学模型。但是，螺旋桨气动噪声在复杂构型（如真实机身）周围非均匀流动（边界层）中的传播机理仍是螺旋桨气动噪声研究的亟待解决的关键问题之一，这是由于存在以下难点：1）机身等固壁对噪声的反射、衍射和遮蔽效应是螺旋桨气动噪声传播的核心研究问题，对螺旋桨气动噪声传播的影响最大。由于固壁的作用，螺旋桨气动噪声在旋转平面具有很强的指向性，声强随不同角度相差约40-50dB,在复杂构型下准确捕捉到这些差异，研究噪声的

6、反射、衍射和遮蔽效应随声波波长的变化机理是具有较高难度；2）流场梯度对螺旋桨噪声的折射和反射作用是螺旋桨气动噪声传播的关键问题，很大程度上改变机身表面的声场分布。机身周围流动的梯度较大（如机身周围边界层），流场梯度对声场的反射和折射作用较强，准确预测声场在流场中的反射和折射机理，具有较高难度；3）多个螺旋桨间气动噪声的干涉作用是螺旋桨气动噪声传播的重要问题，较大程度上改变了机身表面的噪声分布。探索螺旋桨相位差对气动噪声传播的干涉机理对准确预测机身表面的噪声分布和后续的降噪设计具有重要的理论意义。Freqnancy(Hz)图3螺旋桨气动噪声频谱6gp而孟Isun需a)dPirlos1.2.1.

7、螺旋桨气动噪声的传播机理研究螺旋桨气动噪声在复杂构型(如真实机身)周围非均匀流动(边界层)中的传播机理可以通过试验、解析和数值模拟方法开展研究。Hanson7、Salikuddin8等人采用试验方法研究了螺旋桨噪声在机身表面的分布和散射问题。研究表明，机身对单个螺旋桨具有很强的遮蔽效应；当来流马赫数大于0.7时，边界层的折射作用较强地改变了机身表面噪声分布。Hanson9、McAninch10、Lull和Salikuddin8等人采用解析方法对边界层的反射和折射作用开展了研究，通过将机身简化成圆柱和简化、均匀的边界层，证明了折射作用主要来自边界层的速度梯度。试验方法需要在消声室内进行，该方法存

8、在费用高，周期长，滞后性等问题，并且使用的是缩比模型，难以兼顾气动和声学的相似性，需要对试验数据进行修正。解析方法成本低，周期短，通过理论推导得到明确的表达式，能够对各个参数的影响规律作出分析，但考虑的几何构型和流动情况简单，只能在已有的试验数据范围内进行研究，在条件以外公式会失效。当考虑几何构型复杂的工程问题时，数值方法是必需的，且具备下列优势：1) 可以模拟复杂构型，如考虑真实的机身几何外形；2) 背景流动真实，如采用CFD(ComputationalFluidDynamics)模拟的方法模拟真实边界层。近年来，随着计算机技术的飞速发展，计算气动声学(ComputationalAeroAc

9、oustics,CAA)得到了较大发展，并且度过了一个黄金周期12，在求解方法和以高精度差分格式为代表的数值技术方面取得了一系列重要进展12,13。计算气动声学采用数值模拟的方法，结合专门为声波传播定制的控制方程以及高精度、低耗散、低色散数值格式，求解声场的产生、传播和辐射。计算气动声学需要解析复杂构型（如起落架、缝翼等）非定常流动产生的湍流引起的声场，其计算量是巨大的。CAA的成本较低，周期短，能够兼顾气动和声学的相似性，针对复杂构型开展气动噪声研究，是低噪声设计的重要手段之一，是较先进的技术和未来气动噪声研究的发展方向，也是本项目拟使用的研究手段。1.2.2. 数值模拟方法目前，国内外大多

10、使用混合方法开展螺旋桨气动噪声传播的数值模拟研究，使用直接方法的研究很少，这是因为直接方法的计算量太大，计算周期过长。CAA混合方法根据声场的物理机制，将声场的物理过程分成三部分：声场的产生、近场传播和远场辐射，分别采用不同的控制方程和求解方法开展研究，大大减少了计算量和计算时间。针对螺旋桨气动噪声传播的数值模拟研究，混合方法首先开展单独螺旋桨的数值模拟或建立模型提取声源，通过扰动方程和高精度数值格式求解螺旋桨气动噪声的传播。准确预测螺旋桨气动噪声的传播和机体表面的分布，需要采用专门的数值方法，这是由于螺旋桨气动噪声传播的数值模拟技术存在以下难点：1）声波幅度小：气动噪声本质上是气体的压力脉动

11、，声场的幅度通常与流场也相差5到6个数量级14。传统CFD的格式截断误差大，耗散和色散大，不适合开展声传播数值模拟；2）传播距离远：相较于流体动力脉动（如涡脉动），声场的传播距离更远，通常大2到3个数量级，计算气动声学的计算量大；3）传播机理复杂：涡桨飞机巡航阶段的典型马赫数约为0.42，机身附近的流动复杂，流场对声场的反射和折射作用较强，需要准确模拟流场对声场的折射和反射作用；4）计算量大：由于噪声的本质是宽频的，人耳能感觉到的频率从20Hz到20000Hz,全机尺寸下解析高频噪声需要海量的网格量（109级别的网格量），气动噪声的计算量巨大；5）声源模型复杂：螺旋桨桨叶的运动较为复杂，在飞行

12、的过程中，桨叶的运动包含了旋转、向前平移等，声源模型需要准确模拟桨叶的空间位置；此外，包含了多种噪声机理和成分，如桨叶表面压力分布引起的载荷噪声，需要高精度的声源模型来给出；6)边界条件要求高：为了减少计算量，需要将物理域人为截断成尽可能小的物理域，同时会对数值模拟带来较大的影响。因此需要引入对声场反射少或吸收强的边界条件，提高数值模拟的精确性。因此，螺旋桨气动噪声传播的高精度数值模拟技术是当前美国机体噪声计算会议(BenchmarkproblemsofAirframeNoiseComputations,BANC)15和欧洲国际航空研究论坛(InternationalForumforAviat

13、ionResearchNoiseWorkingGroup,IFAR)的研究重点和热点之一，具有较高的难度和技术要求。传统声学使用波动方程来求解均匀来流下声波的反射和衍射现象，难以捕捉流场梯度对声波的反射和折射作用。目前国内外气动声学领域大多使用声扰动方程(AcousticPerturbationEquation,APE)或LEE结合高精度格式来开展声传播的数值模拟研究，包含流场梯度对声场的反射和折射作用。APE是一种声类比方法16，通过对LEE开展特征值分析，保留产生声波的源项，过滤掉产生涡波和熵波的源项，存在一定程度上的近似，能够使计算更加稳定，在德宇航(DLR)、南京航空航天大学和西工大的

14、研究中应用较多。典型工作包括：DLR的Ewert等人17,18采用色散关系保持(DispersionRelationKeeping,DRP)差分格式19离散APE,预测了二维翼型尾缘的噪声，分析了后缘噪声的主要来源，并通过在后缘打孔降低了噪声。南安普顿大学的Ma等20采用预系数化紧致格式21和APE，研究了声衬对前缘缝翼的降噪效果。南京航空航天大学的陈荣钱等人22,采用APE和DRP差分格式，研究了NACA0012尾缘噪声的产生和传播。西工大的白俊强研究组23进一步对30P30N二维增升装置的声传播开展了数值模拟。南京航空航天大学的司海青等人24对APE进行了改进，降低了计算成本，同时保留了计

15、算精度。LEE保留了声波、涡波和熵波的贡献，精度更高，在DLR、法宇航(ONERA)、南安普顿大学的空客噪声技术中心(AirbusNoiseTechnologyCentre,ANTC)、美国国家宇航局(NASA)兰利研究中心(LangleyResearchCentre,LaRC)和英国的声学和振动中心(InstituteofSoundandVibration,ISVR)的声传播的相关研究中使用,产生了大量的研究成果和数值格式，其中与本项目最相关的是螺旋桨气动噪声在非均匀流场中的折射和反射数值模拟研究。徳宇航的Siefert和Delfs25在马赫数为0.7,0.75和0.8时，使用LEE研究了不

16、同频率的静止点声源在机身表面有限区域的噪声分布；研究表明，边界层能够显著减少上游区域机身表面的噪声分布，巡航状态下边界层的折射和反射作用不可忽略，如图4所示。徳宇航的Dierke等人26扩展了之前的研究成果，采用涡扇发动机风扇噪声的模态替换了之前的点声源，使用相同的数值模拟方法在机身表面有限区域内研究了机身表面的噪声分布，得到了相同的结论。之后，Dierke27,28推导并开发了精度较高的螺旋桨声源模型，并采用LEE和DRP差分格式模拟了单个螺旋桨气动噪声的传播。针对声源模型和高精度数值求解方法，对比了试验和数值模拟结果，证明该螺旋桨气动噪声传播数值模拟方法具有较高的精度，如图5和图6所示。法

17、宇航的Garrec和Reboul29使用类似的声源模型和LEE预测了单侧开式转子的桨叶噪声在机尾表面的噪声分布，如图7所示。南安普顿大学的Miao30,31，使用LEE和相同的声源模型研究了单个螺旋桨噪声在圆柱表面的噪声分布，如图8所示。之后，Miao31研究了在起飞阶段，相同相位下，双螺旋桨气动噪声在ATR72机身表面的分布规律，如图9所示。国内关于螺旋桨气动噪声的研究有一些代表工作。乔渭阳等人32对螺旋桨厚度和载荷噪声的时间历程提出了一种新的求解方法，结果精度令人满意。刘沛清研究组33使用Kirchhoff-FWH方法对单独螺旋桨的噪声开展了计算，并对噪声频谱开展了分析。高永卫、乔志德等3

18、4改进了升力面方法，快速求解螺旋桨的远场声学特性，结算结果与试验结果符合较好。毛熙昌等35采用试验方法对比了四叶桨和三叶桨的气动和声学特性，结果表明四叶桨的气动性能更好，噪声更低。王同庆36对螺旋桨飞机缩比模型表面的声压开展了预测和试验对比结果显示预测方法与试验符合较好。刘秋洪37等人深入研究了气动噪声传播的物理模型和数值预测方法，提出了非紧致阻抗边界气动噪声辐射/散射统一积分计算方法,并通过圆柱绕流算例验证了方法的有效性。HvEl-5-ij5ndb)-eT1pIJOJi-4)图41000Hz下点声源在机体表面产生的声压级分布：a)无边界层；b)有边界层。201Load2rzO.fiLd2mi

19、c2PiAHO10100Dp_i窃604090SPL(im),dECarplete2HPF1Hh2Coimplffte2SPF2linsgTMeasurernefiiB-PF1#PaeasurememLBPF1hI!_I5J3O10001500200025003&00J0*1图5螺旋桨声源模型数值解与理论解对比图6螺旋桨气动噪声传播数值解与试验对比a)开式转子压力脉动的空间分布b)开式转子机尾表面声压级分布图7开式转子噪声在机尾表面的噪声分布a)螺旋桨噪声的空间声压分布b)圆柱表面的声压级图8螺旋桨噪声在圆柱表面的分布a)螺旋桨噪声的空间声压分布b)机身表面的声压级图9螺旋桨噪声在机身表面的分

20、布1.2.3.基于噪声干涉机理的降噪研究多个螺旋桨间气动噪声传播中的干涉作用能够在较大程度上改变机身表面的噪声分布，利用干涉现象能够对机体表面开展降噪研究。固定相位方法(Synchrophasingorphase-lockedpropeller)正是利用多个相同转速、不同相位螺旋桨间气动噪声的干涉原理，达到降低机体表面噪声的目的38,39。传统的螺旋桨气动噪声降噪手段，如桨叶外形优化、流动控制、降低螺旋桨和机身间的气动干扰等方法，可以从根本上降低螺旋桨的气动噪声，但是设计难度大，周期长，成本高，并且可能损失飞机整体的气动性能。安装声衬能够降低机身表面和远场观测点的噪声分布，但会增加飞机重量。并

21、且声衬的厚度需要与声波的波长相当才能起到较好的吸声效果38,39，螺旋桨噪声的主要来源是100Hz-250Hz的低频谐波，声衬不适合螺旋桨气动噪声的降低。相关研究表明40，使用固定相位方法对LockheedP-3飞机的试验数据进行分析，在某些机身观测点上噪声可以降低8到14dB。因此，固定相位是降低螺旋桨气动噪声的有效方法之一,对其开展深入研究，量化降噪效果是必要的。目前，国外对固定相位对螺旋桨飞机降噪的研究较多，普遍使用解析和试验两种方法。弗吉尼亚理工学院暨州立大学的Jones和Fuller41采用试验方法，采用圆柱代替机身，单极子声源代替螺旋桨噪声，在消声室内测量了固定相位方法对圆柱表面噪

22、声分布的影响，机身内部某些观测点上最大可降噪30dB。多伦多大学的Richarz和Dale42使用更真实的，缩比后的螺旋桨和简化后的机身，测量了双桨下固定相位对机身表面声压级(SoundPressureLevel,SPL)分布的影响，累计可减少1dB。美国联合技术公司的Magliozzi40对LockheedP-3飞机的试验数据进行分析，表明固定相位方法可以使机身表面SPL平均减少1.5dB。国内对固定相位的研究较少，南京航空航天大学的Huang39,43等采用试验和优化的方法对双桨下固定相位方法进行了研究，并对试验处理算法的精度做出了优化，减少了计算量。采用解析方法研究螺旋桨相位差对声散射影

23、响的工作较少。弗吉尼亚理工学院暨州立大学的Fuller38采用无限长圆柱模拟机身，简单偶极子模拟螺旋桨声源，建立了圆柱表面声压的解析表达式，并研究了固定相位方法的降噪效果，得出了最佳相位差随声波频率和观测点位置变化的结论。相对于试验和解析方法，采用数值模拟研究固定相位方法的研究很少。1.3.亟待解决的关键问题关于螺旋桨气动噪声传播的高精度数值模拟方法和降噪研究，综合国内外研究现状，存在以下5个特点：1) 流场梯度对螺旋桨气动噪声折射作用的数值模拟研究较少，多数研究开展流场对简单声源的折射和反射研究，或是忽略流场梯度对螺旋桨气动噪声传播的影响；2) 螺旋桨气动噪声干涉现象的研究较少，多数螺旋桨气

24、动噪声数值模拟只研究了单桨的声传播，回避了双桨气动噪声的干涉问题；3）基于噪声干涉现象的固定相位方法的降噪研究，多数使用试验和解析方法使用数值模拟方法的研究较罕见；4）在螺旋桨气动噪声传播的数值模拟方法中，关于数值求解方法和数值格式的研究较多，螺旋桨声源模型的研究较少；高精度的、准确的螺旋桨气动声源模型对于整个螺旋桨气动噪声数值模拟方法是关键技术之一；5）关于声传播的数值模拟研究，大多停留在二维构型上的计算，三维复杂构型的数值模拟研究较少。综上所述，虽然螺旋桨气动噪声传播的高精度数值模拟技术在求解方法、数值格式和声源模型等研究点上已经取得了一定的突破，但将以上研究点综合起来，形成准确、可靠的螺

25、旋桨气动噪声传播的数值模拟程序，研究螺旋桨气动噪声在复杂流动中的传播机理，进一步使用干涉现象作为降噪的手段，对三维全尺寸复杂构型开展精细化降噪研究却很罕见。因此，基于高精度螺旋桨气动噪声传播的数值模拟技术是当前涡桨飞机噪声数值模拟和降噪设计中亟待解决的关键问题之一。1.4.本研究的总体构想本项目针对我国涡桨飞机研制中预测和降低机体表面气动噪声亟待解决的技术难题，发展基于高精度数值格式和声源模型的螺旋桨气动噪声传播的高精度数值模拟方法，针对螺旋桨气动噪声在机身周围的传播机理开展大规模、精细化的数值模拟。一方面，针对螺旋桨气动噪声在复杂流动中的传播机理，发展基于LEE、高精度声源模型和高精度差分格

26、式的螺旋桨气动噪声传播数值模拟方法，准确预测机体表面的噪声分布；另一方面，基于螺旋桨噪声传播过程中的干涉机理，采用固定相位方法来降低机身表面的噪声分布，分析机体表面降噪效果随相位角变化的规律，得出最佳固定相位角，量化降噪效果。本研究旨在研究螺旋桨气动噪声在机身周围复杂流动中的传播机理，基于传播机理探索有效的机身表面降噪方法，为涡桨飞机的降噪设计和噪声防护提供理论支持和研究方法。1.4参考文献1 .A.Ghobbar,andC.Friend.EvaluationofForecastingMethodsforIntermittentPartsDemandintheFieldofAviation:a

27、PredictiveModel,Computers&OperationsResearch,vol.30,No.14,pp.2097-2114,2003.2 .N.E.Antoine,andI.M.Kroo.AircraftOptimizationforMinimalEnvironmentalImpact,JournalofAircraft,vol.41,No.4,pp.790-797,2004.3 .P.Arguelles,M.Bischoff,P.Busquin,B.Droste,S.R.Evans,W.Kroll,J.Lagardere,A.Lina,J.Lumsden,andD.Ranq

28、ue.EuropeanAeronautics:AVisionfor2020,EuropeanCommission,2001.4 .Flightpath2050EuropesVisionforAviation,ReportoftheHighLevelGrouponAviationResearch,EuropeanCommission,2011.5 .A.Filippone.AircraftNoisePrediction,ProgressinAerospaceSciences,vol.68,pp.27-63,2014.6 .H.H.Hubbard.AeroacousticsofFlightVehi

29、cles:TheoryandPractice.Volume1.NoiseSources,NASA-RP-1258-VOL-1,NASALangleyResearchCenter;Hampton,VA,1991.7 .D.Hanson,andB.Magliozzi.Propagationofpropellertonenoisethroughafuselageboundarylayer,JournalofAircraft,vol,22,No.1,pp.63-70,1985.8 .S.M,T.A.M.C,B.R,andG.J.Anexperimentalandtheoreticalinvestiga

30、tionofthepropagationofsoundwavesthroughaturbulentboundarylayer,AIAA-1986-1968,1986.9 .D.Hanson.ShieldingofProf-Fancabinnoisebythefuselageboundarylayer,JournalofSoundandVibration,vol,92,No.4,pp.591-598,1984.10 .P.L.Spence.Effectsoffuselageboundarylayeronnoisepropagationfromadvancedpropellers,Journalo

31、fAircraft,vol,29,No.6,pp.1005-1011,1992.11 .H.Y.Lu.Fuselageboundary-layereffectsonsoundpropagationandscattering,AIAAJournal,vol,28,No.7,pp.1180-1186,1990.12 .S.K.Lele,andJ.W.Nichols.ASecondGoldenAgeofAeroacoustics?,PhilosophicalTransactionsoftheRoyalSocietyA:Mathematical,PhysicalandEngineeringScienc

32、es,vol.372,No.2022,2014.13 .X.Li,M.Jiang,J.Gao,D.Lin,L.Liu,andX.Li.RecentAdvancesofComputationalAeroacoustics,AppliedMathematicsandMechanics,vol,36,No.1,pp.131-140,2015.14 .S.W.Rienstra,andA.Hirschberg.AnIntroductiontoAcoustics.EindhovenUniversityofTechnology,2003.15 .E.Manoha,andB.Caruelle.“Summary

33、oftheLAGOONSolutionsfromtheBenchmarkproblemsforAirframeNoiseComputations-IIIWorkshop”,AIAA2015-2846,2015.16 .R.Ewert,andW.Schroder.AcousticPerturbationEquationsBasedonFlowDecompositionviaSourceFiltering,JournalofComputationalPhysics,vol,188,No.2,pp.365-398,2003.17 .R.Ewert,andW.Schroder.OntheSimulat

34、ionofTrailingEdgeNoisewithaHybridLES/APEMethod,JournalofSoundandVibration,vol.270,No.3.pp.509-524,2004.18 .B.W.Fassmann,C.Rautmann,R.Ewert,andJ.Delfs.PredictionofPorousTrailingEdgeNoiseReductionviaAcousticPerturbationEquationsandVolumeAveraging,AIAA2015-2525,2015.19 .C.K.W.Tam,andJ.C.Webb.Dispersion

35、-Relation-PreservingFiniteDifferenceSchemesforComputationalAcoustics,JournalofComputationalPhysics,vol,107,No.2,pp.262-281,1993.20 .Z.Ma,andX.Zhang.NumericalInvestigationofBroadbandSlatNoiseAttenuationwithAcousticLinerTreatment,AIAAJournal,vol.47,No.12.pp.2812-2820,2009.21 .R.Hixon.PrefactoredSmall-

36、StencilCompactSchemes,JournalofComputationalPhysics,vol,165,No.2,pp.522-541,2000.22 .陈荣钱，伍贻兆，夏健基于SNGR方法的后缘噪声数值模拟，”计算物理，2011,28(05):698-704.23 .周旺仪，白俊强，余培汛，韩啸基于SNGR方法的二维增升装置噪声数值模拟,西北工业大学学报，2017，35(02):183-188.24 .司海青，王兵.气动声学计算中一种声扰动方程的改进，航空计算技术，2012，(42)05:1-3+8.25 .S.Malte，D.Jan，andC.Bastien.Refract

37、ionandScatteringinHighMachNumberBoundaryLayers，AIAA2011-2847，2011.26 .D.Juergen，E.Roland，D.Jan，S.Christian，andR.Marco.TheEffectofaBoundaryLayeronEngineNoisepropagatingtotheFuselageatFlightConditions，AIAA2013-2006，2013.27 .A.A.Rinie，P.P.Michael，B.Heino，D.Jan，andA.Daniela.InstallationEffectsofaPropell

38、erMountedonaHigh-LiftWingwithaCoandaFlap.PartI:AeroacousticExperiments，AIAA2014-3191，2014.28 .D.Juergen，A.A.Rinie，D.Jan，andE.Roland.InstallationEffectsofaPropellerMountedonaWingwithCoandaFlap.PartII:NumericalInvestigationandExperimentalValidation，AIAA2014-3189，2014.29 .G.ThomasLe，andR.Gabriel.Comput

39、ationalAeroAcousticsofCounterRotatingOpenRotorModelonRearFullScaleAirplaneinCruiseCondition,AIAA2012-2125,2012.30 .S.Miao,Y.Hou,andX.Zhang.EvaluationandComparisonofLinearizedEulerEquationsandEquivalentSourceMethod,The22ndInternationalCongressonSoundandVibration,Florence,Italy,No.169,2015.31 .S.Miao.

40、GPUBasedAeroacousticComputationwithPrefactoredCompactSchemes,PhDthesis,FacultyofEngineeringandtheEnvironment,UniversityofSouthampton,2017.32.乔渭阳,李文兰.飞机螺旋桨噪声的数值计算,推进技术,1992,05:44-49.33.刘振刚,刘沛清.螺旋桨的噪声计算,北京力学学会第12届学术年会.中国北京,2006,p.2.34.高永卫,乔志德.一种快速计算螺旋桨气动声学特性的数值方法,空气动力学学报,2006,24(02):187-193.35.毛熙昌,王界斌

41、,李红民.两个缩尺模型螺旋桨的气动与声学性能对比实验,航空动力学报,1991,6(04):320-324+375.36.王同庆.螺桨飞机缩比模型表面声压预测的验证试验,航空动力学报,1997,12(04):58-61+103-104.37.王芳,刘秋洪,蔡晋生.非紧致结构气动噪声辐射散射统一积分计算方法.航空学报,2013,34(11):2482-2491.38 .FullerCR.NoiseControlCharacteristicsofSyn-chrophasing,Part1:AnalyticalInvestigationJ.AIAAJOURNAL,1986,24(7):1063-106

42、8.39 .HuangXH,ShengL,WangYY.PropellerSyn-chrophaseAngleOptimizationofTurboprop-DrivenAircraft-AnExperimentalInvestigationJ.JournalofEngineeringforGasTurbinesandPower,2014,136(112606):1-9.40 .MagliozziB.SynchrophasingforCabinNoiseRe-ductionofPropeller-DrivenAirplanes,AIAA-1983-717R.Atlanta:AIAA,1983.

43、41 .FullerCR.NoiseControlCharacteristicsofSyn-chrophasing,Part2:ExperimentalInvestigationJ.AIAAJOURNAL,1986,24(8):1271-1276.42 .RicharzWG,DaleBR.ModelTestsofPhase-LockedPropellerInstallations,AIAA-1985-1R.Reno:AIAA,1985.43.HuangXH,WangYY,ShengL.SynchrophasingControlinaMulti-PropellerDrivenAircraftC.

44、Chicago:2015AmericanControlConference,2015:1836-1841.2项目的研究内容、研究目标，以及拟解决的关键科学问题（此部分为重点阐述内容）；2.1研究内容本项目针对涡桨飞机研制中机身表面螺旋桨气动噪声分布亟待解决的关键问题，通过发展螺旋桨气动噪声传播的高精度数值模拟方法，研究螺旋桨气动噪声传播的物理机制，探索噪声干涉效应对机身表面降噪效果的基本规律，获得双螺旋桨的最佳固定相位角和量化降噪效果。具体研究内容包括以下三个方面：2.1.1. 发展螺旋桨气动噪声传播的高精度数值模拟方法针对该数值方法的特殊性，发展一套完整的、基于高精度差分格式的数值模拟方法，

45、包括：1）LEE声传播模型，用来准确预测流场梯度对声场的反射和折射作用，固壁对声场的反射和衍射作用；2）高精度、低耗散和低色散的数值格式：空间离散使用6阶精度的预系数化紧致格式（PrefactoredCompactScheme,PSC），用来减少耗散、色散和网格量；空间滤波使用10阶显式滤波，用来过滤未能准确解析的高频数值伪波；时间推进采用4阶精度，低耗散、低色散、低存储的显式龙哥库塔格式（LowDissipationLowDispersionRunge-Kutta,LDDRK），增大时间推进步长，提高计算效率，减少耗散和色散；3）吸收边界条件：采用缓冲区边界条件，在显著减小物理域的同时，减少

46、截断边界带来的声反射，提高数值解的精度；4）旋转偶极子声源模型：旋转偶极子随时间做旋转运动，模拟桨叶的运动轨迹，通过相位差控制每个桨叶的相对位置；偶极子幅度代表桨叶的载荷噪声大小，准确描述桨叶的噪声源信息；5）分布式并行计算方法：基于本单位的大规模计算集群，使用消息传递接口（MessagePassingInterface,MPI）开展气动声学并行计算技术研究，提高计算效率，缩短计算时间。2.1.2. 研究螺旋桨气动噪声传播的物理机制基于多种流动和螺旋桨状态的变化，通过精细化数值模拟，分析螺旋桨噪声的传播机制，揭示机身表面螺旋桨气动噪声分布的基本规律。具体包括：1）通过均匀来流、圆柱形机身和单桨

47、噪声模拟，分析固壁对噪声的反射、衍射和遮蔽作用；2）通过真实边界层、圆柱形机身和单桨噪声模拟，揭示流动梯度对噪声的折射和反射作用；3）通过真实边界层、圆柱形机身和双桨噪声模拟，分析相干噪声的干涉作用，揭示不同相位差对机身表面噪声分布的影响。2.1.3. 分析复杂机身构型下双螺旋桨的噪声传播特征，量化揭示噪声干涉现象对机身表面噪声分布的影响和降噪效果。在真实的、较复杂的机身几何外形下，研究降低螺旋桨气动噪声的最优固定相位角，并对降噪效果做出量化评价。具体包括：1）基于圆柱型机身下双桨噪声干涉的研究结果，对ATR72飞机的噪声传播特征开展大规模模拟分析；2）对ATR72的最优固定相位较和降噪效果开

48、展精细研究，如图10所示，揭示噪声干涉现象对机身表面噪声分布的影响；3）建立机身表面降噪效果的量化评价方法，通过对声能量累加的方法对真实机身表面的降噪效果给出较全面的评价。图10机身构型计算示意图2.2研究目标本研究的研究目标为：1）针对涡桨飞机机身表面螺旋桨气动噪声分布的重要问题，开发基于高精度声源模型和数值格式的螺旋桨气动噪声传播数值模拟方法，开发精确、可靠的数值模拟程序，提高噪声传播数值模拟的准确度和可靠性；2）研究螺旋桨气动噪声的传播机理，分析机身对螺旋桨气动噪声的遮蔽和反射效应，流体梯度对噪声的折射和反射作用，以及多个螺旋桨间噪声的干涉效应，总结气动噪声在机身表面分布的基本规律；3）

49、研究复杂机身下双螺旋桨气动噪声的传播特征，基于不同相位差下，多个螺旋桨间气动噪声的干涉原理，找到最优的固定相位角，实现机身表面的最大降噪效果，为后续的降噪设计和防护提供理论与技术支持。2.3拟解决的关键科学问题1）准确可靠的螺旋桨气动噪声传播高精度数值模拟方法的实施；2）螺旋桨气动噪声传播的物理机制分析；3）复杂机身下双螺旋桨气动噪声的传播特征分析；4）最佳固定相位角下机身表面的降噪效果探索。3拟采取的研究方案及可行性分析（包括研究方法、技术路线实验手段、关键技术等说明）；3.1. 研究方案本研究拟采用理论分析和数值模拟结合的方法，从4个方面开展研究。首先采用理论分析推导螺旋桨气动噪声传播的高

50、精度数值模拟方法，包括控制方程空间格式、时间格式和边界条件，以及螺旋桨噪声源模型，包括桨叶运动轨迹和桨叶运动产生的噪声信息，作为本研究的理论基础。其次，以程序的形式实现螺旋桨气动噪声传播的高精度数值模拟方法，作为本研究的研究工具。再次，通过对圆柱形机身下，不同流动状态和螺旋桨气动噪声传播开展精细化数值模拟，研究螺旋桨气动噪声的传播机制。最后，在复杂机身下，基于多个螺旋桨气动噪声的干涉现象，对最优固定相位角和机身表面的降噪效果展开精细研究通过对声能量累加的方法做出较全面的评价。图11是本研究技术路线图。W5ST-R研无曲喜研究目标代功噪.由传播WrftM度數血喷拟方法为唧究挫供理论韮础为研丸提供

51、埋论絶础邮旋桨动嗓卢1专播的禹料度数恒模抵程序沟研咒挝俣数庇枝拟工具救值樓拟利理论分*1嫌琏桨乂厨嗓由隹播的物即机制-妙折怙播机制难确预测机体表由|嗥声仔布数值桎拉和理论労析固定和任厅址注机专衷面的降嗓效黑硏丸荒得阳蝉庞1ET动瞬南的怜播特性堆:干涉现號降咙图11本研究技术研究路线图1）螺旋桨气动噪声传播的高精度数值模拟方法该方法研究包括两个部分：分别为通用的气动噪声传播的高精度数值模拟方法，以及专门的螺旋桨声源模型。a）气动噪声传播的高精度数值模拟方法研究为了满足气动噪声传播的数值模拟技术的特殊性，本研究拟使用有限差分法(FiniteDifferenceMethod,FDM)来开展数值模拟研

52、究。相较于有限体积法或有限元法，基于有限差分法的高精度声传播数值模拟方法是发展最早的，也是分析最透彻和技术最成熟的方法。该方法具有数值格式种类多，数值性能高，计算量小，评价体系成熟和易实现等优点，本项目拟使用的数值技术如下：线化欧拉方程：本项目拟使用LEE来模拟声场在流体中的传播。传统声学使用波动方程对声场开展模拟，可以反映流场对声场的对流作用和固体对声场的反射作用。LEE相较于波动方程，可以考虑各点流场的非均匀性，模拟流场梯度对声波的折射和反射。特别是针对涡桨飞机的巡航阶段，马赫数在0.42左右，此时流场各点的非均匀性较强，边界层内流场梯度对声场的折射和反射作用不可忽略。如图12所示，在均匀

53、来流下，螺旋桨噪声以0角直线传播到机身表面观测点；当边界层存在时，流场梯度的折射作用改变了传播方向，声源以0R角传播到机身表面同一观测点。螺旋桨声源的强度随着角度而变化，因此折射作用会使同一个观测点上的声压级发生变化，需要使用LEE来准确模拟声场的折射。I冷直线時播角I沛m倍播屛、谥动方向b线传播路銭/折肘路统边界泾机身表面图12螺旋桨噪声折射原理图6阶精度PCS格式：该格式具有截断误差小、耗散小、色散小和易执行的特点，可以明显提高声波的解析能力，如图13所示，横轴代表声波真实的波数，竖轴代表数值格式解析的数值波数，红色符号代表各个数值格式对应的数值波数，即解析能力。一维情况下，PSC格式的数

54、值波数为1.38，2阶中心差分的数值波数为0.314，表示每个波长下分别需要3.7(PointsPerWavelength,PPW)和20.0个网格点，网格量约减少5.4倍，三维情况下，可减少约157.9倍，可以大大减少计算网格量。0.5truew购电numberideal 2ndordercentral 尸Cmordercentral*7pointDRPHi富口ncoffn|Mci: AHheroflaptmized4thorder图13典型差分格式的解析波数10阶显式滤波：每个数值格式对声波具有特性的解析能力，如图13所示，超过数值格式数值解析波数的声波(高频)会失真，称为数值伪波(Spu

55、riousWave)，需要使用数值滤波将其过滤或衰减，否则会污染数值解甚至造成求解失败。Hixon的6阶PSC的数值解析波数约为1.38,10阶显式滤波的数值滤波波数约为1.32，如图14所示，可与Hixon的6阶PCS配合4阶LDDRK格式：该格式适用于求解声传播等非定常现象，具有截断误差小、低耗散、低色散和低存储的特点，同时能够使用较大的时间推进步长，可以提高执行效率，减少误差。缓冲区边界条件：当人为将物理域截断成较小的区域时，截断的边界会对数值解带来较强的反射，需要引入吸收边界或无反射边界条件来消除反射本研究拟使用缓冲区(BufferZone)边界条件结合拉伸网格的方法来吸收或衰减声波，

56、如图15所示。缓冲区是吸收边界的一种，其基本原理是进入该区域的声波会被快速衰减成0，从而减少对缓冲区外数值解的反射，该方法具有直观、易执行和吸声性能好的特点。此外，声波的数值模拟需要相应的网格密度来支持；反而言之，当网格较疏时，声波会严重衰减，并且原本在较密网格上被准确解析的低频信号在网格较疏时会变成高频数值伪波，会被数值滤波过滤。因此，结合快速拉伸网格，缓冲区衰减声场的性能会进一步提升。图15缓冲区边界结合拉伸网格b)高精度的螺旋桨噪声源模型本研究拟通过LEE右端源项引入螺旋桨载荷噪声源模型，该模型通过随时间旋转变化的偶极子替代桨叶来实现。螺旋桨气动噪声的主要贡献来自于100Hz-250Hz

57、的低频噪声，声波波长为1.4m-3.4m，通常涡桨飞机的桨叶半径约2m,此时桨叶形状对声传播的影响较小，桨叶可以看作紧致声源，使用旋转偶极子代替桨叶的形状来简化螺旋桨声源模型。旋转偶极子的位置随时间变化，和桨叶的运动轨迹相同；声场幅度由桨叶表面的压力分布决定；引入相位差来控制桨叶间的相对位置，实现多个螺旋桨之间固定相位角的控制。2) 螺旋桨气动噪声传播的高精度数值模拟程序开发本研究拟在已有的方法基础上，开发螺旋桨气动噪声传播的高精度数值模拟程序。一方面，使用高精度程序模拟螺旋桨气动噪声在自由空间的传播，通过与Farassat1A的解析解作对比，验证整个程序和声源模型的正确性。另一方面，开展基于结构化网格的大规模并行计算技术研究，测量并优化程序的并行效率设计的程序流程图如图16所示。图16程序流程图3) 螺旋桨气动噪声传播的物理机制针对圆柱形机身，基于多种流动和螺旋桨状态的变化，通过精细化数值模拟，分析螺旋桨噪声的传播机制，揭示机身表面螺旋桨气动噪声分布的基本规律。在Ma=0.42，高度为8