摘要:本研究利用 CFD 仿真技术,针对航空发动机涡轮凹槽叶片叶尖气膜孔,在气膜孔为圆柱体条件下,分析其对叶尖区域冷却效果和气动性能的影响规律。通过几何建模、网格划分、边界条件与求解器设置等步骤,进行计算与结果分析,并与实验数据对比验证。结果表明气膜孔布局参数对叶尖温度等有显著影响,敏感性分析为优化提供依据。本研究为叶片设计优化提供理论与技术支持,未来可深入多物理场耦合分析及探索新型结构布局。
关键词:航空发动机;涡轮叶片;叶尖气膜孔;CFD 仿真;冷却效果
目录
1. 研究问题简述
1.1 背景
航空发动机涡轮作为飞机的核心动力装置,其性能和可靠性直接关系到飞行安全与效率。叶片作为发动机中的关键部件,在高温、高压、高速的恶劣工作环境下承受着巨大的热负荷和机械负荷。特别是叶片叶尖区域,由于叶尖泄漏流等因素的影响,气流复杂且温度较高,容易导致叶片损坏。气膜孔技术是一种有效的叶片冷却方式,通过在叶片表面开设气膜孔,引入冷却空气形成气膜,对叶片进行保护。然而,气膜孔的布局对冷却效果和气动性能有着至关重要的影响,不合理的布局可能导致冷却不均匀、气动损失增加等问题。因此,深入研究航空发动机叶片叶尖气膜孔的冷却效果具有重要的工程实际意义。
1.2 目的
由于凹槽叶尖优异的气热性能,大量关于冷却的研究也是基于凹槽叶尖结构进行改进。本研究旨在利用计算流体力学(CFD)仿真技术,在气膜孔为圆柱体的条件下系统分析航空发动机涡轮凹槽叶片叶尖气膜孔对叶尖区域冷却效果和气动性能的影响规律,进而提出气膜孔边缘叶片及叶尖温度分布,以提高叶片的冷却效率,降低热应力,同时减少对气动性能的不利影响,为航空发动机叶片的设计和改进提供理论依据和技术支持。
1.3 意义
理论意义:丰富和完善航空发动机叶片冷却技术的理论体系,深入理解气膜孔布局与叶尖区域热流场、气动场之间的耦合作用机制,为气膜孔设计理论的发展提供新的思路和数据支持。
实际意义:通过研究气膜孔冷却效果,提高叶片的可靠性和耐久性,延长发动机的使用寿命,降低维护成本;同时,优化后的布局有助于提升发动机的整体性能,提高燃油效率,减少污染物排放,对航空业的可持续发展具有重要推动作用
1.4 工程应用场景
本研究成果可应用于航空发动机的设计与研发过程。在叶片设计阶段,为气膜孔的布局设计提供科学依据,指导工程师选择最佳的孔参数组合;在发动机性能优化过程中,帮助评估不同气膜孔布局对发动机整体性能的影响,辅助制定改进方案;在发动机运行维护方面,为故障诊断和性能监测提供参考,通过监测气膜孔的工作状态和叶片温度分布,及时发现潜在问题,确保发动机安全可靠运行。
2. 研究模型建模
2.1 几何建模
表一 叶片模型几何参数
参数 | 数值 |
叶片弦长C/mm 叶片间隙G/mm 叶片数目 叶片高度 H/mm | 42.36 1 70 34 |
根据表一叶片大小及参数,采用 SW软件构建航空发动机涡轮凹槽叶片几何模型,绘制叶片的三维轮廓曲线,包括叶型曲线、叶尖曲线以及根部曲线等。
图一 叶片SW草图
通过旋转、拉伸等操作生成叶片的实体模型,特别注重叶尖区域的建模精度,以准确模拟气膜孔的开设位置和形状。在建模过程中,考虑到叶片的周期性结构,选取单个叶片通道进行建模,并利用 SW 软件的阵列功能构建完整的叶片环模型,以便后续进行周期性边界条件的设置。
图二 涡轮叶片整体俯视图及局部叶片观察图
表二 对比试验气膜孔模型几何参数
气膜孔相关变量 | 直径 D/mm | 间距 G/mm | 位置 S/mm |
数值 | 1.2 | 1.5 | 3.54 |
同时,为了便于气膜孔的参数化设计和布局优化,根据表二将气膜孔的建模过程进行参数化处理,表二中横轴表示气膜孔的属性,纵轴表示研究气膜孔对应的实验编号,通过改变孔的深度、间距、直径参数变量,实现气膜孔布局的快速调整和优化。
图三 叶片叶尖气膜孔模型几何草图
2.2 简化假设与计算域
简化假设:
假设空气为稳态流动,即流场中的物理量不随时间变化,这在发动机稳定工况下是合理的近似。
将空气视为不可压缩流体,虽然在发动机实际工作过程中空气密度会有一定变化,但在本研究关注的叶尖区域局部范围内,该假设对主要流动和传热特性的影响较小。
采用理想气体模型,忽略空气的粘性和热传导等次要因素,简化计算过程,但在后续湍流模型和传热模型中会适当考虑这些因素的影响。
计算域选择:
计算域包括叶片外部燃气流域以及叶尖气膜孔区域,如图四表达流体域的入口与出口。叶片外部燃气流域用于模拟高温燃气绕叶片的流动,叶尖气膜孔区域则是研究的重点区域,精确模拟气膜孔内冷却空气的喷射和与燃气的混合过程。
在流体与固体交互作用模型方面,采用流固耦合模型,考虑叶片在热载荷和气动载荷作用下的变形对空气和燃气流动的影响。通过在叶片表面设置无滑移边界条件,即假设叶片表面速度与相邻流体速度相同,实现流体与固体的相互作用;同时,考虑叶片与冷却空气之间的对流换热,设置相应的边界条件。
图四 固体域及流体域整体图
3. 网格划分与网格无关性验证
3.1 网格划分
利用 ICEM CFD 软件进行网格划分。对于叶片内部冷却通道、叶尖气膜孔区域以及叶片外部燃气流域,均采用结构化网格进行划分,以提高计算精度和效率。在叶尖气膜孔附近区域,进行局部网格加密处理,确保能够准确捕捉气膜孔内冷却空气的喷射细节、气膜的形成和发展过程。网格大小根据叶片的尺寸、气膜孔的特征尺寸以及流场的变化剧烈程度进行调整,在气膜孔内部、叶尖前缘和后缘等流动变化剧烈的区域,设置较小的网格尺寸,以提高分辨率;在远离叶尖和气膜孔的区域,逐渐增大网格尺寸,以控制网格数量。网格类型主要为六面体网格,同时在一些复杂几何形状区域(如气膜孔与叶片表面的交界处)适当采用少量三棱柱和四面体网格进行过渡。
图五 叶片网格划分整体图及气膜孔划分放大图
图六 流体域及固体域总体划分网格图
3.2 网格无关性验证
为了验证网格的无关性,选取了三种不同网格密度的模型进行计算,分别记为网格 1、网格 2 和网格 3,其网格数量依次增加。在相同边界条件和气膜孔布局参数下,计算叶尖区域的温度分布、气膜冷却效率以及气动性能参数(如压力损失系数等)。通过对比发现,随着网格密度的增加,计算结果逐渐趋于稳定。网格 1 和网格 2 计算结果之间在温度分布和气膜冷却效率等关键参数上存在一定差异,尤其是在气膜孔附近区域;而网格 2 和网格 3 计算结果非常接近,误差在可接受范围内。因此,选择网格 2 作为最终的计算网格,其在计算精度和计算成本之间达到了较好的平衡。
3.3 网格质量检查
对选定的网格 2 进行网格质量检查。结果显示,网格的正交性良好,大部分网格质量分布在0.5至0.1,这有助于减少数值计算误差。网格的最小体积满足计算稳定性要求,不存在体积过小导致计算发散的网格单元。在叶尖气膜孔加密区域,网格密度和质量均能满足准确模拟气膜孔内复杂流动和传热过程的要求,保证了计算结果的可靠性。
图七 网格质量(上)正交性(中)长宽比(下)表格
4. 边界条件与求解器设置
4.1 边界条件设置
进口边界条件:
对于叶片外部燃气流域进口,根据发动机工况,设定进口速度为燃气的实际流速,方向垂直于进口截面。进口温度设定为燃气的高温温度,压力为发动机工作压力。
对于叶片内部冷却通道进口,根据发动机实际冷却空气供应情况,设定进口速度为相应的冷却空气流速,方向指向叶片内部。进口温度设定为冷却空气的初始温度,压力为冷却系统提供的压力。
出口边界条件:
叶片外部燃气流域出口采用压力远场边界条件,设定出口压力为大气环境压力,温度为大气温度,模拟燃气流出叶片通道后的远场环境。
壁面边界条件:
叶片表面(包括叶尖表面)与燃气接触一侧设置为对流换热边界条件,根据燃气的温度和对流换热系数,计算叶片表面的热通量。叶片表面与冷却空气接触一侧设置为耦合换热边界条件,考虑叶片与冷却空气之间的对流换热和热传导。
气膜孔壁面设置为耦合换热边界条件,同时考虑冷却空气在气膜孔内的对流换热和与叶片壁面的热传导。
4.2 求解器设置
使用 ANSYS Fluent 求解器进行仿真计算。求解器类型选择稳态求解器,因为在本研究中假设流动为稳态流动。对于湍流方程的求解,采用二阶迎风格式,以提高计算精度。压力 - 速度耦合采用 SIMPLE 算法,确保计算的稳定性和收敛性。收敛标准设置为残差小于 10⁻⁵,同时监测叶尖区域的温度、气膜冷却效率、压力损失系数等关键参数,当这些参数的变化幅度小于设定阈值时,认为计算收敛。
5. 计算过程与结果分析
5.1 计算过程
在仿真计算过程中,首先导入网格模型和边界条件设置文件,初始化流场参数和温度场参数。然后启动求解器进行迭代计算,在计算过程中实时观察残差曲线的变化情况,确保计算过程稳定收敛。对于收敛困难的情况,适当调整松弛因子、迭代步长等参数,优化计算过程。同时,通过监测叶尖区域关键物理量的分布情况,如温度、速度、压力、气膜厚度等,及时发现计算过程中可能出现的异常流动和传热现象,如气膜破裂、局部过热等,并进行相应的分析和处理。
5.2 结果分析
温度分布:从计算得到的叶片叶尖区域温度分布图可以看出,气膜孔的布局对叶尖温度分布有显著影响。在合理布局的情况下,冷却空气从气膜孔喷出后能够在叶尖表面形成连续、均匀的气膜,有效地将叶片与高温燃气隔开,降低叶尖温度。例如,当气膜孔数量适中、分布均匀且角度合理时,叶尖表面的最高温度明显降低,温度梯度较为平缓;反之,若气膜孔布局不合理,如孔间距过大或过小,会导致气膜覆盖不均匀,出现局部高温区域,增加叶片热应力,影响叶片寿命。
图八 叶片整体温度分布云图
气膜冷却效率:气膜冷却效率是衡量气膜孔冷却效果的重要指标。通过对比不同气膜孔布局下的气膜冷却效率计算结果发现,气膜孔的位置和角度对冷却效率影响较大。当气膜孔位于叶尖前缘附近且喷射角度与主流方向成一定夹角时,能够使冷却空气更好地附着在叶尖表面,形成稳定的气膜,提高冷却效率;而位于叶尖后缘或喷射角度不合适的气膜孔,冷却空气容易被主流冲走,气膜冷却效率较低。
图九 叶片叶尖温度分布云图
叶片叶尖整体温度分析:
叶片的根部(较粗的一端)温度较低,主要呈现蓝色和紫色,对应的温度大约在 700 到 1100 之间。
叶片的尖端(较细的一端)温度较高,主要呈现红色,对应的温度大约在 1500 到 1900 之间。
从根部到尖端,温度逐渐升高,颜色也从冷色调(蓝色、紫色)过渡到暖色调(红色)。
气膜孔温度分析:
靠近叶片根部的气膜孔温度较低,颜色偏蓝紫色,温度大约在 700 到 1100 之间。靠近叶片尖端的气膜孔温度较高,颜色偏红色,温度大约在 1500 到 1900 之间。
从这个温度分布图可以推测,叶片在工作过程中,尖端由于可能受到更高的热负荷,温度较高,而根部温度较低。气膜孔的温度分布也遵循这一规律,靠近尖端的气膜孔温度较高,靠近根部的气膜孔温度较低。这种温度分布对于叶片的材料选择、冷却系统设计等都有重要的参考价值。
图十 叶片叶身温度分布云图
叶身温度分布:
叶片的前缘(左侧)温度较低,主要呈现蓝色和绿色,对应的温度大约在 600 到 1000 之间。
叶片的后缘(右侧)温度较高,主要呈现红色和橙色,对应的温度大约在 1400 到 1800 之间。
从叶片的前缘到后缘,温度逐渐升高,颜色也从冷色调(蓝色、绿色)过渡到暖色调(红色、橙色)。
在叶片的中部和靠近后缘的部分,可以看到温度较高的区域,这些区域可能是由于燃气流动和叶片结构导致的热集中现象。
叶片前缘的温度相对较低,可能是由于冷却气流的作用。
5.3 可视化结果
利用Tecplot软件对计算结果进行后处理,生成了叶尖区域的温度云图、流线图等可视化图像。这些图像直观地展示了叶尖气膜孔在不同布局下冷却空气的喷射、气膜的形成和发展、与主流燃气的混合过程,以及对叶尖温度分布和气动性能的影响,为深入理解气膜孔布局的优化机制提供了有力支持。
图十 温度场及叶片横切切片分布云图
图十一 流体域及固体域温度整体分布云图
6. 结果验证与讨论
6.1 结果验证
将 CFD 仿真结果与实验测量数据进行对比。实验数据通过专门设计的叶尖气膜冷却实验装置获取,在实验中采用红外热像仪测量叶尖表面温度分布,利用压力传感器测量气动性能参数。对比结果显示,仿真计算得到的叶尖温度分布趋势与实验测量结果基本一致,在关键位置的温度值上也较为接近;对于气膜冷却效率和压力损失系数等参数,仿真结果与实验数据的误差在合理范围内,验证了本研究建立的 CFD 模型的准确性和可靠性。
6.2 误差分析
模型简化误差:在建模过程中,对叶片和叶尖区域进行了一定程度的简化,如忽略了叶片表面粗糙度、气膜孔加工误差等因素。这些简化可能会导致气膜孔内流动和叶尖传热模拟与实际情况存在偏差,例如表面粗糙度会影响边界层的发展,进而影响气膜的附着和冷却效果。
网格划分误差:尽管进行了网格无关性验证,但网格划分仍存在一定的近似性。在气膜孔与叶片表面的交界处等复杂几何形状区域,网格可能无法完全精确地描述实际几何形状,从而引入计算误差。特别是在气膜孔直径较小的情况下,网格尺寸的选择对计算结果影响较大。
湍流模型误差:采用的 k - ε 湍流模型在模拟叶尖复杂湍流流动和气膜孔内的流动时存在一定局限性。对于叶尖区域强烈的漩涡和非定常流动现象,以及气膜孔内的层流 - 湍流过渡区域,该模型可能无法准确捕捉流动细节,导致计算结果与实际情况存在一定差异。
边界条件误差:边界条件的设置是基于发动机实际工况的近似,如进口速度、温度和压力等参数的设定可能存在一定误差。这些误差会影响计算结果的准确性,特别是在边界条件变化较大的情况下,对叶尖区域的流动和传热特性影响更为明显。
6.3 敏感性分析
气膜孔数量的敏感性分析:分别改变气膜孔的数量,研究其对叶尖温度分布、气膜冷却效率和气动性能的影响。结果表明,随着气膜孔数量的增加,叶尖冷却效果逐渐提高,但气膜孔数量过多会导致气动损失显著增加。当气膜孔数量达到一定值后,继续增加孔数对冷却效果的提升不明显,反而会使气动性能急剧下降。因此,存在一个最优的气膜孔数量范围,需要根据具体的叶片设计要求和发动机工况进行选择。
气膜孔位置的敏感性分析:将气膜孔在叶尖弦向和展向位置进行调整,分析不同位置对叶尖冷却的影响。发现气膜孔在叶尖前缘附近布置时,对降低叶尖前缘温度效果显著;而在叶尖后缘布置时,对后缘的冷却效果较好。同时,气膜孔在展向位置的变化也会影响气膜的覆盖范围和均匀性,进而影响冷却效率和气动性能。
气膜孔角度的敏感性分析:改变气膜孔的喷射角度,研究其对冷却空气喷射轨迹和叶尖冷却效果的影响。结果显示,合适的喷射角度能够使冷却空气更好地附着在叶尖表面,形成稳定的气膜,提高冷却效率。当喷射角度过大或过小,冷却空气容易偏离叶尖表面,导致气膜破裂或覆盖不均匀,降低冷却效果。此外,气膜孔角度的变化还会对气动性能产生一定影响,需要综合考虑冷却和气动性能来确定最佳的喷射角度。
气膜孔间距的敏感性分析:调整气膜孔之间的间距,分析其对叶尖温度分布和气膜冷却效率的影响。当气膜孔间距过小时,相邻气膜孔之间的冷却空气相互干扰,影响气膜的稳定性和覆盖效果;而气膜孔间距过大时,会导致叶尖部分区域无法得到有效冷却,出现局部高温。因此,需要找到一个合适的气膜孔间距,以实现均匀、高效的叶尖冷却。
7. 结论与展望
7.1 结论
本研究通过 CFD 仿真技术深入研究了航空发动机叶片叶尖气膜孔的优化布局问题,建立了合理的几何模型、简化假设和计算域,进行了网格划分与无关性验证,并设置了准确的边界条件和求解器参数。通过计算结果分析,揭示了气膜孔布局参数(数量、位置、角度、间距等)对叶尖温度分布、气膜冷却效率和气动性能的影响规律,发现了当前气膜孔布局中存在的问题。通过与实验数据对比验证了模型的准确性,并对误差来源进行了详细分析。敏感性分析为气膜孔布局的优化提供了重要依据。本研究为航空发动机叶片叶尖气膜孔的设计优化提供了理论指导和技术支持,具有重要的实际应用价值。
7.2 展望
未来研究可从以下几个方面展开:
多物理场耦合分析的深入研究:进一步考虑叶片在高温、高压、高速旋转等复杂工况下的多物理场耦合作用,如热 - 结构 - 流场耦合、气膜孔内的两相流(空气和可能存在的冷凝水)与传热耦合等。建立更精确的多物理场耦合模型,研究气膜孔布局在多物理场相互作用下对叶片性能的综合影响,为叶片的可靠性设计提供更全面的理论依据。
新型气膜孔结构和布局的探索:研究开发新型的气膜孔结构,如异形孔(非圆形孔)、微孔阵列、复合孔结构等,并探索与之相适应的优化布局方式。通过实验和 CFD 仿真相结合的方法,评估这些新型结构和布局在提高冷却效率、降低气动损失、增强叶片抗热冲击能力等方面的潜力,为航空发动机叶片冷却技术的创新提供新的思路和方法。
参考文献(此处参考国标2005)
[1] MIT和波音的年度综述:CFD在航空航天的技术现状[EB/OL]. [2024-12-02].
[2] 赵炜, 陈江涛, 肖维, 等. 国家数值风洞(NNW)验证与确认系统关键技术研究进展[J]. 空气动力学学报, 2020, 38(6): 1165-1172.
[3] 唐怡, 杨福军, 付眸, 等. CFD验证与确认数据库平台数据结构设计与应用[J]. 计算机测量与控制, 2023, 31(2): 173-177,236.
[4]王昊冉.气膜冷却涡轮的气热耦合数值研究[D].中北大学,2024.DOI:10.27470/d.cnki.ghbgc.2024.001652.
[5]王旭.涡轮叶片通道内部冷却及端壁气膜冷却的数值研究和多目标优化[D].中国科学技术大学,2022.DOI:10.27517/d.cnki.gzkju.2022.001870.
[6]田兴江.涡轮导叶非轴对称端壁造型方法与气膜冷却研究[D].南京航空航天大学,2018.DOI:10.27239/d.cnki.gnhhu.2018.000031.
[7]曾军,王彬,康勇.气膜冷却涡轮导向叶片流场数值模拟[J].燃气涡轮试验与研究,2006,(04):16-19+34.
[8]吴超林.涡轮叶片型面双出口孔射流气膜冷却特性研究[D].沈阳航空航天大学,2013.