Modeling of Thermal Transport Properties in Metallic and Oxide Fuels

Thermal conductivity is a critical fuel performance property not only for current UO2 oxide fuel based light water reactors but also important for next-generation fast reactors that use U-Zr based metallic fuels. In this work, the thermal transport properties of both UO2 based oxide fuels and U-Zr based metallic fuels have been studied. At first, molecular dynamics (MD) simulations were conducted to study the effect of dispersed Xe fission gas atoms on the UO2 thermal conductivity. Numerous studies have demonstrated that xenon (Xe) fission gas plays a major role on fuel thermal conductivity degradation. Even a very low Xe concentration can cause significant thermal conductivity reduction. In this work, the effect of dispersed Xe gas atoms on UO2 thermal conductivity were studied using three different interatomic potentials. It is found that although these potentials result in significant discrepancies in the absolute thermal conductivity values, their normalized values are very similar at a wide range of temperatures and Xe concentrations. By integrating this unified effect into the experimentally measured thermal conductivities, a new analytical model is developed to predict the realistic thermal conductivities of UO2 at different dispersed Xe concentrations and temperatures. Using this new model, the critical Xe concentration that offsets the grain boundary Kapitza resistance effect on the thermal conductivity in a high burnup structure is studied. Next, the mechanisms on how Xe gas bubbles affect the UO2 thermal conductivity have been studied using MD. At a fixed total porosity, the effective thermal conductivity of the bubble-containing UO2 increases with Xe cluster size, then reaches a nearly saturated value at a cluster radius of 0.6 nm, demonstrating that dispersed Xe atoms result in a lower thermal conductivity than clustering them into bubbles. In comparison with empty voids of the same size, Xe-filled bubbles lead to a lower thermal conductivity when the number ratio of Xe atoms to uranium vacancies (Xe:VU ratio) in bubbles is high. Detailed atomic-level analysis shows that the pressure-induced distortion of atoms at bubble surface causes additional phonon scattering and thus further reduces the thermal conductivity. For metallic fuels, temperature gradient and irradiation induced constituent redistribution in U-Zr based fuels cause the variation in fuel composition and the formation of different phases that have different physical properties such as thermal conductivity. In this work, a semi-empirical model is developed to predict the thermal conductivities of U-Zr alloys for the complete composition range and a wide range of temperatures. The model considers the effects of (a) scattering by defects, (b) electron-phonon scattering, and (c) electron-electron scattering. The electronic thermal resistivity models for the two pure components are empirically determined by fitting to the experimental data. A new mixing rule is proposed to predict the average thermal conductivity in U-Zr alloys based on their nominal composition. The thermal conductivity predictions by the new model show good agreement with many available experimental data. In comparison with previous models, the new model has further improvement, in particular for high-U alloys that are relevant to reactor fuel compositions and at the low-temperature regime for the high-Zr alloys. The average thermal conductivity model for the binary U-Zr fuel is also coupled with finite element-based mesoscale modeling technique to calculate the effective thermal conductivities of the U-Zr heterogeneous microstructures. For a U-10wt.%Zr (U-10Zr) fuel at temperatures below the ɑ phase transition temperature, the dominant microstructures are lamellar δ-UZr2 and ɑ-U. Using the mesoscale modeling, the phase boundary thermal resistance R (Kapitza resistance) between δ-UZr2 and ɑ-U has been determined at different temperatures, which shows a T-3 dependence in the temperature range between 300K and 800K. Besides, the Kapitza resistance exhibits a strong dependence on the aspect ratio of the δ-UZr2 phase in the alloying system. An analytical model is therefore developed to correlate the temperature effect and the aspect ratio effect on the Kapitza resistance. Combining the mesoscale modeling with the newly developed Kapitza resistance model, the effective thermal conductivities of many arbitrary δ-UZr2 + ɑ-U heterogeneous systems can be estimated.