2025-11-20T16:49:15.210607

Non-equilibrium evaporation of Lennard-Jones fluids: Enskog-Vlasov theory and Hertz-Knudsen model

Li, Gibelli, Zhang

Enskog-Vlasov equation is currently the most sophisticated kinetic model for describing non-equilibrium evaporative flows. While it enables more efficient simulations than the molecular dynamics (MD) methods, its accuracy in reproducing the flow properties of real fluids is limited by both the assumptions underlying the Vlasov forcing term and the approximation introduced by the Enskog collision term for short-range molecular interactions. To address this limitation, this work proposes a molecular kinetic model specifically designed for real fluids, with the Lennard-Jones fluids as an example. The model is first applied to evaluate the equilibrium characteristics of a liquid-vapour system, including the liquid-vapour coexistence curve, transport coefficients, vapour pressure, and surface tension coefficient. The results show excellent agreement with the MD simulation and experimental data. Furthermore, the model is used to investigate non-equilibrium evaporation, with a particular focus on the velocity distribution function adjacent to the liquid-vapour interface. The results confirm that deviations from the Maxwellian distribution persist in the vapour region, indicating limitations of the classical Hertz-Knudsen relation under pronounced non-equilibrium conditions. This work represents a critical step towards the development of an accurate and efficient computational framework for modelling non-equilibrium liquid-vapour flows for real fluids, with direct relevance to practical applications such as flow cooling.

academic

Non-equilibrium evaporation of Lennard-Jones fluids: Enskog-Vlasov theory and Hertz-Knudsen model

基本信息

论文ID: 2508.15785
标题: Non-equilibrium evaporation of Lennard-Jones fluids: Enskog-Vlasov theory and Hertz-Knudsen model
作者: Shaokang Li, Livio Gibelli, Yonghao Zhang
分类: cond-mat.soft physics.flu-dyn
发表时间: 2025年10月14日 (arXiv v2)
论文链接: https://arxiv.org/abs/2508.15785

摘要

Enskog-Vlasov方程是目前描述非平衡蒸发流动最精密的动力学模型。虽然它比分子动力学(MD)方法更高效，但其在再现真实流体流动性质方面的准确性受到Vlasov强迫项假设和Enskog碰撞项对短程分子相互作用近似的限制。为解决这一局限性，本工作提出了专门针对真实流体设计的分子动力学模型，以Lennard-Jones流体为例。该模型首先用于评估液-气系统的平衡特性，包括液-气共存曲线、输运系数、蒸气压和表面张力系数。结果显示与MD模拟和实验数据具有优异的一致性。此外，该模型被用于研究非平衡蒸发，特别关注液-气界面附近的速度分布函数。结果证实在气相区域持续存在偏离麦克斯韦分布的现象，表明经典Hertz-Knudsen关系在显著非平衡条件下的局限性。

研究背景与动机

1. 研究问题

本研究要解决的核心问题是在纳米尺度蒸发过程中，现有动力学模型无法准确描述真实流体（如Lennard-Jones流体）的非平衡蒸发行为。

2. 问题重要性

纳米尺度蒸发在多个前沿技术领域具有关键作用：

纳米多孔膜中的物质分离效率提升
高性能电子设备的高效散热
流动冷却等实际应用

3. 现有方法局限性

Hertz-Knudsen关系：假设气相区域局部平衡，忽略了Knudsen层中的动力学结构
Enskog-Vlasov方程：使用Sutherland势能模型，与真实Lennard-Jones势能存在显著差异
分子动力学方法：计算成本过高，难以应用于工程实际

4. 研究动机

发展一个统一的框架，能够：

准确描述液相、气相及其界面
消除对任意蒸发系数的依赖
在保持计算效率的同时提高预测精度

核心贡献

提出了改进的Enskog-Vlasov动力学模型：专门针对Lennard-Jones流体设计，通过修正对相关函数准确再现真实流体的状态方程
建立了精确的热力学参数校准方法：通过同时匹配临界点和三相点参数，确保模型与真实流体性质的一致性
验证了模型的平衡性质预测能力：在液-气共存曲线、输运系数、蒸气压和表面张力等方面与实验数据高度吻合
揭示了非平衡蒸发中的速度分布特征：证实了气相区域速度分布偏离麦克斯韦分布，为Hertz-Knudsen关系的局限性提供了理论依据

方法详解

任务定义

研究Lennard-Jones流体在非平衡条件下的蒸发过程，重点分析液-气界面附近的分子速度分布特征，并评估经典Hertz-Knudsen关系的适用性。

模型架构

1. 基础动力学方程

模型基于修正的Enskog-Vlasov方程：

$\frac{\partial f}{\partial t} + \boldsymbol{\xi} \cdot \frac{\partial f}{\partial \boldsymbol{x}} + \frac{\boldsymbol{F}}{m} \cdot \frac{\partial f}{\partial \boldsymbol{\xi}} = \Omega$

其中：

$f(\boldsymbol{x}, \boldsymbol{\xi}, t)$ ：分子速度分布函数
$\boldsymbol{F}$ ：Sutherland势能产生的自洽力场
$\Omega$ ：简化的碰撞项

2. 简化碰撞项

碰撞项采用Taylor展开近似： $\Omega = J_S + J_e$

其中 $J_S$ 为Shakhov模型项， $J_e$ 为超额项，确保正确的Prandtl数和体粘性。

3. 状态方程修正

采用广义van der Waals形式： $p = nk_BT(1 + \rho b\chi) - an^2$

通过修正对相关函数 $\chi$ 来匹配Lennard-Jones流体的热力学性质。

技术创新点

1. 对相关函数的经验修正

引入经验函数 $\Phi$ ： $\Phi = \rho b + a_1(\rho b)^2 + a_2(\rho b)^3 + a_3(\rho b)^4 + a_4(\rho b)^5$

其中 $\chi = \Phi' = d\Phi/d(\rho b)$

2. 多约束条件校准

同时满足：

临界点条件： $\left(\frac{\partial p}{\partial n}\right)_{T=T_c,n=n_c} = 0$ ， $\left(\frac{\partial^2 p}{\partial n^2}\right)_{T=T_c,n=n_c} = 0$
三相点条件：压力和化学势相等

参考流体：氩气(Argon)作为Lennard-Jones流体的典型代表
实验数据：液-气共存曲线、输运系数、蒸气压、表面张力的实验测量值
MD模拟数据：用于对比验证

评价指标

平衡性质：
- 液-气共存密度
- 剪切粘性和热导率
- 饱和蒸气压
- 表面张力系数
非平衡特征：
- 速度分布函数偏离麦克斯韦分布的程度
- Knudsen层结构

对比方法

原始Enskog-Vlasov方程
实验数据
分子动力学模拟结果

实现细节

计算域长度： $L$ ，采用周期性边界条件（平衡）或吸收边界条件（非平衡）
温度范围：从三相点到临界点
氩气参数： $a = 162.5 \text{ J·l·mol}^{-2}$ ， $n_c = 13.41 \text{ mol·l}^{-1}$ ， $T_c = 150.69 \text{ K}$