涡轮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使用的软件有哪些？
超越伯努利——水力涡轮理论的百年历程
水力涡轮的基础方程来自欧拉的涡轮机械方程（1754年）。这个简单的原理是"流体角动量变化决定涡轮输出"。但实际设计中，需要考虑叶片周围的粘性边界层、转轮的离心力场、尾水管的扩散效率等因素，仅用一维理论无法完成设计。在CFD普及之前，设计人员结合势流理论和经验性损失修正系数（损失图），通过重复手算和模型试验来找到最优叶型。如今将这些经验知识转换为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分析中最难的问题之一是处理旋转转轮与固定导叶之间的流体力学干涉。冻结转子法计算快速但忽视了动静叶干涉的非定常效应。滑动网格法精度高但计算成本是冻结转子的5～10倍。实务判断标准是"所评估物理量的非定常性"——效率图预测用冻结转子通常足够，但压力脉动或疲劳评估需要滑动网格。
涡轮CFD的实务应用
涡轮叶列网格
涡轮叶列的网格与压缩机相同吗？
基本结构相同，但涡轮特有的注意事项。
- 后缘的薄度：涡轮叶片的后缘非常薄（0.3～0.8mm）。O网格后缘周围需要充分的单元
- 冷却孔：L2/L3模型中需要冷却孔周围的局部细化
- 叶片表面跨音速区域：吸力面的超音速补丁和后缘冲击波的分辨
后缘0.3mm的话网格会很细吧。
后缘的O网格在径向至少需要10个单元，后缘下游的尾迹区域网格也要细化。TurboGrid的trailing edge cutoff功能可以控制后缘形状。
跨音速涡轮叶列
涡轮的流动会达到超音速吗？
HP涡轮的叶间马赫数达到1.1～1.3。吸力面加速到超音速后，从后缘发射斜冲击波。这个后缘冲击波系统撞击隣近叶片是CFD精度的关键。
冲击波分辨需要多细的网格？
垂直于冲击波方向的单元尺寸应小于弦长的0.5%，冲击波前后至少需要10个单元。使用自适应网格细化（AMR）将网格集中在冲击波位置也很有效。Fluent和STAR-CCM+的AMR功能可用。
性能预测精度
涡轮CFD的精度如何？
| 指标 | 精度 |
|---|---|
| 段效率（多级） | ±0.5～1.5百分点 |
| 叶片表面压力分布 | 良好（与实验定性一致） |
| 叶片表面热传递系数 | ±10～20%（乱流模型依赖） |
| 后缘冲击波位置 | 弦长的±2% |
奥只见大坝的CFD——国产水力与仿真的历史
日本水力发电CFD的历史与1980年代大型水力开发齐步。东芝、日立、三菱重工各自从势流求解器开始，1990年代迁移到雷诺平均NS方程求解。如今在奥只见发电站（560MW）等大型设备更新中，CFD用于设计最终确认，CFD预测效率与模型试验（1/10缩尺）的差异已控制在±0.3%以内。模型试验与CFD的相关性确立成为按照国际标准（IEC 60041）判断实机设计适用性的依据。
涡轮CFD的软件对比
蒸汽涡轮的CFD
燃气涡轮和蒸汽涡轮的CFD有什么区别？
蒸汽涡轮LP级膨胀过程中蒸汽会越过饱和线，变成湿蒸汽（两相流）。湿蒸汽损失的预测是蒸汽涡轮特有的课题。
湿蒸汽如何处理？
欧拉-拉格朗日法是通用方法。蒸汽（连续相）用欧拉法求解，水滴（分散相）作为拉格朗日粒子追踪。对核生成产生的水滴和水滴的成长进行建模。
CFX的湿蒸汽模型
CFX有蒸汽涡轮专用模型吗？
CFX内置有湿蒸汽模型。可以考虑非平衡凝聚（自发核生成），从过冷度预测威尔逊线。还可以使用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市场由ANSYS和西门子（STAR-CCM+）两家主导，但这个领域也有有趣的利基玩家。瑞士的ANDRI Engineering公司的CFDTURBO专注于水力涡轮专用功能，在欧洲中小水力制造商中根基深厚。OpenFOAM生态中turboBlade和OpenTURNS组合在研究机构中广泛应用。日本国内的电力中央研究所（CRIEPI）和国内涡轮制造商持有专有的分析平台，通常采用商用工具混合运用。
涡轮CFD的前沿研究
叶列干涉
涡轮的动叶和静叶干涉如何评估？
通过滑动网格非定常计算，直接捕捉上游喷嘴尾迹撞击下游转子叶面的过程。尾迹引起的叶面压力脉动是HCF（高周期疲劳）的主要原因。
幽灵效应
幽灵效应是什么？
多级涡轮中前级喷嘴尾迹撞击后级喷嘴前缘或通过叶间时，后级性能发生变化的现象。喷嘴间的周向相对位置（幽灵位置）优化可改善段效率0.3～1百分点。
CFD如何评估幽灵效应？
改变喷嘴的相对位置5～10个阶段，进行滑动网格非定常计算，比较各位置的时间平均效率。计算成本大但设计效果显著。
叶片表面过渡
涡轮叶片表面过渡能用CFD预测吗？
涡轮叶片表面由上游尾迹的周期通过控制"旁路过渡"。Gamma-Theta过渡模型可以预测稳定的旁路过渡，但为了精确捕捉尾迹诱导过渡的非定常效应，需要滑动网格+过渡模型的组合。
涡轮冷却的最前沿
叶片冷却研究的最新趋势？
弗朗西斯水轮机的恶梦——涡绳被CFD解明的日子
弗朗西斯水轮机部分负荷运行时产生的"涡绳"在尾管内形成螺旋空洞，导致振动、噪音和出力波动。1990年代之前涡绳预测依赖经验法则，但2010年代起LES和VOF法组合的CFD普及后，首次实现了定量预测。瑞士EPFL与ABB的联合研究成功将涡绳发生频率的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分析的Tips
CHT（共役热传递）分析的诀窍请指教。
| Tips | 详细 |
|---|---|
| 固体网格匹配 | 流体-固体界面节点一致时精度提高 |
| 固体导热率 | 镍基超合金：11～25 W/(mK)，需要考虑温度依赖 |
| TBC处理 | 用Thin Wall BC模拟薄的隔热涂层 |
| 内部冷却通道 | 可用1D流网模型简化（CFX的Boundary Source Term） |
| 收敛判定 | 叶片表面温度波动在±1K以内时稳定 |
空蚀漏诊——水力涡轮CFD的危险误判
水力涡轮CFD分析中空蚀预测失误最常见的原因是"忽视蒸汽压的温度依赖"。水温10℃和20℃的蒸汽压相差2.3倍（1.2kPa vs 2.8kPa），在20℃夏季设计的转轮在10℃冬季可能出现意外空蚀。另一个易被忽视的是"局部低压峰值"——平均压力系数分布看似没问题，但在叶前缘或吸力面某个仅2～3网格宽的局部区域超过空蚀条件。网格分辨率检验是必须的。
有益
详细
错误