Three-dimensional multi-physics modeling of hydrogen and hydride distribution in zirconium alloy cladding

Localized phenomena within the reactor core, specifically those associated with the nuclear fuel, require high-fidelity simulations to enable accurate physics predictions. One example is the zirconium cladding, which absorbs hydrogen from the light water coolant during normal reactor operation. Absorbed in the cladding, this hydrogen is in solid solution and its distribution is sensitive to temperature and concentration gradients. At high enough concentrations, the hydrogen will precipitate as a hydride. Thus, the hydrogen distribution as a hydride precipitate in cladding has been identified as an important safety concern, and a possible ersatz for validating reactor simulation code temperature models. This study reports development efforts of using high-fidelity multi-physics codes to model temperature, hydrogen, and hydride distributions in three dimensions under realistic operating conditions. The Consortium for the Advanced Simulation of Light Water Reactors multi-physics code, Tiamat, is used to model selected sub-assemblies. Then, a single fuel pin is selected from the sub-assembly and modeled as a three-dimensional BISON problem. The outer cladding temperatures from the Tiamat calculation are used as boundary conditions for the BISON problem in order to obtain hydrogen and hydride distributions. Areas of interest for hydride precipitation include locations along the fuel rod experiencing highest temperatures with significant spatial variation, particularly in the vicinity of the spacer grids and mixing vanes.