Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Spectrum Unfolding”

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 19 records

Spectrum Unfolding with the MC-15

The Multiplicity Counter 15 tube detector or MC-15 is an optimized detector designed for use in the field. It is composed of 15 3 He tubes embedded in high density polyethylene (HDPE). Recent work has explored expanding the use of the MC-15 beyond multiplicity counting to neutron dosimetry applications. Knowledge of the neutron energy spectrum information is required to use a detector as a neutron dosimeter. The MC-15 tube layout is shown in Figure 1. The unique layout makes it possible to use the detector for neutron spectroscopy via spectrum unfolding. Spectrum unfolding requires (1) energy dependence of the detector response, (2) a detector response matrix that precisely quantifies the response to mono-energetic neutrons, (3) an initial guess spectrum, (4) an unfolding algorithm, and (5) measured data (counts in the case of the MC-15). An energy dependent detector response matrix (DRM) can be constructed by considering either each of the three rows of 3 He tubes as a distinct detector or each individual tube as a distinct detector. The HDPE separating the 3 He in the MC-15 provides the distinct energy dependent response for the rows and individual tubes. In this report we detail the development of detector response matrices for the MC-15 and the application of the Los Alamos Unfolding Code (LUC) to both simulated and measured data. Three MC-15 orientations were studied: (1) standard orientation with the MC-15 front facing the source, (2) standard orientation with Cd sheet, (3) 90° orientation with the side of the MC-15 facing the source.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Data Augmentation for Neutron Spectrum Unfolding with Neural Networks

Neural networks require a large quantity of training spectra and detector responses in order to learn to solve the inverse problem of neutron spectrum unfolding. In addition, due to the under-determined nature of unfolding, non-physical spectra which would not be encountered in usage should not be included in the training set. While physically realistic training spectra are commonly determined experimentally or generated through Monte Carlo simulation, this can become prohibitively expensive when considering the quantity of spectra needed to effectively train an unfolding network. In this paper, we present three algorithms for the generation of large quantities of realistic and physically motivated neutron energy spectra. Using an IAEA compendium of 251 spectra, we compare the unfolding performance of neural networks trained on spectra from these algorithms, when unfolding real-world spectra, to two baselines. We also investigate general methods for evaluating the performance of and optimizing feature engineering algorithms.

McGreivy, James (ORCID:0000000321723411)↗

Impact of Source Geometry on Detector Response Matrix Efficiency: Simulations of the PFUNS-MUSiC Experiments

This work looks at two different bare highly enriched uranium (HEU) systems, configuration one from Measurements of Uranium Subcritical and Critical (MUSiC) and Prompt Fission Uranium Neutron Spectrum (PFUNS) to see the effect geometry has on detector efficiency and the response matrix. The distance from the multiplying source, the medium between the source and detector, and the geometry of the source all play a part in the efficiency of the detector. These factors are especially important when performing spectrum unfolding. Spectrum unfolding is a process used to reconstruct a true spectrum from measured detector response data. It involves interpreting a set of measured values, such as a light signal from a scintillator, to recover the original neutron energy spectrum. This leads to the question, how much will the efficiency of a detector change when you measure a point source compared to an extended source? The distance and solid angle from the source to the detector may be different. It may only be a minuscule change, but in the spectrum unfolding process it can have a considerable effect on the unfolded spectrum

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Validation and Verification of Python based Neutron Spectrum Unfolding Software

To validate and verify the python-based code (PySL), designed to replicate the programs used by STAYSL for Beam Correction Factor (BCF) and Self-Shielding Factor (SHIELD), a series of tests were performed. To test BCF a python script was written to generate a random flux history file and both versions of the code processed the data. The test verified matching values up to at least one decimal place, approximately 10,000 tests where run and each one passed. Isotopes began to fail the tests once neutron saturation was reached. To verify this the total time of exposure was varied the isotopes that failed were compared to a list of their half-lives. The test process for SHIELD was very similar but, in this case, the code began by producing an input file with varying thickness and device type/environment for the SHIELD input. The failure condition for this test was if any of the data points for an isotope had a difference above 3%. Approximately 40 of these tests were run and there were only 3 isotopes that had reoccurring failures but only 2% of their points were above the 3% difference. A visual comparison was conducted by plotting the results from both programs. Although the test failed, the differences between their values were minuscule, and the self-shielding factor’s shape was preserved when plotted. Next steps for this project will be validating and verifying the python-based SigPhi code and then reproducing and testing the least squares unfolding performed by STAYSL.

73 - NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Development of a compact fast-neutron spectrometer for nuclear emergency response applications

We have developed a Compact Fast Neutron Spectrometer (CFNS) for passive assay of special nuclear material (SNM) through the observation of fast neutrons. The CFNS consists of eight organic glass scintillators (OGS) coupled to silicon photomultipliers and a waveform digitizer, which are integrated within a human-portable box. The CFNS determines the neutron energy profile by spectrum unfolding using the Maximum-Likelihood Expectation Maximization method. The detector acquisition system was optimized to have a dynamic range of up to 10 MeV neutron energy. Bulk special nuclear material (SNM) measurements from the National Criticality Experiments Research Center were analyzed for SNM validation/examination. Additionally, the results show that the CFNS can be used to distinguish between fission and (α, n) neutron emitters, regardless of intervening material type (Cu and polyethylene) and thickness, by taking the ratio of neutron counts at different regions in the unfolded energy spectrum. Additionally, by fitting an exponential curve to the unfolded energy spectrum of PuO 2 and Pu neutron emitters, the CFNS showed the ability of distinguishing between pure Pu oxide, pure Pu metal and mixed oxide-metal configurations.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Advancements of the nSpec system

For this study, a new characterization of the detector response for the EJ301 liquid scintillator based nSpec system has been conducted at the Los Alamos Neutron Science Center. The new characterization was accomplished in two orders of magnitude less time and with better statistics than the previous characterization. The new and old characterization methods were compared by conducting both types of characterization on two detectors (one EJ301 and one EJ301D). The two methods show consistent results and compare well to results previously published in the literature; however, the results from the characterization of the original nSpec detector show discrepancies that have been identified as being caused by oxygen being absorbed into the detector due to a leak. The new characterized detector response shows improved spectrum unfolding performance than the previous characterization. The cross talk of the detector system was also characterized as a function of incident neutron energy and found to have a negligible impact on spectrum unfolding results. Alternative materials (EJ301D and stilbene) were also characterized in a first step towards an upgraded system. The different materials show little difference in unfolding performance, but the comparison is ongoing. A SiPM based stilbene detector was used to evaluate the use of SiPMs in a future system upgrade. The SiPM provides adequate performance, but the custom SiPM readout results in a much longer pulse width than PMT based systems. The measured light output of the SiPM based stilbene compares well with the results in literature. Future work will look at additional materials (deuterated stilbene and organic glass) and optimizing the SiPM readout to decrease the pulse width.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Neutron Leakage Spectra of the EUCLID Experiment [Abstract]

Integral experiments with sub-critical and critical configurations of special nuclear material are performed in support of nuclear data validation and adjustment. Different nuclear data evaluations may have different values for individual cross sections due to uncertainties in (or lack of) differential experiments, but compensating errors in these data sets can lead to the same k eff results for one application while vastly different for another application. One example of this is 239 Pu, where both ENDF/B-VIII.0 and JEFF-3.3 correctly compute k eff of the Jezebel critical assembly, but the individual contributions from each reaction are vastly different. To help resolve this specific case, the Experiments Underpinned by Computational Learning for Improvements in Nuclear Data (EUCLID) project used machine learning to design a set of experiments to help resolve the compensating errors in 239 Pu. A total of six responses were measured during the experimental campaign, which constrain the data in ways that k eff alone cannot and will be used for adjustment of the nuclear data. One of these responses is the neutron leakage spectrum, which recent work has shown to be useful for constraining the prompt fission neutron spectrum and inelastic scattering. The neutron leakage spectra were measured utilizing a 3 in. right cylinder EJ-301D detector. The measured signal in the detector was deconvoluted using spectrum unfolding techniques, which are presented and compared to simulations.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Quantification of the light output anistropy in deuterated stilbene

Deuterated stilbene is an organic scintillator that is a desirable material for fast neutron spectroscopy using spectrum unfolding techniques without requiring time-of-flight information. Due to the crystal structure of the material, some anisotropy of the light output exists, which is dependent on the direction of heavy charged particle recoil relative to the crystal structure. The anisotropy of trans-stilbene (hereafter referred to as stilbene) has been well characterized in previous work, but for deuterated stilbene, the anisotropy has only been partially characterized along the a and b crystal axes, while the artificial c' axis, which shows the largest anisotropy in stilbene, has not been characterized until this publication. Here, in this work, two deuterated stilbene crystals were characterized with neutron energies up to 35 MeV at the Los Alamos Neutron Science Center. For one of the crystals, the response is characterized along the a, b, and c' axes. This characterization shows a distinct anisotropy along the axes in deuterated stilbene, which is very similar to that found in regular stilbene, such that the a axis is the brightest, while the b' and c' axes are approximately 3% and 20%–35% lower relative to the a axis.

36 MATERIALS SCIENCE↗

Low-energy measurement of the 25 Mg ⁢(𝛼,𝑛)⁢ 28 Si reaction via neutron spectroscopy

During core helium and carbon burning in massive stars, neutrons are produced mainly by the 22 Ne⁢ (𝛼,𝑛) ⁢25 Mg reaction. Some of these released neutrons are captured by heavy seed nuclei from previous nucleosynthesis events, resulting in the slow production of many of the elements between masses 60 ≤ 𝐴 ≤ 90 via the weak 𝑠 process. Determining the overall neutron flux available in this environment is one of the main challenges in modeling its contributions to heavy element production. Not only must the reaction rate be well defined for the primary neutron source, but the rates of neutron poison and recycling reactions must also be well determined. One critical reaction in the simulation network is 25 Mg ⁢(𝛼,𝑛) ⁢ 28 Si. This reaction, together with 25 Mg ⁢(𝑛,𝛾)⁢ 26 Mg and 25 Mg⁢ (𝛼,𝛾) ⁢ 29 Si, determine how 25 Mg nuclei affect the available neutron flux. Past attempts to measure the 25 Mg ⁢(𝛼,𝑛)⁢ 28 Si cross section using neutron moderator counters have been greatly hindered by lower-𝑍 background reactions. Here, in the present work, neutron spectroscopy with deuterated liquid scintillator detectors has been used. The experimental spectra have been analyzed by applying spectrum unfolding techniques to achieve improved background discrimination for the 25 Mg⁢ (𝛼,𝑛)⁢ 28 Si reaction at low energies, down to 𝐸 𝛼 = 1.75 MeV. The separation of the different background contributions gives further insight into the results of previous moderator counter measurements and the measurements lead to a revised and more reliable determination of the reaction rate.

hydrostatic stellar nucleosynthesis↗

Radiological Source Term Estimation and Isotopic Identification with Parallel Log Domain Particle Filters

This paper presents a parallel log-domain particle filtering algorithm combined with gamma spectrum unfolding to perform localization, identification, and evaluation of multiple point sources of various isotopes in an environment with attenuating obstacles. The method uses sets of precomputed attenuation kernels that map the attenuation characteristics of the environment. These kernels are specific to the energy level of a photopeak of interest. The spectral measurements are deconvolved into count measurements of each photopeak. These count measurements are fed into a set of parallel particle filters using attenuation kernels computed for that photopeak’s energy level. The individual regularized particle filters perform all likelihood calculations in the logarithmic domain to mitigate the effects of particle degeneracy. The output of each particle filter is combined to estimate which isotopes are present as well as their positions and strengths. The performance of the algorithm is characterized in a lab-scale environment using a mobile robot equipped with a gamma ray spectrometer in the presence of up to three different radioactive isotopes simultaneously. The sources were localized to within 10 cm, and their strengths were estimated within 10% of their true values. Furthermore, the isotopes were all correctly identified, and no spurious sources were reported.

42 ENGINEERING↗

Neutron Fluence and Kerma Responses to Bare, Steel-Reflected, and Lucite-Reflected Steady-State Operations of the ORNL Health Physics Research Reactor

The experiments evaluated in this report were conducted at the Health Physics Research Reactor (HPRR), also known as the Fast Burst Reactor. The reactor was designed and built at Oak Ridge National Laboratory (ORNL) in 1961. In this evaluation, a total of 12 experiments were evaluated (6 shield configurations and two responses). Six of them are considered acceptable as benchmark experiments and are thoroughly described: the measurands are bare, steel-shielded and Lucite-shielded neutron fluences and neutron KERMA responses to steady-state operations of the HPRR. It is important to mention that there is a general lack of information in the experimental methodology and neutron spectrum unfolding, and in some cases, unexpected and unexplained discrepancies between calculation and experiment. The HPRR’s unique characteristics were that it is a fast reactor, the use of a Lucite shield as the shielding material, and the use of the neutron fluence as the main benchmark measurand. There are no known, significant correlations between benchmarks in this evaluation and benchmarks in other evaluations.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Determination of particle-dependent light response for defining detector response functions

This study comprehensively characterizes two scintillation detectors, an organic deuterated liquid and a stilbene crystal, concerning their response to charged particles such as electrons, protons, alphas, and carbon nuclei. The Birks function was used to determine the detectors’ response to charged particles, enabling the determination of particle-dependent light responses. The particle-dependent light characterization was used to determine detector response functions (DRF) using Geant in combination with Gamma Detector Response and Analysis Software (GADRAS). These DRFs enhance GADRAS by enabling full spectral analysis functionality for previously underrepresented detector materials. Additionally, neutron-specific DRFs are determined for neutron energy spectrum unfolding applications.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Final Measurement of the 235 U Antineutrino Energy Spectrum with the PROSPECT-I Detector at HFIR

This Letter reports one of the most precise measurements to date of the antineutrino spectrum from a purely 235 U-fueled reactor, made with the final dataset from the PROSPECT-I detector at the High Flux Isotope Reactor. By extracting information from previously unused detector segments, this analysis effectively doubles the statistics of the previous PROSPECT measurement. Further, the reconstructed energy spectrum is unfolded into antineutrino energy and compared with both the Huber-Mueller model and a spectrum from a commercial reactor burning multiple fuel isotopes. A local excess over the model is observed in the 5–7 MeV energy region. Comparison of the PROSPECT results with those from commercial reactors provides new constraints on the origin of this excess, disfavoring at 2.0 and 3.7 standard deviations the hypotheses that antineutrinos from 235 U are solely responsible and noncontributors to the excess observed at commercial reactors, respectively.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Deuterated Xylene-Based Neutron Spectrometry

Deuterated scintillators, such as deuterated xylene (C8D10) have been theorized to offer better pulse shape discrimination (PSD) and neutron spectrometry performance than traditional protiated scintillators due to the anisotropy of neutron scattering on deuterium. However, the light output responses of said deuterated scintillators to monoenergetic has not been characterized in its entirety at both high and low energies. Several time-of-flight experiments were conducted at the Idaho national Laboratory to characterize a pair of EJ-301D deuterated xylene detectors. The low energy detector response to quasi-monoenergetic neutrons was measured using a Cf-252 source, and the high energy detector response was measured using a D-T source and angles of 20, 45, and 60 deg, corresponding to theoretical neutron energies of 13.2, 10.4, and 8.2 MeV. During post-processing, the neutron and gamma bands were clearly able to be separated, even at low light outputs, solidifying the superior PSD capability of deuterated xylene. An exponential of fit of LO=0.63E-2.22*(1-exp(-0.25E)), which has good agreement with previously measured data, was found as the light output response from deuteron recoils. Additionally, measurements of AmBe, D-T, and D-D sources were unfolded using the light output response curve, which resulted in favorable agreement with theoretical values. Further testing of the response matrix by unfolding more sources (such as PuBe) needs to be done, but the unfolding accuracy for currently available sources suggests that deuterated xylene can be used to measure and unfold both fission, monoenergetic, and continuous spectra.

scintillator↗

Introduction, Investigation, and Experimental Validation of a Novel Passive Neutron Spectrometer

Unfolding neutron energy spectra are instrumental for determining personal health effects and calculating dose received. This area of study is heavily researched, and Lawrence Livermore National Laboratory (LLNL) is investigating a passive neutron spectrometer for the purpose of acquiring the information needed to determine personnel dose in the event of a criticality accident. A part of this investigation is presented in this article through the examination of four experimental detector responses (DRs). These four DRs were acquired in the presence of 252 Cf, AmBe, GODIVA, and National Ignition Facility (NIF) neutron sources. An algorithm developed at LLNL was used to unfold the neutron fluence from each of the four DRs, and subsequently, fluence-to-dose conversion factors provided by the American National Standards Institute were used to calculate dose. Additionally, a multistep unfolding process was developed and employed to calculate the effects of both direct (from the source) and indirect (from room return) neutrons. The average error when unfolding the direct DR was less than 8%. The dose from 252 Cf was predicted with only 8% error. The multistep approach allowed for the identification of the low-energy neutrons in the 252 Cf, AmBe, and NIF DRs.

Nuclear Criticality Safety Program (NCSP)↗

Fast neutron leakage spectra of the EUCLID experiment

Special nuclear material in sub-critical and critical configurations measured in integral experiments are important for validation and adjustment of nuclear data. Many different evaluations of nuclear data exist, and these different evaluations can provide different values for individual cross sections that vary due to the uncertainties in differential experiments or lack of such data. For integral experiments, differences in these individual cross sections can have compensating errors, which lead to the same answer. One example of this is the Jezebel critical assembly, where k eff of the system is correctly computed by both ENDF/B-VIII.0 and JEFF-3.3, despite having substantially different underlying evaluated values for specific reactions (such as elastic and inelastic cross sections). To reduce compensating errors in nuclear data, the Experiments Underpinned by Computational Learning for Improvements in Nuclear Data (EUCLID) project has utilized machine learning to design a set of sub-critical and critical experiments. These experiments include slab- and cube-like configurations of 239 Pu in the form of the ZPPR plates. Six different responses were measured on a total of thirteen different configurations. One of these responses, the neutron leakage spectrum, was measured using an EJ301D detector. Finally, the results of the neutron leakage spectra show good agreement (within 1–2 σ ) with the expected spectrum from simulations and will be used in the subsequent nuclear data adjustment done by the EUCLID team.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Passive Confirmation of the Presence Of High-Explosive Material Via Neutron Transmission Spectroscopy

The goal of this project was to explore a novel approach for identifying presence and potential types of high explosives (HE) in treaty-controlled items with passive neutron sources by using passive neutron transmission spectroscopy. We predominantly look to measure relative elemental abundance of Carbon, Hydrogen, Oxygen, and Nitrogen (CHON) since most relevant materials are certain mix of these elements, as shown in Table 1. Previously, most studies exploring potential identification of CHON elemental content of targets in the vicinity of passive neutron source were focused solely on gamma spectroscopy using high resolution gamma detectors. Using neutron transmission spectroscopy with pulse-shape discrimination (PSD) capable organic scintillators is not only a novel idea in of itself, but arguably necessary for accurate identification of the CHON elemental content of an interrogated target. While high resolution gamma spectroscopy in principle can determine a presence of hydrogen and nitrogen, it is effectively blind to a quantitatively measuring areal densities of C, O. The initially proposed approach, described in detail in the project proposal, was to leverage previous Monte Carlo study on neutron transmission spectroscopy using slab shaped target interrogated with well-collimated (“pencil”) neutron beam. In this study, a simulated test CO 2 target was interrogated by a pencil neutron beam and transmitted neutron spectra were measured using PSD capable organic scintillator using MLEM based unfolding technique. To establish relative elemental content of carbon and oxygen, the unfolded neutron spectrum was fit with a parametrized combination of their respective elemental spectral templates. The individual spectral template was calculated by convoluting ENDF neutron crosssection with the known resolution of the PSD detector used in the study. This approach in the studied configuration successfully quantitatively established relative elemental composition of carbon and oxygen.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Quantum chaos on edge

Recently, the physics of many-body quantum chaotic systems close to their ground states has come under intensified scrutiny. Such studies are motivated by the emergence of model systems exhibiting chaotic fluctuations throughout the entire spectrum [the Sachdev-Ye-Kitaev (SYK) model being a renowned representative] as well as by the physics of holographic principles, which likewise unfold close to ground states. Interpreting the edge of the spectrum as a quantum critical point, here we combine a wide range of analytical and numerical methods to the identification and comprehensive description of two different universality classes: the near edge physics of “sparse” and the near edge of “dense” chaotic systems. The distinction lies in the ratio between the number of a system's random parameters and its Hilbert space dimension, which is exponentially small or algebraically small in the sparse and dense case, respectively. Notable representatives of the two classes are generic chaotic many-body models (sparse) and invariant random matrix ensembles or chaotic gravitational systems (dense). While the two families share identical spectral correlations at energy scales comparable to the level spacing, the density of states and its fluctuations near the edge are different. Considering the SYK model as a representative of the sparse class, we apply a combination of field theory and exact diagonalization to a detailed discussion of its edge spectrum. Conversely, Jackiw-Teitelboim gravity is our reference model for the dense class, where an analysis of the gravitational path integral and random matrix theory reveal universal differences to the sparse class, whose implications for the construction of holographic principles we discuss. Published by the American Physical Society 2024

Altland, Alexander (ORCID:0000000229914805)↗