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.
In 2020, the HAWC Collaboration announced the discovery of several gamma-ray sources emitting above 56 TeV. This was the highest-energy astrophysical source catalog ever published.
Underground storage tanks at the Hanford Site, in southeastern Washington State, hold radioactive waste generated from four decades of plutonium production. The 149 single-shell tanks and the 28 double-shell tanks have all exceeded their initial design life of approximately 25 years. At least 67 tanks are assumed to have leaked in the past, resulting in radioactive releases into the vadose zone. Gamma ray logging within dry monitoring wells is currently the primary method for tracking the migration of leaked tank waste through the vadose zone. While this approach provides an accurate assessment of radioactive contamination, that information is only provided near (within ~1m) the borehole, leaving most of the vadose zone unmonitored, particularly the important region directly beneath the tank. This report describes a numerical study that investigates the feasibility and performance of time-lapse 3D electrical resistivity tomography (ERT) for long-term monitoring of a hypothetical tank waste location and migration through the vadose zone. ERT is a method of remotely imaging the bulk electrical conductivity (EC) of the subsurface, which is significantly impacted by the presence of conductive solid and liquid tank waste. The release of liquid tank wastes increases subsurface fluid conductivity and saturation over time, creating a target to use time-lapse ERT for long-term monitoring. Although the presence of metallic infrastructure can cause ERT interference, recent advancements in ERT data processing enable the deleterious effects of buried metallic infrastructure (e.g. pipes, wellbore casings, tanks) to be removed to better determine the liquid tank waste migration over time. Three hypothetical realistic scenarios were simulated in the ERT evaluation. The first two scenarios assume the same leak amount and rate (i.e., between 1/1/1951 and 12/31/1951 at the rate of 347 m 3 per year) but different leaky tanks. Scenario 1 assumes leaks under tank B-102, which is located on the edge of the B-tank farm and surrounded by a few metallic infrastructure including cased pipes/wells/tanks. Scenario 2 assumes leaks under tank B-108, which is located near the center of the B-tank farm and surrounded by larger amount of metallic infrastructure than B-102. Scenario 3 assumes the same metallic infrastructure as B-102, with a more recent contaminant leak that was simulated to have occurred between 1/1/2018 and 12/31/2023 at a rate of 1.89 m 3 per year. The leak time in Scenarios 1 and 2 corresponds to a historical overfill event in 1951 and Scenario 3 corresponds to a recent found tank leak in 2019. In each scenario, a “true” bulk EC model vs. time reflecting contaminant migration was generated. ERT data was simulated from these “true” bulk EC models and a time-lapse ERT inversion produced “imaged” bulk EC vs. time. Three electrode configurations in two, four and eight boreholes surrounding the leak tank were used in the ERT simulations in each scenario. These borehole configurations were considered logistically feasible and cost-effective for monitoring. The hypothetical ERT boreholes are assumed to have non-metallic casing. By comparing the “imaged” bulk EC with the “true” bulk EC, it was demonstrated that the three configurations of wells used (two, four, and eight wells) were able to successfully monitor the migration of tank leaks through the vadose zone, with bulk EC resolution increasing with the number of down borehole ERT arrays for all scenarios. Therefore, the use of eight boreholes to perform ERT monitoring beneath the tanks provided the best spatiotemporal information.
The Taylor series (differential operator) perturbation method, as implemented in the PERT capability in MCNP6 (Ref. 1), can be used for first- and second-order sensitivity analyses in fixed-source problems. Unfortunately, the PERT capability does not yet work with point-detector (F5) or pulse-height (F8) tallies. Also, the second-order Taylor series term lacks the cross terms that would allow the calculation of mixed sensitivities. Finally, the PERT capability does not compute any sensitivities involving the source in a fixed-source problem. When neutrons or gamma rays are emitted due to radioactive decay, the source rate density depends intrinsically on the nuclide densities in the source material. When the PERT capability in MCNP6 is used to compute sensitivities of a response to nuclide and material densities, those sensitivities do not include the effect of the source.
Detecting diversion of spent fuel elements in dry storage casks is challenging due to the thick shielding used in cask construction. Measurements on top of the cask to map the underlying arrangement of the fuel elements and looking for anomalous changes over time has proven difficult to achieve using gamma rays due to the high scattering and attenuation from the thick steel structure, weakening information on the present or absence of fuel bundles. Simulations and laboratory experiments suggest that the high-energy neutron flux (>200 keV) measured directly above each fuel bundle is sufficient to produce a position map that enables detection of the present or absence of fuel bundles, and therefore diversion of a spent fuel bundle. A single-detector spent-fuel monitoring technique based on this principle was development at the Lawrence Livermore National Laboratory (LLNL). The INL Cask Farm in the INTEC technical area at Idaho National Laboratory (INL) offers the capability to test this technique on an MC-10 storage cask which has a distribution of full and empty fuel positions. An experimental test plan for the single-detector verification system was developed in consultation with INL personnel to be completed in FY2021. Due to travel advisories related to COVID-19, the experimental test plan was adapted to enable INL personnel to carry out the measurements in consultation with LLNL personnel following shipment of the LLNL system to INL. Field test measurements of the single detector verification system were successfully carried out at the INL cask farm on September 7-9, 2021. This document summarizes results of the field test.
Moving from traditional fixed fuel nuclear reactor systems to a mobile, dynamic fuel system that has dissolved special nuclear material in a molten salt is a paradigm shift in several respects. One major consideration is how to develop effective nuclear material controls and accounting for these novel reactor systems. The Molten Salt Reactor Experiment is one example of a critical molten salt system and provided a significant reference library based on the documented effort. But it was low thermal power that was not intended to reflect a commercial scale electricity production design. Therefore, to facilitate and assist vendors with domestic licensing considerations, research is underway to identify methods that could be used for domestic safeguards approaches for these novel reactor systems. This research presents the results of a signature analysis of data generated from three simulated scenarios using the molten salt demonstration reactor model defined in the Transient Simulation Framework of Reconfigurable Modules. Each scenario provides 1 hour isotope inventories over a 180 day period. The scenarios investigated provide test cases to examine if direct gamma-ray spectroscopy of the fuel can be used to identify changes comparing a base case (no reactivity control and fixed fission contribution) to an insertion of reactivity (10 pcm no change in fission composition) and a change in fission composition. The analysis demonstrates that monitoring the total count rate in a highly collimated high-resolution photon energy spectrum is sensitive to perturbations imposed into the reactor model. The total photon count rate changes ≈2% for the fission composition change and ≈4.5% for the reactivity insertion compared to the base case. However, both scenarios show an increase in the total photon count rate. The total photon count rate can be used to identify changes due to power (number of fissions) but not due to a change in the material undergoing fission. To distinguish between the two cases of increased power, the photon spectrum would require an intensive analysis technique. A photon peak count rate ratio analysis could be used to identify changes in the fissile material fission generation in the core through identification of a static peak that shows little variation to the source of fission and a highly varying peak. The photon peak strength will ultimately be determined by the isotope’s fission yield. A preliminary analysis investigating the coupling of the isotope’s fission yield and its concentration in the fuel salt derived from the modeling has been performed. A ratio analysis of the photon peak count rates of 140 La to 99 Mo demonstrated that the reactivity insertion creates a distinct difference in the ratio compared to the fission composition change scenario.
The evolution of philosophy and computations in the International Data Center (IDC) related to aerosol samples have had profound impacts on the number of recorded detections in the network since routine operations began in 2000. Key decisions from policymakers have been the list of triggering radionuclides, the scheme for categorizing these into interest levels 1-5, and an algorithm for determining when an anthropogenic isotope is seen so often that it is no longer interesting, known as the Exponential Weighted Moving Average (EWMA). These are described in the Operations Manual of the IDC. Key parameters that are controlled by the IDC but for which the IDC receives occasional input from policymakers include the constants in EWMA and the peak significance threshold for individual gamma rays, the latter of which directly leads to determination of the presence or absence of a radionuclide in a sample. There are also changes in computations which the IDC makes and informs policy makers about, such as changes in how background is computed, which could also affect the ease of detecting a peak – real or false. Rather than focus on the quantitative changes due to computation changes, this work records some thinking on how isotopes and peak significance levels were chosen, and the resulting detections seen over 18 years during the buildup of the International Monitoring System (IMS). These detections are considered on a global scale to try to determine the relative impact on monitoring, and in some cases, the nature of their existence. Repeated detections of 131I and 133I are the most troublesome, but they are not so frequent to be a major problem for the Verification Regime. These detections could probably be handled adequately using scientific methods currently under development for xenon backgrounds. It is also somewhat problematic that top-level analysis of aerosol backgrounds has not been reported previous to this. The steep increase in the rate of detections after 2016 are a concern, either in the actual backgrounds or from changes in the calculations methods used to generate the Reviewed Radionuclide Report (RRR.) Final conclusions of the authors are that the computational stability of the RRR is very important. With computational stability, changes can be usefully analyzed as being due to changes in radioactivity in Earth’s atmosphere This report is a distillation into text of a talk given in the Radionuclide Experts Group (RNEG) in Vienna during Working Group B (WGB) in February of 2019. This report does not directly contain any IDC data, only summaries by year, or by isotope, or by location. No specific IDC detection by time, location, or isotope is included.
Magnetic fields can be the dominant component of astrophysical plasmas, so that the magnetic energy density might exceed even the rest-mass energy density of matter. In this extreme (and largely unexplored) regime the magnetic field controls the overall plasma evolution, dissipation, and acceleration of non-thermal particles. This plasma regime, applicable to a variety of astrophysical sources --- magnetars, pulsars and pulsar wind nebulae (PWNe), jets of Active Galactic Nuclei (AGNs) and Gamma-Ray Bursters (GRBs) --- is dramatically different from laboratory plasmas, the magnetospheres of planets, and the interplanetary plasma. Relativistic astrophysical sources then provide an unique opportunity to study the fundamental plasma physics of magnetically-dominated plasmas; a novel and quickly evolving field of theoretical research which, by investigating energy conversion and particle energization processes in plasmas, is of vital importance to the Fusion Energy Sciences DoE program.
The nodal diffusion method is one of the most widely used approaches in modern reactor analysis. In the nodal diffusion method, a coarse multi-group set of “homogenized” parameters is constructed such that the complex geometry of a reactor core along with the energy dependence of neutron and gamma ray cross sections in a nuclear reactor are conserved in the simpler geometry. The homogenization is typically done on a fuel assembly level as is the case in the DIF3D code developed at Argonne National Laboratory. The nodal methodology is used primarily to predict fuel cycle behavior of nuclear systems of which there is a substantial amount of validation in the literature. Another use of the nodal method is to obtain reactivity coefficients and kinetics parameters for use in a safety analysis of a given nuclear reactor. While there are many ways to obtain reactivity worth and kinetics parameters, the work presented in this manuscript is unique as it provides the user with the ability to compute reactivity worths, kinetics parameters, and cross section sensitivities with a Cartesian and hexagonal geometry based transport code.
The nodal diffusion method is one of the most widely used approaches in modern reactor analysis. In the nodal diffusion method, a coarse multi-group set of “homogenized” parameters is constructed such that the complex geometry of a reactor core along with the energy dependence of neutron and gamma ray cross sections in a nuclear reactor are conserved in the simpler geometry. The homogenization is typically done on a fuel assembly level as is the case in the DIF3D code developed at Argonne National Laboratory. The nodal methodology is used primarily to predict fuel cycle behavior of nuclear systems of which there is a substantial amount of validation in the literature. Another use of the nodal method is to obtain reactivity coefficients and kinetics parameters for use in a safety analysis of a given nuclear reactor. While there are many ways to obtain reactivity worth and kinetics parameters, the work presented in this manuscript is unique as it provides the user with the ability to compute reactivity worths, kinetics parameters, and cross section sensitivities with a Cartesian and hexagonal geometry based transport code.
In this lecture, we cover deterministic particle transport. We focus on neutral particles (neutrons, gamma-rays, thermal photons). Our goal is to predictively and accurately model how free particles interact with their environments and properties of systems involving those particles.
A recent journal article by J.A. Gomez et al. measured gamma-ray die-off curves of three subcritical static highly-enriched uranium (HEU) assemblies driven by an external neutron source with a new detector system. This new detector system was developed for being used in dynamically driven subcritical assemblies as part of the Neutron Diagnosed Subcritical Experiments (NDSE) program. Simulations of these die-off curves with various nuclear data and comparison to experimental data indicated that a decrease of the ENDF/B-VII.1 235 U(n,inl) cross section by a factor 0.8 and 0.85 for ENDF/B-VIII.0 would lead to better predictions of experimental data. Here, we test the proposed changes with another type of measurement response, namely neutron-leakage spectra emitted in LLNL pulsed sphere measurements. These spheres were pulsed by 14-MeV neutrons produced via the D+T reaction in their center. The proposed changes in nuclear data have a distinctly smaller impact on predicting pulsed-sphere neutron-leakage spectra than for the die-off curves; they lead to a worsened prediction of the inelastic valley of LLNL pulsed-sphere neutron spectra indicating that the proposed change could constitute a compensating error. Changes in the inelastic angular distributions along with the cross section might be worthwhile to study
We performed a branching ratio measurement of the most intense gamma rays in the decay of 111Ag. The sample was harvested at Argonne National Laboratory and measured at Texas A&M with the same technique as we used in our previous measurements. We made two samples and were able to measure both of them. The data is currently under analysis.
The central importance of the nuclear data curated by the US Nuclear Data Program (USNDP) for clean energy generation, national security, nonproliferation, medical applications, and space exploration as well as basic science was described in a prior report issued by the DOE/NSF Nuclear Science Advisory Committee subcommittee on Nuclear Data (NSAC-ND) in September 2022. In this report, we present a set of fourteen (14) recommendations that will enhance and advance DOE-NP's stewardship of nuclear data. The first three recommendations focus on the existing core USNDP capabilities, namely: 1) Support the nuclear structure evaluation workforce to improve the currency, consistency, and accessibility of the Evaluated Nuclear Structure Data File (ENSDF); 2) Enhance nuclear reaction evaluation within the USNDP in support of the Evaluated Nuclear Data File (ENDF) through expansion of the workforce and integration of high-performance computing, automation, and machine learning and; 3) Continue atomic mass evaluation in support AME and NUBASE databases. This is followed by eight (8) recommendations representing new cross-cutting initiatives involving both measurement and evaluation to address outstanding nuclear data needs. These new initiatives require a highly trained, diverse workforce that includes personnel with expertise from both inside and outside the nuclear physics community from which evaluators have traditionally been recruited. As such, many of these initiatives are accomplished via a Topical Nuclear Data Collaborations (TNDC). A TNDC is made up of domestic and international stakeholders, subject matter and nuclear data experts, and nuclear data evaluators and features a workforce development plan to ensure that nuclear data evaluators maintain currency in the relevant applications and are seen as equity partners in the endeavor. These include: 1) Establish a coordinated effort to improve evaluation and modeling in nuclear astrophysics for stellar dynamics, multi-messenger astronomy and nucleosynthesis; 2) Initiate a TNDC to develop and maintain nuclear structure evaluation beyond discrete states, including nuclear level densities, photon strength functions and photonuclear data for improved reaction modeling, and exploring nuclear structure at finite temperature; 3) Create a TNDC to perform correlated fission data evaluation, including cross sections, fragment yields, v(A), v(E n ) for nuclear energy, national security, nonproliferation and basic science; 4) From a panel of subject matter experts to establish and annually update a roster of key decay data to nurture its accelerated dissemination including both measurement and evaluation for targeted high-value nuclides for national security, nonproliferation and medical applications; 5) Comprehensive, consistent neutron-induced structure and reaction data for nuclear energy, national security, nonproliferation and planetary nuclear spectroscopy; 6) Charged-particle stopping powers for detector design, space effects and ion beam therapy; 7) High-energy reactions for space exploration and medical nuclide production, and; 8) The creation of an infrastructure for open data and data preservation for use by the entire nuclear physics community. All told, these initiatives require approximately $6.5M increase in NP support of the USNDP in fiscal year 2023 dollars and would require at least 3-5 years to carry out due to the length of time needed to recruit and train new nuclear data researchers. This relatively modest investment would help ensure that the fruits of the nuclear data research carried out by DOE-NP and its collaborators would be brought to bear to address some of the most important needs of our nation and the world. To ensure effective execution of this plan, we present an overview of recruitment, training, and retention goals for the USNDP, the centerpiece of which is a mutually agreed upon code of conduct. Finally, we identify the facility and instrumentation needed to perform the recommended experimental activities. This includes a short review of target fabrication capabilities, reactors, neutron beam, light- and heavy-stable ion, gamma-ray, high-energy and radioactive ion beam facilities. Lastly, a more complete appendix of experimental facilities previously compiled is included with new input provided for 6 facilities.
Density and electron fraction of neutron star stuff leftover after a black hole eats the neutron star. Events like this one are believed to power gamma ray bursts, some of the most energetic events in the universe, and be the source of heavy elements like gold and platinum.
Optomechanical strain sensing provides attractive opportunities for novel particle detection schemes, as well as studying stress-induced (i.e. non-radiogenic) phonon bursts, which have been demonstrated to limit the coherence times of superconducting qubits [1,2] and are a suspected culprit in the low energy excesses observed by many dark matter direct detection experiments [3,4]. We are investigating SiN microring optical resonator strain sensors, developed at Purdue University [5], for applications in fundamental particle sensing and QIS. These sensors can be embedded in the substrate upon which superconducting qubits are patterned, providing a handle to distinguish decoherence events of radiogenic origin from those due to crystal stress. In a similar way, these sensors can be operated in conjunction with superconducting detectors (e.g., MKIDs, TES) to enable multi-channel readout of particle interactions in the device substrate or serve as anticoincidence detectors, which may be r equired to identify low-energy interactions from dark matter particles down to the fermionic thermal relic mass limit of a few keV. Such sensors can also be used to directly observe resonant scattering processes of gamma rays (and perhaps neutrinos) where no detectable quanta are produced in the target, via the microscopic strain induced by the momentum transfer to the (fixed-in-place) crystal lattice as a whole. These strain sensors have so far found application in photonics and communications, but have yet to be adopted for HEP uses, where they can provide unique capabilities in the search for dark matter, observing neutrino interactions, and improving the coherence times of superconducting qubits.
In this lecture, we cover deterministic particle transport. We focus on neutral particles (neutrons, gamma-rays, and thermal photons). Our goal is to predictively and accurately model how free particles interact with their environments and properties of systems involving those particles.
This report discussing the March 2022 irradiation of U233 on the Godiva critical assembly using burst mode and subsequent analysis of the fission products starting at 65 minutes post irradiation, and ending 12 days later. The Short-Lived Fission Product Yield (SLFPY) project is one of many efforts under the Nuclear Physics Multi-Lab and related Venture Project. The objective of the SLFPY project is to provide improved measurements of fission product yields for select nuclides. The SLFPY project consists of measurement and analysis campaigns. To produce fission data, samples of actinide material are irradiated in a high-flux neutron environment. These samples are quickly placed in front of a detection system where both gamma singles and coincidence spectra are collected as a function of time. Gamma spectroscopy is performed on the irradiated material on site near the irradiation facility, starting as soon as feasible after the irradiation and continuing for 7 days or longer. Following the data collection, time-dependent fission yield analyses are performed using gamma-ray spectroscopy. This document summarizes the measurement campaign performed March 21 - April 4, 2022 where sub-gram quantities of 233 U nitrate (UO 2 (NO 3 ) 2 ) were irradiated in the Godiva reactor. The Godiva burst #2076 conducted under experiment IER 504 had a temperature change of +156 degrees Celsius, lasted 13.32 microseconds, and had a reactivity of $\$$1.089. The targets consisted of uranium oxide solution deposited in a quartz tube and packed with quartz wool. The flame-sealed quartz ampoules were used to contain the uranium samples for the irradiation and count activities. On March 22, 2022, at 10:27:00 a.m. the three uranium ampoules, along with a witness foil pack, were irradiated in the Godiva critical assembly using burst mode operations.
The first observation of gravitational waves (GW) from a binary neutron star collision, known as GW170817, was detected by Laser Interferometer Gravitational-Wave Observatory (LIGO) and Virgo GW observatories. This event marked a profound discovery of the kilonova, which occurs when two neutron stars collide, resulting in an energetically explosive release of gamma-ray bursts (GRBs) and emitting bright electromagnetic radiation. This project focuses on assessing how well kilonovae will be detected in the Legacy Survey of Space and Time (LSST) at the Vera C. Rubin Observatory. By using Kasen models with varying ejecta mass, ejecta velocity, and lanthanide fraction, we analyze these 329 models and create kilonova light curves to assess the number of visible kilonovae for each band, based on the minimum magnitude at various redshifts.