Molecular dynamics simulation of hydrogen diffusion into brine: Implications for underground hydrogen storage

Underground hydrogen storage has been the subject of extensive research recently. Hydrogen diffusion coefficient into brine critically affects hydrogen loss. However, the available data pertain to hydrogen diffusivity in pure water or NaCl solutions, with rare interpretations of temperature and salinity effects. In this paper, molecular dynamics simulations at realistic conditions were conducted to calculate these coefficients in various brines, particularly those containing divalent salts. The results indicate lower diffusivities at higher salinities and lower temperatures, the value being 7.29 x 10-9 m2/s at 323 K, increasing to 10.2 x 10-9 m2/s at 353 K for 1 molal NaCl solution. The diffusion coefficient decreases up to 38 % as the salinity increases from 1 to 5 molal. Besides, hydrogen diffusivity in MgCl2 solution is up to 60 % smaller than in NaCl solutions at the same molality. Radial distribution functions and hydrogen bond analyses confirm the results. This research provides new data for estimating hydrogen loss in aquifers or water-saturated caprocks.