23ACTRANTM与ACTRANDGM互相验证二维数值模型的研究.docx
- 文档编号:28226941
- 上传时间:2023-07-09
- 格式:DOCX
- 页数:17
- 大小:961.32KB
23ACTRANTM与ACTRANDGM互相验证二维数值模型的研究.docx
《23ACTRANTM与ACTRANDGM互相验证二维数值模型的研究.docx》由会员分享,可在线阅读,更多相关《23ACTRANTM与ACTRANDGM互相验证二维数值模型的研究.docx(17页珍藏版)》请在冰豆网上搜索。
23ACTRANTM与ACTRANDGM互相验证二维数值模型的研究
ACTRANTM与ACTRANDGM的对比验证:
二维模型数
值研究
PaulTraub
MTU航空发动机公司
一、MTU航空发动机公司简介
MTU(Motoren-undIXirbinen-Uiiion)是全球领先的独立MRO供应商。
MTU航空发动机公司的前身为MTU慕尼黑公司。
MTU航空发动机公司是徳国领先的飞机发动机制造商。
它的客户是全球的发动机与工业燃气轮机领域的制造商及运营商。
MTU航空发动机公司主要致力于喷气发动机的子系统(例如压气机、涡轮等)的发展和生产,为军用和民用飞机发动机的提供技术支持。
图1MTU航空发动公司总部,徳国幕尼黑
MTU公司所参与的航空发动机项目主要涉及到航空发动机的压气机、涡轮这两大部件。
采用MTU航空发动机公司的部件的发动机型号如下;
低压级压气机:
JT8D
高圧级压气机:
PW600\1000G\800
高圧级涡轮:
PW2000>CF6,JT8D,GP7000
低压级涡轮:
PW4000Giowth\6000\1000G\300\500\2000,GP7000,JT8D,V2500
FMU
WJXG.
OAtYi
W・
I屮
<.uip
MTUsactivities
MTU同时还承担一些航空发动机的测试及其它部件的研发工作。
图2所示为MTU航空发动机公司的主要产品以及所参与的一些航空发动机项目示意图。
relatedtosub-systems
PtduUui
七刁・.!
0“
图2MTU的主耍产品及所参与航空发动机项目
二、航空发动机的噪声特点及降噪措施
來自于工程上不断增加的对静音航空飞行器的需求,给工程技术上带來很大的挑战。
现代航空发动机的设计需要考虑各个部件的噪声影响。
航空发动机的噪声主要來口于前传的风扇、压气机噪声和后传的风扇、燃烧室、涡轮以及喷流噪声,如下图3所示。
图3航空发动机的噪声示意图
图4为现代典型的涡扇发动机的各部分噪声示意图,从图中可以看到,现代的大涵道比涡轮风扇发动机中风扇的前传/后传噪声和喷流噪声占有主导地位,其次就是发动机内部的涡轮/燃烧室及床气机向外辐射的噪声。
日前在风扇的前传/后传噪声的控制方面有许多比较成熟的技术:
一方面可以从声源着手,根据气动声学的相关理论合理的选择风扇的转子、静子数目,从而使声波在叶片通过频率处满足“截止”条件,有效的降低声波的能量:
另一方面可以采取抑制声波的传播和向外辐射的途径。
在工程上一般采用合理的设计消声短舱來控制噪声的传播。
使用消声短舱来降噪的核心技术是如何设计合理、且有效的消声部件——声衬。
在声衬的研究方面国内外的工程技术人员都做了大量的工作。
同时,大涵道比涡轮风扇发动机的广泛使用有效地降低了喷流速度,使得喷流噪声得到了一定的控制。
因此在航空发动机噪声的诸多成分中,航空发动机排气过程中的涡轮噪声问题越來越受到重视。
本文的工作主要关注涡轮噪声辐射问题,如图5所示。
现代航空发动机设计过程中需要对各种土行工况卜的噪声进行评估,包括降落时进近和起0(例如£越一返场)工况条件下的噪声对机场环境的影响。
本文的工作主要是研究这两种匸况条件下的涡轮噪声的环境影响。
相关降噪技术的成功运用使得航空发动机进一步的降噪变得越來越困难。
对工程设计人员來说,需要有更加精确的声学设计工具来帮助进一步降低航空发动机的噪声水平。
图4现代航空发动机各部分噪声示意图
Numericalstudyfocusonturbinenoiseradiationproblem
ainbiait
图5涡轮噪声辐射问题示意图
Focus
ACTRAN是当今市场上最完善的声学模拟软件。
基于有限元和无限元方法的通用技术,ACTRAN提供了丰富的单元库、材料库、边界条件、求解配置和求解器。
ACTRAN在航空声学领域已经被证明是一种快速、精确的声学模拟软件。
三、ACTRAN两个求解模块结果相互对比的目的
本文将ACTRAN/TM和ACTRAN/DGM模块进行对比主要是基于以下的理由;
(1)ACTRAN/DGM是专门预测短舱排气噪声辐射的求解模块。
它的程序
主要是基于间断伽辽金方法(DiscontinuousGalerkinMethod)和4步
Rung-Kutta时间格式结合的方法求解线化欧拉方程(LEE)。
它可以模拟包含旋流和非等温条件下的短舱排气噪声传播问题。
在排气噪声问题中很大的一个难题是怎么处理声波通过比较大的速度和温度梯度剪切层的问题。
因为声源在这些剪切层的里会激发放大或者中性不稳定波,不稳定波在向下游传播的过程中将会导致声波的产生并辐射到远场。
在DGM模块中求解的是线化欧拉方程,可以很方便的处理这种剪切层的问题。
在排气噪声为题方面,ACTRAN/DGM的计算能力在MESSIAEN>TURNEX等项目上己经得到验证。
但是ACTRAN/DGM求解的线化欧拉方程是一种时域求解的方法,这种方法的最大弱点是计算时间太长。
因此在工程问题上我们希望有一种快速的、冋时求解精度可以接受的方法。
(2)ACTRAN/TM模块是专门求解旋转机械的噪声问题。
在航空发动机上主要是求解风扇的前传/后传声问题。
它可以处理包含平均势流和非平均温度场以及包含声衬边界的问题。
它可以在一次计算中同时加载多个声源激励,这中技术能力在工程上可以大大节省计算时间。
在短舱的排气噪声问题中,国外的研究者使用“Membraneelement”來定义剪切层,这种方法的求解精度取决于“Membraneelement”的形状是否和真实的剪切层一致。
从2009年以后的版本开始,ACTRAN/TM己经包含了“Moehring"求解器。
它在包含非均匀流动和温度场的声传播计算中有完美的表现。
尽管已经有些人使用“Moehring"方法做过一些研究,但都是在一些低马赫数(Ma<0.1)的情况下的
一些验证。
从实际应用的观点來看,如何用“Moehiing”方法來进行声波通过比较大的速度梯度和温度梯度的剪切层的传播模拟,在航空发动机涡轮噪声的研究中有很重要的意义。
本文同时采用ACTRAN/TM和ACTRAN/DGM模块求解排气中的涡轮噪声问题。
以专业求解排气噪声问题的软件——ACTRAN/DGM的结果作为参照,用它来验证ACTRAN/m在计算排气噪声问题的能力。
使用ACTRAN进行数值模拟的几个关键步骤
使用ACTRAN/TM模块中的“Moehring”求解器进行数值模拟的步骤如下;
1、分析模型特征(包括网格数据和基本的气动声学计算参数)
2、在CFD网格上使用RANS模型进行流场计算
3、使用ACTRAN/CFD用户界面将流场结果插值到声学网格上。
——背景流动的速度结果包含在“FLOW,数据块中。
——背景流动的压力和密度场结果包含在“ACOUSHC—HETROGEINITr数据块中。
4、ACTRAN参数设置(声学网格、插值后的速度场、边界条件等)
5、ACTRAN计算机后处理
操作流程如图6所示。
图6ACTRAN的计算流程图
五、模型特征
模型特征主要有以下几个部分:
•儿何及维数特征:
——排气噪声需要考虑现代发动机全尺度几何尺寸的模拟
——本文中的计算模型是一个二维轴对称模型
•工况点、频率及声模态
——考虑真实的降落时进近和起飞时飞越(返场)的工况条件
——在叶片的每个通过频率下,可以解析出有那些声模态是可以传播。
声模态是由最后的一级涡轮转子所产生的。
•平均流参数
——背景流场是采用RANS模型计算。
同时采用非均匀、无旋假设。
•网格参数
1、几何模型及维数
15m
图7计算域模型
图7为计算域的示意图c整个计算域总体长15米、宽5米。
包括内涵和外涵的喷流计算域以及外部声辐射的计算域。
Radiationdistance
图8远场声辐射示意图
图8为远场声辐射示意图。
在ACTRAN中离声源半径46米,0-150°方位内设置测点。
在这些测点处,可以提取计算结果來做声辐射的指向性结果图。
2、发动机的匸况、频率及声模态
Cut-OnSoundModes
Approach
BPFat2575Hz(kR~25)
Cutback
BPFat3683Hz(kR~35)
|m|=3
n=1,2,3
|m|=3
n=1,2,3,4,5
|m|=10
n=1,2,3
|m|=10
n=1,2,3,4
|m|=16
n=1,2
|m|=16
n=1,2,3,4
|m|=23
n=1,2,3
Specifiedforcrossvalidation
|m|=29
n=1
图9进近和返场工况下的主要传播模态及通过频率
根据发动机特定工况下的转速、涡轮叶片数目、流动参数等条件,可以根据气动声学相关理论得到主要传播模态及传播频率。
图9所示为传播模态及频率。
3、平均流参数(尾喷口处的参数,包括进近及返场工况)
MeanFlowProperties(atthenozzleducts5inlet)
Parameter
Approach
Cutback
(ambient)
0.230
0.278
(bypass)
0.395
0.677
Mfrb(core)
0.183
0.328
%[m/s]
346.0
344.9
CFan
350.3
349.2
Cturb
510.0
526.8
P„Ikg/rrP]
1.167
1.147
Pfon
1.226
1.393
Pturb
0.550
Q.526
图10流动参数
4、网格参数
MeshParameters(basics)
TriaB
Topology:
Finenews:
interpolationorder2(quadratic)
>4elementsper^vavelength*
JDopplercorrected
图11网格参数
本文的计算采用二阶的三角形非结构单元。
网格密度的定义根据经验在单位波长上至少要有4个单元。
单位波长是指考虑背景流动的波长,其定义为
1=(Co-U)If
式中Co为当地声速,u为背景流动速度,f为计算频率。
从式中可以看岀在背景流动比较强的地方及计算频率比较人的时候计算网格需要加密。
六、计算结果
计算结果主要包括以下儿个部分:
•流场结果。
采用RANS模型计算,并将其插值到声场网格上。
•声场结果
单模态的结果
——近场的云图结果(声压)
——远场辐射的结果(46米远场处的辐射声压级指向性图)
多模态禮加的结果
场辐射的结果(46米远场处的辐射声压级指向性图)
图12流场结果(插值到声学网格后的结果)
1、流场结果
背景流场计算,依赖于CFD计算网格。
如果流场不考虑粘性的话,可以使用ACTRAN中门带的程序來计算。
更加准确的流场结果可以使用其它主流的CFD商业程序来计算。
本文结果基于RANS模型,并将其插值到声场网格上。
如上图12所示的结果。
2、声场的结果
(一)声源为单模态的数值结果,降落时进近工况。
ACTRAN/™(Moehring)方法和ACTRAN/DGM方法结果对比。
图13ACTRAN/TM与ACTRAN/DGM求解模块的远场指向性结果对比(进近)
图14ACTRAN/TM与ACTRAN/DGM求解模块的公图结果对比(进近)
图13为ACTRAN两个求解模块的远场指向性结果对比图,。
从图中可以看到除了声模态(16,1)从大约70。
的位置开始由于两种方法不同所造成的差异外,其它模态的对比结果都比较好。
图14为ACTRAN两个求解模块的声圧云图结果。
对比的声模态为(10,1),两种求解方法的云图结果对比良好。
(二)声源为单模态的数值结果,起飞后返场工况。
ACTRAN/TM(Moeliring)方法和ACTRAN/DGM方法结果对比。
CotfeacM:
—0.1)
FMs<1loton
图15ACTRAN/TM与ACTRAN/DGM求解模块的远场抬向性结果对比(返场)
图16ACTRAN/TM与ACTRAN/DGM求解模块的云图结果对比(返场)
图15.16为两种结果对比结果,从指向性图及云图结果来看对比结果良好。
(三)声源为单模态的数值结果,起飞后返场工况。
ACTRAN/TM(Moeliring)方法和ACTRAN/DGM方法结果对比。
高阶径向或者周向模态。
*TDGK
图17ACTRAN/TM与ACTRAN/DGM求解模決的远场指向性结果对比(返场)
图17为高阶模态下两种求解模块的远场指向性结果对比。
结果显示从100。
位置开始出现一些结果差异。
在其它位置结果对比良好。
(四)多个模态叠加的结果(10,1-4),起飞后返场工况。
BPFToneatCutback
SuperpositionofModes(tn.n^lO.1-4)
—rrr_D3k
MTUNhebing
w
图18ACTRAN/TM与ACTRAN/DGM求解模決的远场指向性结果对比(返场)
图18为高阶模态下两种求解模块的远场指向性结果对比。
周向模态为10,径向模态从1到4的所有模态叠加。
结果显示从100。
位置开始出现一些结果差异。
在其它位置结果对比良好。
以上所有的计算的频率均为航空发动机在降落时进近及起飞后返场工况下的叶片通过频率。
七、结论
(1)ACTRAN/™(Moeliring)方法在涡轮噪声数值模拟方面的能力可以通过与ACTRAN/DGM的数值结果对比从而成功得到验证。
数值计算的结果达到了我们所预期的目标。
(2)TM/Moeluing方法在排气噪声声辐射问题方面是一个非常有价值的分析工具。
即使是考虑真实的发动机配置和流动参数,这种方法也体现了比较强大的数值模拟能力。
(3)TM/Moeluing方法适合于改进发动机(声学)的整个设计过程。
这种方法的物理模型假设是可以接受的,并且它的计算时间也是合理的。
致谢
(1)本文的数值研究受到(徳国)政府投资项目MASSIFEffekT(大量减少空中交通噪声一有效降低航空发动机噪声研究)的资助。
(2)同样要感谢以下作者:
SteeveLal丙ee(FFD,协助完成ACTRAN/DGM部分的计算匸作。
YvesBurgold(ATENA),协助完成CFD的计算网格并且提供CFD数值计算结果。
Dr.DominikBroszat,协调MASSIF项目下的MTU子任务工作。
Dr.KarlEngelandDr.AiidieasFiala>核査CFD的计算结果。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 23 ACTRANTM ACTRANDGM 互相 验证 二维 数值 模型 研究