边界元法在燃气涡轮叶栅气热耦合计算与应用.pdfVIP

中国工程热物理学会 热机气动热力学 学术会议论文 编号：082108 边界元法在燃气涡轮叶栅气热耦合计算 中的应用1 王振峰，颜培刚，黄洪雁，韩万金 （哈尔滨工业大学能源科学与工程学院 458 信箱，哈尔滨 150001） (Tel:0451 E-mail：power6411715@126.com） 摘要：本文将边界元法引入燃气涡轮的气热耦合计算中，耦合有限差分法和边界元法对某气冷燃气涡 轮叶栅热环境进行耦合分析。流体部分计算代码 (HIT-NS Code) 采用有限差分方法求解 Navier-Stokes 方程。固体部分计算代码 (3D-BEM Code) 应用边界元法求解固体区域热传导方程。文中应用解析解 和商业软件 Fluent 对边界元程序进行准确性验证，边界元法的计算结果与解析解以及商业软件模拟结 果非常的一致。不同于其它计算方法，边界元法避免了计算域的网格生成工作，节省了数值模拟的前 处理和求解时间。另外，边界元法引用了问题的基本解，具有解析与离散相结合的特点，计算精度高。 应用边界元法耦合计算的结果与实验数据较为接近，说明采用边界元法的气热耦合程序计算结果更接 近于叶栅内真实的热环境。另外，本文应用商业软件 CFX 11.0 对实验叶栅进行三维数值模拟，不同的 湍流模型所得的计算结果不同，讨论了不同湍流模型对涡轮叶栅热环境准确性计算的影响。 关键词：气热耦合，边界元法，有限差分法，燃气涡轮，湍流模型 0 前言 当前，随着对涡轮发动机性能越来越高的要求，使得涡轮叶栅内热环境分析变得非 常重要。因此，在涡轮发动机设计过程中，除了要考虑气动效率外，涡轮叶栅内部高温 部件的温度分布成为另一个重要的因素。近几年，气热耦合[1,2]计算方法成为透平机械领 域中热环境分析的有效工具，耦合计算中热边界条件由热量交换过程动态地加以决定而 不是如传统计算方法一样预先规定热边界条件，因此耦合分析更加真实可信。文献[3]借 助数值模拟技术分析了三维跨音速涡轮导向叶栅在绝热、等温、耦合换热三种不同壁面 边界条件下的温度场分布，结果表明耦合换热计算的精确性。文献[4]也指出采用和不采 用耦合计算，得出的叶片的温度场分布并不一样，采用耦合换热计算时叶片温度场更为 真实可靠。准确的温度分布可以指导设计人员对高温部件进行有效的冷却，提高涡轮发 动机性能。一些研究人员已经开始应用耦合热分析方法研究涡轮叶栅的气热耦合问题 [5-7] 。大多数研究方法基于有限差分、有限元或有限体积法，这些计算方法需要同时计 算固体和流体区域的网格单元，在现代的涡轮设计领域中为了提高气动性能涡轮叶片通 常具有弯、掠等几何特征，叶片特殊的几何轮廓对其网格划分工作带来很大的挑战，网 格质量的好坏对耦合计算精度有很大的影响。然而，Li和Kassab[8-11]发展了一种不同的 耦合热分析方法。这种耦合热分析方法采用边界元计算方法求解固体区域拉普拉斯方 资助项目：中国自然科学基金资助项目(No 程。区别于有限元和有限体积等计算方法，边界元法不需要计算固体内部网格单元，而 仅仅计算固体表面的网格，这部分网格可由流体侧网格提供，因此在耦合换热计算过程 中避免了叶片等复杂几何体的网格生成工作，节约了数值模拟前处理和计算的时间。同 时，边界元法引用了问题的基本解，具有解析与离散相结合的特点，计算精度高。本文 对三维边界元程序进行了详细验证，并且应用边界元法和有限差分法对某气冷涡轮进行 了气热耦合计算。着重分析了采用边界元法进行耦合热分析计算的优势，并将耦合计算 结果与实验结果比较，结合商业软件CFX 的不同湍流模型数值模拟结果，讨论了湍流模 型对计算结果的影响。 1 控制方程 1.1 流体区域 流体区域计算程序基于可压 Navier-Stokes 方程进行编写，采用有限差分方法求解， 差分格式为具有 TVD 性质的三阶精度 Godunov 格式，湍流模型为 B-L 模型，湍流普朗 特数为 0.9，层流普朗特数为 0.