涡轮CFD分析
涡轮CFD的理论基础
概览
涡轮CFD与压缩机有什么区别？
涡轮是从流体中提取能量的一方。流动加速，因此不易发生像压缩机那样的大规模分离。相反，叶片面冷却、二次流损失、跨音速冲击波成为主要课题。
段功和等熵效率
涡轮功率如何表示？
从欧拉方程推导输出和效率。
$h_{02s}$ 是等熵膨胀后的焓。燃气轮机HP段 $\eta_{is}=90\sim92\%$、LP段 $88\sim90\%$ 是现代设计水准。
叶片负荷系数
叶片负荷大小如何评估？
Zweifel叶片负荷系数是标准。
$s$：节距、$c_x$：轴向弦长。$Z_w \approx 0.8$ 是传统最优值，但近年高负荷设计中 $Z_w > 1.0$ 也在研究。
软件选择
涡轮CFD使用什么软件？
涡轮效率理论的基础——兰金循环与热效率（1859年）
蒸汽涡轮热效率定量化的兰金循环理论，由苏格兰工程师William Rankine（1820-1872）确立。相对于卡诺效率，实际蒸汽涡轮系统能达到的效率用兰金循环计算求得的方法至今仍是发电厂设计的基础。Rankine在1859年发表蒸汽机热力学分析论文至今165年后的今天，最新的燃气轮机联合循环（GTCC）发电已达到64%的热效率。这一提升的约1/3来自CFD叶型优化（包括高温域冷却技术），160年前的理论框架与现代CFD结合后继续发挥作用。
涡轮CFD的数值计算方法
叶片面冷却的重要性
涡轮叶片冷却在CFD中如何处理？
HP涡轮入口气体温度达1500～1800℃，远超叶片材料耐热极限（镍基超合金约1000℃）。通过内部冷却通道和膜冷却降低叶片面温度。
冷却模型的阶层
冷却如何融入CFD？
精度与成本的权衡中有多层次方案。
| 等级 | 模型 | 计算成本 | 精度 |
|---|---|---|---|
| L0 | 无冷却流（绝热壁） | 最低 | 无冷却基准评估 |
| L1 | 源项（质量/能量注入） | 低 | 膜冷却的粗略估算 |
| L2 | 离散孔（单独冷却孔BC） | 中 | 膜效率的定量评估 |
| L3 | 分辨冷却孔（孔网格化） | 高 | 最高精度但工作量大 |
| L4 | CHT（流体+固体耦合） | 最高 | 预测叶片内温度分布 |
L3或L4现实可行吗？
单叶片的L3/L4已作为单独计算（Singleton）实现工程化。STAR-CCM+的CHT在该用途中评价高。多级的L3/L4目前仍处研究阶段。
膜冷却效率
膜冷却效果如何评估？
用绝热膜冷却效率定义。
$T_g$：主流气体温度、$T_{aw}$：绝热壁面温度、$T_c$：冷却空气温度。$\eta_f = 0$ 为无冷却，$\eta_f = 1$ 为完全冷却。CFD中计算叶片面绝热壁面温度进而得出。
燃气涡轮叶片CFD数值设置——高温燃烧气体物性与传热模型选择
燃气涡轮叶片CFD分析中，作动气体（高温燃烧气体）物性的准确设置决定传热预测精度。燃烧后气体是CO2、H2O、N2、O2混合物，需用多项式近似或WSGG模型设置温度相关的比热Cp(T)、粘性mu(T)、热导率k(T)。常见的"取1500K空气物性值为常数"近似在高温侧粘性被低估10～15%，这是叶片表面边界层厚度和HTC（传热系数）预测误差的根源。NASA的CHT（共轭传热）基准实验表明，物性用多项式合理近似的CFD能控制叶片面HTC分布与实验误差在±8%内，验证了物性设置精度对高温叶片CFD的重要性。
涡轮CFD的实务应用
涡轮叶列的网格
涡轮叶列网格与压缩机相同吗？
基本结构相同，但涡轮特有的注意事项。
- 后缘细薄：涡轮叶片后缘非常细薄（0.3～0.8mm）。O-grid后缘周围需充分的网格单元
- 冷却孔：L2/L3模型中冷却孔周围需局部细化网格
- 叶片面跨音速区域：吸压面超音速区域与后缘冲击波的分辨
后缘0.3mm，那网格需要非常细密吧。
后缘O-grid半径方向至少10个单元，后缘下游尾迹区域网格也要细化。TurboGrid的trailing edge cutoff功能可控制后缘几何。
跨音速涡轮叶列
涡轮流动会超音速吗？
HP涡轮叶间马赫数达1.1～1.3。吸压面超音速加速后，后缘射出斜冲击波。该冲击波入射到相邻叶片的Trailing Edge Shock System的精准预测决定CFD精度。
冲击波分辨需多细的网格？
垂直冲击波方向的单元尺寸应为弦长0.5%以下，冲击波前后至少10个单元。自适应网格细化（AMR）可将网格集中在冲击波位置。Fluent和STAR-CCM+的AMR功能可用。
性能预测精度
涡轮CFD精度如何？
| 指标 | 精度 |
|---|---|
| 段效率（多级） | ±0.5～1.5百分点 |
| 叶片面压力分布 | 良好（与实验定性一致） |
| 叶片面传热系数 | ±10～20%（湍流模型相关） |
| 后缘冲击波位置 | 弦长±2% |
工业用蒸汽涡轮CFD定期诊断——密封磨损与间隙变化的性能影响
工业蒸汽涡轮在10～20年长期运行中会经历密封磨损、叶片侵蚀、叶顶间隙扩大，导致段效率经年衰减0.5～2百分点。用CFD定量评估该性能衰减并决定最优检修时机的"CFD诊断"已在大型火力电厂实现工程化应用。分析中将定期检查数据（密封间隙实测值、叶片形状变化）反映到CFD模型，与设计值对比特定损失增加。某日本电力公司的事例中，CFD诊断预测出某一级密封间隙已达设计值的3倍，经实机确认后仅进行部分密封更换，相比全面检修削减成本60%，该案例已在电力中央研究所报告中公开。
涡轮CFD的软件比较
蒸汽涡轮的CFD
燃气涡轮和蒸汽涡轮的CFD有什么不同？
蒸汽涡轮LP段膨胀过程中蒸汽超过饱和线，变成湿蒸汽（两相流）。湿蒸汽损失预测是蒸汽涡轮特有课题。
湿蒸汽怎么处理？
欧拉-拉格朗日法很常见。蒸汽（连续相）用欧拉法求解，水滴（分散相）作拉格朗日颗粒追踪。模型化成核生成产生水滴和水滴长大。
CFX的湿蒸汽模型
CFX有蒸汽涡轮专用模型吗？
CFX搭载了湿蒸汽模型。考虑非平衡凝结（自发成核），可从过冷度预测Wilson线。还可利用IAPWS-IF97水蒸汽物性表。
| 功能 | CFX | STAR-CCM+ | OpenFOAM |
|---|---|---|---|
| 湿蒸汽模型 | 湿蒸汽模型 | 拉格朗日+自定义 | wetSteamFoam（研究用） |
| 物性表 | IAPWS-IF97 | IAPWS-IF97 | janaf或自定义 |
| 水滴追踪 | 欧拉多相 | 拉格朗日 | 受限 |
长叶片的课题
蒸汽涡轮LP最后一级叶片很长吧？
现代蒸汽涡轮LP最后一级叶片超过1.5m，叶尖周速超过马赫1.5。从叶根到叶尖流动状态完全不同，3D设计必不可少。叶尖处超音速+湿蒸汽，叶根处亚音速+干蒸汽。
沿展向相差这么大？
沿高度方向在50～100个展向截面改变叶型的3D堆叠是标准设计手法。CFD需精密评估展向全长的效率分布。
燃气涡轮CFD工具——商业工具与企业内部代码的产业生态
全球燃气涡轮制造商拥有独特的内部CFD代码和商业工具相结合的生态。Rolls-Royce的"Hydra"、GE Aerospace的"WIND"、MTU的"TRACE"各有50年以上的开发历史，作为从叶型设计到认证试验都信赖的企业内部工具。商业工具方面，ANSYS CFX是航空发动机叶型分析的事实标准，NUMECA FINE/Turbo在欧洲航空航天企业中常用。学术研究机构还使用SU2（斯坦福大学）和OpenFOAM，但实机设计认证文件用的几乎都是内部代码或商业验证代码。工具的信赖性取决于数十年的实绩验证——再高功能的新工具若缺乏充分实绩就不会被采用，这是现实。
涡轮CFD的前沿研究
叶列干涉
涡轮动叶与静叶的干涉如何评估？
用滑动网格非定常计算，上游喷嘴的尾迹冲击下游转子叶片的过程直接捕捉。尾迹引发的叶片面压力波动是HCF（高周期疲劳）的主因。
间歇效应
间歇效应是什么？
多级涡轮中前级喷嘴尾迹在后级喷嘴前缘打中或穿过叶间，导致后级性能变化的现象。喷嘴之间周向相对位置（间歇位置）优化可改善段效率0.3～1百分点。
用CFD评估间歇效应怎么做？
改变喷嘴相对位置分别进行5～10次滑动网格非定常计算，比较各位置的时间平均效率。计算成本大但设计效果显著。
叶片面转捩
涡轮叶片面转捩能用CFD预测吗？
上游尾迹周期性通过时，涡轮叶片面呈现"旁路转捩"为主。Gamma-Theta转捩模型可预测定常旁路转捩，但要精确捕捉尾迹诱导转捩非定常效应需要滑动网格+转捩模型的组合。
涡轮冷却的最前线
叶片冷却研究的新进展？
涡轮叶片最先进冷却——AI与深度强化学习的膜孔配置优化与CFD验证
为最大化燃气涡轮叶片冷却效率，用AI与深度强化学习（DRL）对膜冷却孔配置进行优化成为研究热点。数百个孔的位置、直径、倾斜角的组合形成事实上无限的设计空间，传统参数化优化难以应对。GE与NASA合作研究表明，强化学习代理利用CFD评估的奖励信号进行10万余次迭代，找到了冷却空气用量削减20%同时维持冷却效果的孔位方案。该设计的试验叶片通过3D打印试制并验证，CFD预测与实验偏差控制在±5%以内。
涡轮CFD故障排除
传热系数的不匹配
CFD的叶片面传热系数经常与实验不符…
叶片面传热预测是CFD最困难的项目之一。梳理典型不匹配要因。
| 要因 | 影响 | 对策 |
|---|---|---|
| 湍流模型 | SST误差±15～25% | 添加Gamma-Theta转捩模型 |
| y+的管理 | y+ > 2时传热被低估 | 确保y+ < 1 |
| 入口湍流强度 | 燃烧器出口TI: 10～20% | 反映实验值，默认5%过小 |
| 自由流湍流衰减 | 到叶前缘的TI衰减 | 正确设置湍流长度尺度 |
入口湍流强度竟然有10～20%？
燃烧器出口还存有涡流，湍流强度较高。用5%计算会导致前缘附近传热严重低估。
后缘传热
后缘附近的传热特别合不上。
后缘是尾迹区与叶片面边界层交叉的复杂流场。RANS往往难以精确再现后缘近处湍流结构。用SAS或SDES解析该区域非定常涡可改善。
CHT分析的技巧
CHT（共轭传热）分析有什么诀窍？
| 技巧 | 详情 |
|---|---|
| 固体网格的对应 | 流体-固体界面节点一致能提高精度 |
| 固体热导率 | 镍基超合金：11～25 W/(m·K)，考虑温度相关性 |
| TBC的处理 | 用薄壁BC模拟薄隔热涂层 |
| 内部冷却通路 | 可用1D流网模型简化（CFX的Boundary Source Term） |
| 收敛判定 | 叶片面温度波动±1K以内为稳定 |
涡轮CFD叶片面传热为实验2倍——膜冷却孔周围网格解析度不足
燃气涡轮叶片CFD分析中"叶片面传热系数（HTC）比实验高2倍左右"的问题，多数源于膜冷却孔周边局部网格解析度不足。冷却孔直径D（典型0.5～1mm）附近，孔出口下游10D以内的单元尺寸需为孔径的至少1/20，不满足此条件的粗网格会严重误算孔出口速度和温度分布。此外孔周生成的肾脏形涡对（CRVP）是绝热膜冷却效率的主要控制因素，但涡心径仅D/2左右，解析度不足时涡会消失。实用对策是在孔周围应用AMR（自适应网格细化），并对孔间距和喷嘴叶片阵列都进行热流收敛确认。
相关主题
了
更详细
错误