14亿求解点!E3高压涡轮全环瞬态LES仿真

2026-01-30

航空发动机领域,高压涡轮作为核心动力部件,其流场特性的精准仿真直接影响着发动机的效率与可靠性的提升。NASA“高效节能发动机(E3)”项目系统公开了涡轮转子/导叶的详细几何、网格及试验数据,为全球发动机仿真技术提供了权威验证标准

长期以来,公开文献中对E3高压涡轮的数值验证,多局限于单通道稳态RANS仿真,而全环整周非稳态大规模计算,因求解规模与精度要求极高,始终是行业技术难题。

针对这一技术挑战,DIMAXER软件突破性地实现了双级高压涡轮(含冷气边界)14亿求解点全环大规模瞬态LES(大涡模拟),验证了其在旋转机械仿真领域的顶尖解决方案能力。这一突破依托DIMAXER四大关键技术

• 时空一体化高阶求解算法:兼顾高解析度与计算经济性,实现可负担高解析度求解;

• 先进GPU异构计算平台:高效并行计算架构,轻松处理十亿、百亿级求解规模,大幅缩短超大规模仿真周期和成本;

• 动静交界面高精度无损传递:消除冻结转子法、掺混面法的模型误差,还原动静部件间真实流场交互;

• 旋转机械专用前后处理模块:量身定制的专业化功能,提升复杂仿真全流程效率。

01 几何模型

E3双级高压涡轮参数(双级四排叶轮):

文字文稿2503 

 

02 网格划分

网格采用O4H拓扑进行正六面体结构网格划分,并对壁面附近网格进行加密,而后生成涡轮的全环整周网格。其中,壁面网格第一层高度3e-06m,y+≈20。DIMAXER基于局部通量重构算法能以低网格量配合高阶格式实现精细流场模拟:本次E3高压涡轮计算中,计算总网格量2251万,经四阶精度重构,等效求解点14亿,足以保证精度,同时减少计算资源消耗。

1

2

3

 

 03 计算设置

求解器:

DIMAXER软件 NS可压缩 K3求解器(四阶精度)

边界条件:

入口(total pressure total temperature inlet):总压344740Pa,总温709K

出口(pressure outlet):60980Pa;

壁面:绝热、无滑移

转子转速(rpm):8283r/min

本次计算,按照NASA文献位置在叶片表面分离出了叶片冷却孔结构,保证冷气孔在叶片的相对位置上与文献中给出的位置一一对应,并在计算中将这些冷气结构按照源项进行处理,根据试验测试数据给定其质量流量入口和入射角度。

文字文稿21034 

第一级静子叶片冷气孔位置

文字文稿21049 

静子叶片冷气源项位置

文字文稿21062 

第一级转子叶片冷气孔位置

文字文稿21077 

转子叶片冷气源项位置

03 计算结果

本算例采用40张NVIDIA 4090D显卡并行运算,3天时间即可完成高压涡轮转动一整圈的完整计算。

文字文稿21151 

1. 总体气动参数对比

DIMAXER软件的POST后处理内置了旋转机械专用模块,可快速提取计算得到的边界信息物理量,同时进行性能参数计算。

以下是DIMAXER计算的E3高压涡轮气动性能参数与试验数据的对比:

文字文稿21259 

偏差相对较小,其中流量偏差1%,总功率偏差2%,性能参数基本吻合。

2. E3涡轮内流场参数对比分析

  • DIMAXER计算的E3高压涡轮沿流向的气动参数分布情况与试验对比,二者吻合较好。

文字文稿21356 

沿流向温度分布

文字文稿21366 

沿流向静压分布

文字文稿21376 

沿流向流量系数分布

文字文稿21388 

沿流向总压分布

  • E3涡轮第一级叶轮和第二级叶轮沿展向的反动度分布情况与试验的对比:

文字文稿21434 

1级反动度沿展向分布对比

文字文稿21449 

2级反动度沿展向分布对比

DIMAXER计算的反动度在展向方向的变化趋势与试验基本符合,第二级叶轮的反动度计算值相对较低,整体偏差不大。

  • E3双级高压涡轮各排叶片进出口的相对气流角沿展向的分布趋势对比(静子进出口的气流角为绝对值,转子叶片进出口的气流角为相对值):

文字文稿21586 

第一级静子出口相对气流角沿展向分布

文字文稿21606 

第一级转子入口相对气流角沿展向分布

文字文稿21626 

第一级转子出口相对气流角沿展向分布

文字文稿21646 

第二级静子入口相对气流角沿展向分布

文字文稿21666 

第二级静子出口相对气流角沿展向分布

文字文稿21686 

第二级转子入口相对气流角沿展向分布

 

在第二级转子入口截面的气流角分布略微偏大,其余进出口截面的气流角分布与试验值较为吻合。

  • 叶轮进出口展向的相对马赫数分布情况:

文字文稿21779 

第一级静子出口

相对马赫数沿展向分布

文字文稿21800 

第一级转子出口

相对马赫数沿展向分布

文字文稿21821 

第二级静子出口

相对马赫数沿展向分布

文字文稿21842 

第二级转子出口

相对马赫数沿展向分布

 

可以看出,在第一级静子出口截面,相对马赫数略微偏大,其余截面的相对马赫数分布较为接近,主要的偏差位置发生于端壁处,可能是端壁二次流脉动影响了计算稳定导致。

以上分析表明,在E3高压涡轮计算流场中,DIMAXER结果与试验结果中的各个气动参数的分布基本吻合,验证了其高精度求解计算的准确性。

3. E3涡轮内流场分布云图分析

DIMAXER软件POST后处理内置了和截面功能,可快速对流场进行S1和S2流面的分析。以下是10%叶高和90%叶高的相对马赫数云图和静熵云图:

文字文稿22134 

B2B截面相对马赫数分布

文字文稿22149 

B2B截面静熵分布

Blade to Blade截面中马赫数流场显示,最高马赫数分布在流道喉部位置,分布与实际相符。静熵云图清晰展示了流道内的能量损失分布,包括叶片尾迹区域的流动损失,以及气流进入叶轮时产生的冲击损失。与时均结果相比,DIMAXER计算的流场显示出了更多的流动细节,如涡系的发展、范围,尾迹的传递等信息。

 

此外,从叶片表面温度、叶片表面压力及两级转子的叶顶泄漏流的流场变化情况,可清晰地看到高压涡轮转动时的卸压过程,以及叶顶间隙泄漏流随叶片转动形成的局部高速射流与叶片尾迹的相互作用。

  涡轮逐级压降过程与动叶叶顶间隙流动

 

涡轮叶片表面温度变化

从涡轮叶片表面的温度变化情况可以看出,冷气由冷气口流出后,叶片被冷却的范围随着叶轮转动而变化。

在旋转机械仿真领域,从单通道到全环整周、稳态到瞬态、千万级到十亿级求解规模,DIMAXER充分展现了大规模、高精度、高效率的技术优势。无论是应对湍流发展、动静交界等复杂场景,还是满足性能优化、精度验证等实际工程需求,DIMAXER都能提供稳定高效的精准解决方案,成为驱动航空发动机等高端装备研发迭代的核心技术工具。

在此过程中,我们感谢东方汽轮机有限公司的相关领导与技术人员给予的专业支持与协作。通过双方的紧密合作,DIMAXER的技术能力在工程实践中得到充分验证,也为后续旋转机械领域的技术攻关与产品迭代奠定了坚实的合作基础。期待未来持续携手,共同推动高端动力装备的技术突破与产业升级。

 

推荐文章

引领 CAE 进入高解析度系统级工业仿真新时代

预约试用

Copyright © 2026 深圳秩益科技有限责任公司苏ICP备2025208929号

网站地图 | 隐私政策

回到顶部