Numerical Study of the Hydride Embrittlement in Zirconium Alloy using XFEM

A numerical model for hydride embrittlement in Zirconium alloy (Zr-2.5Nb) is developed utilizing the extended finite element method (XFEM). Hydride embrittlement reduces the ductility and failure time of a metal/alloy. During hydride embrittlement, stress-directed hydrogen diffusion, metal-hydride phase transformation, mechanical deformation, and hydride precipitation occur simultaneously. The present model incorporates all these processes and is able to predict the hydrogen concentration and the hydride fraction distribution under any externally applied stress field. In this work, both the steady and transient hydrogen diffusion cases are evaluated. Further, the XFEM is utilized to develop a model of hydride embrittlement in the presence of a crack. The first step of the hydride embrittlement process is the diffusion of hydrogen. According to Fick's law of diffusion, hydrogen diffusion is directly dependent on hydrostatic stresses and hydrogen concentration gradient under external stresses. The next step is the hydride precipitation in hydride embrittlement, where the expansion of material takes place that changes the hydrostatic stress field. Thus, studying the effect of precipitation of hydride on hydrostatic stresses is essential. Moreover, the process of hydride embrittlement is highly influenced by residual stresses in the structure. Hence, the effect of residual stress present in the zirconium alloy pressure tube (PT) is also evaluated. The results indicate that the residual tensile stresses contribute to the growth of hydride, which will reduce the material failure time.