涡轮叶栅气热耦合数值方法和验证研究.pdfVIP

中国工程热物理学会 热机气动热力学 学术会议论文 编号： 072032 涡轮叶栅气热耦合数值方法与验证研究 颜培刚 董平周驰 冯国泰 哈尔滨工业大学推进理论与技术研究所 (Tel：045l peigang—y@sina．com．cn) 摘要：应用三阶精度TvD格式和自由型曲面网格技术，对涡轮内冷叶片进行气热耦合数值模拟。结合 商用软件采用多种湍流模型计算方法的结果进行实验验证和对比分析，结果表明，商用软件采用传统 湍流模型对温度边界层的模拟是不合理的，采用考虑转捩的湍流模型，调整转捩起始雷诺数，可以较 好的模拟温度边界层的热传导。本文对非转捩区边界层温度的计算结果与实验结果吻合较好。 关键词：涡轮：气热耦合；热传导；数值验证 l、引言： 在现代高性能燃气涡轮的设计中，随着高压涡轮的气动负荷与热负荷的进一步增加， 撞击、对流、气膜结合的高效气冷设计与气热多场耦合研究已经成为了及其关键的研究 内容。特别是对于当今高负荷高温涡轮，其流动、热传导与热应力是互相制约的，单物 理场的计算误差会造成整体设计的严重偏差，因此必须进行准确而有效的耦合研究…。 国外学者进行了很多气热耦合的定常与非定常数值计算，但是计算涡轮叶片表面以及叶 片固体内部温度方面，对充分发展的湍流壁面流动区域预测较好，而对射流区和转捩区， 传统的湍流模型都不够准确心‘5J。文献[6]对典型的Mark工I叶片进行了验证，在参数光 滑区预测结果与实验吻合较好，在吸力面激波后区域，误差仍旧偏大。利用边界元方法 进行计算的CHT方法H1也受限于较低的湍流模拟精度。耦合换热研究目前仍然是学术界 公认的有待解决的一个问题。 本文的研究作为高压高负荷涡轮的气热弹耦合设计体系的一部分，通过生成较高质 量的多块网格结构，以及利用高精度的数值计算格式，发展计算程序并进行实验验证， 以获得热端部件数值研究的可靠性。结果表明，在计算壁面流体附面层内热传导问题时， 湍流模型对计算结果的影响是及其显著的。对商用软件以及国外公开文献发表的计算综 合分析结果表明，采尉考虑转捩影响的湍流模型对MarkII叶片的计算可以获得最接近 实验的温度和热传导系数值，对C3X叶片较难得到满意的预测结果。本文采用了低雷诺 数q—w湍流模型与高精度的计算格式，获得了与CFX转捩模型最接近的计算结果，在非 转捩区具有较高的可信度。准确预测气热耦合在研究热端问题方面，特别是对近来迅速 发展的高压、高温、高负荷叶片的研究是必不可少的。 2、数值方法： 2．1控制方程及离散 流场基本控制方程采用可压缩雷诺平均Ns方程，离散方法采用具有TVD性质的三阶 精度格式和隐式时间推进格式佃’。湍流模型采用q—w模型旧1。 a71 温度场控制方程的形式为：÷=V．(rV丁)，其中，=旯／心，详细见文献[10]。 刁f 2．2耦合方法： 42 热传导边界条件：Ar(a丁／a船)，=办(L—c) 采用预处理技术n¨，将温度作为守恒量来处理。采用内部边界条件直接耦合方法¨引： ， (旯。／见，)·(△胛，／△玎，)疋+弓 ” 1+(旯，／A，)·(△，?r／△刀，) 其中w表示叶片表面，s、f分别表示固体和流体，n表示到壁面的法向距离。在计 算过程中，可以在保证计算精度的同时，加快耦合速度。本文使用绝热壁面条件下的流 场作为初始值进行耦合计算，流场和叶片温度场只需耦合10次左右即可收敛。 2．3计算网格： 结合AutoGrid网格生成器与自由型曲面网格生成技术”引，生成具有较高网格质量的 多块结构化网格结构。流场区域采用H一0一H型网格拓扑结构，叶片区域采用外部0型网 格，内部H型网格的拓扑结构。 图1计算网格 为精确模拟固体场与流场交接面的温度梯度，叶片表面第一层网格的无量纲y+距