关于“具有均匀加热侧壁的矩形腔内瞬态热毛细对流流动”研究的学术报告

本文旨在向广大科研工作者介绍一篇关于微重力条件下热毛细对流（Thermocapillary Convection）瞬态演化过程的重要研究。该工作由北京交通大学土木工程学院的陈恩惠（Enhui Chen）与徐丰（Feng Xu）教授（通讯作者）共同完成，发表于Physics of Fluids期刊第33卷，论文于2020年10月22日提交，2020年12月14日被接受，并于2021年1月6日正式在线发表，文章识别号为013602。

一、 研究的学术背景

本研究的核心科学领域是流体力学，特别是微重力流体物理与传热传质学中的热毛细对流现象。热毛细对流，又称Marangoni对流，是指由于液体自由表面（或界面）上表面张力随温度变化而产生的梯度所驱动的流动。当流体表面存在温度不均匀性时，表面张力差异会牵引流体从低表面张力区域（通常对应高温区）流向高表面张力区域（通常对应低温区），从而形成对流循环。这一现象在空间材料加工（如晶体生长、合金凝固）、微流体器件、3D打印以及许多地球上的微尺度流动过程中至关重要。

尽管针对热毛细对流稳态流动及其失稳特性的研究已有数十年历史，但对其瞬态动力学演化过程的理解仍不够清晰。在实际应用中，例如空间材料制备的初始加热阶段，流动处于瞬态发展过程，此阶段的流动与传热特性对最终产品质量有决定性影响。因此，深入研究瞬态热毛细对流，揭示其在不同参数下的演化规律和标度关系，具有重要的理论意义和应用价值。本研究正是在此背景下展开，其核心目标是：针对零重力条件下，一侧壁均匀加热的矩形腔内瞬态层流热毛细对流，通过标度分析和数值模拟相结合的方法，系统揭示流动随时间的演化规律，识别不同的流动区域和动力学演化阶段，并建立预测流动关键特征（如表面流速度、热层厚度）的标度律，最终通过数值模拟验证这些理论预测。

二、 研究的详细工作流程

本研究采用了理论分析与数值验证相结合的研究范式，主要包含两个核心步骤：系统的标度分析和详尽的数值模拟验证。

第一步：系统的标度分析 标度分析是本研究的理论核心，旨在通过物理量级的比较，推导出不同流动阶段控制方程中主导项的平衡关系，从而获得关键物理量（如速度、厚度、时间尺度）的幂律表达式。研究团队针对一个高宽比A<1、普朗特数Pr>1的二维矩形腔模型（如图1所示）进行了深入分析。其具体流程如下：

1. 建立物理模型与控制方程：首先明确了研究问题：一个初始静止、温度为T0的流体，在零重力条件下，其右侧壁面突然施加均匀热流q，左侧壁面保持恒温T0，底部绝热，顶部为可施加热毛细应力的平坦自由表面（视为刚性、不可渗透）。控制方程为纳维-斯托克斯方程和能量方程，边界条件包括热毛细应力条件（公式39）。
2. 识别关键物理过程与时间尺度：分析从加热瞬间开始，流动的演化涉及多个相互耦合的物理过程：a) 靠近热壁的垂直热层的形成与增厚；b) 由表面张力梯度驱动的表面流的产生与发展；c) 表面流通过粘性效应带动下方流体形成的水平粘性层的扩散；d) 表面流引起的压力梯度驱动的垂直回流。研究首先推导了这些特征层的厚度随时间演化的标度，例如垂直热层厚度δ_t ~ κ^{¹⁄₂}t^{¹⁄₂}，水平粘性层厚度δ_ν ~ (Pr κ)^{¹⁄₂}t^{¹⁄₂}，并定义了它们达到稳态的特征时间尺度，如粘性层扩散至整个腔高的时间t_ν，热层达到稳态厚度的时间t_κ。
3. 划分流动区域与动力学演化阶段：这是本研究最具创新性的部分。研究团队深入分析了表面流动量方程中非定常项、粘性项、热毛细力项之间的主导平衡关系，这些平衡关系取决于时间以及三个关键无量纲参数：马兰戈尼数Ma、普朗特数Pr和腔体高宽比A。通过细致的量级比较，研究预测了六种不同的动力学演化路径，在Ma-Pr-A参数空间中形成了六个区域（见图2）。每条路径都包含一系列按时间顺序出现的流动区域，例如：
   • 非定常惯性-热毛细平衡区域：在极早期，表面流加速度项与热毛细力项平衡。
   • 非定常粘性-热毛细平衡区域：随后，粘性力项与热毛细力项平衡，此区域可能根据温度梯度特征长度尺度的不同（是热扩散尺度δ_t还是表面流平流尺度u*t）进一步细分。
   • 稳态粘性-热毛细平衡区域：最终，流动达到稳态，粘性力与热毛细力平衡。 研究详细推导了每个区域对应的表面流速度标度律（例如，公式21的u_2，公式27的u_5，公式33的u_9等）以及区域之间的转换时间（如t_1, t_2, t_3, … t_8）。这些标度律清晰地展示了速度如何随时间（t的幂次）和系统参数（Ma, Pr, A）变化。
4. 推导稳态与传热标度：对于充分发展后的稳态，研究推导了稳态表面流速度u_s、垂直回流速度v_s、稳态热层厚度δ_ts和粘性层厚度δ_νs的标度律（公式13），以及努塞尔数Nu的标度律（公式17），并给出了垂直热层显著存在的判据Ma > A^{-3}（公式14）。

第二步：二维数值模拟验证 为了验证标度分析的理论预测，研究团队进行了系统的二维数值模拟。

1. 数值方法与模型设置：采用有限体积法求解无量纲化的控制方程（公式35-38）。压力-速度耦合采用SIMPLE算法，对流项采用QUICK格式离散，扩散项采用中心差分，时间项采用二阶格式。热毛细效应通过自由表面的Marangoni应力边界条件（∂u/∂y = Ma^{¹⁄₃} ∂T*/∂x*）实现。计算域为矩形腔，网格在边界处加密以确保精度。
2. 参数范围与网格无关性验证：数值模拟覆盖了广泛的参数空间：Ma从10^2到10^5，Pr从10到100，高宽比A从0.14到0.33（对应1/7到1/3）。针对不同的A，进行了细致的网格和时间步长依赖性测试（见表II），确保计算结果与网格和时间步长无关，最大速度差异小于1%。例如，对于A=1/6，采用100×300的网格和0.04的无量纲时间步长。
3. 方法验证：为了确认数值方法的可靠性，研究者将本文的模拟结果（Ma=6330, A=0.18, Pr=8.4）与Gillon和Homsy (1996)的经典实验结果进行了对比（见图5）。两者在速度剖面上表现出良好的一致性，从而验证了当前数值模型模拟热毛细对流的准确性。
4. 数据提取与分析流程：数值模拟提供了完整的瞬态流场和温度场数据。为了验证标度律，研究者从模拟结果中提取了关键物理量：
   • 垂直热层厚度：定义为从加热壁面到温度T* = 0.1的等温线（在y* = 0.25A高度处）的距离。
   • 表面流速度：追踪自由表面（y* = A）上特定等温线（如T* = 0.05）的水平运动速度来代表表面流特征速度。
   • 稳态表面流速度与努塞尔数：在充分发展的稳态（t* = 1000）下，测量表面固定点的速度以及加热壁面的平均传热系数（Nu）。

三、 研究的主要结果

数值模拟的结果全面支持了标度分析的理论预测，具体如下：

1. 流场结构可视化验证：图6展示了在不同动力学演化区域（I至VI）下的典型流场（等温线和流线）。结果直观地证实了理论预测：当Ma < A^{-3}（区域I）时，流动很弱，没有明显的垂直热层，传热以导热主导；随着Ma增大，对流增强，等温线扭曲，流场对称性被打破，垂直热层变得明显且越来越薄（区域III-VI），特别是在Ma非常大的区域VI，热层极薄，流动高度集中在加热壁角附近。这些流场结构的变化与理论划分的流动区域特征完全吻合。
2. 瞬态垂直热层厚度标度律验证：图7展示了在所有验证的动力学演化路径（II-VI）中，无量纲化的垂直热层厚度δ_t* Ma^{¹⁄₃}与无量纲时间的平方根t^{¹⁄₂}呈完美的线性关系。这与标度律δ_t ~ κ^{¹⁄₂}t^{¹⁄₂}（公式8）的预测一致。通过线性拟合得到比例系数约为1.95，将标度关系定量化为公式：δ_t Ma^{¹⁄₃} ≈ 1.95 t*^{¹⁄₂}。
3. 瞬态表面流速度标度律验证：
   • 对于u_2标度律（公式21，预测u_2 ~ Ma Pr^{¹⁄₂} κ^{³⁄₂} t^{¹⁄₂} / l^2），图8显示，从多个演化路径（II-VI）的早期数据中提取的u_2* Pr^{-¹⁄₂}与t*^{¹⁄₂}呈线性关系，验证了该标度。
   • 对于u_5标度律（公式27，预测u_5 ~ Ma^{¹⁄₃} κ^{¹⁄₂} A^{¹⁄₂} / t^{¹⁄₂}），图9展示了在动力学演化V中，u_5* A^{-¹⁄₂}与t*^{-¹⁄₂}的线性关系，证实了该阶段速度随时间衰减的规律。
   • 对于u_9标度律（公式33，预测u_9 ~ Ma^{¹⁄₃} Pr^{¹⁄₄} κ^{³⁄₄} / (l^{¹⁄₂} t^{¹⁄₄})），图10展示了在演化V和VI中，u_9* Ma^{¹⁄₆} Pr^{-¹⁄₄}与t*^{-¹⁄₄}的线性关系，验证了另一种速度衰减模式。
4. 稳态表面流速度标度律验证：图11汇总了所有达到稳态的演化路径（II-VI）中，稳态表面流速度u_s*与高宽比A的关系。数据点清晰地落在一条通过原点的直线上，与标度律u_s ~ Ma^{²⁄₃} κ A / l（公式13）的预测（即u_s应与A成正比）高度一致。线性拟合得到u_s ≈ 0.85A。
5. 稳态传热标度律验证：图12显示了稳态努塞尔数Nu与Ma^{¹⁄₃}A的线性关系，完美验证了标度律Nu ~ Ma^{¹⁄₃}A（公式17）。拟合结果为Nu ≈ 160 Ma^{¹⁄₃}A，为预测矩形腔在给定参数下的传热能力提供了定量公式。

这些结果层层递进：首先，瞬态热层和速度的演化标度得到了验证，说明了理论对瞬态过程描述的准确性；其次，稳态速度和传热标度的验证，表明理论能够正确预测流动的最终状态；最后，不同参数下流场结构的差异，直观地印证了理论所划分的不同流动区域是真实存在的。

四、 研究的结论与价值

本研究的核心结论是：通过系统的标度分析，成功预测并数值验证了零重力条件下矩形腔内瞬态热毛细对流存在多种不同的动力学演化路径和流动区域。具体而言，当Ma > A^{-3}时，根据Ma、Pr、A的不同组合，表面流的瞬态演化可划分为六种典型路径。在每条路径中，流动会依次经历由不同物理力主导的平衡阶段（如非定常惯性-热毛细平衡、非定常粘性-热毛细平衡），最终达到稳态的粘性-热毛细平衡。研究给出了每个阶段的表面流速度标度律、特征层厚度标度律以及阶段转换的时间尺度标度律。

本研究的科学价值在于：首次对均匀热流边界条件下矩形腔热毛细对流的完整瞬态演化过程进行了系统的标度分析，构建了一个清晰的理论框架，将复杂的瞬态流动行为分类并参数化。所推导出的标度律为理解和预测此类流动的关键物理量（速度、厚度、时间尺度、传热率）提供了简洁而有力的理论工具。其应用价值体现在：可为空间材料加工、微重力流体管理、微流体系统设计等涉及瞬态热毛细对流的工程问题提供理论指导和快速评估依据，例如预测加热初始阶段的流动发展时间、表面速度峰值以及稳态传热效率。

五、 研究的亮点

1. 理论创新性：研究超越了传统的稳态分析，聚焦于瞬态演化这一相对空白但至关重要的领域。提出的基于Ma-Pr-A参数空间的六路径动力学演化分类图（图2）是该研究最突出的理论贡献，为理解复杂参数依赖的瞬态行为提供了“地图”。
2. 系统性与完整性：标度分析极为详尽，不仅考虑了表面流，还耦合分析了垂直热层、水平粘性层、垂直回流及其相互作用，推导出了一整套覆盖从初始瞬态到最终稳态的完整标度律体系。
3. 严谨的验证：通过覆盖广泛参数空间的精细数值模拟，对多个关键标度律（包括瞬态和稳态）进行了交叉验证，数据与理论预测符合度极高，增强了结论的可靠性。
4. 明确的物理图像：研究将复杂的数学推导转化为清晰的物理图像（如不同阶段的力平衡）和直观的演化序列图（图3），有助于读者理解和应用其理论成果。

六、 其他有价值的内容

研究在附录（Nomenclature）中提供了完整的符号列表，便于读者查阅。此外，文中简要回顾了热毛细对流领域的重要前期工作（如Ostrach, Cowley & Davis, Birikh等人的贡献），并将自己的研究置于该学术脉络中，体现了研究的继承性与发展性。作者也明确指出了模型的适用范围（零重力或浮力效应可忽略的条件α_t/(ρgβl^2) >> 1），并提及了在中国博士后科学基金和国家111计划资助下完成，体现了研究的资助背景。