Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “radiation physics”

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

Efficacy of the symmetry-adapted basis for ab initio nucleon-nucleus interactions for light- and intermediate-mass nuclei

We study the efficacy of a new ab initio framework that combines the symmetry-adapted (SA) no-core shell-model approach with the resonating group method (RGM) for unified descriptions of nuclear structure and reactions. We obtain ab initio neutron-nucleus interactions for 4 He, 16 O, and 20 Ne targets, starting with realistic nucleon-nucleon potentials. We discuss the effect of increasing model space sizes and symmetry-based selections on the SA-RGM norm and direct potential kernels, as well as on phase shifts, which are the input to calculations of cross sections. We demonstrate the efficacy of the SA basis and its scalability with particle numbers and model space dimensions, with a view toward ab initio descriptions of nucleon scattering and capture reactions up through the medium-mass region.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

FITEVT: A FORTRAN program for arrival-time analysis of nuclear-decay events

FITEVT is a FORTRAN program designed to analyze recorded arrival times of nuclear-decay events affected by dead-time. Here, the primary goal of the program is precise and accurate determination of half-lives. Its method involves imposition of a known sufficiently long extending dead-time to the recorded event sequence, so that the original dead-time effects are completely obliterated and the remaining live-times of the survived events are known exactly. Upon completion of the analysis, the arrival-time spectrum and the live-time spectrum are predicted and compared to those constructed from the survived events.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Methodology for physics-informed generation of synthetic neutron time-of-flight measurement data

Accurate neutron cross section data are a vital input to the simulation of nuclear systems for a wide range of applications from energy production to national security. The evaluation of experimental data is a key step in producing accurate cross sections. There is a widely recognized lack of reproducibility in the evaluation process due to its artisanal nature and therefore there is a call for improvement within the nuclear data community. This can be realized by automating/standardizing viable parts of the process, namely, parameter estimation by fitting theoretical models to experimental data. This automation effort could greatly benefit from a synthetic data resource. This work leverages problem-specific physics, Monte Carlo sampling, and a general methodology for data synthesis to generate unlimited, labelled experimental cross-section data that is statistically indistinguishable to the observed data. Heuristic and, where applicable, rigorous statistical comparisons to observed data support this claim. The demonstration is based on/limited to transmission measurements at Rensselaer Polytechnic Institute (RPI) and energy-differential cross sections in the resolved resonance region (RRR). An open-source software is published alongside this article that executes the complete methodology to produce high-utility synthetic datasets. The goal of this work is to provide an approach and corresponding tool that will allow the evaluation community to begin exploring more data-driven, ML-based solutions to long-standing challenges in the field.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

High-order finite element method for atomic structure calculations

We introduce featom, an open source code that implements a high-order finite element solver for the radial Schrödinger, Dirac, and Kohn-Sham equations. The formulation accommodates various mesh types, such as uniform or exponential, and the convergence can be systematically controlled by increasing the number and/or polynomial order of the finite element basis functions. The Dirac equation is solved using a squared Hamiltonian approach to eliminate spurious states. Here, to address the slow convergence of the $κ=±1$ states due to divergent derivatives at the origin, we incorporate known asymptotic forms into the solutions. We achieve a high level of accuracy (10 -8 Hartree) for total energies and eigenvalues of heavy atoms such as uranium in both Schrödinger and Dirac Kohn-Sham solutions. We provide detailed convergence studies and computational parameters required to attain commonly required accuracies. Finally, we compare our results with known analytic results as well as the results of other methods. In particular, we calculate benchmark results for atomic numbers (Z) from 1 to 92, verifying current benchmarks. We demonstrate significant speedup compared to the state-of-the-art shooting solver dftatom. An efficient, modular Fortran 2008 implementation, is provided under an open source, permissive license, including examples and tests, wherein particular emphasis is placed on the independence (no global variables), reusability, and generality of the individual routines.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

SIMULATeQCD: A simple multi-GPU lattice code for QCD calculations

The rise of exascale supercomputers has fueled competition among GPU vendors, driving lattice QCD developers to write code that supports multiple APIs. Moreover, new developments in algorithms and physics research require frequent updates to existing software. These challenges have to be balanced against constantly changing personnel. At the same time, there is a wide range of applications for HISQ fermions in QCD studies. This situation encourages the development of software featuring a HISQ action that is flexible, high-performing, open source, easy to use, and easy to adapt. In this technical paper, we explain the design strategy, provide implementation details, list available algorithms and modules, and show key performance indicators for SIMULATeQCD, a simple multi-GPU lattice code for large-scale QCD calculations, mainly developed and used by the HotQCD collaboration. The code is publicly available on GitHub.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

The TDHF code Sky3D version 1.2

The Sky3D code has been widely used to describe nuclear ground states, collective vibrational excitations, and heavy-ion collisions. The approach is based on Skyrme forces or related energy density functionals. The static and dynamic equations are solved on a three-dimensional grid, and pairing is been implemented in the BCS approximation. This updated version of the code aims to facilitate the calculation of nuclear strength functions in the regime of linear response theory, while retaining all existing functionality and use cases. The strength functions are benchmarked against available RPA codes, and the user has the freedom of choice when selecting the nature of external excitation (from monopole to hexadecapole and more). Some utility programs are also provided that calculate the strength function from the time-dependent output of the dynamic calculations of the Sky3D code.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Computing the QRPA level density with the finite amplitude method

Here, we describe a new algorithm to calculate the vibrational nuclear level density of an atomic nucleus. Fictitious perturbation operators that probe the response of the system are generated by drawing their matrix elements from some probability distribution function. We use the Finite Amplitude Method to explicitly compute the response for each such sample. With the help of the Kernel Polynomial Method, we build an estimator of the vibrational level density and provide the upper bound of the relative error in the limit of infinitely many random samples. The new algorithm can give accurate estimates of the vibrational level density. Since it is based on drawing multiple samples of perturbation operators, its computational implementation is naturally parallel and scales like the number of available processing units.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Machine learning-assisted identification of potential sources of bias in measurements of prompt-fission neutron spectra

Unrecognized sources of uncertainty (USU) can bias the reported mean and/or covariance of experimental nuclear data. These biases, in turn, can propagate through evaluated nuclear data to application simulations or may poorly inform nuclear theory that is fitted to the experimental data. Such unknown sources of bias must be tied to the inherent physical constituents of the measurements such as the characteristics of a detector response or a background reduction technique. Here, in this article, a sparse Bayesian learning model is used to support experts in their efforts to identify and characterize USU in experimental prompt fission neutron spectra (PFNS) for spontaneous fissioning of 252 Cf by linking observed biases to features of the measurement system. Three different bias components were found. The first acts as a verification case for the algorithm as it identifies a bias coming from a well-known source related to the use of 6 Li in the neutron detection system. The second two cases demonstrate how this method can benefit the evaluation of experimental nuclear data by identifying, quantifying, and relating unknown biases to potential causes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

High-performance data format for scientific data storage and analysis

Here, in this article, we present the High-Performance Output (HiPO) data format developed at Jefferson Laboratory for storing and analyzing data from Nuclear Physics experiments. The format was designed to efficiently store large amounts of experimental data, utilizing modern fast compression algorithms. The purpose of this development was to provide organized data in the output, facilitating access to relevant information within the large data files. The HiPO data format has features that are suited for storing raw detector data, reconstruction data, and the final physics analysis data efficiently, eliminating the need to do data conversions through the lifecycle of experimental data. The HiPO data format is implemented in C++ and JAVA, and provides bindings to FORTRAN, Python, and Julia, providing users with the choice of data analysis frameworks to use. In this paper, we will present the general design and functionalities of the HiPO library and compare the performance of the library with more established data formats used in data analysis in High Energy and Nuclear Physics (such as ROOT and Parquete). In columnar data analysis, HiPO surpasses established data formats in performance and can be effectively applied to data analysis in other scientific fields.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

HFBTHO-AD: Differentiation of a nuclear energy density functional code

The HFBTHO code implements a nuclear energy density functional solver to model the structure of atomic nuclei. HFBTHO has previously been used to calibrate energy functionals and perform sensitivity analysis by using derivative-free methods. To enable derivative-based optimization and uncertainty quantification approaches, we must compute the derivatives of HFBTHO outputs with respect to the parameters of the energy functional, which are a subset of all input parameters of the code. Here, we use the algorithmic/automatic differentiation (AD) tool Tapenade to differentiate HFBTHO. We compare the derivatives obtained using AD against finite-difference approximation and examine the performance of the derivative computation.

Algorithmic differentiation↗
DOE OSTI2025

A new database website for nuclear level densities

We introduce a new open-access, web-based database (http://nld.ascsn.net), Current Archive of Nuclear Density of Levels (CANDL), that hosts experimental nuclear level density (NLD) datasets from a variety of techniques and energy ranges. Built using the Dash framework in Python, the database is designed to be interactive and user-friendly, allowing researchers to search, visualize, fit, and export NLD data with minimal effort. This resource includes data extracted from evaporation spectra, Oslo method variants, and other experimental techniques that cover excitation energies beyond the neutron resonance region. The database supports on-the-fly fitting with two widely-used phenomenological models—the Constant Temperature (CT) model and the Back-Shifted Fermi Gas (BSFG) model—selected for their simplicity and computational efficiency. Future versions aim to include additional datasets and model types, as well as easy-to-use interfaces to data science techniques. Here, this platform offers a vital tool for the nuclear physics, astrophysics, medicine, and reactor design communities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Data on Cu- and Ni-Si-Mn-rich solute clustering in a neutron irradiated austenitic stainless steel

The data presented in this article is supplementary to the research article “Phase instabilities in austenitic steels during particle bombardment at high and low dose rates” (Levine et al.). Needle-shaped samples were prepared with focused ion beam milling from a 304L stainless steel that was irradiated with fast neutrons (E 0.1 MeV) in the BOR-60 reactor at 318 °C to 47.5 dpa. Atom probe tomography (APT) experiments in voltage mode were then conducted on a Cameca LEAP 5000X HR. Atom position, range, and mass spectrum files after reconstruction with Cameca’s IVAS software are included. Cu- and Ni-Si-Mn-rich solute nanoclusters were identified and analyzed using the Open Source Characterization of APT Reconstructions (OSCAR) program. Python code for OSCAR, information on the program’s underlying algorithm, and sample output files are provided. A proximity histogram of a Ni-Si-Mn-rich cluster and a 1D density/solute concentration profile of a Cu-rich cluster are given to demonstrate OSCAR’s analytical functionalities. The provided APT dataset is valuable for benchmarking phase instabilities in neutron-irradiated austenitic stainless steels that occur at high doses. The OSCAR program can be reused to process other APT data sets where solute nanoclustering is of interest.

42 ENGINEERING↗
DOE OSTI2020

Preliminary results in using Deep Learning to emulate BLOB, a nuclear interaction model

Purpose: A reliable model to simulate nuclear interactions is fundamental for Ion-therapy. We already showed how BLOB (“Boltzmann-Langevin One Body”), a model developed to simulate heavy ion interactions up to few hundreds of MeV/u, could simulate also 12 C reactions in the same energy domain. However, its computation time is too long for any medical application. For this reason we present the possibility of emulating it with a Deep Learning algorithm. Methods: The BLOB final state is a Probability Density Function (PDF) of finding a nucleon in a position of the phase space. We discretised this PDF and trained a Variational Auto-Encoder (VAE) to reproduce such a discrete PDF. As a proof of concept, we developed and trained a VAE to emulate BLOB in simulating the interactions of 12 C with 12 C at 62 MeV/u. To have more control on the generation, we forced the VAE latent space to be organised with respect to the impact parameter (b) training a classifier of b jointly with the VAE. Results: In this work, the distributions obtained from the VAE are similar to the input ones and the computation time needed to use the VAE as a generator is negligible. Conclusions: We show that it is possible to use a Deep Learning approach to emulate a model developed to simulate nuclear reactions in the energy range of interest for Ion-therapy. We foresee the implementation of the generation part in C++ and to interface it with the most used Monte Carlo toolkit: Geant4.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Origin of water masses in Floridan Aquifer System revealed by 81 Kr

Groundwater in coastal aquifers serves as an essential resource in densely populated areas throughout the world. However, hydrological connection to the ocean leaves it at risk of salinization via changes in climate and hydrological cycles. Therefore, an accurate hydrological characterization of coastal aquifers is of the utmost importance. Although localized salinization of shallow aquifers has commonly been studied, geochemical constraints on the time scale of regional freshwater and seawater in coastal aquifers are limited especially in carbonate aquifers. Such information is not only crucial for water resource management, but also for understanding solute flux over the land-ocean boundary and for interpreting the geochemical signatures that each water mass carries. Toward this goal, we report tracer-based, subsurface residence times of groundwater in the Floridan Aquifer System in southern Florida, which is one of the most productive aquifers on Earth. The first application of Kr 81-an emerging tracer of old groundwater-to this aquifer has identified freshwater recharge during the last glacial period as anticipated from previous studies using other geochemical tools including radiocarbon, stable isotopes and noble gas concentrations. Moreover, freshwater of Holocene age was also detected at 100 km-distance from the recharge area, suggesting the possibility of an active flow system in the upgradient region. A contribution of fossil seawater that predates the last glacial maximum was locally identified, suggesting slow seawater circulation and limited but significant solute flux from the ocean into the aquifer for dolomitization. This also implies that the deep saline coastal groundwater can potentially serve as a paleo-seawater archive.

58 GEOSCIENCES↗
DOE OSTI2020

Nuclear response of additive manufactured GRCop-84 copper for use in Lower hybrid launchers in a fusion environment

Recent advances in selective laser melting (SLM) 3D printing technology allow additive manufacture of lower hybrid current drive (LHCD) RF launchers from a new material, Glenn Research Copper 84 (GRCop-84), a Cr 2 Nb (8 at. % Cr, 4 at. % Nb) precipitation hardened alloy, in configurations unachievable with conventional machining. Tensile strength and thermal conductivity are compared to other competing high conductivity copper alloys for fusion use. Swelling and conductivity at 100 DPA when exposed to the Hanford Fast Flux Test Facility neutron spectrum is predicted in the 1–2 % and 65-70% IACS range, respectively, by comparison to copper alloys with similar precipitate and grain size. Finally, self-similar ion irradiation of a GRCop-84 target at 430 °C to 20 DPA resulted in no noticeable void formation in FIB cross sections or TEM lamella.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

3D modelling of MHD mixed convection flow in a vertical duct with transverse magnetic field and volumetric or surface heating

Understanding phenomena associated with the multiple effects/interactions of the fusion nuclear environment on liquid metal flow is required to correctly design liquid metal (LM) blankets for fusion facilities. These effects are investigated in the present work by numerically simulating 3D LM MHD flow. The simulated geometry consists of a straight, vertical duct which runs perpendicular to a strong, fringing applied magnetic field. There is also a region of applied heating as the primary goal is to explore buoyancy effects in MHD duct flows. Results are presented for both buoyancy assisted (upwards) and buoyancy opposed (downwards) flows in conducting and insulating ducts for a range of Hartmann numbers (Ha) up to 100, Reynolds numbers (Re) from 10 3 to 10 4 and Grashof (Gr) numbers from 10 7 to 10 8 . While increasing Gr or decreasing Re increases buoyancy effects, increasing Ha was shown to increase maximum temperature through turbulence reduction. The extent to which the MHD mixed convection flows are quasi-2D is analyzed and buoyant effects, in competition with electromagnetic forces, are shown to bring about 3D flow features not seen in purely MHD flows. Volumetric nuclear heating with steep gradients is applied to the vertical MHD flows for comparison to flows with surface heating only. Finally, surface heating generates stronger buoyancy effects than volumetric heating of the same total power; however, many of the same phenomena occur. Therefore, surface heating, the only option for lab experiments, can provide indication of the effects of volumetric heating in MHD flows.

42 ENGINEERING↗
DOE OSTI2021

Introduction to MDSplus using Docker

With increased use of MDSplus comes an influx of new users. With them comes a need for more and better ways to learn the suite of tools that is MDSplus. Therefore, the MDSplus team continually evaluates new technologies to improve our software and user experience. To this end, we investigated Docker to determine if and how it could help new users understand MDSplus, make MDSplus easier to install, and allow us to easily test old/new versions. To achieve this, a set of Docker images and instructions have been developed. This paper will provide an overview of MDSplus, and detail the methods to create and use the Docker images. Additionally, we will explore the limitations of such an approach, and the recommended applications. The project where these Docker Images were built, along with the Demo is here: https://github.com/WhoBrokeTheBuild/DockerizedMDSplus https://hub.docker.com/r/whobrokethebuild/mdsplus The now official Docker Images are available here: https://github.com/MDSplus/Docker https://hub.docker.com/r/mdsplus/mdsplus

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗