SEARCH · Engineering Papers

Results for “thermal system modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 361 records · Page 20

All sources NASA NTRS DOE OSTI

DOE OSTI2025

A hybrid surrogate modeling framework for the Digital Twin of a Fluoride-salt-cooled High-temperature Reactor (FHR)

While nuclear energy is a non-greenhouse-gas emitting energy source, expensive operational costs due to the high-level of safety requirements decreases their competitiveness in the sustainable energy market. Advanced reactor concepts paired with Digital Twins aim to increase the commercialization gains of nuclear energy by reducing operational costs, increasing reactor reliability and enhancing power generation. To support Digital Twin tasks such as real-time autonomous control, proactive maintenance monitoring or optimizing power demand operations, a fast and accurate virtual representation of the Nuclear Power Plant (NPP) is required. The computational cost of high-fidelity, physics-based models are unsuitable for real-time analysis or scalability. Here, in this work, a hybrid surrogate modeling framework is developed fora Fluoride-salt-cooled High-temperature Reactor (FHR) that leverages physics-inspired models for key reactor components and uses data-driven methods for rapid system state space prediction. The Xenon reactivity feedback model is integrated to inform the surrogate model about the reactor core and the homologous pump theory model is the basis for representing pump degradation. Using a detailed, two dimensional thermal hydraulics model to generate data on the FHR, we train a network of Vectorized Autoregressive Moving-Average with eXogenous input (VARMAX) models to predict the remaining state values. The result is a surrogate model that provides a detailed reactor state representation of 41 system states and a pump degradation analysis. The framework is applied to Load Follows profiles, yielding high accuracy and a speedup that is more than 4000x faster compared to the higher- fidelity thermal hydraulics model, enabling real-time operational intelligence and applications in long horizon predictions. While the surrogate model framework is demonstrated for the particular case of FHR, the hybrid physical/data-driven modeling approach including the network of surrogates and the underlying modularity has the potential to be applied to other physical asset systems.

Digital Twins↗

DOE OSTI2024

Sensitivity Analyses of Natural Convection in the AP1000 Passive Containment Cooling System Following LBLOCA Using CFD

The AP1000 power plant has implemented multiple layers of passive safety systems to enhance reactor safety and ensure system integrity during postulated scenarios, such as loss of coolant accident (LOCA) and main steam line break (MSLB). Among these, the passive containment cooling system (PCCS) is a safety-related system designed to prevent the containment from exceeding the design limits of both pressure and temperature following a postulated design basis accident, by transferring heat from the steel containment vessel to the atmosphere. During LOCA, natural air convection plays a vital role in removing heat from the steel containment vessel to the atmosphere. This is integrated with the condensate film from the Passive Containment Cooling Water Storage Tank (PCCWST). Air natural convection is also the only heat removal mechanism after the dry-out of the PCCWST and during the loss of shutdown decay heat removal (LOSDHR) event. Therefore, it is essential to investigate the thermal hydraulics behavior of the containment under transient conditions to ensure its safety and integrity. This paper presents a simplified CFD/ANSYS FLUENT model developed to analyze the impact of different parameters, such as air relative humidity, air temperature, and steel containment temperature, on the natural convection capability to remove the decay heat following LOCA. The model aims to examine the thermal behavior of the containment and provides a comprehensive understanding of the system's natural convection capability during postulated accidents. The results obtained from the model can be used to improve the safety and integrity of the containment system and provide valuable insights for future research in the field.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2020

Development and Validation of a Simulation Testbed for the Intelligent Building Agents Laboratory (IBAL) Using TRNSYS

This paper documents the development and validation of a dynamic primary cooling and thermal storage system simulation testbed. The system simulation testbed, sIBAL, is based on the Intelligent Building Agents Laboratory (IBAL) at the National Institute of Standards and Technology (NIST), which is a research infrastructure and testbed for the development, evaluation, and demonstration of intelligent control algorithms. The sIBAL testbed developed in this project will serve as a virtual twin of the real facility for future control algorithm development. The details of the methodologies used to develop and validate the simulation testbed, which replicates the dynamic behaviors of the primary cooling and ice storage system in the IBAL facility, are presented. The simulation testbed was developed in TRNSYS using built-in component models and MATLAB functions to replicate the two water-cooled chillers, a thermal storage tank, pumps, valves and other components for four different operation modes. Experiments on IBAL components were designed and executed to generate experimental data for model development and verification of the simulation platform. The validation of the simulation results was carried out in two phases: 1) independent component simulations for the chillers and thermal storage tank, and 2) a combined testbed simulation of the entire hydronic system. Comparison of simulation results to the experimental data obtained from the IBAL facility showed errors within 1 °C for the temperatures outputs of both the chiller and the thermal storage model. The error is within an acceptable range for further intelligent control algorithms development. The findings from the study are summarized and presented along with areas where additional research is needed. In addition, data filtering procedures and model refinement measures utilized to improve the accuracy and accelerate the computation time of the simulation are presented.

cooling systems↗

DOE OSTI2025

Anomalous thermal transport in Eshelby twisted van der Waals nanowires

Dislocations in van der Waals (vdW) layered nanomaterials induce strain and structural changes that substantially impact thermal transport. Understanding these effects could enable the manipulation of dislocations for improved thermoelectric and optoelectronic applications, but experimental insights remain limited. Here, in this study, we use synthetic Eshelby twisted vdW GeS nanowires (NWs) with single screw dislocations as a model system to explore the interplay between dislocation-induced structural modifications and lattice thermal conductivity. Our measurements reveal a monoclinic structure stabilized by the dislocation, leading to a substantial drop in thermal conductivity for larger-diameter NWs (70% at room temperature), supported by first-principles calculations. Interestingly, we also find an anomalous enhancement of thermal conductivity with decreasing diameter in twisted NWs, contrary to typical trends in non-twisted GeS NWs. This is attributed to increased conductivity near the NW cores due to compressive strain around the central dislocations, and aligns with a density-functional-theory-informed core–shell model. Our results highlight the critical role of dislocations in thermal conduction, providing fundamental insights for defect and strain engineering in advanced thermal applications.

36 MATERIALS SCIENCE↗

DOE OSTI2021

A measure of active interfaces in supported catalysts for high-temperature reactions

Formulating knowledge of structure-function relationships in heterogeneous catalysis is central to the design of efficient catalysts; yet, the elucidation of dominant reaction sites has remained as a challenge. Here, in this paper, we present a methodology that can be used to visualize metal-gas and metal-oxide-gas interfaces in three dimensions and to quantify their catalytic activity levels. As a case study, CH 4 oxidation occurring in a Pt/CeO 2 system is chosen. By employing thermally robust Pt@CeO 2 model catalysts with size-tunable and monodisperse cores, and gas-permeable shells, we reconstruct a series of structures in 3D via electron tomography and match the information to activity data and theoretical calculations. This strategy reveals that two different interfaces catalyze the CH 4 oxidation and that their contribution to the overall rate changes with the Pt size, temperature, and gas atmosphere. Our results provide an analytic platform on which to explore reaction pathways and mechanisms applicable to multiple reactions and materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Comparison of carrier localization effects between InAs quantum dashes and quantum dots in a DWELL (dashes- or dots-in-a-well) configuration

The optical properties of InAs quantum dashes (QDashes) grown on InP and InAs quantum dots (QDots) grown on GaAs in a dashes- or dots-in-a-well (DWELL) configuration are comparatively investigated using temperature-dependent photoluminescence (PL) measurements. The trends in PL characteristics such as exciton energy, spectral bandwidth and integrated intensity with respect to temperature are found to be distinctly dissimilar between the two systems. A rate-equation model involving exciton recombination and thermal transfer in a localized-state ensemble is used to quantitively interpret the experimental data. The results of this study suggest that QDashes in this configuration exhibit PL properties more consistent with a lower degree of carrier localization compared to QDots. A preliminary structural analysis highlighting the shape/size differences between the two nanostructures is also presented.

carrier localization↗

DOE OSTI2024

Critical slowing of the spin and charge density wave order in thin film Cr following photoexcitation

We report on the evolution of the charge density wave (CDW) and spin density wave (SDW) orders of a chromium film following photoexcitation with an ultrafast optical laser pulse. The CDW is measured by ultrafast time-resolved x-ray diffraction of the CDW satellite that tracks the suppression and recovery of the CDW following photoexcitation. We find that as the temperature of the film approaches a discontinuous phase transition in the CDW and SDW orders, the time scales of recovery increase exponentially from the expected thermal time scales. We extend a Landau model for SDW systems to account for this critical slowing with the appropriate boundary conditions imposed by the geometry of the thin film system. This model allows us to assess the energy barrier between the available CDW/SDW states with different spatial periodicities.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Thermal Performance of Spandrel Assemblies in Glazed Wall Systems: Laboratory Test Design – Challenges and Test Results

Accurate thermal performance calculation procedures for opaque spandrel areas in curtain wall and window wall systems are essential for rating systems when comparing spandrel systems. However, there is a lack of consensus in thermal modeling needed for accurately characterizing heat transfer through spandrel assemblies due to the complex arrangement of materials and structural components. Several studies indicate that conventional 2D thermal simulations may overestimate R-values by 30% compared to physical testing and 3D simulations. Detailed simulations and well-curated laboratory test data are necessary to build confidence in simulation models, which will later be used to develop correlations to improve widely used conventional 2D thermal simulations. This study aims to experimentally test heat transfer through various spandrel assemblies to validate 3D simulation models. Also, the challenges of conducting a thorough testing design along with the solutions would be documented. The team developed a design for testing spandrel assemblies, making appropriate modifications to the existing heat, air, and moisture (HAM) chamber to accommodate the testing needs. Two moveable baffles were designed and fabricated to guide airflow direction parallel to the test article surface. The data acquisition capabilities in the chamber were upgraded to add more than two hundred sensors to the climate and indoor side of the chamber. The goal is to provide a quality dataset for validating complex 3D modeling simulations, which will be used to develop improved thermal simulation techniques that more accurately represent the thermal behavior of spandrel assemblies and their integration within the building envelope. This paper will summarize the results for the boundary conditions of the testing and the temperature variation across different locations of the spandrel assemblies.

Kunwar, Niraj [ORNL] (ORCID:0000000263457652)↗

DOE OSTI2024

An analysis of parameters affecting ampacity in aircraft bipolar MVDC power cables via coupled electrical, thermal, and computational fluid dynamic modelling

Abstract The next generation of aircraft, including more electric aircraft and all‐electric aircraft (AEA), requires electric power systems with high power density and low system mass specifications. Increasing the voltage of the system to the range of a few kV, medium voltage (MV), is a reasonable approach to achieving high‐power‐density and low‐system‐mass EPSs for aircraft applications. Higher voltages, however, pose many challenges for aviation MV power cables such as arcs and arc tracking, partial discharges (PDs), and thermal management. In this regard, thermal management is more challenging since heat transfer by convection is greatly reduced at wide‐body aircraft's cruising altitudes due to the reduced air pressure. In this paper, a finite element method (FEM) model is developed in COMSOL Multiphysics for an aircraft bipolar MVDC (±5 kV) power cable. Using the model, the maximum permissible cable current at a low pressure of 18.8 kPa (at an altitude of 12.2 km from sea level, the usual cruising altitude for wide‐body aircraft) is calculated. Also, an analytical model is developed based on analytical and proven empirical correlations governing conductive, radiative, and convective heat transfers at the steady state to estimate the ampacity of the bipolar cable system at reduced pressure. It was shown that the proposed analytical model can be used for atmospheric pressure and systems with a larger number of poles, expanding its range of applications. The results of the FEM and analytical models correlate at wide ranges of parameters such as ambient temperature, duct size, distance between the positive and negative pole cables, and the overall diameter of the cables. The influence of horizontal and vertical arrangement of poles is included in the analytical model. The results of this study can be used to design bipolar MVDC power cable systems for the envisaged wide‐body AEA.

Azizi, Arian↗

DOE OSTI2025

Gas-Cooled High-Temperature Pebble-Bed Reactor Reference Plant Model Updates

This work presents the latest improvements to, and investigations performed with, the pebble-bed high-temperature gas-cooled reactor (PB-HTGR) reference plant models for the United States Nuclear Regulatory Commission. These models serve as the foundation for the future development of detailed design evaluation models based on license applications. The reference plant models have been developed with the Comprehensive Reactor Analysis Bundle, or BlueCRAB, which is the code suite proposed for non-light-water reactor systems safety analysis. It incorporates various simulation tools developed by the Nuclear Energy Advanced Modeling and Simulation program, including the Griffin code for reactor physics, the Pronghorn and SAM codes for core thermal fluids, the BISON code for solid conduction and fuel performance, and the SAM code for system analysis. The primary objective of the work that was performed was to assess BlueCRAB’s level of readiness for modeling a PB-HTGR. To do so, we first developed numerical models in BlueCRAB that include the key physics for this technology to ensure an adequate level of fidelity for modeling PB-HTGR core performance and for performing multiphysics simulations for equilibrium core conditions and different accident scenarios. Then we simulated transient scenarios, including depressurized and pressurized loss of forced cooling accidents, over-cooling, and control rod withdrawal events with delayed and prompt supercritical reactivity insertions. The analysis in this report includes comparisons of the 2D thermal fluid porous media models in Pronghorn and SAM, and comparisons of coupled SAM/Griffin/SAM and coupled Pronghorn/Griffin for depressurized and pressurized loss of forced cooling, over-cooling, and control rod withdrawal events. In addition, we compare 3D, 2D, and 0D/PKE neutronic models for the two control rod withdrawal scenarios with coupled Pronghorn/Griffin. The comparisons show that the BlueCRAB models lead to physically intuitive solutions for the scenarios examined. The changes in the various scalar and vector fields, such as neutron flux, power, temperature, density, pressure, and velocity, are within the expected ranges, and their distributions can be explained by the system response of the transients and the geometric and material variations. Several comparisons suggest that the porous media models in Pronghorn and SAM can lead to similar solutions, even though they are based on different methodologies. This work further highlights the need for flexible tools with various levels of fidelity to cover the breadth and depth of needs that may arise in future technical evaluations of the PB-HTGR. We believe that the BlueCRAB capabilities will be a significant asset for confirmatory analyses in order to resolve important safety questions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2023

Modeling phonons in nanomaterials

The vibrational properties of nanomaterials play fundamentally important roles in determining their stability and various functionalities, from rippling of otherwise flat two-dimensional (2D) monolayers to the extreme thermal transport behaviors observed in one-dimensional (1D) and 2D carbon-based systems. Here we discuss modeling efforts to understand how phonons in lower dimensional materials differ from those in bulk, and how these behaviors manifest in their unusual properties. In particular, this chapter gives a heavy focus to calculations of phonon dynamics, interactions, and lattice thermal transport in 1D chains and nanotubes, 2D monolayers, and van der Waals coupled bulk materials from first principles methods. This chapter will discuss important phonon calculation details related to nanomaterials with examples and discussions providing both historical and modern perspectives.

Lindsay, Lucas↗

DOE OSTI2020

Microcanonical rates from ring-polymer molecular dynamics: Direct-shooting, stationary-phase, and maximum-entropy approaches

We address the calculation of microcanonical reaction rates for processes involving significant nuclear quantum effects using ring-polymer molecular dynamics (RPMD), both with and without electronically non-adiabatic transitions. After illustrating the shortcomings of the naive free-particle direct-shooting method, in which the temperature of the internal ring-polymer modes is set to the translational energy scale, we investigate alternative strategies based on the expression for the microcanonical rate in terms of the inverse Laplace transform of the thermal reaction rate. It is shown that simple application of the stationary-phase approximation (SPA) dramatically improves the performance of the microcanonical rates using RPMD, particularly in the low-energy region where tunneling dominates. Using the SPA as a Bayesian prior, numerically exact RPMD microcanonical rates are then obtained using maximum entropy inversion of the thermal reaction rates for both electronically adiabatic and non-adiabatic model systems. Finally, the direct-shooting method is revisited using the SPA-determined temperature for the internal ring-polymer modes, leading to a simple, direct-simulation method with improved accuracy in the tunneling regime. This work suggests a general strategy for the extraction of microcanonical dynamical quantities from RPMD (or other approximate thermal) simulations.

Tao, Xuecheng↗

DOE OSTI2023

Thermal acclimation of plant photosynthesis and autotrophic respiration in a northern peatland

Peatlands contain one-third of global soil carbon (C), but the responses of peatland ecosystems to long-term warming are not well understood. Here, we pursue an emergent understanding of warming effects on ecosystem C fluxes at peatlands by constraining a process-oriented model, the terrestrial ECOsystem model, with observational data from a long-term warming experiment at the Spruce and Peatland Responses Under Changing Environments site. Model-based assessments show that ecosystem-level photosynthesis and autotrophic respiration exhibited significant thermal acclimation, with temperature sensitivities being linearly decreased with warming. Using the thermal-acclimated parameter values, simulated gross primary production, net primary production, and plant autotrophic respiration (R a ), were all lower than those simulated with non-thermal acclimated parameter values. In contrast, ecosystem respiration simulated with thermal acclimated parameter values was higher than that simulated with non-thermal acclimated parameter values. Net ecosystem CO 2 exchange was much higher after constraining model parameters with observational data from the warming treatments, releasing C at a rate of 28.3 g C m -2 yr -1 °C -1 . Our data-model integration study suggests that peatlands are likely to release more C than previously estimated. Earth system models may overestimate C uptake by peatlands under warming if physiological thermal acclimation of plants is not incorporated. Thus, it is critical to consider the long-term physiological thermal acclimation of plants in the models to better predict global C dynamics under future climate and their feedback to climate change.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Monsoonal MCS Initiation, Rainfall, and Diurnal Gravity Waves over the Bay of Bengal: Observation and a Linear Model

Abstract Previous observational studies have indicated that mesoscale convective systems (MCSs) contribute the majority of precipitation over the Bay of Bengal (BoB) during the summer monsoon season, yet their initiation and propagation remain incompletely understood. To fill this knowledge gap, we conducted a comprehensive study using a combination of 20-yr satellite observations, MCS tracking, reanalysis data, and a theoretical linear model. Satellite observations reveal clear diurnal propagation signals of MCS initiation frequency and rainfall from the west coast of the BoB toward the central BoB, with the MCS rainfall propagating slightly slower than the MCS initiation frequency. Global reanalysis data indicate a strong association between the offshore-propagating MCS initiation frequency/rainfall and diurnal low-level wind perturbations, implying the potential role of gravity waves. To verify the hypothesis, we developed a 2D linear model that can be driven by realistic meteorological fields from reanalysis. The linear model realistically reproduces the characteristics of offshore-propagating diurnal wind perturbations. The wind perturbations, as well as the offshore propagation signals of MCS initiation frequency and rainfall, are associated with diurnal gravity waves emitted from the coastal regions, which in turn are caused by the diurnal land–sea thermal contrast. The ambient wind speed and vertical wind shear play crucial roles in modulating the timing, propagation, and amplitude of diurnal gravity waves. Using the linear model and satellite observations, we further show that the stronger monsoonal flows lead to faster offshore propagation of diurnal gravity waves, which subsequently control the offshore propagation signals of MCS initiation and rainfall. Significance Statement Rainfall over the Bay of Bengal (BoB) is primarily contributed by large and organized rainfall systems in the summer monsoon season. During this season, these systems are commonly observed over the east coast of India around midnight, the western BoB in the morning, and the central BoB in the afternoon. This eastward rainfall propagation is confirmed by observations, reanalysis data, and a theoretical model to have a strong association with atmospheric diurnal gravity waves. These waves are caused by land–sea thermal contrast and can trigger rainfall systems over the offshore regions. We also found that the diurnal gravity waves, as well as their triggered rainfall systems, can be greatly modulated by the large-scale monsoonal flows. All the above findings improve our understanding of diurnal rainfall cycle over the BoB and may contribute to the future improvement of rainfall forecast over the region.

Meteorology & Atmospheric Sciences↗

DOE OSTI2025

Application of Ocean Thermal Energy Conversion (OTEC) Systems for Powering Safety Monitoring Systems of Offshore Oil and Gas Operations. OESI 2.0 M-1 T-1-P1.1 – Objective 1 Interim Report: Summary of Offshore Oil and Gas Well Monitoring Systems and Process Power Requirements

This report is part of a comprehensive research initiative aimed at evaluating the technoeconomic feasibility of deploying greenhouse gas (GHG) emissions monitoring systems powered by small-scale renewable energy sources on the U.S. Outer Continental Shelf (OCS). The study has three primary objectives: summarizing available GHG monitoring systems for the OCS, modeling a renewable marine energy source (specifically ocean thermal energy conversion, or OTEC), and modeling the collocation of OTEC power sources with offshore oil and gas activities to offset power demand.

02 PETROLEUM↗

DOE OSTI2021

Ambiguities in the hadro-chemical freeze-out of Au+Au collisions at SIS18 energies and how to resolve them

The thermal fit to preliminary HADES data of Au+Au collisions at $\sqrt{s_{NN}}$=2.4 GeV shows two degenerate solutions at T≈50 MeV and T≈70 MeV. The analysis of the same particle yields in a transport simulation of the UrQMD model yields the same features, i.e. two distinct temperatures for the chemical freeze-out. While both solutions yield the same number of hadrons after resonance decays, the feeddown contribution is very different for both cases. This highlights that two systems with different chemical composition can yield the same multiplicities after resonance decays. The nature of these two minima is further investigated by studying the time-dependent particle yields and extracted thermodynamic properties of the UrQMD model. It is confirmed, that the evolution of the high temperature solution resembles cooling and expansion of a hot and dense fireball. The low temperature solution displays an unphysical evolution: heating and compression of matter with a decrease of entropy. These results imply that the thermal model analysis of systems produced in low energy nuclear collisions is ambiguous but can be interpreted by taking also the time evolution and resonance contributions into account.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2026

Project Phase 1 Report: Reducing Data Center Peak Cooling Demand and Energy Costs With Cold Underground Thermal Energy Storage (Cold UTES)

Cold Underground Thermal Energy Storage (Cold UTES) is an ultra-long duration grid energy storage technology. With Cold-UTES, low-cost grid power is converted to cold thermal energy and stored in the native subsurface rock at the point of use. Cold UTES is one approach within the general category of engineered geothermal systems. Cold UTES for peak-hour cooling of data centers (DCs) was studied for deployment in Maricopa County, Arizona and Loudoun County, Virgina using thermal storage capacities from 4 GWh-th to over 1,000 GWh-th (>1 Terawatt-hour). The two sites have different power and transmission systems, available grid energy resources, daily and seasonal load profiles, weather conditions, and grid regulatory requirements. The study results indicate high value for both locations and because of this, likely indicates value across most of the US and the world. The basis of the study was a 1,000 MW-e hourly electric use of DC computing and auxiliary loads, which was modeled as 1 GW-th of thermal load to a dry cooled heat rejection system - i.e. a cooling system that does not consume water. The electric power required for cooling the DC varies as the air temperature changes. In cold weather only the dry-coolers are used, with an electrical load for cooling load as low as 10 MW-e. In hot summer hours, chillers and dry- coolers are required, which raises the electrical load for cooling load to as much as 300 MW-e. The continuous and peak cooling electrical loads result in a grid interconnection requirement of no less than 1,300 MW-e. From both a grid and thermal design modeling perspective the 1.3 GW-e could either be a single facility or result from the total load at multiple sites.

15 GEOTHERMAL ENERGY↗

DOE OSTI2022

A combined theoretical-experimental investigation of thermal transport in low-dose irradiated thorium dioxide

During reactor operation, nuclear fuels are subject to extreme temperature and irradiation conditions which can significantly degrade the fuel's thermal transport properties. The reduction in thermal conductivity of the fuel as a result of irradiation-induced lattice defects is arguably the most important fuel performance metric in regard to reactor efficiency and safety. Because thorium dioxide (ThO2) is suitable as a model system for more complex materials such as UO2 and its mixed oxides, we present a theoretical investigation of thermal conductivity reduction seen in defect-bearing thorium dioxide and compare directly to experimental measurements. Phonon-mediated thermal transport of the fuel is modeled by a solution to the Boltzmann transport equation (BTE) for phonons. A cluster dynamics (CD) model for lattice defect evolution during irradiation predicts defect densities which are used as input to the BTE for modeling phonon-defect scatterings. Phonon scatterings by lattice defects include those from point defects and vacancy clusters and interstitial clusters of various sizes. The CD model is benchmarked against structural defect characterization of irradiated thorium dioxide using electron microscopy. Thermal conductivity predicted by the BTE model is compared to measured values for irradiated thorium dioxide specimens below room temperature to isolate effects of phonon-defect scattering from intrinsic 3-phonon processes, which dominate at higher temperatures. The computed conductivity values are in partial agreement at temperatures close to room temperature while slight deviations are observed at the lowest measured temperatures, suggesting that implemented phonon-defect scattering cross-section expressions may not be adequate for low temperatures. The presented work provides a necessary investigation of the influence of irradiation induced defects on fuel performance and represents a first step toward a full characterization of phonon mediated thermal transport in irradiated materials with complex defect microstructure.

36 MATERIALS SCIENCE↗