风摩损失
风摩损失的理论基础
概述
什么是风摩损失？这是一个不太常见的术语。
旋转体与静止壁面间的空气（或流体）摩擦导致的动力损失。高速旋转机械中不可忽视的损失源。特别重要应用于涡轮发电机、飞轮、高速电动机等。
旋转圆盘的摩擦损失
具体的损失有多大？
闭空间内旋转圆盘的摩擦损失表示如下。
$C_M$是力矩系数，是旋转雷诺数 $Re_\theta = \rho \omega r^2 / \mu$ 和间隙比 $s/r$ 的函数。
与$\omega^3$成比例？旋转数的影响很大啊。
是的。转速增加一倍，风摩损失增加8倍。因此在高速旋转机械中，它可能成为主要的损失来源。
流动体制
间隙流动中也有层流和湍流吗？
有。Daily & Nece (1960) 的分类被广泛使用。
| 流动体制 | $Re_\theta$ | 间隙比 | 特征 |
|---|---|---|---|
| I: 层流·合并BL | 低 | 小 | 两侧边界层结合 |
| II: 层流·分离BL | 低 | 大 | 边界层独立，核心旋转 |
| III: 湍流·合并BL | 高 | 小 | 最常见的工业条件 |
| IV: 湍流·分离BL | 高 | 大 | 核心旋转存在 |
工业涡轮机械几乎全部处于流动体制III或IV。
$C_M$ 的经验式
力矩系数有经验公式吗？
流动体制III的情况下，Daily的公式最具代表性。
通过将CFD结果与该经验式进行比较，可以简便地验证计算的合理性。
旋转圆盘的流体力学——von Kármán的圆盘分析（1921年）与Stokes的开创
对静止流体中旋转圆盘流动进行分析的是Theodore von Kármán（1921年）。他不仅分析了转子的Kármán涡，还推导了旋转圆盘边界层（von Kármán边界层）的严格解析解，明确定义了圆盘力矩与涡粘性的关系。Cochran(1934)以更高的精度计算了这个解，Daily & Nece(1960)通过实验扩展了模型，奠定了现代相关式的基础。当代的涡轮机械圆盘腔CFD是建立在100年前von Kármán的分析基础之上，并在此基础上通过数值方法添加了多级、多体旋转、热气侵入等复杂的三维效应。这一点值得特别指出：一位数学物理学家在一个世纪前做出的革命性成果，如今直接应用于喷气发动机冷却设计中。
风摩损失的数值计算方法
CFD模型的构成
如何用CFD计算风摩损失？
将旋转圆盘与静止壁间的间隙空间建模为流体区域。
- 旋转壁面：圆盘表面（无滑移、指定旋转速度）
- 静止壁面：外壳内表面（无滑移、速度为零）
- 入口/出口：如果存在向间隙泄漏的流动，需要设置
如果是轴对称，2D轴对称模型就足够了。如果存在3D效应（螺栓头、冷却孔等），使用扇区模型。
湍流模型的选择
哪种湍流模型比较适用？
推荐SST k-omega。旋转壁面附近的边界层很关键，需要低Reynolds数解法来确保y+ < 1。k-epsilon模型不能准确捕捉壁面附近的旋转效应，不推荐。
| 湍流模型 | 风摩预测精度 | 备注 |
|---|---|---|
| SST k-omega (Low-Re) | ±5～10% | 推荐 |
| k-epsilon + 壁函数 | ±15～25% | 壁面附近精度不足 |
| RSM (Reynolds Stress) | ±3～8% | 捕捉旋转效应，计算成本大 |
| LES | ±2～5% | 最高精度，计算成本大 |
RSM效果更好吗？
在旋转圆盘流中，周向和径向Reynolds应力的各向异性很强。RSM（如BSL-RSM）直接模化这种各向异性，所以在某些情况下精度会优于SST。但计算成本增加1.5～2倍。
网格要求
间隙内需要多少网格？
间隙方向（轴向）40～60个网格单元，径向100～200个。在两侧壁面配置棱柱层，确保y+ < 1。2D轴对称模型只需数万个网格单元。
风摩损失的计算——Daily-Nece相关式与CFD精度对比
旋转圆盘的风摩损失力矩M可以从Daily & Nece(1960)的实验相关式计算得出。四种流动体制（G-Reomega的组合）各有不同的系数式，在湍流域（体制III、IV）中，M与转速的平方、半径的5次方成正比，与Re的1/5次方成反比。这个相关式在发动机和涡轮设计中至今仍被广泛用于初步估算，但当工况条件（G=0.02～0.2、Reomega=10^4～10^7）超出适用范围时，误差可能超过30%。多项CFD对比研究表明，在间隙很小（G<0.02）或圆盘有凹陷、孔洞等实际形状的情况下，Daily-Nece式与CFD结果相比往往低估损失10～20%。对于形状复杂的发动机圆盘，标准做法是先用相关式进行初步估算，然后用CFD进行精密验证。
风摩损失的实务应用
转子系的温度上升
风摩损失会导致温度升高吗？
损失全部转化为热。如果密闭空间冷却不足，温度会持续上升。这在蒸汽涡轮汽室或电机气隙中是个问题。
冷却流量 $\dot{m}_{cooling}$ 越小，温度上升越大。
CFD可以预测温度分布吗？
启用能量方程并进行共轭传热（CHT）计算可以预测转子表面温度分布。但在高温涡轮中，还需考虑辐射传热，不仅仅是对流。
密封流动的耦合
迷宫密封的泄漏流会影响风摩吗？
影响很大。通过迷宫密封的泄漏流进入圆盘间隙时，会改变核心旋转速率，从而改变风摩损失。CFD最好将密封和间隙一体化建模。
与实验的验证
如何验证风摩损失的CFD结果？
| 验证方法 | 精度 |
|---|---|
| Daily & Nece经验式 | ±10%（简单圆盘） |
| 扭矩测量计实测 | 最精确，但实验成本高 |
| 从温度上升反算 | ±15%（热损失估算的不确定性） |
对于简单圆盘，CFD精度为±5～10%，与经验式相符。含有凸起或螺栓头的情况下，如果形状建模准确，可获得±10～15%的精度。
航空发动机的圆盘腔——二次流冷却系统的CFD热设计
喷气发动机高压涡轮中，高温燃气可能侵入旋转圆盘与外壳之间的"圆盘腔"，导致圆盘过热和破损危险。为防止这一点，采用"二次流冷却系统（Secondary Air System: SAS）"，从压气机抽气注入腔内，形成防止热气侵入的清洗流。过多的清洗流会降低发动机效率，过少则圆盘达到危险温度。通过CFD（旋转圆盘+可压缩RANS）来优化圆盘腔内的压力、温度、流量分布，可在保证热安全系数的前提下最小化清洗流量，这已成为航空发动机制造商CFD应用的一个重要领域。
风摩损失的软件比较
损失降低的方法
如何降低风摩损失？
有以下方法。
| 方法 | 原理 | 效果 |
|---|---|---|
| 减小间隙 | 流体体积减少 | 流动体制III→I转变，$C_M$下降 |
| 填充低密度气体 | $P_w \propto \rho$ | 用He替代空气，密度降低1/7 |
| 真空化 | 密度降低的极限 | 大幅降低，但需要密封 |
| 降低表面粗糙度 | 摩擦系数降低 | 研磨表面可降低5～15% |
| 去除转子表面突起 | 压力阻力降低 | 用螺栓头盖降低20～30% |
飞轮应该都是真空的吧？
是的。飞轮储能装置在10 Pa以下真空中运转，使风摩损失近乎为零。但真空密封设计很困难。
CFD的参数化评估
用CFD评估损失降低方案的步骤是？
2D轴对称模型每个工况只需几分钟计算，很适合参数化研究。通过改变间隙比、转速、气体种类，生成$C_M$的特征图，找到最优设计条件。
求解器的选择
哪种求解器适合风摩分析？
圆盘腔CFD的专用工具——国际SAS联盟与CFX的应用实绩
航空发动机的二次流系统（SAS）分析中，产业界与学术界联合开展的国际研究项目对分析方法标准的建立作出了贡献。欧洲的RFICAS（Internal Cooling Flows and Turbine Airfoils Study）、美国的NASA Turbofan SAS研究等。在工具方面，ANSYS CFX在劳斯莱斯、GE、普惠等公司的应用实绩最多，在涡轮机械领域几乎成为事实标准。尤其是CFX的SST湍流模型在圆盘腔内旋流（Swirl）的预测中，相比其他模型精度更为稳定，圆盘腔内旋流系数（Swirl Ratio）预测精度的实绩数据累计到±5%。OpenFOAM的rotating Frame功能也在研究用途中使用，但在需要商业许可认证的设计应用中，CFX仍然占主导地位。
风摩损失的先端研究
泰勒-库埃特流
旋转圆盘间的流动会产生涡吗？
同心圆筒间的流动（泰勒-库埃特流）会产生泰勒涡。当旋转Reynolds数超过泰勒数的临界值时，轴对称的环形涡会自发形成。
泰勒涡如何影响风摩损失？
涡引起的动量传输大幅增加，$C_M$急剧上升。准确预测泰勒涡的转变需要细致的网格和谨慎的湍流模型选择。LES或DNS可以直接分辨涡结构，但RANS虽可捕捉涡的产生，但涡结构的细节会丢失。
轴向流动的泰勒涡
当密封有泄漏流时，泰勒涡如何变化？
轴向流加入后会形成螺旋泰勒涡。流动体制由轴向Reynolds数与旋转Reynolds数的比决定。CFD通过非定常计算追踪这一转变，可以揭示密封内部流动与损失的关系。
最近的研究
风摩损失研究的前沿是什么？
风摩损失的前沿——超高速飞轮与真空外壳设计
飞轮储能（FES: Flywheel Energy Storage）系统中，转子在每分钟40,000～100,000转运转，周速可达500～1,000 m/s，风摩损失成为主要损失部分。真空外壳（0.1～1 Pa）内转体的空气动力损失属于稀薄气体（分子流），普通连续体CFD（Knudsen数Kn>0.1）不适用。需采用DSMC模拟（Direct Simulation Monte Carlo）分子论分析或格子Boltzmann方法（LBM）。飞轮损失最小化涉及超高速转体形状（辐条配置、轮毂形状）的CFD优化，已有报告证实可将风摩损失降低50%的超导磁轴承飞轮研究论文。
风摩损失的故障排除
2D vs 3D模型
风摩分析用2D就够了吗？
看情况。
| 条件 | 推荐 |
|---|---|
| 轴对称转子（光滑表面） | 2D轴对称足够 |
| 有螺栓头或冷却孔 | 3D扇区模型 |
| 附加迷宫密封 | 2D轴对称（密封齿轴对称时） |
| 非轴对称形状（带翼转子） | 3D全模型或扇区 |
2D和3D结果差异有多大？
光滑圆盘的2D和3D差异在5%以内。如果有突出的螺栓头，3D是必需的，否则阻力被严重低估。单个螺栓头可导致$C_M$增加10～30%。
收敛判定
风摩分析的收敛如何判定？
监测转子表面的扭矩。定常计算中，当扭矩波动在0.1%以内时即收敛。
常见的错误
风摩分析常见的错误是？
2. 层流计算：工业机械的Re_theta通常 > $10^6$，肯定是湍流。层流设定会使$C_M$偏低好几个数量级
3. 忽略核心旋转速率：间隙内流体不会与圆盘同速旋转。核心旋转速率（通常0.3～0.5）必须准确预测，需要充分的网格分辨率
4. 省略突起物**：螺栓、导线等突起被忽略会严重低估损失
核心旋转速率是什么？
圆盘间中央（核心）的流体以圆盘角速度的多少比例旋转的指标。间隙越大，核心旋转速率越低（约0.3～0.4），即圆盘的一半以下速度。这个值决定了风摩损失的大小。
风摩损失CFD预测实测的2倍——圆盘间隙流动体制选择错误
涡轮机械的圆盘与外壳间"风摩损失"的CFD预测高于实测2倍的情况，多源于旋转流动体制判别失误。圆盘间流动分为4种体制（Daily & Nece 1960），由间隙比G = s/R（s:间隙、R:半径）与旋转Reynolds数Reomega = omega*R^2/v决定。当G<0.01的狭窄间隙、Reomega>10^5的湍流域（体制IV）时，边界层会合流，速度分布大幅改变。这个体制判别出错会导致风摩力矩预测偏差2～3倍。CFD的基本条件是：几何上准确表现间隙（不能省略或简化），壁面处理采用y+<1的低Reynolds数方法。
相关主题
价值
详细
错误