Molecular dynamic simulations of selective self-diffusion of CH4/CO2/H2O/N2 in coal

Abstract

The self-diffusion coefficients (D) of CH4/CO2/H2O/N2 at a relatively broad range of temperatures(298.15∼ 458.15K)and pressures (1∼6MPa) under the NPT, NPH, NVE, and NVT ensembles were obtained after the calculations of molecular mechanics(MM), annealing kinetics(AK), giant canonical Monte Carlo(GCMC), and molecular dynamics (MD) based on Wiser bituminous coal model (WM). The Ds of the adsorbates at the saturated adsorption configurations are D CH4 < D N2 (NVT)≈D N2 (NPH)≈D N2 (NPT) (T<418K) and D N2 (NVE) is remarkable higher than other ensembles when T>418K. The average swelling ratios manifest as H2O (14.7∼35.18%)>CO2 (13.38∼32.25%)>CH4 (15.35∼23.71%)> N2 (11.47∼22.14%) (NPH, 1∼6MPa). There exits differences in D, swelling ratios and E among various ensembles, indicating that the selection of ensembles has an important influence on the MD calculations for self-diffusion coefficients.