Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results 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.

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

At least 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Measuring and simulating ice–ablator mix in inertial confinement fusion

Fuel–ablator mix has been established as a major performance degrading effect in the burning plasma regime of recent inertial confinement fusion (ICF) experiments. As such, the study of fuel–ablator mix with experiments and simulations can provide valuable insight for our understanding of these experiments and establish a path for even higher yields and increased robustness. Here we present a novel high-yield experimental ICF design that is motivated by recent experiments measuring ice–ablator mix with a CH ablator instead of a high-density carbon (HDC) ablator. Here we review these experiments in more detail and describe the modeling assumptions and parameters used to obtain agreement with the data from implosion and burn simulations with mix. Using this mix model calibrated a posteriori to the experimental data, we design an implosion that uses a CH ablator that is predicted to achieve better performance than a recent experiment that achieved net target gain of 1.5 in HDC. Because hydrodynamic instabilities are greatly reduced with this new design, we also expect a high reproducibility at the same implosion adiabat as current record yield experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Supervised machine learning-based multivariate regression of parallel closures for a high-collisionality deuterium-carbon plasma

Many plasmas of interest in laboratory experiments and space consist of multiple ion species. In tokamak edge plasmas, for instance, ionized impurities expelled from the vessel wall influence plasma transport. When describing multi-species plasmas using fluid equations, we need accurate closure relations to close the set of fluid equations. In this study, we introduce the development of fitting formulas for parallel closures using supervised machine learning, in conjunction with the recent closure theory, considering multi-ion collisions and arbitrary ion temperatures. We apply this approach to a high-collisionality deuterium-carbon plasma and demonstrate its effectiveness. As a result, the machine learning-based method for developing practical and accurate closures can be extended to a wider range of plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Solid, structured composite neutron detectors with high dynamic range capability

Neutron detectors are essential across disciplines such as fundamental science, nuclear security, safeguards, and civilian applications. While 3 He-filled gas proportional counters have long been revered for their efficacy in detecting thermal neutrons and praised for their efficiency, neutron/gamma discrimination, and stability, the scarcity of 3He has spurred a search for alternatives. Here, we explore a solid structured scintillating particle composite (SPC) consisting of 6 Li-containing scintillating glass particles within an acrylic matrix as a neutron detector for high dynamic range applications. We show for the first time that an SPC neutron detector can boast an intrinsic detection efficiency of 0.261% for pure 252 Cf fission neutrons and an overall neutron detection efficiency of (0.546 ± 0.003)% at the Neutron Free-in-Air facility while being able to function in an intense gamma-ray environment. We also show that the SPC neutron detector supports fast neutron capture times and enables a dual-readout scheme that extends the detector dynamic range to high incident neutron fluxes. A scalable fabrication process allows for tailoring the SPC detector properties to the requirements of specific applications. Good agreement is found between the experimental results taken with a National Institute of Standards and Technology traceable 252 Cf source and the coupled MCNP6 and optical-ray-tracing simulations.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Sequential Bayesian Experimental Design for Calibration of Expensive Simulation Models

Simulation models of critical systems often have parameters that need to be calibrated using observed data. For expensive simulation models, calibration is done using an emulator of the simulation model built on simulation output at different parameter settings. Using intelligent and adaptive selection of parameters to build the emulator can drastically improve the efficiency of the calibration process. The article proposes a sequential framework with a novel criterion for parameter selection that targets learning the posterior density of the parameters. The emergent behavior from this criterion is that exploration happens by selecting parameters in uncertain posterior regions while simultaneously exploitation happens by selecting parameters in regions of high posterior density. Furthermore, the advantages of the proposed method are illustrated using several simulation experiments and a nuclear physics reaction model.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Advanced Isotope Separation Technology for Fusion Fuel

Deuterium-tritium fusion is the easiest nuclear fusion reaction among known fusion reactions. Since tritium is extremely rare, it is artificially produced by irradiating lithium metal. The separation, isolation, and storage of the tritium isotope has been a major focus of the Savannah River Site (SRS) for many decades. Thermal diffusion, fractional absorption, and cryogenic distillation have all been used in the past, and each has significant operational and safety challenges. A process known as the Thermal Cycling Absorption Process (TCAP) was invented at SRS, and because of its overwhelming advantages in safety, efficiency, size, and reduced tritium inventory, it has replaced all other hydrogen isotope separation processes at SRS. Here, the working principles and current development of hydrogen isotope separation using TCAP at SRS are explained as a potential advanced isotope separation process for the fusion fuel cycle.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Get on the BAND Wagon: a Bayesian framework for quantifying model uncertainties in nuclear dynamics

In this paper, we describe the Bayesian Analysis of Nuclear Dynamics (BAND) framework, a cyberinfrastructure that we are developing which will unify the treatment of nuclear models, experimental data, and associated uncertainties. We overview the statistical principles and nuclear-physics contexts underlying the BAND toolset, with an emphasis on Bayesian methodology's ability to leverage insight from multiple models. In order to facilitate understanding of these tools we provide a simple and accessible example of the BAND framework's application. Four case studies are presented to highlight how elements of the framework will enable progress on complex, far-ranging problems in nuclear physics. By collecting notation and terminology, providing illustrative examples, and giving an overview of the associated techniques, this paper aims to open paths through which the nuclear physics and statistics communities can contribute to and build upon the BAND framework.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Variational inference of effective range parameters for 3 He− 4 He scattering

We use two different methods, Monte Carlo sampling and variational inference (VI), to perform a Bayesian calibration of the effective-range parameters in 3 He– 4 He elastic scattering. The parameters are calibrated to data from a recent set of 3 He– 4 He elastic scattering differential cross section measurements. Analysis of these data for E lab ≤ 4.3 MeV yields a unimodal posterior for which both methods obtain the same structure. However, the effective-range expansion amplitude does not account for the 7/2 − state of 7 Be so, even after calibration, the description of data at the upper end of this energy range is poor. The data up to E lab = 2.6 MeV can be well described, but calibration to this lower-energy subset of the data yields a bimodal posterior. After adapting VI to treat such a multi-modal posterior we find good agreement between the VI results and those obtained with parallel-tempered Monte Carlo sampling.

effective field theory↗
DOE OSTI2025

Emulation of the calculations of final r -process abundance patterns with a neural network

This work explores the construction of a fast emulator for the calculation of the final pattern of nucleosynthesis in the rapid neutron capture process (the r-process). An emulator is built using a feed-forward artificial neural network (ANN). We train the ANN with nuclear data and relative abundance patterns. We take as input the β-decay half-lives and the one-neutron separation energy of the nuclei in the rare-earth region. The output is the final isotopic abundance pattern. In this work, we focus on the nuclear data and abundance patterns in the rare-earth region to reduce the dimension of the input and output space. We show that the ANN can capture the effect of the changes in the nuclear physics inputs on the final r-process abundance pattern in the adopted astrophysical conditions. We employ the deep ensemble method to quantify the prediction uncertainty of the neural network emulator. The emulator achieves a speed-up by a factor of about 20 000 in obtaining a final abundance pattern in the rare-earth region. The emulator may be utilized in statistical analyses such as uncertainty quantification, inverse problems, and sensitivity analysis.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Connecting ground-state properties of 6 Li to each other and to scattering data

We examine the relationship between the asymptotic normalization coefficient (ANC) of 6 Li and other low-energy observables in the α–deuteron system. Our analysis uses a set of calculations carried out within the ab initio no core shell model with continuum (NCSMC) using a variety of inter-nucleon interactions and basis sizes, and yielding 6 Li deuteron separation energies between 1.3 and 1.8 MeV (Hebborn et al 2022 Phys. Rev. Lett. 129 042503). These NCSMC calculations show that the square of the ANC is strongly correlated with the separation energy over this range. In this work, we investigate the origin of this correlation using the phenomenological R-matrix, a single-channel potential and a perturbative approach. We show that this correlation occurs because the depth of the α–deuteron central potential changes by only a small relative amount as the separation energy varies. We then investigate if the ANC can be accurately extracted from α–deuteron phase shifts in an ideal case in which low-energy data are available and there are no experimental errors. We find that both R-matrix and Coulomb-modified effective-range theory (CM-ERE) yield extracted ANCs close to, although not exactly equal to, the NCSMC value, provided the extrapolation is constrained by the known position of the bound-state pole and at least three terms are included in the fit function. The R-matrix approach converges faster than the CM-ERE as the number of parameters increases and is also more robust against the inclusion of low-energy and high-energy phase shift data. Finally, our study also shows that a naive quantification of uncertainties by comparing different truncations used in both theories is not accurate, and suggests the accuracy of ANCs extracted from phase shift data needs further investigation.

R-matrix↗
DOE OSTI2021

Radioactive nuclei in the early Solar system: analysis of the 15 isotopes produced by core-collapse supernovae

Short-lived radioactive isotopes (SLRs) with half-lives between 0.1 and 100 Myr can be used to probe the origin of the Solar system. In this work, we examine the core-collapse supernovae production of the 15 SLRs produced: 26 Al, 36 Cl, 41 Ca, 53 Mn, 60 Fe, 92 Nb, 97 Tc, 98 Tc, 107 Pd, 126 Sn, 129 I, 135 Cs, 146 Sm, 182 Hf, and 205 Pb. Here, we probe the impact of the uncertainties of the core-collapse explosion mechanism by examining a collection of 62 core-collapse models with initial masses of 15, 20, and 25 M ⊙ , explosion energies between 3.4 × 10 50 and 1.8 × 10 52 erg and compact remnant masses between 1.5 and 4.89 M ⊙ . We identify the impact of both explosion energy and remnant mass on the final yields of the SLRs. Isotopes produced within the innermost regions of the star, such as 92 Nb and 97 Tc, are the most affected by the remnant mass, 92 Nb varying by five orders of magnitude. Isotopes synthesized primarily in explosive C-burning and explosive He-burning, such as 60 Fe, are most affected by explosion energies. 60 Fe increases by two orders of magnitude from the lowest to the highest explosion energy in the 15 M ⊙ model. The final yield of each examined SLR is used to compare to literature models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

On the synthesis of heavy nuclei in protomagnetar outflows and implications for ultra-high energy cosmic rays

ABSTRACT It has been suggested that strongly magnetized and rapidly rotating protoneutron stars (PNSs) may produce long duration gamma-ray bursts (GRBs) originating from stellar core collapse. We explore the steady-state properties and heavy element nucleosynthesis in neutrino-driven winds from such PNSs whose magnetic axis is generally misaligned with the axis of rotation. We consider a wide variety of central engine properties such as surface dipole field strength, initial rotation period, and magnetic obliquity to show that heavy element nuclei can be synthesized in the radially expanding wind. This process is facilitated provided the outflow is Poynting-flux dominated such that its low entropy and fast expansion time-scale enables heavy nuclei to form in a more efficient manner as compared to the equivalent thermal GRB outflows. We also examine the acceleration and survival of these heavy nuclei and show that they can reach sufficiently high energies ≳ 1020 eV within the same physical regions that are also responsible for powering gamma-ray emission, primarily through magnetic dissipation processes. Although these magnetized outflows generally fail to achieve the production of elements heavier than lanthanides for our explored electron fraction range 0.4–0.6, we show that they are more than capable of synthesizing nuclei near and beyond iron peak elements.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Numerical relativity simulations of the neutron star merger GW190425: microphysics and mass ratio effects

ABSTRACT GW190425 was the second gravitational wave (GW) signal compatible with a binary neutron star (BNS) merger detected by the Advanced LIGO and Advanced Virgo detectors. Since no electromagnetic counterpart was identified, whether the associated kilonova was too dim or the localization area too broad is still an open question. We simulate 28 BNS mergers with the chirp mass of GW190425 and mass ratio 1 ≤ q ≤ 1.67, using numerical-relativity simulations with finite-temperature, composition dependent equations of state (EOS) and neutrino radiation. The energy emitted in GWs is $\lesssim 0.083\mathrm{\, M_\odot }c^2$ with peak luminosity of 1.1–$2.4\times ~10^{58}/(1+q)^2\, {\rm {erg \, s^{-1}}}$. Dynamical ejecta and disc mass range between 5 × 10−6–10−3 and 10−5–$0.1 \mathrm{\, M_\odot }$, respectively. Asymmetric mergers, especially with stiff EOSs, unbind more matter and form heavier discs compared to equal mass binaries. The angular momentum of the disc is 8–$10\mathrm{\, M_\odot }~GM_{\rm {disc}}/c$ over three orders of magnitude in Mdisc. While the nucleosynthesis shows no peculiarity, the simulated kilonovae are relatively dim compared with GW170817. For distances compatible with GW190425, AB magnitudes are always dimmer than ∼20 mag for the B, r, and K bands, with brighter kilonovae associated to more asymmetric binaries and stiffer EOSs. We suggest that, even assuming a good coverage of GW190425’s sky location, the kilonova could hardly have been detected by present wide-field surveys and no firm constraints on the binary parameters or EOS can be argued from the lack of the detection.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Progenitors and explosion properties of supernova remnants hosting central compact objects: II. A global systematic study with a comparison to nucleosynthesis models

ABSTRACT Core-collapse explosions of massive stars leave behind neutron stars, with a known diversity that includes the ‘Central Compact Objects’ (CCOs). Typified by the neutron star discovered near the centre of the Cas A supernova remnant (SNR), CCOs have been observed to shine only in X-rays. To address their supernova progenitors, we perform a systematic study of SNRs that contain a CCO and display X-ray emission from their shock-heated ejecta. We make use of X-ray data primarily using the Chandra X-ray observatory, complemented with XMM–Newton. This study uses a systematic approach to the analysis of each SNR aimed at addressing the supernova progenitor as well as the explosion properties (energy and ambient density). After fitting for the ejecta abundances estimated from a spatially resolved spectroscopic study, we compare the data to six nucleosynthesis models making predictions on supernova ejecta yields in core-collapse explosions. We find that the explosion models commonly used by the astrophysics community do not match the ejecta yields for any of the SNRs, suggesting additional physics, for example multidimensional explosion models or updated progenitor structures, are required. Overall we find low-mass (≤25 solar masses) progenitors among the massive stars population and low-energy explosions (<1051 ergs). We discuss degeneracies in our model fitting, particularly how altering the explosion energy affects the estimate of the progenitor mass. Our systematic study highlights the need for improving on the theoretical models for nucleosynthesis predictions as well as for sensitive, high-resolution spectroscopy observations to be acquired with next generation X-ray missions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Secular outflows from 3D MHD hypermassive neutron star accretion disc systems

Magnetized hypermassive neutron stars (HMNSs) have been proposed as a way for neutron star mergers to produce high electron fraction, high-velocity ejecta, as required by kilonova models to explain the observed light curve of GW170817. The HMNS drives outflows through neutrino energy deposition and mechanical oscillations, and raises the electron fraction of outflows through neutrino interactions before collapsing to a black hole (BH). Here, in this study, we perform 3D numerical simulations of HMNS–torus systems in ideal magnetohydrodynamics, using a leakage/absorption scheme for neutrino transport, the nuclear APR equation of state, and Newtonian self-gravity, with a pseudo-Newtonian potential added after BH formation. Due to the uncertainty in the HMNS collapse time, we choose two different parametrized times to induce collapse. We also explore two initial magnetic field geometries in the torus, and evolve the systems until the outflows diminish significantly ($\sim\!\! 1\!\! - \!\!2\ \mathrm{s}$). We find bluer, faster outflows as compared to equivalent BH–torus systems, producing M ∼ 10 −3 M ⊙ of ejecta with Y e ≥ 0.25 and v ≥ 0.25c by the simulation end. Approximately half the outflows are launched in disc winds at times $t\lesssim 500 \ \mathrm{ms}$, with a broad distribution of electron fractions and velocities, depending on the initial condition. The remaining outflows are thermally driven, characterized by lower velocities and electron fractions. Nucleosynthesis with tracer particles shows patterns resembling solar abundances in all models. Although outflows from our simulations do not match those inferred from two-component modelling of the GW170817 kilonova, self-consistent multidimensional detailed kilonova models are required to determine whether our outflows can power the blue kilonova.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Proto-neutron star convection and the neutrino-driven wind: implications for the r-process

The neutrino-driven wind from proto-neutron stars is a proposed site for r-process nucleosynthesis, although most previous work has found that a wind heated only by neutrinos cannot produce the third r-process peak. However, several groups have noted that introducing a secondary heating source within the wind can change the hydrodynamic conditions sufficiently for a strong r-process to proceed. One possible secondary heating source is gravito-acoustic waves, generated by convection inside the proto-neutron star. As these waves propagate into the wind, they can both accelerate the wind and shock and deposit energy into the wind. Additionally, the acceleration of the wind by these waves can reduce the total number of neutrino captures and thereby reduce the final electron fraction of the wind. In neutron rich conditions, all of these effects can make conditions more favourable for r-process nucleosynthesis. Here, we present a systematic investigation of the impact of these convection-generated gravito-acoustic waves within the wind on potential nucleosynthesis. We find that wave effects in the wind can generate conditions favourable for a strong r-process, even when the energy flux in the waves is a factor of 10 –4 smaller than the total neutrino energy flux and the wind is marginally neutron rich. Nevertheless, this depends strongly on the radius at which the waves become non-linear and form shocks. We also find that both entropy production after shock formation and the acceleration of the wind due to stresses produced by the waves prior to shock formation impact the structure and nucleosynthesis of these winds.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Long-lived neutron-star remnants from asymmetric binary neutron star mergers: element formation, kilonova signals and gravitational waves

We present 3D general-relativistic neutrino-radiation hydrodynamics simulations of two asymmetric binary neutron star mergers producing long-lived neutron stars remnants and spanning a fraction of their cooling time scale. The mergers are characterized by significant tidal disruption with neutron rich material forming a massive disc around the remnant. The latter develops one-armed dynamics that is imprinted in the emitted kilo-Hertz gravitational waves. Angular momentum transport to the disc is initially driven by spiral-density waves and enhanced by turbulent viscosity and neutrino heating on longer timescales. The mass outflows are composed by neutron-rich dynamical ejecta of mass ∼10 −3 –10 −2 M ⊙ followed by a persistent spiral-wave/neutrino-driven wind of ≳ 10 −2 M ⊙ with material spanning a wide range of electron fractions, ∼0.1–0.55. Dynamical ejecta (winds) have fast velocity tails up to ∼0.8 (∼0.4) c. The outflows are further evolved to days timescale using 2D ray-by-ray radiation-hydrodynamics simulations that include an online nuclear network. We find complete r-process yields and identify the production of 56 Ni and the subsequent decay chain to 56 Co and 56 Fe. Synthetic kilonova light curves predict an extended (near-) infrared peak a few days postmerger originating from r-process in the neutron-rich/high-opacity ejecta and UV/optical peaks at a few hours (ten minutes) postmerger originating from weak r-process (free-neutron decay) in the faster ejecta components. Additionally, the fast tail of tidal origin generates kilonova afterglows potentially detectable in radio and X band on a few to ten years time scale. Quantitative effects originating from the tidal disruption merger dynamics are reflected in the multimessenger emissions.

abundances↗
DOE OSTI2025

2D end-to-end modelling of kilonovae from binary neutron star merger remnants

We investigate the kilonova emission resulting from outflows produced in a 3D general-relativistic magnetohydrodynamic (GRMHD) simulation of a hypermassive neutron star (HMNS) remnant. We map the outflows into the flash hydrodynamics code to model their expansion in axisymmetry, and study the effects of employing different r-process heating rates. Except for the highest heating rate prescription, we find no significant differences with respect to overall ejecta dynamics and morphology compared to the simulation without heating. Once homologous expansion is attained, typically after ~2s for these ejecta, we map the outflows to the sedona radiative transfer code and compute the spectral evolution of the kilonova and broad-band light curves in various Legacy Survey of Space and Time (LSST) bands. The kilonova properties depend on the remnant lifetime, with peak luminosities and peak time-scales increasing for longer lived remnants that produce more massive ejecta. For all models, there is a strong dependence of both the bolometric and broad-band light curves on the viewing angle. While the short-lived (12 ms) remnant produces higher luminosities when viewed from angles closer to the pole, longer lived remnants (240 ms and 2.5 s) are more luminous when viewed from angles closer to the equator. Our results highlight the importance of self-consistent, long-term modelling of merger ejecta, and taking viewing-angle dependence into account when interpreting observed kilonova light curves. We find that magnetized outflows from an HMNS – if it survives long enough – could explain blue kilonovae, such as the blue emission seen in AT2017gfo.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Optical follow-up observation for GW event S190510g using Subaru/Hyper Suprime-Cam

Abstract A gravitational wave event, S190510g, which was classified as a binary-neutron-star coalescence at the time of preliminary alert, was detected by LIGO/Virgo collaboration on 2019 May 10. At 1.7 hours after the issue of its preliminary alert, we started a target-of-opportunity imaging observation in the Y band to search for its optical counterpart using the Hyper Suprime-Cam (HSC) on the Subaru Telescope. The observation covers a 118.8 deg2 sky area corresponding to $11.6\%$ confidence in the localization skymap released in the preliminary alert and $1.2\%$ in the updated skymap. We divided the observed area into two fields based on the availability of HSC reference images. For the fields with the HSC reference images, we applied an image subtraction technique; for the fields without the HSC reference images, we sought individual HSC images by matching a catalog of observed objects with the PS1 catalog. The search depth is 22.28 mag in the former method and the limit of search depth is 21.3 mag in the latter method. Subsequently, we performed visual inspection and obtained 83 candidates using the former method and 50 candidates using the latter method. Since we only have the one-day photometric data, we evaluated the probability of candidates being located inside the 3D skymap by estimating their distances with photometry of associated extended objects. We found three candidates are likely located inside the 3D skymap and concluded they could be a counterpart of S190510g, while most of the 133 candidates were likely to be supernovae because the number density of candidates was consistent with the expected number of supernova detections. By comparing our observational depth with a light-curve model of such a kilonova reproducing AT2017gfo, we show that early deep observations with the Subaru/HSC can capture the rising phase of the blue component of a kilonova at the estimated distance of S190510g (∼230 Mpc).

79 ASTRONOMY AND ASTROPHYSICS↗