热交换器CFD分析

手工计算法(LMTD、ε-NTU)对整体性能估算有效，但管阵列偏流、挡板周边的死水域、温度分布的不均匀性等局部现象只有CFD才能评估。

需要联合求解流动和传热。Navier-Stokes + 能量方程。

整体性能评估的基本式是LMTD法和ε-NTU法。

对。从CFD求得的入出口温度计算 $Q$，再从中反推 $U$ 与设计值比较。

差异很大。壳管式的管外流动很复杂，挡板的影响很大。板式则是通道内的流动和传热占主导。

数百到数千根管一根根网格化是现实不可行的。实务中采用将管束作为多孔质介质+分布阻力(Distributed Resistance)来建模。

管束的压力损失用与管排列相应的相关式表示。正方排列的情况：

广泛使用Zukauskas(1972)的相关式。将其转换为多孔质介质的阻力参数输入到CFD中。

根据建模方法分别说明。

双管或管数在20根以下的小型热交换器可以实现全几何建模。将管壁作为固体区域网格化，用Conjugate Heat Transfer(CHT)联合流体和固体。

设置步骤：

使用Fluent的Heat Exchanger模型(Macro Model)。壳侧流动用3D CFD求解，管内侧用1D温度分布处理。管外侧的局部h值由子模型自动计算后求得传热量。

所需的输入数据：

Fluent的Heat Exchanger Model在内部基于Bell-Delaware法计算修正系数。使用者只需指定挡板的几何形状即可。

除了残差外，还需确认以下物理量已趋于稳定。

正是。能量平衡的确认是热交换器CFD最重要的检查项目。

壳侧从喷嘴导入流体时，会在挡板间发生偏流，部分管流量集中。通过评估这种偏流，优化喷嘴配置和挡板设计。

常见问题：

对。低流速区域颗粒容易沉降堆积，高流速区则易发生侵蚀。用CFD可视化流速分布来评估两种风险。

对鱼骨纹角度(人字形图案)差异的传热和压损权衡进行评估。以1通道周期模型进行参数研究。

常用的是Goodness Factor($j/f^{1/3}$、$j$是Colburn j-factor)。CFD计算各人字纹角度的j/f比进行对比，导出最优设计。

空调用翅片管的翅片表面温度分布和传热量用CFD评估。通过CHT求解翅片内部温度分布，得出翅片效率。

CFD不是用上式的理论值，而是从翅片表面温度分布直接计算翅片效率。特别是通道间距很小或翅片形状复杂时，理论式不适用，CFD必需。

除了3D CFD，热交换器专用设计工具也介绍一下。

HTRI是基于Bell-Delaware法扩展的1D/准2D设计计算工具，不是3D CFD。但基于实验数据库的相关式很充实，设计阶段的性能预测往往比HTRI可信性更高。

使用分工是这样的。

在Cell Zone条件中有Heat Exchanger Model，输入以下参数。

该模型对壳侧流动用3D CFD求解，管内侧作为1D温度分布处理。计算管外侧局部h值后自动求得传热量。

对。管内有复杂插件(涡流器等)或管内二相流时需要完整几何的CHT。

主要趋势为以下几个。

通过3D打印，可实现了传统上不可能的复杂流路形状(TPMS: Triply Periodic Minimal Surface、陀螺体结构等)。用CFD评估、优化这些新形状的传热性能的研究很活跃。

自支持周期性曲面结构，表面积/体积比很大。用式表示为：

此结构传热面积大，流路蛇形布置促进传热效应高。CFD中用OpenFOAM或COMSOL直接网格化STL数据。

用VOF(Volume of Fluid)法或欧拉多相模型求解气液二相。使用Lee模型等相变化模型来建模沸腾、凝聚。

正是。Lee模型的 $r_l$ 是经验值($0.1$～$10^6$ s⁻¹)，值的不同结果差异很大。必须与实验数据比较。

以传热量最大化、压损最小化为目标函数，通过拓扑优化自动设计流路形状的研究在进展。

目标函数例：

对。用Adjoint法高效计算感度，用Level Set法或Density法更新形状。

仍在研究阶段，但对离散设计变量(挡板数、截切方向)的优化已证明强化学习有效。

检查点：

对。Pr数大于1时(水: Pr≒7、油: Pr≒100～1000)，温度边界层比速度边界层薄，壁面附近网格分解不足会传热量估低。

症状: 高温侧放热量与低温侧受热量差5%以上。

对策：

对策：

症状: 管出口处管中心与壁面温度倒转。

原因: 管内壁面网格不足。在管截面配置最低15～20个网格。用O-grid网格保证壁面附近分解。

对策：

对。气泡的生成、成长、脱离精确建模目前仍是研究课题。工程上用RPI(Rensselaer Polytechnic Institute)模型的欧拉二相流模型较为实用。