Multiple Coulomb scattering of muons in lithium hydride
Not Available
46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
SEARCH · Engineering Papers
Results for “nuclear lattice simulations”
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.
Recovering high-purity uranyl nitrate from simulated used nuclear fuel dissolver solutions by crystallization: rejecting technetium
The separation of U from Tc and other problematic fission product elements like Mo and Ru, along with Sr, Zr, Cs, and Nd, has been achieved via the crystallization of uranyl nitrate hexahydrate (UNH). Rejection of technetium as pertechnetate anion ( 99 TcO 4 – ) is an especially important feature of this system, as it otherwise tends to follow U (VI) in extractive separations. It also raises the salient question regarding why this oxoanion cannot replace nitrate within the crystalline lattice of UNH. Results showed high-yield (>90 %), high-purity (>99 %) recovery of U as UNH from solutions containing 99 TcO 4 – by simple reduction of temperature from 60°C to 20°C. There was no observable interaction of 99 TcO 4 – with UO 2 2+ . The addition of other cations like, Sr 2+ , Zr 4+ , Cs + , and Nd 3+ , also did not form secondary, contaminant solid phases, leaving the > 99 % of the fission product elements in the mother liquor, while the U was recovered at > 90 %. Similarly, Mo and Ru, when added to the mixture, were shown to behave as the other fission-product elements, remaining in the mother liquor during crystallization. As a result, DFT calculations showed that, despite the higher binding strength of TcO 4 – , HMoO 4 – , and BiO 3 – with the UO 2 2+ cation compared to NO 3 – , the hydrogen-bonding network of the two coordinated ions and four waters of hydration in the UNH crystal structure is the driving force for the high specificity of this separation.
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
Detecting the critical point through entanglement in the Schwinger model
Not Available
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Nuclear two point correlation functions on a quantum computer
The calculation of dynamic response functions is expected to be an early application benefiting from rapidly developing quantum hardware resources. The ability to calculate real-time quantities of strongly correlated quantum systems is one of the most exciting applications that can easily reach beyond the capabilities of traditional classical hardware. Response functions of fermionic systems at moderate momenta and energies corresponding roughly to the Fermi energy of the system are a potential early application because the relevant operators are nearly local, and the energies can be resolved in moderately short real time, reducing the spatial resolution and gate depth required. This is particularly the case in quasielastic electron and neutrino scattering from nuclei, a topic of great interest in the nuclear and particle physics communities and directly related to experiments designed to probe neutrino properties. In this paper we use current quantum hardware and error mitigation protocols to calculate response functions for a highly simplified nuclear model through calculations of a 2-point real time correlation function for a modified Fermi-Hubbard model in two dimensions with three distinguishable nucleons on four lattice sites.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
KinCat v.1.0
SAND2024-02099O The software is designed to allow researchers to perform kinetic Monte Carlo (KMC) simulations of catalytic reactions on a 2D lattice. The code is written efficiently to run on a variety of shared memory computing architectures (e.g. GPU, multi-core) and to natively express the full complexity of lateral interactions on reaction rates. The software allows researchers to perform KMC simulations of catalytic reactions on a 2D lattice. It uses parallel shared-memory computing architectures to reduce run-times and allows for simultaneous simulation of multiple independent runs. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.
Najm, Habib↗
Modeling and Simulation of Fuel Burnup in Pebble Bed Reactors
Modeling and simulation of fuel burnup plays important roles in reactor design, operation, safety, and security as well as nuclear material control and accounting (MC&A) [1]. This task is uniquely challenging for pebble bed reactors (PBR) because the pebbles are continuously added and recycled into the reactor core, and their paths through the core are random. To address this problem, we present two simulation models in this paper. Brookhaven National Laboratory (BNL) developed a simple lattice model of a PBR in Serpent software to generate used pebble isotopic concentrations. The benefit of using Serpent software in this specific application is that it helps streamline the data generation process without having to use too many independent software codes in combination to achieve a simple task. For example, transport, burnup and zero power decay can be implemented in a single pass. Three-dimensional core models were developed using Serpent to simulate the burnup process of 5 subject pebbles starting from fresh till they reach nearly target burnup, with each pebble placed in one of the five artificially designated radial channels to capture the changing neutron spectra along the core radius. Equilibrium isotopic concentrations were assumed in all other pebbles in the core. To provide a verification for the Serpent isotope transmutation and decay results, Oak Ridge National Laboratory (ORNL) performed simple SCALE/ORIGEN calculations using the average neutron spectra calculated by Serpent for each of the 5 pebbles. The 252-group neutron spectra from Serpent were then used by ORIGEN to produce the one-group library for depletion and decay calculations. The isotopic concentrations of a few nuclides of interest and neutron and gamma source terms produced from the ORIGEN calculations were compared with the ones from the Serpent calculations. The model simulated in this work was based on the Pebble Bed Modular Reactor (PBMR)-400 design because many data needed for the simulation such as core power profiles, fuel and reflector temperatures, and equilibrium core composition are publicly available. In this paper, we will compare the results between these two approaches and benchmark the results against a set of well-established simulation results for PBMR-400.
Dim, Odera↗
Tensor network simulations of quasi-GPDs in the massive Schwinger model
Generalized parton distribution functions (GPDs) are off-diagonal light-cone matrix elements that encode the internal structure of hadrons in terms of quark and gluon degrees of freedom. In this work, we present the first nonperturbative study of quasi-GPDs in the massive Schwinger model, quantum electrodynamics in 1+1 dimensions (QED 2 ), within the Hamiltonian formulation of lattice field theory. Quasidistributions are spatial correlation functions of boosted states, which approach the relevant light-cone distributions in the luminal limit. Using tensor networks, we prepare the first excited state in the strongly coupled regime and boost it to close to the light-cone on lattices of up to 400 lattice sites. We compute both quasiparton distribution functions and, for the first time, quasi-GPDs, and study their convergence for increasingly boosted states. In addition, we perform analytic calculations of GPDs in the two-particle Fock-space approximation and in the Reggeized limit, providing qualitative benchmarks for the tensor network results. Our analysis establishes computational benchmarks for accessing partonic observables in low-dimensional gauge theories, offering a starting point for future extensions to higher dimensions, non-Abelian theories, and quantum simulations.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Conjectures about the Chiral Phase Transition in QCD from Anomalous Multi-Instanton Interactions
To date numerical simulations of lattice QCD have not found a chiral phase transition of first order that is expected to occur for sufficiently light pions. We show how the restoration of an exact global chiral symmetry can strongly decrease the breaking of the approximate, anomalous U A ( 1 ) symmetry. This is testable on the lattice through simulations for one through four flavors. In QCD a small breaking of the U A ( 1 ) symmetry in the chirally symmetric phase generates novel experimental signals. Published by the American Physical Society 2024
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NuMI Beam Monitoring Simulation and Data Analysis Status and Progress
With the Main Injector Neutrino Oscillation Search (MINOS) experiment decommissioned, muon and hadron monitors became an important diagnostic tool for the NuMI O↵- axis vμ Appearance (NOvA) experiment at Fermilab to monitor the Neutrinos at the Main Injector (NuMI) beam. The goal of this study is to maintain the quality of the monitor signals and to establish correlations with the neutrino beam profile. And we carry out a systematic study of the response of the muon monitors to the changes in the parameters of the proton beam and lattice parameters. We report here on the progress of the beam data analysis and comparison with the simulation results.
43 PARTICLE ACCELERATORS↗
Liouvillian dynamics of the open Schwinger model: String breaking and kinetic dissipation in a thermal medium
Understanding the dynamics of bound state formation is one of the fundamental questions in confining quantum field theories such as Quantum Chromodynamics (QCD). One hadronization mechanism that has garnered significant attention is the breaking of a string initially connecting a fermion and an antifermion. Deepening our understanding of real-time string-breaking dynamics with simpler, lower dimensional models like the Schwinger model can improve our understanding of the hadronization process in QCD and other confining systems found in condensed matter and statistical systems. Here in this paper, we consider the string-breaking dynamics within the Schwinger model and investigate its modification inside a thermal medium, treating the Schwinger model as an open quantum system coupled to a thermal environment. Within the regime of weak coupling between the system and environment, the real-time evolution of the system can be described by a Lindblad evolution equation. We analyze the Liouvillian gaps of this Lindblad equation and the time dependence of the system’s von Neumann entropy. We observe that the late-time relaxation rate decreases as the environment correlation length increases. Moreover, when the environment correlation length is infinite, the system exhibits two steady states, one in each of the sectors with definite charge-conjugation-parity (CP) quantum numbers. For parameter regimes where an initial string breaks in vacuum, we observe a delay of the string breaking in the medium, due to kinetic dissipation effects. Conversely, in regimes where an initial string remains intact in vacuum time evolution, we observe string breaking (melting) in the thermal medium. We further discuss how the Liouvillian dynamics of the open Schwinger model can be simulated on quantum computers and provide an estimate of the associated Trotter errors.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Real-time lattice gauge theory actions: Unitarity, convergence, and path integral contour deformations
Not Available
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
High-temperature structure, elasticity, and thermal expansion of ε-ZrH 1.8
Zirconium hydride is a promising candidate material for nuclear microreactor applications as a solid-state moderator component, owing to its favorable neutronics properties and good thermal stability over other metal hydrides. Here, in the present work, the crystal structure, thermal expansion, and elastic properties of the hydrogen-rich ε phase hydride were measured at elevated temperatures in the range 300–900 K. Samples were prepared by direct hydriding Zircaloy-4 metal – a nuclear-grade zirconium alloy. Room-temperature lattice parameters agree well with those reported from literature for unalloyed zirconium hydride and fall within an observed quadratic H-content dependence. The coefficients of thermal expansion, determined from lattice expansion and dilatometry, agree well within our work but were about 30 % lower than those reported by others for unalloyed hydrides. Density functional theory-based molecular dynamics simulations were used to compare with thermal expansion and elasticity measurements. Results showed lattice parameter temperature dependence and slope of thermal expansion align with those from measurements. Based on diffraction scans at select temperatures, ε phase remained stable in air up to at least 770 K. Likewise, dilatometry showed smooth thermal expansion up to the thermal decomposition temperature around 950 K. The precise decomposition temperature was not determined via diffraction due to sparse scanning. The complete elastic property measurements were gathered for ε-phase Ziracloy-4 hydride for the first time. Young's modulus was lower compared to the metal and δ hydride phases. High-temperature elasticity measurements were limited to <350 K due to acoustic dissipation effects.
36 MATERIALS SCIENCE↗
On the uncertainties in helium generation predictions for fission and fusion alloys
With ongoing advances in fusion and advanced fission reactors, quantifying irradiation effects in materials is critical. Transmutation-induced helium in cladding and structural materials can drive swelling and embrittlement, thereby reducing these components’ lifespans. Yet most studies ignore the considerable uncertainties in predicting helium generation rates. In this work, we created a code wrapper, F-SCATTER, that automatically performs simulations in FISPACT-II. We used this tool to investigate potential variance in helium generation rate, or He/dpa, calculations based on deviations in alloy composition, irradiating neutron flux spectrum, computational methodology, and nuclear data sources. We used 12 wt% Cr HT9 steel as the reference case and observed a 6.5%–98.3% He/dpa spread based on compositional variation within a single chemical specification, a 1.8%–11.5% He/dpa variation upon the incorporation of a 15% artificial uncertainty in flux at each energy, and a He/dpa difference as high as 231% when using ENDF/B-VIII.0 versus TENDL-2021 data libraries. Similar results were found for other prominent iron-based alloys, including Grade 91, castable nano-structured alloy, and 316H—where additional variations exist based on reactor type (e.g. thermal, fast, or fusion) and alloying elements such as carbon, nitrogen, and nickel. Based on the simulated results, we conclude that a significant part of the heat-to-heat variability in swelling responses of Fe-based alloys can be driven by impurity content in alloy compositions, and, therefore, chemical control should be a key element in supply chain design for advanced nuclear energy systems. Furthermore, we provide critical recommendations on best practices for evaluating and reporting helium production and lattice damage rates when computing predictions with multiphysics programs such as FISPACT-II.
FISPACT-II↗
Radiation Response of Ga 2 O 3 MOSFETs Probed via Focused Particle Beams
Gallium Oxide (Ga 2 O 3 ), particularly in its β-phase, is attracting lots of interest for high-power and high-voltage electronics due to its wide bandgap, high breakdown field, and thermal stability. This study investigates the radiation response of Ga 2 O 3 Metal-Oxide-Semiconductor Field-Effect Transistors (MOSFETs) under Total Ionizing Dose (TID) and Displacement Damage (DD) conditions, which are critical for applications in radiation environments. Utilizing a dual-beam Focused Ion Beam and Scanning Electron Microscope setup, microscale analysis of radiation effects on individual devices is performed. The findings reveal distinct behaviors under TID and DD + TID conditions, with TID leading to threshold voltage shifts due to trapped charges, while DD results in decreased drive current attributed to increased carrier scattering from lattice defects. Notably, it is demonstrated that the TID effect can be mitigated through dynamic threshold voltage adjustments and that the predicted TID from ions calculated by Monte Carlo simulations overestimates actual TID due to unaccounted charge yield effects. In conclusion, this research enhances the understanding of Ga2O3 MOSFETs' performance in harsh radiation environments, providing insights for the design of robust electronic devices for space and nuclear applications.
MOSFET↗
Semiclassical treatment of bottomonium suppression and regeneration in 𝑝 + Pb collisions
Here, we study bottomonium suppression in 𝑝 + Pb relative to 𝑝 + 𝑝 collisions at center-of-mass energies of $\sqrt{s_{NN}}$ = 5.02 and 8.16 TeV. Specifically, we combine cold nuclear matter effects (nuclear modifications of the parton densities, energy loss, and momentum broadening) with those from hot nuclear matter (suppression and regeneration) by implementing the formation of a quark-gluon plasma in hydrodynamic simulations. Bottomonium transport in the quark-gluon plasma is evaluated semiclassically, employing two different reaction rates. The first includes quasifree inelastic scattering and gluodissociation employing a perturbative coupling to the medium. The second is based on in-medium 𝑇-matrix calculations where the input potential is constrained by lattice quantum chromodynamics to extract the bottomonium masses and dissociation rates. These semiclassical results are compared to previous calculations in an open quantum system approach and to the experimental data. Predictions for 𝜒 𝑏 suppression at $\sqrt{s_{NN}}$ = 8.16 TeV are also presented.
Physics - Physics of elementary particles and fiel↗
VERA BWR progression problems
During the first phase of the Consortium for Advanced Simulation of Light Water Reactors (CASL) program, the Virtual Environment for Reactor Applications (VERA) was developed with a focus on capabilities for high-fidelity, multiphysics simulation of pressurized water reactors (PWRs). During this development effort, a set of progression problems was created ranging from smaller pin cell calculations to larger 3D full-core calculations. These progression problems helped to guide the development of the software and served as benchmarks against which to test VERA. Since 2019, efforts have been made to extend the capabilities of VERA to model boiling water reactors (BWRs). Because BWR simulations come with many unique challenges, a set of BWR progression problems was developed to aid in this new effort. The BWR progression problems range from 2D lattice calculations to 3D mini-core problems, and reference neutronic solutions were computed using continuous-energy Monte Carlo codes. MPACT, one of the neutronics code in VERA, was benchmarked using the BWR progression problems. The code is capable of computing solutions to all problems. The eigenvalues computed by MPACT agree well with the Monte Carlo reference solutions. (authors)
21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
Fundamental Studies of Tritium Diffusivity in Pure and Defective Zircaloy-4 Getters (Tritium Science Program FY2020 Report)
Zirconium (Zr) alloys are extensively used as tritium (T) getters in nuclear reactors due to their low absorption cross section to thermal neutron, good mechanical properties at high temperature and resistant to corrosion in different environmental conditions. Zr alloys have better thermal properties than many other refractory alloys including stainless steel. The nuclear characteristic is simulated, and T is produced in tritium-producing burnable absorber rods (TPBAR) by the irradiation of neutron flux in pressurized water reactor (PWR). T thus produced diffuses through the pellets and is captured by the getter where it chemically reacts with Zr to form metal hydride (ZrTx). These hydrides are brittle and affect adversely the mechanical properties of alloys. In addition to that, mismatch in the lattice structures after hydride formation creates a stress which needs to be considered in component design and their life evaluation. Therefore, understanding the behavior of T and its species becomes significant as fuel burnup is increased, which leads to increase in hydrogen pickup and oxide formation. The transport of T and its species in pure Zr alloys, and in alloys with certain alloying elements and defects are important to understand in order to enhance the performance of the materials.
36 MATERIALS SCIENCE↗
Shedding light on thermal photon and dilepton production
Electromagnetic radiation from the quark-gluon plasma (QGP) is an important observable to be considered in heavy ion collision experiments. I will provide an update on recent advancements from perturbation theory and quenched lattice simulations. The resummed next-to-leading order (NLO) emission rate has recently been decomposed into transverse and longitudinal components, and extended to non-zero baryon chemical potential. The associated spectral function has also been tested against the Euclidean correlator, for continuum-extrapolated lattice data (at μ B = 0).
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗