本文研究由来自中国科学院力学研究所国家微重力实验室的Tang Zemei与Hu Wenrui共同完成，发表于2007年12月的《Science in China Series G: Physics, Mechanics & Astronomy》期刊上。该研究得到国家自然科学基金的资助。

一、 学术背景 本研究隶属于微重力流体物理领域，具体聚焦于热毛细对流（Thermocapillary convection）及其失稳机制。随着载人航天和空间科学的发展，微重力环境下的流体行为，尤其是热毛细对流，因其在空间材料加工等应用中的重要性，以及在基础流体物理研究中的独特价值，成为过去几十年的研究热点。热毛细对流由流体自由表面（气液界面）的表面张力梯度驱动，而表面张力通常对温度敏感。在存在温度梯度的条件下，会产生沿表面的切向应力，驱动流体运动，即马兰戈尼对流（Marangoni convection）。

当施加的温度差超过某个临界值时，稳态的热毛细对流会失稳，转变为振荡流。其中一种重要的失稳模式被称为“热液波”（Hydrothermal Wave, HTW）。Smith和Davis在1983年的理论分析中指出，对于大普朗特数（Prandtl number, Pr）流体，在平行于自由表面存在温度梯度的无限延伸液层中，可能产生一种向上游（即逆着表面主流方向）传播的振荡波，即热液波，其能量来源于液层内的垂直温度梯度。然而，在地面重力环境中，热毛细对流常与浮力对流耦合，使得现象复杂化，难以观察到纯粹的热液波。此外，对于大普朗特数流体（如本研究使用的10厘斯托克硅油，Pr=105.6），其高粘性、低导热特性可能导致其振荡对流行为与常用流体（Pr < 20）不同。因此，本研究旨在通过数值模拟，在微重力条件下，探究大普朗特数流体在浅槽内的热毛细对流失稳，特别是热液波的产生机制、传播特性及其激发机理，为空间实验的开展提供理论依据和前期准备。

二、 详细研究流程 本研究主要采用二维数值模拟方法，具体工作流程如下：

1. 物理模型与数学建模：

   • 研究对象：研究一个充满10 cst硅油（Pr = 105.6）的矩形浅槽。液层高度（H）为4毫米，长度（L）为52毫米，长高比（Aspect Ratio, L/H）为13。
   • 边界条件：左侧壁面为高温（Th），右侧壁面为低温（Tc），温差为ΔT = Th - Tc。底部为刚性绝热壁面，顶部为可变形自由表面（与被动气体接触），且假设为绝热。表面张力σ随温度线性变化。
   • 控制方程：在引入涡量（ω）和流函数（ψ）将连续性方程自动满足后，将动量方程和能量方程转化为以涡量和流函数表示的非定常、二维、不可压缩形式。控制方程中考虑了马兰戈尼数（Ma）和雷诺数（Re）等无量纲参数。自由表面的边界条件包括切向应力平衡（引入表面张力梯度效应）和法向应力平衡（用于确定自由表面形状）。
   • 初始条件：从没有温差、无对流的静止状态开始。

2. 数值方法与程序验证：

   • 求解方法：采用分数步杂交有限元法求解控制方程及其边界条件。计算域在x和z方向分别划分为121和21个网格，共生成4800个三角形单元，2541个节点。
   • 求解流程：初始时刻，自由表面设为平面，两侧壁温度相同（T0）。随后，左侧壁温度以0.5℃/s的速率从T0升高至目标值（如T0+25℃）并保持恒定，右侧壁温度保持为T0。程序迭代求解流场、温度场，并耦合法向应力平衡方程更新自由表面形状。
   • 程序验证：为确保数值方法的可靠性，作者将程序计算结果与Carpenter和Homsy（1990）关于方腔内热毛细对流的研究结果进行了对比。对比参数包括流函数最小值、涡量核心值、表面速度以及努塞尔数等。对于Pr=1, 10, 30, 50等多种情况，本文结果与文献结果吻合良好。此外，还将计算得到的零重力下液层中截面（ξ=0.5）的平均速度分布与Smith和Davis提出的基本回流状态理论解进行了比较，结果高度一致。这些验证工作充分证明了所开发数值程序的准确性。

3. 数值模拟与分析过程：

   • 参数设定：重点模拟微重力（g=0）环境，以排除浮力效应的影响，专注于热毛细效应。研究的关键变量是施加的侧壁温差ΔT。
   • 数据处理：将随时间变化的瞬时物理量（如速度u、温度T、流函数ψ、压力p）分解为时间平均的稳态部分（下标0）和振荡的扰动部分（上标‘）。例如，u = u0(x, z) + u’(x, z, t)。通过分析扰动量的演变来识别振荡对流的发生。
   • 临界现象识别：逐步增加ΔT，观察流场和温度场从稳态到振荡态的转变。由于扰动非常微弱，无法直接从瞬时总场的流线或等温线图中直观判断振荡开始，因此通过监测自由表面上几个典型位置（如x1=8.233 mm, x2=21.233 mm, x3=34.233 mm）的扰动速度u’的频谱来精确判断振荡频率和临界温差（ΔTc）。
   • 详细分析：在确定发生振荡（如ΔT=25℃）后，详细分析以下内容：
     • 平均场（u0, T0, p0）在自由表面的分布。
     • 扰动场（u’, T’, p’, 自由面变形h’）在自由表面的时空演化。
     • 平均流函数ψ0和平均温度T0在液层内的分布。
     • 扰动流函数ψ’和扰动温度T’在液层内的时空演化。
     • 瞬时总场（ψ, T）在液层内的形态。
     • 扰动速度、扰动温度、扰动压力在自由表面上的相位关系。
     • 计算各扰动量的幅值与对应平均量的比值，量化扰动的强度。

三、 主要研究结果 1. 稳态对流与失稳发生：当ΔT大于零但低于临界值时（如ΔT=7℃），流动为稳定的单胞结构：流体在表面张力驱动下从热壁流向冷壁，在底部区域则从冷壁流回热壁，温度沿垂直方向（z轴）为正梯度。当ΔT超过约9℃时，首先在靠近热壁的区域出现扰动胞。继续增大ΔT至25℃时，扰动胞扩展到整个液层，此温差被定义为临界温差ΔTc。对应的振荡频率为0.045144 Hz。

1. 热液波的确认与特征：

   • 传播方向：扰动速度u’、温度T’、压力p’和自由面变形h’在自由表面上的演化图显示，扰动波以行波形式传播，且传播方向与自由表面的主流平均速度方向相反，即向上游（热壁方向）传播，这符合Smith和Davis理论预测的热液波特征。
   • 波参数：在自由表面温度近似线性分布的区域（x = 10.4 mm 至 40.3 mm），计算得到的无量纲波数和相速度分别为2.417和0.0666，与Smith和Davis针对Pr→∞流体线性稳定性分析得到的结果（0.47和0.0622）接近。
   • 扰动强度：扰动非常微弱。在ΔT=25℃时，最大扰动速度、扰动温度和扰动压力与各自平均最大值的比值分别为0.51%、0.05%和0.38%。因此，尽管扰动存在，但瞬时总流场和温度分布图看起来与稳态场几乎相同（见图12），无法直接观察到明显的振荡多胞流结构。这解释了为何在实验中直接观测大Pr数流体的热液波较为困难。

2. 平均场与扰动场结构：

   • 平均流函数ψ0显示，在热壁附近区域出现了两个稳态胞。平均温度T0呈现典型的大Pr数热毛细对流温度分布，垂直方向为正梯度。
   • 扰动流函数ψ’和扰动温度T’的分布图（图11）揭示了热液波在液层内部的时空结构。扰动胞大致起源于x=40-45 mm区域，随时间向上游（逆主流方向）迁移，然后在靠近热壁区域下沉。流动胞的中心大致位于z=2H/3的高度线上。扰动温度模式由平均流速和扰动流速共同作用形成。

3. 相位关系与激发机制：

   • 这是本研究的关键发现。图13清晰地展示了自由表面上扰动温度T’、扰动压力p’和扰动速度u’之间的相位关系。扰动温度T’和扰动压力p’几乎同相，并且它们在波的传播方向上，领先于扰动速度u’约1/4个周期。
   • 具体表现为：T’和p’的最大值点、u’的零点以及扰动流函数ψ’的垂直零线（分隔正负扰动涡胞）在x坐标上位置一致。负ψ’胞的中心对应u’的最大值，正ψ’胞的中心对应u’的最小值。更重要的是，自由表面上的扰动速度u’是从较高的扰动温度区域流向较低的扰动温度区域。
   • 这一相位关系直接证明了驱动扰动波的动力来源于温度扰动。温度扰动导致表面张力分布发生变化，产生附加的切向应力（马兰戈尼应力），从而驱动了速度扰动。这种由温度扰动驱动、并与流动耦合的波，正是热液波。

4. 激发机制的物理解释：结合图13©中倾斜于自由表面的扰动等温线（虚线为负值，实线为正值）可以理解其物理图像。在自由表面附近的下方区域，扰动温度场存在负的温度梯度。在自由表面的高温扰动点，热量被带走，这部分热量由来自液层内部的垂直扰动速度所补充。这种热量的对流输运与传导之间的平衡，耦合了马兰戈尼效应和热毛细对流，形成并维持了热液波的传播。

四、 结论与意义 本研究通过数值模拟，成功揭示了大普朗特数流体在微重力浅槽中热毛细对流失稳后产生的热液波特性。主要结论如下： 1. 当侧壁温差超过临界值时，液层中会激发向上游（热壁方向）传播的行波，即热液波。 2. 该热液波由温度扰动驱动，扰动温度与扰动压力同相，并领先于扰动速度1/4周期。自由表面上的扰动速度从高温扰动区流向低温扰动区。 3. 对于Pr=105.6的硅油，热液波引起的流动和温度扰动幅值非常小（不到平均场的1%），导致瞬时总场与平均场视觉上几乎无差别，无法直接观察到振荡的多胞流结构。 4. 因此，在地面或空间的小尺度实验中，要判断热液波是否出现，不能依赖对总流场的直接目视观察，而需要采用更灵敏的检测手段，例如监测液层中某点速度或温度的频谱，或者利用光学方法高灵敏度地探测自由表面的振荡。

科学价值： 1. 机制阐释：研究通过详尽的数值结果，特别是首次清晰展示了热液波扰动场中速度、温度、压力三者的精确相位关系，为理解热液波的激发机制提供了直观、确凿的物理图像和证据，深化了对热毛细对流失稳物理本质的认识。 2. 方法指导：研究明确指出了对于大Pr数流体热液波实验观测的特殊性（扰动微弱，总场变化不明显），提出了有效的检测方法（频谱分析、光学测振），对后续的空间实验设计和数据分析具有重要的指导意义。 3. 区分不同失稳模式：研究结果暗示，在空间半浮区液桥实验中观察到的强烈振荡对流，其激发机制可能与浅液层中热液波的机制有所不同，这为进一步研究不同几何配置下的热毛细振荡提供了对比和思考。

五、 研究亮点 1. 研究对象特殊：专注于高普朗特数（Pr=105.6）流体，其行为与常用中低Pr数流体差异显著，扰动极其微弱，揭示了热液波在强粘性流体中的独特表现。 2. 揭示关键相位关系：研究的核心亮点在于细致刻画并明确了热液波中扰动温度、压力、速度三者之间精确的相位关系（T’与p’同相，并领先u’ 1/4周期），这为“温度扰动驱动波”的理论预测提供了直接的、定量的数值证据，是理解激发机制的关键。 3. 澄清观测误区：明确指出并解释了对于大Pr数流体，热液波的存在可能不会导致瞬时流场出现肉眼可见的多胞振荡结构，纠正了可能存在的观测预期误区，对实验研究具有重要警示和指导价值。 4. 数值方法可靠：采用了经过严格验证的分数步杂交有限元法，并与多个经典理论解和数值结果进行了对比，确保了研究结论的可信度。

六、 其他有价值内容 研究还简要讨论了热液波与贝纳德对流（浮力失稳）和皮尔逊不稳定性（纯马兰戈尼失稳）的区别。后两者涉及的是从静态到流动的转变，且通常需要垂直方向负的温度梯度（底部加热）。而本研究中的热液波发生在垂直方向为正温度梯度的稳定层结背景下，涉及的是从稳态对流到振荡对流的转变，其不稳定性机制与液层内热对流和热传导的平衡相关。这一对比有助于将热液波置于更广阔的流体失稳研究谱系中进行定位。