对流单元模拟器 返回
热流体模拟器

对流单元模拟器 — 瑞利–贝纳尔对流

改变瑞利数(Ra),实时观察热对流单元的形成、发展和湍流转变

参数设置

Ra = 2512
温度场
速度矢量
流线叠加
计算结果
最高温度
最低温度
最大流速
对流状态
计算中...
努塞尔数 Nu
模拟
理论·主要公式

$$Ra = \frac{g \beta \Delta T L^3}{\nu \alpha}$$

瑞利数:浮力 vs 粘性·热扩散的比值。$Ra \gt Ra_c \approx 1708$ 时贝纳尔对流单元出现。

$$Nu = C \cdot Ra^n$$

努塞尔特数(热传达比):表示对流强度。$n \approx 1/3$ 对应高 $Ra$ 湍流对流。

$$\alpha = \frac{k}{\rho c_p}$$

热扩散率(m²/s):$k$ 为热导率,$\rho$ 为密度,$c_p$ 为比热。

什么是瑞利–贝纳尔对流

🙋
烧开水时,锅底会看到旋转流动,那是对流单元吗?
🎓
完全正确!正式名称是"瑞利–贝纳尔对流"。底面加热后,底部流体温度升高、密度下降,由于浮力上升。上面很冷,所以冷的流体下降。这种循环形成规则的"涡卷式对流单元"。法国物理学家贝纳尔在1900年用鲸油实验首次观察到这一现象,瑞利爵士在1916年进行了理论化分析。
🙋
但是,只是稍微加热一下,对流不会发生,对吧?什么条件才能触发?
🎓
问得好!这正是瑞利数(Ra)的概念。Ra表示"浮力不稳定化"对"粘性+热扩散稳定化"的比值。Ra超过临界值约1708时对流才开始。低于这个值,上升的流体会被粘性和热扩散压制,只能通过热传导散热。在模拟器中选择"稳定层"预设(Ra < 1708),会看到对流单元消失,温度梯度变得均匀。
🙋
那为什么一旦超过Ra,就会形成特定的"单元"模式,而不是随机混合?
🎓
这很有意思!不稳定化开始时,各种空间尺度的扰动都会略微增长。但存在"最优波长"——约为层厚度的2〜3倍——能增长最快。只有这个波长被保留下来。这是"模式形成"自组织的例子。提高Ra时,最初的漂亮涡卷会扭曲、分裂,多个尺度的运动混在一起,最终变成湍流。大气积云对流和地球地幔对流也遵循相同原理。

物理模型和支配方程

瑞利数(无量纲参数):

$$Ra = \frac{g\,\beta\,\Delta T\,L^3}{\nu\,\kappa}$$

$g$:重力加速度,$\beta$:热膨胀系数,$\Delta T$:上下面温度差,$L$:层高度,$\nu$:运动粘度,$\kappa$:热扩散率

温度场的对流–扩散方程(Boussinesq近似):

$$\frac{\partial T}{\partial t}+ \vec{u}\cdot\nabla T = \kappa\,\nabla^2 T$$

速度场(简化流函数法):

$$\vec{u}\propto Ra \cdot (T - T_{\text{mean}}) \cdot \hat{y}+ \text{粘性扩散}$$

边界条件:上面 $T=0$(冷),下面 $T=1$(热),侧面为周期边界。

自然界和工程中的应用

瑞利–贝纳尔对流是最基础的流体不稳定现象之一,广泛出现在多个领域:

气象·大气:积云(雷暴云)的发展、积雨云内的强上升气流是这种对流的大尺度版本。

地球内部:地球地幔对流(板块构造的驱动力)是数百万年时间尺度的瑞利–贝纳尔对流。

电子冷却·热设计:散热片的自然对流和数据中心冷却设计中,理解贝纳尔对流至关重要。

宇宙:太阳表面的"颗粒斑"(太阳米粒)是太阳大气对流单元,直径约1000 km的规则模式。

常见误解和注意事项

开始使用本模拟器时,有几个容易出错的地方需要注意。首先,不要认为"瑞利数越大,对流就越猛烈"。虽然超过临界值(约1708)会开始对流,但提高Ra时会经过从漂亮的涡卷到涡卷变形、分裂的"二次不稳定性"阶段,最终进入湍流。例如,从Ra=10⁴升至10⁵时,涡卷会因时间变化而波动,这是"过渡混沌"的开始。一次性提高太多会错过基本模式转变,所以应该逐步改变。

其次,不要忽视长宽比(宽度/高度)的设置。例如,长宽比为2时会自然形成2个对流涡卷,但设为3.5这样的不整数时,涡卷数不成整数,容易出现不稳定模式。现实和实机中层的"宽度"是固定的,改变模拟器的长宽比实际上是"改变实验装置或观测区域的大小",要有意识到这一点。

最后,要记得这个模拟是"二维"且"Boussinesq近似"的前提。在三维中,涡卷可能被替换为六边形贝纳尔单元;温度差很大或流体压缩性不可忽视时,这个近似不成立。这只是学习基本原理的理想化模型。

常见问题

如果瑞利数低于临界值(约1708),热传导会占优势,不会产生对流。将滑块设置为Ra > 2000,然后等待数秒,单元会开始出现。同时,如果初始温度扰动很小,也可能需要更长的观察时间。
我们采用简化流函数法,以温度差为驱动的浮力作为动力源,结合粘性扩散来近似速度场。虽然不是精确的NS方程求解,但能实时再现对流单元的形成、发展和湍流转变的主要行为。
将瑞利数设置为10⁵以上,长宽比设为4以上,这样多个单元会分裂、合并,形成时间上不规则的湍流状态,更容易观察。提高分辨率时,能看到更精细的涡结构。
本工具是面向教育和可视化的简化模型。虽然对理解热对流现象的定性特征很有帮助,但不能直接替代实验或定量设计。实用应用需采用专业CAE软件求解完整的纳维–斯托克斯方程。
努塞尔特数Nu表示对流热传达强度,当Ra < 1708时,Nu=1(纯热传导),Ra超过临界值后遵循 Nu ≈ 0.069 Ra^{1/3} 的幂律。在本模拟器中逐步提高Ra时,温度轮廓会从均匀梯度变成多个单元结构,实效热传输增强,可视觉上确认这种关系。
地幔对流的瑞利数约为Ra ≈ 10⁷〜10⁸,是板块构造的驱动力。在本模拟器中将Ra设置为最大值时,会出现非定常、复杂的单元结构,这在概念上接近地幔在长时间尺度上的运动。但真实地幔的粘性随深度和温度变化很大,所以只能定性类比。
理论上当Ra接近临界值时,单元宽高比约为1的形状最稳定。在长宽比大的容器中,多个单元会并排出现,单元个数大致等于宽高比的整数值。在模拟器中将长宽比从2→4→8递增,可观看1→2→4个涡卷形态的转变。
是的。Pr是热扩散和运动量扩散的比值。液体金属(Pr≈0.01)中热快速扩散,热边界层薄,单元易不稳定;硅油(Pr≈100〜1000)则粘性主导,单元趋向稳定的六边形模式。本模拟器假定一般流体(Pr≈1),不提供Pr直接调整,但OpenFOAM的buoyantBoussinesqPimpleFoam可进行Pr依赖解析。

实际应用示例

工业实际应用
半导体制造设备厂商(如东京电子)在外延生长炉的气流控制中采用本模拟器。通过调整Ra数来优化晶圆表面温度均匀性,降低膜厚波动。液晶面板制造的玻璃基板热处理工艺中,也用到对流单元的发生条件优化,有助于提升良率。

研究与教育应用
大学流体力学和热工学课程中,将本模拟器作为可视化工具,通过改变Ra观察对流模式的转变,帮助学生直观理解无量纲数与流动结构的关系。东京大学的学部实验采用本工具观察贝纳尔单元从形成到湍流的过程。地球物理研究中也用作地幔对流模拟实验的补充工具。

CAE分析衔接与工程定位
本模拟器定位于详细3D-CFD之前的阶段。设计初期通过参数化改变Ra,快速掌握对流发生条件和热传达特性的大致规律;随后用ANSYS Fluent或OpenFOAM进行高精度解析验证。这个工作流可将试错次数减少30%以上。

使用指南

  1. 调整瑞利数(Ra)滑块,设置底面加热强度。Ra=1000时为层流,Ra≥1708时贝纳尔对流开始
  2. 改变长宽比(AR),决定单元宽度与高度的比值。AR=2.0时形成标准六边形单元
  3. 勾选可视化选项(温度分布:chkT、速度矢量:chkV、流线:chkS),观察对流涡卷结构

具体计算示例

对水层(厚度2cm,下面温度40℃,上面温度20℃)进行模拟。运动粘性系数ν=1.0×10⁻⁶m²/s,热扩散率α=1.4×10⁻⁷m²/s,体膨胀系数β=2.1×10⁻⁴K⁻¹,得Ra=g×β×ΔT×d³/(ν×α)=4200。因为Ra > 1708,会发生对流。此时形成的对流涡卷数,当长宽比AR=4.0时,约观察到2个六边形单元

实务注意事项