[论文解读] Nucleation of spherical shell-like interfaces by second gradient theory: numerical simulations

Francesco dell’Isola, Henri Gouin+1 more|ArXiv.org|Jun 10, 2009

nanoparticles nucleation surface interactions参考文献 28被引用 79

一句话总结

本文利用二阶梯度理论对流体相中的球形壳状界面进行建模，实现了表面张力、界面厚度及平衡半径随化学势变化的数值模拟。该理论预测了微观气泡的最小成核半径——即使在分子尺度下也成立，为经典Laplace-Tolman模型提供了一种物理解释一致的替代方案。

ABSTRACT

The theory of second gradient fluids (which are able to exert shear stresses also in equilibrium conditions) allows us: (i) to describe both the thermodynamical and the mechanical behavior of systems in which an interface is present; (ii) to express the surface tension and the radius of microscopic bubbles in terms of a functional of the chemical potential; (iii) to predict the existence of a (minimal) nucleation radius for bubbles. Moreover, the above theory supplies a 3D-continuum model which is endowed with sufficient structure to allow the construction of a 2D-shell-like continuum representing a consistent approximate 2D-model for the interface between phases. In this paper we use numerical simulations in the framework of second gradient theory to obtain explicit relationships for the surface quantities typical of 2D-models. In particular, for some of the most general two-parameter equations of state, it is possible to obtain the curves describing the relationship between the surface tension, the thickness, the surface mass density and the radius of the spherical interfaces between fluid phases of the same substance. These results allow us to predict the (minimal) nucleation radii for this class of equations of state.

研究动机与目标

• 基于二阶梯度理论构建一个三维连续体模型，以捕捉具有剪切应力能力的流体界面的力学与热力学行为。
• 通过从三维二阶梯度理论推导出的类壳二维模型，消除对Gibbs过剩量的依赖。
• 通过数值方法确定表面张力、界面厚度、表面质量密度及平衡半径对化学势的依赖关系。
• 预测球形气泡在平衡状态下的最小成核半径，尤其在经典理论失效的微观尺度下。
• 为临界点附近相变过程中的表面张力与曲率效应提供一个物理解释一致的框架。

提出的方法

• 构建一个包含质量密度二阶梯度非局部效应的三维二阶梯度流体界面模型。
• 利用力平衡定律推导出在球对称条件下控制平衡质量密度分布的偏微分方程。
• 应用“等效Laplace气泡”概念，基于三维场解定义表面张力与半径。
• 针对一般双参数状态方程（如范德瓦尔斯型）数值求解所得常微分方程，改变化学势参数。
• 从数值获得的质量密度分布中提取界面厚度、表面质量密度及曲率依赖的表面张力。
• 区分由自由能定义的热力学压强与由应力张量定义的机械压强，以提高物理一致性。

实验结果

研究问题

• RQ1 二阶梯度理论能否提供一个一致的三维描述，用于预测球形界面的最小成核半径？
• RQ2 在微观尺度下，表面张力如何随气泡半径变化？这种变化能否表示为化学势的函数？
• RQ3 在二阶梯度模型中，界面厚度、表面质量密度与平衡半径之间存在何种关系？
• RQ4 二阶梯度理论的预测与关于成核半径和表面张力的实验数据相比如何？
• RQ5 从二阶梯度理论推导出的类壳二维模型能否准确表征流体界面的力学与热力学行为？

主要发现

• 对于环己烷（C6H12），在 t ≈ 0.5（还原温度）时，预测的最小成核半径约为 9.8 Å，高于Fisher和Israelachvili的实验值，提示实验数据在临界点附近可能与经典Laplace假设不一致。
• 对于水，在 t ≈ 0.5 时，预测的最小成核半径为 R_m ≈ 8.3 Å，与Lamer和Pound的实验数据（8.0–8.9 Å）高度吻合。
• 对于CCl4，在 t ≈ 0.5 时，预测的平衡半径范围为12–14 Å，与Kumar的实验数据（12–17 Å）非常接近。
• 在最小成核半径处，界面厚度约为该半径的40%（比值 ≈ 0.4），表明在小尺度下界面占据气泡的显著部分。
• 计算得到的表面质量密度数量级与Alts和Hutter理论对水的预测一致，验证了该模型的物理合理性。
• 表面张力随气泡半径减小而降低，且其对化学势的函数依赖关系已明确推导并针对多种状态方程进行了数值验证。