Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Decay Curve”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Inverse energy transfer in decaying, three-dimensional, non-helical magnetic turbulence due to magnetic reconnection

ABSTRACT It has been recently shown numerically that there exists an inverse transfer of magnetic energy in decaying, non-helical, magnetically dominated, magnetohydrodynamic turbulence in three dimensions (3D). We suggest that magnetic reconnection is the underlying physical mechanism responsible for this inverse transfer. In the two-dimensional (2D) case, the inverse transfer is easily inferred to be due to smaller magnetic islands merging to form larger ones via reconnection. We find that the scaling behaviour is similar between the 2D and 3D cases, i.e. the magnetic energy evolves as t−1, and the magnetic power spectrum follows a slope of k−2. We show that on normalizing time by the magnetic reconnection time-scale, the evolution curves of the magnetic field in systems with different Lundquist numbers collapse on to one another. Furthermore, transfer function plots show signatures of magnetic reconnection driving the inverse transfer. We also discuss the conserved quantities in the system and show that the behaviour of these quantities is similar between the 2D and 3D simulations, thus making the case that the dynamics in 3D could be approximately explained by what we understand in 2D. Lastly, we also conduct simulations where the magnetic field is subdominant to the flow. Here, too, we find an inverse transfer of magnetic energy in 3D. In these simulations, the magnetic energy evolves as t−1.4 and, interestingly, a dynamo effect is observed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

A Stochastic Approach to Reconstruct Gamma-Ray-burst Light Curves

Gamma-ray bursts (GRBs), as they are observed at high redshift ( z = 9.4), are vital to cosmological studies and investigating Population III stars. To tackle these studies, we need correlations among relevant GRB variables with the requirement of small uncertainties on their variables. Thus, we must have good coverage of GRB light curves (LCs). However, gaps in the LC hinder the precise determination of GRB properties and are often unavoidable. Therefore, extensive categorization of GRB LCs remains a hurdle. We address LC gaps using a stochastic reconstruction, wherein we fit two preexisting models (the Willingale model; W07; and a broken power law; BPL) to the observed LC, then use the distribution of flux residuals from the original data to generate data to fill in the temporal gaps. We also demonstrate a model-independent LC reconstruction via Gaussian processes. At 10% noise, the uncertainty of the end time of the plateau, its correspondent flux, and the temporal decay index after the plateau decreases by 33.3%, 35.03%, and 43.32% on average for the W07, and by 33.3%, 30.78%, 43.9% for the BPL, respectively. The uncertainty of the slope of the plateau decreases by 14.76% in the BPL. After using the Gaussian process technique, we see similar trends of a decrease in uncertainty for all model parameters for both the W07 and BPL models. These improvements are essential for the application of GRBs as standard candles in cosmology, for the investigation of theoretical models, and for inferring the redshift of GRBs with future machine-learning analyses.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Biochemical inhibition of acid phosphatase activity in two mountain spruce forest soils

The product inhibition of the 4-Methylumbelliferyl phosphate (MUB-P) decay rate, a measure of potential acid phosphatase activity, has not been considered in most of the published kinetic studies. The aim of this study was to determine the type and strength of the product inhibition in order to better define reaction conditions at which the acid phosphatase activity assay produces unbiased results. The MUB-P decay rate was measured in the forest floor and topsoil organic horizons of two spruce forest catchments at different initial MUB-P and P-PO4 3- concentrations. Furthermore, the type and strength of the inhibition was analyzed by non-linear regression. We found that MUB-P decay was competitively inhibited by the P-PO4 3- but also by dissolved organic P. By using estimated kinetic parameters, we calculated an underestimation of potential acid phosphatase activity by up to ~20% at concentrations of added MUB-P as high as 300 µmol g-1 dry weight and an incubation time one hour due to background inorganic and organic P concentrations in range of units of µmol per gram of dry soil. We discuss approaches that can be used to minimize the underestimation and show that the previously defined recommendations might not be always sufficient to avoid the bias. Therefore, we recommend to analyze progress curves at a wide range of initial MUB-P concentrations using non-linear methods when the enzyme assay is optimized for a given soil. When changes in inorganic and/or organic P concentrations are expected, we further recommend measuring background inorganic and organic P concentrations.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

End-to-end jet classification of boosted top quarks with the CMS open data

Here we describe a novel application of the end-to-end deep learning technique to the task of discriminating top quark-initiated jets from those originating from the hadronization of a light quark or a gluon. The end-to-end deep learning technique uses low-level detector representation of high-energy collision event as inputs to deep learning algorithms. In this study, we use low-level detector information from the simulated Compact Muon Solenoid (CMS) open data samples to construct the top jet classifiers. To optimize classifier performance we progressively add low-level information from the CMS tracking detector, including pixel detector reconstructed hits and impact parameters, and demonstrate the value of additional tracking information even when no new spatial structures are added. Relying only on calorimeter energy deposits and reconstructed pixel detector hits, the end-to-end classifier achieves an area under the receiver operator curve (AUC) score of 0.975 ± 0.002 for the task of classifying boosted top quark jets. After adding derived track quantities, the classifier AUC score increases to 0.9824 ± 0.0013, serving as the first performance benchmark for these CMS open data samples.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

The impact of r -process heating on the dynamics of neutron star merger accretion disc winds and their electromagnetic radiation

Neutron star merger accretion discs can launch neutron-rich winds of >10 -2 M⊙. This ejecta is a prime site for r-process nucleosynthesis, which will produce a range of radioactive heavy nuclei. The decay of these nuclei releases enough energy to accelerate portions of the wind by ~0.1c. Here, we investigate the effect of r-process heating on the dynamical evolution of disc winds. We extract the wind from a 3D general relativistic magnetohydrodynamic simulation of a disc from a post-merger system. This is used to create inner boundary conditions for 2D hydrodynamic simulations that continue the original 3D simulation. We perform two such simulations: one that includes the r-process heating, and another one that does not. We follow the hydrodynamic simulations until the winds reach homology (60 s). Using time-dependent multifrequency multidimensional Monte Carlo radiation transport simulations, we then calculate the kilonova light curves from the winds with and without dynamical r-process heating. We find that the r-process heating can substantially alter the velocity distribution of the wind, shifting the mass-weighted median velocity from 0.06c to 0.12c. The inclusion of the dynamical r-process heating makes the light curve brighter and bluer at $\sim 1\, \mathrm{d}$ post-merger. However, the high-velocity tail of the ejecta distribution and the early ($\lesssim 1\, \mathrm{d}$) light curves are largely unaffected.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Modeling Kilonova Light Curves: Dependence on Nuclear Inputs

The mergers of binary neutron stars, as well as black hole–neutron star systems, are expected to produce an electromagnetic counterpart that can be analyzed to infer the element synthesis that occurred in these events. We investigate one source of uncertainties pertinent to lanthanide-rich outflows: the nuclear inputs to rapid neutron capture nucleosynthesis calculations. We begin by examining 32 different combinations of nuclear inputs: eight mass models, two types of spontaneous fission rates, and two types of fission daughter product distributions. We find that such nuclear physics uncertainties typically generate at least one order of magnitude uncertainty in key quantities such as the nuclear heating (one and a half orders of magnitude at 1 day post-merger), the bolometric luminosity (one order of magnitude at 5 days post-merger), and the inferred mass of material from the bolometric luminosity (factor of 8 when considering the 8–10 day region). Since particular nuclear processes are critical for determining the electromagnetic signal, we provide tables of key nuclei undergoing β-decay, α-decay, and spontaneous fission important for heating at different times, identifying decays that are common among the many nuclear input combinations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Primary $\gamma$-ray intensities and $\gamma$-strength functions from discrete two-step $\gamma$-ray cascades in radiative proton-capture experiments

Background: Reaction rates of radiative capture reactions can play a crucial role in the nucleosynthesis of heavy nuclei in explosive stellar environments. These reaction rates depend strongly on $\gamma$-ray decay widths in the reaction products, which are, for nonresonant capture reactions at high excitation energies, derived from the $\gamma$-ray strength function and the nuclear level density. Recently, the ratio method was applied to primary $\gamma$ rays observed from ($\textit{d, p}$) reactions and nuclear resonance fluorescence measurements to extract the dipole strength in atomic nuclei and to test the generalized Brink-Axel hypothesis. Purpose: The purpose of this work is to apply the ratio method to primary $\gamma$-ray intensities of the 63,65 Cu($\textit{p},\gamma$) reactions to extract $\gamma$-ray strength information on the nuclei 64,66 Zn. Here, the impact of spin distribution, total $\gamma$-ray decay widths, level densities, and width fluctuations on the application of the ratio method will be discussed. Additionally, by comparing the relative $\gamma$-ray strength at different excitation energies, conclusions on the validity of the generalized Brink-Axel hypothesis can be made. Method: The radiative proton capture reaction measurements have been performed at the HORUS $\gamma$-ray spectrometer of the University of Cologne at one excitation energy for each reaction. Primary $\gamma$-ray intensities have been determined by normalizing secondary $\gamma$-ray transitions in two-step cascades using their absolute branching ratio. The ratio method was applied to the measured primary $\gamma$-ray intensities as well as to previous measurements by Erlandsson et al. at different excitation energies. Results: The relative strength function curve for 64 Zn from our measurement shows no significant deviation from the previous measurement at a different excitation energy. The same is true for 66 Zn where both measurements were at almost the same excitation energy. Absolute $\gamma$-strength function values have been obtained by normalizing the relative curves to quasiparticle random phase approximation calculations because of the absence of experimental data in the respective energy region. Conclusion: The generalized Brink-Axel hypothesis, i.e., the independence of the strength function on the excitation energy, seems to hold in the studied energy region and nuclei. The method to obtain primary $\gamma$-ray intensities from two-step cascade spectra was shown to be a valuable and sensitive tool although its uncertainties are connected to the knowledge of the low-energy level scheme of the investigated nucleus. The scaling in the ratio method should be taken with care, because the relative strength is not a simple sum of $f_{E1}$ and $f_{M1}$ but a somewhat complex linear combination dependent on the excitation energy of the nucleus.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

The 100-month Swift catalogue of supergiant fast X–ray transients

Supergiant fast X-ray transients (SFXTs) are high mass X-ray binaries (HMXBs) displaying X-ray outbursts that can reach peak luminosities up to 10 38 erg s -1 and spend most of their lives in more quiescent states with luminosities as low as 10 32 -10 33 erg s -1 . During the quiescent states, less luminous flares are also frequently observed with luminosities of 10 34 -10 35 erg s -1 . The main goal of the comprehensive and uniform analysis of the SFXT Swift triggers presented in this paper is to provide tools to predict whether a transient that has no known X-ray counterpart may be an SFXT candidate. These tools can be exploited for the development of future missions exploring the variable X-ray sky through large field-of-view instruments. We examined all available data on outbursts of SFXTs that triggered the Swift/Burst Alert Telescope (BAT) collected between 2005 August 30 and 2014 December 31, in particular those for which broad-band data, including the Swift/X-ray Telescope (XRT) data, are also available. This work complements and extends our previous catalogue of SFXT flares detected by BAT from 2005 February 12 to 2013 May 31, since we now include the additional BAT triggers recorded until the end of 2014 (i.e. beyond the formal first 100 months of the Swift mission). Due to a change in the mission’s observational strategy, virtually no SFXT triggers obtained a broad-band response after 2014. We processed all BAT and XRT data uniformly by using the Swift Burst Analyser to produce spectral evolution dependent flux light curves for each outburst in the sample. The BAT data allowed us to infer useful diagnostics to set SFXT triggers apart from the general γ-ray burst population, showing that SFXTs uniquely give rise to image triggers and are simultaneously very long, faint, and ‘soft’ hard-X-ray transients. We find that the BAT data alone can discriminate very well the SFXTs from other classes of fast transients, such as anomalous X-ray pulsars and soft gamma repeaters. On the contrary, the XRT data collected around the time of the BAT triggers are shown to be decisive for distinguishing SFXTs from, for instance, accreting millisecond X-ray pulsars and jetted tidal disruption events. The XRT observations of 35 (out of 52 in total) SFXT BAT triggers show that in the soft X-ray energy band, SFXTs display a decay in flux from the peak of the outburst of at least three orders of magnitude within a day and rarely undergo large re-brightening episodes, favouring in most cases a rapid decay down to the quiescent level within three to five days (at most).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Public Reference Data for Megawatt-Scale Hydrogen Electrolysis - NLR Historical Solar PV

The U.S. Department of Energy and National Laboratory of the Rockies (NLR) demonstrate hydrogen electrolysis from variable sources, hydrogen compression and storage, and hydrogen fuel cell power production using megawatt-scale equipment at NLR’s Flatirons Campus as part of the Advanced Research on Integrated Energy Systems (ARIES) research platform. This dataset represents part of that effort and is intended for academic, national laboratory, industrial, and other stakeholders to plan, design, and validate models of megawatt-scale hydrogen technologies and diverse energy infrastructure nationwide. These data provide a baseline for how existing hydrogen electrolysis technologies perform when coupled with various energy technologies. Future datasets will demonstrate how existing hydrogen fuel cell technologies can provide controllable, dispatchable, and variable power output for artificial intelligence data centers and other variable loads. This dataset entry describes the behavior of a 1.25-MW proton exchange membrane MC250 electrolyzer system, manufactured by Nel Hydrogen , [1] when fed historical data generated by the 430-kW, fixed-axis solar photovoltaic (PV) array located at NLR’s Flatirons Campus. (While the electrolyzer balance of plant supports up to 2.5 MW of electrolysis, NLR only has a single 1.25-MW electrolysis stack.) Solar PV power output data for the 2020 calendar year were categorized on a daily basis by total energy generation and standard deviation. Each day was then ranked by these metrics, and the 25th, 50th, and 100th percentiles were selected. The 75th percentile day did not exhibit sufficient variability to make for a valuable experiment. A similar process was used for the related historical wind dataset . [2] The historical days in 2020 that represented these percentiles are Dec. 19, March 29, and May 4, respectively. The entire solar day’s power profile was then fed through the MC250 electrolyzer. Due to its length, the 100th percentile day experiment was split into two parts, and the final 3 hours of the solar day were not captured. These final 3 hours contained no spikes or dips of interest and simply represented a slow decay of input solar power. Also, a single timestamp (13:13:47 on Jan. 14, 2026) was lost in the hydrogen system supervisory control and data acquisition. Finally, during the 25th percentile experiment (solar day Dec. 19, 2020) data recording was lost from 11:00:13 to 11:14:45. The roughly 15 minutes of the solar profile were rerun at the end of the experiment and spliced into this time slot during post-processing. The electrolysis system controls hydrogen production by varying direct current applied to the stack, from a maximum of 3,000 A to a minimum safe operation of 300 A, or 10%. Because the current–voltage characteristic changes as the stack ages and efficiency degrades, the actual minimum safe operating power changes over time. The historical solar profiles were translated from power (kilowatts) to current (amperes) using a curve fit with calibration data and sent to the electrolyzer power supply at 1-Hz frequency. For more details on the statistical analysis process, see the slide deck “Public Reference Data for Megawatt-Scale Hydrogen Electrolysis: NLR Historical Solar PV Analysis and Profile Generation” accessible with this data entry. These datasets report relevant hydrogen balance-of-plant and system data, all captured at 1 Hz, including hydrogen mass production measured with an Emerson Coriolis flow meter. Each .zip file represents a single solar PV electrolysis experiment and is formatted as: {technology}_{percentile}_{scaling factor} For instance, “solarPV-430kW_25_2x.zip” reports the experiment using the 25th percentile solar data from the historical 2020 solar PV dataset, scaled to 200%. Scaling factors were applied to the generated solar PV power output files to more closely match the 1.25-MW capacity of the electrolyzer. Each .zip folder contains the following files: A .csv file containing raw data. An .xlsx file explaining all the fields in the raw data. A .png plot showing the time series of hydrogen production, electrolysis power consumption, and solar power input. A PDF file detailing the historical solar data statistical analysis used to generate the solar profile. An experiment labeled “characterization_200.zip” demonstrates the MC250 electrolyzer steady-state response with 30-minute load steps for a total duration of 5 hours. Finally, a .csv file is provided with all experiments combined into one dataset labeled "combined_solarPV_experiments.csv". [1] nelhydrogen.com/product/mc-series-electrolyser . [2] data.nlr.gov/submissions/316 .

08 HYDROGEN↗
DOE OSTI2024

Vortex bursting and associated twist dynamics on helical vortex tubes and vortex rings

The interaction of opposite-signed twist waves on vortex tubes can lead to vortex bursting, a process where the core expands into a double ring-like structure with strong swirling flows. Previous works have studied vortex bursting on rectilinear vortices by axially perturbing the initial core size to generate the twist waves, and observed largely axisymmetric bursting dynamics. In this work, we numerically study bursting on vortical structures with curved centrelines, analysing the interaction between the centreline dynamics, twist wave generation and propagation, and vortex bursting. We focus on axially perturbed helical vortex tubes with small radius-to-pitch ratios up to$0.0625$, as well as vortex rings with a large radius-to-core size ratio$10$, both at a circulation-based Reynolds number$5000$. The results show that though the initial twist wave propagation speeds are relatively unaffected by the curvature and torsion of the centreline, the bursting process is altered significantly compared with rectilinear vortices. The self-induced rotation of the centreline of the helical tube induces a non-axisymmetric distortion of the bursting structure, which rapidly breaks up the vortex core into small-scale helical structures. A similar destabilization of the bursting structure also occurs on vortex rings. The enstrophy increase and accelerated energy decay associated with bursting are predominantly determined by the twist wave strength, rather than the curvature and torsion of the centreline. Combined, our findings imply that bursting could play an important role in transferring and dissipating energy of vortical structures in wakes, and turbulent flows in general.

Mechanics↗
DOE OSTI2023

Exploring unsupervised top tagging using Bayesian inference

Recognizing hadronically decaying top-quark jets in a sample of jets, or even its total fraction in the sample, is an important step in many LHC searches for Standard Model and Beyond Standard Model physics as well. Although there exists outstanding top-tagger algorithms, their construction and their expected performance rely on Montecarlo simulations, which may induce potential biases. For these reasons we develop two simple unsupervised top-tagger algorithms based on performing Bayesian inference on a mixture model. In one of them we use as the observed variable a new geometrically-based observable \tilde{A}_{3} A ̃ 3 , and in the other we consider the more traditional \tau_{3}/\tau_{2} τ 3 / τ 2 N N -subjettiness ratio, which yields a better performance. As expected, we find that the unsupervised tagger performance is below existing supervised taggers, reaching expected Area Under Curve AUC \sim 0.80-0.81 ∼ 0.80 − 0.81 and accuracies of about 69% - − 75% in a full range of sample purity. However, these performances are more robust to possible biases in the Montecarlo that their supervised counterparts. Our findings are a step towards exploring and considering simpler and unbiased taggers.

Alvarez, Ezequiel↗
DOE OSTI2021

Cumulative short-lived photofission product yields for nuclear forensic application

Experiments have long been carried out with linear accelerator (linac) produced bremsstrahlung X-rays to study high-energy photon-induced fission, known as “photofission”. One benefit of photofission, when compared with neutron-induced fission, is the ability to investigate fundamental physics of fission at excitation energies lower than what is possible with neutrons [1]. Absorption of an incident neutron results in a compound nucleus with an excitation energy at least equal to the neutron binding energy; thus, the lowest energy region of the fission barrier cannot be studied. Fission product yield distributions are dependent on the incident particle energy and the target nucleus mass. They are generally represented by a double-humped curve, with low and high mass “peaks” and an intermediate mass “valley”. As incident particle energy increases, the excitation energy of the nucleus increases, thereby increasing the overall yield of isotopes in the valley with respect to the two maxima or peaks of the fission product curve [3]. Photofission is most probable for incident photons around 14 MeV due to the giant dipole resonance. Induced photofission around this energy results in a fission product yield distribution that resembles 14 MeV DT neutron fission more than thermal or fast (500 keV) neutron fission. In post-detonation nuclear forensics, the peak-to-valley ratio of the fission product yield curve could be used for characterizing device type by providing an indication of the neutron energy spectrum that had induced fission (e.g., 500 keV vs. 14 MeV). The goal of this work is to improve nuclear data for isotopes of interest in nuclear forensics and special nuclear material (SNM) detection, specifically the fission product yield data via photofission experiments for the isotopes of 238U and 232Th. Photofission has been proposed as a less costly and more flexible production method for the valley isotopes that are representative of DT neutron fission product yield distributions, which would benefit nuclear forensics exercises in testing for analysis techniques. Natural uranium and thorium targets were irradiated with bremsstrahlung photons (nominal endpoint energies of 8, 14, and 20 MeV), generated via impingement of high-energy electrons on a tungsten electron-to-photon radiator. Gamma-ray spectroscopy was used to observe and identify the decay of the resulting fission products. The photofission yields of 102 fission products were measured for the set of the 238U and 232Th targets at the three bremsstrahlung X-ray endpoint energies 8, 14, and 20 MeV. These fragments included masses ranging from A=84 to A=144, with half-lives as short as 1.07 seconds.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2024

Quantification of Residual Water in Spent Fuel Dry Storage Canisters Using Mass Spectrometry

The purpose of this report is to document updates on the apparatus to simulate commercial vacuum drying procedures at the Nuclear Energy Work Complex at Sandia National Laboratories. Validation of the extent of water removal in a dry spent nuclear fuel storage system based on drying procedures used at nuclear power plants is needed to close existing technical gaps. Operational conditions leading to incomplete drying may have potential impacts on the fuel, cladding, and other components in the system during subsequent storage and disposal, such as fuel degradation; cladding corrosion, embrittlement, or breaching; and the creation of a flammable environment via radiolysis of water. A general lack of data suitable for model validation of commercial nuclear canister drying processes necessitates well-designed investigations of drying process efficacy and water retention. Scaled tests that incorporate relevant physics and well-controlled boundary conditions are essential to provide insight and guidance to the simulation of prototypic systems undergoing drying processes. This report documents details on the quantification of residual water in the Advanced Drying Cycle Simulator (ADCS), an apparatus built to simulate commercial drying procedures and quantify the amount of residual water remaining in a pressurized water reactor (PWR) fuel assembly after drying. The ADCS was constructed with a prototypic 17×17 PWR fuel skeleton and waterproof heater rods to simulate decay heat. The ADCS is outfitted with thermocouples to measure the thermal response of the ADCS to simulated decay heats and internal helium fill pressures relevant to commercial drying procedures. The ADCS is also instrumented with pressure transducers to measure the pressures and vacuum levels observed during simulated commercial drying. The most unique instrumentation used for quantifying residual water in the ADCS is a Hiden Analytical HPR-30 mass spectrometer (MS), which measures gas compositions of the ADCS internal free volume, based on partial pressures calculated from relative proportions of gas molecules detected by the MS. This report details the methodology used to implement MS measurements in quantifying residual water in the ADCS. This methodology includes the calibration of the HPR-30 MS to a Buck Research Instruments CR-4 chilled mirror hygrometer, which itself is calibrated to a NIST-traceable standard. Data collected by both the MS and the chilled mirror hygrometer from water/helium mixtures ranging from 150 to 500,000 ppmv water in helium were used to generate calibration curves, establishing a source of verification of MS measured water contents. Details regarding water content measurement uncertainties are included in this report, defining the accuracy and verifiability of the HPR-30 MS in measuring residual water content in simulated dry storage canister environments.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

PORFLOW modeling of Vadose zone flow and transport for the E-area intermediate level vault

In support of the E-Area Performance Assessment, a two-dimensional model of water flow and radionuclide transport through the E-Area Intermediate Level Vault (ILV) and local vadose zone has been developed using the PORFLOWTM software. The purpose of the model is to calculate flux to the water table for radionuclides eluted from the ILV during its operational life, the period of institutional control, and times following site closure. Results of model calculations will be used by a three-dimensional PORFLOW model of transport through the aquifer to determine radionuclide concentrations at a hypothetical 100 meter well and at the site boundary where contaminated groundwater could be assessed by members of the public following site closure. While newly developed, the model structure closely follows that used in the 2008 PA while incorporating a refined computational mesh, updated material properties for the vadose zone soil and vault concrete, revised infiltration rates, and new solid-liquid distribution coefficients. The model also addresses degradation of concrete hydraulic properties by blending soil and concrete water retention curves over a 500 year period. That is, the hydraulic properties of ILV concrete start out as fresh concrete at the time of site closure and degrade to soil properties over a 500 year time period. This approach has not been used previously in E-Area PA’s. The model has been used to calculate water flow through the ILV vadose zone and, on a trial basis, transport of a limited number of radionuclides from the waste region inside the ILV to the water table. Figure ES-1 provides a view of the model and an example of water flow after the concrete has substantially degraded. Figure dimensions are in centimeters. The red area in the left hand side figure is the waste disposal region within the ILV. Colored bands in the figures indicate different soil and concrete regions. Figure ES-2 shows a calculation of flux to the water table (mol/yr) for U-234 and its one-year half-life decay chain (U-234 → Th-230 → Ra-226 → Pb-210). Pb-210 flux is too small to show on the same scale used for Th-230 and Ra-226 in Figure ES-2. The model has also been used to make a preliminary evaluation of the release of tritium from TPBAR disposal containers that will be place in the ILV. A separately calculated tritium release rate as a function of time was used as a source in the PORFLOW model. The TPBAR disposal schedule and tritium release rate differ from the values used in the 2008 PA. Results from these initial trial runs are compared to calculations made for the 2008 PA ILV analysis. In brief, a new PORFLOW model of flow and transport through the ILV vadose zone has been developed and tested. The model is ready to be used for PA calculations.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Preliminary Design of a New Pulsed-Neutron Die--Away Experiments with Absorbers (IER-552 CED-1)

This report presents the preliminary design (CED-1) of IER-552, pulsed-neutron die-away (PNDA) experiments with absorbers, to be conducted by Lawrence Livermore National Laboratory (LLNL). The goal of the PNDA experiments with absorbing materials is to produce high-quality, low-cost integral benchmarks that can be used to validate cross sections that are high priority nuclear data for the Department of Energy’s Nuclear Criticality Safety Program (NCSP). PNDA experiments can be used to generate nuclear data validation benchmarks quickly and cheaply, as they do not require fissile material or a nuclear facility to conduct them. The original PNDA experiments were set up at LLNL under IER-501. Extending the PNDA experiments to include absorbing materials leverages existing facilities and resources at LLNL to validate cross sections, with a basic experiment design that uses a neutron generator to impinge a short, mono-energetic neutron pulse on a poisoned, moderating target sample. The neutron pulse lasts 100 µs, after which the neutron population in the sample reaches thermal equilibrium and later spatial equilibrium. In this fully equilibrized state, the neutron population has a characteristic decay-time eigenvalue. The eigenvalue can be used as an integral parameter to validate nuclear data involved with neutron scattering and absorption. Two observables can be measured: the neutron population and the population of secondary γ-rays produced by inelastic and capture reactions. This report explores the feasibility of both approaches. For γ-ray measurements, this allows for direct measurement of reactions of interest for validation. During the neutron burst time, inelastic (prompt) γ-rays dominate, and these data can be isolated for validation of neutron inelastic scattering data. After the pulse ends, neutron capture (delayed) γ-rays dominate and are produced by many reaction channels as neutrons slow down to thermal energies. Two sets of γ-ray spectra (inelastic, capture) can be produced for validation. Additionally, the time-dependent decay of the capture γ-ray population when the neutron population is thermalized and in its fundamental mode can be used for validation. Two experiments can be imagined where pulsed neutrons are used to validate the nuclear data of absorbing materials. The first integrates well into the existing PNDA testbed at LLNL. The absorber material would be introduced into a moderating target and the same neutron generator, neutron detectors, electronics, and shielding box used for IER-501 would be employed. From the same experimental setup, many absorbing materials can be validated by simply changing the neutron poison in the moderator. The only significant modification would be the introduction of a γ-ray detector. Such an experiment would primarily serve to validate nuclear data in the thermal energy range. The second kind of experiment would validate the fast and epithermal energy range. It would involve pulsing neutron into a bare assembly of the absorber material and measuring the die-away of neutrons. The fast leakage of such an assembly requires nanosecond neutron-generator pulse widths and fast detectors/electronics that could measure ~200 ns die-away curves. This kind of experiment would ideally be performed at a cyclotron, like that available at Lawrence Berkeley National Laboratory.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Advances in phosphor two-color ratio method thermography for full-field surface temperature measurements

We report thermographic phosphors can be employed for optical sensing of surface, gas phase, and bulk material temperatures through different strategies including the time-decay method, time-integrated method, and frequency-domain method. We focus on the time-integrated method, also known as the ratio method, as it can be more practical in many situations. This work advances the ratio method using two machine vision cameras with CMOS detectors for full-field temperature measurements of a solid surface. A phosphor calibration coupon is fabricated using aerosol deposition and employed for in situ determination of the temperature-versus-intensity ratio relationship. Algorithms from digital image correlation are employed to determine the stereoscopic imaging system intrinsic and extrinsic parameters, and accurately register material points on the sample to subpixel locations in each image with 0.07 px or better accuracy. Detector nonlinearity is carefully characterized and corrected. Temperature-dependent, spatial non-uniformity of the full-field intensity ratio—posited to be caused by a blue-shift effect of the bandpass filter for non-collimated light and/or a wavelength-dependent transmission efficiency of the lens—is assessed and treated for cases where a standard flat-field correction fails to correct the non-uniformity. In sum, pixel-wise calibration curves relating the computed intensity ratio to temperature in the range of T = 300–430 K are generated, with an embedded error of less than 3 K. This work offers a full calibration methodology and several improvements on two-color phosphor thermography, opening the door for full-field temperature measurements in dynamic tests with deforming test articles.

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

Constraining Inputs to Realistic Kilonova Simulations through Comparison to Observed r -process Abundances

Kilonovae, one source of electromagnetic emission associated with neutron star mergers, are powered by the decay of radioactive isotopes in the neutron-rich merger ejecta. Models for kilonova emission consistent with the electromagnetic counterpart to GW170817 predict characteristic abundance patterns, determined by the relative balance of different types of material in the outflow. Assuming that the observed source is prototypical, this inferred abundance pattern in turn must match r -process abundances deduced by other means, such as what is observed in the solar system. We report on analysis comparing the input mass-weighted elemental compositions adopted in our radiative transfer simulations to the mass fractions of elements in the Sun, as a practical prototype for the potentially universal abundance signature from neutron star mergers. We characterize the extent to which our parameter inference results depend on our assumed composition for the dynamical and wind ejecta and examine how the new results compare to previous work. We find that a dynamical ejecta composition calculated using the FRDM2012 nuclear mass and FRLDM fission models with extremely neutron-rich ejecta ($Y_e$ = 0.035) along with moderately neutron-rich ($Y_e$ = 0.27) wind ejecta composition yields a wind-to-dynamical mass ratio of $M_w$/$M_d$ = 0.47, which best matches the observed AT2017gfo kilonova light curves while also producing the bestmatching abundance of neutron capture elements in the solar system, though, allowing for systematics, the ratio may be as high as of order unity.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Leveraging Machine Learning Capabilities for the Characterization of Irradiated Uranium: A Case Study of Analysis Methods for Nuclear Safeguards and Nuclear Forensics

Nondestructively determining the initial enrichment of irradiated uranium is a complex and laborious multivariable problem due to the presence of fission products. This work demonstrates the capabilities of machine learning to analyze gamma-ray spectral data to determine initial enrichment without knowledge of the decay time of the sample. The approach developed is agnostic to the particular scenario and is applicable to a wide variety of applications in nuclear forensics and nuclear safeguards. We irradiated 5 mg uranium standard reference materials at discrete enrichment values ranging from 0.02% to 97% 235 U (weight percent) in UT Austin’s Nuclear Engineering Teaching Laboratory TRIGA Mark II 1.1 MW research reactor, allowed each to decay for 8 hours, and then measured each sample via gamma-ray spectrometry for 50 hours post-irradiation yielding 1,400 individual gamma-ray spectra discretized into 8,192 energy bins. We then trained decision trees models to analyze individual gamma-ray spectra and estimate the associated initial enrichment without knowledge of the time since end of irradiation. We evaluated the performance of the models with a reserved test set not used for training or calibrating the model. A decision tree model constructed with this procedure achieved a mean absolute error in initial enrichment determination of 2.3% (weight percent 235 U). Next, we implemented a principal component analysis pre-processing routine of the gamma-ray spectrometry data to reduce the dimensionality of the dataset from 8,192 channels in the spectrum to 10 principal components while retaining over 99% of the inherent variance in the data. Decision tree models constructed with these data demonstrated decreased mean absolute error in enrichment determination, reduced computational time, and decreased complexity. A single decision tree model constructed with this procedure achieved a mean absolute error in initial enrichment determination of 0.05% (weight percent 235 U). Furthermore, we analyzed these models with learning curves to ensure that overfitting did not occur. The capabilities provided by these models can be naturally extended to other application-focused measurements in the fields of nuclear safeguards, nuclear forensics, and nuclear non-proliferation.

Drescher, Adam↗