A COMSOL framework for predicting hydrogen embrittlement -- Part I: coupled hydrogen transport

Hydrogen threatens the structural integrity of metals and thus predicting hydrogen-material interactions is key to unlocking the role of hydrogen in the energy transition. Quantifying the interplay between material deformation and hydrogen diffusion ahead of cracks and other stress concentrators is key to the prediction and prevention of hydrogen-assisted failures. In this work, a generalised theoretical and computational framework is presented that for the first time encompasses: (i) stress-assisted diffusion, (ii) hydrogen trapping due to multiple trap types, rigorously accounting for the rate of creation of dislocation trap sites, (iii) hydrogen transport through dislocations, (iv) equilibrium (Oriani) and non-equilibrium (McNabb-Foster) trapping kinetics, (v) hydrogen-induced softening, and (vi) hydrogen uptake, considering the role of hydrostatic stresses and local electrochemistry. Particular emphasis is placed on the numerical implementation in COMSOL Multiphysics, releasing the relevant models and discussing stability, discretisation and solver details. Each of the elements of the framework is independently benchmarked against results from the literature and implications for the prediction of hydrogen-assisted fractures are discussed. The second part of this work (Part II) shows how these crack tip predictions can be combined with crack growth simulations.