空气贴近高温壁面后吸热,温度升高,密度随之降低。在重力场中,这部分较轻的空气受到浮升作用而向上运动。周围温度较低、密度较大的空气随后补入,空间内逐渐形成连续的循环流动。这就是自然对流的基本过程。
因此,自然对流至少需要三个条件:流体内部存在温度差,密度能够随温度变化,同时计算域处于非零重力场中。能量方程负责计算温度变化,却不会单独产生浮升力。如果空气仍采用恒定密度,温度即使已经升高,动量方程也无法识别冷热空气之间的密度差。
在封闭腔体中,热流体上升后会沿顶部流向较冷区域,再沿冷壁下降,最终回到底部,形成内部循环。开放空间还会伴随空气从开口进入和排出,流动路径同时受开口位置、压力边界及外部环境影响。

图1 热壁面附近空气上升,冷壁面附近空气下降,形成封闭循环(原理示意)
Fluent中的重力加速度默认值为零。建立自然对流模型时,需要先确认几何模型的实际竖直方向,再把重力分量设置到对应坐标轴。例如,模型的竖直方向为Y轴向上,重力通常应设置在Y轴负方向,而不是凭界面习惯直接填写某个分量。
重力方向设置反了,浮升方向也会随之反转。发热元件上方本应出现的热羽流可能转向下方,冷空气补充路径也会颠倒。求解器并不知道哪个方向才是现实中的上下,只会按照输入的坐标和方程计算,所以错误模型仍可能得到平滑云图和较低残差。
检查时不要只看重力数值,还要在前处理或后处理中确认坐标轴。模型经过旋转、镜像或从其他软件导入后,竖直方向尤其容易与最初设想不一致。

图2 重力方向反转会使计算得到的浮升方向随之反转(原理示意)
恒定密度适合不需要考虑热浮升的不可压缩流动。用于自然对流时,恒定密度会切断温度与浮升力之间的联系。此时即使发热体附近形成明显温升,也不会因为这部分空气变轻而产生正确的有组织上升流。
Boussinesq近似是一种常用处理方法。它把密度在大部分控制方程中视为常数,只在重力产生的浮升项中考虑密度随温度的线性变化。这样既保留了自然对流的主要驱动力,又通常比完整的变密度模型更容易收敛。
这项近似要求实际密度变化较小,判断依据应回到流体的热膨胀系数、参考温度和计算温差,而不是套用一个适用于所有问题的温差界限。温度范围较宽或物性呈明显非线性时,可以考虑不可压缩理想气体、温度函数或其他合适的密度关系,并同步检查工作压力与参考密度的处理。

图3 恒定密度无法建立温度差驱动的有组织浮升流(原理示意)
采用Boussinesq近似时,参考温度是密度线性化的基准。参考温度应能代表计算域内的特征温度,热膨胀系数和参考密度也应与该温度附近的物性相匹配。参考值选得不合理,会改变浮升项的计算基础。
边界条件决定热量从哪里进入、空气能否交换以及热量最终从哪里离开。封闭腔体依靠内部循环传热,系统总质量保持不变;带开口的机柜或房间则可能通过压力入口、压力出口等边界与外界交换空气。两类模型的物理含义不同,不能只因为几何上都有一个孔,就使用相同的处理方式。
外壁是否绝热、是否给定对流换热系数、是否考虑辐射,也会影响温度水平和循环强度。自然对流往往速度较低,辐射换热在高温表面附近可能占据不可忽略的比例,需要根据温度和表面性质判断。

图4 自然对流模型的五项设置与结果检查顺序(原理示意)
第一步是看总体方向。高温流体应在重力反方向产生主要浮升趋势,低温流体从周围补入。这里说的是整体趋势,局部几何阻挡、分离和回流仍可能改变某一小片区域的流向,不能把所有局部回流都判为错误。
第二步是看循环路径和速度量级。速度很小并不自动等于死区,自然对流本身可能就是低速流动。应结合温差、特征长度、流体物性以及瑞利数判断流动强弱,并观察速度监测量和温度监测量是否稳定。某些高瑞利数流动可能具有明显非定常特征,稳态求解未必能够给出稳定解。
最后检查质量守恒和能量平衡。开放系统需要核对各开口的质量流量是否闭合;传热计算需要比较热源输入、壁面散热和流体带走的热量。云图颜色符合直觉,只能作为初步检查,不能代替守恒关系和关键监测量。
提交计算前,可以依次复核以下五项:
① 确认几何坐标和实际竖直方向,检查重力分量与单位;
② 确认密度能够随温度影响浮升项,并核对所选模型的适用范围;
③ Boussinesq模型下检查参考温度、参考密度和热膨胀系数;
④ 区分封闭空间与开放空间,检查热边界、开口边界和必要的辐射换热;
⑤ 检查总体流向、循环路径、监测量、质量守恒和能量平衡。
自然对流模型最常见的问题并不是某个求解按钮没有打开,而是物理条件之间没有形成闭环。把重力、密度、参考状态和边界条件逐项对应起来,结果才具有可解释性。
如果模型已经收敛,但热羽流方向、温度水平或开口流量仍不合理,可以结合具体几何、物性和边界条件进一步排查。需要时也可提供模型排错、结果分析和案例复现辅导。
• Ansys Fluent User's Guide Operating Conditions https://ansyshelp.ansys.com/public/Views/Secured/corp/v251/en/flu_ug/flu_ug_bcs_sec_operating.html
• Ansys Fluent User's Guide Density https://ansyshelp.ansys.com/public/Views/Secured/corp/v251/en/flu_ug/flu_ug_sec_density.html
• Ansys Innovation Space How to Model Natural Convection in Ansys Fluent
https://innovationspace.ansys.com/courses/courses/heat-transfer-modeling-in-ansys-fluent/lessons/how-to-model-natural-convection-in-ansys-fluent-lesson-4