燃气轮机的 “降温神器”,原来这么设计-夜雨聆风

燃气轮机的 “降温神器”,原来这么设计

导读：燃气轮机作为能源发电、航空航天等领域的核心设备，工作时内部温度常突破 1400℃，极端的高温环境会直接导致叶片材料老化、蠕变甚至断裂，严重影响设备的运行安全和效率。为了解决这一问题，高效的冷却技术成为燃气轮机设计的核心环节，而内部冷却通道则是叶片冷却的关键结构。本次研究借助计算流体动力学（CFD） 技术，重点分析了不同冷却空气流速和湍流模型对燃气轮机叶片冷却效果的影响，最终为实际工程中优化叶片冷却设计、选择合适的仿真模型提供了实用参考，既能提升冷却效率，又能兼顾计算成本和设备运行经济性。

<一、研究背景：为何燃气轮机叶片冷却如此重要？>

燃气轮机的热效率和工作温度密切相关，温度越高，燃料的能量利用率就越高，但这也对叶片材料提出了严苛要求。即便采用镍基高温合金这类耐高温材料，若无有效冷却，叶片仍会因热应力过大出现损坏。

目前内部冷却通道是燃气轮机叶片的主流冷却方式，通过让冷却空气在叶片内部流动，带走热量降低叶片温度。但冷却通道内的湍流流动和热量传递相互作用，过程十分复杂，想要精准预测冷却效果，湍流模型的选择和冷却风速的控制成为两大关键难题。

以往的研究要么只关注湍流模型，要么固定冷却风速分析，而本次研究首次将三种主流湍流模型与三种不同风速结合，全面分析其对叶片温度分布、换热效率和压力损失的影响，填补了相关研究的空白，也让结论更具工程参考价值。

<二、研究方法：从模型搭建到仿真分析，步步严谨>

本次研究以燃气轮机第一级转子叶片为研究对象，通过数值建模、网格划分、设定边界条件，结合三种湍流模型开展 CFD 仿真，全程保证分析的准确性和可靠性。

1. 叶片几何模型与网格划分

研究的转子叶片为空心结构，包含凹面压力侧、凸面吸力侧和支撑平台，内部设计了三条独立的冷却流道，能让加压冷却空气充分流经叶片内部。叶片根部因热影响小、主要起机械连接作用，未纳入仿真分析。

为了让仿真结果精准，研究对叶片进行了高质量网格划分，最终采用的网格包含 1447157 个节点和 7385005 个单元。网格的最大偏斜率为 0.24178（属于优秀等级），最小正交质量为 0.76005（属于非常好等级），能有效捕捉冷却通道内的复杂流动和温度梯度。

同时研究开展了网格无关性验证，分别测试了 350 万、740 万和 1520 万单元的网格，发现 740 万单元后，继续增加网格数量，叶片平均温度仅轻微变化，而计算时间却大幅增加。因此740 万单元的网格是精度和计算效率的最优选择。

图 1 燃气轮机转子叶片的几何结构：(a) 燃气轮机叶片外部视图；(b) 转子叶片整体结构；(c) 叶片内部冷却通道；(d) 冷却气流路径横截面；(e) 燃气轮机叶片翼型形状

2. 核心控制方程

仿真基于流体动力学和传热学的基本原理，采用质量守恒、动量守恒和能量守恒方程，同时忽略瞬态效应，按稳态流场分析（聚焦长期冷却性能，瞬态波动对整体冷却效率影响极小）。

为了描述冷却通道内的湍流流动，分别引入了三种湍流模型的专属传输方程，精准捕捉湍流动能、耗散率等关键参数，这也是不同湍流模型仿真结果存在差异的核心原因。

3. 湍流模型与边界条件

本次研究选用工程中最常用的三种湍流模型，各有其适用特点：

k-ω 剪切应力传输模型（k-ω SST）

擅长捕捉壁面附近的湍流和逆压梯度流动，结合了壁面区 k-ω 模型和远场区 k-ε 模型的优势；
k-ε 模型

结构简单、计算效率高，适用于完全湍流流动，是工业仿真的常用选择；
雷诺应力模型（RSM）

高阶湍流模型，能直接求解雷诺应力，精准捕捉各向异性的湍流效应，适合复杂流动场景。

仿真的边界条件严格贴合燃气轮机实际工作状态，关键参数设定如下：

冷却空气入口风速：30m/s、60m/s、90m/s（分别代表低、中、高风速），入口温度 20℃；
燃气主流温度 1150℃，主流换热系数 400W/(m²・℃)；
叶片材料为镍基高温合金 CMSX-4，热导率 23W/(m・℃)，具备优异的高温强度和抗腐蚀性；
壁面采用无滑移条件，非冷却表面为绝热面，入口湍流强度 5%。

仿真采用 SIMPLEC 算法求解，当连续性方程残差低于 10⁻⁶、能量方程残差低于 10⁻⁹时，判定为计算收敛，保证结果的可靠性。

图 2 燃气轮机转子叶片的网格结构注：网格单位为毫米，核心区域做精细化处理，壁面设置膨胀层提升近壁面仿真精度

图 3 基于燃气轮机叶片平均温度的网格无关性研究注：横坐标为网格单元数量，纵坐标为叶片平均温度（℃），网格单元达到 740 万后，温度结果趋于稳定

<三、研究结果：风速和湍流模型，如何影响冷却效果？>

本次研究从温度分布、热通量、湍流动能、流速流场、压力分布五个维度，分析了不同风速和湍流模型的作用效果，核心结论清晰展现了冷却性能的变化规律，同时也揭示了各模型的优劣。

1. 叶片温度分布：风速越高冷却越好，k-ω SST 模型精度最优

冷却风速是影响叶片温度的关键因素，风速越高，叶片降温效果越显著：30m/s 低风速下，叶片整体温度偏高，前缘和中弦区域高温尤为明显；60m/s 中风速时，冷却效果大幅提升，叶片温度均匀性显著改善；90m/s 高风速时，叶片温度降至最低，相比 30m/s 工况，最高温度降幅达 7%，核心高温区得到有效控制。

而不同湍流模型的温度预测效果差异明显，k-ω SST 模型表现最优：在所有风速下，该模型都能实现最低的叶片最高温度和最均匀的温度分布，温度预测偏差低于 3%；雷诺应力模型（RSM）次之，能精准捕捉湍流的各向异性效应，但低温区略少于 k-ω SST 模型；k-ε 模型表现最差，各风速下均出现更高的局部高温，对换热效率的预测存在偏差，不适合高精度冷却仿真。

图 4 燃气轮机叶片的温度（℃）分布注：(a) 30m/s 风速下 k-ω SST 模型；(b) 30m/s 风速下 k-ε 模型；(c) 30m/s 风速下 RSM 模型；(d) 60m/s 风速下 k-ω SST 模型；(e) 60m/s 风速下 k-ε 模型；(f) 60m/s 风速下 RSM 模型；(g) 90m/s 风速下 k-ω SST 模型；(h) 90m/s 风速下 k-ε 模型；(i) 90m/s 风速下 RSM 模型

2. 壁面热通量：风速提升强化换热，RSM 模型高风速下换热最显著

热通量代表叶片表面的热量传递速率，冷却风速越高，壁面热通量越大，热量带走越迅速：从 30m/s 提升至 90m/s，所有模型的叶片最大热通量均显著增加，前缘、冷却孔等核心区域的换热效率大幅提升，部分区域热通量翻倍。

值得注意的是，部分区域出现负热通量，代表这些区域叶片未实现散热，反而吸收热量，主要原因是局部流场回流或冷却效率不足，而高风速能有效减少这类区域的面积。

三种模型中，高风速下 RSM 模型的热通量最高，能捕捉复杂湍流带来的强化换热效应，前缘和冷却孔附近的换热效果最显著；k-ω SST 模型的热通量分布最均匀，无明显的换热薄弱区；k-ε 模型的热通量整体偏低，即便是 90m/s 高风速，最大热通量仍低于前两者，换热预测偏保守。

图 5 燃气轮机叶片的热通量（W/m²）分布注：(a) 30m/s 风速下 k-ω SST 模型；(b) 30m/s 风速下 k-ε 模型；(c) 30m/s 风速下 RSM 模型；(d) 60m/s 风速下 k-ω SST 模型；(e) 60m/s 风速下 k-ε 模型；(f) 60m/s 风速下 RSM 模型；(g) 90m/s 风速下 k-ω SST 模型；(h) 90m/s 风速下 k-ε 模型；(i) 90m/s 风速下 RSM 模型

3. 湍流动能：集中于流道弯管区，高风速加剧湍流程度

湍流动能是描述湍流强度的核心指标，本次研究发现，湍流动能主要集中在冷却通道的弯管、分叉处，这些区域气流流动方向突变，易出现流场分离和再附，湍流混合剧烈；而流道直管段的湍流动能较低，气流流动稳定。

风速对湍流动能的影响显著，风速越高，湍流动能越大：30m/s 低风速下，各模型的湍流动能整体偏低，仅弯管区有轻微湍流；60m/s 和 90m/s 风速下，湍流动能大幅提升，90m/s 时部分弯管区湍流动能是 30m/s 的数倍，湍流混合更剧烈，换热效率也随之提升。

三种模型中，RSM 模型能捕捉到最显著的湍流动能变化，高风速下弯管区湍流动能最高；k-ω SST 模型的湍流动能分布更合理，贴合实际流动状态；k-ε 模型则低估了湍流强度，各风速下的湍流动能均为三者最低，这也是其换热预测偏差的重要原因。

图 6 燃气轮机叶片的湍流动能（m²/s²）分布注：(a) 30m/s 风速下 k-ω SST 模型；(b) 30m/s 风速下 k-ε 模型；(c) 30m/s 风速下 RSM 模型；(d) 60m/s 风速下 k-ω SST 模型；(e) 60m/s 风速下 k-ε 模型；(f) 60m/s 风速下 RSM 模型；(g) 90m/s 风速下 k-ω SST 模型；(h) 90m/s 风速下 k-ε 模型；(i) 90m/s 风速下 RSM 模型

4. 流速流场：高风速易引发流场分离，k-ω SST 模型流场最稳定

冷却空气的流速流场直接决定冷却效果的均匀性，30m/s 低风速下，所有模型的气流流场均较为稳定，流线平滑，无明显的流场分离和回流；而当风速提升至 60m/s 和 90m/s 时，气流在冷却通道弯管区易出现流场分离和旋流，湍流程度加剧。

三种模型的流场预测差异明显：k-ω SST 模型的流场最稳定，即便 90m/s 高风速，流场分离区域也较小，气流流动更可控，能在强化换热的同时，减少压力损失；RSM 模型能捕捉到复杂的旋流和分离流场，流线最复杂，湍流混合最剧烈；k-ε 模型的流场则相对平缓，对高风速下的流场分离预测不足。

图 7 燃气轮机叶片的速度（m/s）流线分布注：(a) 30m/s 风速下 k-ω SST 模型；(b) 30m/s 风速下 k-ε 模型；(c) 30m/s 风速下 RSM 模型；(d) 60m/s 风速下 k-ω SST 模型；(e) 60m/s 风速下 k-ε 模型；(f) 60m/s 风速下 RSM 模型；(g) 90m/s 风速下 k-ω SST 模型；(h) 90m/s 风速下 k-ε 模型；(i) 90m/s 风速下 RSM 模型

5. 叶片展向温度与压力：中段冷却最优，高风速伴随压力损失增加

研究将叶片从根部到叶尖的展向位置归一化为 0-1（0 为根部，1 为叶尖），分析了展向的温度和压力变化规律：温度方面，所有模型和风速下，叶片展向 0.4-0.6 中段区域温度最低，冷却效果最优，而根部和叶尖区域温度偏高，主要原因是中段冷却流道分布更密集，气流流动更充分；高风速能让全展向温度整体降低，其中 k-ω SST 模型在全展向的温度均为最低，RSM 模型叶尖区域温度略高。

压力方面，风速越高，叶片整体压力越大，90m/s 高风速下的压力远高于 30m/s 低风速；展向 0.5 位置出现明显的压力骤降，代表该区域气流流动受阻，易出现流场分离；RSM 模型的压力分布变化最显著，高风速下部分区域压力偏高，而 k-ω SST 和 k-ε 模型的压力分布更平稳。

同时研究发现，风速提升带来冷却效果提升的同时，会伴随压力损失的增加，90m/s 高风速下，冷却通道的压力损失相比 30m/s 增加了 15%，这意味着需要消耗更多的辅助动力来推动冷却空气，存在冷却效率和运行成本的权衡。

图 8 不同风速和湍流模型下，无量纲温度 T/Tin 沿归一化展向的变化注：横坐标为归一化展向位置，纵坐标为无量纲温度，包含 30m/s、60m/s、90m/s 三种风速下 k-ω SST、k-ε、RSM 模型的变化曲线。

<四、关键结论：模型和风速怎么选？这份工程建议请收好>

本次研究通过系统的 CFD 仿真，明确了冷却风速和湍流模型对燃气轮机叶片冷却性能的影响规律，同时对比了三种湍流模型的优劣，为工程实际提供了清晰、实用的参考建议，核心结论和建议如下：

1. 冷却风速的选择：优先选 60m/s 中风速，兼顾冷却和经济性

高风速（90m/s） 能让叶片温度降幅达 7%，冷却效果最优，但会带来 15% 的压力损失增加，导致辅助动力消耗上升，运行成本增加；低风速（30m/s） 压力损失小，但冷却效果差，叶片高温区明显，易引发材料损坏；中风速（60m/s） 是最优选择，能在实现较好冷却效果的同时，将压力损失控制在合理范围，兼顾冷却效率和设备运行经济性。

2. 湍流模型的选择：常规仿真选 k-ω SST，复杂流场选 RSM

三种湍流模型的性能和适用场景差异显著，工程中需根据仿真需求选择，表 1 燃气轮机冷却分析的湍流模型对比清晰展现了各模型的特点：

k-ω SST 模型

工程常规仿真的首选，壁面预测精度高，能有效捕捉逆压梯度流动，温度预测偏差低于 3%，同时计算成本适中，在精度和效率之间实现了最优平衡，适合内部冷却通道的主流仿真分析；
k-ε 模型

适合快速设计迭代的粗略仿真，计算效率最高，但壁面预测精度低，换热效率低估达 10%，不适合高精度的冷却性能分析，仅可用于完全湍流的外部流场仿真；
RSM 模型

适合复杂冷却结构的精细化仿真，能精准捕捉各向异性湍流和流场分离效应，高风速下换热预测最精准，但计算成本最高，比 k-ω SST 模型高 5%，收敛时间长 40%-50%，仅在流道结构复杂、对仿真精度要求极高时使用。

表 1 燃气轮机冷却分析的湍流模型对比

湍流模型	核心优势	主要局限	最佳适用场景
k-ω SST	壁面预测精度高，适配逆压梯度流动，兼顾计算成本和精度	计算成本略高于 k-ε 模型，无法完全捕捉湍流的各向异性效应	存在强近壁面效应、计算资源适中的内部冷却通道仿真
k-ε	计算效率高，工业应用广泛，对高雷诺数完全湍流流动仿真稳定	壁面预测精度低，换热效率低估达 10%	完全湍流外部流场、近壁面精度要求低的快速设计迭代
RSM	湍流表征最详细，精准捕捉各向异性效应，适配分离流场	计算成本最高，需精细化网格保证仿真稳定性	存在回流区、湍流各向异性显著的复杂冷却几何结构

3. 冷却通道设计建议：优化流道结构，减少压力损失

冷却通道的弯管、分叉处是湍流最剧烈的区域，也是压力损失的主要来源，设计时应采用平滑过渡的弯管结构，减少尖锐转角，在强化湍流换热的同时，降低流场分离和压力损失；同时可优化冷却流道的分布，增加叶片中段流道密度，保持根部和叶尖的流道通畅，提升全叶片冷却的均匀性。

<五、研究展望：这些方向，未来值得深入探索>

本次研究为燃气轮机叶片冷却设计提供了重要参考，但仍存在一些未涉及的内容，未来可从以下方向进一步深入，让仿真结果更贴合燃气轮机的实际工作状态：

引入旋转效应

实际燃气轮机工作时，叶片处于高速旋转状态，会产生科里奥利力和二次流，可能改变冷却通道内的流场和换热特性，未来仿真需纳入旋转效应，提升结果的真实性；
开展瞬态仿真

本次研究为稳态仿真，而实际燃气轮机的温度、压力存在动态变化，瞬态仿真能捕捉短期波动对冷却性能的影响，更贴合实际运行工况；
结合实验验证

本次研究为数值仿真，未来可开展物理实验，通过实测数据验证 CFD 仿真结果，进一步优化湍流模型的参数，提升仿真精度；
探索复合冷却技术

将内部冷却通道与气膜冷却、冲击冷却等技术结合，分析复合冷却方式的性能，为更高温度下的叶片冷却提供解决方案。