本文介绍由中国科学院力学研究所微重力重点实验室及中国科学院大学工程科学学院的姜欢、段俐与康琦合作完成的研究。该研究成果以题为“Instabilities of thermocapillary–buoyancy convection in open rectangular liquid layers”的论文形式，发表于2017年的《Chinese Physics B》期刊第26卷第11期。
一、 学术背景
本研究属于流体力学，特别是热对流稳定性研究领域，聚焦于热毛细-浮力耦合对流的不稳定性问题。在微重力或常规重力环境下，具有自由表面的液层在水平温度梯度驱动下，会产生由表面张力梯度驱动的热毛细对流（Thermocapillary convection）以及由密度差驱动的浮力对流（Buoyancy convection）。这两种对流机制的耦合作用，会导致复杂的流动失稳现象，产生诸如 hydrothermal waves（热毛细波）、纵向对流卷（longitudinal rolls）等时空结构。理解这些不稳定性的产生条件、形态特征及转变规律，对于空间流体管理、晶体生长、薄膜涂覆等工程应用以及基础流体物理研究都具有重要意义。
研究的直接背景源于Smith和Davis在1983年通过线性稳定性分析预测，在水平温度梯度作用下，无限延伸的薄液层中会出现从冷壁斜向热壁传播的hydrothermal waves。随后，Riley和Neitzel等人在实验中观察到了此类波动。然而，流动失稳的模式强烈依赖于液层厚度、普朗特数（Pr）、马兰戈尼数（Ma）和动态邦德数（Bd）等参数。特别是在常规重力下，浮力效应不可忽略，其与热毛细效应的竞争（由Bd数衡量，Bd>1时浮力主导，Bd时热毛细主导）会显著改变失稳的形态。此前的研究虽多有涉及，但本文作者指出，这是首次基于其矩形液层实验模型，对不同对流失稳模式进行的全面总结性研究。本研究旨在通过实验，系统观测开式矩形液层中热毛细-浮力对流在转捩过程中的不稳定性，详细区分不同流动模式，并利用创新的时空演化分析方法直观展示行波传播特性。
二、 详细工作流程
本研究为实验研究，主要流程包括实验系统搭建、实验参数设置、流动可视化观测与数据采集、以及后续的图像与数据分析。
1. 实验系统与对象： 研究团队构建了一套专门的实验系统，核心是一个开式矩形腔体。腔体尺寸为52 mm（长）× 108 mm（宽）× 8 mm（高），两侧壁为铜板，分别作为冷壁和热壁，前后壁及底部为K9光学玻璃以便观测。温度控制系统包含用于加热的电热膜、用于冷却的半导体制冷器、直流电源以及两个Eurotherm温度控制器（精度0.1°C），以实现对两侧壁温的精确控制和高差（ΔT）的线性增加。
研究的对象是不同厚度的硅油液层。工作流体为1 cSt, 1.5 cSt和2 cSt的Shin-Etsu KF96硅油，对应的普朗特数（Pr）分别为16, 25和29。通过改变液层深度（d），可以改变纵横比（γ = L/d，其中L=52 mm为冷热壁间距）和动态邦德数（Bd = ρgβd²/(∂σ/∂T)），从而系统研究浮力与热毛细力的相对重要性。表征热毛细力强度的无量纲数为马兰戈尼数（Ma = (∂σ/∂T)(ΔT/L)d²/(ρνκ)）。
2. 流动观测与数据采集： 本研究采用了一种创新的观测方法：红外热成像技术。使用配备红外微距镜头的FLIR T420热像仪（分辨率320×240像素，热灵敏度0.05°C，采样频率15 fps）直接观测并记录液层自由表面的温度场（即表面波）。微距镜头将观测区域聚焦在液层自由表面中部一个32 mm × 24 mm的范围内，以减少固体边界的干扰。这是一种非侵入式测量方法，不会干扰流场运动。
3. 实验步骤与数据处理： 对于每组给定的Pr和γ（通过设定液深d实现），实验以一个“实验周期”运行。在此周期内，采用线性增加温度差的模式，逐步增大Ma数。热像仪持续记录自由表面的红外图像视频序列。 获得红外图像序列后，研究进行了两方面的分析： * 流动模式识别与分类： 直接观察不同Ma数（即不同ΔT）下自由表面温度场的形态变化，识别并分类出现的对流结构，如稳态二维流、对流纵向卷、hydrothermal waves、表面流（surface flow）等。 * 时空演化分析： 这是本研究采用的一种特色分析方法。在自由表面观测区域内，分别选取一条平行于温度梯度方向（x方向，从冷壁到热壁）的竖直线和一条垂直于温度梯度方向（y方向）的水平线。将图像序列沿时间轴堆叠，即可得到沿该线温度随时间变化的时空演化图。通过分析这些图中斜条纹的特征，可以直观地展示波动的传播过程、规律，并定量计算行波的波速和传播角。具体计算方法是：从x方向和y方向的时空图中分别得到波速分量v_x和v_y，进而合成总波速v和传播角θ（θ为波矢与温度梯度方向，即x方向的夹角）。
三、 主要结果
实验系统观测并区分了热毛细-浮力对流转捩过程中出现的三种主要不稳定性模式，并利用时空演化图分析了其动态特性。
1. 对流纵向卷： 这是在整个对流转捩过程中广泛存在的一种结构。当γ ≥ 34.7时，在矩形几何结构中可以观察到这种由表面张力梯度产生的卷状结构。它们独立于Pr数，但随液层纵横比γ变化。研究发现，纵向卷的无量纲波数α = 2πd/λ（λ为波长）随着γ的增大而呈增长趋势，表明不同厚度液层具有不同的流动构型。随着Ma数的增加，这些纵向卷会发生前后振荡。当温度差足够大时，对流卷的结构会被表面波动破坏，流动向振荡对流转变。
2. 行波： 当Bd < 1（热毛细效应主导）且20.8 ≤ γ ≤ 34.7时，可以观察到明显的斜向行波，即hydrothermal waves。 * 特征： 这些波从冷壁斜向热壁传播，与Smith和Davis的线性稳定性理论预测一致。它们“骑”在对流纵向卷之上出现。例如，在γ=34.7， Pr=16的条件下，当Ma较小时（Ma=871），流场仅呈现二维的对流纵向卷；当Ma增至临界值（Ma=1574）时，流动失稳为三维流，出现行波；当Ma更大时（Ma=1807），会出现行波叠加和振荡对流。 * 波源： 在γ=26，不同Pr的条件下，观测到对称分布的两列波在自由表面斜向相交重叠，对称轴即为波源线。波源线从容器的这一端延伸到另一端，产生两列波。 * 相互作用： 在转捩过程中，存在行波与对流纵向卷的明显相互作用。例如在γ=20.8， Pr=25， Ma=2475时，行波叠加在纵向卷结构上；当Ma增至3079时，强烈的表面波动和振荡对流打断了纵向卷结构，造成大量缺陷，形成非线性相互作用，这属于热毛细主导的I型振荡流。在此条件下（Bd≈1），还会观察到由浮力效应引起的、源自热边界层的微弱表面扰动叠加在行波之上。
3. 表面流与扰动波： 当γ ≤ 20.8，即Bd ≥ 1（浮力效应主导）时，浮力对流主导，在自由表面出现体流（bulk fluid flow），即表面流。 * 特征： 表面流从热壁流向冷壁，并随着Ma数增大而增强。其成因是热边界层底部的流体受热膨胀，快速上升并在自由表面形成从热侧到冷侧的流动。 * 扰动波： 当Ma数足够大时，不稳定的热边界层还会导致扰动波的出现。这种波类似于Garcimartín等人观察到的边界层失稳波，源于热壁附近边界层内的小扰动被向下游拖曳并放大。当Ma数非常大时，强烈的浮力效应扩展到整个流场，产生显著的振荡流。表面流与各种扰动波的重叠和相互作用构成了浮力主导的II型振荡流。在此条件下，hydrothermal waves虽然仍存在，但被浮力对流的不稳定性所掩盖。
4. 稳定性相图： 基于实验结果，研究者绘制了以纵横比γ和马兰戈尼数Ma为参数的稳定性相图。相图清晰展示了流动随Ma数增加而失稳的路径：从小Ma下的稳态流，到中等Ma下的三维流（分为Bd<1时的hydrothermal waves区和Bd>1时的表面流区），再到高Ma下的振荡流（分为I型和II型）。对于不同的Pr数（16, 25, 29），流动特性相似，但失稳的临界Ma数有所不同。
5. 时空演化分析结果： 时空演化图清晰地展示了波动随Ma数演变的三个阶段：无波阶段、规则行波阶段、以及伴随振荡对流的行波阶段。图12（γ=26, Pr=25）和图13（γ=17.3, Pr=25）分别展示了热毛细主导和浮力主导条件下的典型演化过程。通过该方法，研究者定量测量了不同实验条件下hydrothermal waves的传播速度（v）和传播角（θ），数据列于文中表2。结果表明，在相同γ、不同Pr条件下，波速大小差异不大；波速和传播角在一定程度上更依赖于γ的变化，这与不同流场结构有关。
四、 结论
本研究通过红外热成像技术，对开式矩形液层中热毛细-浮力对流的不稳定性进行了详细的实验观测和分析。主要结论如下： 1. 对流纵向卷存在于整个对流转捩过程中。 2. 当Bd < 1（热毛细主导）时，自由表面出现斜向传播的hydrothermal waves；当Bd ≥ 1（浮力主导）时，流场中出现表面流和扰动波。 3. 时空演化分析能有效展示对流转捩过程中波动的特征，并可定量测量行波的波速和传播角。Hydrothermal waves的传播特性更依赖于纵横比γ，因为γ决定了不同的流场结构。 4. 不同类型的对流不稳定性取决于不同的流场结构，而这在常规重力环境下普遍受到浮力效应的影响。
五、 研究亮点
六、 其他有价值内容
研究还提及了部分流动结构的PIV（粒子图像测速）截面速度场结果（图3b），与红外表面温度场观测形成了互补，有助于理解三维流动结构。此外，文中引用了大量该领域的经典与近期文献，显示了研究团队对学术背景的深刻把握，并将自身工作置于更广阔的学术对话之中。