强制对流 — CAE术语解说

具体来说，风扇或泵，以及物体本身的运动是典型的外力。例如，在服务器机箱的冷却中，直径120mm的DC风扇以每分钟1500转的速度送风，产生强制对流。与自然对流的决定性区别是Reynolds数变大，容易形成乱流。

数量级不同。自然对流的热传递率最多5~25 W/(m²·K)，而空气的强制对流可达25~250 W/(m²·K)，水则可达500~10,000 W/(m²·K)。这是因为乱流使热边界层保持很薄，热传递更高效。代表性的关联式Dittus-Boelert公式为，

是对流项。特别是能量方程中的

典型的设置有两种。一种是「速度入口」，根据风扇性能曲线得到的均匀流速，例如2.5 m/s。另一种是「质量流量入口」，例如指定冷却空气为0.05 kg/s。实务中也常直接内置风扇模型，Ansys Fluent的「Fan Boundary Condition」可输入压力-流量曲线。

如果要准确捕捉壁面附近的热传递，SST k-ω模型是标准选择。它与壁面函数结合，能在不求解粘性底层的情况下估计热流。但对于内部流动，如果没有复杂的分离，实际上Standard k-ε模型加「Enhanced Wall Treatment」也很常用。无论如何，壁面的y+值应保持在1左右，或者应用壁面函数时保持在30以上，这很重要。

对于强制对流，当对流占主导且流道形状相对简单时，SIMPLE系列算法（SIMPLEC或PISO）就足够了。计算稳定性好，内存消耗少。但对于复杂的三维形状，流向变化急剧的情况，Ansys Fluent的「Coupled」求解器可能收敛更快。代价是内存使用增加1.5~2倍。

首先是Reynolds数的概算。例如，流道水力直径D_h=0.01m、流速u=1m/s、空气动力粘度ν=1.5e-5 m²/s，则，

必须在壁面垂直方向贴边界层网格。乱流情况下，根据y+值的目标确定第一层厚度。以y+≈1为目标时，第一层厚度Δy按照

即使没有实验数据，也可以做以下检查。1) 质量守恒：入口和出口的质量流量差在0.5%以内吗。2) 热平衡：固体向流体传递的热量与流体的焓增是否一致（允许差5%以内）。3) 无量纲数对照：圆管内流，从分析结果求出的Nusselt数与Gnielinski关联式等文献值比较。这些都过关的话，一定程度上可以认为有信度。

核心物理模型没有很大差异，但工作流程和前处理差别很大。Fluent是「几何→网格（Ansys Meshing/Fluent Meshing）→设置→求解」的分离式工作流。而Star-CCM+是在单一GUI内完成所有工序的一体化方式。强制对流中常见的「参数化研究」，比如风扇转速从500rpm变到2000rpm，Star-CCM+的「Field Function」和「Table」配置直观高效。

这类「共轭传热」问题中，各软件的「简化模型」有无很关键。例如，Fluent的「PCB模型」把基板均质化成块，可设定发热量和等效导热系数。Star-CCM+也有类似的「Thin Conductor」模型。另一方面，COMSOL Multiphysics在「热阻网络」直接建模方面很强，能考虑IC芯片包JEDEC规范的热阻（Θ_ja等），这方面很有优势。

OpenFOAM有强制对流基础方程的求解功能。标准求解器「buoyantSimpleFoam」或「chtMultiRegionFoam」可用。但商业软件那样高级的风扇模型或GUI界面的直观边界条件设置没有。全部要用文本形式编辑「fvSchemes」「fvSolution」等配置文件。乱流模型选择多，但这意味着用户需要深入理解详细内容。大规模并行计算是它的优势。

最常见原因是不恰当的初始条件和求解器方案。强制对流中，流速越高对流项越强，越容易发散。对策按这个顺序试：1) 不要用均匀值（流速0、压力0）初始化所有变量，而要用接近入口值的值。2) 降低求解器的「欠松弛系数」。特别是动量和压力降到0.3~0.5。3) 从「一阶迎风」等稳定的一阶精度方案开始，收敛后切换到「二阶迎风」。

几乎肯定是网格，特别是壁面的y+值。y+值在1~5范围（解析粘性底层）却用了「壁面函数」，或反过来y+>30还用「Enhanced Wall Treatment」，这样剪切应力与热流的关系就不能正确计算。首先要可视化壁面y+分布，确认与使用的乱流壁面处理前提条件是否一致。别忘了在Ansys Fluent的「Report → Reference Values」中设定正确的代表长度和速度。

出口反向流是出口边界离流向太近，或出口静压设置过高引起的。对策有两个。第一，如果可能的话延长流道，把出口设在完全展开流的位置。第二，在Fluent的「Pressure Outlet」设置中，把「Backflow Direction Specification」从「Normal to Boundary」改成「From Neighboring Cell」，设置反向流时更现实的条件（比如入口温度和乱流强度）。根本上，改用「Mass Flow Outlet」固定质量流量最有保障。

应该连续的界面温度出现不连续，是网格不连续或界面边界条件设置错误。首先检查固体和流体网格在界面是否一致（conforming mesh）。如果是不一致的非conformal网格，在Ansys Fluent中要正确定义「Interface」对，检查「Mesh Interface」设置中是否选了「Couple」。COMSOL会自动生成「Identity Pair」，但热流连续性的维持可能需要调整「对的扩展」。