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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Enhancing Monte Carlo Workflows for Nuclear Reactor Analysis with Metamodel-Driven Modeling

Monte Carlo codes are essential components of many reactor physics simulation workflows as high-fidelity continuous-energy neutron transport solvers. Among Monte Carlo radiation transport codes, MCNP is particularly notable due to its diverse simulation capabilities, large user base, and long validation history. Despite being a powerful simulation tool, MCNP provides limited capabilities to allow automated execution, model transformation, or support for user-defined logic and abstractions that limit its compatibility with modern workflows. Here, to better integrate MCNP into a modern scientific workflow, we have developed an intuitive yet full-featured MCNP Application Program Interface (API) in Python, named MCNPy, which provides a specialized set of classes for MCNP input development. Moreover, to guarantee that our reading, writing, and modeling capabilities remain self-consistent (and to render the huge scope of the MCNP API manageable), we have adopted a strategy of model-driven software development in which a generalized model of the MCNP input format has been created. From this generalized model, or “metamodel,” problem-specific implementations such as an engine for input validation or a codebase for programmatic operations may be automatically generated. Since MCNPy primarily acts as a Python front-end to the underlying Java API that directly interfaces with the metamodel, it is intrinsically linked to the metamodel and thus remains maintainable. With MCNPy, users can programmatically read, write, and modify any syntactically valid MCNP input file regardless of its origin. These capabilities allow users to automate complicated tasks like design optimization and model translation for nuclear systems. As examples, this work demonstrates the use of MCNPy to find the critical radius of a plutonium sphere and to translate a 9000+ line MCNP input file into a corresponding OpenMC model.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Modification of Quark-Gluon Distributions in Nuclei by Correlated Nucleon Pairs

We extend the QCD Parton Model analysis using a factorized nuclear structure model incorporating individual nucleons and pairs of correlated nucleons. Our analysis of high-energy data from lepton deep-inelastic scattering, Drell-Yan, and W and Z boson production simultaneously extracts the universal effective distribution of quarks and gluons inside correlated nucleon pairs, and their nucleus-specific fractions. Such successful extraction of these universal distributions marks a significant advance in our understanding of nuclear structure properties connecting nucleon- and parton-level quantities. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Ruthenium solubility and its impact on the crystallization behavior and electrical conductivity of MoO 3 -containing borosilicate-based model high-level nuclear waste glasses

The present study focuses on investigating the solubility of RuO 2 in a borosilicate-based model high-level waste glass and understanding its impact on the crystallization behavior and electrical conductivity of the resulting vitrified waste forms. The solubility limit of RuO 2 in the investigated glass composition has been determined to be 460 ppm by weight. The higher concentration of RuO 2 results in sub-micron sized Ru-rich inclusions in the glassy matrix, which eventually agglomerate to form needle-like and polyhedral RuO 2 crystals. It is observed that RuO 2 selectively promotes the crystallization of the rare-earth apatite phase over the powellite phase. The as-synthesized RuO 2 -containing glasses exhibit semiconducting behavior with a similar level of electrical conductivity below the glass transition. Here, the percolation of non-uniformly distributed RuO 2 inclusions may result in a formation of short-range low-resistive conduction pathways in the host glass matrix leading to an apparent metallic-like behavior of selected thin samples with the highest ruthenium content.

36 MATERIALS SCIENCE↗
DOE OSTI2023

CoSMo - Shell Model Code

CoSMo is a comprehensive shell model code suite designed for nuclear shell model and configuration interaction calculations. Its highly structured and templated nature allows for flexibility and ease in applications, including those to open quantum systems with non-Hermitian Hamiltonians, clustering, and time-dependent dynamics.

Volya, Alexander↗
DOE OSTI2025

SPARC - Plans for a New Critical Experiment Facility with a Horizontal Split Table

Several critical experiment facilities, sometimes referred to as zero power reactor facilities, have provided crucial data to aid understanding and validate nuclear-physics models since the beginning of nuclear technology. Indeed, the first man-made reactor, Chicago Pile-1, was essentially this type of reactor. However, there was a downturn in nuclear technology development toward the turn of the millennium, and the need for these specialized research facilities waned. Now there are few of these experimental facilities operational in the world and those that remain have relatively small critical assembly machines. The need for criticality safety benchmark experiments at intermediate neutron energy levels and the modern resurgence of interest in advanced reactors designs, many of which do not have historical precedents in terms of nuclear fuel composition, moderator, and coolant combinations, all combine to create a substantial need for a critical experiment facility with a large horizontal split-table (HST) machine. A HST machine is used to arrange two separate and subcritical parts of a core assembly, bring them together in a precise manner to achieve criticality using remote controls, and separate them to achieve a subcritical configuration again. A new effort was recently performed to develop user needs for a HST, assess candidate locations at the Idaho National Laboratory (INL), and develop a plan for deployment. This project is referred to as the System Physics Advanced Reactor Critical facility (SPARC). A few months after this assessment began, and shortly after as a viable pathway was emerging, a series of important presidential executive orders were issued to revitalize nuclear energy in the United States (U.S.). The relevance of SPARC to these executive orders was immediately apparent. The far-reaching potential of SPARC to these executive orders will reside in its ability to produce data which facilitates licensing of advanced nuclear reactor designs while reducing uncertainties to help increase energy production alongside new criticality safety data to enable more efficient nuclear fuel manufacture, transport, and storage.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Gas jet targets for direct reaction studies

The study of direct reactions is of broad interest in nuclear physics, providing constraint to models of nuclear structure evolution and data to better understand the creation of the elements. In many cases, however, the data of interest are hindered by backgrounds and poor resolution from contaminants in either the beam, the target, or both. The use of a gas jet can overcome some of these issues through clever engineering, providing a reaction target that is chemically pure and thin enough to significantly reduce the impact on experimental resolution. This Perspective will discuss the effort to design, construct, and operate gas jet targets for direct reaction studies in the rare isotope era.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2003

Development of A Monte Carlo Radiation Transport Code System For HEDS: Status Update

Modifications of the Monte Carlo radiation transport code HETC are underway to extend the code to include transport of energetic heavy ions, such as are found in the galactic cosmic ray spectrum in space. The new HETC code will be available for use in radiation shielding applications associated with missions, such as the proposed manned mission to Mars. In this work the current status of code modification is described. Methods used to develop the required nuclear reaction models, including total, elastic and nuclear breakup processes, and their associated databases are also presented. Finally, plans for future work on the extended HETC code system and for its validation are described.

Townsend, Lawrence W.↗
DOE OSTI2024

Comparing modelling approaches for a generic nuclear waste repository in salt

This paper contains a comparison of five modelling approaches for a simplified nuclear waste repository in a domal salt formation. It is the result of a four-year collaboration between five international teams on Task F of the DECOVALEX-2023 project on performance assessment modelling. The primary objectives of Task F are to build confidence in the models, methods, and software used for performance assessment (PA) of deep geologic nuclear waste repositories, and/or to bring to the fore additional research and development needed to improve PA methodologies. This work demonstrates how these objectives are accomplished through staged development and comparison of the models and methods used by participating teams in their PA frameworks. Participating teams made a wide range of model assumptions, ranging from compartmentalized networks to full 3D models of the salt formation and repository. Despite differences in the modelling strategies, all models indicate that salt compaction and diffusion of radionuclides in brine are key processes in the repository. For the isothermal spent nuclear fuel and vitrified waste scenario with multiple early failures considered, all models indicate little of the disposed radionuclides will migrate beyond the repository seal over the 100,000-year simulations. In general, the model output quantities have the largest differences over the short term and near the waste. Disparities between the models are believed to be due to differing simplifications from the conceptual model.

DECOVALEX↗
DOE OSTI2022

Simultaneous Optimization of Nuclear–Electronic Orbitals

Accurate modeling of important nuclear quantum effects, such as nuclear delocalization, zero-point energy, and tunneling, as well as non-Born-Oppenheimer effects, requires treatment of both nuclei and electrons quantum mechanically. The nuclear–electronic orbital (NEO) method provides an elegant framework to treat specified nuclei, typically protons, on the same level as the electrons. In conventional electronic structure theory, finding a converged ground state can be a computationally demanding task; converging NEO wavefunctions, due to their coupled electronic and nuclear nature, is even more demanding. Herein, we present an efficient simultaneous optimization method that uses the direct inversion in the iterative subspace method to simultaneously converge wavefunctions for both the electrons and quantum nuclei. In conclusion, benchmark studies show that the simultaneous optimization method can significantly reduce the computational cost compared to the conventional stepwise method for optimizing NEO wavefunctions for multicomponent systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

System-size dependence of charged-particle suppression in ultrarelativistic nucleus-nucleus collisions

High-energy partons lose energy while propagating through the hot, strongly interacting medium produced in ultrarelativistic nucleus-nucleus collisions, leading to a suppression of particle production at high transverse momentum (p T ). The dependence of this energy loss on the size of the colliding nuclear system has yet to be firmly established experimentally. This Letter presents a systematic study of charged-particle suppression across four different nucleus-nucleus collision systems using nuclear modification factors (R AA ) measured by the CMS Collaboration at the CERN LHC. Previous CMS measurements of R AA in oxygen-oxygen, xenon-xenon, and lead-lead collisions are recast with identical p T intervals and are complemented by the first measurement of the charged-particle R AA in neon-neon collisions at $\sqrt{S_{NN}}$ = 5.36 Te V . The neon-neon data correspond to an integrated luminosity of 0.76 nb −1 . The R AA in all collision systems examined show similar qualitative trends as a function of p T , but have a magnitude which is ordered with the nucleon number A. The R AA feature a downward slope at low p T , a local minimum at around 5–7 Ge V , and an upward slope with increasing p T . The R AA are also compared in terms of A 1/3 , which is proportional to the nuclear radius. Models including only initial-state nuclear effects fail to reproduce the observed trends, whereas energy loss models reproduce the trends in the region p T > 9.6 Ge V .

CMS↗
DOE OSTI2024

Strong Coupling of Hydrodynamics and Reactions in Nuclear Statistical Equilibrium for Modeling Convection in Massive Stars

We build on the simplified spectral deferred corrections (SDC) coupling of hydrodynamics and reactions to handle the case of nuclear statistical equilibrium (NSE) and electron/positron captures/decays in the cores of massive stars. Our approach blends a traditional reaction network on the grid with a tabulated NSE state from a very large, ${\mathcal O }(100)$ nuclei network. We demonstrate how to achieve second-order accuracy in the simplified-SDC framework when coupling NSE to hydrodynamics, with the ability to evolve the star on the hydrodynamics time step. We discuss the application of this method to convection in massive stars leading up to core collapse. We also show how to initialize the initial convective state from a 1D model in a self-consistent fashion. All of these developments are done in the publicly available Castro simulation code and the entire simulation methodology is fully GPU-accelerated.

Explosive nucleosynthesis↗
DOE OSTI2025

Air Force Institute of Technology Dust Cloud Model Vertical Dispersion

The Air Force Institute of Technology (AFIT) has developed and refined a nuclear dust cloud model for modeling nuclear fallout, and the dust mass and radiation encountered by aircraft flying through a nuclear dust cloud. Recently, a question was raised requesting an explanation of why the vertical mass and activity distributions had an anomalous shape for yields above 1 Mt. This report explains that the anomalous shape occurs because the dust cloud model is a discrete numeric model, and not a continuous model.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Decoding Zeolite Crystallization and Stage III in Nuclear Waste Glasses by Coupled Modeling and Experiments

Under specific conditions of pH and temperature, nuclear waste immobilization borosilicate glasses may exhibit a sudden acceleration in their corrosion kinetics (stage III)—a behavior that has been associated with the formation of zeolite crystals. Such accelerated dissolution may compromise the integrity of nuclear wasteforms placed in geological depositories. However, thus far, none of the available models is able to predict the thermodynamic propensity and kinetics of zeolite precipitation as a function of the solution conditions due to (i) a lack of fundamental knowledge regarding the nucleation & growth mechanisms of zeolitic phases, (ii) uncertainty regarding the compositions (types) of zeolites that may form and the rate-limiting step in their precipitation as a function of the solution conditions, and (iii) the complexities that arise due to the vast parametric space (i.e., solution chemistry, temperature, number of secondary phases, etc.) that encompass these systems under conditions of environmental exposure. To resolve these challenges, this project aimed to unambiguously identify the thermodynamic propensity for zeolite precipitation and the kinetics thereof as a function of the solution conditions (composition, pH, and temperature). To achieve this goal: 1) We identified the solution conditions and zeolite phases relevant to nuclear glass dissolution. 2) We performed a series of ab initio molecular dynamics (AIMD) simulations to compute the thermodynamic properties of a group of characteristic zeolites that features a large range of compositions, various hydration levels, a wide range of framework structures, and partial atomic site occupancies. 3) We released a first-of-a-kind self-consistent thermodynamic database that can be used to assess the kinetics and the stability fields of zeolitic phases within a Gibbs energy minimization (GEM) framework. 4) We developed a robust geochemical modeling method allowing us to predict the stability of secondary phases (including zeolites, calcium–silicate–hydrate gels, and clays) upon the dissolution of nuclear waste immobilization glasses. 5) We introduced a model that predicts the dissolution kinetics of a series of borosilicate nuclear waste immobilization glasses in terms of the topology of their atomic network. 6) We investigated the roles of the solution composition on the crystallization kinetics of phillipsite zeolites and tobermorite silicate hydrates. Via PNNL’s collaboration and engagement, this project directly supports DOE’s nuclear waste immobilization activities by offering a technical, science-based foundation that will (i) facilitate predictions of the long-term corrosion rates and extents of existing nuclear waste immobilization glasses to help ensure safe and successful vitrification operations, (ii) inform the development of advanced glass formulations with enhanced durability, and, (iii) enable cost-savings that result from making more decisive and hence less conservative predictions while offering higher levels of nuclear waste embedment in smaller, more compact glass volumes.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Modeling Kilonova Light Curves: Dependence on Nuclear Inputs

The mergers of binary neutron stars, as well as black hole–neutron star systems, are expected to produce an electromagnetic counterpart that can be analyzed to infer the element synthesis that occurred in these events. We investigate one source of uncertainties pertinent to lanthanide-rich outflows: the nuclear inputs to rapid neutron capture nucleosynthesis calculations. We begin by examining 32 different combinations of nuclear inputs: eight mass models, two types of spontaneous fission rates, and two types of fission daughter product distributions. We find that such nuclear physics uncertainties typically generate at least one order of magnitude uncertainty in key quantities such as the nuclear heating (one and a half orders of magnitude at 1 day post-merger), the bolometric luminosity (one order of magnitude at 5 days post-merger), and the inferred mass of material from the bolometric luminosity (factor of 8 when considering the 8–10 day region). Since particular nuclear processes are critical for determining the electromagnetic signal, we provide tables of key nuclei undergoing β-decay, α-decay, and spontaneous fission important for heating at different times, identifying decays that are common among the many nuclear input combinations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Search for a Non-Zero Value of the Electric Dipole Moment of the Neutron (Final Technical Report)

The research grant addressed in this final report began its most recent 3-year segment June 1, 2022 and terminated May 31, 2025. The work funded by the grant supported our group at the Massachusetts Institute of Technology (MIT) in its collaboration on the nEDM@SNS (SNS refers to the Spallation Neutron Source) project at the Oak Ridge National Laboratory (ORNL). This project had as its goal the measurement of a non-zero value of the electric dipole moment of the neutron (nEDM). If indeed the neutron has a non-zero electric dipole moment, it will tell us very important information about the Standard Model of Nuclear and Particle Physics. The Standard Model is remarkably successful in describing much of the world we live in, but we know it is not complete. A non-zero nEDM would indicate the existence of strong CP violation, which would be a major addition to the Standard Model.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Multilevel Analysis, Design, and Modeling of Coupling Advanced Nuclear Reactors and Thermal Energy Storage in an Integrated Energy System

This report discusses the different options for coupling thermal energy storage (TES) systems to advanced nuclear power plants (A-NPPs) in order to enable flexible and hybrid plant operation. An advanced light-water reactor (A LWR), a high-temperature gas-cooled reactor (HTGR) and a liquid-metal fast reactor (LMFR) were selected as the initial use cases for demonstrating a thermally balanced energy storage coupling design for thermal power extraction. The models presented herein showcase several design considerations, focusing on optimal deployment methodologies for achieving steady-state and transient-state operation with minimum disruption to the nuclear power cycle. This first part of the study presents steady-state models developed using Aspen HYSYS®, with the thermal energy bypass for NPP-TES coupling being varied at up to 50%. The various components were sized using the Aspen Process Economic Analyzer (APEA) and Aspen Exchanger Design and Rating (EDR), when applicable. Cost functions from these models were developed using the latest publicly available data obtained from APEA V11. The TES-coupled A-NPP steady-state models and cost functions then provided a baseline for dynamic operation and process optimization by using Idaho National Laboratory (INL)’s Framework for Optimization of Resources and Economics (FORCE) tools. A stochastic optimization of the various energy storage systems coupled to the A-NPPs was then performed using the Risk Analysis Virtual Environment (RAVEN) and its dispatch optimization plugin, the Holistic Energy Resource Optimization Network (HERON). The signal processing and synthetic history capabilities of RAVEN were used to account for the unpredictable behavior of electricity markets. An autoregressive moving average (ARMA) model was used to analyze price signals from the Pennsylvania-New Jersey-Maryland (PJM) market and were applied to the HERON analysis in order to optimize a system with the best economics. Transient modeling evaluation was then performed using Modelica models within the HYBRID repository, which was developed at INL for the Department of Energy Integrated Energy Systems program for the characterization of dynamic integrated system behavior and feedback. This includes evaluation of the TES-coupled A-LWR systems’ impact on physical and thermal system response during imposed system demands. Additional TES-coupled reactor types, coupling approaches, markets, and TES technologies will be evaluated in future work.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Impact of the MARATHON data on F 2⁢n⁡ /F 2⁢p and off-shell effects in light nuclei

The neutron structure function, F 2⁢n⁡ , has historically been extracted from measurements of the deuteron structure function, which is limited by our understanding of the nuclear effects on the bound proton and neutron. The MARATHON Collaboration recently extracted F 2⁢n⁡ from the comparison of 3 H and 3 He targets, where the nuclear effects are larger but nearly identical, yielding a precise extraction of F 2⁢n⁡ /F 2⁢p . Here, to ensure that this comparison is not biased by the specific model of nuclear effects used by MARATHON, we examine a range of models of the nuclear effects to obtain a more conservative, but more model-independent, extraction of F 2⁢n⁡ /F 2⁢p for comparison with deuteron extractions. Even with this approach, the results are precise enough to compare to deuteron extractions, providing new constraints on the nuclear effects in the deuteron and indicating the need for larger nuclear corrections than those included in most models.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗