首页/文章/ 详情

涡轮丨金属所:高推重比下航空发动机涡轮叶片单晶高温合金/热障涂层体系的应用研究进展

4小时前浏览2

随着先进航空发动机对于高推重比和高服役可靠性的追求,不断对涡轮叶片材料工艺选用提出更高的要求,采用以镍基单晶高温合金作为基体,叶身外表面涂覆热障涂层的涡轮叶片已成为先进航空发动机的典型特征。镍基单晶高温合金和热障涂层材料的研发、镍基单晶高温合金/热障涂层体系的考核评价是保障涡轮叶片工作安全性的关键,是当前航空发动机领域的研究热点之一。本文对涡轮叶片用镍基单晶高温合金、热障涂层材料的研究及应用情况进行了详细介绍,简要介绍了镍基单晶高温合金/热障涂层体系的考核评价及失效机理研究进展,并从单晶基体/金属底层界面匹配性、新型金属底层和陶瓷面层研发以及考核评价等方面阐述了涡轮叶片用镍基单晶高温合金/热障涂层体系的研究重点,以期为充分发挥镍基单晶高温合金/热障涂层体系使用潜力提供参考。

关键词:镍基单晶高温合金;热障涂层;评价;界面;失效机理

涡轮叶片是航空发动机的核心部件,其承温能力是提升发动机性能、效率、可靠性的关键技术指标,一般认为提高涡轮进口温度是提高航空发动机效率的有效方法,目前推重比为10的发动机涡轮进口温度已达1580~1650℃,这就导致涡轮叶片实际使用温度越来越高。涡轮叶片所需要的高温性能70%~80%依靠材料的改进,目前采用以镍基单晶高温合金作为基体,叶身外表面涂覆热障涂层的涡轮叶片已成为先进航空发动机的典型特征。镍基单晶高温合金因消除了晶界高温下的薄弱环节,且随着 Re元素含量的不断提高,工程应用的镍基单晶合金最高使用温度可达1100℃,热障涂层(thermal barrier coating,TBC)涂覆于涡轮叶片表面,能使涂层表面与涡轮叶片基体间产生超过100~170K的温降,同时具有优异的抗氧化腐蚀性能,为进一步提升涡轮叶片的承温能力,具有更高使用温度的镍基单晶高温合金和新型热障涂层材料的设计、制备是当前航空发动机领域的研究热点。值得一提的是,为了获取镍基单晶高温合金/热障涂层体系的相关性能数据以支撑其在航空发动机上的工程应用,尤其是在高代次单晶合金或新型热障涂层试制阶段,热震实验、抗CMAS腐蚀实验、燃气热冲击实验等试样级、零件级考核评价是必不可少的环节。除此之外,随着航空发动机试车状态的增大和试车时长的累计,涂覆热障涂层的涡轮叶片存在局部剥落和表面附着CMAS(CaO-MgO-Al₂O₃-SiO₂混合物)等问题,国内外相关研究人员从单晶基体/金属底层界面、金属底层/陶瓷面层界面以及CMAS侵蚀等方面开展了失效机理研究,为镍基单晶高温合金/热障涂层体系匹配性优化、性能提升设计提供指导,进而保证其工作可靠性。


本文详细介绍了镍基单晶高温合金和热障涂层的研究及应用,综述了镍基单晶高温合金/热障涂层体系考核评价方法以及失效机理研究进展,并对今后涡轮叶片用镍基单晶高温合金/热障涂层体系的研究重点进行了展望。

1

镍基单晶高温合金的发展及应用

镍基单晶高温合金具有优异的蠕变、疲劳、氧化及腐蚀抗力等综合性能,是目前制造先进航空发动机涡轮叶片的主要材料。为了满足高性能先进航空发动机的使用需求,各国十分重视镍基单晶高温合金的设计和开发,形成了各自的单晶合金体系,极大地促进了镍基单晶高温合金的发展。


自20世纪70年代成功研制出第一代单晶高温合金PWA1480,并成功应用于PW2037,JT9D-7R4等军民用发动机的涡轮叶片,引发了开发单晶高温合金的热潮,相继出现了与PWA1480性能水平相当的CMSX-2,SRR99,René N4,KC-26Y等多个第一代单晶高温合金,第一代单晶高温合金之间的主要区别在W,Mo,Ta等元素含量的不同,并且几乎完全弃用C,B,Hf等降低合金初熔点的晶界强化元素。目前,国外应用最成熟的是以CMSX-4,René N5和PWA 1484为代表的第二代镍基单晶高温合金,第三代镍基单晶高温合金René N6,CMSX-10,TMS-75等,其中报道较多的是CMSX-10(RR 3000),20世纪90年代即在国外多型发动机上进行了实验验证,已应用于航空发动机涡轮叶片。21世纪初,随着合金设计水平的提高和生产工艺的改进,日本国立材料研究所(NIMS)通过调整合金中Re和Ru含量,研发出TMS-138(第四代)、TMS-196(第五代)和TMS-238(第六代)等一系列综合性能优异的单晶合金,但仍处于实验室阶段,尚需解决铸造性能、组织稳定性、成本等一系列问题才可实现工程化应用。


值得注意的是,在单晶合金的发展中,Re是一个典型的合金元素,第二、三代单晶高温合金Re含量分别达到3%(质量分数,下同)和6%,显著提高了合金的高温强度和抗氧化性,但因为Re是强TCP相形成的元素,所以降低了合金组织的稳定性,从而影响合金的高温持久性能;另一方面,Re是一种非常昂贵的合金元素,无Re合金的成本至少要比含3%Re合金低50%。Ru是第四代单晶合金的代表性元素,其主要作用是抑制TCP相形成,提高组织稳定性,进而提高其他强化元素含量的调整空间以获得更高的蠕变强度。除此之外,早期的单晶合金去除C,B,Hf等晶界强化元素,以增大热处理温度范围,但随着涡轮叶片尺寸的增大以及内腔结构的复杂化,涡轮叶片精铸件的单晶完整性控制变得十分困难,往往存在小角度晶界等缺陷,降低单晶合金力学性能。为保证涡轮叶片工作可靠性及单晶涡轮叶片精铸件合格率,高代次单晶高温合金中开始限 量使用C,B,Hf等元素,并对C,B,Hf的作用机制开展了深入研究。


我国单晶高温合金的发展基本沿用国外研制思路,通过开展测仿和自主创新,成功研制出与国外PWA 1480力学性能水平相当的第一代单晶高温合金DD3,DD8等,并实现工程应用。20世纪90年代末,成功研制出自主知识产权为2%Re低成本的第二代单晶高温合金DD6,并试制了仿René N5的、DD5、仿PWA1484的DD412,仿ЖC32的DD32等第二代单晶高温合金。与国外应用情况相似,目前国内应用最广泛的是以DD6,DD5等为代表的第二代镍基单晶高温合金,DD9,DD10,DD90等第三代单晶高温合金也已开展涡轮叶片试制及实验验证,同时开展了含Ru四代单晶合金DD15,DD91的研制。


可以肯定的是,现阶段及未来一段时间内,镍基单晶合金仍将是制备先进航空发动机涡轮叶片的首选材料,应继续加强对高性能镍基单晶高温合金的设计研发,获取合金元素的相互作用机理,研发低成本、高性能新一代镍基单晶高温合金。

2

热障涂层的发展及应用

涡轮叶片热障涂层是将耐高温、低导热、抗腐蚀的陶瓷材料通过特殊工艺与镍基单晶高温合金基体相复合,可降低涡轮叶片基体温度,提高耐冲刷、抗氧化及热疲劳等性能。涡轮叶片热障涂层一般包括金属底层与陶瓷面层,其中金属底层主要作用为增强陶瓷面层和金属基体结合力,提高二者热膨胀系数匹配性,并提升基体抗氧化性,陶瓷面层起隔热和抗高温冲蚀作用。自20世纪70年代初进入工业应用以来,热障涂层受到广泛重视并得到迅速发展。


目前,国内外航空发动机涡轮叶片应用最广泛的热障涂层是MCrAlY金属底层+Y2O3⋅ZrO2陶瓷面层的双层结构热障涂层(如图1所示)。金属底层材料MCrAlY与镍基单晶高温合金具有良好的界面结合能力,并且具有良好的抗氧化和抗热腐蚀能力,组元中M通常为Ni和Co,Ni具有优良的缓和热应力的能力,Co和Cr具有优异的抗氧化和抗腐蚀性能,Al可以和陶瓷面层中扩散来的氧反应形成致密Al2O3膜(thermally grown oxide,TGO)。此外,在底层中加入少量的Y可以起到氧化物钉扎和细化晶粒的作用,以提高热循环条件下Al₂O₃膜与基体的结合力,改善涂层的抗热震性能。陶瓷面层材料Y2O3⋅ZrO2应用最广泛的为6%~8%Y₂O₃稳定的ZrO₂,由于ZrO₂在1100℃左右的高温下会发生由单斜相(M相)向四方相(T相)的转变,伴随着3%~10%的体积变化,进而影响其抗热震性能,因此通常需要往ZrO₂中加入一些稳定组元,使相结构完全稳定或部分稳定,早期采用CaO,MgO作为稳定剂(国外对应材料牌号为Metco 201/202,CaO,MgO稳定ZrO₂涂层若长期或周期性地暴露于1039℃以上的环境中,CaO,MgO有扩散出稳定化ZrO2晶体之外的倾向,从而限制了该体系的最高使用温度,研究发现Y2O3作为稳定剂,且含量为6%~8%时,一般可获得不可转变的四方相。

图片  

图1 涡轮叶片热障涂层结构示意图

但当涡轮叶片在1200℃以上长期服役时,Y2O3⋅ZrO2陶瓷面层中仍会发生单斜相与四方相之间的转变,从而影响涂层的使用寿命。此外,MCrAlY金属底层仅可以在1100℃以下的环境中长期稳定地工作,当温度高于1150℃时,MCrAlY金属底层的氧化速率加快,TGO层变厚且内部产生的热应力和生长应力变大,最终导致氧化膜开裂和陶瓷面层剥落。因此,MCrAlY金属底层+Y2O3⋅ZrO2面层的最高服役温度在1150~1200℃左右,难以满足更高服役温度的需求。同时,随着航空发动机涡轮进口温度的提升,CMAS 在涡轮叶片表面的附着程度愈加严重,当Y2O3⋅ZrO2表面温度达到CMAS熔点时,致使CMAS熔化并在Y2O3⋅ZrO2表面润湿,然后通过毛细作用渗入到多孔的Y2O3⋅ZrO2涂层内部,在冷却的过程中CMAS凝固造成涂层杨氏模量增加、应变容限下降,并且CMAS的热膨胀系数比陶瓷面层和高温合金基体都小,在冷却过程中会产生很大的热应力,导致涂层中产生裂纹,在循环过程中裂纹逐步扩张,最终使涂层脱落。


为进一步提升热障涂层的使用温度(>1200℃)和抗CMAS腐蚀能力,国内外学者均开展了大量研究。关于金属底层材料,添加少量Re元素可大幅提升MCrAlY金属底层的高温抗氧化性能,Hf元素掺杂可以在TGO表面形成稳定的HfO2,使涂层的热循环寿命延长10倍,在MCrAlY金属底层表面沉积或者预制一层具有抗氧化性的薄层,比如在MCrAlY层上再喷涂上2~5μm的Al2O3阻挡层,可以进一步降低金属底层的氧化。除了对MCrAlY金属底层开展优化改进外,还发展了Pt-Al,NiAl合金、功能梯度材料等新型金属底层,其中Pt-Al底层已在单晶涡轮叶片热障涂层中应用。关于陶瓷面层材料,主要有对Y2O3⋅ZrO2进行掺杂改性和开发新型陶瓷面层材料两条研究思路,针对Y2O3⋅ZrO2的掺杂改性,单元素掺杂改变Y2O3⋅ZrO2热物理性能的好坏取决于掺杂元素的离子半径、含量以及掺杂元素能否形成高温稳定相,其中离子半径是最关键的因素,但是任何一种单元素的掺杂都很难获得比较理想的热物理性能,因此很多学者进行了多元素掺杂研究。郭洪波等针对一定浓度稀土元素氧化物共同掺杂的Y2O3⋅ZrO2性能进行研究,发现Gd2O3和Yb2O3共掺杂可明显提高相稳定性,1500℃热处理10h材料中M相增量仅有10%。其中,3%Gd2O3和3%Yb2O3共掺杂材料在25~1350℃范围内依然保持良好的相稳定,在25~1200℃范围内其热导率降至1.18~1.25W·m-1·K-1,明显低于Y2O3⋅ZrO2,热膨胀系数为(9.67~13)×10-6K-1,与Y2O3⋅ZrO2相当。针对新型陶瓷面层材料,开展了系列稀土锆酸盐热障涂层研究,LaZrCeO和Gd2Zr2O7是目前研究最广泛、也被认为最有前景的热障涂层材料,其在室温至熔点温度范围内,其热稳定性好,抗烧结性能好,氧扩散率低,此外新型热障涂层中的稀土元素可与CMAS附着物反应,迅速形成稳定的“牺牲层”,在涂层表面形成的致密连续膜层,阻碍反应进一步发生,并抑制CMAS向涂层内扩散,从而提高涂层抗CMAS腐蚀性能,保护涂层不产生剥落失效,稀土锆酸盐热障涂层已在国内外航空发动机单晶涡轮叶片上得到应用。

3

镍基单晶高温合金/热障涂层体系的考核评价

根据航空发动机涡轮叶片工作环境,对镍基单晶高温合金/热障涂层体系进行考核评价,是其工程应用前必须开展的工作。目前国内涡轮叶片热障涂层验收规范中对外观、厚度、气膜孔通畅性(缩孔率)、结合强度、相结构、显微组织、表面粗糙度、弯曲性能、热导率、抗热震性能、抗氧化性能以及对基体力学性能影响等提出了要求,并已经积累了大量数据,支撑了镍基单晶高温合金/热障涂层体系在航空发动机中的应用。但相关数据仅可用于对涂层质量进行评价,无法判断镍基单晶高温合金/热障涂层体系的使用寿命,以抗热震性能为例,其实验条件通常为将涂覆热障涂层的单晶涡轮叶片在马弗炉中加热至服役温度,保温5~10min后取出迅速放置常温水中(水冷)或通过压缩空气(空冷)进行冷却,此为一个循环,通过涂层不发生剥落的循环次数来判定其抗热震性能,在该过程中涡轮叶片的温度始终是均匀的,这与涡轮叶片在工作过程中的温度分布是不一致的,而涡轮叶片基体或热障涂层存在温度梯度可能会导致提前失效。


为更真实模拟涡轮叶片工作条件,目前国内外开始采用高温燃气火焰热冲击方式进行镍基单晶高温合金/热障涂层体系寿命考核。火焰热冲击是通过高温燃气对涂层表面加热,同时高温合金基材使用压缩空气进行冷却,进而在单晶基体及热障涂层内部形成温度梯度,这与涡轮叶片工作条件高度一致,燃气加热结束后采用压缩空气对涂层正面进行强制冷却,如此反复进行冷热循环,模拟发动机启停造成的热循环工况,可以说燃气火焰热冲击考核能够很好地模拟镍基单晶高温合金/热障涂层体系在发动机真实工作条件下的寿命。根据调研国外相关研究机构的先进模拟工况考核设备具备至少1500℃长期稳定考核能力,采用精准燃气流量控制技术实现温度控制的高稳定性。为充分评价镍基单晶高温合金/热障涂层体系模拟工况燃气火焰热冲击性能和使用寿命,必须建立相应的模拟工作环境工况考核平台,而考核平台建设的主要难点是超高温下焰流稳定控制和温度精准调控,国内也已开展了模拟服役环境的考核设备开发,具备开展高温考核评价的能力。


随着第三、四代航空发动机试车状态的增大和试车时间的延长,CMAS在涡轮叶片表面的附着程度愈加严重,因此抗CMAS腐蚀性能成为镍基单晶高温合金/热障涂层体系的重要评价标准之一。国内多家单位设计开展了抗CMAS腐蚀性能实验,主要思路为对使用后涡轮叶片表面黏附的CMAS进行成分分析,采用CaO,SiO₂,MgO,Al₂O₃等粉末材料,通过烧结、破碎、研磨等过程制备CMAS粉末,CMAS粉末添加至实验件表面的方式一般有两种,一是将CMAS粉末均匀涂抹在实验件表面,随后开展热循环实验,实验温度通常设置为1250℃,以保证CMAS粉末变为熔融状态,CMAS在热障涂层表面的润湿角测试即通过这种方式;另一种是将CMAS粉末添加在热冲击的火焰中,通过高温高速的火焰将CMAS喷涂在实验件表面,并进行热冲击循环实验。两种方式各有优劣,前者便于评估镍基单晶高温合金/热障涂层体系的抗CMAS腐蚀能力,后者更接近真实工作环境,便于获取镍基单晶高温合金/热障涂层体系的综合性能。


本团队采用上述方法,通过静态抗CMAS腐蚀性能实验、带CMAS燃气热冲击实验以及多因素耦合实验等对镍基单晶高温合金/热障涂层体系进行了考核评价,对比分析了Y2O3⋅ZrO2热障涂层和新型热障涂层的抗热冲击性能和抗CMAS腐蚀性能,其中DD5单晶/热障涂层体系静态抗CMAS腐蚀实验后宏观形貌如图2所示,可以看出新型热障涂层抗CMAS腐蚀性能优于Y2O3⋅ZrO2热障涂层。

图片  

图2 DD5单晶高温合金试片静态抗CMAS腐蚀后宏观形貌(a)Y2O3⋅ZrO2热障涂层2个循环;(b)新型热障涂层10个循环

4

镍基单晶高温合金/热障涂层体系的失效机理

研究表明,镍基单晶高温合金/热障涂层体系失效的主要原因是涡轮叶片工作过程中应力的产生,包括金属底层/陶瓷面层界面TGO生长应力、陶瓷面层相变应力、温度梯度分布引起的热应力、单晶基体与热障涂层热膨胀不匹配引起的热应力,其中最主要的是TGO生长应力。近些年,由CMAS侵蚀导致的镍基单晶高温合金/热障涂层体系失效越来越受到人们的重视。除此之外,单晶基体/金属底层界面的互扩散对基体性能的影响也应引起注意,在单晶合金和金属底层成分设计时加以考虑。


金属底层/陶瓷面层界面的TGO层在镍基单晶高温合金/热障涂层体系制备及使用过程中均可产生,其主要成分为α-Al₂O₃。一方面,TGO层由于其致密的结构可以保护金属底层免受进一步的氧化,另一方面,涡轮叶片工作过程中TGO的过度和不均匀生长会在界面处产生巨大的应力,在高温环境下,随着金属底层中Al元素的消耗,导致Ni(Al,Cr)₂O₄尖晶石的形成,而这种快速生长的氧化物将会促进TGO层的增厚,同时使垂直于金属底层/陶瓷面层界面的外平面生长应力也随着增大,从而引起界面的开裂和陶瓷面层内裂纹的生成,这将引起失效。


镍基单晶高温合金与金属底层在化学成分和相结构上存在显著差异,金属底层制备及应用过程中与单晶合金基体间会不可避免地出现互扩散行为。在高温保温后,在单晶基体/金属底层界面通常会形成元素互扩散区(IDZ)和二次扩散区(SDZ),其中IDZ是以β-NiAl相为基体,同时含有大量的TCP相,而SDZ仍然保留了合金基体中的大部分γ'相以及单晶取向,但由于合金中Mo,Re,W等难熔元素溶解度的下降会导致大尺寸的针状TCP相析出。IDZ和SRZ的形成都是单晶合金表面层的γ/γ'两相组织被破环的同时形成新相,导致单晶基体/金属底层界面组织失稳进而诱发金属底层失效。单晶基体/金属底层界面在热力耦合条件下,界面元素互扩散会导致单晶基体表面再结晶、界面相变、单晶γ'相反常筱化等现象,造成单晶合金基体力学性能下降。


CMAS对镍基单晶高温合金/热障涂层体系的影响取决于涡轮叶片的工作温度,低温时(<735℃)CMAS颗粒撞击表面涂层,会引起腐蚀磨损、气膜孔堵塞和涂层局部剥落,在1200℃左右,CMAS会在涂层表面变为熔融态侵蚀涂层。高温下CMAS侵蚀热障涂层的失效机制包括热化学与热力学两方面,热化学方面CMAS会与陶瓷面层发生反应,导致陶瓷面层分解并丧失性能;在热力学方面,腐蚀产物的生成以及ZrO₂的相变会使涂层内部出现大量的体积膨胀,造成涂层内部产生拉应力,多次热循环后,导致涂层内部产生横向裂纹,造成分层剥落。本团队对航空发动机DD5单晶导向叶片热障涂层剥落开展失效分析时,发现存在CMAS侵蚀热障涂层形貌,如图3所示,表明CMAS侵蚀是导致该叶片热障涂层提前失效的原因之一。

图片  

图3 DD5单晶导向叶片热障涂层剥落区域面扫描结

5

结束语

随着先进航空发动机涡轮进口的不断提高,涡轮叶片镍基单晶高温合金/热障涂层体系的应用研究必然受到国内外研究人员的持续关注,结合相关调研以及本团队在工程应用中遇到的问题,认为镍基单晶高温合金/热障涂层体系未来的应用研究重点应集中在以下几个方面:


(1) 目前国内在镍基单晶高温合金成分设计、组织性能评价等方面已具备丰富经验,单晶合金/金属底层界面的互扩散行为及其对单晶基体力学性能的影响的应用基础研究应进一步重视,从单晶合金与金属底层成分协同设计等方面开展研究,提升单晶合金/金属底层界面匹配性。


(2) 稀土改性ZrO₂和Gd₂Zr₂O₇,LaZrCeO等新型热障涂层材料可以提高陶瓷面层的耐温能力和隔热效果,已开展涡轮叶片热障涂层试制及应用,但其工作可靠性仍需进一步提升,在新型陶瓷面层材料开发时,应综合考虑材料的力学性能、热物理性能及其微观结构,使三者达到最优匹配。


(3) 金属底层起到了重要的抗高温氧化和抗脱落作用,但是相比于镍基单晶高温合金和陶瓷面层而言,金属底层材料的研制和应用研究力度不足,功能梯度底层具有广阔的应用前景,需开展制备工艺优化研究;另外,在涂覆陶瓷面层前对金属底层进行后处理也可以提高其综合性能,应开展系统研究。


(4) 在高代次单晶合金及新型热障涂层材料研制阶段,对镍基单晶高温合金/热障涂层体系进行考核评价,可以降低整机验证风险,考核评价的实验条件应尽可能的接近涡轮叶片工作环境,应研究构建可同时施加机械应力、高温燃气和CMAS颗粒的多因素耦合实验装置,并基于实时监测技术,综合掌握镍基单晶高温合金/热障涂层体系状态衰减及表面状态变化。

声明: 本文来源于网络, 仅供交流分享, 若涉及版权等问题请留言, 我们会及时处理


来源:两机动力先行
疲劳化学航空核能铸造裂纹材料控制
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2025-11-09
最近编辑:4小时前
两机动力先行
其它 聚焦航空发动机/燃气轮机关键技术...
获赞 99粉丝 121文章 421课程 0
点赞
收藏
作者推荐

涡轮丨上交大:增材制造燃气轮机涡轮叶片的扰流柱三维形貌对流动换热影响机制研究

扰流柱单元是涡轮叶片内部冷却中常见的增强换热的手段。随着先进加工方法的成熟发展,部分涡轮叶片逐步开始采用增材制造的方式进行加工。为了研究增材制造表面特征对扰流柱流动换热的影响,本文开展计算机断层扫描(CT)实验,测得铸造和增材制造加工的涡轮叶片扰流柱单元的表面特征,表面特征显示增材制造相较于铸造平均增加了85%的粗糙度,且圆形度较低。采用数值模拟的方式依次探索了包括光滑扰流柱在内的三种扰流柱在不同雷诺数下的流动换热特性。仿真结果表明:仅从增加粗糙度的角度考虑,加工方式引入的表面粗糙度变化能够增强换热,相对于光滑壁面设置,增材制造的提升换热效果更加明显,最大分别可以提升21%和10%的上、下端壁换热,同时增材制造工艺也会导致流动损失增加;进一步考虑加工方式引入的表面形貌的变化,增材制造的扰流柱表面起伏较大,引起流动在扰流柱表面过早分离,对换热起到了削弱的作用。关键词:燃气轮机;涡轮叶片冷却;扰流柱;强化换热;增材制造;CT扫描1引言为了满足工程应用中对燃气轮机更高效率的追求,涡轮前入口温度在不断提高,目前已经远远超过金属所能承受的极限温度。为了避免涡轮叶片热烧蚀,冷却技术的引入是极其关键的。常用的冷却技术包括外部冷却和内部冷却,其中内部冷却通过把冷却空气通入叶片内部的专门通道进行冷却,其冷却性能依赖于内部通道的几何结构设计。扰流柱排列是增强内部冷却较为常见的结构布置。在实际应用过程中,叶片加工所引入的表面结构特征对扰流柱的换热性能影响往往是不能忽视的,因此厘清加工制造引入的表面特征在涡轮叶片内部扰流柱冷却过程中流动换热的影响一直是相关领域人员的关注方向。作为典型的涡轮叶片内部换热单元,扰流柱冷却通过增大换热面积,增强流动掺混来增强内部表面的对流换热。国内外学者开展了大量关于扰流柱表面特征对其流动换热影响的相关研究。Pham等通过k-ω湍流模型仿真分析了C型端壁对于圆柱扰流柱换热方面的影响。Luo等提出了一种弯曲扰流柱的设计,并通过realizable k-ε数值模拟方法验证了其优于传统的右上扰流柱的换热性能。Zhou等采用SST k-ω模型研究了不同形状的扰流柱横截面(方形、圆形、椭圆形、NACA叶型轮廓和水滴形)对扰流柱换热的影响,发现椭圆形和NACA叶型轮廓的扰流柱换热最好。Yeranee等同样研究了不同形状的扰流柱(圆柱、方菱形、三角形和优化型)的流动换热性能,数值模拟发现优化型结果从增强流动掺混和增加换热面积两个方面起到了对换热的增强作用。除了几何形状,一部分学者研究了扰流柱表面粗糙度对流动换热的影响。Buresti通过开放式喷气风洞实验,研究表面粗糙度对圆柱绕流中亚临界、临界、超临界和过临界状态转变的影响。雷诺数在2.6×104~2.8×105,相对粗糙度在1.0×10−3~12×10−3,实验证实了粗糙度对扰流流态的显著影响,特别是粗糙度可以缩短过渡区域,在高粗糙度的圆柱中甚至可能消失。李广超等模拟并分析了在不同粗糙度下,雷诺数在Re=1×104~5×104的局部和平均努塞尔数分布,发现粗糙度的增加使扰流柱表面平均换热系数增大;雷诺数越大时,粗糙度对换热系数的影响越明显。根据以上的研究可以发现,扰流柱及其端壁的表面形貌特征对于其本身的流动换热特性影响是显著的。目前,随着金属增材制造技术的成熟发展,新的加工方式引入的特殊表面结构对加工件应用过程中流热方面的影响亟需开展相关的评估研究。Kaur等通过X射线计算机断层扫描(CT)扫描和光学轮廓仪对增材制造部件进行了表面形貌获取,并以此进行数值模拟建模,研究了其对十四面体晶格的流动换热的影响。增材制造方法在涡轮叶片冷却领域同样开始普及,Snyder等采用CT扫描技术获得了不同加工角度的气膜孔样,扫描的结果表面气膜孔内部流道的加工偏差受到加工角度明显的影响。Furgeson等同样采用CT扫描技术获得了某种增材制造型孔的内部流道特征,发现气膜孔的实际直径要大于设计值,而气膜孔壁面厚度则略小于设计值。可以发现,增材制造对叶片气膜冷却孔的影响是比较显著的,但是目前有关涡轮叶片扰流柱的增材制造表面形貌表征以及相关的流动换热影响研究较少。因此,本文研究的出发点正是基于CT扫描技术分布,获取增材叶片和传统铸造叶片内部扰流柱几何结构与表面特征,通过直接对加工的扰流柱进行计算域的建模,采用数值模拟的方式探讨两种加工方式的表面结构对扰流柱单元流动换热的影响。研究结果首先将与典型圆柱绕流努塞尔数实验关联式进行验证,后续从努塞尔数、流动损失等方面进行讨论分析。2方法2.1 叶片加工与扰流柱形貌提取本文首先通过金属铸造和增材制造分别获得了两种加工方式下的某实验型号涡轮叶片。为了获取叶片内部冷却通道内的扰流柱及其两端壁面的表面特征,常规的光学形貌测量方法是无法实现的。参考Kaur等通过X射线计算机断层扫描(CT)法对燃气轮机叶片尾缘针翅结构表面形貌扫描的方法,本文采用相同的方法对本实验获得的真实叶片进行扫描,具体测量方法如图1(a)所示,通过把测试叶片的尾缘部分放置在旋转载物台上环形扫描即可获得其内部的形貌特征。对叶片扫描后,实际加工获得的叶片模型如图1(b)和(c)所示。 图1 计算机断层扫描(CT)法获得的叶片模型本文的研究目标是某实验型号涡轮叶片内部的扰流柱单元,采用CT扫描的方法分别扫描了叶片内多个扰流柱的表面形貌,获取了两种加工制造方式下的扰流柱单元模型各个位置的平均表面粗糙度。随后测量了同一排6个扰流柱的直径,选取了其中具有代表性的扰流柱单元开展后续的研究,该扰流柱的表面形貌特征与未选择的扰流柱是相似的,且直径最接近6个扰流柱的平均直径,如图2所示。具体的扰流柱形貌如图3所示,选取了从扰流柱根部的向上的7个截面的数据进行两种加工方式下典型扰流柱形貌的分析,其中各处的截面的水利直径dₙ与圆形度φ的计算公式如下: 式中A为各个截面的面积,C为各个截面的周长,光滑扰流柱在叶片内的设计直径为1mm。通常,圆形度φ的取值在0~1之间,等于1则表明是标准的圆,越接近1表示该截面接近圆形,远小于1时表示该截面形状偏离圆形,可能是不规则或者是椭圆等。 图2 扰流柱单元提取图3表明铸造扰流柱的直径变化较大,扰流柱与端壁连接处呈现扩张的喇叭状,但是圆形度比较接近1;增材制造扰流柱的直径变化较小,但是存在很多不规则的截面,这是因为在金属粉末固化过程中温度控制不均匀,扰流柱表面存在更多的加工缺陷,也表现得更加粗糙。表1总结了7个截面上扰流柱直径与圆形度的统计数据,其结果与图3观察到的结果是一致的,并且增材制造的扰流柱平均直径要大于铸造扰流柱。此外,根据CT扫描获得的上下端壁及扰流柱形貌,计算了三个表面上的算术平均粗糙度Ra,具体结果如表2所示。可以发现铸造所引入的表面粗糙度要小于金属增材制造所引入的表面粗糙度,其中金属增材制造的上壁面粗糙度要远大于铸造所引入的粗糙度。相较于铸造,增材制造引入的粗糙度在三个面上平均增加了85%的粗糙度。后续的研究正是基于实际加工获得的扰流柱形貌进行计算域的简化建模,同时将光滑扰流柱作为对比基准。 图3 两种加工方式扰流柱形貌对比 表1 不同加工方式下扰流柱形貌统计量 表2 不同加工方式下的扰流柱及上下端壁算术平均粗糙度Ra(单位:μm)2.2 数值模拟验证为了匹配雷诺数以及后续圆柱绕流实验关联式的验证,本文将光滑扰流柱放大8倍(D=8mm)进行计算。由于对于单个扰流柱单元的附近端壁曲率变化不大,计算域模型的上下端壁简化为平直壁面,夹角与扰流柱所处位置的上下壁面夹角一致,经过测量为12.6°。整个计算域几何参数如图4所示,除了光滑壁面的扰流柱外,还建立了铸造和金属增材制造获得的扰流柱计算域模型,铸造和增材制造扰流柱根据表1中的平均直径放大到8mm。 图4 扰流柱单元计算域根据以前学者对扰流柱研究采用的湍流模型,本文采用SSTk-ω湍流模型来计算扰流柱附近的流动与换热。在进行本文的数值模拟研究之前,首先需要对湍流模型进行验证。模型验证的计算域如下图5所示,速度入口按照不同的雷诺数(Re=2917、5835、9725、14587)分别设置,雷诺数定义见式(3),来流温度为318.15K;出口静压为0.1MPa,温度为318.15K;上下壁面与圆柱壁面设置为恒温壁面,温度为298.15K;前后两边则为周期性边界条件。圆柱中截面的热流密度通过数值模拟获得,温差已知,因此数值模拟的圆柱中间截面上的平均努塞尔数可由式(4)获得。 从表3可以看出,随着雷诺数的增大,努塞尔数的增强表面圆柱表面的换热是增强的。同时雷诺数越大,仿真获得的圆柱绕流平均努塞尔数与实验关联式所计算的相对误差越小。当雷诺数大于9000左右时,预测误差在3%以内。同时本文对比了在高雷诺数条件下圆柱表面随角度θ变化的局部努塞尔数分布,从图6可知,仿真结果与实验结果的量值与趋势是一致的,预测的回升点均在80°左右,因此湍流模型的选择是合理的。 图5 湍流模型验证计算域模型 图6 局部努塞尔数分布(Re=70800) 表3 不同雷诺数卜仿真与实验关联式获得的圆柱周向平均努塞尔数对比对于带有增材制造特征的扰流柱单元计算域,采用Fluent meshing进行多面体网格的划分,网格的细节如图7(a)所示,在扰流柱的壁面上对网格进行了加密,同时在所有壁面上附加了棱柱层网格,第一层网格厚度为了保证复杂形状扰流柱表面的计算精度,保证y+在1左右。为了验证网格无关性,本文计算了平均努塞尔数随网格数目的变化和扰流柱前缘驻点位置沿着z方向上压力随网格数目的变化曲线,分别如图 7(b)和(c)所示。从结果可以看出,随着网格数目的增加,换热和压强结果随着网格数目改变而变化的幅度在逐渐减小,当网格数目大于15万时,计算结果可以认为已经趋于稳定,考虑到计算效率与计算准确度,后续的光滑扰流柱、铸造扰流柱和金属增材制造扰流柱均按照15.5万网格的生成方式进行构建。 图7 计算域多面体网格划分与网格独立性验证3结果与讨论在验证完湍流模型的合理性和网格数目的无关性后,将与验证湍流模型时采用的边界条件设置赋予后续的扰流柱单元计算域。其中扰流柱本身保留CT扫描获得的表面形貌特征。上下端壁则根据FLUENT内部常用的沙砾模型赋予粗糙度,其中等效沙砾直径Ks按照式(6)的变换关系分别计算并设置,粗糙度均匀度Cs默认为0.5。为了验证此类壁面粗糙度添加方法的可靠性,根据粗糙壁面圆柱绕流的换热实验结果,首先在图5的计算域中的圆柱表面赋予了Ks/D=900×10−5量级的粗糙度,随后进行了两个不同来流雷诺数下的验证,将仿真获得的平均努塞尔数与实验经验关联式结果进行对比,其中光滑壁面采用式(5)进行计算,粗糙壁面采用式(7)进行计算,表4中结果显示仿真计算的换热努塞尔数与实验结果的相对误差均在5%以内,并且在壁面添加均匀粗糙度能够增强壁面的换热。 表4 圆柱周向平均努塞尔数对比验证后续正式计算的工况除了加工方式不同的扰流柱计算域以外,本文还研究了不同来流雷诺数(Re=8536、17071、25606)下对扰流柱流动换热的影响,雷诺数范围涵盖了叶片内部扰流柱内的大多数来流情况,本节后续将讨论数值模拟结果。3.1 不同加工方式下的扰流柱流场分析为了进一步分析增材制造加工的扰流柱对内部流动的影响,本文提取了z/D=0.5平面上的流场进行分析。图8显示了Re=17071时z/D=0.5平面上的流线以及标准化速度云图,其中光滑扰流柱和铸造扰流柱在该平面上的流场是比较对称的,而金属增材制造获得的扰流柱因其表面突起的不对称性,从扰流柱右侧经过的流体出现分离壁面的现象要早于光滑扰流柱,流体分离后的扰流柱壁面努塞尔数将因此大幅减小。图9绘制了Re=17071下的z/D=0.5平面上的湍动能云图,其中湍动能使用入口速度的平方进行了归一化处理。可以看到金属增材制造扰流柱下游的湍动能水平最高,光滑扰流柱的最低。湍动能越大,说明该区域的湍流水平越高,能量掺混传递的越快,越有利于该区域对应端壁的换热。此外,由于增材制造引起的扰流柱不对称性在湍动能分布上的表现更加明显,铸造和金属增材制造扰流柱均表现出了较为明显的湍动能不对称分布。最后,图9对比了不同雷诺数下各种加工方式的扰流柱的阻力因子大小。阻力因子的定义如下式: 式中𝑝in和𝑝out分别为流体域入口和出口的压力,𝜌in为入口处来流的密度,𝐷H为入口水力直径,𝐿为整个计算域长度。对于所有工况,随着雷诺数的增大,阻力因子不断增大,其中铸造的扰流柱在低雷诺数下阻力因子要略小于光滑扰流柱,这主要是由于铸造扰流柱表面在靠近中部区域的直径较小,使得其对流道的阻塞面积小导致的;铸造扰流柱与端壁之间为渐扩连接,衔接更加光顺,能够减小马蹄涡等二次流的流动损失。而在高雷诺数下,由于壁面粗糙度导致摩擦引起的损失占主导,因而铸造扰流柱的阻力因子大于光滑扰流柱的。对于金属增材制造加工获得的扰流柱,由于其表面粗糙度远大于另外两种构型,其阻力因子在不同工况下都是最大的,并且在高雷诺数下损失系数要明显高于光滑和铸造扰流柱扰流柱。其中,在本文雷诺数研究范围内流经金属增材制造获得的扰流柱时阻力因子最大,最大流动损失增加了1.6%。 图8 三种扰流柱的绕流流场(Re=17071,z/D=0.5) 图9 三种扰流柱下游的湍动能分布(Re=17071,z/D=0.5) 图10 扰流柱单元的阻力因子𝑓随雷诺数的变化 3.2 不同加工方式和雷诺数对扰流柱表面努塞尔数的影响 图11 不同加工方式和雷诺数下的扰流柱表面努塞尔数分布 图12 扰流柱的面平均努塞尔数随雷诺数的变化3.3 不同加工方式和雷诺数对上下端壁表面努塞尔数的影响 属增材制造粗糙度下的下端壁换热仅仅略强于铸造粗糙度的,并且二者也是要大于光滑下端壁的换热。粗糙度对换热的增强随着雷诺数的增加也愈发明显。需要注意的是,不同于上节中扰流柱表面的换热结果,由于端壁并未引入真正的形貌而是仅仅设置了粗糙度数值变化,所以相较于光滑端壁,铸造和金属增材制造的上下端壁的换热均得到了增强。 图13 不同加工方式和雷诺数下的上下端壁表面努塞尔数分布 图14 上下端壁的面平均努塞尔数随雷诺数的变化4结论本文通过研究,得到如下结论:(1)相较于铸造工艺,增材制造加工方式引入的粗糙度在三个面上平均增加了85%的粗糙度,其扰流柱的平均圆形度为0.832,铸造的为0.995,同时增材制造扰流柱沿径向的直径变化要大于铸造扰流柱;(2)扰流柱能够增强下游的湍动能,从而增强掺混,提高扰流柱端壁的换热,并且端壁的粗糙度越大,换热越强。其中增材制造方式相较于光滑的扰流柱单元壁面能够提高上端壁7%~21%的换热,提高下端壁1.3%~10%的换热;铸造的上下端壁提升量分别为2.5%~12%和0.2%~9.8%。(3)除了能够增强换热,增材制造引入的表面形貌以及端壁粗糙度会增大流动损失,其中金属增材制造的扰流柱最大增加1.6%的流动损失。(4)本研究范围内,考虑加工引入的表面特征,相较于光滑扰流柱,铸造工艺下引入的表面形貌能够略微提高扰流柱表面的努塞尔数;而增材制造工艺下引入的表面形貌会导致扰流柱表面流动死区较多,绕流流动过早分离,进一步地将轻微削弱扰流柱本身的换热。声明: 本文来源于网络, 仅供交流分享, 若涉及版权等问题请留言, 我们会及时处理 来源:两机动力先行

未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习计划 福利任务
下载APP
联系我们
帮助与反馈