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 505 records · Page 28

Bayesian estimation of the $S$ factor and thermonuclear reaction rate for 16 O(p, γ) 17 F

The 16 O(p, γ) 17 F reaction is the slowest hydrogen-burning process in the CNO mass region. Its thermonuclear rate sensitively impacts predictions of oxygen isotopic ratios in a number of astrophysical sites, including AGB stars. The reaction has been measured several times at low bombarding energies using a variety of techniques. The most recent evaluated experimental rates have a reported uncertainty of about 7.5% below 1 GK. However, the previous rate estimate represents a best guess only and was not based on rigorous statistical methods. We apply a Bayesian model to fit all reliable 16 O(p, γ) 17 F cross section data, and take into account independent contributions of statistical and systematic uncertainties. The nuclear reaction model employed is a single-particle potential model involving a Woods-Saxon potential for generating the radial bound state wave function. The model has three physical parameters, the radius and diffuseness of the Woods-Saxon potential, and the asymptotic normalization coefficients (ANCs) of the final bound state in 17 F. Here, we find that performing the Bayesian S -factor fit using ANCs as scaling parameters has a distinct advantage over adopting spectroscopic factors instead. Based on these results, we present the first statistically rigorous estimation of experimental 16 O(p, γ) 17 F reaction rates, with uncertainties (±4.2%) of about half the previously reported values.

6 ≤ A ≤ 19↗

SSA Monte Carlo and Master Equation Modeling of Neutron Leakage Distributions

Strongly stochastic neutron populations in highly subcritical multiplying media containing weak spontaneous fission sources are unique indicators of the presence of SNM. For instance, in neutron multiplicity counting, an extended tail in the count number probability distribution is observed that is a characteristic of and sensitive to the amount of SNM present. Theoretical modeling and efficient computation of stochastic neutron distributions are complementary to experimental measurements and moreover shed light on the fundamental interactions and transport processes that underlie nuclear signatures. Modeling approaches that are intermediate between gross-system level models and high-fidelity but computationally demanding Monte Carlo simulations are of particular interest in this regard. Furthermore we present three advances on existing work in point kinetic or lumped geometry: (i) implementation of the Stochastic Simulation Algorithm (SSA), a system state-updating Monte Carlo method that is more e cient than single-event Monte Carlo in lumped geometries, (ii) application of the backward Master equation to obtain time-gated count moments, and (iii) reconstructing the count number distribution from the deterministically computed moments using a generalized Laguerre polynomial expansion.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A parameter sensitivity study for simulating DNA damage after proton irradiation using TOPAS-nBio

Monte Carlo (MC) track structure simulation tools are commonly used for predicting radiation induced DNA damage by modeling the physical and chemical reactions at the nanometer scale. However, the outcome of these MC simulations is particularly sensitive to the adopted parameters which vary significantly across studies. In this study, a previously developed full model of nuclear DNA was used to describe the DNA geometry. The TOPAS-nBio MC toolkit was used to investigate the impact of physics and chemistry models as well as three key parameters (the energy threshold for direct damage, the chemical stage time length, and the probability of damage between hydroxyl radical reactions with DNA) on the induction of DNA damage. Our results show that the difference in physics and chemistry models alone can cause differences up to 34% and 16% in the DNA double strand break (DSB) yield, respectively. Additionally, changing the direct damage threshold, chemical stage length, and hydroxyl damage probability can cause differences of up to 28%, 51%, and 71% in predicted DSB yields, respectively, for the configurations in this study.

60 APPLIED LIFE SCIENCES↗

Assessment and Requirements of Nuclear Reaction Databases for GCR Transport in the Atmosphere and Structures

The transport properties of galactic cosmic rays (GCR) in the atmosphere, material structures, and human body (self-shielding) am of interest in risk assessment for supersonic and subsonic aircraft and for space travel in low-Earth orbit and on interplanetary missions. Nuclear reactions, such as knockout and fragmentation, present large modifications of particle type and energies of the galactic cosmic rays in penetrating materials. We make an assessment of the current nuclear reaction models and improvements in these model for developing required transport code data bases. A new fragmentation data base (QMSFRG) based on microscopic models is compared to the NUCFRG2 model and implications for shield assessment made using the HZETRN radiation transport code. For deep penetration problems, the build-up of light particles, such as nucleons, light clusters and mesons from nuclear reactions in conjunction with the absorption of the heavy ions, leads to the dominance of the charge Z = 0, 1, and 2 hadrons in the exposures at large penetration depths. Light particles are produced through nuclear or cluster knockout and in evaporation events with characteristically distinct spectra which play unique roles in the build-up of secondary radiation's in shielding. We describe models of light particle production in nucleon and heavy ion induced reactions and make an assessment of the importance of light particle multiplicity and spectral parameters in these exposures.

Cucinotta, F. A.↗

Coherent photoproduction of ρ0 vector mesons in ultra-peripheral Pb-Pb collisions at $ \sqrt{{\mathrm{s}}_{\mathrm{NN}}} $ = 5.02 TeV

Cross sections for the coherent photoproduction of ρ0 vector mesons in ultra-peripheral Pb-Pb collisions at $ \sqrt{{\mathrm{s}}_{\mathrm{NN}}} $ = 5.02 TeV are reported. The measurements, which rely on the π+π− decay channel, are presented in three regions of rapidity covering the range |y| < 0.8. For each rapidity interval, cross sections are shown for different nuclear-breakup classes defined according to the presence of neutrons measured in the zero-degree calorimeters. The results are compared with predictions based on different models of nuclear shadowing. Finally, the observation of a coherently produced resonance-like structure with a mass around 1.7 GeV/c2 and a width of about 140 MeV/c2 is reported and compared with similar observations from other experiments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Loosely bound composite dark matter

We investigate loosely bound composite states made of dark matter, where the binding energy for constituent particles is considerably less than the constituent mass. We focus on models of nuclear and molecular dark matter, where constituents are separated by length scales larger than the inverse constituent mass, just like nuclei and atoms in the Standard Model. The cosmology, structure, and interactions at underground experiments are described. We find that loosely bound composites can have a very large cross section for scattering with nuclei that scales with nucleon number like ∼ A 4 . For some couplings, these composites produce extremely soft (≪ keV) individual atomic recoils while depositing a large amount of total recoil energy (≫ keV) in a single passage through a detector, implying an interesting new class of signatures for low threshold direct detection.

79 ASTRONOMY AND ASTROPHYSICS↗

Measurement of the azimuthal anisotropy of charged particles in sNN=5.36TeV O16+O16 and Ne20+Ne20 collisions with the ATLAS detector

This paper presents the first measurements of the azimuthal anisotropy coefficients , which quantify the -order Fourier modulation of charged-particle azimuthal distributions, for in and collisions recorded with the ATLAS detector at the CERN Large Hadron Collider in 2025. The coefficients are measured as a function of transverse momentum ( ), collision centrality, and event multiplicity. They are extracted using two complementary methods: two-particle correlations with a template-fit subtraction of short-range nonflow contributions, and four-particle subevent cumulants, which intrinsically suppress nonflow effects and provide sensitivity to flow fluctuations. The results show a clear hierarchy and a nonmonotonic dependence on , reaching a maximum around , consistent with trends observed in heavy-ion collisions. Detailed comparisons between the two collision systems reveal an enhanced in central collisions, consistent with theory expectations based on the predicted prolate deformation of neon nuclei, in contrast to the slightly tetrahedral structure predicted for oxygen. The four-particle cumulant results highlight strong event-by-event fluctuations and provide the greatest sensitivity to nuclear shape effects. These measurements can place new constraints on the initial geometry and the hydrodynamic response in light-ion collisions, offering valuable input for models of nuclear structure.

Aad, G↗

Measurement of the azimuthal anisotropy of charged particles in $\sqrt{s_{\mathrm{NN}}}=5.36$ TeV $^{16}$O$+^{16}$O and $^{20}$Ne$+^{20}$Ne collisions with the ATLAS detector

This paper presents the first measurements of the azimuthal anisotropy coefficients $v_{n}$, which quantify the $n^{\mathrm{th}}$-order Fourier modulation of charged-particle azimuthal distributions, for $n=2$-4 in $\sqrt{s_{\mathrm{NN}}}=5.36$ TeV $\mathrm{^{16}O}+\mathrm{^{16}O}$ and $\mathrm{^{20}Ne}+\mathrm{^{20}Ne}$ collisions recorded with the ATLAS detector at the Large Hadron Collider in 2025. The $v_{n}$ coefficients are measured as a function of transverse momentum ($p_{\mathrm{T}}$), collision centrality, and event multiplicity. They are extracted using two complementary methods: two-particle correlations with a template-fit subtraction of short-range non-flow contributions, and four-particle subevent cumulants, which intrinsically suppress non-flow effects and provide sensitivity to flow fluctuations. The results show a clear hierarchy $v_{2} > v_{3} > v_{4}$ and a non-monotonic dependence on $p_{\mathrm{T}}$, reaching a maximum around 2 GeV, consistent with trends observed in heavy-ion collisions. Detailed comparisons between the two collision systems reveal an enhanced $v_{2}$ in central $\mathrm{^{20}Ne}+\mathrm{^{20}Ne}$ collisions, consistent with theory expectations based on the predicted prolate deformation of neon nuclei, in contrast to the slightly tetrahedral structure predicted for oxygen. The four-particle cumulant results highlight strong event-by-event fluctuations and provide the greatest sensitivity to nuclear shape effects. These measurements can place new constraints on the initial geometry and the hydrodynamic response in light-ion collisions, offering valuable input for models of nuclear structure.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Improved Nuclear Reactor and Shield Mass Model for Space Applications

New technologies are being developed to explore the distant reaches of the solar system. Beyond Mars, solar energy is inadequate to power advanced scientific instruments. One technology that can meet the energy requirements is the space nuclear reactor. The nuclear reactor is used as a heat source for which a heat-to-electricity conversion system is needed. Examples of such conversion systems are the Brayton, Rankine, and Stirling cycles. Since launch cost is proportional to the amount of mass to lift, mass is always a concern in designing spacecraft. Estimations of system masses are an important part in determining the feasibility of a design. I worked under Michael Barrett in the Thermal Energy Conversion Branch of the Power & Electric Propulsion Division. An in-house Closed Cycle Engine Program (CCEP) is used for the design and performance analysis of closed-Brayton-cycle energy conversion systems for space applications. This program also calculates the system mass including the heat source. CCEP uses the subroutine RSMASS, which has been updated to RSMASS-D, to estimate the mass of the reactor. RSMASS was developed in 1986 at Sandia National Laboratories to quickly estimate the mass of multi-megawatt nuclear reactors for space applications. In response to an emphasis for lower power reactors, RSMASS-D was developed in 1997 and is based off of the SP-100 liquid metal cooled reactor. The subroutine calculates the mass of reactor components such as the safety systems, instrumentation and control, radiation shield, structure, reflector, and core. The major improvements in RSMASS-D are that it uses higher fidelity calculations, is easier to use, and automatically optimizes the systems mass. RSMASS-D is accurate within 15% of actual data while RSMASS is only accurate within 50%. My goal this summer was to learn FORTRAN 77 programming language and update the CCEP program with the RSMASS-D model.

Robb, Kevin↗

Advancing Multi-Hazard Risk and Safety Considerations for Aging Nuclear Facilities (Revision 1)

This project will demonstrate a multi-hazard time-dependent probabilistic risk assessment (PRA) approach for nuclear facilities considering aging-related deterioration of structures. A generic pressurized water reactor (PWR) reactor subjected to seismic mainshock-aftershock sequences considering the aging of the containment structure will be used as a case study to demonstrate the multi-hazard PRA approach. Using advanced modeling and simulation, seismic mainshock-aftershock fragility functions will be simulated for the containment structure considering aging effects. A multi-hazard PRA model for a generic PWR reactor will be built to quantify the multi-hazard core damage frequency (CDF) and large early release frequency (LERF) with explicit time-dependent modeling of event sequences. To date, the cascading impacts of multi-hazards are not adequately accounted for in the PRA models for nuclear facilities. In addition, for both initial and periodic evaluation of facilities to withstand natural phenomena hazards (NPH), deterioration of the SSCs due to aging and other effects is not adequately considered. By advancing multi-hazard considerations accounting for aging effects, this project will contribute to improved understanding of the safety of aging nuclear facilities. Project outcomes such as the multi-hazard CDF and LERF will also allow facility owners to optimize upgrade and retrofit protocols. This project is divided into two components: (1) advanced modeling and simulations; and (2) multi-hazard PRA. For the first component, Idaho National Laboratory’s (INL) Multi-hazard Analysis for STOchastic time-DOmaiN phenomena (MASTODON) and BlackBear codes will be used to simulate the seismic response and damage of a containment structure under cascading mainshock-aftershock sequences with aging effects. Uncertainties related to the seismic inputs, material parameters, and environmental factors such as temperature and humidity will be identified and propagated to the fragility functions. These fragility functions, which consider aging effects, will be time dependent. The second component will take the fragility information from the first component to build a multi-hazard time-dependent PRA model starting from a generic PWR model to evaluate the multi-hazard CDF and LERF and characterize the associated consequences. For this component, OpenPRA Web Application and SAPHIRE code will be used. Events such as the Fukushima Daiichi accident have highlighted the importance of considering the cascading impacts of multi-hazards for PRA. Moreover, many of the reactors in the current nuclear fleet in the United States (US) are already operating well beyond their initially planned design life, and applications for further extensions to operating licenses are being considered. Therefore, considering multi-hazard effects and aging deterioration in the NPH risk assessment process will contribute to the safety of both existing and future nuclear facilities. Outcomes of this project will thus directly benefit the standards DOE-STD-1020-16 and DOE-HDBK-1224-18.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Selection and Ranking of Experiments from the Halden Database in support of Multiscale Model Validation

Validating fuel performance codes, such as BISON, requires an extensive amount of experiments covering a wide range of operating conditions and fuel types. Within the light-water reactor (LWR) space, there have been several international experimental programs that have contributed to the wealth of available experimental data available for use. One of those international programs, the Halden Reactor Project (HRP), began in 1958 and utilized the Halden Boiling Water Reactor to conduct many highly instrumented experiments until the reactor closed in 2018. Idaho National Laboratory, through the U.S. Department of Energy, has utilized several Halden experiments to perform the initial validation of the BISON code based upon their inclusion in international modeling and simulation benchmarks. Recently, the HRP has provided member organizations a complete copy of all available data, reports, and presentations since the HRP began. This report provides an initial exploration of the data available in the database for use in validating the multiscale models under development in the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program for LWR applications. Ranking tables that identify potential validation cases are provided for the high priority models of interest. It was found that some of the recommended high priority experiments correspond to additional rods in existing assemblies already available in the BISON validation suite. It is expected that several of these cases will be incorporated into future NEAMS milestones in the fuels technical area for increased validation of BISON.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Towards the Development of Molten Salt Reactor Multiphysics Modeling Capabilities: Depletion, Thermochemistry, Reactor Dynamics, and Species Tracking

This report documents the simulation and modeling work performed under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, utilizing the tools and codes developed under the program. The focus in this Fiscal Year (FY) is on the depletion, thermochemistry, reactor dynamics, and species transport in MSRs. The fission products travelling in the entire primary loop is unique in MSRs. There are many physical/chemical mechanisms that can significantly affect the fission products concentration in the fuel salt, e.g., diffusion into gas bubbles, diffusion in the core moderators, deposition on the surface of structures/pipes, etc. The fission product concentration in the fuel salt can potentially affect the neutronics, system dynamics, chemical potentials, etc., so accurately predicting the fission products distribution in the primary loop is important. Under the main objective, different activities are pursued to assess the NEAMS tools applicability for a variety of areas on species transport, including the depletion analyses using Griffin (fission product source), the redox potential assessment at different depletion points, the assessment of xenon transport effect on the neutronics and reactor dynamics, a survey of correlations for calculating the mass transfer coefficient of different fission products, etc.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Industry Level Integrated Fire Modeling Using Fire Risk Investigation in 3D (FRI3D)

The software Fire Risk Investigation in 3D (FRI3D) has been developed over the last 2 years to integrate 3D spatial modeling with existing fire probabilistic risk assessment (PRA) models and fire simulation codes. The goal of this research and development is to automate many of the fire analysis manual tasks to reduce industry efforts in the initial fire modeling and operational costs for the model maintenance and evaluations required during normal plant operations. The tasks for Fiscal Year (FY) 2021 include first testing the FRI3D modeling capabilities by importing an industry fire model into FRI3D and making a 3D model of a complex/high-risk significant area. (For this work, the switchgear room was chosen.) Then, the second task of FY 2021 is to develop a dynamic fire PRA process that can help optimize traditional fire PRA models. The switchgear room model will be used for the dynamic fire PRA work. This report describes the work and insights learned when using FRI3D software to model both a Nuclear Regulatory Report (NUREG) example models and a full industry switchgear room.

97 MATHEMATICS AND COMPUTING↗

Quantifying model prediction sensitivity to model-form uncertainty

Computational and mathematical models are essential to understanding complex systems and phenomena. However, when developing such models, limited knowledge and/or resources necessitates the use of simplifying assumptions. It is therefore crucial to quantify the impact of such simplifying assumptions on the reliability and accuracy of resulting model predictions. This work develops a first-of-its-kind approach to quantify the impact of physics modeling assumptions on predictions. Here, we leverage the emerging field of model-form uncertainty (MFU) representations, which are parameterized modifications to modeling assumptions, in combination with grouped Sobol’ indices to quantitatively measure an assumption’s importance. Specifically, we compute the grouped Sobol’ index for the MFU representation’s parameters as a single importance measure of the assumption for which the MFU representation characterizes uncertainty. To ensure this approach is robust to the subjective choice of how to parameterize a MFU representation, we establish bounds for the difference between sensitivity results for two different MFU representations based on differences in model prediction statistics. The capabilities associated with this approach are demonstrated on three exemplar problems: an upscaled subsurface contaminant transport problem, ablation modeling for hypersonic flight, and nuclear waste repository modeling. We found that our grouped approach is able to assess the impact of modeling assumptions on predictions and offers computational advantages over classical Sobol’ index computation while providing more interpretable results.

97 MATHEMATICS AND COMPUTING↗

Leveraging differentiable programming in the inverse problem of neutron stars

Neutron stars (NSs) probe the high-density regime of the nuclear equation of state (EOS). However, inferring the EOS from observations of NSs is a computationally challenging task. Here, in this work, we efficiently solve this inverse problem by leveraging differential programming in two ways. First, we enable full Bayesian inference in under one hour of wall time on a GPU by using gradient-based samplers, without requiring pretrained machine learning emulators. Moreover, we demonstrate efficient scaling to high-dimensional parameter spaces. Second, we introduce a novel gradient-based optimization scheme that recovers the EOS of a given NS mass-radius curve. We demonstrate how our framework can reveal consistencies or tensions between nuclear physics and astrophysics. First, we show how the breakdown density of a metamodel description of the EOS can be determined from NS observations. Second, we demonstrate how degeneracies in EOS modeling using nuclear empirical parameters can influence the inverse problem during gradient-based optimization. Looking ahead, our approach opens up new theoretical studies of the relation between NS properties and the EOS, while effectively tackling the data analysis challenges brought by future detectors.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Nucleonic shells and nuclear masses

The binding energy of an isotope is a sensitive indicator of the underlying shell structure as it reflects the net energy content of a nucleus. Since magic nuclei are significantly lighter, or more bound, compared to their neighbors, the presence of nucleonic shell structure makes an imprint on nuclear masses. In this study, using a carefully designed binding-energy indicator, we catalog the appearance of spherical and deformed shell and subshell closures throughout the nuclear landscape. After presenting experimental evidence for shell and subshell closures as seen through the lens of nuclear masses, we study the ability of global nuclear mass models to predict local binding-energy variations related to shell effects.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Nuclear Propulsion Technical Interchange Meeting, volume 2

The purpose of the meeting was to review the work performed in fiscal year 1992 in the areas of nuclear thermal and nuclear electric propulsion technology development. These proceedings are an accumulation of the presentations provided at the meeting along with annotations provided by authors. The proceedings cover system concepts, technology development, and system modeling for nuclear thermal propulsion (NTP) and nuclear electric propulsion (NEP). The test facilities required for the development of the nuclear propulsion systems are also discussed.

Source record↗

The Virtual Test Bed (VTB) Repository: A Library of Multiphysics Reference Reactor Models using NEAMS Tools

With the next generation of nuclear reactors under development, modeling and simulation (M&S) tools are being developed by the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program in order to support their design, licensing, and future operation. Mirroring the physical test beds currently under construction (i.e., EBR-II and ZPPR), the Virtual Test Bed (VTB) was launched by the National Reactor Innovation Center (NRIC) in collaboration with NEAMS to support the advanced reactor community. This collaborative effort, which involves multiple teams at both Idaho National Laboratory and Argonne National Laboratory aims to use NEAMS tools to model a wide range of reactor designs. Those models are automatically tested to ensure their continued functionality as the tools are further developed. Examples are extensively documented, each acting as a tutorial for applying the relevant NEAMS tools to that reactor design. Currently, five advanced reactor types (with a total of eight specific design variants) are simulated by a variety of different models. These models range from steady-state, core multiphysics simulations to integrated plant analysis during loss-of-flow transients.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗