一维、二维、三维和四维中的热对流

《Journal of Fluid Mechanics》：Thermal convection in one, two, three and four dimensions

　　

摘要：研究人员通过直接数值模拟（DNS）研究了对流维度对流动性质的影响。研究人员通过综合考虑一维、二维、三维和四维的结果，提请关注几个一般性原则。具体而言，研究人员探讨了两个实际方面：（i）瞬态时间，即流动达到稳态所需的时间；以及（ii）对所谓的最终状态的可能启示。

一、研究背景与问题

湍流热对流的直接数值模拟（DNS）在基础研究中占据日益重要的地位。随着这一发展，数值收敛性、分辨率、达到稳态所需的瞬态时间等问题愈发受到关注。以往研究主要聚焦于二维和三维瑞利-贝纳德对流（RBC），但对维度如何系统性地影响对流动力学特性，特别是瞬态时间与最终状态（ultimate state）的判定，仍缺乏整体性认识。Pandey和Sreenivasan（简称PS）此前在二维RBC研究中发现，伴随大尺度脉动热输运的长时间瞬态过程使最终状态的观测变得相当不确定。最终状态指努塞尔数（Nusselt number, Nu）标度超越经典的1/3指数并趋近1/2幂次（可能带有对数修正）的渐近高瑞利数（Rayleigh number, Ra）区间。然而，三维与四维情况下的瞬态行为是否不同、更高维度能否在更低瑞利数下趋近最终状态、以及一维情形是否存在对流，均是尚未解决的关键问题。

二、研究目标与意义

本研究旨在通过系统比较一维至四维热对流的动力学行为，阐明维度对对流特性的影响规律。研究重点关注两方面：一是瞬态时间（即流动达到统计稳态所需的时间）随维度的变化趋势；二是高维度对流是否能为最终状态的研究提供新视角。物理上，四维虽在真实宇宙中不存在，但统计力学中伊辛（Ising）模型的上临界维度恰为4——四维及以上的临界行为可由平均场理论描述。相应地，湍流在四维中可能呈现无异常的正常标度行为。因此，研究四维对流有助于检验这一推测，并为三维最终状态问题提供可能的参考框架。该工作发表于《Journal of Fluid Mechanics》。

三、主要关键技术方法

研究人员利用直接数值模拟（DNS）方法，在二维、三维和四维空间中求解奥伯贝克-布辛涅斯克（Oberbeck–Boussinesq）方程组。二维方形封闭域模拟采用谱元法求解器Nek5000，水平周期边界情形及三维、四维模拟采用有限差分求解器Dhara（基于交错网格的二阶空间离散与三阶龙格-库塔时间推进）。四维模拟在周期超立方体中进行，网格规模最高达256?，在ALCF的Polaris超算上使用128块NVIDIA A100 GPU执行。计算中通过全域平均动能E的时间演化判断瞬态过程，以努塞尔数（Nu）和雷诺数（Reynolds number, Re）量化热与动量输运。研究覆盖Ra范围因维度而异：二维为10?至1012，三维为10?至101?，四维为10?至2×10?，Pr数范围为0.021至10（表1和表2给出了详细参数）。2.1节进一步说明二维模拟采用的谱元法Nek5000中，计算域分解为N?个单元并使用N阶勒让德多项式离散；三、四维模拟采用一致的有限差分框架以保证可比性（表2）。

四、研究结果

1. 基本响应量

研究中以全域平均动能E=〈|u|2/2〉_V/u_f2描述对流的瞬态演化。稳态时通过雷诺数Re = u_rms H/ν表征流动强度，其中u_rms为均方根速度。努塞尔数通过Nu = 1 + H/(κΔT)〈u_zT〉_V,t计算热输运效率。不同方法计算的努塞尔数一致性在百分之几以内，验证了空间和时间收敛性。

2. 一维流动

在一维不可压缩情形中，热传导是上下壁面间唯一的输运机制。连续性方程与上下不可穿透边界在数学上严格排除了产生经典对流所需的定常分岔，线性稳定性分析和DNS均证实了这一约束。可压缩情形中压力扰动可传播，但纯纵向压缩不产生单向热输运，因此一维条件下在任何Ra和Pr下都不会出现对流，瞬态时间可视为无限长。

3. 二维对流

在封闭方形域中，瞬态时间t_trns随Ra增大而增长，标度指数在0.59至0.71之间（取决于Pr），与Lindborg理论一致，但比例系数取决于Pr。研究还发现t_trns与Re近似线性相关，综合所有Pr数据的最佳拟合为t_trns = 0.0044 Re^1.02。能量演化经历初始指数增长阶段和随后的缓慢指数增长阶段，达到名义稳态后仍伴随巨大脉动。瞬态时间随Re的线性增长对高Re二维对流的探索构成严重计算瓶颈。

4. 三维对流

在封闭立方体域中，能量E在初始阶段快速增长并开始脉动，约50至70个自由落体时间t_f后达到合理的E_av。初始增长同样具有指数特性。不同Pr下的演化特征定性相似，在长方形长方体（2.4H×0.8H）中结果亦然。与二维不同，三维中E_av随Ra变化较弱：低Pr流体中E_av随Ra略微下降，高Pr流体中则略微上升。瞬态时间随Pr增大而略有延长，即三维中瞬态时间随Re减小而微弱增加。总体而言，三维情形对Ra或Re的依赖趋势远不如二维强烈。

5. 四维对流

四维模拟在周期超立方体中进行，Pr=1，Ra范围10?至2×10?。能量E的演化与三维相似：初始指数增长后进入统计稳态，瞬态时间约40至50 t_f，略小于三维但总体相当。关于全局响应量，在给定Ra下Nu随维度D增加而增大，但三维与四维之间的差异小于二维与三维之间的差异。Ra≥5×10?时，三维标度为Nu=0.22Ra^0.27，四维为Nu=0.21Ra^0.28，指数与前因子均非常接近；二维则为Nu=0.063Ra^0.29，且伴随巨大时间脉动。Ra=2×10?时，三维的Nu约为二维的2.4倍，而四维仅比三维高15%。Re方面，二维（Re=0.071Ra^0.59）明显高于三维和四维（Re≈0.36Ra^0.46），后两者几乎不可区分。二维高Re源于涡量守恒禁止涡拉伸、能量逆级联导致大尺度环流增强均方根速度；三维和四维中涡拉伸可实现能量正向级联，大尺度结构减弱。二维周期盒数据在Ra≈10?附近出现明显间断，归因于强带状流的出现。四维结果在Ra^0.28的准经典标度下持续至2×10?，表明空间维度的增加并未绕过经典热输运机制。

五、总结与讨论

研究揭示了维度对热对流动力学的根本性影响。一维中严格约束在物理和数学上排除了对流。二维瞬态时间异常漫长且近似线性正比于Re，构成探索极高Ra区间的严重计算瓶颈。随着维度增至三维和四维，瞬态时间对Re的依赖显著减弱。四维动能演化与三维高度相似，表明充足的空间自由度使稳态过渡更高效、对热强制的依赖更弱。研究人员基于一至四维的瞬态时间数据构建了经验关系式：瞬态时间/t_f = βRe^η，其中η=(D-3)(D-4)/[2(D-1)]，三维和四维时β≈40，二维时β≈4.4×10?3。

研究结论可归纳如下：第一，鉴于二维中的长瞬态时间，在二维中开展决定性的最终状态研究需要过长计算时间，在极高Re下可能永远无法实现。第二，三维对流的瞬态时间约为50个自由落体时间量级。第三，若能视三维对流为四维对流的微扰，研究四维对流可能有助于解决包括最终状态在内的各种标度问题。然而，即使在瞬态时间可确定的精度水平下，三维与四维情形之间预计不会出现巨大差异。此外，三维和四维中Nu与Re标度指数的稳健性表明，空间维度不会改变极端高热强制下对流的基本经典机制——所有维度趋近最终状态可能都需要同样高的热强制水平。