Residual trapping is a critical mechanism influencing the efficiency of Underground Hydrogen Storage (UHS). This study investigates the underlying processes of residual trapping by bypassing, through bifurcating geometries, focusing on how geometrical parameters and flow characteristics affect the trapping process. We develop a dynamic simulation framework based on the lattice Boltzmann method (LBM) to simulate full drainage/imbibition cycles. Various geometries, based on the pore doublet model, were investigated and supported by theoretical analysis. In addition, trapping behavior of hydrogen was compared to that of CO₂ and CH₄. It is found that the channel width ratio, specially across the local bifurcating geometries, and the roundness of the grains, are among the key factors which control hydrogen trapping. Results indicate that the suited reservoirs for underground hydrogen storage have narrower channel-size ratios and smoother edges at micro-scale. Operational conditions also play a significant role. Lower flow rates enhance bypassing, which increases trapping.

Large-scale geological storages of hydrogen (H₂) and carbon dioxide (CO₂) in saline aquifers present feasible options for a sustainable energy future. We compared the plume migration of CO₂ and H₂ in aquifers using the FluidFlower benchmark, incorporating the state-of-the-art thermophysical and petrophysical properties. The H₂ plume, with its higher buoyancy and mobility compared to CO₂, remains predominantly in the gas phase due to its lower solubility, increasing the chances of escaping through fractures or migration to distant regions. This additionally leads to a higher pressurized reservoir, which, along with higher buoyancy, increases the chance of caprock penetration. Dissolution trapping of CO₂ into brine increases over time due to its fingering, while H₂ does not show fingering. Our findings show that while geological carbon storage (GCS) benefits significantly from all structural, dissolution, and residual trapping, underground hydrogen storage (UHS) relies mainly on structural trapping, making the integrity of sealing elements of the system a key factor in its performance.

Data for several key thermodynamic and transport properties needed for technologies using hydrogen (H₂), such as underground H₂ storage and H₂O electrolysis are scarce or completely missing. Force field-based Molecular Dynamics (MD) and Continuous Fractional Component Monte Carlo (CFCMC) simulations are carried out in this work to cover this gap. Extensive new data sets are provided for (a) interfacial tensions of H₂ gas in contact with aqueous NaCl solutions for temperatures of (298 to 523) K, pressures of (1 to 600) bar, and molalities of (0 to 6) mol NaCl/kg H₂O, (b) self-diffusivities of infinitely diluted H₂ in aqueous NaCl solutions for temperatures of (298 to 723) K, pressures of (1 to 1000) bar, and molalities of (0 to 6) mol NaCl/kg H₂O, and (c) solubilities of H₂ in aqueous NaCl solutions for temperatures of (298 to 363) K, pressures of (1 to 1000) bar, and molalities of (0 to 6) mol NaCl/kg H₂O. The force fields used are the TIP4P/2005 for H₂O, the Madrid-2019 and the Madrid-Transport for NaCl, and the Vrabec and Marx for H₂. Excellent agreement between the simulation results and available experimental data is found with average deviations lower than 10%.

Underground Hydrogen Storage (UHS) is an attractive technology for large-scale (TWh) renewable energy storage. To ensure the safety and efficiency of the UHS, it is crucial to quantify the H₂ interactions with the reservoir fluids and rocks across scales, including the micro scale. This paper reports the experimental measurements of advancing and receding contact angles for different channel widths for a H₂/water system at P = 10 bar and T = 20 °C using a microfluidic chip. To analyse the characteristics of the H₂ flow in straight pore throats, the network is designed such that it holds several straight channels. More specifically, the width of the microchannels range between 50 μm and 130 μm. For the drainage experiments, H₂ is injected into a fully water saturated system, while for the imbibition tests, water is injected into a fully H₂-saturated system. For both scenarios, high-resolution images are captured starting the introduction of the new phase into the system, allowing for fully-dynamic transport analyses. For better insights, N₂/water and CO₂/water flows were also analysed and compared with H₂/water. Results indicate strong water-wet conditions with H₂/water advancing and receding contact angles of, respectively, 13°–39°, and 6°–23°. It was found that the contact angles decrease with increasing channel widths. The receding contact angle measured in the 50 μm channel agrees well with the results presented in the literature by conducting a core-flood test for a sandstone rock. Furthermore, the N₂/water and CO₂/water systems showed similar characteristics as the H₂/water system. In addition to the important characterization of the dynamic wettability, the results are also crucially important for accurate construction of pore-scale simulators.