涡轮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使用什么软件？
朗肯循环和CFD——热力学与流体力学的接点
蒸汽涡轮的CFD分析中，不仅需要纳维-斯托克斯方程，还需要热力学状态方程。理想气体近似在高压阶段完全失效，必须引入IAPWS-IF97（水和水蒸气的国际标准状态方程）。该状态方程由数十个系数组成的复杂多项式构成，集成到CFD求解器后计算成本增加约20%。朗肯循环热效率提高0.1%在大型发电站相当于年度数亿日元的燃料节约，因此状态方程精度与经济效益直接相关。
涡轮CFD的数值计算方法
叶片冷却的重要性
涡轮叶片的冷却在CFD中如何处理？
HP涡轮进口燃气温度达到1500～1800℃，远超叶片材料耐热极限（镍基超合金约1000℃）。通过内部冷却通道和膜冷却来降低叶片表面温度。
冷却模型的层级
如何在CFD中引入冷却？
精度和成本的权衡中存在多个层级。
| 级别 | 模型 | 计算成本 | 精度 |
|---|---|---|---|
| L0 | 无冷却流（绝热壁） | 最低 | 无冷却的基准评价 |
| L1 | Source Term（质量/能量注入） | 低 | 膜冷却概估 |
| L2 | Discrete Hole（独立冷却孔BC） | 中 | 膜冷却效率定量评价 |
| L3 | Resolved Cooling Holes（孔网格化） | 高 | 最高精度但工作量大 |
| 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的诞生——雾中的数值求解法
蒸汽涡轮最后阶段的蒸汽凝聚产生液滴，冲击叶片表面引起冲蚀。这种"湿蒸汽"的数值处理方法在1990年代才确立，需要结合经典均质凝聚模型和纳米级液滴核生成理论。即使在今天，湿蒸汽CFD的计算成本仍比单相气体高3～5倍，液滴粒径分布的假设改变就可能导致效率预测偏差超过1%。
涡轮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% |
高温高压的洗礼——蒸汽涡轮实机计测的困难
蒸汽涡轮高压阶段常见温度600℃、压力240bar以上的环境，用皮托管计测动叶内部压力分布等"常规上不可能"的情况屡见不鲜。从1980年代开始尝试精密内部流动计测，但实用性仅限于出口总压、总温的级间计测。CFD验证的现场中，计测精度±0.5%相对于CFD误差1～2%，"计测和模型都存在误差"成为前提，统计性验证方法成为实务智慧。
涡轮CFD的软件对比
蒸汽涡轮的CFD
燃气涡轮和蒸汽涡轮的CFD有什么区别？
蒸汽涡轮LP阶段膨胀过程中蒸汽超过饱和线，成为湿蒸汽（二相流）。湿蒸汽损失的预测是蒸汽涡轮特有的课题。
湿蒸汽怎么处理？
欧拉-拉格朗日法最为普遍。蒸汽（连续相）欧拉法求解，水滴（分散相）作为拉格朗日粒子追踪。需要建模核生成产生的水滴和水滴增长。
CFX的湿蒸汽模型
CFX有蒸汽涡轮专用模型吗？
CFX搭载Wet Steam Model。考虑非平衡凝聚（Spontaneous Nucleation），能从过冷度预测Wilson线。也可利用IAPWS-IF97水蒸汽物性表。
| 功能 | CFX | STAR-CCM+ | OpenFOAM |
|---|---|---|---|
| 湿蒸汽模型 | Wet Steam Model | 拉格朗日+自定义 | wetSteamFoam（研究用） |
| 物性表 | IAPWS-IF97 | IAPWS-IF97 | janaf或自定义 |
| 水滴追踪 | 欧拉多相 | 拉格朗日 | 受限 |
长叶片的课题
蒸汽涡轮LP最后阶段叶片很长吧？
最新蒸汽涡轮LP最后阶段叶片长度超过1.5m，叶端周速超过马赫1.5。从叶根到叶端的流动条件完全不同，3D设计不可或缺。叶端为超音速+湿蒸汽，叶根为亚音速+干燥蒸汽。
沿展向差异这么大吗。
沿叶高变更50～100个展向断面的翼型是标准设计方法（3D Stacking）。需要CFD精密评估整个展向的效率分布。
蒸汽涡轮CFD工具选择——ANSYS vs CFX vs OpenFOAM
蒸汽涡轮专家长期支持ANSYS CFX的理由是转动机械专用Mixing Plane接口和湿蒸汽模型多年前就被集成。Fluent在2010年代后期开始追赶，现在功能上基本相当。另一方面OpenFOAM的WetSteamFoam求解器存在，但IAPWS状态方程的完全实现在许多情况下需要用户补充。考虑企业许可证成本（年费5000～2000万日元），越来越多机构在考虑向开源迁移，但支持体系的差异成为迁移壁垒。
涡轮CFD的前沿研究
叶列干涉
涡轮动叶和静叶的干涉如何评估？
Sliding Mesh非定常计算直接捕捉上游喷嘴尾流冲击下游转子叶片的过程。尾流引起的叶片表面压力波动是HCF（高周期疲劳）的主要原因。
时钟效应
时钟效应是什么？
多级涡轮中前级喷嘴尾流冲击后级喷嘴前缘还是叶间，取决于喷嘴之间的周向相对位置（时钟位置），后级性能会改变。喷嘴相对位置优化可使级效率改善0.3～1个百分点。
CFD如何评估时钟效应？
改变喷嘴相对位置5～10个阶段，对各位置做Sliding Mesh非定常计算，比较时间平均效率。计算成本很大，但设计效果也显著。
叶片表面转捩
涡轮叶片表面转捩能用CFD预测吗？
上游尾流周期通过导致涡轮叶片表面"旁路转捩"占主导。Gamma-Theta转捩模型可预测稳态旁路转捩，但要精确捕捉尾流诱导转捩的非定常效应，需要Sliding Mesh加转捩模型的组合。
涡轮冷却的最前线
叶片冷却研究的发展趋势是？
蒸汽涡轮叶型与AI——1亿通设计搜索
2020年代，GE和Siemens Energy使用AI代理模型，以比传统快100倍以上的速度搜索蒸汽涡轮叶型。每次CFD计算需数小时，学习型神经网络在0.01秒内近似，使1亿通设计候选现实可评估。这提升动叶热效率较传统设计增加0.3～0.5个百分点。AI并非取代CFD工程师，而是大幅扩展可探索的设计空间。
涡轮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/(mK)，需考虑温度相关性 |
| TBC处理 | 用Thin Wall BC模拟薄的隔热涂层 |
| 内部冷却通道 | 用1D流动网络模型简化（CFX的Boundary Source Term） |
| 收敛判定 | 叶片表面温度变动在±1K以内稳定 |
蒸汽涡轮CFD的陷阱——"收敛不了"的背后
蒸汽涡轮CFD残差不下降的典型原因之一是"堵塞附近的计算"。高压阶段叶间流路中蒸汽达到马赫1的喉部，非物理衝击波引起残差振荡。此时下游压力边界条件稍降就导致超音速区扩大，收敛破裂。资深工程师的经验法则是"先做堵塞判定"——用1D计算事先确认范廷系数和马赫数分布，再启动3D分析才是关键。
细节
错误