雷诺输运定理
雷诺输运定理的理论基础
概要
老师，什么是雷诺输运定理？它与NS方程有什么关系？
定理的记述
物质体积 $V_m(t)$（随流体一起运动的体积）内任意示量性物理量 $B$ 的时间变化率可以表示为：
其中 $b = B/m$ 是单位质量的量。RTT将其转换为固定检查体积 $V_{cv}$ 的表达式。
右边的两项分别是什么意思？
第一项是检查体积内 $B$ 的时间变化率（蓄积率）。第二项是通过检查体积表面的 $B$ 的正净流出率（通量）。两者之和等于整个系统 $B$ 的时间变化率。
物质导数与的关系
在微分形式中，RTT对应物质导数。对于任意标量场 $\phi$，
其中物质导数为 $\frac{D\phi}{Dt} = \frac{\partial\phi}{\partial t} + \mathbf{u}\cdot\nabla\phi$。
守恒律的导出
根据代入 $b$ 的不同选择，可以得到各种守恒律。
| $b$ 的选择 | 导出的守恒律 | 得到的方程 | ||
|---|---|---|---|---|
| $b = 1$ | 质量守恒 | 连续方程 | ||
| $b = \mathbf{u}$（速度） | 动量守恒 | NS方程（动量方程） | ||
| $b = e + \frac{1}{2}$ | \mathbf{u} | ^2$（全能量） | 能量守恒 | 能量方程 |
| $b = \mathbf{r}\times\mathbf{u}$（角动量） | 角动量守恒 | 涡轮机械欧拉式 |
所有这些都来自一个定理。太了不起了。
是的。RTT可以说是流体力学守恒律的大一统原理。CFD求解的所有方程的出发点都在这里。
雷诺输运定理的起源——奥斯博恩·雷诺的多彩成就（1842～1912年）
以"雷诺输运定理"命名的奥斯博恩·雷诺（Osborne Reynolds）不仅发现了湍流（1883年），还在流体力学整体上做出了革命性贡献。1895年，他导出了时间平均化的纳维-斯托克斯方程——"雷诺平均纳维-斯托克斯（RANS）方程"，为工程上处理湍流奠定了基础。雷诺输运定理（物质系与控制体积的关系）是连续体力学的普遍变换原理，适用于所有质量、动量、能量和角动量的守恒律。时间快进140年到现在，CFD的控制方程（纳维-斯托克斯）仍基于雷诺导出的形式，他的成就在计算机模拟时代依然闪闪发光。
雷诺输运定理的数值计算方法
RTT与有限体积法的关系
RTT与CFD离散化是怎样连接起来的？
有限体积法就是直接离散化RTT的积分形式。将每个计算单元作为检查体积，计算通过单元面的通量收支。
其中 $F_f$ 是通过面 $f$ 的对流通量，$D_f$ 是扩散通量，$S_i$ 是源项。这就是RTT的离散版本。
移动检查体积（ALE法）
当检查体积本身运动时怎么办？
在ALE（任意拉格朗日-欧拉）方法中，检查体积以速度 $\mathbf{u}_g$ 运动。RTT相应修改为：
关键是对流速度变成 $\mathbf{u} - \mathbf{u}_g$（流体速度与网格速度的差）。动态网格（活塞运动、旋转机械等）就使用这种形式。
几何守恒律（GCL）
在移动网格中必须满足几何守恒律（Geometric Conservation Law）。要求在代入 $\phi = 1$（均匀场）时严格守恒。
如果GCL不满足会怎样？
均匀流无法保持均匀，会出现人为的质量生成或消耗。特别在网格大幅变形的问题（流固耦合、自由表面等）中会变得很严重。商用求解器通常会自动满足GCL，但用户定义的动态网格需要特别注意。
检查体积分析的实务应用
RTT的积分形式也可直接用于从CFD结果计算力和力矩。
除了物体表面的压力和摩擦力积分外，还有"远场力算法"——从围绕物体的检查体积表面的通量计算力。这在航空器阻力分解（压力阻力、摩擦阻力、诱导阻力）中特别有用。
RTT不仅是理论，还直接连接到实际CFD后处理。
控制体积的"切割方法"使计算效率改变10倍
将雷诺输运定理落实到数值解法时，控制体积的设置方式会大幅影响计算成本。在涡轮机械的叶间分析中，每级之间设置"混合平面"来处理静翼-动翼界面的方法，与一次性非定常求解所有级数的方法，计算时间可相差10～100倍。有限体积法教科书上说"处处相同"，但在实际工业设计中，控制体积的切割方法本身就成了技术差异。
雷诺输运定理的实务应用
实践指南
怎样在实务中运用RTT的思想？
RTT是"选定检查体积，应用守恒律"这一强大思维工具。可用于CFD结果的验证和快速估算。
管路的压力损失
以管路入口（截面1）和出口（截面2）为检查面，应用动量RTT。
其中 $F_{\text{wall}}$ 是壁面摩擦力。从CFD结果取得各截面的平均速度和压力，检验是否与此式一致。
喷气的推力
对包围喷气发动机的检查体积应用动量RTT，得到：
用CFD求解喷嘴流动后，从出口截面的流量、速度和压力用这个式子计算推力。
CFD后处理的应用
在CFD后处理中怎样运用RTT？
可以在以下场景活用：
| 应用场景 | 方法 | 求解器操作 |
|---|---|---|
| 力的计算 | 动量通量的面积分 | Fluent: Report > Forces |
| 质量流量确认 | 质量通量的面积分 | Fluent: Report > Fluxes |
| 混合程度评估 | 标量通量的面积分 | 任意截面的标量平均值 |
| 能量收支 | 焓通量的面积分 | 入出口Total Enthalpy差 |
动态网格问题检查清单
使用移动网格时的要点整理如下：
- GCL合规性：用均匀流测试验证质量守恒
- 网格品质监控：在每个时间步验证变形过程中品质未劣化
- 重新网格阈值：斜度超过0.9时重新网格
- 时间步长：考虑网格移动速度的CFL条件
我还以为RTT是抽象的定理，没想到它这么贴近实务。
要有"验证CFD结果的眼光"，非得理解RTT的思想不可。数值出来后，首先检查入出口的通量收支。那就是RTT的实践啊。
流量计校准实际上就是输运定理的教科书问题
现场常用的超声波流量计和科里奥利流量计的校准工作，其实就是将雷诺输运定理的质量守恒应用到现实中。当管路内流速分布不均匀时，截面平均流速与实际体积流量会有偏差。大型水力发电站会测量入口的流速分布，用雷诺输运定理积分，以0.1%的精度管理进入涡轮的流量。"使用定理"会进化成"用定理赚钱"，这就是实践阶段。
雷诺输运定理的软件比较
RTT与动态网格的商用工具支持
移动网格的处理在各求解器中有区别吗？
移动检查体积版RTT（ALE法）的实现在各工具中不同。
动态网格方法的比较
| 方法 | Fluent | CFX | STAR-CCM+ | OpenFOAM |
|---|---|---|---|---|
| 弹簧平滑 | 支持 | 支持 | 支持 | displacementLaplacian |
| 扩散基平滑 | 支持 | 支持 | 支持 | displacementMotionSolver |
| 重新网格 | 局部重新网格 | 无（仅网格变形） | 支持 | 有限制 |
| 重叠（Chimera） | Fluent 2020+ | 不支持 | 支持 | overPimpleDyMFoam |
| 滑移网格 | MRF / Sliding | GGI interface | Sliding Interface | cyclicAMI |
什么是重叠网格？
背景网格与物体周围的移动网格相互重叠的方法。无需受网格变形的约束，适合物体大幅移动、转动的问题（阀门开关、多体运动等）。但重叠区域的插值会导致守恒性略有劣化，需要注意。
旋转机械的处理
旋转机械（泵、涡轮等）作为RTT旋转检查体积版来处理。
| 方法 | 精度 | 成本 | 用途 |
|---|---|---|---|
| MRF （多参考系） | 定常近似 | 低 | 设计探索、初步评估 |
| 混合平面 | 定常近似 | 低 | 多级涡轮机械 |
| 滑移网格 | 非定常、严格 | 高 | 详细设计、动静叶干涉 |
MRF与滑移网格的选用有什么区别？
MRF把旋转域当定常来处理，无法捕捉叶片数×通过频率的非定常效应。当压力脉动和噪音很重要时，非得用滑移网格。计算成本会增加10～100倍，所以通常先用MRF掌握概貌，必要时再用滑移网格做详细分析。
力和力矩的计算
整理各求解器中基于RTT的力算法。
| 计算量 | Fluent | STAR-CCM+ | OpenFOAM |
|---|---|---|---|
| 表面力（压力+摩擦） | Report > Forces | Report > Force | forces functionObject |
| 力矩 | Report > Moments | Report > Moment | forceCoeffs |
| 远场力 | Custom Field Function | Custom Report | 自定义后处理 |
无论用哪个求解器，只要理解RTT的物理，输出值的意义就清楚了。
CFD的质量、动量守恒实现差异——Fluent与OpenFOAM的守恒性方案比较
CFD求解器的数值守恒性（Conservative Property）在各工具的实现中有所不同，影响长时间非定常计算和多相流计算的误差累积。ANSYS Fluent的非结构网格求解器采用有限体积法设计，力求满足离散守恒性（Discrete Conservation），但二阶以上精度的方案会引入非守恒补间。OpenFOAM采用单元中心有限体积法，质量守恒的严格性较高，rhoPimpleFoam通过压缩性质量通量修正妥善处理密度与速度耦合。在多相流（VOF法）中，界面处的质量守恒误差会作为"体积误差"积累，需要根据工具的体积守恒特性规格进行确认，并定期监控长时间计算的体积分数累积误差。
雷诺输运定理的先进研究
先进主题
RTT还有哪些先进应用？
RTT本身是古典的，但其扩展和应用仍在发展。
流固耦合（FSI）
在FSI中，流体与结构的界面进行动量与能量交换。用RTT计算流体侧的力，传递给结构分析侧。
分离耦合（Partitioned）：用不同的求解器分别求解流体和结构，交换界面数据。为了稳定，对界面力施加欠松弛或Aitken加速。
单体耦合（Monolithic）：流体和结构作为同一方程组求解。收敛性好但实现复杂。
商用工具中哪种比较多？
伴随法形状优化
伴随法用RTT定义的目的函数（阻力、压降等）的形状敏感度进行高效计算。
其中 $\mathbf{u}^\dagger, p^\dagger$ 是伴随变量。计算成本不依赖设计变量数，可进行数千个设计变量的形状优化。
声学类比
Ffowcs Williams-Hawkings音响类比是RTT的扩展。从移动控制面上的通量计算远场声压。
用于汽车风切声和喷气噪声预测。Fluent、STAR-CCM+都标配FW-H类比。
能量收支的损失分析
利用RTT的能量形式，可求得检查体积内熵生成率，定量化损失的空间分布。
这种方法被广泛用于涡轮机械效率改进。
RTT是现代CFD所有应用的基础，这一点真清楚了。
雷诺输运定理在太空飞行器燃料管理中的应用
在宇宙飞行器的轨道控制中，无重力环境下燃料箱内液体推进剂如何运动是生死攸关的问题。将雷诺输运定理扩展到微重力环境下进行的"晃动分析"对人造卫星姿态控制设计至关重要。NASA和JAXA会精密跟踪燃料在何时何处，否则推进器会吸入气体，推力瞬间变零。在任意检查体积中整理质量、动量、能量——这个定理是先端宇宙工学的根基。
雷诺输运定理的故障排除
故障排除
RTT相关的实务故障有哪些？
通量收支不匹配、动态网格守恒性、力计算误差是主要问题。
1. 入出口质量收支不符
症状：入口和出口的质量流量之差超过总流量的1%。
检查步骤：
1. 用Flux Report输出各边界的质量流量
2. 确认内部面（internal face）的截面流量
3. 检查残差的收敛水平
对策：收敛条件改为$10^{-5}$以下。若仍未解决，检查网格的非正交性。
2. 动态网格中质量泄漏
活塞运动分析中，质量逐步增加……
原因：几何守恒律（GCL）未满足。
检查方法：用网格移动但速度与压力均匀（$\mathbf{u} = 0$, $p = \text{const}$）的静止流体问题进行计算，验证速度和压力是否保持均匀。
对策：
- Fluent中使用In-Cylinder模型自动满足GCL
- OpenFOAM中设置
correctPhi yes; - 减小时间步长
3. 力计算结果异常
症状：升力系数是文献值的2倍以上。
| 检查项目 | 检查方法 |
|---|---|
| 参考面积是否正确 | $C_L = F_L / (0.5\rho U^2 A_{\text{ref}})$ 中的 $A_{\text{ref}}$ |
| 力的方向是否正确 | 升力垂直于流向，阻力平行于流向 |
| 压力基准值 | Operating Pressure设置是否正确 |
| 定常解是否收敛 | 力的监测数据是否振动 |
4. 滑移网格界面处不连续
症状：旋转体与静止部分的界面处速度或压力出现跳跃。
原因：
- 界面设置错误（配对错误）
- 界面网格不匹配
- 补间方法不当
对策：
- Fluent中正确定义Interface > Mesh Interface
- 旋转域和静止域的边界网格密度保持一致
- 确认GGI（General Grid Interface）设置
意识到RTT的守恒原理，可以及早发现问题。
没错。"进入的东西必然流出"、"力等于动量变化"——用这样的基本原则来验证结果，养成习惯。
质量守恒错误——CFD中"出口流量与入口不符"意味着什么？
在CFD中，入口质量流量与出口质量流量相差超过1%时，说明有问题。常见原因：①用不可压缩求解器处理可压缩流（密度=常数假设破裂）、②对称边界条件误设为墙面导致质量"泄漏"、③周期边界条件的压力差设置错导致逆流、④网格存在孔洞，物理上是非闭合区域。其中④经常由CAD导入时曲面法向不一致引发，常出现在"Negative Volume Cell"错误附近。收敛后必须输出全边界质量流量检查报告，验证入出平衡——这是CFD质量管理的基本功。
相关主题
更详细
错误