SEARCH · Search NASA

Results for “TRANSPORT EQUATION”

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 55 records · Page 3

All sources NASA NTRS DOE OSTI

DOE OSTI2025

A geometric perspective on kinetic matter–radiation interaction and moment systems

Here, in this paper, we provide a geometric perspective on the kinetic interaction of matter and radiation, based on a pair bracket approach. We discuss the interaction of kinetic theories via dissipative brackets, with our fundamental example being the coupling of matter, described by the Boltzmann equation, and radiation, described by the radiation transport equation. We explore the transition from kinetic systems to their corresponding moment systems, provide a Hamiltonian description of such moment systems, and give a geometric interpretation of the moment closure problem for kinetic theories. As an application, we discuss in detail diffusion radiation hydrodynamics as an example of a pair bracket formulation on a space of moments corresponding to kinetic matter–radiation interaction. Additionally, using the variable moment closure framework of Burby [J. W. Burby, Variable-moment fluid closures with Hamiltonian structure, Sci. Rep. 13 (2023) 18286, doi:10.1038/s41598-023-45416-5], we show how to construct Hamiltonian moment closures for kinetic transport equations with arbitrary Hamiltonian. Using this general construction, we derive novel Hamiltonian moment closures for pure radiation transport.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

On the evolution of scalar iso-surface area density in a turbulent mixing layer

In this paper, the transport equation for scalar iso-surface area density (Σ) in a turbulent, temporally developing mixing layer is examined. Exploring the spatial and temporal evolution of the terms in the Σ transport equation is vital to improving our understanding of turbulent flows characterized by distinct interfaces, e.g. the flame surface or the turbulent/non-turbulent interface. Previous work reported by the authors identified that Σ exhibits self-similar behaviour consistent with the development of the temporal mixing layer. Accordingly, each of the terms in the Σ transport equation is found to behave in a self-similar manner, though there are notable differences in the self-similar behaviours for each term. Based on the results presented herein, it is suggested that the rate of change of Σ (∂Σ/∂t) and the advection term scale with hλ Φ /ΔU, where h is the width of the mixing layer, λ Φ is the scalar Taylor length scale and ΔU is the velocity difference. The production and destruction terms are found to scale with an additional factor (ReSc) 1/2 . In contrast, the molecular diffusion term is found to scale with a factor (ReSc) –1/2 compared to ∂Σ/∂t. Importantly, it is found that the difference between the production and destruction terms, or net surface ‘stretch’, scales with the same factor as ∂Σ/∂t and the advection term, which may have a significant impact on how the evolution of Σ is understood and modelled in turbulent flows.

42 ENGINEERING↗

DOE OSTI2024

A Second Moment Method for k -Eigenvalue Acceleration with Continuous Diffusion and Discontinuous Transport Discretizations

The second moment method is a linear acceleration technique that couples the transport equation to a diffusion equation with transport-dependent additive closures. The resulting low-order diffusion equation can be discretized independent of the transport discretization, unlike diffusion synthetic acceleration, and is symmetric positive definite, unlike quasidiffusion. While this method has been shown to be comparable to quasidiffusion in iterative performance for fixed source and time-dependent problems, it is largely unexplored as an eigenvalue problem acceleration scheme due to the belief that the resulting inhomogeneous source makes the problem ill posed. Recently, a preliminary feasibility study was performed on the second moment method for eigenvalue problems. The results suggested comparable performance to quasidiffusion and more robust performance than diffusion synthetic acceleration. This work extends the initial study to more realistic reactor problems using state-of-the-art discretization techniques. Finally, the results in this paper show that the second moment method is more computationally efficient than its alternatives on complex reactor problems with unstructured meshes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2026

Offline Maximizing Minimally Invasive Proper Orthogonal Decomposition for Reduced-Order Modeling of S n Radiation Transport

Deterministic solutions to the Sn radiation transport equation can be computationally expensive to calculate. Reduced-order modeling enables efficient approximation of the full-order model (FOM) solution. We propose a novel method for constructing reduced-order models (ROMs) of the S n radiation transport equation, offline maximizing minimally invasive (OMMI) proper orthogonal decomposition (POD). POD uses the method of snapshots to create a reduced-order basis for constructing an ROM. Minimally invasive POD leverages the sweep infrastructure existing in deterministic transport codes to create a POD-based ROM, even when infeasible by traditional methods. Offline maximizing minimally invasive proper orthogonal decomposition (OMMI-POD) extends minimally invasive POD by performing sweeps offline, therefore maximizing the potential speedup. OMMI-POD does so by creating a library of reduced systems from a training set. This library of reduced systems is then interpolated to provide a rapid approximate solution of the S n radiation transport equation. The model is evaluated on a set of test problems, achieving a low error with a 466 times speedup over the FOM. Also presented is a study of the effect of sampling method on the performance of OMMI-POD, specifically comparing naive uniform sampling to the more accurate and computationally expensive greedy sampling.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Assessment of diffuse-interface methods for compressible multiphase fluid flows and elastic-plastic deformation in solids

This work describes three diffuse-interface methods for the simulation of immiscible, compressible multiphase fluid flows and elastic-plastic deformation in solids. The first method is the localized-artificial-diffusivity approach of Cook, Subramaniam et al., and Adler and Lele, in which artificial diffusion terms are added to the individual phase mass fraction transport equations and are coupled with the other conservation equations. The second method is the gradient-form approach that is based on the quasi-conservative method of Shukla et al., in which the diffusion and sharpening terms (together called regularization terms) are added to the individual phase volume fraction transport equations and are coupled with the other conservation equations. The third approach is the divergence-form approach that is based on the fully conservative method of Jain et al., in which the regularization terms are added to the individual phase volume fraction transport equations and are coupled with the other conservation equations. In the present study, all three diffuse-interface methods are used in conjunction with a four-equation, multicomponent mixture model, in which pressure and temperature equilibria are assumed among the various phases. The primary objective of this work is to compare these three methods in terms of their ability to: maintain constant interface thickness throughout the simulation; conserve mass, momentum, and energy; and maintain accurate interface shape for long-time integration. The second objective of this work is to consistently extend these methods to model interfaces between solid materials with strength. To assess and compare the methods, they are used to simulate a wide variety of problems, including (1) advection of an air bubble in water, (2) shock interaction with a helium bubble in air, (3) shock interaction and the collapse of an air bubble in water, and (4) Richtmyer–Meshkov instability of a copper–aluminum interface. The current work focuses on comparing these methods in the limit of relatively coarse grid resolution, which illustrates the true performance of these methods. In conclusion, this is because it is rarely practical to use hundreds of grid points to resolve a single bubble or drop in large-scale simulations of engineering interest.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Implicit and coupled fluid plasma solver with adaptive Cartesian mesh and its applications to non-equilibrium gas discharges

In this work, we present a new fluid plasma solver with adaptive Cartesian mesh (ACM) based on a full-Newton (nonlinear, implicit) scheme for non-equilibrium gas discharge plasma. The electrons and ions are described using drift-diffusion approximation coupled to Poisson equation for the electric field. The electron-energy transport equation is solved to account for electron thermal conductivity, Joule heating, and energy loss of electrons in collisions with neutral species. The rate of electron-induced ionization is a function of electron temperature and could also depend on electron density (important for plasma stratification). The ion and gas temperature are kept constant. The transport equations are discretized using a non-isothermal Scharfetter-Gummel scheme to resolve possible large temperature gradients in the sheaths. We demonstrate the new solver for simulations of direct current (DC) and radiofrequency (RF) discharges. The implicit treatment of the coupled equations allows using large time steps. The full-Newton method (FNM) enables fast nonlinear convergence at each time step, offering significantly improved simulation efficiency. We discuss the selection of time steps for solving different plasma problems. The new solver enables solving several problems we could not solve before with existing software: two- and three-dimensional structures of the entire DC discharges including cathode and anode regions, electric field reversals and double-layer formation, the normal cathode spot and an anode ring, moving striations in diffuse and constricted DC discharges, and standing striations in RF discharges. The developed FNM-ACM technique offers many benefits for tackling the disparity of gas discharge plasma systems' time scales and nonlinearity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Correlation Function Approach for Diffusion in Confined Geometries

This paper describes a formalism for extracting spatially varying transport coefficients from simulations of a molecular fluid in a nano channel. This approach is applied to self-diffusion of a Lennard-Jones fluid confined between two parallel surfaces. A numerical grid is laid over the domain confining the fluid, and fluid properties are projected onto the grid cells. The time correlation functions between properties in different grid cells are calculated and can be used as the basis for a fitting procedure for extracting spatially varying diffusion coefficients from the simulation. Results for the Lennard-Jones system show that transport behavior varies sharply near the liquid-solid boundary and that the changes depend on the details of the liquid-solid interaction. A quantitative difference between the reduced and detailed models is discussed. It is found that the difference could be associated with assumptions about the form of the transport equations at molecular scales in lieu of problems with the method itself. The study suggests that this approach to fitting molecular simulations to continuum equations may guide the development of appropriate coarse-grained equations to model transport phenomena at nanometer scales.

Nanoscale flow, Transport, molecular simulation↗

DOE OSTI2020

Benchmark of interfacial area concentration approaches for the two-fluid model in gas-dispersed condensing flow

Accurate predictions of interfacial area concentration (IAC) are required for closure of the condensation phase change term. The interfacial area transport equation (IATE) can dynamically predict interfacial area concentration and is increasingly being used in the benchmarking of phase change flows over static based correlations which require a flow regime map. For higher flow regimes, the two-group two-fluid IATE model can separately model the transport of spherical bubbles and cap/slug bubbles. Recently, a modified formulation of the one-group and two-group two-fluid model with interfacial area transport equation was proposed for gas-dispersed condensing flows. Three datasets have been used for validation consisting of condensing two-phase flow data in the unheated section of a vertical annulus with conditions spanning a range of pressures, inlet subcoolings and mass fluxes for benchmarking the various IAC approaches. For bubbly flows, the coupled void transport - IATE simulation results in good agreement in the predictions of the void fraction and the interfacial area concentration under nearly all sixty conditions. Considering the proposed two-group two-fluid IATE model formulation, the group-2 Nusselt number correlation is investigated to account for the appropriate heat transfer length scale for group-2 bubbles. Validation of the proposed two-group two-fluid model against the new dataset shows good agreement with experimental data with significant improvement over the one-group IATE model for conditions with substantial group-2 void fraction throughout the simulation domain.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2022

A Forward Analytic Model of Neutron Time-of-Flight Signals for Inferring Ion Temperatures from MagLIF Experiments

A forward analytic model is required to rapidly simulate the neutron time-of-flight (nToF) signals that result from magnetized liner inertial fusion (MagLIF) experiments at Sandia’s Z Pulsed Power Facility. Various experimental parameters, such as the burn-weighted fuel-ion temperature and liner areal density, determine the shape of the nToF signal and are important for characterizing any given MagLIF experiment. Extracting these parameters from measured nToF signals requires an appropriate analytic model that includes the primary deuterium-deuterium neutron peak, once-scattered neutrons in the beryllium liner of the MagLIF target, and direct beamline attenuation. Here, mathematical expressions for this model were derived from the general-geometry time- and energy-dependent neutron transport equation with anisotropic scattering. Assumptions consistent with the time-of-flight technique were used to simplify this linear Boltzmann transport equation into a more tractable form. Models of the uncollided and once-collided neutron scalar fluxes were developed for one of the five nToF detector locations at the Z-Machine. Numerical results from these models were produced for a representative MagLIF problem and found to be in good agreement with similar neutron transport simulations. Twenty experimental MagLIF data sets were analyzed using the forward models, which were determined to only be significantly sensitive to the ion temperature. The results of this work were also found to agree with values obtained separately using a zero scatter analytic model and a high-fidelity Monte Carlo simulation. Finally, inherent difficulties in this and similar techniques are identified, and a new approach forward is suggested.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Monte-Carlo Neutrino Transport in Neutron Star Merger Simulations

Gravitational waves and electromagnetic signals from merging neutron star binaries provide valuable information about the the properties of dense matter, the formation of heavy elements, and high-energy astrophysics. To fully leverage observations of these systems, we need numerical simulations that provide reliable predictions for the properties of the matter unbound in these mergers. An important limitation of current simulations is the use of approximate methods for neutrino transport that do not converge to a solution of the transport equations as numerical resolution increases, and thus have errors that are impossible to quantify. Here, we report on a first simulation of a binary neutron star merger that directly solves the transport equations using Monte-Carlo techniques. We demonstrate that a computationally cheap Monte-Carlo transport can be implemented without ill effects from sampling noise. We apply Monte-Carlo transport to an unequal mass neutron star binary merger up to 5ms past merger, and report on the properties of the matter and neutrino outflows. Finally, we compare our results to the output of our best approximate `M1' transport scheme, demonstrating that an M1 scheme that carefully approximates the neutrino energy spectrum only leads to ~10% uncertainty in the composition and velocity of the ejecta, and ~20% uncertainty in the electron (anti)neutrino luminosities and energies. Here, the most significant disagreement found between M1 and Monte-Carlo results is a factor of ~2 difference in the luminosity of heavy-lepton neutrinos.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Simultaneous velocity and density measurements of fully developed Rayleigh-Taylor mixing

Here, the dynamics of molecular mixing and the energy transfer process in the Rayleigh-Taylor instability (RTI) are studied through the collection of simultaneous density-velocity measurements. These experiments provide simultaneous density-velocity field measurements, in contrast to previous point measurements. Statistically stationary experiments are performed in a “convective-type” gas tunnel facility, with density contrast achieved through the injection of helium into the bottom stream. Three experiments at Atwood number ≈ 0.1 are captured at three outer scale Reynolds numbers Re = 520, 2260, and 4050. Particle image velocimetry and laser induced fluorescence are employed simultaneously. Statistics of the density and velocity show self-similar collapse of RTI profiles at large Reynolds number Re >2000. Flat velocity profiles indicate homogeneous turbulence characteristics in the core of the mixing region. Significant anisotropy develops in the flow, with horizontal velocity fluctuations being only 60% of the vertical velocity fluctuations. The turbulent mass flux is found to be asymmetric about the centerline, with increased peak towards the spike. Measurements of the molecular mixing show that mixing is maximized at the core of the flow and increases with increased Reynolds number. The transport equation of density-specific-volume correlation b shows that it is mostly produced in the core of the mixing region, and that the spatial evolution of its profile is the result of transport by bulk motion of the bubble and spike. Energy transfer from gravitational potential energy to turbulent kinetic energy and viscous dissipation is observed to occur in the experiment with a ratio of dissipated energy to potential energy released of 38%. The analysis of the turbulent kinetic energy transport equation budget reveals that production is the primary mechanism towards the growth of turbulent kinetic energy in the core of the flow, and is asymmetrically slightly skewed towards the spike. However, it is through the transport that the strong advection at the edges of the mixing region is maintained.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

Griffin: A MOOSE-based reactor physics application for multiphysics simulation of advanced nuclear reactors

Griffin is a Multiphysics Object-Oriented Simulation Environment (MOOSE) based reactor physics application for multiphysics simulations of advanced reactor designs jointly developed by Idaho National Laboratory and Argonne National Laboratory. This paper summarizes the motivation, significance, architecture, design, and features of Griffin. Griffin offers flexible and extensible features to address the challenges associated with advanced reactor designs. These features range from fundamental particle transport to specific reactor physics tasks. The features cover a wide range including on-the-fly and traditional two-step cross-section generation methods, steady-state and transient transport solvers suitable for both heterogeneous and homogeneous models, high-fidelity depletion where thousands of isotopes can be tracked and low-fidelity depletion characterized by burnup, etc. The most fundamental aspect that sets Griffin apart from other reactor analysis codes is that it is developed based on the MOOSE framework. A modular development approach is strongly enforced, with multiphysics being an essential element considered since the beginning of Griffin’s development. Griffin links various MOOSE physics modules and couples to other MOOSE-based applications and non-MOOSE-based applications for multiphyiscs simulations. Griffin includes three modules: ISOXML for preparing and managing multigroup cross sections, radiation transport for solving the neutron transport equation, and reactor analysis for user-oriented reactor physics analysis functionalities. Griffin uses various finite element methods for spatial discretization, multigroup approximation for energy discretization and discrete ordinates method, spherical harmonics expansion method, and diffusion approximation for streaming direction discretization to solve the neutron transport equation. Griffin’s flexibility is evidenced through Griffin’s various applications to fast reactor, high-temperature reactor, pebble bed reactor, molten salt reactor, and microreactor designs. Griffin development follows the software quality assurance procedure for MOOSE-based applications and with software requirements consistent with the ASME NQA-1 standard. Griffin has been adopted into the reactor analysis system for the U.S. NRC and is in use at U.S. companies, universities and national laboratories.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Real-time estimation of the electron temperature profile in DIII-D by leveraging neural-network surrogate models

Control of both the magnitude and the shape of tokamak profiles will be necessary to achieve stable, high-performance plasmas. In order to reject disturbances in real time, feedback-control algorithms rely on accurate real-time knowledge of the plasma state. When diagnostics alone are insufficient, either because they are limited in number or their measurements are too noisy, observers can be used to combine diagnostic data with a response model to provide a better estimation of different plasma properties. An observer has been developed to estimate the electron temperature profile in real time using both diagnostic data from the Thomson scattering system and a model based on the electron heat transport equation describing the evolution of the electron temperature profile. Here, neural network surrogate models are leveraged to help improve the overall model prediction while staying within computation time constraints for real-time use. The observer algorithm is shown in offline tests to produce smooth profiles that are consistent with both the diagnostic data and the electron heat transport equation. When implemented into the real-time plasma control system, this observer will provide valuable information on the electron temperature profile to many potential feedback-control applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Neutron transport methods for multiphysics heterogeneous reactor core simulation in Griffin

Griffin is a reactor physics application based on the Multiphysics Object-Oriented Simulation Environment (MOOSE). This work discloses the methods, algorithms, and implementation for simulating heterogeneous reactor dynamics models. Griffin utilizes a discontinuous finite-element method with discrete ordinates (DFEM-S ) to discretize the field variable of the multigroup neutron transport equation. Multiphysics feedback is handled using two-step tabulated cross-section methodology. Feedback quantities are evaluated using the MOOSE-MultiApp system to couple various engineering phenomena, such as heat conduction and thermal fluids. The multiphysics DFEM-S system is solved using fixed-point iteration with a fully asynchronous parallel sweeper, unstructured coarse-mesh finite difference acceleration, and a multi-timescale improved quasi-static method scheme. The implementation is applied to a multiphysics microreactor model, with two transients: one initiated by a single heat-pipe failure and another by control drum rotation. Importantly, these examples demonstrate the ability of Griffin to tractably solve the neutron transport equation considering seven independent variables and feedback.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Large-eddy simulation and Reynolds-averaged Navier-Stokes modeling of three Rayleigh-Taylor mixing configurations with gravity reversal

In this work, high-fidelity large-eddy simulation (LES) is performed of Rayleigh-Taylor (RT) mixing in three different configurations involving gravity reversal. In each configuration, LES results are compared with one-dimensional Reynolds-averaged Navier-Stokes (RANS) results, and a deficiency in a commonly used transport equation for the mass-flux velocity, $a_j$, is identified. In the first configuration, a classical two-component RT mixing layer is allowed to develop before it is subjected to rapid acceleration reversal. In the second configuration, a three-component RT mixing layer with an intermediate density layer is allowed to develop before being subjected to rapid acceleration reversal. Finally, in the third configuration, a light layer is interposed between two heavy layers; in this configuration, only one interface is RT-unstable at a time as it undergoes rapid acceleration reversal. In all cases, a commonly used buoyancy production closure in the $a_j$ transport equation is shown to lead to significant over-prediction of mixing layer growth after gravity reversal. An alternative formulation for this closure is then presented which is shown to more accurately capture the stabilization effect of gravity reversal.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

A novel energy balance approach for a verifiable and accurate solution of radiation extinction in purely absorbing particle clouds

Here we consider the problem of radiation transport through purely absorbing particle clouds. The gold-standard solution of particle-resolved Monte Carlo ray-tracing method to this problem is computationally expensive and therefore solving the radiation transport equation (specifically, the Beer-Bouguer law (BB-law)) on a Eulerian mesh is often preferred. While the absorption coefficient in the real problem is infinite, the BB-law approximates it to be a finite number in the form of number density through a set of assumptions. For particle clouds that do not obey these assumptions, the BB-law predicts an incorrect exponential decay. Also, when the number density is computed using the nearest-neighbor approach, the BB-law solution diverges when the Eulerian mesh size becomes closer to or smaller than the particle size. This numerical divergence is due to the homogenization error. Although the filtering strategy for number density minimizes the homogenization error, it still converges to an incorrect exponential decay for particle clouds that break the BB-law constraints and the cost of the filtering is equivalent to that of the gold-standard solution. In this study, we develop a novel, highly accurate, verifiable and cost-effective solution to the radiation transport equation on a Eulerian domain using an energy balance approach where we derive an expression for the absorption coefficient as a function of particle and Eulerian mesh sizes. We apply our new method to Poisson and turbulent particle clouds that violate all the BB-law constraints and show that the solution is converging upon the mesh refinement, and eventually, we recover the same gold-standard solution for a much cheaper computational cost.

42 ENGINEERING↗

DOE OSTI2024

A robust spectral element implementation of the $k - τ$ RANS model in Nek5000/NekRS

The $k - ω$ Reynolds Averaged Navier Stokes (RANS) model is one of the industry standard approaches for modeling of turbulent flows. It performs better than the $k - ϵ$ model for low Reynolds number flows and is also more suitable for boundary layers with adverse pressure gradients. Major drawback of the model, however, is that the asymptotic value of $ω$ at the walls is singular, necessitating the use of a contrived “sufficiently” large value for $ω$ as the boundary condition for its transport equation. Here, this invariably leads to the solution being sensitive to near wall grid spacing. While an acceptable solution for low order (finite volume) methods, the excessive near wall gradients lead to persistent numerical stability issues in high order codes. To alleviate the problem, specifically in the context of the high order spectral element code Nek5000, a regularized $k - ω$ approach was formulated in our prior work (Tomboulides et al., 2018). The formulation, however, relies on the use of wall distance and its gradients for modeling the closure terms and can pose problems for simulations in complex geometries. This work presents a novel implementation of the $k - τ$ RANS model in Nek5000, where $τ = 1/ω$, eliminating the need for regularization, owing to the asymptotically bounded behavior of the source terms in the $τ$ transport equation, and also eliminating dependence on wall distance. Robustness and stability of the $k - τ$ model is ensured through implicit treatment of the source terms and their careful numerical implementation and demonstrated through several cases aimed at verification and validation. Studies include both canonical and engineering relevant problems, viz., turbulent channel flow, pipe flow, backward facing step, flow over NACA 0012 airfoil and flow in a T-junction. Results from the $k - τ$ model are shown to be consistent with regularized $k - ω$ model and also with the $k - ω$ SST model in OpenFOAM (for select studies). Comparison with experimental data is also shown, where available, to bolster validation efforts for the $k - τ$ model implementation through prediction of key turbulent quantities of interest.

Nek5000↗

DOE OSTI2021

Assessment of empirical interatomic potential to predict thermal conductivity in ThO 2 and UO 2

Computing vibrational properties of crystals in the presence of complex defects often necessitates the use of (semi-)empirical potentials, which are typically not well characterized for perfect crystals. In this study we explore the efficacy of a commonly used embedded-atomempirical interatomic potential for the U x Th 1- x O 2 system, to compute phonon dispersion, lifetime, and branch specific thermal conductivity. Our approach for ThO 2 involves using lattice dynamics and the linearized Boltzmann transport equation to calculate phonon transport properties based on second and third order force constants derived from the empirical potential and from first-principles calculations. For UO 2 , to circumvent the accuracy issues associated with first-principles treatments of strong electronic correlations, we compare results derived from the empirical interatomic potential to previous experimental results. It is found that the empirical potential can reasonably capture the dispersion of acoustic branches, but exhibits significant discrepancies for the optical branches, leading to overestimation of phonon lifetime and thermal conductivity. The branch specific conductivity also differs significantly with either first-principles based results (ThO 2 ) or experimental measurements (UO 2 ). These findings suggest that the empirical potential needs to be further optimized for robust prediction of thermal conductivity both in perfect crystals and in the presence of complex defects.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗