Modeling thermochemistry of fuel and coupling to fuel performance codes

A thermochemical library, Thermochimica, utilizes new algorithms and solvers designed for calculating equilibria of multicomponent and multiphase systems and has been coupled to the nuclear fuel performance code BISON. Thermochimica utilizes a generalized thermochemical database for urania-based fuels that has been under development for a number of years and continues to be expanded. One of the key motivations for use of thermodynamic models is that they can better represent the physics of oxygen migration. The coupling of the BISON code with Thermochimica and examples of its use in representing fuel phenomena are described. In the current implementation in BISON, fluxes are driven by the elemental concentration gradient. However, equilibrium is characterized by the absence of spatial variations of chemical potentials of the system species/components. The driving force to equilibrium (species flux) is proportional to the deviation from equilibrium, i.e., the gradient in chemical potential. It is therefore necessary to integrate chemical potential driven diffusion in the BISON representation of oxygen transport. The diffusion kernels based on thermodynamic models have been developed to simulate multicomponent transport in light water reactor (LWR) materials. The Thermochimica library was used to calculate various thermodynamics properties needed for transport calculations. The library also can calculate properties that are specific to various transport mechanisms for different fuel materials. For example, it can calculate defect site fractions among the sublattices in oxide fuels as a function of temperature, burnup, and stoichiometry. The transport and thermodynamic models are being integrated with other simulation efforts such as calculation of mobility functions that couple fluxes and chemical gradients.