Division of Mechanical and Space Engineering, Hokkaido University, Sapporo, Hokkaido, Japan.
PLoS One. 2021 Mar 16;16(3):e0248660. doi: 10.1371/journal.pone.0248660. eCollection 2021.
The kinetic boundary condition (KBC) represents the evaporation or condensation of molecules at the vapor-liquid interface for molecular gas dynamics (MGD). When constructing the KBC, it is necessary to classify molecular motions into evaporation, condensation, and reflection in molecular-scale simulation methods. Recently, a method that involves setting the vapor boundary and liquid boundary has been used for classifying molecules. The position of the vapor boundary is related to the position where the KBC is applied in MGD analyses, whereas that of the liquid boundary has not been uniquely determined. Therefore, in this study, we conducted molecular dynamics simulations to discuss the position of the liquid boundary for the construction of KBCs. We obtained some variables that characterize molecular motions such as the positions that the molecules reached and the time they stayed in the vicinity of the interface. Based on the characteristics of the molecules found from these variables, we investigated the valid position of the liquid boundary. We also conducted an investigation on the relationship between the condensation coefficient and the molecular incident velocity from the vapor phase to the liquid phase. The dependence of the condensation coefficient on the incident velocity of molecules was confirmed, and the value of the condensation coefficient becomes small in the low-incident-velocity range. Furthermore, we found that the condensation coefficient in the non-equilibrium state shows almost the same value as that in the equilibrium state, although the corresponding velocity distribution functions of the incident velocity significantly differ from each other.
动力学边界条件 (KBC) 代表了分子在气-液界面处的蒸发或凝结,这对于分子气体动力学 (MGD) 来说非常重要。在构建 KBC 时,有必要将分子运动分为蒸发、凝结和反射,这需要在分子尺度模拟方法中进行分类。最近,已经有一种方法涉及设置蒸汽边界和液体边界来对分子进行分类。蒸汽边界的位置与 MGD 分析中应用 KBC 的位置有关,而液体边界的位置尚未被唯一确定。因此,在本研究中,我们进行了分子动力学模拟,以讨论 KBC 构建中液体边界的位置。我们获得了一些可以描述分子运动的变量,例如分子到达的位置和它们在界面附近停留的时间。基于这些变量中发现的分子特征,我们研究了液体边界的有效位置。我们还研究了从气相到液相的分子入射速度与凝结系数之间的关系。证实了凝结系数与分子入射速度的依赖性,并且在低入射速度范围内凝结系数变小。此外,我们发现非平衡态的凝结系数几乎与平衡态的凝结系数相同,尽管对应的入射速度分布函数存在显著差异。