热交换器CFD分析
理论与物理
概述
老师！进行换热器CFD分析的目的是什么？用NTU法或LMTD法手工计算不够吗？
手工计算法（LMTD、ε-NTU）对于整体性能估算有效，但管排偏流、折流板周围的死区、温度分布不均匀等局部现象，不用CFD则无法评估。
控制方程
请告诉我换热器CFD的基本方程。
流动与传热耦合求解。是Navier-Stokes方程 + 能量方程。
整体性能评估的基本公式是LMTD法和ε-NTU法。
$U$ 是总传热系数，$A$ 是传热面积对吧。CFD中是从计算结果反算这些值吗？
是的。根据CFD求得的进出口温度计算出 $Q$，然后反算出 $U$ 与设计值进行比较。
按换热器类型的建模
管壳式和板式换热器的分析方法不同吗？
有很大不同。管壳式换热器的壳程流动复杂，折流板影响大。板式换热器则主要是通道内的流动与传热。
| 类型 | CFD模型 | 主要评估项目 |
|---|---|---|
| 管壳式 | 管束 + 折流板 | 壳程偏流、死区 |
| 板式 | 单通道周期模型 | 波纹板的传热强化效果 |
| 翅片管（空调） | 单排周期模型 | 翅片效率、结露 |
| 套管式 | 全3D模型 | 流动模式（逆流/顺流） |
管壳式换热器的多孔介质方法
管壳式换热器的管束全部建模很麻烦吧？
对数百到数千根管子分别进行网格划分是不现实的。实际工作中，将管束建模为多孔介质 + 分布阻力。
管束的压降用根据管排布置的相关式表示。正方形排列时：
$N_r$ 是管排数，$\chi$ 是排列修正系数，$f$ 是摩擦系数对吧。$V_{max}$ 是最小截面处的流速。
广泛使用的是Zukauskas (1972) 的相关式。将其转换为多孔介质的阻力参数输入到CFD中。
实际注意事项
换热器CFD需要特别注意什么？
换热器理论的确立——LMTD法与ε-NTU法（1940年代）的诞生
作为换热器热设计理论至今仍在使用的“LMTD（对数平均温差）法”和“ε-NTU（效率-传热单元数）法”，是由1940年代的Mason(1954)和Kays & London(1964)等人系统建立的。特别是Kays & London的名著《Compact Heat Exchangers》，整理了数百种散热片/翅片形状的St（斯坦顿数）和f（摩擦系数）的实验数据，成为紧凑型换热器设计的圣经。这些实验关联式在50年后的现代仍被用于CFD结果的合理性检查。如果CFD计算出的Nu和f与Kays-London关联式偏差超过20%，通常是网格或边界条件存在问题的信号。
各项的物理意义
- 时间项 $\partial(\rho\phi)/\partial t$：想象一下拧开水龙头的瞬间。最初水会不稳定地哗啦哗啦流出，过一会儿就变成稳定的水流了吧？描述这个“正在变化的过程”的就是时间项。心脏搏动导致血流脉动，发动机阀门每次开闭导致流动变化，这些都是非定常现象。那么定常分析是什么？就是只看“经过足够时间流动稳定之后”——也就是将此项设为零。计算成本大幅降低，因此先用定常求解是CFD的基本策略。
- 对流项 $\nabla \cdot (\rho \mathbf{u} \phi)$：把落叶丢进河里会怎样？会被水流带着运往下游对吧。这就是“对流”——流体运动搬运物质的效果。暖风的暖气能到达房间的另一端，也是因为空气这个“搬运工”通过对流输送热量。这里有趣的是——这项包含“速度×速度”，因此是非线性的。也就是说，流速变快这项会急剧增强，变得难以控制。这就是湍流的根本原因。常见的误解：“对流和传导差不多”→ 完全不一样！对流是流动搬运，传导是分子传递。效率有天壤之别。
- 扩散项 $\nabla \cdot (\Gamma \nabla \phi)$：有过在咖啡里倒入牛奶后放置的经历吗？即使不搅拌，过一会儿也会自然混合对吧。那就是分子扩散。那么下一个问题——蜂蜜和水，哪个更容易流动？当然是水对吧。因为蜂蜜的粘性（$\mu$）高，所以不易流动。粘性越大扩散项越强，流体的运动就变得“粘稠”。雷诺数小的流动（缓慢、粘稠）中扩散占主导。相反，Re数大的流动中对流占压倒性优势，扩散则成为配角。
- 压力项 $-\nabla p$：注射器的活塞一推，液体就从针头有力地喷出对吧？为什么呢？因为活塞侧压力高，针头侧压力低——这个压力差成为推动流体的力。水坝放水也是同样原理。天气图上等压线密集的地方会怎样？没错，会刮强风。“有压力差的地方就会产生流动”——这就是纳维-斯托克斯方程压力项的物理意义。这里的误解点：CFD的“压力”多为表压而非绝对压力。切换到可压缩分析时结果突然变得奇怪，原因可能就是混淆了绝对压力/表压。
- 源项 $S_\phi$：被加热的空气会上升——为什么呢？因为比周围空气轻（密度低），被浮力推上去了。这个浮力作为源项添加到方程中。此外，燃气灶火焰产生化学反应热、工厂电磁泵对金属熔液施加洛伦兹力……这些都是“从外部向流体注入能量或力”的作用，用源项表示。忘记源项会怎样？自然对流分析中忘记加入浮力，流体就完全不动——冬天房间里开了暖气但热空气不上升，得到这种物理上不可能的结果。
假设条件与适用范围
- 连续介质假设：克努森数 Kn < 0.01（分子平均自由程 ≪ 特征长度）时成立
- 牛顿流体假设：剪切应力与应变速率呈线性关系（非牛顿流体需要粘度模型）
- 不可压缩假设（Ma < 0.3时）：将密度视为常数处理。马赫数0.3以上需考虑可压缩性效应
- Boussinesq近似（自然对流）：仅在浮力项中考虑密度变化，其他项使用恒定密度
- 不适用的情形：稀薄气体（Kn > 0.1）、超音速/高超音速流动（需要捕捉激波）、自由表面流动（需要VOF/Level Set等方法）
量纲分析与单位制
| 变量 | SI单位 | 注意事项·换算备忘 |
|---|---|---|
| 速度 $u$ | m/s | 入口条件从体积流量换算时，注意截面面积的单位 |
| 压力 $p$ | Pa | 区分表压与绝对压力。可压缩分析使用绝对压力 |
| 密度 $\rho$ | kg/m³ | 空气: 约1.225 kg/m³@20°C、水: 约998 kg/m³@20°C |
| 粘性系数 $\mu$ | Pa·s | 注意与运动粘性系数 $\nu = \mu/\rho$ [m²/s] 的混淆 |
| 雷诺数 $Re$ | 无量纲 | $Re = \rho u L / \mu$。层流/湍流转换的判断指标 |
| CFL数 | 无量纲 | $CFL = u \Delta t / \Delta x$。直接关系到时间步长的稳定性 |
数值解法与实现
数值方法的细节
请告诉我换热器CFD的具体设置方法。
我来解释根据建模方法的不同设置。
方法1: 全几何模型（小规模·详细评估）
是对每根管子进行网格划分的方法吧。
对于套管式或管子数量在20根以下的小型换热器，可以采用全几何模型。将管壁作为固体区域进行网格划分，通过共轭传热耦合流体与固体。
设置步骤：
1. 分别定义流体域（管内、管外）和固体域（管壁）
2. 在流体-固体界面设置耦合壁面条件
3. 湍流模型: Realizable k-ε + 增强壁面处理（管外）
4. 能量方程: 开启
5. 物性值: 温度相关（特别是粘度）
| 参数 | 管内侧 | 管外侧（壳程） |
|---|---|---|
| 湍流模型 | Realizable k-ε | SST k-ω（存在折流板分离时） |
| y+ | ≈ 1 | ≈ 1〜30 |
| 网格 | O型网格 + 棱柱层 | 四面体 + 棱柱层 |
方法2: 多孔介质模型（大规模·整体评估）
是将管束视为多孔体的方法吧。
使用Fluent的Heat Exchanger模型（宏观模型）。将壳程视为多孔介质，通过子模型计算管束的阻力和传热。
需要的输入数据：
- 管外径、管间距、管排列模式（正方形、三角形）
- 管内流量和入口温度
- 管外侧的传热关联式（Zukauskas, Kern, Bell-Delaware等）
- 折流板数量、间距、切割率
可以直接输入Bell-Delaware法的参数吗？
Fluent的Heat Exchanger Model内部会计算基于Bell-Delaware法的修正系数。用户只需指定折流板的几何形状即可。
共轭传热的设置
CHT计算需要注意什么？
| 材料 | 热导率 [W/(m K)] | 固体网格 |
|---|---|---|
| 铜 | 385 | 可以粗（3层） |
| 铝 | 205 | 可以粗（3层） |
| 不锈钢 | 16 | 稍细（5〜8层） |
| 钛 | 22 | 稍细（5〜8层） |
收敛判定
换热器CFD的收敛判定看什么好？
除了残差之外，以
この記事の評価ご回答ありがとうございます！参考に
なったもっと
詳しく誤りを
報告
换热器CFD的收敛判定看什么好？
除了残差之外，以
なった
詳しく
報告