Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nuclear structure data”

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 865 records · Page 48

Thermal Scattering Law Data Development for Paraffin Wax

Paraffin wax is often used as a nuclear moderator to slow down the fast neutrons in experimental critical assemblies [1]. It is a colorless and soft solid material that consists primarily of straight-chain alkanes (n-alkanes), which are hydrocarbons with the general formula CnH2n+2 [2-3]. The length of the hydrocarbon chain ranges from C20 to C30 and higher [2]. It is distinguished by its solid state at room temperature and begins to melt above approximately 310 K [4]. Paraffin wax is a commonly employed substance in the manufacture of shielding. One of its noteworthy characteristics is its ability to effectively absorb the neutrons. Also, it possesses a high macroscopic cross section, which enables it to efficiently moderate neutrons. As a result, paraffin wax is extensively utilized in various applications where moderation and shielding of neutrons are needed. For simulations, it is necessary to evaluate its thermal scattering law (TSL) and cross sections. Computationally, classical molecular dynamics (CMD) simulations provide the capability of simulating atomic details. For example, several unary, binary, and few multi component mixtures have been investigated of the paraffin model by using molecular dynamics simulations [5-12]. An assessment of thermal neutron scattering in a heavy paraffinic oil treated both as a solid and a viscous fluid containing 25% linear branched paraffin (C30H62), 35% one ring cycloalkane (C30H60), 15% two rings cycloalkane (C30H58), and 25% aromatic (C30H60) chains has been studied using CMD simulations for producing TSL data [13]. Nevertheless, there is lack of TSL and cross section data for paraffin wax as most of the reported analyses focus on the unary and binary mixture of n-alkanes, which is not consistent with actual paraffin wax [2]. In this work, we applied the equilibrium CMD simulations technique to explore the structure and dynamical properties of wax, which are fundamental input to calculate the TSL. A paraffin wax system was modeled using the CMD code LAMMPS (Large-scale Atomic/Molecular Massively Parallel Simulator) [14-15] with the semi-empirical COMPASS [16] force field. The density of state (DOS) was calculated from the normalized velocity autocorrelation function (VACF), which is the Fourier transform of the normalized VACF. The DOS was used for the calculation of the TSL and thermal scattering cross sections. The paraffin wax atomic system was constructed by using the MedeA material design platform [17], and was benchmarked using available properties (i.e., density, bond lengths, angles, diffusivity, and viscosity).

Nuclear Criticality Safety Program (NCSP)↗

Impact of structural biology and the protein data bank on us fda new drug approvals of low molecular weight antineoplastic agents 2019–2023

Abstract Open access to three-dimensional atomic-level biostructure information from the Protein Data Bank (PDB) facilitated discovery/development of 100% of the 34 new low molecular weight, protein-targeted, antineoplastic agents approved by the US FDA 2019–2023. Analyses of PDB holdings, the scientific literature, and related documents for each drug-target combination revealed that the impact of structural biologists and public-domain 3D biostructure data was broad and substantial, ranging from understanding target biology (100% of all drug targets), to identifying a given target as likely druggable (100% of all targets), to structure-guided drug discovery (>80% of all new small-molecule drugs, made up of 50% confirmed and >30% probable cases). In addition to aggregate impact assessments, illustrative case studies are presented for six first-in-class small-molecule anti-cancer drugs, including a selective inhibitor of nuclear export targeting Exportin 1 (selinexor, Xpovio), an ATP-competitive CSF-1R receptor tyrosine kinase inhibitor (pexidartinib,Turalia), a non-ATP-competitive inhibitor of the BCR-Abl fusion protein targeting the myristoyl binding pocket within the kinase catalytic domain of Abl (asciminib, Scemblix), a covalently-acting G12C KRAS inhibitor (sotorasib, Lumakras or Lumykras), an EZH2 methyltransferase inhibitor (tazemostat, Tazverik), and an agent targeting the basic-Helix-Loop-Helix transcription factor HIF-2α (belzutifan, Welireg).

60 APPLIED LIFE SCIENCES↗

Probing elemental speciation in hydrochar produced from hydrothermal liquefaction of anaerobic digestates using quantitative X-ray diffraction

Valorization of hydrochar, a solid byproduct from hydrothermal liquefaction (HTL) of anaerobically-digested agriculture wastes (digestates), requires fundamental knowledge of elemental speciation. This study investigated the effects of reaction temperatures (320–360 °C), digestate pH (3.5–8), and digestate cellulose-to-lignin ratios (0.2–1.8) on the speciation (chemical form) and composition of organics and inorganics in hydrochars produced during hydrothermal treatment. Quantitative X-ray diffraction (XRD) method was the primary technique used to characterize hydrochars. The comprehensive XRD pattern processing including the Rietveld refinement protocols demonstrated that the organic phase was comprised of mostly crystalline monocyclic, heterocyclic, and polycyclic aromatics with diverse aliphatic and aromatic substituents, while the inorganic mineral phase consisted of calcium-phosphates, magnesium-phosphates, calcium-carbonates, and magnesium-carbonates. Further, XRD results were validated by the elemental yields of products and the distribution of chemical functionalities measured using solid-state nuclear magnetic resonance (NMR) spectroscopy. The characterization data were used to evaluate proposed mechanistic pathways using compositional analysis of biocrude and aqueous-phase coproducts. Mechanistic pathways developed in the study suggested that benzoic acids, phenols, benzaldehydes, phenolic aldehydes, α-dicarbonyls, and α-hydroxycarbonyls were responsible for the precipitation of organics through various reactions depending on operating conditions. Meanwhile, the formation of inorganic compounds appeared to be consistently represented by reactions including dehydration, hydrolysis, endergonic reduction, and structure rearrangement of native minerals in the digestates. This study provides basic knowledge needed to create and assess potential elemental speciation pathways. In addition, the results of the study facilitate the specification of process conditions to optimize targeted utilization routes of hydrochar for more economically-feasible and sustainable HTL processing.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Floating zone crystal growth, structure, and properties of a cubic Li 5.5 La 3 Nb 1.5 Zr 0.5 O 12 garnet-type lithium-ion conductor

As a promising candidate for solid-state electrolytes in Li-ion batteries, the garnet-type Li-ion conductor series Li 5+x La 3 Nb 2-x Zr x O 12 (LLNZO) (0 ≤ x ≤ 2) exhibits high ionic conductivity at room temperature. However, no previous single-crystal growth or characterization has been reported for LLNZO compositions 0 ≤ x ≤ 1. To obtain a complete understanding of the trend in the structure–property relationship in this class of materials, we used the floating zone (FZ) method to grow a single crystal of Li 5.5 La 3 Nb 1.5 Zr 0.5 O 12 that was 4 mm in diameter and 10 mm in length. Using Laue neutron single-crystal diffraction, two distinct Li sites were observed: tetrahedral 24d and octahedral 96h sites. The maximum entropy method (MEM) based on neutron single-crystal diffraction data was used to map Li nuclear density and estimate that the bottleneck of Li transport exists between neighboring tetrahedral and octahedral sites, and that Li is delocalized between split octahedral sites. Room-temperature Li-ion conductivity in Li 5.5 La 3 Nb 1.5 Zr 0.5 O 12 measured with electrochemical impedance spectroscopy (EIS) was 1.37 × 10 -4 S cm -1 . The Li migration activation energy was estimated to be 0.50 eV from EIS and 0.47 eV from dielectric relaxation measurements. The Li-ion jump attempt rate was estimated to be 1.47 × 10 12 Hz while the time scale of successful migration is 10 -7 to 10 -6 s.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Degradation of REBCO coated conductors due to a combination of epoxy impregnation, thermal cycles, and quench: Characteristics and a method of alleviation

REBa2Cu3Ox (REBCO) coated conductors are an enabling technology for >1 GHz nuclear magnetic resonance spectrometers that need >23.5 T solenoids, compact high-field fusion reactors beyond the International Thermonuclear Experimental Reactor, and 20 T class accelerator magnets. An in-depth understanding of stress and strain management is critically needed for these applications. Here, we report experimental data and an analysis of damages of REBCO-coated conductors in several scenarios relevant for their applications in high-field magnets and unique to their multi-layer thin film structure. The examined scenarios include failure induced by quench, failure induced by epoxy impregnation, and failure induced by a combination of epoxy impregnation and quench. The root cause is film blistering and peeling induced by thermally interfacial stresses characteristic of multilayered heteroepitaxial thin film structures on thick substrates utilized. The quantitative stress analysis is introduced using simple and easy to use formulas and clearly indicates how materials and structural characteristics affect the magnitude and the distribution of stresses and deflections. Finally, a simple method is introduced for removing degradation of REBCO coated conductors due to epoxy impregnation. Our data and analysis help our understanding of the practical limit of REBCO-coated conductors and can be utilized as a guidance for the physical design of REBCO coated conductors and high-current cables assembled from them and for minimizing damages.

Yin, Shijian (ORCID:0000000224312274)↗

Spatially-Resolved Characterization Techniques and Their Implications for Nuclear Debris Formation.

Debris from nuclear tests has a complex formation process and can inform multiple fields of study such as geology, atmospheric science, shock physics, and chemistry under extreme conditions. Macroscopic nuclear debris forms as the result of fireball interaction with surrounding materials (e.g., structural and environmental component), which rapidly undergo melting and vaporization followed by condensation, convective/diffusive mixing, and solidification over the course of seconds. In atmospheric events, some of this material may disperse over long distances, but much of the material is deposited close-in to ground zero, often in multicomponent, partially (or entirely) amorphous debris formations. During the U.S. nuclear testing program, this material was collected and analyzed, particularly for radionuclide composition. Fallout formation models were developed from nuclear test data based on such radionuclide compositional analyses. These models were not just used to understand fallout dispersion (i.e. the spread of radioactivity over geographical regions) but also guide radiochemical interpretations of historical nuclear tests.6 Developments in analytical techniques over the past several decades have made it possible to make new analyses on historical debris, some of which may be more than a half-century old. For example, inductively-coupled plasma mass spectrometry (ICP-MS) has been used to look at the trace elements in glassy fallout material from aboveground nuclear tests (including the Trinity test). Advanced analyses using X-ray absorption allowed for the measurement of oxidation state of fallout constituents, including actinides. The ability to measure trace actinides and their subsequent oxidation state is important to understanding how the surrounding environment may have influenced the resultant fallout composition and may have implications for why certain fractionation trends have been observed. Here we present data from historic nuclear test debris illustrating the power of spatially resolved methods to connect interaction of the near-field environment with the explosion and provide new insights into nuclear debris formation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

β decay of Ba 141

The β-decay strength function of nuclides produced in fission is important as it dictates the distribution of decay energy between electrons, neutrinos, and γ rays and so is critical for calculating decay heat in reactors and for estimating the reactor antineutrino spectrum. Several experimental techniques are available to determine this strength function, including electron spectroscopy, γ-ray calorimetry (TAGS spectroscopy), and detailed, high-resolution spectroscopy with modern large high-purity germanium arrays. This work investigates the decay of the well-known and strongly produced fission fragment 141 Ba. A beam of 141 Cs was implanted at the target position of the Gammasphere and the subsequent decay of the daughter 141 Ba was studied. Extensive decay spectroscopy was possible up to the decay Q value of 3.197(7) MeV, including a significant extension of the level scheme and detailed angular correlation measurements for all levels with greater than 0.25% β feeding. The distribution of the β-decay strength was then inferred and compared to previous calorimetric studies. In conclusion, the agreement was excellent and provides a benchmark for comparing strength function methods and data for a more detailed understanding of the structure of 141 La.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

An accelerated framework for predicting creep rupture lifetimes in engineering alloys

Confidently predicting high-temperature deformation, including creep and creep rupture, is paramount for the design and commercialization of candidate materials for advanced nuclear energy systems. To accelerate creep quantification, we introduce a framework that enables rapid, cost-effective, and reliable prediction of creep rupture lifetimes, minimizing reliance on time-intensive bulk creep testing. Unlike conventional creep analysis, which requires extensive time and resources, our method leverages a maximum of four short-term bulk creep tests as training data for prediction. This framework combines high-throughput nanoindentation up to 700 °C with these targeted bulk tests to inform our creep rupture model in order to predict rupture lifetimes. The strong agreement between our predictions and conventional experimental data demonstrates the effectiveness of our approach for accelerated creep analysis and lifetime prediction of structural components in high-temperature applications. Our multi-pronged approach motivates further integration of computational tools and advanced instrumentation to establish a universal framework for understanding high-temperature material responses.

36 MATERIALS SCIENCE↗

Nonmonotonic-potential description of polarization effects, fusion, and nuclear rainbows in elastic scattering of 6 Li + 12 C at 4.5–600 MeV

The experimental differential cross-section (CS) and analyzing power (AP) data of the 6 Li + 12 C elastic scattering over a wide laboratory energy scale (4.5MeV≤𝐸 lab ≤600MeV) are analyzed within the framework of the optical model (OPM) using nonmonotonic (NM) nucleus–nucleus potentials. The real part of the NM potentials is derived from the Pauli-embodied energy density-functional (EDF) formalism with the sudden approximation. The real part of the noncentral spin-orbit and tensor terms, as well as the imaginary parts, are treated phenomenologically. The effect of the radius of sensitivity on the CS and AP data is found to be more important at lower energies. The diffractive and refractive scattering with Airy structures in the whole angular region of the elastic scattering across the studied energy range is successfully described within the OPM using the NM 6 Li + 12 C potential. The near- and far-side (N and F) decomposition of the total elastic-scattering amplitudes has also been studied using our NM potentials. The evolution of the Airy minima in the angular distributions, coupled with the fitting of the AP data, provides an accurate description of Airy minima of different orders. The OPM calculations with the NM potentials describe exceptionally well the CS, vector analyzing power (VAP), and tensor analyzing power data at 𝐸 lab =9.0,19.24,20,30, and 50MeV. In agreement with our past successful descriptions of CS and the opposite signs of the VAP data for the 6 Li and 7 Li elastic scattering using NM potentials in OPM, the present results appear to provide a better fit, so far, than those obtained from the coupled-channels method. The fusion cross sections of 6 Li + 12 C have been predicted in the energy range (4.5MeV≤𝐸 lab ≤20MeV), fitting the experimental data well in the range 𝐸 lab =2.97–11.77MeV. The EDF potential without any energy dependence and renormalization is also found to describe satisfactorily the experimental CS and AP data at energies up to several hundreds of MeV.

6 ≤ A ≤ 19↗

Demonstration of NEAMS Multiphysics Tools for Fast Reactor Applications

The SHARP toolkit is a high-fidelity reactor simulation tool developed under the U.S. Department of Energy, Office of Nuclear Energy Advanced Modeling and Simulation (NEAMS) Campaign. SHARP toolkit is comprised of the neutronics module PROTEUS thermal hydraulics module Nek5000, and structural mechanics module Diablo. During FY17 and FY18, the PROTEUS and Nek5000 components of SHARP were applied to solve challenging sodium-cooled fast reactor (SFR) problems. In particular, selected hot channel factors (HCF) for a prototype metal-fueled SFR design (the AFR-100) were analyzed in high fidelity, and the “SHARP zooming capability” for SFRs was developed and demonstrated to reduce computational expense for full core problems in cases where detailed data is needed in selected fuel assemblies. After the previous success applying SHARP to challenging SFR problems, the focus in FY19 and FY20 expanded to additional fast reactor applications including lead cooled fast reactors (LFR) and sodium cooled fast reactors (SFR). The specific technical tasks were (1) assessment of hot channel factors for LFR, for which no data currently exists, and (2) demonstration of zooming capability in assemblies of the Versatile Test Reactor (VTR). First-of-a-kind hot channel factor (HCF) estimation for LFR with high fidelity codes (PROTEUS/Nek5000) was successfully demonstrated in this study which began in FY19 and continued in FY20. Selected HCF were computed and compared with SFR data (AFR-100, EBR-II). The findings confirm that different reactor types, design parameters and uncertainties lead to different HCFs. Careful estimation of HCF for a specific design is necessary to obtain appropriate HCFs. In addition to improvement in HCF accuracy, high fidelity tools generate data to help the designer better understand the mechanism of the impact from these uncertainties. For example, the impact of cladding thickness manufacturing tolerance resulted in non-intuitive effects in the corner pins of the LFR assembly. This procedure of computing HCF using high fidelity models shows promise and flexibility for being repeated for any arbitrary reactor of choice. Along with the application on SFR and LFR, the capability of the tools has also been matured to deal with different reactor types and designs. Progress was made towards extending the previously demonstrated SHARP zooming capability to non-fueled SFR assemblies. In particular, in FY19 a gamma transport capability was implemented in both high fidelity PROTEUS solvers in order to accurately account for heat deposition caused by gamma particles, which accounts for ~10% of total core power. Neutronics verification cases were carried out for a candidate Versatile Test Reactor (VTR) design using the new gamma transport capability in PROTEUS. Comparisons were made with continuous energy MCNP calculations and shown to agree well. The models for the full core design with heterogeneous control and fuel assemblies is in progress for PROTEUS-SN and completed with MCNP. The MCNP power distributions were transferred to Nek5000 to perform thermal hydraulic calculations of the control and fuel assembly.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

30 cm Drop Modeling

This report documents modeling and analysis of the 30 cm cask drop scenario for spent nuclear fuel (SNF) casks. The 30 cm drop is an important hypothetical dynamic loading scenario for SNF because 10 CFR 71.71 defines testing requirements for normal conditions of transport (NCT) of SNF packages, and free drop (of 30 cm for the heaviest packages) tests are required. Other NCT conditions specified in 10 CFR 71.71 include vibration normally incident to transport, and DOE’s multimodal transportation test (MMTT) and related analysis concluded that mechanical shock and vibration incident to transportation is negligible for SNF structural integrity because the loads are so small. The SNF cask 30 cm drop has a lot more potential to cause damage to SNF because the kinetic energy and momentum of the system is stopped abruptly during the cask impact event. Test data and model predictions both agree that the SNF assemblies are subjected to significantly higher loads in a cask drop event than normal transportation shock and vibration, but currently there is no reason to expect gross structural damage of the fuel assembly or widespread SNF cladding failure as a result of a 30 cm cask drop. The ultimate goal of this work is to close the knowledge gap of stress profiles on SNF rods and fuel assembly components during NCT, and this report describes the progress toward reaching that goal as well as a plan for finishing this work over the next year or two. Key remaining tasks include refining the fuel assembly model, calculating a selection of specific fuel assembly cases to populate the damage model, and finally running the damage model to predict the full general response range of 30 cm drop loads on SNF.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DUNE and MINERvA Flux Studies and a Measurement of the Charged-Current Quasielastic Antineutrino Scattering Cross Section with $<E_{\nu}> \sim 6 $ GeV on a CH Target

One of the biggest open questions in physics is the prevalence of matter over anti-matter in the present universe. One of possible answers lies in the violation of charge and parity symmetry in the lepton sector that would favor matter over anti-matter so that the universe becomes dominated by matter over time. Recent results from neutrino experiments like T2K and NOvA have indicated that there is CP violation in the leptonic sector \citep{Abe:2019vii}. However, the precise measurement of CP violation will require an unprecedented level of accuracy in these oscillation experiments. One of the challenges lies with the systematic uncertainties that come with the measurements, which include uncertainties related to the neutrino flux modeling, interaction modeling and detector related systematic uncertainties.\\ The first part of this thesis will go through the challenges related to flux uncertainties. It will discuss the flux uncertainties related to hadron production in the LBN F beamli ne for the DUNE experiment. Hadron production models used in flux simulations vary widely depending on the simulation and choice of physics models used to simulate the hadron production. This thesis will explain the method of using existing hadron production data to constrain these uncertainties, a method which has been used in MINERvA experiment.\\ The second part of the thesis goes through the effect of possible flux mismodeling and MINERvA's approach to address this effect. In doing so, we discovered that the suspected mismodeling could also result from incorrect estimation of the energy scale of the muons. This project not only demonstrated the challenges of flux modeling but also led to a novel approach of using neutrino energy spectra to understand the correlation between detector and neutrino beamline parameters. \\ A high statistics measurement of the anti-neutrino scattering cross section in the charged current quasielastic (CCQE) channel is the final part of the th esis. Th e 2-body interactions in this channel have the advantage of allowing reconstruction of the interaction kinematics from the outgoing muon trajectory which can be measured very well. Double differential cross sections as a function of muon kinematics, one of the deliverables of this analysis, will help future oscillation experiments understand their data. The cross-section as a function of the four-momentum transferred from the leptonic system to the hadronic system ($Q^{2}_{QE}$) can be used to test models used in simulating antineutrino interactions and the nuclear effects that can modify the predicted cross sections. Nuclear effects arise from the complex nuclear environment and both modify the initial scattering and change the fate of the particles produced from in neutrino-nucleus interactions. Study of nuclear effects will help to understand the structure of the atomic nucleus and its impact on neutrino interactions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Machine learning-enabled prediction of chemical durability of A 2 B 2 O 7 pyrochlore and fluorite

Pyrochlore-structure type and its derivative in a general formula A 2 B 2 O 7 (A = rare earth elements and actinides; B = Ti, Sn, Zr, Hf, Pb, Si, etc.) display excellent structural flexibility and rich crystal chemistry as promising nuclear waste form materials capable of immobilizing actinides and fission products. It is essential to understand these materials’ chemical durability and element release of radionuclides in order to evaluate their performance in near-field environment. However, it is a formidable grand technological challenge to experimentally perform durability testing across hundreds of thousands of possibilities resulting from their extreme compositional complexities due to cation substitutions at both A and B-sites. In this work, we demonstrate a machine learning approach to determine the key materials parameters and structural characteristics governing the leaching behaviors from a small set of selected compositions as model systems, enabling a science-based prediction of their chemical durability that can be extended to a wide range of chemical compositions. The combination of four key structural characteristics and materials parameters, including ionic radius size difference , ionic potential difference , electronegativity difference , and lattice parameter , creates features an optimized prediction of the chemical durability. Two machine learning models, linear regression and Kernel ridge regression models, are trained on the randomly-split training dataset derived from the experimentally-determined elemental release rates, and subsequently tested on the testing dataset. The predicted leaching rates from both machine learning models show an excellent agreement with the experimental data, demonstrating the feasibility of rapidly evaluating the material properties of new compositions. These results highlight the immense potential of synergizing informatics through machine learning-based models and well-controlled experiments of selected model systems to accelerate materials design and discovery with optimized compositions and performance of promising materials for effective nuclear waste management.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Assessment and Improvement of Fission Product Transport Predictions of Particle Fuel in BISON

The U.S. Department of Energy’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) program aims to develop predictive capabilities by applying computational methods to the analysis and design of advanced reactor and fuel cycle systems. This program has been providing engineering-scale support for the development of BISON, a high-fidelity and high-resolution fuel performance tool. This study was motivated by the need to incorporate more physics-based models in BISON in order to foster tri-structural isotropic (TRISO) applications. This document details the integration of new modeling capabilities in BISON, including (1) development of pyrolytic carbon (PyC) and silicon carbide (SiC) layer anisotropic thermal and mass transport capabilities, (2) verification of the mass diffusion solution in TRISO modeling, (3) calibration of fission product diffusivity using Advanced Gas Reactor (AGR) experiments, (4) improved fission product release modeling by developing compact diffusion modeling capabilities, and (5) documentation of accelerated failure analysis on the BISON website. Improvements made to the diffusion models and parameters were documented and validated against AGR-1 and -2 experiment data.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Summary of US DOE R&D Activities on Graphite Oxidation (2006–2021)

The objective of the international collaboration between United States Department of Energy (U.S.-DOE) and Generation IV International Forum (GIF) is the development of the next generation of nuclear energy systems. The current GIF Project Arrangement (PA) on Materials (2018-2022) was revised in 2019 and extended for another 10 years (2020-2030). The Work Package 1 (“Graphite”) of the extended Project Plan (PP) on Materials specifies technical tasks and High Level Deliverables for research and development (R&D) activities related to using graphite in fuel elements, reflectors, and support structures of Very High Temperature Reactors (VHTR). The graphite tasks include specification and acquisition, qualification and development of new grades, characterization of properties, and development of behavior models. Specifically, Task 1.4 (“Graphite Oxidation Behavior”) outlines planned activities related to acute oxidation by air and chronic oxidation by impurities in the helium coolant. A final report on experimental data regarding graphite oxidation behavior is scheduled for 2022 (deliverable 3.1.1.4.a). In preparation of this deliverable, this document summarizes the R&D activities funded by U.S.-DOE from 2006 (the inception of the VHTR system arrangement) through present (2021). This report is being submitted to the GIF Graphite Working Group (GWG) to serve as input for the GWG high-level deliverable to the Project Management Board (PMB) of PA on Materials. Besides U.S.-DOE, other organizations participating to Task 1.4 of the current PA on Materials are: European Commission’s Joint Research Center (JRC), Korea Atomic Energy Research Institute (KAERI), and Japan Atomic Energy Agency (JAEA). U.S.-DOE is the main contributor on graphite oxidation R&D, with 85 % commitment of total funding during 2018-2022.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Machine learning molecular dynamics simulations toward exploration of high-temperature properties of nuclear fuel materials: case study of thorium dioxide

Predicting materials properties of nuclear fuel compounds is a challenging task in materials science. Their thermodynamical behaviors around and above the operational temperature are essential for the design of nuclear reactors. However, they are not easy to measure, because the target temperature range is too high to perform various standard experiments safely and accurately. Moreover, theoretical methods such as first-principles calculations also suffer from the computational limitations in calculating thermodynamical properties due to their high calculation-costs and complicated electronic structures stemming from f-orbital occupations of valence electrons in actinide elements. Here, we demonstrate, for the first time, machine-learning molecular-dynamics to theoretically explore high-temperature thermodynamical properties of a nuclear fuel material, thorium dioxide. The target compound satisfies first-principles calculation accuracy because f-electron occupation coincidentally diminishes and the scheme meets sampling sufficiency because it works at the computational cost of classical molecular-dynamics levels. We prepare a set of training data using first-principles molecular dynamics with small number of atoms, which cannot directly evaluate thermodynamical properties but captures essential atomistic dynamics at the high temperature range. Then, we construct a machine-learning molecular-dynamics potential and carry out large-scale molecular-dynamics calculations. Consequently, we successfully access two kinds of thermodynamic phase transitions, namely the melting and the anomalous λ transition induced by large diffusions of oxygen atoms. Furthermore, we quantitatively reproduce various experimental data in the best agreement manner by selecting a density functional scheme known as SCAN. Our results suggest that the present scale-up simulation-scheme using machine-learning techniques opens up a new pathway on theoretical studies of not only nuclear fuel compounds, but also a variety of similar materials that contain both heavy and light elements, like thorium dioxide.

74 ATOMIC AND MOLECULAR PHYSICS↗

Independent normalization for γ -ray strength functions: The shape method

Here, the shape method, a novel approach to obtain the functional form of the γ-ray strength function (γSF), is introduced. In connection with the Oslo method the slope of the nuclear level density (NLD) and γSF can be obtained simultaneously even in the absence of neutron resonance spacing data. The foundation of the shape method lies in the primary γ-ray transitions which preserve information on the functional form of the γSF. The shape method has been applied to 56 Fe, 92 Zr, and 164 Dy, which are representative cases for the variety of situations encountered in typical NLD and γSF studies. The comparisons of results from the shape method to those from the Oslo method demonstrate that the functional form of the γSF is retained regardless of nuclear structure details or J π values of the states fed by the primary transitions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Quantum Monte Carlo Calculations of Magnetic Form Factors in Light Nuclei

Here, we present Quantum Monte Carlo calculations of magnetic form factors in A = 6-10 nuclei, based on Norfolk two- and three-nucleon interactions, and associated one- and two-body electromagnetic currents. Agreement with the available experimental data for 6 Li, 7 Li, 9 Be and 10 B up to values of momentum transfer q ~ 3 fm -1 is achieved when two-nucleon currents are accounted for. We present a set of predictions for the magnetic form factors of 7 Be, 8 Li, 9 Li, and 9 C. In these systems, two body currents account for ~ 40-60% of the total magnetic strength. Measurements in any of these radioactive systems would provide valuable insights on the nuclear magnetic structure emerging from the underlying many-nucleon dynamics. A particularly interesting case is that of 7 Be, as it would enable investigations of the magnetic structure of mirror nuclei

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗