SEARCH · Engineering Papers

Results for “Models & methods for nuclear reactions”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 289 records · Page 16

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Modelling populations of kilonovae

Abstract The 2017 detection of a kilonova coincident with gravitational-wave emission has identified neutron star mergers as the major source of the heaviest elements and dramatically constrained alternative theories of gravity. Observing a population of such sources has the potential to transform cosmology, nuclear physics, and astrophysics. However, with only one confident multi-messenger detection currently available, modelling the diversity of signals expected from such a population requires improved theoretical understanding. In particular, models that are quick to evaluate and are calibrated with more detailed multi-physics simulations are needed to design observational strategies for kilonovae detection and to obtain rapid-response interpretations of new observations. We use grey-opacity models to construct populations of kilonovae, spanning ejecta parameters predicted by numerical simulations. Our modelling focuses on wavelengths relevant for upcoming optical surveys, such as the Rubin Observatory Legacy Survey of Space and Time (LSST). In these simulations, we implement heating rates that are based on nuclear reaction network calculations. We create a Gaussian-process emulator for kilonova grey opacities, calibrated with detailed radiative transfer simulations. Using recent fits to numerical relativity simulations, we predict how the ejecta parameters from binary neutron star (BNS) mergers shape the population of kilonovae, accounting for the viewing-angle dependence. Our simulated population of BNS mergers produce peak i-band absolute magnitudes of −20 ≤ Mi ≤ −11. A comparison with detailed radiative transfer calculations indicates that further improvements are needed to accurately reproduce spectral shapes over the full light curve evolution.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Neural Networks for Nuclear Reactions in MAESTROeX

We demonstrate the use of neural networks to accelerate the reaction steps in the MAESTROeX stellar hydrodynamics code. A traditional MAESTROeX simulation uses a stiff ODE integrator for the reactions; here, we employ a ResNet architecture and describe details relating to the architecture, training, and validation of our networks. Our customized approach includes options for the form of the loss functions, a demonstration that the use of parallel neural networks leads to increased accuracy, and a description of a perturbational approach in the training step that robustifies the model. We test our approach on millimeter-scale flames using a single-step, 3-isotope network describing the first stages of carbon fusion occurring in Type Ia supernovae. We train the neural networks using simulation data from a standard MAESTROeX simulation, and show that the resulting model can be effectively applied to different flame configurations. This work lays the groundwork for more complex networks, and iterative time-integration strategies that can leverage the efficiency of the neural networks.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Thermonuclear Reaction Rates and Primordial Nucleosynthesis

Assuming the best numerical value for the cosmic baryonic density and the existence of three neutrino flavors, standard Big Bang nucleosynthesis is a parameter-free model. It is important to assess if the observed primordial abundances can be reproduced by simulations. Numerous studies have shown that the simulations overpredict the primordial 7 Li abundance by a factor of ≈3 compared to the observations. The discrepancy may be caused by unknown systematics in 7 Li observations, poorly understood depletion of lithium in stars, errors in thermonuclear rates that take part in the lithium and beryllium synthesis, or physics beyond the standard model. Here, we focus on the likelihood of a nuclear physics solution. The status of the key nuclear reaction rates is summarized. Big Bang nucleosynthesis simulations are performed with the most recent reaction rates, and the uncertainties of the predicted abundances are established using a Monte Carlo technique. Correlations between abundances and reaction rates are investigated based on the metric of mutual information. The rates of four reactions impact the primordial 7 Li abundance: 3 He(α,γ) 7 Be, d(p,γ) 3 He, 7 Be(d,p)2α, and 7 Be(n,p) 7 Li. We employ a genetic algorithm to search for simultaneous rate changes in these four reactions that may account for all observed primordial abundances. Furthermore, when the search is performed for reaction rate ranges that are much wider than recently reported uncertainties, no acceptable solutions are found. Based on the currently available evidence, we conclude that it is highly unlikely for the cosmological lithium problem to have a nuclear physics solution.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Disruption prediction with artificial intelligence techniques in tokamak plasmas

In nuclear fusion reactors, plasmas are heated to very high temperatures of more than 100 million kelvin and, in so-called tokamaks, they are confined by magnetic fields in the shape of a torus. Light nuclei, such as deuterium and tritium, undergo a fusion reaction that releases energy, making fusion a promising option for a sustainable and clean energy source. Tokamak plasmas, however, are prone to disruptions as a result of a sudden collapse of the system terminating the fusion reactions. As disruptions lead to an abrupt loss of confinement, they can cause irreversible damage to present-day fusion devices and are expected to have a more devastating effect in future devices. Disruptions expected in the next-generation tokamak, ITER, for example, could cause electromagnetic forces larger than the weight of an Airbus A380. Furthermore, the thermal loads in such an event could exceed the melting threshold of the most resistant state-of-the-art materials by more than an order of magnitude. In this work, to prevent disruptions or at least mitigate their detrimental effects, empirical models obtained with artificial intelligence methods, of which an overview is given here, are commonly employed to predict their occurrence—and ideally give enough time to introduce counteracting measures.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Fermion Pair Dynamics in Open Quantum Systems

Three-body decay is a rare decay mode observed in a handful of unbound rare isotopes. The angular and energy correlations between emitted nucleons are of particular interest, as they provide invaluable information on the interplay between structure and reaction aspects of the nuclear open quantum system. To study the mechanism of two-nucleon emission, we developed a time-dependent approach that allows us to probe emitted nucleons at long times and large distances. We successfully benchmarked the new method against the Green’s function approach and applied it to low-energy two-proton and two-neutron decays. In particular, we studied the interplay between initial-state nucleon-nucleon correlations and final-state interaction. We demonstrated that the time evolution of the two-nucleon wave function is strongly impacted by the diproton/dineutron dynamics and that the correlations between emitted nucleons provide invaluable information on the dinucleon structure in the initial state.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

Diffusion quantum Monte Carlo approach to the polaritonic ground state

Making and using polaritonic states (i.e., hybrid electron-photon states) for chemical applications has recently become one of the most prominent and active fields that connects the communities of chemistry and quantum optics. Modeling of such polaritonic phenomena using ab initio approaches calls for new methodologies, leading to the reinvention of many commonly used electronic structure methods, such as Hartree-Fock, density functional, and coupled cluster theories. Here in this work, we explore the formally exact diffusion quantum Monte Carlo approach to obtain numerical solutions to the polaritonic ground state during the dissociation of the H 2 molecular system. Here, we examine various electron-nuclear-photon properties throughout the dissociation, such as changes to the minimum of the cavity Born-Oppenheimer surface, the localization of the electronic wave function, and the average mode occupation. Finally, we directly compare our results to that obtained with state-of-the-art, yet approximate, polaritonic coupled cluster approaches.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2022

Development of a three-dimensional APOLLO3 neutrons deterministic scheme for the CABRI reactor

CABRI is an experimental reactor to study the fuel behavior during reactivity injection transients. These transients being highly multiphysics, the development of suitable modeling and simulation tools to simulate them is important for the optimization of the tests and the control of the experimental conditions. This paper focuses on the development of an APOLLO3 deterministic core calculation tool dedicated to the CABRI transient analysis. It represents the first stage of the incremental process for the implementation of a multiphysics time-dependent modeling of the CABRI transient. The neutron calculation scheme is based on a classical two-step approach. The first step consists of a 281-energy group calculation flux with the TDT-MOC (Method Of Characteristics) solver for cross-section space and energy (23 groups) collapsing for the CABRI different assembly clusters. The bias on a 2D core neutron calculation due to the self-shielding calculation and collapsing on a restricted pattern are investigated thanks to a comparison with a direct full 2D calculation on a quarter of core. The second step relies on a pin-resolved transport 3D transport core calculation with the SN solver MINARET. A progressive numerical validation process is followed to quantify the calculation biases on reactivity and reaction rates at each step using reference calculations with the stochastic code TRIPOLI4. The next development stage toward a multiphysics scheme will be the implementation of the 3D-kinetics equation resolution and the coupling with a core thermal-hydraulics model. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2021

Search for in-band transitions in the candidate superdeformed band in Si 28

Background: Superdeformed (SD) bands are suggested by theory around Ca 40 and in lighter alpha-conjugate nuclei such as Mg 24 , Si 28 , and S 32 . Such predictions originate from a number of theoretical models including mean-field models and antisymmetrized molecular dynamics (AMD) calculations. While SD bands have been identified in Ca 40 and its near neighbors, evidence of their existence in the lighter, midshell nuclei is circumstantial at best. Additionally, the key evidence of superdeformation would be the observation of transitions with high B ( E 2 ) transition strengths connecting states in a rotational sequence. This is challenging information to obtain since the bands lie at a high excitation energy and competition from out-of-band decay is dominant. Purpose: The purpose of the present study is to establish a new methodology to circumvent the difficulties in identifying and quantifying in-band transitions through directly populating candidate states in the SD band in Si 28 through inelastic alpha scattering, selecting such states with a spectrometer, and measuring their gamma-ray decay with a large array of high-purity germanium detectors, allowing direct access to electromagnetic transition strengths. Methods: Excited states in Si 28 were populated in the Si 28 ( α , α ' ) reaction using a 130-MeV He 4 beam from the K140 AVF cyclotron at the Research Center for Nuclear Physics. Outgoing alpha particles were analyzed using the Grand Raiden spectrometer positioned at an angle of 9 . 1 ° to favor the population of states with J ≈ 4 . Coincident gamma rays were detected with the CAGRA array of 12 HPGe clover detectors augmented by a set of four large LaBr 3 detectors. Results: Data analysis showed that it was possible to identify additional low-energy transitions in competition with high-energy decays from excited states in Si 28 in the vicinity of 10 MeV. However, while the candidate 4 + SD state at 10.944 MeV was populated, a 1148-keV transition to the candidate 2 + SD state at 9.796 MeV was not observed, and only an upper limit for its transition strength of B ( E 2 ) < 43 W.u. could be established. This contradicts AMD predictions of ≈ 200 W.u. for such a transition. Conclusion: The present study strongly rejects the hypothesis that the candidate set of states identified in Si 28 represents an SD band, which demonstrates the potential of the methodology devised here.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS2015

Investigation of HZETRN 2010 as a Tool for Single Event Effect Qualification of Avionics Systems - Part II

An accurate prediction of spacecraft avionics single event effect (SEE) radiation susceptibility is key to ensuring a safe and reliable vehicle. This is particularly important for long-duration deep space missions for human exploration where there is little or no chance for a quick emergency return to Earth. Monte Carlo nuclear reaction and transport codes such as FLUKA can be used to generate very accurate models of the expected in-flight radiation environment for SEE analyses. A major downside to using a Monte Carlo-based code is that the run times can be very long (on the order of days). A more popular choice for SEE calculations is the CREME96 deterministic code, which offers significantly shorter run times (on the order of seconds). However, CREME96, though fast and easy to use, has not been updated in several years and underestimates secondary particle shower effects in spacecraft structural shielding mass. Another modeling option to consider is the deterministic code HZETRN 20104, which includes updates to address secondary particle shower effects more accurately. This paper builds on previous work by Rojdev, et al. to compare the use of HZETRN 2010 against CREME96 as a tool to verify spacecraft avionics system reliability in a space flight SEE environment. This paper will discuss modifications made to HZETRN 2010 to improve its performance for calculating SEE rates and compare results with both in-flight SEE rates and other calculation methods.

Rojdev, Kristina↗

DOE OSTI2021

Development of New Experimental Methods for Correlated Operando Surface/Gas Characterization

The predictive understanding of catalytic surface reactions requires accurate microkinetic models, and while decades of work has been devoted to the elucidation of the reaction steps in these models, many open questions remain. One key issue is a lack of approaches enabling the local spatially resolved assessment of catalytic activity over a surface. In this report, we detail efforts to develop a new diagnostic approach to solve this problem. The approach is based upon laser resonance enhanced multiphoton ionization of reaction products emitted into the gas phase followed by spatially resolved imaging of the resultant ions or electrons. Ion imaging is pursued with a velocity-selected spatially resolved ion imaging microscope, while electron imaging was attempted in a low energy electron microscope. Successful demonstration of the ion imaging microscope coupled with the development of transport simulations shows promise for a revolutionary new tool to assess local catalytic activity

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS2002

Hierarchies of Models: Toward Understanding Planetary Nebulae

Stars like our sun (initial masses between 0.8 to 8 solar masses) end their lives as swollen red giants surrounded by cool extended atmospheres. The nuclear reactions in their cores create carbon, nitrogen and oxygen, which are transported by convection to the outer envelope of the stellar atmosphere. As the star finally collapses to become a white dwarf, this envelope is expelled from the star to form a planetary nebula (PN) rich in organic molecules. The physics, dynamics, and chemistry of these nebulae are poorly understood and have implications not only for our understanding of the stellar life cycle but also for organic astrochemistry and the creation of prebiotic molecules in interstellar space. We are working toward generating three-dimensional models of planetary nebulae (PNe), which include the size, orientation, shape, expansion rate and mass distribution of the nebula. Such a reconstruction of a PN is a challenging problem for several reasons. First, the data consist of images obtained over time from the Hubble Space Telescope (HST) and spectra obtained from Kitt Peak National Observatory (KPNO) and Cerro Tololo Inter-American Observatory (CTIO). These images are of course taken from a single viewpoint in space, which amounts to a very challenging tomographic reconstruction. Second, the fact that we have two disparate and orthogonal data types requires that we utilize a method that allows these data to be used together to obtain a solution. To address these first two challenges we employ Bayesian model estimation using a parameterized physical model that incorporates much prior information about the known physics of the PN. In our previous works we have found that the forward problem of the comprehensive model is extremely time consuming. To address this challenge, we explore the use of a set of hierarchical models, which allow us to estimate increasingly more detailed sets of model parameters. These hierarchical models of increasing complexity are akin to scientific theories of increasing sophistication, with each new model/theory being a refinement of a previous one by either incorporating additional prior information or by introducing a new set of parameters to model an entirely new phenomenon. We apply these models to both a simulated and a real ellipsoidal PN to initially estimate the position, angular size, and orientation of the nebula as a two-dimensional object and use these estimates to later examine its three-dimensional properties. The efficiency/accuracy tradeoffs of the techniques are studied to determine the advantages and disadvantages of employing a set of hierarchical models over a single comprehensive model.

Knuth, Kevin H.↗

NASA NTRS2003

Hierarchies of Models: Toward Understanding Planetary Nebulae

Knuth, Kevin H.↗

DOE OSTI2020

Recent advances in the quantification of uncertainties in reaction theory

Uncertainty quantification has become increasingly more prominent in nuclear physics over the past several years. In few-body reaction theory, there are four main sources that contribute to the uncertainties in the calculated observables: the effective potentials, approximations made to the few-body problem, structure functions, and degrees of freedom left out of the model space. In this work, we illustrate some of the features that can be obtained when modern statistical tools are applied in the context of nuclear reactions. This work consists of a summary of the progress that has been made in quantifying theoretical uncertainties in this domain, focusing primarily on those uncertainties coming from the effective optical potential as well as their propagation within various reaction theories. We use, as the central example, reactions on the doubly-magic stable nucleus 40 Ca, namely neutron and proton elastic scattering and single-nucleon transfer 40 Ca(d,p) 41 Ca. First, we show different optimization schemes used to constrain the optical potential from differential cross sections and other experimental constraints; we then discuss how these uncertainties propagate to the transfer cross section, comparing two reaction theories. Finally, we finish by laying out our future plans.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2020

Electrochemical Impedance Spectroscopy and Cyclic Voltammetry Methods for Monitoring SmCl 3 Concentration in Molten Eutectic LiCl-KCl

Molten salt solutions consisting of eutectic LiCl-KCl and concentrations of samarium chloride (0.5 to 3.0 wt%) at 500 C were analyzed using both cyclic voltammetry (CV) and electrochemical impedance spectroscopy (EIS). The CV technique gave the average diffusion coefficient for Sm 3+ over the concentration range. Equipped with Sm 3+ diffusion coefficient, the Randles-Sevcik equation predicted Sm 3+ concentration values that agree with the given experimental values. From CV measurements; the anodic, cathodic, and half-peak potentials were identified and subsequently used as a parameter to acquire EIS spectra. A six-element Voigt model was used to model the EIS data in terms of resistance-time constant pairs. The lowest resistances were observed at the half-peak potential with the associated resistance-time constant pairs characterizing the reversible reaction between Sm 3+ and Sm 2+ . By extrapolation, the Voigt model estimated the polarization resistance and established a polarization resistance-concentration relationship.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2022

Self-consistent description of the halo nature of 31 Ne with continuum and pairing correlations

Here a relativistic structure model has previously been used to predict a halo structure for 31 Ne, consistent with halo signatures from measured reaction cross sections of Ne isotopes bombarding Carbon targets. However, previous attempts to calculate those cross sections with reaction models were missing contributions from resonances and pairing correlations in their structure input. Use a reaction model with our relativistic fully microscopic structure model input to predict these cross sections and momentum distributions and analyze for possible halo signatures. Structure input for exotic Ne isotopes were obtained via the analytical continuation of the coupling constant (ACCC) method based on the relativistic mean field (RMF) theory with Bardeen - Cooper - Schrieffer (BCS) pairing approximation, the RAB approach. Total reaction cross sections, one-neutron removal cross sections, and momentum distributions of breakup reaction products were calculated with a Glauber model using our relativistic structure input. Our predictions of total reaction and one-neutron removal cross sections of 31 Ne on a Carbon target were significantly enhanced compared with those of neighboring Neon isotopes, agreeing with measurements at 240 MeV/nucleon and consistent with a single neutron halo. Furthermore, our calculations of the inclusive longitudinal momentum distribution of the 30 Ne and valence neutron residues from the 31 Ne breakup reaction indicate a dilute density distribution in coordinate space, another halo signature. We give a full description of the halo nature of 31 Ne that includes a self-consistent use of pairing and continuum contributions that makes predictions consistent with reaction cross section measurements. This approach can be utilized to determine the halo structure of other exotic nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

H 2 production in the 10 B(n,α) 7 Li reaction in water

Herein we demonstrate a method for measuring the H2 produced in water from the 10 B(n,α) 7 Li fission reaction. Low energy neutrons from the NIST Center for Neutron Research interact with borate-containing water in a temperature-controlled high pressure cell made from titanium. After exposure for one to several hours, the water is extracted and sparged with argon. H 2 entrained in the sparging gas is sampled with a small mass spectrometer. To determine the neutron exposure, a small amount of sodium is included in the borate solution. The water is collected and 24 Na activation is measured in a counting apparatus on the following day. The G-value for H 2 at room temperature is found to be (1.18 ± 0.10) molecules H 2 /100eV, in good agreement with previous estimates and recent modeling calculations.

10B↗

DOE OSTI2020

Low-lying level structure of the neutron-unbound $\textit{N}$ = 7 isotones

This article reports on results of an experimental search for ground state shell inversion in the neutron-unbound $\textit{N}$ = 7 isotones 9 He and 10 Li. Two different radioactive ion beams ( 11 Be, E lab = 44 MeV/u and 12 B, E lab = 45 MeV/u) impinging on a beryllium target were used to directly and selectively populate unbound states of a given ℓ in the nuclides of interest. Be( 11 Be, 8 He + $\textit{n}$) and Be( 12 B, 8 He + $\textit{n}$) reactions populated unbound states in 9 He. Fragments and neutrons were detected in coincidence to reconstruct the decay energy of 9 He using invariant mass spectroscopy. Similarly, Be( 11 Be, 9 Li + $\textit{n}$) and Be( 12 B, 9 Li + $\textit{n}$) reactions were used to populate unbound states in 10 Li, and the time-of-flight method was used to determine the relative velocity of the 9 Li fragments and neutrons in coincidence. Various states in both 9 He and 10 Li were observed and characterized. Here, the data indicate possible ground state shell inversion in both cases.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2023

Exascale Multiphysics Nuclear Reactor Simulations for Advanced Designs

ENRICO is a coupled application developed under the U.S. Department of Energy's Exascale Computing Project (ECP) targeting the modeling of advanced nuclear reactors. It couples radiation transport with heat and fluid simulation, including the high-fidelity, highresolution Monte-Carlo code Shift and the Computational fluid dynamics code NekRS. NekRS is a highly-performant open-source code for simulation of incompressible and low-Mach fluid flow, heat transfer, and combustion with a particular focus on turbulent flows in complex domains. It is based on rapidly convergent high-order spectral element discretizations that feature minimal numerical dissipation and dispersion. State-of-the-art multilevel preconditioners, efficient high-order time-splitting methods, and runtime-adaptive communication strategies are built on a fast OCCA-based kernel library, libParanumal, to provide scalability and portability across the spectrum of current and future high-performance computing platforms. On Frontier, Nek5000/RS has recently achieved an unprecedented milestone in breaching over 1 billion spectral elements and 350 billion degrees of freedom. Shift has demonstrated the capability to transport upwards of 1 billion particles per second in full core nuclear reactor simulations featuring complete temperature-dependent, continuous-energy physics on Frontier. Shift achieved a weak-scaling efficiency of 97.8% on 8192 nodes of Frontier and calculated 6 reactions in 214,896 fuel pin regions below 1% statistical error yielding first-of-a-kind resolution for a Monte Carlo transport application.

Hamilton, Steven P.↗