Numerical Solution of Natural Convection on a Heated Vertical Plate at Low Diffusivity Ratios

技术摘要：低扩散率比条件下加热垂直平板的自然对流数值解

问题陈述
本研究调查了等温垂直平板在低扩散率比条件下的瞬态自然对流流动，特别针对热扩散率显著高于动量扩散率（低普朗特数，$Pr \ll 1$）的流体。虽然垂直平板上的自然对流已是成熟的研究基准，但在低 $Pr$数状态下（这类流体常见于液态金属和熔盐），其瞬态演化过程和边界层结构与中高$Pr$ 数流体（如空气或油）存在显著差异。本文旨在填补文献中关于在这些高扩散条件下，热边界层和传热机制的具体瞬态行为的研究空白；在这些条件下，热扩散主导了由浮力驱动的动量输运。

研究方法
本研究采用基于有限差分法（FDM）的数值方法，求解基于 Boussinesq 近似推导出的无量纲控制方程组。该方程组包括连续性方程、动量方程和能量方程。

• 离散化： 空间域采用结构化均匀笛卡尔网格。扩散项使用二阶中心有限差分进行近似，以保持空间精度；而非线性对流项则使用一阶上风格式进行离散，以抑制具有陡峭梯度的区域中可能出现的数值振荡。
• 时间积分： 采用显式欧拉格式进行时间推进。为确保数值稳定性，时间步长（$\Delta\tau$）严格遵循库朗-弗里德里希斯-列维（CFL）条件，由于增强的热扩散作用，该条件对于低 $Pr$ 流体而言尤为苛刻。
• 验证： 通过网格无关性研究（收敛于 $100 \times 100$ 网格）以及与经典相关性及近期文献数据进行对比，验证了数值方法的有效性。模拟涵盖了 $Pr$数为$0.01$至$0.71$的范围，并固定雷诺数（$Re=10^2$）和瑞利数（$Ra=10^5$）。

主要结果
数值模拟揭示了低 $Pr$ 数机制下与常规对流截然不同的物理行为：

1. 边界层结构： 随着普朗特数降低，热边界层变得显著厚于流体动力学（速度）边界层。对于极低 $Pr$数（例如$Pr=0.01$），温度场向流体域深处渗透，导致热衰减过程非常缓慢。
2. 速度与温度分布： 在低 $Pr$流体中，峰值速度向远离加热壁面的方向偏移，且动量边界层变宽。这是因为浮力作用分布在更大的热区域内，导致靠近壁面的局部加速减弱。相反，高$Pr$ 数流体表现出更陡峭的热梯度和更靠近表面的速度峰值。
3. 传热速率： 随着普朗特数下降，局部努塞尔数（$Nu_x$）降低。热扩散的主导作用将热能快速扩散至流体中，从而降低了法向壁面温度梯度，进而降低了对流传热率。
4. 瞬态演化： 低 $Pr$数对流表现出延迟的热响应。$Pr=0.01$的流体表现出向热平衡趋近的渐近收敛过程较为缓慢，而中高$Pr$ 数流体则稳定得更快。
5. 定标律： 结果证实了热边界层厚度（$\delta_t$）与普朗特数的平方根成反比（$\delta_t \sim Pr^{-1/2}$）的定标关系，这与扩散主导机制下的理论预测一致。

意义与主张
本文声称提供了一个稳健的数值框架，能够准确捕捉低 $Pr$ 数自然对流中固有的强扩散热输运机制。通过利用具有特定稳定策略（上风格式和严格遵守 CFL 条件）的显式有限差分法，研究表明可以解析这些机制中速度场与温度场之间复杂的耦合关系。

作者断言，其研究结果为涉及液态金属冷却系统、核反应堆热管理及先进热能技术的工程应用提供了物理见解。本工作强调，理解低扩散率流体中热场的瞬态演化，对于设计以热扩散决定输运行为的被动冷却系统和换热器至关重要。研究结论指出，所提出的数值模型是可靠的，这从其与基准数据的吻合程度，以及能够重现预期的物理趋势（如热层增厚和低 $Pr$ 数条件下传热系数降低）得到了证实。