Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nuclear Modeling”

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

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

At least 865 records · Page 48

Hardware, methodology and applications of 2 + D backscatter Mössbauer spectroscopy with simultaneous x-ray and γ-ray detection

Abstract A unique method is presented for the acquisition and analysis of 57 Fe backscatter Mössbauer spectra with simultaneous detection of the resonant 14.4 keV γ -rays and the characteristic 6.4 keV x-rays, using a custom-built multi-parameter analyser constructed on the basis of commercial analogue to digital converters and high-speed digital latches. The system allows for the simultaneous registration of Doppler-modulation velocities and photon energies, with up to 4096 and 8192 digital channels respectively. This arrangement is in contrast to most related systems, which detect at a single narrow energy window per detector. Samples of arbitrary atomic structure, morphology and surface topography can be studied without altering the setup or the analysis procedure, provided that the samples are at least micrometre sized. The hardware and software that are used to acquire data with minimal dead time are described and the custom and self-contained methods for post-measurement energy discrimination, background correction and velocity-axis folding are discussed. The data are fit using a general Hamiltonian model for the nuclear energy levels of 57 Fe and a quantum mechanical description of the angular momentum coupling is utilised, with consideration of the crystalline and chemical disorder of the sample under examination. Three examples of distinct magnetic systems, with thicknesses ranging from 5 μ m to 6 mm, that were studied using this method are presented, these are: an amorphous CoFeB-based ribbon with ultra-soft coercivity for high-frequency applications, magnetically hard Nd-Fe-B thick films on Si substrates, examined in both as-deposited and annealed states, and a sample from the nickel-rich iron meteorite NWA 6259 that contains the atomically ordered, elevated coercivity, L 1 0 phase of FeNi, tetrataenite. The wide applicability and usefulness of this method is thus demonstrated on three distinct sample morphologies that required little to no surface preparation prior to examination.

Engineering↗

Assessing the difficulty of capturing the distribution function of neutrinos in neutron star merger simulations

The collision of two neutron stars is a rich source of information about nuclear physics. In particular, the kilonova signal following a merger can help us elucidate the role of neutron stars in nucleosynthesis, and informs us about the properties of matter above nuclear saturation. Approximate modeling of neutrinos remains an important limitation to our ability to make predictions for these observables. Part of the problem is the fermionic nature of neutrinos. By the exclusion principle, the expected value 𝑓 𝜈 for the number of neutrinos in a quantum state is at most 1. Any process producing neutrinos is suppressed by a blocking factor (1 −𝑓 𝜈 ). Recent simulations focused on neutrino physics mostly use a gray two-moment scheme to evolve neutrinos. This evolves integrals of 𝑓 𝜈 over momentum space, preventing direct calculations of blocking factors. Monte Carlo methods may be an attractive alternative, providing access to the full distribution of neutrinos. Their current implementation is, however, inadequate to estimate 𝑓 𝜈 : in our most recent simulations, a single Monte Carlo packet causes, in the worst cases, estimates of 𝑓 𝜈 to jump from 𝑓 𝜈 =0 to 𝑓 𝜈 ∼10 5 . While this is concerning, this brazen violation of the fermionic nature of neutrinos has been largely inconsequential, as the interactions used in simulations avoid direct calculations of 𝑓 𝜈 . We are, however, reaching a level of modeling at which this problem can no longer be ignored. Here, we discuss the relatively simple origin of this issue. We then show that very rough estimates of 𝑓 𝜈 can in theory be obtained in merger simulations, but that they will require a combination of unintuitive weighting schemes for Monte Carlo packets and smoothing of the neutrino distribution at coarser resolution than what the merger simulation uses.

79 ASTRONOMY AND ASTROPHYSICS↗

Inverse Method for Noninvasive Material Discrimination with Multienergy X-Ray Radiography

Active x-ray inspections have great potential for noninvasive detection of illicit materials. Dual-energy x-ray radiography, in particular, has been used to determine material composition by utilizing the energy-dependence of the x-ray attenuation coefficients. However, current implementations of this method are limited in their ability to determine the material composition of composite objects. Here, in this paper, we present an inverse algorithm that uses multienergy x-ray radiography data to noninvasively reconstruct the contents of a container. A critical feature of the current contribution is that material identification can be performed using conventional detectors and x-ray spectra with different endpoint energies produced by varying nothing more than the tube voltage. Adaptive regularization is used to increase the accuracy of material estimations from multienergy data sets. The utility of these methods is demonstrated with experimentally acquired radiographs obtained using a tunable x-ray source that produces spectra with endpoint energies of 100–450 keV. The object inspected is a scale model of a nuclear materials storage container composed of three-dimensional printed plastic and stainless-steel spheres inside a thin-walled steel container. Reconstructions of the steel sphere thicknesses are within a root-mean-square error of 0.37 cm.

47 OTHER INSTRUMENTATION↗

Bulk viscosity from Urca processes: npeμ matter in the neutrino-trapped regime

In this work, we extend our previous study of the bulk viscosity of hot and dense npe matter induced by the Urca process in the neutrino trapped regime to npeμ matter by adding the muonic Urca processes as well as the purely leptonic electroweak processes involving electron-muon transition. Here, the nuclear matter is modeled in a relativistic density functional approach with two different parametrizations which predict neutrino dominated matter (DDME2 model) and anti-neutrino dominated matter (NL3 model) at temperatures for which neutrinos/anti-neutrinos are trapped. In the case of neutrino-dominated matter, the main equilibration mechanism is lepton capture, whereas in the case of antineutrino-dominated matter this is due to neutron decay. We find that the equilibration rates of Urca processes are higher than that of the pure leptonic processes, which implies that the Urca bulk viscosity can be computed with the leptonic reactions assumed to be frozen. We find that the bulk viscosity decreases with temperature as ζ~T –2 at moderate temperatures. At high temperatures this scaling breaks down by sharp drops of the bulk viscosity close to the temperature where the proton fraction is density-independent and the matter becomes scale-invariant. This occurs also when the matter undergoes a transition from the antineutrino-dominated regime to the neutrino-dominated regime where the bulk viscosity attains a local maximum. We also estimate the bulk viscous dissipation timescales and find that these are in the range ≳ 1 s for temperatures above the neutrino trapping temperature. These timescales would be relevant only for long-lived objects formed in binary neutron star mergers and hot proto-neutron stars formed in core-collapse supernovas.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

First-principles determination of the phonon-point defect scattering and thermal transport due to fission products in ThO 2

This work presents the first-principles calculations of the lattice thermal conductivity degradation due to point defects in thorium dioxide using an iterative solution of the Pierels-Boltzmann transport equation. We have used the non perturbative Green’s function methodology to compute the phonon- point defect scattering rates that consider the local distortion around the point defect, including the mass difference changes, interatomic force constants, and structural relaxation near the point defects. The point defects considered in this work include the vacancy of thorium ($V$ $TH$ ) and oxygen ($V$ $O$ ), substitutions of helium ($He$ $Th$ ), krypton ($Kr$ $Th$ ), zirconium ($Zr$ $Th$ ), iodine ($I$ $Th$ ) and xenon ($Xe$ $Th$ ) in the thorium site, and the three different configurations of the Schottky defects. The results of the phonon-defect scattering rate reveal that, among all the considered intrinsic defects, the thorium vacancy and helium substitution in the thorium site scatter the phonon most due to the substantial changes in the force constant and structural distortions. The scattering of phonons due to the substitutional defects unveils that the zirconium atom scatters phonons the least, followed by xenon, iodine, krypton, and helium. This is contrary to the intuition that the scattering strength follows $He$ $Th$ > $Kr$ $Th$ > $Zr$ $Th$ > $I$ $Th$ > $Xe$ $Th$ based on the mass difference. This striking difference in the zirconium phonon scattering is due to the local chemical environment changes. Zirconium is an electropositive element with valency similar to thorium and, therefore, can bond with the oxygen atoms, thus creating less force constant variance compared to iodine, an electronegative element noble gas helium, xenon, and krypton. In conclusion, these results can serve as the benchmark for the analytical models and help the engineering-scale modeling effort for nuclear design.

36 MATERIALS SCIENCE↗

Light element ($$Z=1,2$$) production from spontaneous ternary fission of $$^{252}$$Cf

Abstract The yields of light elements ( $$Z=1,2$$ Z = 1 , 2 ) obtained from spontaneous ternary fission of $$^{252}$$ 252 Cf are treated within a nonequilibrium approach, and the contribution of unstable nuclei and excited bound states is taken into account. These light cluster yields may be used to probe dense matter, and to infer in-medium corrections such as Pauli blocking which is determined by the nucleon density. Continuum correlations are calculated from scattering phase shifts using the Beth-Uhlenbeck formula, and the effect of medium modification is estimated. The relevant distribution is reconstructed from the measured yields of isotopes. This describes the state of the nucleon system at scission and cluster formation, using only three Lagrange parameters which are the nonequilibrium counterparts of the temperature and chemical potentials, as defined in thermodynamic equilibrium. We concluded that a simple nuclear statistical equilibrium model neglecting continuum correlations and medium effects is not able to describe the measured distribution of H and He isotopes. Moreover, the freeze-out concept may serve as an important ingredient to the nonequilibrium approach using the relevant statistical operator concept.

Röpke, G.↗

Microstructure study of U–35 wt.% Zr alloy after quick annealing at 650 °C

Uranium-35wt.%zirconium (U-35wt.%Zr) alloy was annealed for 1 hour and 24 hours at 650 ºC and characterized to understand the early stage microstructure evolution. Dendritic microstructure with fine (~300 nm in length) α-U precipitates clustered between dendrite branches were observed in the 1-hour annealed sample. After 24-hour annealing 650 ºC, the α-U precipitates coarsened, and the dendritic microstructure disappeared due to microstructure homogenization. Furthermore, microchemical homogenization observed with Energy-Dispersive X-ray Spectroscopy analysis suggests that α-U precipitates are approaching thermodynamic equilibrium in the 24-hour annealed sample. The findings from this study have potential impacts on the manufacturing and computer modeling of metallic nuclear fuel.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Exploring Uncertainty in Moment Estimation for Small Earthquakes in Southern Nevada Using the Coda Envelope Method

Compiling source parameter estimates for small earthquakes is important both for our understanding of earthquake physics and for accurately assessing earthquake hazard. Reliable source parameter estimates are difficult to achieve for small earthquakes, in part due to our inability to accurately model the relevant physical processes at high frequencies. The coda envelope methodology developed by Mayeda and Walter (1996) and Mayeda et al. (2003) can mitigate this concern and estimate the moment of small earthquakes by determining the parameters that control the shape of the S-wave coda envelope while eliminating path effects by minimizing the scatter between seismic stations. Here, we use an open-source implementation of this technique called the Coda Calibration Tool (CCT; Barno, 2017) to calculate CCT-based moment magnitude estimates of small earthquakes (M L 0–3) in the Rock Valley, Nevada, region within the Nevada National Security Site. The Rock Valley data set is of particular interest because it allows us to explore the changes in uncertainties of the coda calibration method with earthquake size and depth. We found that a consistent linear relationship exists between the local magnitude M L and our coda-derived M w estimates for earthquakes as small as M L 0–3, but that current CCT workflows do not accurately characterize very shallow events. We also demonstrate that the epistemic uncertainty in the apparent stress value assumed by the CCT algorithm can influence magnitude estimates of small earthquakes. In conclusion, these results provide valuable insight into the seismicity of this region, and inform future analysis and modeling efforts for nuclear monitoring and seismic hazard.

58 GEOSCIENCES↗

Status of the NEAMS and ARC neutronic fast reactor tools integration to the NEAMS Workbench

The Workbench initiative was launched in FY-2017 within the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program to facilitate the transition from conventional tools to high-fidelity tools. The Workbench provides a common user interface for model creation, real-time validation, execution, output processing, and visualization for integrated codes. The integration of the Argonne Reactor Computation (ARC) suite of codes into the NEAMS Workbench was initiated in FY-2017.

97 MATHEMATICS AND COMPUTING↗

RAVEN Theory Manual

RAVEN is a software framework able to perform parametric and stochastic analysis based on the response of complex system codes. The initial development was aimed at providing dynamic risk analysis capabilities to the thermohydraulic code RELAP-7, currently under development at Idaho National Laboratory (INL). Although the initial goal has been fully accomplished, RAVEN is now a multi-purpose stochastic and uncertainty quantification platform, capable of communicating with any system code. In fact, the provided Application Programming Interfaces (APIs) allow RAVEN to interact with any code as long as all the parameters that need to be perturbed are accessible by input files or via python interfaces. RAVEN is capable of investigating system response and explore input space using various sampling schemes such as Monte Carlo, grid, or Latin hypercube. However, RAVEN strength lies in its system feature discovery capabilities such as: constructing limit surfaces, separating regions of the input space leading to system failure, and using dynamic supervised learning techniques. The development of RAVEN started in 2012 when, within the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, the need to provide a modern risk evaluation framework arose. RAVEN’s principal assignment is to provide the necessary software and algorithms in order to employ the concepts developed by the Risk Informed Safety Margin Characterization (RISMC) program. RISMC is one of the pathways defined within the Light Water Reactor Sustainability (LWRS) program. In the RISMC approach, the goal is not just to identify the frequency of an event potentially leading to a system failure, but the proximity (or lack thereof) to key safety-related events. Hence, the approach is interested in identifying and increasing the safety margins related to those events. A safety margin is a numerical value quantifying the probability that a safety metric (e.g. peak pressure in a pipe) is exceeded under certain conditions. Most of the capabilities, implemented having RELAP-7 as a principal focus, are easily deployable to other system codes. For this reason, several side activates have been employed (e.g. RELAP5-3D, any MOOSE-based App, etc.) or are currently ongoing for coupling RAVEN with several different software. The aim of this document is to provide a set of commented examples that can help the user to become familiar with the RAVEN code usage.

97 MATHEMATICS AND COMPUTING↗

Assessment of Fast Reactor Hot Channel Factor Calculation Capability in Griffin and NekRS

The DOE-NE Advanced Modeling and Simulation (NEAMS) program develops, demonstrates, and deploys predictive computational modeling tools for nuclear reactor design and analysis. Under the Fast Reactor Applications activity within the Application Drivers Technical Area of NEAMS, the ability to perform high fidelity “hot channel factor” (HCF) coupled multiphysics simulations has been identified as a valuable simulation capability.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Request for Information: Access to Quantum Systems

Quantum science and technology is a focal point of research at Lawrence Livermore National Laboratory (LLNL). Quantum-coherent devices offer the potential for unprecedented precision in sensing and the ability to directly simulate complex quantum phenomena that have no known efficient classical algorithms. Thus, development and implementation of quantum technologies is expected to have a significant impact on our ability to address some of the most complex national security problems. LLNL maintains projects involving a broad spectrum of Quantum Information Science (QIS) research activities including investigating sources of decoherence, the 3DQ microscope, studies into quantum interconnects, 3-D printed ion traps, fusion energy and nuclear physics simulations, modeling for data-starved environments, and quantum sensing, including accelerometry and timing.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Recent Updates in ETOE-2 and MC 2 -3

The MC 2 -3 code is a multigroup cross section generation code for fast reactor applications, developed by Argonne National Laboratory (ANL) under the DOE Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. In this work, first, the cross section libraries for ENDF/B-VII.1 and ENDF/B-VIII.0 data were generated using the updated ETOE-2 code. Preliminary verification and validation tests of the ENDF/B-VII.1 MC 2 -3 library were performed with the selected fast reactor benchmark problems. Processing and verifying the libraries suggested that significant efforts would be required for thoroughly verifying the ENDF/B-VII.1 MC 2 -3 library and successfully processing the ENDF/B-VIII.0 MC 2 -3 library. Secondly, the cross section generation capability of MC 2 -3 was updated with the intermediate group lattice calculation and the equivalent Dancoff-factor cell (EDC) method to significantly improve the performance of a twodimensional assembly calculation using the method of characteristics (MOC). This effort was made useful in implementing and verifying the EDC method in the Griffin cross section API. Finally, we analyzed a SFR problem, for which unusually large deviations in core eigenvalues from Monte Carlo solutions were reported, by generating multigroup cross sections with MC 2 -3, performing core calculations with DIF3D-VARIANT, and analyzing cross section and eigenvalue results against Monte Carlo solutions. The analysis confirmed that the MC 2 -3/DIF3D solutions were in good agreement with Monte Carlo solutions, providing the correct process of accurately generating broad-group cross sections with MC 2 -3

97 MATHEMATICS AND COMPUTING↗

NEAMS-Multiphysics Technical Assistance in FY21

The Multiphysics Object Oriented Simulation Environment (MOOSE) [1] is a massively parallel finite-element/volume package for multiphysics simulation in science and engineering. The package focuses on providing rapid-development capabilities for engineering applications by leveraging well-built features from libMesh [2] and the Portable Extensible Toolkit for Scientific Computation (PETSc) [3]. Fiscal year 2021 (FY-21) was the first year with funding dedicated to supporting MOOSE-derived applications relevant to the Nuclear Engineering Advanced Modeling and Simulation (NEAMS) program. In this report we outline the work done to support NEAMS applications such as BISON, Griffin, Pronghorn, and System Analysis Module (SAM).

42 ENGINEERING↗

Advanced Measurement and Visualization Techniques for High-Temperature Heat Pipe Experiments

This report conducts a preliminary investigation of the available literature and summarizes some potential advanced measurement and visualization techniques for operating high-temperature heat pipe (HP) that can be implemented to support the United States (U.S.) Department of Energy (DOE) Microreactor Program (MRP). The primary objective for the extensive literature research is to investigate the possibilities and feasibility for the design and construction of an advanced experimental test facility with the aim of producing high-fidelity high-resolution heat pipe data during its operation. Based on the existing HP test facilities at Idaho National Laboratory (INL)—including the Single Primary Heat Extraction and Removal Emulator (SPHERE) and the Microreactor Agile Non-Nuclear Experimental Test Bed (MAGNET)—further efforts will be made to examine the technical feasibility to HP measurements and support the research, development, and demonstration (RD&D) process for an HP-cooled microreactor. Continuing collaborations with the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, other research institutes and universities are considered to be extremely important so that the experimental facility and resultant database can satisfy the validation needs of advanced heat pipe modeling codes being developed under DOE NEAMS program.

42 ENGINEERING↗

Materials and Fuels Complex Operations Management Improvement Strategy for FY2022

The Materials and Fuels Complex (MFC) has experienced substantial growth in terms of staff, research, and production in recent years. MFC operational performance has effectively kept pace with this growth. However, to capture the continuous improvement actions needed to improve effectiveness and increase the efficiency of our management systems, a broad operations management strategy is necessary. The MFC Operations Management Improvement (OMI) Strategy is complementary to the MFC Five-Year Mission and Investment Strategies and the MFC Management Plan. The OMI strategy is structured to address the management systems outlined in the Nuclear Facility Management Standard Operations Model. Selected management systems are evaluated independently in chapters that describe the prior 5 years of performance improvement, a description of improvement actions for the MFC staff that directly perform within or contribute to the management system, process improvements, and any needed equipment improvements. Chapter selection is based on a management system’s need for improvement. In some cases, management systems may be combined in a single chapter. Additional chapters are added to address subject areas not formally described in the standard operations model.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Materials and Fuels Complex Operations Management Improvement Strategy for Fiscal Year 2023

The Materials and Fuels Complex (MFC) has experienced substantial growth in terms of staff, research, and production in recent years. MFC operational performance has effectively kept pace with this growth. However, to capture the continuous improvement actions needed to improve effectiveness and increase the efficiency of our management systems, a broad operations management strategy is necessary. The MFC Operations Management Improvement (OMI) Strategy is complementary to the MFC Five-Year Mission and Investment Strategies and the MFC Management Plan. The OMI strategy is structured to address the management systems outlined in the Nuclear Facility Management Standard Operations Model. Selected management systems are evaluated independently in chapters that describe the prior 5 years of performance improvement, a description of improvement actions for the MFC staff that directly perform within or contribute to the management system, process improvements, and any needed equipment improvements. Chapter selection is based on a management system’s need for improvement. In some cases, management systems may be combined in a single chapter. Additional chapters are added to address subject areas not formally described in the standard operations model.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

NEAMS-Multiphysics MOOSE End Year Framework Activities FY20

The Multiphysics Object Oriented Simulation Environment (MOOSE) is a general finite element package meant for high performance solution of multiphysics problems in science and engineering. In this report, we present additions and enhancements to MOOSE funded by the Nuclear Engineering Advanced Modeling and Simulation (NEAMS) program. NEAMS-funded MOOSE framework library improvements include: expansion of support for multi-level multi-application restart; enhancement of coupling between native MOOSE applications and external libraries; addition of a sparse automatic differentiation (AD) container enabling non-local degree of freedom coupling; block-specific quadrature rules; further development of an eigenvalue executioner; faster setup of periodic boundary conditions; and creation of a mesh meta-data system streamlining simulation startup. Besides developing the framework library, NEAMS funds were used to overhaul a swath of important MOOSE infrastructure in order to substantially improve user experience. These critical infrastructure changes included adapting MOOSE to python 3, improving the test harness to catch race conditions, and most importantly transitioning MOOSE to use Conda, a globally known packaging system, that greatly eases software adoption.

97 MATHEMATICS AND COMPUTING↗