[论文解读] Molecular dynamic simulation of water vapor interaction with blind pore of dead-end and saccate type

E. G. Nikonov, Miron Pavluš+1 more|arXiv (Cornell University)|Jul 28, 2017

Nanopore and Nanochannel Transport Studies被引用 4

一句话总结

本研究结合分子动力学（MD）模拟与宏观扩散模型，研究了死端型和囊状型盲孔中的水蒸气传输。通过从MD模拟中推导出扩散系数，并将其应用于三维扩散方程，作者证明该混合方法能准确捕捉干燥与湿润的动力学过程，微尺度与宏观尺度对水蒸气浓度随时间变化的预测结果高度一致。

ABSTRACT

One of the varieties of pores, often found in natural or artificial building materials, are the so-called blind pores of dead-end or saccate type. Three-dimensional model of such kind of pore has been developed in this work. This model has been used for simulation of water vapor interaction with individual pore by molecular dynamics in combination with the diffusion equation method. Special investigations have been done to find dependencies between thermostats implementations and conservation of thermodynamic and statistical values of water vapor - pore system. The two types of evolution of water-pore system have been investigated: drying and wetting of the pore. Full research of diffusion coefficient, diffusion velocity and other diffusion parameters has been made.

研究动机与目标

• 使用分子动力学方法模拟水蒸气与死端型和囊状型盲孔的相互作用。
• 研究恒温器参数对水蒸气-孔隙系统热力学稳定性的影响。
• 测定干燥与湿润过程中扩散系数和速度。
• 基于MD模拟得到的扩散系数，建立并验证宏观扩散模型。
• 将微尺度MD结果与基于MD-derived D值的宏观扩散模型预测结果进行比较，分析水蒸气浓度。

提出的方法

• 构建了尺寸为500×50×50 nm的三维盲孔模型，用于MD模拟。
• 采用Lennard-Jones势能，参数为σ=3.17 Å，ε=6.74×10⁻³ eV，描述水分子间的相互作用。
• 使用Velocity Verlet积分法，结合Berendsen恒温器（τB ∈ [0.5,2] ps），将温度维持在25 °C。
• 模拟运行时长为65.3 ns，时间步长Δt=0.016 ps，记录孔内1000–1800个H₂O分子，孔外1800–9000个H₂O分子。
• 从MD数据中计算干燥与湿润过程的扩散系数。
• 基于偏微分方程∂wv/∂t = D(∂²wv/∂x² + ∂²wv/∂y² + ∂²wv/∂z²)构建三维宏观扩散模型，边界条件考虑了孔口处的蒸气传输。

实验结果

研究问题

• RQ1 不同恒温器实现方式如何影响水蒸气-孔隙系统中热力学与统计性质的守恒？
• RQ2 在盲孔中，水蒸气在干燥与湿润过程中的时间依赖扩散系数为何值？
• RQ3 基于MD模拟得到的D值，微尺度MD结果与宏观扩散模型预测结果在水蒸气浓度上的差异如何？
• RQ4 孔隙几何形状与边界条件在水蒸气传输动力学中起何种作用？
• RQ5 混合MD-扩散模型能否准确描述盲孔中瞬态水蒸气浓度分布？

主要发现

• MD模拟表明，Berendsen恒温器在65.3 ns内有效维持系统温度为25 °C，且热力学性质稳定。
• 干燥过程的扩散系数约为1779.1 nm²/ps，湿润过程的扩散系数为536.33 nm²/ps。
• 宏观扩散模型通过变量分离法求解，可获得水蒸气浓度随时间变化的精确解析解。
• 宏观模型的空间平均水蒸气浓度与MD模拟结果高度一致，如图3所示。
• 该混合方法——即在宏观扩散模型中使用MD推导的D值——被验证为模拟水-孔隙相互作用的可靠方法。
• 该模型成功捕捉了干燥过程（初始浓度高）与湿润过程（初始浓度低）的动力学行为，且在不同时间尺度下浓度分布匹配良好。