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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Extension of Virtual Test Bed Advanced Modeling and Simulation Capabilities for Fusion Energy

The National Reactor Innovation Center (NRIC) was established to accelerate the deployment of novel reactor concepts. This is achieved by providing physical and virtual spaces for building and testing reactor experiments. The Virtual Test Bed (VTB) represents the virtual counterpart to the physical test beds. It is a collaboration with the Department of Energy’s (DOE) Nuclear Energy Advanced Modeling and Simulation (NEAMS) program with the mission to accelerate the deployment of advanced reactors by facilitating the adoption of advanced modeling and simulation (M&S) tools developed by the DOE. This mission has been carried out by the VTB since 2020 by hosting and featuring dozens of advanced fission nuclear reactor models developed by national laboratories and academia. The charter of the NRIC’s definition of advanced reactors also includes fusion nuclear reactors. As the tools developed by the NEAMS program are increasingly used for modeling fusion energy devices, there is an increasing need to host fusion reactor models on the VTB repository. The VTB will be extended in 2024 to support fusion energy modeling and simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

NucMesh: nuclear reactor geometry creation and mesh generation module in NEMoSys

NucMesh is a parameterized geometry and mesh generator for nuclear reactors developed within the Nuclear Energy Modeling System NEMoSys at Illinois Rocstar. NEMoSys is a platform developed for mesh generation, adaptive refinement, and solution verification. NucMesh is implemented to be generalized and extensible with a robust computer aided design engine and multiple mesh generation algorithms for unstructured triangular, quad-dominant, and structured quadrilateral meshing. In this paper, we present the geometric and meshing features of NucMesh. Geometrically objects are constructed bottom-to-top and overlaps are addressed automatically. A sophisticated object tracking algorithm prevents data from being lost for segmented objects. We discuss the primitive objects of circle and polygons that constitute the module and show how they are used with example inputs. Arrays of primitives and arrays of arrays are utilized to build large assemblies of objects. The concept of saved objects is discussed to demonstrate how repetitive objects can be reused easily and augmented in place. Three dimensional meshes can be obtained through mesh extrusion where all materials and side sets are extended to three dimensions. We show that side sets can be defined nearly anywhere within the geometry and can then be applied to the mesh. Finally, example reactor meshes are demonstrated for the Idaho National Laboratory Advanced Test Reactor and Los Alamos National Lab Empire reactor, both of which use control drums that NucMesh handles easily. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Understanding the role of negative charge in the scaffold of an artificial enzyme for CO 2 hydrogenation on catalysis

Here, we have approached the construction of an artificial enzyme by employing a robust protein scaffold, lactococcal multidrug resistance regulator, LmrR, providing a structured secondary and outer coordination spheres around a molecular rhodium complex, [Rh I (P Et2 N gly P Et2 ) 2 ] - . Previously, we demonstrated a 2–3 fold increase in activity for one Rh-LmrR construct by introducing positive charge in the secondary coordination sphere. In this study, a series of variants was made through site-directed mutagenesis where the negative charge is located in the secondary sphere or outer coordination sphere, with additional variants made with increasingly negative charge in the outer coordination sphere while keeping a positive charge in the secondary sphere. Placing a negative charge in the secondary or outer coordination sphere demonstrates decreased activity by a factor of two compared to the wild-type Rh-LmrR. Interestingly, addition of positive charge in the secondary sphere, with the negatively charged outer coordination sphere restores activity. Vibrational and NMR spectroscopy suggest minimal changes to the electronic density at the rhodium center, regardless of inclusion of a negative or positive charge in the secondary sphere, suggesting another mechanism is impacting catalytic activity, explored in the discussion.

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

Design of a supervisory control system for autonomous operation of advanced reactors

Advanced reactors to be deployed in the coming decades will face deregulated energy markets, and may adopt flexible operation to boost profitability. To aid in the transition from baseload to flexible operation paradigm, autonomous operation is sought. This work focuses on the control aspect of autonomous operation. Specifically, a hierarchical control system is designed to support constraint enforcement during routine operational transients. Within the system, data-driven modeling, physics-based state observation, and classical control algorithms are integrated to provide an adaptable and robust solution. A 320 MW Fluoride-cooled High-temperature Pebble-bed Reactor is the design basis for demonstrating the proposed control system. The hierarchical control system consists of a supervisory layer and low-level layer. The supervisory layer receives requests to change the system's operating conditions (e.g., the current reactor power to meet a load -follow), and accepts or rejects them based on constraints that have been assigned. Constraints are issued to keep the plant within an optimal operating region. The low-level layer interfaces with the actuators of the system to fulfill requested changes, while maintaining tracking and regulation duties. Further, to accept requests at the supervisory layer, the Reference Governor algorithm was adopted. To model the dynamics of the reactor, a system identification algorithm, Dynamic Mode Decomposition, was utilized. To estimate the evolution of process variables that cannot be directly measured (e.g., the propagation of delayed neutron precursors), the Unscented Kalman Filter, incorporating a nonlinear model of nuclear dynamics, was adopted. The composition of these algorithms led to a numerical demonstration of constraint enforcement during a 40% power drop transient (at a rate of 5 %/min). Uncontrolled secondary-side temperatures were successfully constrained. Adaptability of the proposed system was demonstrated by modifying the constraint values, and enforcing them during the transient. Robustness was also demonstrated by enforcing constraints under noisy environments.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2025

Butterfly valve performance factors using the multiphysics object oriented simulation environment

Butterfly valves are typically used in nuclear reactors to control incompressible fluid flow with high inlet velocities. Performance factors for butterfly valves include the pressure drop across the valve and the loss coefficient from which hydrodynamic torque and flow coefficients can be computed. This work explores a computational fluid dynamics approach for butterfly valve performance factors using the open-source Multiphysics Object Oriented Simulation Environment (MOOSE) framework. While MOOSE is often used in the nuclear energy modeling and simulation community for simulations ranging from fuel characterization to heat pipe simulation, this work employs the MOOSE open-source Navier–Stokes solver capability for simulating butterfly valve performance factors and compares those to experimentally measured results from the Advanced Test Reactor at Idaho National Laboratory at Reynolds numbers in the order of 10 6 for the partially opened configuration. The MOOSE framework results are compared against experimentally measured butterfly valve performance factors across five valve opening angles using meshes with order 10 4 – 10 5 elements. This validation serves to enable MOOSE-based multiphysics simulations incorporating the open-source Navier–Stokes module.

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Beyond-mean-field approaches for nuclear neutrinoless double beta decay in the standard mechanism

Nuclear weak decays provide important probes to fundamental symmetries in nature. A precise description of these processes in atomic nuclei requires comprehensive knowledge on both the strong and weak interactions in the nuclear medium and on the dynamics of quantum many-body systems. In particular, an observation of the hypothetical double beta decay without emission of neutrinos (0νββ) would unambiguously demonstrate the Majorana nature of neutrinos and the existence of the lepton-number-violation process. It would also provide unique information on the ordering and absolute scale of neutrino masses. The next-generation tonne-scale experiments with sensitivity up to 10 28 years after a few years of running will probably provide a definite answer to these fundamental questions based on our current knowledge on the nuclear matrix element (NME), the precise determination of which is a challenge to nuclear theory. Beyond-mean-field approaches have been frequently adapted for the study of nuclear structure and decay throughout the nuclear chart for several decades. In this review, we summarize the status of beyond-mean-field calculations of the NMEs of 0νββ decay assuming the standard mechanism of an exchange of light Majorana neutrinos. In conclusion, the challenges and prospects in the extension and application of beyond-mean-field approaches for 0νββ decay are discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

3.0 - MOOSE: Enabling massively parallel multiphysics simulations

The development of MOOSE has kept accelerating since the last release, with over 2,100 pull requests merged over the last 30 months that involved nearly fifty contributors across close to a dozen institutions internationally. The growth in MOOSE's capabilities and downstream applications is reflected in the growth of the community. User support provided on the GitHub discussions forum has steadily increased to nearly 50 daily interactions. New simulation projects, notably to model advanced nuclear reactor and fusion devices, are driving a significant expansion of the capabilities. This paper reports on these developments, with several major released features, new physics modules, and key improvements to the user experience and simulation workflow.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Tracing baryon and electric charge transport in isobar collisions

It is of fundamental interest to understand the carrier of conserved quantum charges within protons and nuclei at high energy. Preliminary data from isobar collisions at RHIC reveal a scaled net-baryon to net-electric charge ratio (B/ΔQ × ΔZ/A) at midrapidity between 1.2 and 2, consistent with string junction model predictions. Here, we compute the initial stage scaled net baryon to net-electric charge ratio for isobar collisions. Our model incorporates a realization of the string junction model and models the nuclear structure. Our predictions identify the baseline expectations for such measurement and quantify the impact of the nuclear structure.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Mixed quantum–classical approach to model non-adiabatic electron–nuclear dynamics: Detailed balance and improved surface hopping method

We develop a density matrix formalism to describe coupled electron–nuclear dynamics. To this end, we introduce an effective Hamiltonian formalism that describes electronic transitions and small (quantum) nuclear fluctuations along a classical trajectory of the nuclei. Using this Hamiltonian, we derive equations of motion for the electronic occupation numbers and for the nuclear coordinates and momenta. We show that, in the limit, when the number of nuclear degrees of freedom coupled to a given electronic transition is sufficiently high (i.e., the strong decoherence limit), the equations of motion for the electronic occupation numbers become Markovian. Furthermore, the transition rates in these (rate) equations are asymmetric with respect to the lower-to-higher energy transitions and vice versa. In thermal equilibrium, such asymmetry corresponds to the detailed balance condition. We also study the equations for the electronic occupations in the non-Markovian regime and develop a surface hopping algorithm based on our formalism. To treat the decoherence effects, we introduce additional “virtual” nuclear wave packets whose interference with the “real” (physical) wave packets leads to the reduction in coupling between the electronic states (i.e., decoherence) as well as to the phase shifts that improve the accuracy of the numerical approach. Remarkably, the same phase shifts lead to the detailed balance condition in the strong decoherence limit.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Structure of beta-decaying states in the deformed, neutron-rich nucleus 104 Nb

Excited structures in 104 Mo were populated by β decays of the ground and isomeric states in the neutron-rich nucleus 104 Nb. The beams were produced by the CARIBU facility at Argonne National Laboratory, re-accelerated by the ATLAS accelerator and implanted on a moving-tape system in the middle of the GAMMASPHERE array. Separate decay schemes for the two β-decaying states in 104 Nb were constructed for the first time. The structure of the isomers are discussed in the framework of the deformed Nilsson model and systematics of known quasiparticle structures in neighboring nuclei.

Nilsson model↗
DOE OSTI2026

Evolution of chirality from transverse wobbling in 135 Pr

Chirality is a distinct signature that characterizes triaxial shapes in nuclei. We report the first observation of chirality in the nucleus 135 Pr using a high-statistics Gammasphere experiment with the 123 Sb( 16 O, 4n) 135 Pr reaction. Two chiral-partner bands with the configuration π(1h11/2 )1 ⊗ ν(1h 11/2 ) −2 have been identified in this nucleus. Angular distribution analyses of the ΔI = 1 transitions connecting the two bands reveal a dominant dipole character, and quasiparticle triaxial rotor model calculations show good agreement with the data. Since the simultaneous observation of chirality and transverse wobbling in 135 Pr relies critically on these angular distribution results, we also address and refute the experimental and theoretical criticisms raised in a recent work by Lv et al., presenting additional evidence that further strengthens our interpretation. Furthermore, this marks the first observation of both hallmarks of triaxiality—chirality and wobbling—in the same nucleus.

90 ≤ A ≤ 149↗
DOE OSTI2022

Alternative approach to quantum imaginary time evolution

There is increasing interest in quantum algorithms (QAs) that are based on the imaginary time evolution (ITE), a successful classical numerical approach to obtain ground states. However, most of the proposals so far require heavy postprocessing computational steps on a classical computer, such as solving linear equations. Here we provide an alternative approach to implement ITE. A key feature in our approach is the use of an orthogonal basis set: the propagated state is efficiently expressed in terms of orthogonal basis states at every step of the evolution. We argue that the number of basis states needed at those steps to achieve an accurate solution can be kept on the order of n , the number of qubits, by controlling the precision (number of significant digits) and the imaginary time increment. The number of quantum gates per imaginary time step is estimated to be polynomial in n . Additionally, while in many QAs the locality of the Hamiltonian is a key assumption, in our algorithm this restriction is not required. This characteristic of our algorithm renders it useful for studying highly nonlocal systems, such as the occupation-representation nuclear shell model. Here, we illustrate our algorithm through numerical implementation on an IBM quantum simulator.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Quasielastic lepton scattering and back-to-back nucleons in the short-time approximation

Understanding quasielastic electron- and neutrino-scattering from nuclei has taken on new urgency with current and planned neutrino oscillation experiments, and with electron scattering experiments measuring specific final states, such as those involving nucleon pairs in “back-to-back” configurations. Accurate many-body methods are available for calculating the response of light (A ≤ 12) nuclei to electromagnetic and weak probes, but they are computationally intensive and only applicable to the inclusive response. In the present work we introduce a novel approach, based on realistic models of nuclear interactions and currents, to evaluate the short-time (high-energy) inclusive and exclusive response of nuclei. The approach accounts reliably for crucial two-nucleon dynamics, including correlations and currents, and provides information on back-to-back nucleons observed in electron and neutrino scattering experiments. We demonstrate that in the quasielastic regime and at moderate momentum transfers both initial- and final-state correlations, and two-nucleon currents are important for a quantitatively successful description of the inclusive response and final state nucleons. Finally, the approach can be extended to include relativistic—kinematical and dynamical—effects, at least approximately in the two-nucleon sector, and to describe the response in the resonance-excitation region.

6 ≤ A ≤ 19↗
DOE OSTI2024

Effect of causality constraints on Bayesian analyses of heavy-ion collisions

There have long been questions about the limits to the validity of relativistic fluid dynamics and whether it is being used outside its regime of validity in modern simulations of relativistic heavy-ion collisions. An important new tool for answering this question is a causality analysis in the nonlinear regime—if the solutions of the evolution equations do not respect relativistic causality, then they are not a faithful representation of the underlying relativistic theory (in this case, quantum chromodynamics). Using this nonlinear criterion, it has recently been shown that hydrodynamics is indeed being used outside its regime of validity in simulations, at least sometimes. Here we explore the phenomenological implications, particularly the quantitative effects of demanding limits on acausality in modern Bayesian parameter estimation. We find that, while typically only a small fraction of the system's energy is initially in an acausal regime, placing strict limits on the allowed energy fraction significantly changes the preferred properties of the initial condition, which in turn alters the extracted medium properties such as bulk viscosity, where large values are no longer favored. Furthermore, these findings highlight the importance of developing better theoretical descriptions of the early-time, out-of-equilibrium dynamics of relativistic heavy-ion collisions.

Bayesian methods↗
DOE OSTI2021

Speed of Sound and Baryon Cumulants in Heavy-Ion Collisions

We present a method that may allow an estimate of the value of the speed of sound as well as its logarithmic derivative with respect to the baryon number density in matter created in heavy-ion collisions. To this end, we utilize well-known observables: cumulants of the baryon number distribution. In analyses aimed at uncovering the phase diagram of strongly interacting matter, cumulants gather considerable attention as their qualitative behaviour along the explored range of collision energies is expected to aid in detecting the QCD critical point. We show that the cumulants may also reveal the behavior of the speed of sound in the temperature and baryon chemical potential plane. We demonstrate the applicability of such estimates within two models of nuclear matter, and explore what might be understood from known experimental data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

First nuclear ultraheavy dark matter search in argon time projection chambers with the DarkSide-50 experiment

We report on the first search for nuclear ultraheavy dark matter (UHDM) in a dual-phase liquid argon time projection chamber using the DarkSide-50 experiment. Unlike conventional weakly interacting massive particles, a nuclear UHDM model may be composed of many dark nucleons and scatter numerous times while passing through the detector. Accounting for energy loss through the Earth’s overburden, we apply selection criteria optimized for multiscatter event topologies using the 532-day low-radiation campaign of the DarkSide-50 detector. Excluded limits on the UHDM-nucleon scattering cross section for dark nucleon masses of m χ = 10 , 50 , 100 , 500 GeV / c 2 are presented.

Agnes, P. [GSSI, Aquila; Gran Sasso]↗
DOE OSTI2025

Microscopic optical potentials from a Green's function approach

Optical potentials are a standard tool in the study of nuclear reactions, as they describe the interaction between a target nucleus and a projectile. The use of phenomenological optical potentials built using experimental data on stable isotopes is widespread. Although successful in their dedicated domain, it is unclear whether these phenomenological potentials can provide reliable predictions for unstable isotopes. To address this problem, optical potentials based on microscopic nuclear structure input calculations prove to be crucial and are an important current line of research. In this work we present an explicit implementation of the Feshbach formalism for the systematic derivation of optical potentials using input from nuclear structure models. Numerical tools for the derivation of Green's functions associated with nonlocal potentials are presented. In conclusion, the new optical potential, based on the valence shell model, is applied to the calculations of 𝑛 + 24 Mg elastic scattering and yields a close agreement with the experimental data.

Direct reactions↗
DOE OSTI2022

Radiation effects on the structure and alteration behavior of an SiO 2 –Al 2 O 3 –B 2 O 3 –Na 2 O glass

As borosilicate glasses are used in many countries to immobilize fission products and minor actinides after spent fuel reprocessing before storage in a deep geological repository, assessing that their chemical durability is of paramount importance. Here, pristine and preirradiated (952 MeV, 136Xe) SiO 2 –B2O 3 –Al2O 3 –Na 2 O glasses with the same molar ratios as in the French SON68 and ISG glasses have been subjected to aqueous corrosion in deionized water and in silica-saturated solution to measure the initial and longer term alteration rates. Pristine and preirradiated glasses corrode following the same mechanisms, but the preirradiation has a strong impact on the initial dissolution rate (increase by a factor of 5.6), and on the alteration layer depth in silica-saturated conditions (by two- to threefolds). The later result is related to the formation of a more porous, less passivating gel on the preirradiated glass specimen. Using both experimental spectroscopies (NMR, IR, and SFG) and classical molecular dynamics, the radiation effects on the glass structure and water diffusion have been assessed. After preirradiation, the density and the polymerization degree of the glass decrease, whereas the topological disorder increases. In consequence, water diffusion accelerates. These observations allow to correlate the radiation impact on the alteration behavior to the structural changes.

36 MATERIALS SCIENCE↗