涨落流体动力学与稀有事件方法揭示沸腾的介观物理

《Journal of Fluid Mechanics》：Fluctuating hydrodynamics and rare-event techniques unveil the mesoscale physics of boiling

　　

解读过热流体中相变的初始阶段及其后续演化仍是一项重大挑战，因为它需要桥接微观成核（nucleation）与宏观气泡动力学（bubble dynamics）和传热（heat transfer）。这一长期存在的差距阻碍了预测性沸腾模型的发展，并限制了新兴热技术的进展。在此，研究人员利用大规模的涨落流体动力学（fluctuating hydrodynamics, FHD）模拟，结合稀有事件技术（rare-event technique），以合理描述从成核到气泡生长的液-气转变。研究人员量化了成核气泡和流体动力学场的统计量、起始温度以及润湿性（wettability）对成核路径的非平凡影响。尽管在介观尺度（mesoscale）上运行，这些模拟再现了原子尺度研究中观察到的微观成核特征，包括沸腾起始温度和均相/异相成核路径之间的交叉对表面润湿性的依赖性。在较大尺度上，它们再现了实验观察到的从早期气泡接触面积的指数增长到后临界幂律标度的转变。这些发现建立了从涨落驱动成核到跨尺度涌现沸腾动力学之间的定量桥梁。

**论文解读文章：涨落流体动力学揭示沸腾的介观物理**

**研究背景与问题**
沸腾现象在核电站、微电子热管理等工程领域具有关键技术地位，但其初始相变成核阶段的物理机制长期未被完全理解。传统连续介质模型依赖经验阈值并预设已存在的相变胚体，无法预测成核发生的热力学条件；而分子动力学虽能捕捉原子级成核特征，却受限于计算规模，无法探索介于纳米与微米之间的介观尺度（mesoscale）窗口。这一尺度空白阻碍了从微观成核到宏观气泡动力学及传热之间的定量桥接，导致沸腾模型的预测能力受限。为此，研究人员尝试发展介观框架，将热涨落（thermal fluctuations）自洽地纳入流体动力学描述，以揭示非平衡相变的完整路径。

**研究内容与结论**
研究人员结合大规模涨落流体动力学（FHD）直接模拟与稀有事件（rare-event）分析，系统研究了从成核到气泡生长的完整液-气转变过程。结果表明：尽管模拟仅运行在介观尺度（约2.25 μm³），却能再现原子尺度研究中观察到的微观成核特征，包括亲水壁面附近液体纳米层（liquid nanolayer）的形成，以及表面润湿性（wettability）对沸腾起始温度和均相/异相成核路径交叉的调控。在更大尺度上，模拟定量复现了实验观察到的早期气泡接触面积指数分布向后临界幂律标度（power-law scaling）的转变，临界指数与宏观测量值吻合。这些发现建立了从涨落驱动成核到跨尺度涌现沸腾动力学之间的定量桥梁，为多尺度沸腾模拟和热系统设计提供了原理性策略。论文发表在《Journal of Fluid Mechanics》。

**关键技术方法**
（1）**涨落流体动力学（FHD）**：基于朗道-利夫希茨-纳维-斯托克斯（Landau-Lifshitz–Navier–Stokes）方程与扩散界面热力学（diffuse-interface thermodynamics）耦合，在质量、动量和能量守恒方程中引入高斯随机涨落通量，其振幅由涨落-耗散定理（fluctuation–dissipation theorem）确定，以确保与爱因斯坦-玻尔兹曼分布的一致性。
（2）**最小自由能路径与字符串方法**：基于范德瓦尔斯平方梯度近似构建朗道自由能泛函，通过字符串方法（string method）计算不同润湿性和过热度下的最小自由能路径（MFEP），识别临界核结构及成核能垒。
（3）**大规模数值模拟**：采用交错网格有限差分与二阶龙格-库塔时间积分，在3D立方域（L=2.25 μm，网格数500³）中求解FHD方程，壁面施加恒定热流和润湿性边界条件，顶部保持饱和压力与温度。参考水物性（T_c=647 K, p_c=22 MPa），将结果归一化为约化变量。

**研究结果**

**3. Results**

**模拟快照与成核机制**
图2展示了不同接触角下的沸腾快照。对于完全浸润壁面（φ=0°），气泡在近壁液体体内成核，而非直接接触壁面，呈现压缩液体纳米层结构，该现象与分子动力学（MD）模拟一致。对于中性壁面（φ=90°），成核直接发生在壁面。中间接触角（30°–75°）时，壁面成核气泡数随浸润性减弱而增多，系统逐渐从均相成核过渡到异相成核。

**壁面温度、成核时间与气泡统计**
图3(a)显示平均约化壁温?Θ_wall?随时间演化：达到最大值（onset，沸腾起始）后因潜热提取而下降，再因蒸汽膜形成而回升。图3(b)表明：起始温度Θ_ons随接触角减小而升高，但当φ≤30°时趋于平台（Θ_ons≈0.95）；起始时间t_ons在φ≤60°时也出现平台，表明强亲水性下均相成核主导，与经典成核理论（CNT）的单调趋势不同。图3(c)给出气泡数目随距壁面距离z的分布：φ=90°时气泡紧贴壁面；φ=45°时分布展宽；φ=0°时成核位置更远离壁面，总气泡数减少。图3(d)通过沃罗诺伊镶嵌（Voronoi tessellation）分析壁面气泡空间分布，发现归一化面积α的概率分布与随机泊松点过程（RPPP）完全吻合，方差接近理论值0.286，表明成核事件空间均匀且不相关。图3(e)显示起始温度随热流密度Q的增加而升高（因加热越快，稀有事件可用时间越短），起始时间则随Q增加而缩短，体现了沸腾的非平衡本质。

**场统计与涨落校正**
图4(a)展示了壁面温度与密度概率分布的时间演化：初始为高斯分布，后因蒸汽形成而出现左侧重尾，最终趋于双峰再恢复高斯。图4(b)将平均壁温与平均壁密度的关系与状态方程预测对比，发现考虑涨落校正（自洽二阶泰勒展开）后，理论曲线与FHD模拟结果吻合，证明了宏观层面考虑涨落的必要性。

**最小自由能路径与成核机制过渡**
图5给出了MFEP对应的密度剖面：φ=0°时临界核出现在近壁液体中，存在纳米层；φ=60°时仍有纳米层；φ=75°和90°时成核直接发生在壁面。图6为归一化成核能垒ΔΩ^?随过热度Θ和接触角φ的等值线图。与CNT相比，扩散界面模型自然包含曲率依赖性表面张力，能垒更低且更准。通过比较均相/异相成核能垒，识别出一条过渡边界（白色虚线）：在FHD模拟的典型过热度（Θ?0.7）下，φ?60°时均相成核主导，这与图3(b)中t_ons的平台行为一致。

**气泡接触面积统计**
图7展示了气泡-壁接触面积PDF的时间演化：成核阶段呈指数分布（p(R)~exp(-ΔΩ/k_BT)），随后进入幂律阶段（p~A^-γ），指数γ随时间变化，从初期的-2.81~-2.09演变为后临界平均-1.87，与实验报告值（-1.85至-1.50）吻合。后期幂律区出现凸起，对应大蒸汽簇形成，标志着渗流驱动的近临界团簇凝聚。

**4. Discussion**
研究人员指出，FHD结合扩散界面描述提供了捕获完整成核路径的一致框架，桥接了介观统计热力学与流体动力学。模拟结果在微观上再现了MD中的纳米层现象与润湿性效应，在宏观上复现了接触面积的指数-幂律转变，从而建立了跨尺度定量桥接。所计算的起始过热度和气泡空间统计（可描述为RPPP）可直接输入宏观求解器，开启将涨落驱动成核与连续尺度沸腾模拟耦合的路径。当前工作聚焦于近临界点（约化温度T≈0.95–0.96），该条件在加压水反应堆等工程中具有相关性，且框架可推广至更低温区（需结合稀有事件技术）。总体而言，这一工作将沸腾视为热激活且由流体动力学介导的非平衡过程，推动了相变非平衡统计力学理论的发展。

**研究结论翻译**
总体而言，当前结果表明，FHD结合扩散界面描述提供了一个一致框架，用于捕获流体的完整成核路径，在统一描述中桥接了介观统计热力学与流体动力学。尽管本研究聚焦于临界点附近的区域（此方法自然适用），但这些热力学条件并非纯学术兴趣，例如在约化温度T≈0.9和约化压力p≈0.8的加压水反应堆以及低温高压系统中具有工程相关性。更广泛而言，当辅以稀有事件技术以触及原本不可及的活化区间时，该框架预期在远离临界点时仍然有效。尽管此处针对理想范德瓦尔斯型流体进行演示，该框架通过采用精确状态方程可自然扩展至更真实的热力学描述。结合现代高性能计算能力，该方法有望在介观尺度上实现直接模拟，并逐步捕获从成核到宏观气泡动力学的完整沸腾物理，为相变现象的建模开辟新视角。从基础角度看，这些发现也激励非平衡统计力学的理论发展——通过将沸腾视作由流体动力学介导的涨落激活过程，本文指出了超越经典平衡热力学极限理解相变的新途径。