流体的能量方程
流体能量方程的理论基础
概述
老师，流体的能量方程在什么情况下需要用到呢？
当需要求解温度场时必不可少。热交换器的设计、电子设备冷却、燃烧、自然对流等都需要求解能量方程。等温流动仅需NS方程和连续方程就能闭合，但涉及温度时需要加入能量方程。
能量方程（温度形式）
不可压缩流动的温度场能量方程如下所示。
其中 $c_p$ 是定压比热，$k$ 是热导率，$\Phi$ 是粘性耗散函数，$\dot{q}$ 是内部发热率。
粘性耗散是什么意思？
粘性耗散是指由于粘性作用，动能转化为热的项。对于不可压缩流动，
通常在普通工程流动中粘性耗散极为微小可以忽略，但在高粘流体（聚合物熔融等）或超高速流中就不能忽视。
埃克特数
粘性耗散的重要性用埃克特数（Eckert number）来评估。
当 $Ec \ll 1$ 时粘性耗散可以忽略。例如空气流（$U = 50$ m/s, $\Delta T = 20$ K）时 $Ec \approx 0.12$ 较小。而聚合物挤出（$U = 0.1$ m/s, $\mu = 1000$ Pa·s）中耗散成为温度上升的主要因素。
能量方程（焓形式）
可压缩流动中如何处理？
使用全焓 $h_0 = h + \frac{1}{2}|\mathbf{u}|^2$ 的能量方程。
可压缩流动中加入状态方程 $p = \rho R T$（理想气体），方程组才能闭合。密度、速度、压力、温度四个场完全耦合。
无量纲参数
| 参数 | 定义 | 物理意义 |
|---|---|---|
| 普朗特数 Pr | $\nu/\alpha = \mu c_p/k$ | 动量扩散与热扩散的比值 |
| 努塞尔特数 Nu | $hL/k$ | 对流/传导的传热比 |
| 佩克莱特数 Pe | $Re \cdot Pr$ | 对流/扩散的比值 |
| 埃克特数 Ec | $U^2/(c_p\Delta T)$ | 动能/热能 |
Pr数由流体本身性质确定，而Nu是计算结果呢。
完全正确。空气 $Pr \approx 0.71$，水 $Pr \approx 7$，发动机油 $Pr \approx 100$～1000。Pr数越大，热边界层越薄，为精确解析壁面附近的温度梯度需要更细的网格。
能量方程的历史——从焦耳的热功当量实验（1843年）到流体力学的应用
流体能量方程中包含的热功当量概念源自英国科学家詹姆斯·焦耳（James Joule）1843年的搅拌实验。焦耳常数1cal=4.186J至今仍是热力学基本常数。将此纳入流体力学是通过傅里叶热传导方程（1822年）与纳维-斯托克斯方程的结合，能量方程的标准形式在1850年代确立。特别是"粘性耗散项"——速度梯度转化为热的效应——在超音速流和高粘流体中不可忽视，现代CFD在高精度分析时必须启用这一项。低速低粘流体的工程CFD中常被省略，但其适用范围（Br=ηU²/(kΔT)>0.1）必须认识清楚。
流体能量方程的数值计算方法
能量方程的离散化
能量方程怎样进行数值求解呢？
用有限体积法时，对流项和扩散项作为单元面上的通量进行离散化。
其特点是能量方程是标量输运方程，只要速度场已知，就可以作为线性问题求解。可以先解NS方程，然后作为后处理单独求解能量方程。
对流格式的选择
温度场的对流格式能否和速度场相同？
取决于佩克莱特数（Pe数）。Pe数大时对流主导，中心差分会产生数值振荡。一般规则如下。
| Pe数范围 | 推荐格式 | 备注 |
|---|---|---|
| Pe < 2 | 中心差分（CD） | 扩散主导，稳定 |
| Pe > 2 | 迎风差分（Upwind） | 有数值扩散 |
| 高精度 | QUICK、TVD（MUSCL等） | 兼顾精度和稳定性 |
热边界条件
壁面的热边界条件主要有三种。
| 边界条件 | 数学表达 | 应用场景 |
|---|---|---|
| 固定温度（第一类） | $T_{wall} = T_0$ | 冷却水壁面、恒温槽 |
| 固定热流（第二类） | $-k\frac{\partial T}{\partial n} = q_w$ | 加热器、发热面 |
| 对流传热（第三类） | $-k\frac{\partial T}{\partial n} = h(T - T_\infty)$ | 外界环境热交换 |
湍流中的温度场
湍流情况下温度场也需要模型吗？
RANS中需要对湍流热通量 $\overline{u_i'T'}$ 进行建模。最常用的是涡粘数模型。
$Pr_t$（湍流普朗特数）通常设为0.85～0.9。壁面附近的温度分布（如Jayatilleke壁函数等）也很重要。
湍流Pr数对结果影响有多大？
对壁面传热系数约有10～20%的影响。特别是液态金属（$Pr \ll 1$，$Pr_t \approx 1$～4）中，用标准的0.85不精确，需要专门模型。
能量方程的离散化——"隐式还是显式"由温度的变化速率决定
在数值求解流体能量方程时，显式还是隐式的选择很大程度影响计算效率。显式法易于实现但受热扩散稳定条件约束 Δt≤ρcΔx²/(2k)。在混有绝热材料和金属的模型中，热导率可能相差100倍，稳定的Δt受最薄金属网格限制，整体计算时间大幅增长。此时转向隐式法，Δt可大幅增加。"方程形式相同"却因离散化策略不同导致计算时间相差数十倍——这是数值解法的现实。
流体能量方程的实际应用
流体能量方程的实际应用
共轭传热（CHT）分析最近听得比较多，是什么意思？
这是同时求解流体对流和固体热传导的方法。电子基板的冷却设计、发动机缸体、燃气轮机叶片冷却等，流体-固体界面的温度分布很重要，必须用这个方法。
共轭传热分析的要点
流体-固体界面处需满足以下连续条件。
| 求解器 | CHT方法 | 设置方式 |
|---|---|---|
| Fluent | Coupled Wall | 用Coupled Wall连接固体和流体区域 |
| STAR-CCM+ | Multi-Region | 区域间用interface连接 |
| OpenFOAM | chtMultiRegionFoam | 各区域切换对应求解器 |
Nusselt数的计算
怎样从CFD结果计算努塞尔特数？
从壁面传热系数 $h$ 计算。
$T_{\text{ref}}$ 根据问题选择入口温度或平均温度。圆管内流可与Dittus-Boelter式 $Nu = 0.023 Re^{0.8} Pr^{0.4}$ 比较，这是标准做法。
热分析的网格要求
温度场的精度直接取决于网格质量。
| 分析对象 | 额外要求 | 原因 |
|---|---|---|
| 自然对流 | 壁面附近网格密集 | 浮力驱动的温度边界层 |
| 高Pr流体 | $y^+ \approx 1$ | 热边界层比速度边界层薄 |
| CHT | 固液界面网格整合 | 界面温度插值精度 |
| 辐射 | 视角分辨率 | 取决于S2S、DOM等模型 |
高Pr流体时，速度场用壁函数够了，温度场却需要Low-Re是吗？
正是如此。水（$Pr \approx 7$）的热边界层厚度约为速度边界层的 $Pr^{-1/3} \approx 0.52$ 倍。发动机油（$Pr \approx 500$）更薄，所以壁面第一层网格的分辨率极其重要。
自然对流的注意事项
Boussinesq近似将密度变化仅在浮力项中反映。
瑞利数 $Ra = Gr \cdot Pr = g\beta\Delta T L^3/(\nu\alpha)$ 超过 $10^8$ 时开始湍流转变。$\beta\Delta T < 0.1$ 是Boussinesq近似的判据。
燃气轮机燃烧室的全焓守恒——CFD入口条件设置的实际应用
燃气轮机燃烧室的CFD分析中，压气机出口的高压高温气体作为入口边界条件需要精确设置。能量方程的实施中"静温度（Static Temperature）"与"全温度（Total Temperature）"的混淆是最常见的设置错误，高速流（Mach>0.3）中二者差异可达几十℃。燃料的化学焓值（天然气Lower Heating Value: 43.1MJ/kg）作为热源项的正确处理也至关重要，否则燃烧后温度会被低估500～800K。在实际设计中，出口条件与入口焓守恒的质量-能量平衡验证，误差控制在1%以内是国内航空发动机制造厂的标准CFD验证步骤。
流体能量方程的软件对比
能量方程相关的商用工具对比
热分析功能在各求解器间是否有差异？
基本能量方程各求解器都能处理，但辐射模型、CHT、相变的差异较大。
辐射模型的对比
| 模型 | Fluent | CFX | STAR-CCM+ | OpenFOAM |
|---|---|---|---|---|
| S2S（面到面） | 支持 | 支持 | 支持 | viewFactor模型 |
| DOM（离散坐标法） | 支持 | 支持 | 支持 | fvDOM |
| P-1（球谐近似） | 支持 | 支持 | 支持 | P1 |
| 蒙特卡洛 | Fluent 2024之后 | 不支持 | 支持 | 非标准 |
选择哪种辐射模型比较好？
用光学厚度 $\tau = \kappa L$ 判断。
- $\tau > 3$（光学厚）：P-1足够
- $\tau < 0.1$（光学薄）：S2S最高效
- 其他情况：DOM通用
相变与沸腾模型
| 功能 | Fluent | STAR-CCM+ | OpenFOAM |
|---|---|---|---|
| 凝固融化 | 焓-孔隙度法 | 类似 | solidificationMelting |
| 蒸发凝结 | Lee模型、VOF | VOF+相变 | 自定义UDF |
| 沸腾 | RPI沸腾模型 | RPI沸腾模型 | 研究代码 |
各求解器的设置技巧
记住各求解器启用能量方程的步骤。
- Fluent：Models > Energy 打开。CHT时在Cell Zone Conditions > Solid设置固体材料
- CFX：Domain > Heat Transfer > Thermal Energy或Total Energy
- STAR-CCM+：Physics > Energy 打开。Multi-Region CHT在Region Interface设置
- OpenFOAM：选择buoyantSimpleFoam（有浮力定常）或chtMultiRegionFoam（CHT）
Fluent中打开Energy模型计算时间增加多少？
温度作为一个标量方程加入，计算时间仅增加10～20%。加辐射模型再增20～50%。CHT情况（固体区域增加）则根据单元数增加，但固体无对流项，计算比流体轻。
能量方程实施的工具差异——Fluent与OpenFOAM壁面热流计算的区别
不同CFD求解器对能量方程的离散化和壁面热流计算方法存在细微差异。ANSYS Fluent计算壁面热流qw=h_wall × (T_wall - T_fluid_ref)时，T_fluid_ref的定义（壁面单元中心值 vs 对数律修正值）可通过模型设置切换。而OpenFOAM的buoyantSimpleFoam则在壁面边界条件中用fixedHeatFlux或externalWallHeatFluxTemperature区别处理，设置错误会导致局部热流量偏差超过2倍。特别是壁函数（Wall Function）与低Re处理混合使用时，自动化工具与手工设置的结果会产生较大偏差。为确保可靠的热流预测，必须在简单验证算例（平板强制对流等）上先检验该工具的默认设置。
流体能量方程的前沿研究
前沿话题
能量方程还有哪些前沿研究课题？
温度场建模的深化在产业和研究中都非常活跃。
湍流热通量的高级模型
标准涡粘模型（$Pr_t = 0.85$ 固定）在以下情况不够。
- 液态金属冷却（$Pr \sim 0.01$）：$Pr_t$ 取决于 $Pr$
- 强逆压力梯度：湍流热通量与温度梯度非平行
- 喷流冷却：滞止点附近传热系数预测
高级模型包括代数热通量模型（AFM）和完整两方程热通量模型（$\overline{T'^2}$-$\varepsilon_\theta$ 模型）。
LES中的温度场
LES中温度的亚网格尺度(SGS)通量需要模型。
$Pr_{\text{sgs}} \approx 0.4$～0.6为典型值。动态SGS模型可动态计算 $Pr_{\text{sgs}}$。
超临界流体的传热
超临界CO2循环最近受关注，CFD能处理吗？
超临界状态中，在拟临界点附近物性值（$c_p$、$\rho$、$k$、$\mu$）急剧变化。需与NIST REFPROP或CoolProp物性库联接，使用温度-压力相关的物性表。
标准RANS湍流模型对超临界传热预测精度较低。壁面附近密度变化影响湍流结构。
燃烧与化学反应
燃烧中能量方程加入化学反应发热项。
$h_k^0$ 是物种 $k$ 的生成焓，$\dot{\omega}_k$ 是反应速率。Flamelet模型、EDC（涡耗散概念）、PDF输运方程等多种燃烧模型可用。
能量方程涉及的问题果然非常广泛呀。
不错。不涉及温度的工程问题反而罕见。因此深刻理解能量方程的物理本质至关重要。
数据中心的"液体浸没冷却"——流体能量方程再度成为主角
近年快速普及的液体浸没冷却（immersion cooling）中，服务器主板直接浸入绝缘油或冷却液。此时流体能量方程的对流项与热传导项的平衡决定了CPU温度。相比风冷，传热系数增加20～50倍，因此能量方程的系数在完全不同的量级。最前沿的两相液体浸没（相变冷却）中，蒸发凝结的潜热项加入，能量方程的形式大幅改变。流体能量方程的最前沿正与AI芯片冷却设计紧密相连。
流体能量方程的故障排除
流体能量方程的故障排除
温度场计算常见的问题有哪些？
温度发散、壁面传热系数不匹配，是常见问题。
1. 温度发散为非物理值
现象：温度变成 $10^{10}$ K等异常值。
| 原因 | 对策 |
|---|---|
| 能量残差未收敛 | 降低松弛因子（0.8→0.5） |
| 启用粘性耗散但速度场不准 | 先收敛速度场再启用能量 |
| 辐射模型吸收系数不合理 | 重新检查光学性质 |
| 温度相关物性表范围超出 | 扩展物性表的温度范围 |
2. 努塞尔特数与实验相关式不符
圆管内流的Nu与Dittus-Boelter式相差30%……
| 原因 | 诊断 | 对策 |
|---|---|---|
| y+不合适 | 检查壁面y+ | 热传达用 $y^+ < 5$ 推荐 |
| 湍流Pr数 | 检查是否默认0.85 | Pr>10液体需做敏感性分析 |
| 助流区 | 是否排除助流区 | 确保 $L/D > 10$ |
| 参考温度 | 检查 $T_{\text{ref}}$ 定义 | 用平均温度 |
Dittus-Boelter式本身有 $\pm 25$%的误差。CFD与其吻合在 $\pm 10$%以内算是良好。
3. 自然对流不发生
现象：有加热壁面却完全没有流动。
原因：
- 重力未设定（$g = 0$）
- Boussinesq近似未启用（密度保持常数）
- Operating Density设置不对
对策：Fluent中设置Operating Conditions > Gravity，材料物性中Density选Boussinesq，Operating Density设为参考温度下的密度。
4. 共轭传热界面温度不连续
固体和流体的边界处温度出现跳跃……
原因：界面网格未对齐（非conformal），或interface条件设置错误。
对策：
- 流体和固体网格界面对齐（conformal mesh）
- 非conformal情况检查interface的Mapped设置
- 确认壁面Thermal Condition设为Coupled
温度场的调试基础是速度场计算正确呢。
完全同意。温度场依附于速度场，必须先收敛速度和压力场，再调试能量方程。
"能量不守恒"的Bug曾是电力业的大问题
1970年代，核电站冷却系统的CFD分析出现"计算能量平衡偏离现实超过5%"的问题频繁发生。大多数原因是粘性耗散项的符号错误或对流项与扩散项离散化不一致。5%的误差会导致冷却材温度评估偏离，直接影响安全评估，当时的监管部门因此要求对计算代码进行验证。流体能量方程的故障排除经验告诉我们：理论错误远少于实现错误——这个教训至今仍然适用。
更详细
错误