Search NASA⌕ Search

SEARCH · Search NASA

Results for “neutron software”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

The study of C-13(a, n)O-16 reactions in the Majorana Demonstrator calibration data

Neutron captures and delayed decays of reaction products are common sources of backgrounds in ultra-rare event searches. In this work, we studied 13C(a, n)1640 O reactions induced by a-particles emitted within the calibration sources of the Majorana Demonstrator. These sources are thorium-based calibration standards enclosed in carbon-rich materials. The reaction rate was estimated by using the 6129 keV ?-rays emitted from the excited 1643 O states that are populated when the incoming a-particles have enough energy. Thanks to the excellent energy performance of the Demonstrator’s germanium detectors, these characteristic photons can be clearly observed in the calibration data. Facilitated by Geant4 simulations, a comparison between the observed 6129 keV photon rates and predictions by a TALYS-based software was performed. The measurements and predictions were found to be consistent with each other, albeit with large statistical uncertainties. This agreement provides supports for further background predictions for (a, n)-reactions in future double-beta decay search efforts.

6 ≤ A ≤ 19↗
NASA NTRS2021

A NICER View of the Massive Pulsar PSR J0740+6620 Informed by Radio Timing and XMM-Newton Spectroscopy

We report on Bayesian estimation of the radius, mass, and hot surface regions of the massive millisecond pulsar PSR J0740+6620, conditional on pulse-profile modeling of Neutron Star Interior Composition Explorer X-ray Timing Instrument event data. We condition on informative pulsar mass, distance, and orbital inclination priors derived from the joint North American Nanohertz Observatory for Gravitational Waves and Canadian Hydrogen Intensity Mapping Experiment/Pulsar wideband radio timing measurements of Fonseca et al. We use XMM-Newton European Photon Imaging Camera spectroscopic event data to inform our X-ray likelihood function. The prior support of the pulsar radius is truncated at 16 km to ensure coverage of current dense matter models. We assume conservative priors on instrument calibration uncertainty. We constrain the equatorial radius and mass of PSR J0740+6620 to be-+12.390.981.30km and-+2.0720.0660.067Me respectively, each reported as the posterior credible interval bounded by the 16% and 84% quantiles, conditional on surface hot regions that are non-overlapping spherical caps of fully ionized hydrogen atmosphere with uniform effective temperature; a posteriori, the temperature is=-+TlogK5.99100.060.05([])for each hot region. All software for the X-ray modeling framework is open-source and all data, model, and sample information is publicly available, including analysis notebooks and model modules in the Python language. Our marginal likelihood function of mass and equatorial radius is proportional to the marginal joint posterior density of those parameters(within the prior support)and can thus be computed from the posterior samples.

Millisecond pulsars↗
DOE OSTI2020

Benchmark Modeling and Simulation of the FFTF LOFWOS Test #13 Using SAM

The Fast Flux Test Facility (FFTF) was a 400 MW thermal powered, oxide-fueled, liquid sodium cooled test reactor, built to assist development and testing of advanced fuels and materials for fast breeder reactors. In July 1986, a series of unprotected Loss of Flow Without Scram (LOFWOS) transients were performed in FFTF as part of the Passive Safety Testing (PST) program. The LOFWOS Test #13, which was initiated at 50% power and 100% flow with the pump pony motors left off, has been chosen as a benchmark case by IAEA to support collaborative efforts within international partnerships on the validation of simulation tools and models in the area of sodium fast reactor passive safety in an IAEA Coordinated Research Project (CRP), launched in October 2018. The System Analysis Module (SAM) is an advanced and modern system analysis tool under development at Argonne National Laboratory for advanced non-LWR safety analysis. It utilizes the object-oriented application framework MOOSE to leverage the modern software environment and advanced numerical methods. The capabilities of SAM are being extended to enable the transient modeling, analysis, and design of various advanced nuclear reactor systems. To participate the IAEA CRP and enhance the SAM validation base for advanced reactor transient safety analysis, benchmark simulations of the FFTF LOFWOS Test #13 are performed using the SAM code. In this first phase of the validation effort, the thermal-hydraulic behavior of the reactor system is the focus and the reactor kinetics is not considered in the SAM FFTF model. Instead, the results of Argonne’s neutronics calculations are directly used, including the power shape of the active core region and the power history during the transient. The simulation results of FFTF at steady state agreed well with the measured data from the test. During the transient, reasonably good agreement were also obtained. Future work to improve the model will focus on introducing the reactivity predictions into the model, as well as better understanding or resolving the current discrepancies with the measured data.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Validation of the Public Radiation Exposure Calculation for the Incident at the National Institute of Standards and Technology Center for Neutron Research on February 3, 2021

The Department of Energy/National Nuclear Security Administration Consequence Management Program was contacted by the Health Physics Chief of the National Institute of Standards and Technology Center for Neutron Research (NCNR) to review public radiation exposure calculations for an event that occurred on its Gaithersburg, Maryland, campus on February 3, 2021. Subject matter experts from the Nuclear Emergency Support Team (NEST) assets, specifically the Consequence Management Home Team (CMHT) and the National Atmospheric Release Advisory Center (NARAC), were selected to provide support. CMHT used three separate modeling codes to validate the results the scientists at NCNR calculated using the HotSpot model. The analyses were performed using NARAC’s in-house Lagrangian dispersion codes known as LODI and Aeolus, as well as the Turbo FRMAC software from Sandia National Laboratories. The team used parameters provided by the NCNR scientists regarding the site, applicable observable meteorological data, and environmental survey and sampling data to estimate public exposure. Each model estimated public dose at much less than 0.5 mrem. CMHT concurs with the NCNR public radiation exposure calculations which state that members of the public at the 400-meter boundary would have received a radiological dose of less than 0.5 mrem.

61 RADIATION PROTECTION AND DOSIMETRY↗
NASA NTRS2016

SEXTANT X-Ray Pulsar Navigation Demonstration: Flight System and Test Results

The Station Explorer for X-ray Timing and Navigation Technology (SEXTANT) is a technology demonstration enhancement to the Neutron-star Interior Composition Explorer (NICER) mission. NICER is a NASA Explorer Mission of Opportunity that will be hosted on the International Space Station (ISS). SEXTANT will, for the first time, demonstrate real-time, on-board X-ray Pulsar Navigation (XNAV), a significant milestone in the quest to establish a GPS-like navigation capability available throughout our Solar System and beyond. This paper gives an overview of the SEXTANT system architecture and describes progress prior to environmental testing of the NICER flight instrument. It provides descriptions and development status of the SEXTANT flight software and ground system, as well as detailed description and results from the flight software functional and performance testing within the highfidelity Goddard Space Flight Center (GSFC) X-ray Navigation Laboratory Testbed (GXLT) software and hardware simulation environment. Hardware-in-the-loop simulation results are presented, using the engineering model of the NICER timing electronics and the GXLT pulsar simulator-the GXLT precisely controls NASA GSFC's unique Modulated X-ray Source to produce X-rays that make the NICER detector electronics appear as if they were aboard the ISS viewing a sequence of millisecond pulsars.

SEXTANT↗
NASA NTRS2016

SEXTANT X-Ray Pulsar Navigation Demonstration: Flight System and Test Results

The Station Explorer for X-ray Timing and Navigation Technology (SEXTANT) is a technology demonstration enhancement to the Neutron-star Interior Composition Explorer (NICER) mission. NICER is a NASA Explorer Mission of Opportunity that will be hosted on the International Space Station (ISS). SEXTANT will, for the first time, demonstrate real-time, on-board X-ray Pulsar Navigation (XNAV), a significant milestone in the quest to establish a GPS-like navigation capability available throughout our Solar System and beyond. This paper gives an overview of the SEXTANT system architecture and describes progress prior to environmental testing of the NICER flight instrument. It provides descriptions and development status of the SEXTANT flight software and ground system, as well as detailed description and results from the flight software functional and performance testing within the high-fidelity Goddard Space Flight Center (GSFC) X-ray Navigation Laboratory Testbed (GXLT) software and hardware simulation environment. Hardware-in-the-loop simulation results are presented, using the engineering model of the NICER timing electronics and the GXLT pulsar simulator-the GXLT precisely controls NASA GSFC's unique Modulated X-ray Source to produce X-rays that make the NICER detector electronics appear as if they were aboard the ISS viewing a sequence of millisecond pulsars

Pulsar↗
DOE OSTI2024

Investigating resource-efficient neutron/gamma classification ML models targeting eFPGAs

There has been considerable interest and resulting progress in implementing machine learning (ML) models in hardware over the last several years from the particle and nuclear physics communities. A big driver has been the release of the Python package, hls4ml, which has enabled porting models specified and trained using Python ML libraries to register transfer level (RTL) code. So far, the primary end targets have been commercial field-programmable gate arrays (FPGAs) or synthesized custom blocks on application specific integrated circuits (ASICs). However, recent developments in open-source embedded FPGA (eFPGA) frameworks now provide an alternate, more flexible pathway for implementing ML models in hardware. These customized eFPGA fabrics can be integrated as part of an overall chip design. In general, the decision between a fully custom, eFPGA, or commercial FPGA ML implementation will depend on the details of the end-use application. In this work, we explored the parameter space for eFPGA implementations of fully-connected neural network (fcNN) and boosted decision tree (BDT) models using the task of neutron/gamma classification with a specific focus on resource efficiency. We used data collected using an AmBe sealed source incident on Stilbene, which was optically coupled to an OnSemi J-series silicon photomultiplier (SiPM) to generate training and test data for this study. We investigated relevant input features and the effects of bit-resolution and sampling rate as well as trade-offs in hyperparameters for both ML architectures while tracking total resource usage. The performance metric used to track model performance was the calculated neutron efficiency at a gamma leakage of 10 -3 . The results of the study will be used to aid the specification of an eFPGA fabric, which will be integrated as part of a test chip.

47 OTHER INSTRUMENTATION↗
NASA NTRS2007

High Fidelity, Fuel-Like Thermal Simulators for Non-Nuclear Testing: Analysis and Initial Test Results

Non-nuclear testing can be a valuable tool in the development of a space nuclear power system, providing system characterization data and allowing one to work through various fabrication, assembly and integration issues without the cost and time associated with a full ground nuclear test. In a non-nuclear test bed, electric heaters are used to simulate the heat from nuclear fuel. Testing with non-optimized heater elements allows one to assess thermal, heat transfer. and stress related attributes of a given system, but fails to demonstrate the dynamic response that would be present in an integrated, fueled reactor system. High fidelity thermal simulators that match both the static and the dynamic fuel pin performance that would be observed in an operating, fueled nuclear reactor can vastly increase the value of non-nuclear test results. With optimized simulators, the integration of thermal hydraulic hardware tests with simulated neutronic response provides a bridge between electrically heated testing and fueled nuclear testing. By implementing a neutronic response model to simulate the dynamic response that would be expected in a fueled reactor system, one can better understand system integration issues, characterize integrated system response times and response characteristics and assess potential design improvements at relatively small fiscal investment. Initial conceptual thermal simulator designs are determined by simple one-dimensional analysis at a single axial location and at steady state conditions; feasible concepts are then input into a detailed three-dimensional model for comparison to expected fuel pin performance. Static and dynamic fuel pin performance for a proposed reactor design is determined using SINDA/FLUINT thermal analysis software, and comparison is made between the expected nuclear performance and the performance of conceptual thermal simulator designs. Through a series of iterative analyses, a conceptual high fidelity design is developed: this is followed by engineering design, fabrication, and testing to validate the overall design process. Test results presented in this paper correspond to a "first cut" simulator design for a potential liquid metal (NaK) cooled reactor design that could be applied for Lunar surface power. Proposed refinements to this simulator design are also presented.

Bragg-Sitton, Shannon M.↗
DOE OSTI2021

A NICER View of the Massive Pulsar PSR J0740+6620 Informed by Radio Timing and XMM-Newton Spectroscopy

We report on Bayesian estimation of the radius, mass, and hot surface regions of the massive millisecond pulsar PSR J0740+6620, conditional on pulse-profile modeling of Neutron Star Interior Composition Explorer X-ray Timing Instrument event data. We condition on informative pulsar mass, distance, and orbital inclination priors derived from the joint North American Nanohertz Observatory for Gravitational Waves and Canadian Hydrogen Intensity Mapping Experiment/Pulsar wideband radio timing measurements of Fonseca et al. We use XMM-Newton European Photon Imaging Camera spectroscopic event data to inform our X-ray likelihood function. The prior support of the pulsar radius is truncated at 16 km to ensure coverage of current dense matter models. We assume conservative priors on instrument calibration uncertainty. We constrain the equatorial radius and mass of PSR J0740+6620 to be ${12.39}_{-0.98}^{+1.30}$ km and ${2.072}_{-0.066}^{+0.067}$ M ⊙ respectively, each reported as the posterior credible interval bounded by the 16% and 84% quantiles, conditional on surface hot regions that are non-overlapping spherical caps of fully ionized hydrogen atmosphere with uniform effective temperature; a posteriori, the temperature is ${\mathrm{log}}_{10}(T\,[{\rm{K}}])={5.99}_{-0.06}^{+0.05}$ for each hot region. All software for the X-ray modeling framework is open-source and all data, model, and sample information is publicly available, including analysis notebooks and model modules in the Python language. Furthermore, our marginal likelihood function of mass and equatorial radius is proportional to the marginal joint posterior density of those parameters (within the prior support) and can thus be computed from the posterior samples.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Verification of Spent Fuel Inside Dry Storage Casks by Cask Top Fast Neutron Mapping (FY2023 Mid-Year Report)

This project is developing a prototype scanner array verification system for detection of missing fuel assemblies in spent-fuel storage casks. The prototype consists of six fast-neutron scintillator detectors mounted to a linear actuator frame that is placed on the top of a spent fuel cask to scan across all fuel assembly positions. The scanner array was assembled and tested at LLNL in FY2022. A field test schedule has been requested at the Idaho National Laboratory (INL) Cask Farm site for FY2023. Note that the Cask Farm contractor determines this scheduling and not INL directly. Further system automation will be designed and implemented with the goal of obtaining a level of system operation that meets IAEA needs. This includes integration of the scanner array and data-acquisition control software into a single interface for operator use. In addition, commercial operators and the IAEA may have special requirements for portability, shipping, lifting, and installation. Prior to the Field Test at INL, the system will be operated at LLNL to exercise lifting procedure and linear actuators, monitor stability of detector energy and pulse-shape discrimination calibration, and test system software integration efforts. Following the Field Test, we will present results and discuss the technology with the IAEA. We will incorporate additional improvements to the system based on lessons learned from the field test and feedback from the IAEA. If successful, the technology can be transferred to the IAEA or other stakeholders for assessment.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRSDate not supplied

Pragmatic Stress Prediction on Additively Manufactured Coupons

Prediction of residual stresses from process parameters for additively manufactured large metal parts is computationally expensive. NASA is currently developing meter-scale parts with direct energy deposition. Practically, the predictive computational methods need to efficiently scale-up to meter-scale parts. Coupled thermal-mechanical multi-physics simulations have been developed with the pragmatic method using ABAQUS, COMSOL Multiphysics, ALE3D software. The residual stresses are a result of the manufacturing process which creates thermal cycling of the build layers. The pragmatic method uses lumped thermal layers for stress predictions to reduce computational costs. The stress predictions as well as deformations of the different codes are compared with each other and with ANSYS Additive using identical material models, boundary and initial conditions. The codes were used to simulate three different geometries: a thin wall, hollow cylinder and twin-cantilever part. The coupon parts were then manufactured with Inconel-625. The residual stresses in these parts were measured using X-ray diffraction as well as neutron beam diffraction at NIST. The stress measurements for the two technologies are compared. The pragmatic stress prediction method enabled predictions of the multi-centimeter scale parts using desktop computer workstations in only a few hours for each coupon. The results of the simulated stress predictions compared favorably with the measured stresses even though thermally lumped layers were employed. Finally, a two-meter scale nozzle was simulated using ANSYS Additive. The simulations were used to examine the build orientation trade-space with respect to resulting geometric deformation. The predicted deformations were compared to measurements of an actual subscale part manufactured with direct energy deposition.

Additive Manufacturing↗
NASA NTRSDate not supplied

Pragmatic Stress Prediction on Additively Manufactured Coupons

Prediction of residual stresses from process parameters for additively manufactured large metal parts is computationally expensive. NASA is currently developing meter-scale parts with direct energy deposition. Practically, the predictive computational methods need to efficiently scale-up to meter-scale parts. Coupled thermal-mechanical multi-physics simulations have been developed with the pragmatic method using ABAQUS, COMSOL Multiphysics, ALE3D software. The residual stresses are a result of the manufacturing process which creates thermal cycling of the build layers. The pragmatic method uses lumped thermal layers for stress predictions to reduce computational costs. The stress predictions as well as deformations of the different codes are compared with each other and with ANSYS Additive using identical material models, boundary and initial conditions. The codes were used to simulate three different geometries: a thin wall, hollow cylinder and twin-cantilever part. The coupon parts were then manufactured with Inconel-625. The residual stresses in these parts were measured using X-ray diffraction as well as neutron beam diffraction at NIST. The stress measurements for the two technologies are compared. The pragmatic stress prediction method enabled predictions of the multi-centimeter scale parts using desktop computer workstations in only a few hours for each coupon. The results of the simulated stress predictions compared favorably with the measured stresses even though thermally lumped layers were employed. Finally, a two-meter scale nozzle was simulated using ANSYS Additive. The simulations were used to examine the build orientation trade-space with respect to resulting geometric deformation. The predicted deformations were compared to measurements of an actual subscale part manufactured with direct energy deposition.

pragmatic method↗
NASA NTRS2013

MESSENGER at Mercury: Early Orbital Operations

The MErcury Surface, Space ENvironment, GEochemistry, and Ranging (MESSENGER) spacecraft, launched in August 2004 under NASA's Discovery Program, was inserted into orbit about the planet Mercury in March 2011. MESSENGER's three flybys of Mercury in 2008-2009 marked the first spacecraft visits to the innermost planet since the Mariner 10 flybys in 1974-1975. The unprecedented orbital operations are yielding new insights into the nature and evolution of Mercury. The scientific questions that frame the MESSENGER mission led to the mission measurement objectives to be achieved by the seven payload instruments and the radio science experiment. Interweaving the full set of required orbital observations in a manner that maximizes the opportunity to satisfy all mission objectives and yet meet stringent spacecraft pointing and thermal constraints was a complex optimization problem that was solved with a software tool that simulates science observations and tracks progress toward meeting each objective. The final orbital observation plan, the outcome of that optimization process, meets all mission objectives. MESSENGER's Mercury Dual Imaging System is acquiring a global monochromatic image mosaic at better than 90% coverage and at least 250 m average resolution, a global color image mosaic at better than 90% coverage and at least 1 km average resolution, and global stereo imaging at better than 80% coverage and at least 250 m average resolution. Higher-resolution images are also being acquired of targeted areas. The elemental remote sensing instruments, including the Gamma-Ray and Neutron Spectrometer and the X-Ray Spectrometer, are being operated nearly continuously and will establish the average surface abundances of most major elements. The Visible and Infrared Spectrograph channel of MESSENGER's Mercury Atmospheric and Surface Composition Spectrometer is acquiring a global map of spectral reflectance from 300 to 1450 nm wavelength at a range of incidence and emission angles. Targeted areas have been selected for spectral coverage into the ultraviolet with the Ultraviolet and Visible Spectrometer (UVVS). MESSENGER's Mercury Laser Altimeter is acquiring topographic profiles when the slant range to Mercury's surface is less than 1800 km, encompassing latitudes from 20 deg. S to the north pole. Topography over the remainder of the southern hemisphere will be derived from stereo imaging, radio occultations, and limb profiles. MESSENGER's radio science experiment is determining Mercury's gravity field from Doppler signals acquired during frequent downlinks. MESSENGER's Magnetometer is measuring the vector magnetic field both within Mercury's magnetosphere and in Mercury's solar wind environment at an instrument sampling rate of up to 20 samples/s. The UVVS is determining the three-dimensional, time-dependent distribution of Mercury's exospheric neutral and ionic species via their emission lines. During each spacecraft orbit, the Energetic Particle Spectrometer measures energetic electrons and ions, and the Fast Imaging Plasma Spectrometer measures the energies and mass per charge of thermal plasma components, both within Mercury's magnetosphere and in Mercury's solar-wind environment. The primary mission observation sequence will continue for one Earth year, until March 2012. An extended mission, currently under discussion with NASA, would add a second year of orbital observations targeting a set of focused follow-on questions that build on observations to date and take advantage of the more active Sun expected during 2012-2013. MESSENGER's total primary mission cost, projected at $446 M in real-year dollars, is comparable to that of Mariner 10 after adjustment for inflation.

mercury↗
NASA NTRS2012

MESSENGER at Mercury: Early Orbital Operations

The MErcury Surface, Space ENvironment, GEochemistry, and Ranging (MESSENGER) spacecraft, launched in August 2004 under NASA's Discovery Program, was inserted into orbit about the planet Mercury in March 2011. MESSENGER's three flybys of Mercury in 2008-2009 marked the first spacecraft visits to the innermost planet since the Mariner 10 flybys in 1974-1975. The unprecedented orbital operations are yielding new insights into the nature and evolution of Mercury. The scientific questions that frame the MESSENGER mission led to the mission measurement objectives to be achieved by the seven payload instruments and the radio science experiment. Interweaving the full set of required orbital observations in a manner that maximizes the opportunity to satisfy all mission objectives and yet meet stringent spacecraft pointing and thermal constraints was a complex optimization problem that was solved with a software tool that simulates science observations and tracks progress toward meeting each objective. The final orbital observation plan, the outcome of that optimization process, meets all mission objectives. MESSENGER's Mercury Dual Imaging System is acquiring a global monochromatic image mosaic at better than 90%coverage and at least 250 m average resolution, a global color image mosaic at better than 90%coverage and at least 1 km average resolution, and global stereo imaging at better than 80%coverage and at least 250 m average resolution. Higher-resolution images are also being acquired of targeted areas. The elemental remote sensing instruments, including the Gamma-Ray and Neutron Spectrometer and the X-Ray Spectrometer, are being operated nearly continuously and will establish the average surface abundances of most major elements. The Visible and Infrared Spectrograph channel of MESSENGER's Mercury Atmospheric and Surface Composition Spectrometer is acquiring a global map of spectral reflectance from 300 to 1450 nm wavelength at a range of incidence and emission angles. Targeted areas have been selected for spectral coverage into the ultraviolet with the Ultraviolet and Visible Spectrometer (UVVS). MESSENGER's Mercury Laser Altimeter is acquiring topographic profiles when the slant range to Mercury's surface is less than 1800 km, encompassing latitudes from 201S to the north pole. Topography over the remainder of the southern hemisphere will be derived from stereo imaging, radio occultations, and limb profiles. MESSENGER's radio science experiment is determining Mercury's gravity field from Doppler signals acquired during frequent downlinks. MESSENGER's Magnetometer is measuring the vector magnetic field both within Mercury's magnetosphere and in Mercury's solar wind environment at an instrument sampling rate of up to 20 samples/s. The UVVS is determining the three-dimensional, time-dependent distribution of Mercury's exospheric neutral and ionic species via their emission lines. During each spacecraft orbit, the Energetic Particle Spectrometer measures energetic electrons and ions, and the Fast Imaging Plasma Spectrometer measures the energies and mass per charge of thermal plasma components, both within Mercury's magnetosphere and in Mercury's solar-wind environment. The primary mission observation sequence will continue for one Earth year, until March 2012. An extended mission, currently under discussion with NASA, would add a second year of orbital observations targeting a set of focused follow-on questions that build on observations to date and take advantage of the more active Sun expected during 2012-2013. MESSENGER's total primary mission cost, projected at $446 M in real-year dollars, is comparable to that of Mariner 10 after adjustment for inflation.

Planetary Science↗
DOE OSTI2021

Fast Emulation of Expensive Simulations using Approximate Gaussian Processes [Slides]

Nuclear Computational Low-Energy Initiative (NUCLEI) collaboration uses Density Functional Theory (DFT) simulations to predict the structure and binding energies of nuclei over a wide range of proton (Z) and neutron (N) numbers. The DFT simulations utilize a particular parameterization of a Skyrme energy density functional called UNEDF1 which depends on 12 free parameters that must be fit to data (M Kortelainen et al 2014). Fitting involves comparing (e.g.) predicted binding energies of nuclei to experimentally measured values. We use only binding energies as observables, but DFT with UNEDF1 will predict structure (shape) observables as well. In this work, assessing the capability of approximate GP emulators to balance emulator accuracy with computational speed to facilitate improved UNEDF1 calibration. Sparse GPs are straightforward to train and accurate. Calibration is not straightforward with MCMC (using MH or HMC/NUTS). We produced reusable software for continuing and building on this work as well as accessing and using Darwin cluster compute resources

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Coupling between TOUGH3 and FLAC3D

Sandia National Laboratories has hired Itasca Consulting Group, Inc., the authors of the FLAC3D geomechanics software, to couple FLAC3D with TOUGH3, the porous media flow solver. The work is being done to enable a coupled mechanical-thermal-hydraulic analysis of a potential criticality event in a dual purpose cannister (DPC). The U.S. Department of Energy Office of Spent Fuel and Waste Science & Technology is investigating the performance of DPCs for direct geological disposal of spent nuclear fuel. Post closure criticality control is an important aspect of this investigation. Over geological timescales, it is envisioned that the canister and canister overpack will develop fractures due to stress corrosion processes. A breach in the canister could allow groundwater to fill the canister. Fresh water is a neutron moderator; thus, if the canister internals and fuel assemblies have been sufficiently degraded, a criticality event could occur. Such an event would release enough energy to boil the water between the fuel rods and pressurize the cannister. This internal pressurization may cause the initial fractures in the canister and overpack to grow. It is important to understand the change in hydraulic transmissivity between the canister and surroundings for two reasons: first, because it may control the potential for and frequency of subsequent criticality events; second, because it will control the release of radionuclides from the canister. The motivation for this work is to better understand the potential for periodic criticality events, cannister damage, and release of radionuclides during a criticality event in a DPC.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Sierra/SolidMechanics for Target Pulse Analysis

Sierra is a Multiphysics engineering simulation software suite developed and maintained by Sandia National Laboratories. Created for nuclear weapon needs, it is available to national laboratories with a foreign national restriction. The STS Target Systems analysis team is interested in using Sierra to perform explicit dynamics simulations of the tantalum clad tungsten target blocks as a response to sudden temperature changes induced by proton beam pulses. The capability already exists with Abaqus|Explicit, however, there are time and cost limitations due to licensing. The following document serves as a “quick-start” guide to the successful setup, simulation, and post-processing of a target model in Sierra on the OZ3 computing cluster at the Spallation Neutron Source.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

A Full-Induction Magnetohydrodynamics Solver for Liquid Metal Fusion Blankets in Vertex-CFD

Multiphysics modeling of liquid metal fusion blankets, which produce tritium and convert energy of neutrons created via fusion reactions into heat, is crucial for predicting performance, ensuring structural integrity, and optimizing energy production. While traditional blanket modeling of liquid metal flows during normal steady operating conditions commonly employs the inductionless approximation of the magnetohydrodynamics (MHD) equations, transient scenarios, when the plasma-confining magnetic field varies on millisecond time scales, require a full-induction MHD approach that dynamically evolves the magnetic field via the time-dependent induction equation. This paper presents the formulation, implementation, and initial verification of a full-induction MHD solver integrated within the open-source Vertex-CFD framework, which aims to achieve tight multiphysics coupling, a flexible software design enabling easy extension and addition of physics models, and performance portability across computing platforms. The solver utilizes finite element spatial discretization, implicit Runge–Kutta time integration, and an inexact Newton method to solve the resulting discrete nonlinear system, leveraging Trilinos packages for efficient computation. Verification against selected benchmark problems demonstrates accuracy and robustness of the solver. Furthermore, when the solver is applied to an idealized blanket model in 2.5D and full 3D, results obtained with Vertex-CFD are in good agreement with recently published quasi-2D simulations. These findings establish a computational foundation for future simulations of transient MHD phenomena in liquid metal blankets with Vertex-CFD, and open avenues for future extensions and performance optimizations.

Endeve, Eirik [ORNL] (ORCID:0000000312519507)↗