SEARCH · Engineering Papers
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.
Model Description • ALGE is a 3-dimensional hydrodynamic code to predict the movement of a dissolved tracer throughout a water body. • It has been recoded as an emergency response tool. • Can model large-scale bays and estuaries (Chesapeake Bay, Lake Michigan, etc.) • Model includes suspended sediment transport, dissolved and particulate chemical transport, radioactive decay, and dynamic salinity.
Molten salts play a crucial role in advancing Generation IV nuclear energy technology, with chloride-based salts like NaCl-UCl 3 garnering significant attention due to their distinctive properties. However, the corrosive nature of molten salts causes the dissolution of chromium (Cr), leading to the formation of CrCl 2 and CrCl 3 species in molten chlorides. Moreover, the radioactive decay of nuclear fuel gives rise to fission products, including Cs, Sr, and I chlorides. The first part of this report presents a comprehensive study utilizing ab initio molecular dynamics (AIMD) simulations to investigate the properties of eutectic NaCl-UCl 3 molten salt in the presence of corrosion products (CrCl 2 and CrCl 3 ) and fission product (CsCl). The study focuses on essential structural and thermophysical properties such as density, mixing energy, coordination numbers (CN), and Radial Distribution Functions (RDF) of the salts with varying compositions of corrosion products and fission product, covering a range from 0% to 13.5%. The results offer valuable insights into the behavior of corrosion and fission products in uranium-based molten salts, providing essential data that can be used as input to the MTDB-TC and MSTDB-TP property databases being developed by the NEAMS program. Due to their favorable characteristics such as low melting points, high boiling points, and low costs, MgCl 2 +NaCl+KCl (MNK) eutectic salts have recently attracted attention as high-temperature heat transfer fluids. The incorporation of LiCl into MNK salts can further reduce their melting points and increase their specific heat capacities, which is de- sirable for high-temperature heat transfer applications. The second part of this report presents the development and validation of a new shell-model potential for the MgCl 2 +NaCl+KCl+LiCl system, which captures the polarization of Cl anions. The extensive comparison with experimental data and AIMD simulations demonstrates the accuracy and reliability of the potential. Furthermore, using this potential, we elucidate the intricate network structure in MNK eutectic salts that is formed through polyhedron sharing. This research contributes to a fundamental understanding of the atomic structures and thermophysical properties of multi-component molten salts, which is crucial for future development of heat transfer fluids for applications in molten salt reactors.
LArTPCs provide sensitivity to GeV signals, such as accelerator neutrinos and part of the supernova neutrino spectrum. TinyTPC is a LArTPC test stand for R&D of LAr doping to expand the reach of LArTPCs down to the 1-10 MeV range, which would substantially enhance the flagship analyses of experiments like DUNE, while enabling low energy analyses. We aim to dope LAr with Xe and pho- tosensitive dopants to expand the LArTPC range by converting hard-to-detect scintillation light to efficiently detected ionization charge. A critical element of the data analysis in TinyTPC is calibrating the readout. This poster will cover the calibration of TinyTPC, a pixelated liquid argon detector, where we find the distance a muon travels through each pixel. We then calculate the energy loss of muons traveling through the detector from cosmic data. We can reconstruct the path of the particles through the TPC using a density-based clustering algorithm designed to sort straight cosmic muon tracks from low energy radioactive decay curled paths and electronic noise
Geant4 offers advanced event biasing techniques to significantly accelerate simulations involving rare events. Various biasing methods, such as leading particle selection, cross-section biasing, radioactive decay enhancement, and bremsstrahlung splitting, enable efficient event sampling, though they require careful handling. Additionally, Geant4 provides a Fast Simulation Interface, allowing the replacement of standard processes in specific region and for selected particles, enabling faster execution or external code integration. Applications of fast simulation include electromagnetic shower modeling in calorimeters, machine learning inference, and offloading tasks to specialized hardware like GPUs, making Geant4 a powerful tool for computationally demanding simulations.
Spectroscopic data for all nuclei with mass number A=47 have been evaluated and the corresponding level schemes from radioactive decay and reaction studies are presented. Highlights from this evaluation include first observation of gamma-ray transitions in 47 Ar, a new high precision measurement of 47 K β decay, and first observation of excited levels and gamma-ray transitions in 47 Mn.
In advanced reactor (AR) designs, a common feature is continuous chemical processing and circulation of the nuclear material. This work bridges a significant measurement gap in safeguarding reactors with circulating fuel or continuous refueling by leveraging and building on the isotope ratio method first developed by our team under an FY21 Advanced Reactors International Safeguards Engagement (ARISE) project (Uribe et al. 2021). In circulating fuel designs, the radioisotope inventory changes from traditional effects (e.g., radioactive decay, fission) but also includes material transport due to pressure and temperature gradients. Such designs may also require regular or continuous additions or removals during operation, which significantly increases the rate of inventory change compared to traditional pressurized water reactor (PWR)s. Thus, directly tracking the nuclear inventory is ineffective since the isotopes are continuously added and removed. The isotope ratio method instead focuses on detecting changes to the input and output flows of radioisotopes. Previous work showed that for well-chosen pairs of isotopes, the isotopic ratio provides a sensitive and lasting indicator of deviation from normal conditions (e.g., startup, shutdown, diversion). The isotope ratio method is a process monitoring method with potential for application in for forward-looking approaches to International Atomic Energy Agency (IAEA) safeguards. The original process monitoring method was developed for a specific case—the decay tank of a thorium-fueled molten salt breeder reactor. In this expanded work, we explored other types of reactors and processes with nonstationary (e.g., flowing) nuclear material, which are difficult to safeguard with traditional methods because of the transient nature of the systems. The goal of the present work is to generalize the isotope ratio method for use in processes with continuously flowing nuclear material. All continuous processes have an average time for isotopes to be replaced in the system. The isotope ratio method works by choosing radioisotopes with half-lives both above and below the processing time. The present work seeks to explore which isotopes are suitable for the method by simulating the nuclear inventory, radioisotope emissions, and detector responses for several classes of advanced reactors. While the method can in principle be applied to other processes (e.g., enrichment or reprocessing facilities), the present work limits scope to ARs with continuously flowing fuel. Section 2 details the mathematics supporting the isotope ratio method, and Section 3 introduces the representative ARs selected for this work. Section 4 discusses how each reactor was analyzed, and Section 5 showcases the results for each representative reactor. Finally, Section 6 provides concluding remarks and suggests pathways for further analysis.
Betavoltaics (BV) are long-life power sources that typically convert beta particle radiation into electricity. Largely, the radioactive decays within the source go unharvested by the device. This work seeks to augment the power generation of BV devices by integration of scintillating particles within the radiative getter to convert beta emission which would otherwise not leave the getter into usable light for power generation. Silica-covered barium fluoride scintillating particles were integrated into a tritiated water getter. Power generation was increased from 100s of nW to µW levels with the addition of 0.2 wt. % particles into the getter. Nanowatt-scale sensor platforms were demonstrated with the µW BV devices and the maximum possible lifetime of such platforms was estimated. As this technique enables higher power density BV devices from conventional Si semiconductors (compared to wider bandgap BVs), the further implementation may lower the barrier-to-deployment of these long-life power/sensor platforms.
This study measured radiolytic gas generation and pressure buildup in a sealed stainless-steel system containing alkaline Mo-99 solution with added sodium nitrate as a hydrogen suppressant. In each of two experiments, approximately 200 Ci of Mo-99 solution was transferred into a closed experimental vessel inside a hot cell, isolated, and monitored for pressure rise caused by gas generation during radioactive decay. After pressure buildup, headspace gas samples were collected and analyzed by mass spectrometry to determine hydrogen and oxygen concentrations. The purpose was to quantify the magnitude of pressure buildup in a closed system and to characterize the gas composition produced by radiolysis of the Mo-99 target solution under representative handling and storage conditions. The two experiments used similar total Mo-99 activity but differed in solution volume, headspace volume, and leak integrity. In the first experiment, 197.5 Ci of Mo-99 in 16.92 mL of solution was loaded into a vessel with a 31.1-mL headspace; a small leak was later identified, and the measured peak pressure of 46 psig was extrapolated to about 70 psig in the absence of leakage. The headspace gas from this experiment contained about 21.4% H2 and 63.4% O2, but the composition was influenced by preferential hydrogen loss through the leak. In the second experiment, 193 Ci of Mo-99 in 5.41 mL of solution with 0.46 g NaNO3 was loaded into a vessel with a larger 42.61 mL headspace, and no detectable leak was observed. This experiment reached a peak pressure of 38 psig, and the measured gas composition was approximately 28.2% H2 and 51.5% O2. Because the second experiment was leak-free, it is considered the more reliable indicator of the true pressure buildup and intrinsic radiolytic gas composition of the Mo-99 solution.
Deep space exploration requires specialized sources for both thermal and power applications. Radioactive decay heat of plutonium-238 (238Pu) provides these sources in the form of radioisotope thermoelectric generators (RTGs). The 238 Pu is produced via neutron capture reaction involving neptunium-237 ( 237 Np) target material. Continual optimization of 237 Np target materials and evaluation of potential alternative targets for production of 238 Pu RTGs are advantageous for meeting ongoing space power system resource requirements. Current production of 238 Pu for RTGs for the United States space program utilizes neptunium dioxide ( 237 NpO 2 ) targets; however, the use of neptunium mononitride ( 237 NpN) presents an opportunity to increase the mass of 237 Np per target compared to the dioxide form, as well as increase the thermal conductivity of the target. To assess the viability of a 237 NpN target material, the material chemistry must be thoroughly evaluated, including synthesis methods and dissolution and reprocessing schemes. This review presents a summary of synthesis pathways for 237 NpN based on published literature on actinide mononitrides. Specific literature on 237 NpN is limited, necessitating evaluation of other actinide systems to gather parallels. This suggests a need for additional experimental studies on 237 NpN. A particular limitation in the existing literature is a lack of information on the differences in material characteristics, such as morphology, particle size, and trace chemical impurities, as a function of synthesis method. These parameters may affect subsequent reactor performance or dissolution of irradiated targets. The evaluation of existing literature is presented with a focus on the efficacy of 237 NpN targets for 238 Pu production.
The merger of neutron star binaries is believed to eject a wide range of heavy elements into the universe. By observing the emission from this ejecta, scientists can probe the ejecta properties (mass, velocity, and composition distributions). The emission (a.k.a. kilonova) is powered by the radioactive decay of the heavy isotopes produced in the merger and this emission is reprocessed by atomic opacities to optical and infrared wavelengths. Understanding the ejecta properties requires calculating the dependence of this emission on these opacities. The strong lines in the optical and infrared in lanthanide opacities have been shown to significantly alter the light curves and spectra in these wavelength bands, arguing that the emission in these wavelengths can probe the composition of this ejecta. In this work, we study variations in the kilonova emission by varying individual lanthanide (and the actinide uranium) concentrations in the ejecta. The broad forest of lanthanide lines makes it difficult to determine the exact fraction of individual lanthanides. Nd is an exception. Its opacities above 1 μm are higher than other lanthanides and observations of kilonovae can potentially probe increased abundances of Nd. Similarly, at early times when the ejecta is still hot (first day), the U opacity is strong in the 0.2–1 μm wavelength range and kilonova observations may also be able to constrain these abundances.
The recent detection of the gravitational-wave event GW170817, produced by the coalescence of two neutron stars, and of its optical–infrared counterpart, powered by the radioactive decay of r-process elements, has opened a new window onto gamma-ray astronomy: the direct detection of photons coming from such decays. In this work, we calculate the contribution of kilonovae to the diffuse gamma-ray background in the MeV range, using recent results on the spectra of the gamma-rays emitted in individual events, and we compare it with that from other sources. We find that the contribution from kilonovae is not dominant in such an energy range, but within current uncertainties, and its addition to other sources might help to fit the observational data.
Here, we present optical, radio, and X-ray observations of SN 2020bvc (=ASASSN20bs; ZTF20aalxlis), a nearby (z = 0.0252; d = 114 Mpc) broad-lined (BL) Type Ic supernova (SN) and the first double-peaked Ic-BL discovered without a gamma-ray burst (GRB) trigger. Our observations show that SN 2020bvc shares several properties in common with the Ic-BL SN 2006aj, which was associated with the low-luminosity gamma-ray burst (LLGRB) 060218. First, the 10 GHz radio luminosity (L radio ≈ 10 37 erg s -1 ) is brighter than ordinary core-collapse SNe but fainter than LLGRB-SNe such as SN 1998bw (associated with LLGRB 980425). We model our VLA observations (spanning 13–43 d) as synchrotron emission from a mildly relativistic (v ≳ 0.3c) forward shock. Second, with Swift and Chandra we detect X-ray emission (L X ≈ 10 41 erg s -1 ) that is not naturally explained as inverse Compton emission or as part of the same synchrotron spectrum as the radio emission. Third, high-cadence (6×/night -1 ) data from the Zwicky Transient Facility (ZTF) shows a double-peaked optical light curve, the first peak from shock-cooling of extended low-mass material (mass M e < 10 -2 M ⊙ at radius R e > 10 12 cm) and the second peak from the radioactive decay of 56 Ni. SN 2020bvc is the first double-peaked Ic-BL SN discovered without a GRB trigger, so it is noteworthy that it shows X-ray and radio emission similar to LLGRB-SNe. For four of the five other nearby (z . 0.05) Ic-BL SNe with ZTF high-cadence data, we rule out a first peak like that seen in SN 2006aj and SN 2020bvc, i.e. that lasts ≈ 1 d and reaches a peak luminosity M ≈ -18. X-ray and radio follow-up observations of Ic-BL SNe with well-sampled early optical light curves will establish whether double-peaked optical light curves are indeed predictive of LLGRB-like X-ray and radio emission.
The Type Ia supernova (SN Ia) LSQ14fmg exhibits exaggerated properties that may help to reveal the origin of the "super-Chandrasekhar"(or 03fg-like) group. The optical spectrum is typical of a 03fg-like SN Ia, but the light curves are unlike those of any SNe Ia observed. The light curves of LSQ14fmg rise extremely slowly. At -23 rest-frame days relative to B-band maximum, LSQ14fmg is already brighter than M V = -19 mag before host extinction correction. Here, the observed color curves show a flat evolution from the earliest observation to approximately 1 week after maximum. The near-infrared light curves peak brighter than -20.5 mag in the J and H bands, far more luminous than any 03fg-like SNe Ia with near-infrared observations. At 1 month past maximum, the optical light curves decline rapidly. The early, slow rise and flat color evolution are interpreted to result from an additional excess flux from a power source other than the radioactive decay of the synthesized 56 Ni. The excess flux matches the interaction with a typical superwind of an asymptotic giant branch (AGB) star in density structure, mass-loss rate, and duration. The rapid decline starting at around 1 month past B-band maximum may be an indication of rapid cooling by active carbon monoxide (CO) formation, which requires a low-temperature and high-density environment. These peculiarities point to an AGB progenitor near the end of its evolution and the core degenerate scenario as the likely explosion mechanism for LSQ14fmg.
Abstract Kilonovae are ultraviolet, optical, and infrared transients powered by the radioactive decay of heavy elements following a neutron star merger. Joint observations of kilonovae and gravitational waves can offer key constraints on the source of Galactic r -process enrichment, among other astrophysical topics. However, robust constraints on heavy element production require rapid kilonova detection (within ∼1 day of merger) as well as multiwavelength observations across multiple epochs. In this study, we quantify the ability of 13 wide-field-of-view instruments to detect kilonovae, leveraging a large grid of over 900 radiative transfer simulations with 54 viewing angles per simulation. We consider both current and upcoming instruments, collectively spanning the full kilonova spectrum. The Roman Space Telescope has the highest redshift reach of any instrument in the study, observing kilonovae out to z ∼ 1 within the first day post-merger. We demonstrate that BlackGEM, DECam, GOTO, the Vera C. Rubin Observatory’s LSST, ULTRASAT, VISTA, and WINTER can observe some kilonovae out to z ∼ 0.1 (∼475 Mpc), while DDOTI, MeerLICHT, PRIME, Swift/UVOT, and ZTF are confined to more nearby observations. Furthermore, we provide a framework to infer kilonova ejecta properties following nondetections and explore variation in detectability with these ejecta parameters.
SN 2018aoz is a Type Ia SN with a B-band plateau and excess emission in infant-phase light curves ≲1 day after the first light, evidencing an over-density of surface iron-peak elements as shown in our previous study. Here, we advance the constraints on the nature and origin of SN 2018aoz based on its evolution until the nebular phase. Near-peak spectroscopic features show that the SN is intermediate between two subtypes of normal Type Ia: core normal and broad line. The excess emission may be attributable to the radioactive decay of surface iron-peak elements as well as the interaction of ejecta with either the binary companion or a small torus of circumstellar material. Nebular-phase limits on Hα and He ι favor a white dwarf companion, consistent with the small companion size constrained by the low early SN luminosity, while the absence of [O ι] and He ι disfavors a violent merger of the progenitor. Of the two main explosion mechanisms proposed to explain the distribution of surface iron-peak elements in SN 2018aoz, the asymmetric Chandrasekhar-mass explosion is less consistent with the progenitor constraints and the observed blueshifts of nebular-phase [Fe ιι] and [Ni ιι]. The helium-shell double-detonation explosion is compatible with the observed lack of C spectral features, but current 1D models are incompatible with the infant-phase excess emission, $B_{max}$--$V_{max}$ color, and weak strength of nebular-phase [Ca ιι]. Although the explosion processes of SN 2018aoz still need to be more precisely understood, the same processes could produce a significant fraction of Type Ia SNe that appear to be normal after ~1 day.
GRB 171205A is a low-luminosity, long-duration gamma-ray burst (GRB) associated with SN 2017iuk, a broad-line type Ic supernova (SN). It is consistent with having been formed in the core collapse of a widely separated binary, which we have called the binary-driven hypernova of type III. The core collapse of the CO star forms a newborn NS (νNS) and the SN explosion. Fallback accretion transfers mass and angular momentum to the νNS, here assumed to be born non-rotating. The accretion energy injected into the expanding stellar layers powers the prompt emission. The multiwavelength power-law afterglow is explained by the synchrotron radiation of electrons in the SN ejecta, powered by energy injected by the spinning νNS. We calculate the amount of mass and angular momentum gained by the νNS, as well as the νNS rotational evolution. The νNS spins up to a period of 47 ms, then releases its rotational energy powering the synchrotron emission of the afterglow. The paucity of the νNS spin explains the low-luminosity characteristic and that the optical emission of the SN from the nickel radioactive decay outshines the optical emission from the synchrotron radiation. From the νNS evolution, we infer that the SN explosion had to occur at most 7.36 h before the GRB trigger. Therefore, for the first time, the analysis of the GRB data leads to the time of occurrence of the CO core collapse leading to the SN explosion and the electromagnetic emission of the GRB event.
Abstract Type IIP supernovae (SNe IIP) are thought to originate from the explosion of massive stars >10M ⊙ . Their luminosity is primarily powered by the explosion energy and the radioactive decay energy of 56 Co, with the photosphere location regulated by hydrogen recombination. However, the physical connections between SNe IIP and their progenitor stars remain unclear. This paper presents a comprehensive study of SNe IIP and their progenitor stars by using the one-dimensional stellar evolution code,MESA. Our model grids consider the effects of stellar metallicity, mass, and rotation in the evolution of massive stars, as well as the explosion energy and 56 Ni production in modeling supernovae. To elucidate the observed SNe IIP and their origins, we compare their light curves (LCs) with our models. Furthermore, we investigate the impact of stellar parameters on LCs by considering stellar mass, metallicity, rotation, explosion energy, and 56 Ni production. We find that more massive stars exhibit longer plateaus due to increased photon diffusion time caused by massive ejecta. Higher metallicity leads to increased opacity and mass loss of progenitor stars. Rapid rotation affects internal stellar structures, enhancing convective mixing and mass loss, potentially affecting the plateau’s brightness and duration. Higher explosion energy results in brighter but shorter plateaus due to faster-moving ejecta. 56 Ni mass affects late-time luminosity and plateau duration, with larger masses leading to slower declines.
The precise atomic structure and therefore the wavelength-dependent opacities of lanthanides are highly uncertain. This uncertainty introduces systematic errors in modeling transients like kilonovae and estimating key properties such as mass, characteristic velocity, and heavy metal content. Here, we quantify how atomic data from across the literature as well as choices of thermalization efficiency of r-process radioactive decay heating impact the light curve and spectra of kilonovae. Specifically, we analyze the spectra of a grid of models produced by the radiative transfer code Sedona that span the expected range of kilonova properties to identify regions with the highest systematic uncertainty. Our findings indicate that differences in atomic data have a substantial impact on estimates of lanthanide mass fraction, spanning approximately 1 order of magnitude for lanthanide-rich ejecta, and demonstrate the difficulty in precisely measuring the lanthanide fraction in lanthanide-poor ejecta. Mass estimates vary typically by 25%–40% for differing atomic data. Similarly, the choice of thermalization efficiency can affect mass estimates by 20%–50%. Observational properties such as color and decay rate are highly model dependent. Velocity estimation, when fitting solely based on the light curve, can have a typical error of ∼100%. Atomic data of light r-process elements can strongly affect blue emission. Even for well-observed events like GW170817, the total lanthanide production estimated using different atomic data sets can vary by a factor of ∼6.