Search NASA⌕ Search

SEARCH · Search NASA

Results for “RADIATIVE FORCE”

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.

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

At least 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

First β -delayed neutron spectroscopy of 24 O

The β decay of 24 O was investigated at NSCL/MSU using a combination of neutron and γ-ray spectroscopy. For the first time, the β-delayed neutron energy spectrum of 24 O was measured, revealing three intensely populated, isolated neutron-unbound states in 24 F. This allowed for the extraction of the decay strength in 24 F up to 6.2 MeV. A comprehensive comparison of the experimental results with various nuclear theories, ranging from the empirical shell model to the most advanced ab initio calculations, was conducted. While most theoretical predictions align with the experimental data for low-lying states, discrepancies arise at higher excitation energies. Finally, in the transition from 24 O to 24 F, shell model calculations using the empirical USDB interaction predicted the structure of both nuclei without invoking the need for a stronger proton-neutron tensor force, which was postulated for the neighboring isotone 25 F.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Investigating the Effect of Water on the Mechanical Properties of Cellulose from Multiscale Molecular Dynamics Simulations

Classical molecular dynamics (MD) simulations provide insight into the structure and physicochemical properties of materials with atomic resolution. However, the length and time scales accessible to atomistic MD are orders of magnitude smaller than many relevant processes such as the response of a bulk material to experimentally accessible strain rates, which presents challenges when comparing models to experimental measurements. Bottom-up coarse-graining provides a means for systematically mapping atomistic information to lower resolution models to increase the length and time scales achievable by simulation. Cellulose is an abundant carbohydrate biopolymer with applications to many fields of research, such as materials science and renewable energy, due to its desirable mechanical properties and viability for conversion into biofuel. The effect of moisture content on the Young's modulus of cellulose is of special interest due to its native environment often being in the hydrated secondary plant cell wall and the grinding energy requirements for biomass feedstock preprocessing. The current work investigates the effects of water solvent on the Young's modulus of cellulose calculated from coarse-grained MD mechanical stress simulations. The coarse-grained model was parametrized from atomistic MD calculations of cellulose-cellulose potentials of mean force using umbrella sampling techniques under vacuum and solvated conditions. The Young's moduli of the coarse-grained cellulose assemblies parametrized from cellulose in vacuum or solvated in water were computed via mechanical stress simulations to highlight the importance of capturing solvent interactions for modeling the mechanical behavior of cellulose.

BASIC BIOLOGICAL SCIENCES,RADIATION PROTECTION AND↗
DOE OSTI2024

Benchmarking MCNP models to support gamma irradiation experiments

Seismic isolators and damping devices have been used for 30+ years to reduce seismic forces and deployed in 10,000+ structures, including bridges, buildings, and mission-critical infrastructure. Studies have shown that seismic isolation of nuclear power plants can reduce the seismic risk by several orders of magnitude and substantially reduce the overnight capital cost. Advanced reactor developers are considering seismic isolation as an integral design feature in their reactor design. It can either be deployed at the building level, isolating the entire reactor building, or at the component level, isolating individual pieces of safety-class equipment, such as a reactor vessel or a steam generator. For equipment isolation, located close to the reactor vessel, the seismic isolators and/or damping devices may be exposed to gamma and neutron radiation, with a possible effect on their mechanical properties which would affect the response of the isolated equipment under earthquake shaking. The Department of Energy, as part of its Nuclear Energy University Program (NEUP), in collaboration with the Idaho National Laboratory (INL), is funding a project at the University at Buffalo (UB) to investigate the effect of gamma radiation on the mechanical properties of various seismic protective devices. (When isolating equipment, the seismic protective devices will be located outside the reactor vessel where the neutron exposure will be insignificant.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Logistics for the 2026 International Nuclear Accident Dosimetry Intercomparison Exercise

The purpose of this exercise is to test participants’ nuclear accident dosimeters (NADs) and personnel to the requirements established by American National Standards Institute (ANSI), Health Physics Society (HPS) standard ANSI/HPS-N13.3-2013 (R2025) Dosimetry for Criticality Accidents. IER-634 is the second nuclear accident dosimetry exercise taking place at the Armed Forces Radiobiology Research Institute (AFRRI) in Bethesda, MD, a Department of Defense (DOD) facility. This effort continues use of the high-dose neutron exposure test bed maintained by the DOD and represents the continued collaboration between the DOD and Department of Energy (DOE). Previous work at AFRRI included a reactor characterization (IER 484) in 2023 and a nuclear accident intercomparison exercise (IER-602) in 2024.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2025

High-resolution 3D-printed plastic scintillators with tertiary dye

Additive manufacturing offers efficient production of plastic scintillators with nontrivial geometries using vat polymerization, allowing fabrication of geometries which would be difficult or even impossible to produce using conventional subtractive manufacturing. This work presents a novel photocurable scintillator formula that includes coumarin 450 as a tertiary dye to enable high-resolution 3D printing via the manipulation of the 405 nm cure light. Bulk photocured and 3D printed (with and without tertiary dye) samples were compared through observational assessment and spectral response. All samples showed pulse shape discrimination between neutron and gamma events. Inclusion of the tertiary dye has minimal impact on emission spectrum and light output, but significant impact on print resolution as shown by comparison of printed high-complexity geometries and feature resolution test objects. With the use of a cure-limiting dye, unsupported features — such as freestanding pillars — were resolvable down to 0.7 mm. Even finer resolution at or below 0.1 mm was achieved in fully supported, integrated structures printed with off-the-shelf 405 nm desktop 3D printer. Scintillators demonstrated a light output up to 50% of EJ-200 with a PSD figure of merit up to 1.35 at 0.9–1.1 MeVee.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Direct comparison of gyrokinetic and fluid scrape-off layer simulations

Typically, fluid simulations are used for tokamak divertor design. However, fluid models are only valid if the scrape-off layer (SOL) is highly collisional. This assumption is valid in many present-day experiments but is questionable in the upstream SOL of some high-power scenarios envisioned for burning plasmas and fusion pilot plants. This paper reports on comparisons between fluid and kinetic simulations of the SOL for upstream parameters and geometry representative of the Spherical Tokamak for Energy Production fusion pilot plant. The SOLPS-ITER (fluid) and Gkeyll (gyrokinetic) codes are operated in a two-dimensional axisymmetric mode, which replaces turbulence with ad-hoc diffusivities. In kinetic simulations, we observe that the ions in the upstream SOL experience significant mirror trapping. This substantially increases the upstream temperature and has important implications for impurity dynamics. We show that the mirror force, which is excluded in SOLPS’s fluid equations, enhances the electrostatic potential drop along the field line in the SOL. We also show that the assumption of equal main ion and impurity temperatures, which is made in commonly used fluid codes, is invalid for the regimes explored here. The combination of these effects results in superior confinement of impurities to the divertor region in kinetic simulations, consistent with our earlier predictions [Kotschenreuther et al., in 29th IAEA 29 Fusion Energy Conference (IAEA, London, UK, 2023)]. This effect can be dramatic, reducing the midplane impurity density by orders of magnitude. These results indicate that in lower collisionality SOL’s the tolerable downstream impurity densities may be higher than would be predicted by fluid simulations, allowing for higher radiated power while avoiding unacceptable core contamination. Our results highlight the importance of kinetic simulations for divertor design and optimization for fusion pilot plants.

Computational fluid dynamics↗
DOE OSTI2024

Atmospheric Variability Drives Anomalies in the Bering Sea Air–Sea Heat Exchange

Abstract High latitudes, including the Bering Sea, are experiencing unprecedented rates of change. Long-term Bering Sea warming trends have been identified, and marine heatwaves (MHWs), event-scale elevated sea surface temperature (SST) extremes, have also increased in frequency and longevity in recent years. Recent work has shown that variability in air–sea coupling plays a dominant role in driving Bering Sea upper-ocean thermal variability and that surface forcing has driven an increase in the occurrence of positive ocean temperature anomalies since 2010. In this work, we characterize the drivers of the anomalous surface air–sea heat fluxes in the Bering Sea over the period 2010–22 using ERA5 fields. We show that the surface turbulent heat flux dominates the net surface heat flux variability from September to April and is primarily a result of near-surface air temperature and specific humidity anomalies. The airmass anomalies that account for the majority of the turbulent heat flux variability are a function of wind direction, with southerly (northerly) wind advecting anomalously warm (cool), moist (dry) air over the Bering Sea, resulting in positive (negative) surface turbulent flux anomalies. During the remaining months of the year, anomalies in the surface radiative fluxes account for the majority of the net surface heat flux variability and are a result of anomalous cloud coverage, anomalous lower-tropospheric virtual temperature, and sea ice coverage variability. Our results indicate that atmospheric variability drives much of the Bering Sea upper-ocean temperature variability through the mediation of the surface heat fluxes during the analysis period. Significance Statement A long-term ocean warming trend and a recent increase in marine heatwaves in the Bering Sea have been identified. Previous work showed that anomalies in the exchange of heat between the ocean and the atmosphere were the primary driver of Bering Sea temperature variability, but the processes responsible for the heat exchange anomalies were unknown. In this work, we show that the atmosphere is the primary driver of anomalies in the Bering Sea air–sea heat exchange and therefore plays an important role in altering the thermal state of the Bering Sea. Our results highlight the importance of understanding more about how the ocean and the atmosphere interact at high latitudes and how this relationship will be affected by future climate change.

Hayden, Emily E. [a College of Earth, Ocean, and A↗
DOE OSTI2024

Quantum Thermodynamics of Nonequilibrium Processes in Lattice Gauge Theories

A key objective in nuclear and high-energy physics is to describe nonequilibrium dynamics of matter, e.g., in the early Universe and in particle colliders, starting from the standard model of particle physics. Classical computing methods, via the framework of lattice gauge theory, have experienced limited success in this mission. Quantum simulation of lattice gauge theories holds promise for overcoming computational limitations. Because of local constraints (Gauss’s laws), lattice gauge theories have an intricate Hilbert-space structure. This structure complicates the definition of thermodynamic properties of systems coupled to reservoirs during equilibrium and nonequilibrium processes. We show how to define thermodynamic quantities such as work and heat using strong-coupling thermodynamics, a framework that has recently burgeoned within the field of quantum thermodynamics. Our definitions suit instantaneous quenches, simple nonequilibrium processes undertaken in quantum simulators. To illustrate our framework, we compute the work and heat exchanged during a quench in a Z 2 lattice gauge theory coupled to matter in 1+1 dimensions. Here, the thermodynamic quantities, as functions of the quench parameter, evidence a phase transition. For general thermal states, we derive a simple relation between a quantum many-body system’s entanglement Hamiltonian, measurable with quantum-information-processing tools, and the Hamiltonian of mean force, used to define strong-coupling thermodynamic quantities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Overview of ASDEX upgrade results in view of ITER and DEMO

Experiments on ASDEX Upgrade (AUG) in 2021 and 2022 have addressed a number of critical issues for ITER and EU DEMO. A major objective of the AUG programme is to shed light on the underlying physics of confinement, stability, and plasma exhaust in order to allow reliable extrapolation of results obtained on present day machines to these reactor-grade devices. Concerning pedestal physics, the mitigation of edge localised modes (ELMs) using resonant magnetic perturbations (RMPs) was found to be consistent with a reduction of the linear peeling-ballooning stability threshold due to the helical deformation of the plasma. Conversely, ELM suppression by RMPs is ascribed to an increased pedestal transport that keeps the plasma away from this boundary. Candidates for this increased transport are locally enhanced turbulence and a locked magnetic island in the pedestal. The enhanced D-alpha (EDA) and quasi-continuous exhaust (QCE) regimes have been established as promising ELM-free scenarios. Here, the pressure gradient at the foot of the H-mode pedestal is reduced by a quasi-coherent mode, consistent with violation of the high-n ballooning mode stability limit there. This is suggestive that the EDA and QCE regimes have a common underlying physics origin. In the area of transport physics, full radius models for both L- and H-modes have been developed. These models predict energy confinement in AUG better than the commonly used global scaling laws, representing a large step towards the goal of predictive capability. A new momentum transport analysis framework has been developed that provides access to the intrinsic torque in the plasma core. In the field of exhaust, the X-Point Radiator (XPR), a cold and dense plasma region on closed flux surfaces close to the X-point, was described by an analytical model that provides an understanding of its formation as well as its stability, i.e., the conditions under which it transitions into a deleterious MARFE with the potential to result in a disruptive termination. With the XPR close to the divertor target, a new detached divertor concept, the compact radiative divertor, was developed. Here, the exhaust power is radiated before reaching the target, allowing close proximity of the X-point to the target. No limitations by the shallow field line angle due to the large flux expansion were observed, and sufficient compression of neutral density was demonstrated. With respect to the pumping of non-recycling impurities, the divertor enrichment was found to mainly depend on the ionisation energy of the impurity under consideration. In the area of MHD physics, analysis of the hot plasma core motion in sawtooth crashes showed good agreement with nonlinear 2-fluid simulations. This indicates that the fast reconnection observed in these events is adequately described including the pressure gradient and the electron inertia in the parallel Ohm’s law. Concerning disruption physics, a shattered pellet injection system was installed in collaboration with the ITER International Organisation. Thanks to the ability to vary the shard size distribution independently of the injection velocity, as well as its impurity admixture, it was possible to tailor the current quench rate, which is an important requirement for future large devices such as ITER. Progress was also made modelling the force reduction of VDEs induced by massive gas injection on AUG. The H-mode density limit was characterised in terms of safe operational space with a newly developed active feedback control method that allowed the stability boundary to be probed several times within a single discharge without inducing a disruptive termination. Regarding integrated operation scenarios, the role of density peaking in the confinement of the ITER baseline scenario (high plasma current) was clarified. The usual energy confinement scaling ITER98(p,y) does not capture this effect, but the more recent H20 scaling does, highlighting again the importance of developing adequate physics based models. Advanced tokamak scenarios, aiming at large non-inductive current fraction due to non-standard profiles of the safety factor in combination with high normalised plasma pressure were studied with a focus on their access conditions. A method to guide the approach of the targeted safety factor profiles was developed, and the conditions for achieving good confinement were clarified. Based on this, two types of advanced scenarios (‘hybrid’ and ‘elevated’

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Synthesis and characterization of electron beam irradiation-induced damage in polycrystalline metal thin films

High-energy physics research, industrial sterilizing, and material processing depend extensively on electron beam accelerators. Exit windows are crucial components of such electron accelerator systems, maintaining vacuum integrity inside the machine while providing mechanical strength, thermal stability, and radiation resistance at the beam-target interface. In this study, thin metallic films of Ni, Ti, Cr, and V were explored for use in electron beam exit windows, and their properties were compared with the properties of their bulk counterparts. Simulation results of metal foils predicted Ti to exhibit less beam power dissipation compared to Ni. However, Ni possesses superior mechanical and structural properties compared to Ti. The performance of these films under electron beam irradiation was examined by depositing thin layers of these films on silicon and metallic substrates using magnetron sputtering and exposing them to e-beam irradiation in a controlled setup. The deposited films were subjected to a dose of approximately 66 kGy at a beam energy of 10 MeV and characterized prior to and postexposure to the beam using field emission scanning electron microscopy, atomic force microscopy, x-ray diffraction (XRD), and nanoindentation. Particular emphasis was given to characteristics like the grain structure, surface morphology, dislocation density, and hardness. XRD patterns revealed irradiation-induced changes in peak intensities, while the crystallinity remained largely unchanged. Nanoindentation results showed that the pristine and irradiated Ti and Ni films were twice as hard when compared to bulk Ti and Ni samples, regardless of the substrate type used (Si or bulk Ti, or Ni). These results emphasize the superior mechanical properties of thin metal films compared to their bulk counterparts. In conclusion, this study advances the optimization of thin film materials for robust and efficient e-beam applications, ensuring improved durability and operational reliability.

AFM↗
DOE OSTI2025

Investigating Characteristic Droplet Size Distributions in Large Eddy Simulations of Stratocumulus Clouds

Cloud processes relevant to radiative and precipitation properties depend on the shape of the cloud droplet size distribution. Recent holographic observations revealed that cloud droplet populations do not have the same size distribution shapes throughout but form regions of characteristic distributions with similar microphysical properties. We investigate the existence and properties of these characteristic distributions within Large‐Eddy Simulations of stratocumulus clouds using Lagrangian and bin microphysics schemes. Distribution types are identified, revealing localized characteristic distributions that vary on the scale of the largest convective cell for simulations with bin microphysics. The results from the Lagrangian microphysics scheme hint at similar behavior. Compared to observations, the simulated clouds are much more uniform. Analysis of the LES results suggests a connection to the local entrainment rate, so the poorly resolved entrainment interface in LES may be a cause of the uniformity. The uniformity of the large‐scale forcing could also be a factor.

cloud droplet size distributions↗
DOE OSTI2025

Spatial and Temporal Variability of Vertical Velocity under Shallow Cumulus

Vertical velocity distribution below cloud is one of the key determinants of cloud life cycle, but observations of this variable are extremely sparse in space. Doppler lidar retrievals and large-eddy simulations at the U.S. Department of Energy’s Atmospheric Radiation Measurement User Facility Southern Great Plains site are used to determine whether vertical velocity statistics from temporally dense profiles at a single location can be substituted for spatial vertical velocity statistics. We show that even a small number (five) of widely distributed [ O (1°) latitude/longitude spacing] lidars is sufficient sampling to reconstruct domainwide spatial vertical velocity variance, but not higher moments of the vertical velocity distribution. Spatial and temporal vertical velocity variances in the Doppler lidar observations are nearly interchangeable as long as the spatial variance is temporally averaged and the temporal variance is averaged across lidars. This is true even though the dominant spatial scales of vertical velocity variability are ≲ 3 km, more than an order of magnitude smaller than the spacing between the lidars. Further, in the limit where the temporal variance does not vary across a spatial domain (e.g., if the meteorological and surface forcing of the atmospheric turbulence is homogeneous across the domain) and the domain-mean vertical velocity is zero, the commonly available retrieval of temporal vertical velocity variance at one site is equivalent to the spatial variance over the domain. We use an updraft parcel model to show that substituting temporal for spatial vertical velocity statistics will have a relatively minor effect on cloud droplet number concentrations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

An–imidophosphorane (An = U–Pu) bond covalency and proton-coupled electron transfer thermodynamics driven by orbital energy matching

A series of mid-actinide (An = U–Pu) tetrahomoleptic complexes supported by highly electron-donating imidophosphorane ligands, NPC ([NP t Bu(pyrr) 2 ] − , where t Bu = C(CH 3 ) 3 ; pyrr = pyrrolidinyl = N(C 4 H 8 )), are systematically investigated computationally and experimentally to elucidate the nature of actinide–ligand (An–L) covalency across the An 3+/4+/5+ oxidation states. Trends in An–L bonding and redox properties for these complexes, together with their protonated counterparts, are examined using orbital-, electron density-, and energy-decomposition-based methods. This integrated approach reveals progressively improved energy matching between α-spin An 5f and N im 2p orbitals with increasing atomic number and oxidation state, becoming particularly pronounced in the ligand-dominant π-bonding orbitals of An 4+ and An 5+ . In contrast to the An 3+ species, the enhanced An 5f π contributions in the higher-valent counterparts drive the increase in An–N im covalency for later An, thereby inverting the covalency trend to U < Np < Pu. Redistribution of electron density towards the An and N im atomic basins due to the growing energy-matching assisted covalency correlates with higher pKa values and increased N im –H bond dissociation free energies in protonated An 4+ complexes. Electron density at Nim in An 4+ shows a linear correlation with the p K a values calculated via the Bordwell equation. Calculations predict a cathodic shift of 0.84–1.00 V in the redox couples upon protonation, a trend validated when experimentally accessible. These findings demonstrate an increasing role of covalency driven by orbital energy matching from U to Pu in tuning the thermodynamic driving force for proton-coupled electron transfer in the An 5+ species.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2024

Nuclear-matter saturation and symmetry energy within Δ -full chiral effective field theory

Nuclear saturation and the symmetry energy are key properties of low-energy nuclear physics that depend on fine details of the nuclear interaction. The equation of state around saturation is also an important anchor for extrapolations to higher densities and studies of neutron stars. Here we develop a unified statistical framework that uses realistic nuclear forces to link the theoretical modeling of finite nuclei and infinite nuclear matter. We construct fast and accurate emulators for nuclear-matter observables and employ an iterative history-matching approach to explore and reduce the enormous parameter domain of Δ -full chiral interactions. We perform rigorous uncertainty quantification and find that model calibration including O 16 observables gives saturation predictions that are more precise than those that only use few-body data. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Resolving Low Cloud Feedbacks Globally With E3SM High‐Res MMF: Agreement With LES but Stronger Shortwave Effects

This study investigates low cloud feedback in a warmer climate using global simulations from the High-Resolution Multi-scale Modeling Framework (HR-MMF), which explicitly simulates small-scale eddies globally. Two 5-year simulations—one with present-day sea surface temperatures (SSTs) and a second with SSTs warmed uniformly by 4 K—reveal a positive global shortwave cloud radiative effect (SWCRE = 0.3 W/m 2 /K), comparable to estimates from CMIP models. As the climate warms, significant reductions in low cloud cover occur over stratocumulus regions. This study is the first attempt to compare HR-MMF results with predictions from idealized large-eddy simulations from the CGILS intercomparison. Despite different underlying assumptions, we find qualitative agreement in SWCRE and inversion height changes between HR-MMF and CGILS predictions. This suggests reasonable credibility for the CGILS framework in predicting cloud responses under the out-of-sample conditions found in HR-MMF. However, the HR-MMF exhibits stronger SWCRE changes than predicted by CGILS. We explore potential causes for this discrepancy, examining variations in cloud-controlling factors (CCFs) and cloud conditions. Our results show a fairly homogeneous SWCRE response, with little systematic variation tied to the variations in CCFs. This reveals a dominant role for SST forcing in modulating SWCRE.

boundary-layer clouds↗
DOE OSTI2024

Higher-order tails and RG flows due to scattering of gravitational radiation from binary inspirals

Abstract We establish and develop a novel methodology to treat higher-order non-linear effects of gravitational radiation that is scattered from binary inspirals, which employs modern scattering-amplitudes methods on the effective picture of the binary as a composite particle. We spell out our procedure to study such effects: assembling tree amplitudes via generalized-unitarity methods and employing the closed-time-path formalism to derive the causal effective actions, which encompass the full conservative and dissipative dynamics. We push through to a new state of the art for these higher-order effects, up to the third subleading tail effect, at order$$ {G}_N^5 $$ G N 5 and the 5-loop level, which corresponds to the 8.5PN order. We formulate the consequent dissipated energy for these higher-order corrections, and carry out a renormalization analysis, where we uncover new subleading RG flow of the quadrupole coupling. For all higher-order tail effects we find perfect agreement with partial observable results in PN and self-force theories, where available.

Physics↗
DOE OSTI2024

Investigation of Point-Contact Strategies for CFD Simulations of Pebble-Bed Reactor Cores

This study numerically investigated the effects of various contact strategies on the thermal hydraulic behavior within a structured bed of 100 explicitly modeled pebbles. Four contact strategies and two thermal hydraulic conditions were considered. The strategies to avoid contact singularities include decreasing the pebble diameter, increasing the pebble diameter, bridging the pebble surfaces near the contact region, and capping the pebble surfaces near the contact region. One strategy, Strategy 3a, which involves bridging with a cylinder equal to 10% of the pebble diameter, was selected as the baseline strategy because it addressed the contact singularity while minimizing the geometric changes that affect the bed porosity. The two thermal hydraulic conditions were full-power operation (Case 1) and pressurized loss of forced cooling or PLOFC (Case 2). Simulations of the conjugate heat transfer within the structured bed were performed using the Reynolds-averaged Navier–Stokes approach with the realizable k-ϵ turbulence model and two-layer all y + wall treatment. The thermal-fluid quantities of interest were compared between the contact strategies for each case. In Case 1, the hydraulic behavior was sensitive to the contact strategy, with large differences in the pressure drop (30%) and volume-average velocity (4%). The thermal behavior was not sensitive, with less than a 0.5% difference across the strategies. To better understand the separate effects of each heat transfer mode, Case 2 was divided into the following subcases: conduction (Case 2a); conduction/radiation (Case 2b); and conduction/radiation/convection (Case 2c). Case 2a represents an early phase of the PLOFC transient. Case 2b represents an intermediate phase of the PLOFC transient, with the pebble temperatures sufficiently high for the radiative heat transfer to be non-negligible. Case 2c represents a late phase of the PLOFC transient after the establishment of the natural circulation of the heat transfer fluid. For Case 2, large differences in the contact strategy were observed only in Case 2a with only conduction. The difference in the maximum pebble temperature was 23% in Case 2a, 2% in Case 2b, and 0.3% in Case 2c.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Defect-induced phonon-resonant scattering and its influence on thermal transport of irradiated thorium-dioxide

Thermal transport in proton irradiated thorium-dioxide (ThO 2 ) is investigated. Using a combination of experiments and first-principles computational framework, the role of lattice defects on thermal conductivity is analyzed. A resonant-phonon scattering mechanism beyond the traditionally considered Rayleigh scattering is found to significantly influence low-temperature thermal transport in the presence of irradiation-induced point defects. The existence of localized phonon modes associated with irradiation-induced defects is suggested by the inability of the first-principles based thermal conductivity model—which considers only three-phonon interactions and phonon-defects scattering using the Tamura formalism—to predict the experimental results, unless a resonant scattering mechanism is included. The emergence of additional peaks in the Raman spectra in the proximity of phonon-resonant frequency provides further evidence for the existence of localized modes. Coupled with a microstructure evolution model, this analysis enables more accurate analysis for contrasting the contributions of different phonon scattering mechanisms across all irradiation doses and temperatures.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗