Molecular Perspectives of Interfacial Properties in the Water+Hydrogen System in Contact with Silica or Kerogen

TL;DR

This study uses molecular dynamics simulations to explore interfacial properties of water, hydrogen, silica, and kerogen systems relevant to underground hydrogen storage, revealing how temperature and pressure influence interfacial tension and contact angles.

Contribution

It provides detailed molecular insights into interfacial behaviors of H2O+H2 with silica and kerogen, including effects of temperature and pressure on interfacial tension and contact angles, validated by experimental and theoretical data.

Findings

Interfacial tensions decrease with pressure and temperature, but increase at high T and P.

Contact angles vary with pressure and temperature, differing between silica and kerogen.

H2 enrichment at interfaces is minimal, but H2O shows strong adsorption on silica.

Abstract

Interfacial behaviours in multiphase systems containing H2 are crucial to underground H2 storage but are not well understood. Molecular dynamics simulations were conducted to study interfacial properties of the H2O+H2 and H2O+H2+silica/kerogen systems over a wide range of temperatures (298 - 523 K) and pressures (1 - 160 MPa). The combination of the H2 model with the INTERFACE force field and TIP4P/2005 H2O model can accurately predict the interfacial tensions (IFTs) from the experiment. The IFTs from simulations are also in good agreement with those from the density gradient theory coupled to the PC-SAFT equation of state. Generally, the IFTs decrease with pressure and temperature. However, at relatively high temperatures and pressures, the IFTs increase with pressure. The opposite pressure effect on IFTs can be explained by the inversion of the sign of the relative adsorption of H2. The enrichment of H2 in the interfacial regions was observed in density profiles. Meanwhile, the behaviours of contact angles (CAs) in the H2O+H2+silica system are noticeably different from those in the H2O+H2+kerogen system. The H2O CAs for the H2O+H2+silica and H2O+H2+kerogen systems increase with pressure and decrease with temperature. However, the effect of temperature and pressure on these CAs is less pronounced for the H2O+H2+silica system at low temperatures. The behaviours of CAs were understood based on the variations of IFTs in the H2O+H2 system (fluid-fluid interaction) and adhesion tensions (fluid-solid interaction). Furthermore, the analysis of the atomic density profiles shows that the presence of H2 in between the H2O droplet and the silica/kerogen surface is almost negligible. Nevertheless, the adsorption of H2O on the silica surface outside the H2O droplet is strong, while less H2O adsorption is seen on the kerogen surface.