Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Neutron source localization”

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 37 records · Page 2

Gamma Ray Source Localization for Time Projection Chamber Telescopes Using Convolutional Neural Networks

Diverse phenomena such as positron annihilation in the Milky Way, merging binary neutron stars, and dark matter can be better understood by studying their gamma ray emission. Despite their importance, MeV gamma rays have been poorly explored at sensitivities that would allow for deeper insight into the nature of the gamma emitting objects. In response, a liquid argon time projection chamber (TPC) gamma ray instrument concept called GammaTPC has been proposed and promises exploration of the entire sky with a large field of view, large effective area, and high polarization sensitivity. Optimizing the pointing capability of this instrument is crucial and can be accomplished by leveraging convolutional neural networks to reconstruct electron recoil paths from Compton scattering events within the detector. In this investigation, we develop a machine learning model architecture to accommodate a large data set of high fidelity simulated electron tracks and reconstruct paths. We create two model architectures: one to predict the electron recoil track origin and one for the initial scattering direction. We find that these models predict the true origin and direction with extremely high accuracy, thereby optimizing the observatory’s estimates of the sky location of gamma ray sources.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Remote Radiation Sensing Using Aerial and Ground Platforms

Remote sensing of ionizing radiation has a significant role in waste management, nuclear material management and nonproliferation, and radiation safety. Robotic platforms can surpass the number of tasks that are achieved by humans. With this technique, the operator's radiation exposure can be decreased. Remote sensing allows for the evaluation and monitoring of radiological contamination. Gamma-ray and neutron sensors were integrated onto the robotic platforms. This approach allows for the radiation sensor data to be dynamically tracked and mapped thus enabling further analysis of the radiation flux in temporal and spatial domains. The goal is to complete scheduled tasks while the robot is being irradiated. To achieve this, electronic components must be shielded and radiation hardened. CZT Detector: Cadmium Zinc Telluride (CZT) detector technology has been a promising solution for gamma-ray and x-ray measurements. Detector data is transferred to the Odroid minicomputer that controls and powers the module via the USB. Robot Operating System (ROS) was utilized for data acquisition and data fusion. The Mariscotti method was employed for the spectrum analysis. A function was programmed in ROS for the automatic identification of photopeaks. CLYC Detector: A Cs{sub 2}LiYCl{sub 6}:Ce{sup 3+} (CLYC) detector was used for simultaneous medium-resolution gamma-ray measurements and neutron counting. A 2.54 cm diameter photomultiplier tube (PMT) was equipped with a high voltage supply and a miniature digitizer. Gamma-ray excitation: fast core-to-valence luminescence (CVL) with 1 ns decay constant, and prompt Ce{sup 3+} emission with 50 ns decay constant. Neutron excitation: slow cerium self-trapped excitation (Ce{sup 3+} STE), 1000 ns decay constant. Radiation Source Localization: Maximum Likelihood Estimation (MLE) and gradient-based methods were used to locate the position of a radiation source based on measured radiation intensities. Multi-Particle Transport Code FLUKA: Estimation of radiation damage of the electronic components is important in order to optimize the robot's operational time while it is irradiated. Displacement per atom (DPA) represents the radiation damage in materials exposed to the ionizing radiation. Various shielding layers of different thickness t were analyzed (< 5% statistical error). The model of the controller of the UAS was designed in FLUKA. Conclusion: CZT and CLYC detectors were integrated onto the robotic platforms. Radiation source localization and contour mapping using robotic platforms were studied. Functions for data analysis and fusion were developed in ROS. FLUKA code was utilized to analyze DPA values. Layers of low-density and high-density materials were used to shield the UAS electronics.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Dual particle imaging: Applications in security and environmental imaging

Radiation imaging detectors have been applied in a variety of fields for detection, localization and characterization of radioactive sources. In this work, we describe a handheld detection system (H2DPI) that can simultaneously image neutrons and gamma rays, and reconstruct the energy spectra of the incident particles. This system was applied to the detection and characterization of a plutonium metal source in experiments performed at the National Criticality Experiments Research Center, where the 4.5 kg plutonium source was placed approximately 60 cm from the H2DPI. Neutron and gamma-ray images were measured that show good agreement with the known location of the source. Furthermore, the H2DPI was also used to measure the neutron spectrum, which shows the expected shape of a Watt fission spectrum, and the gamma-ray energy spectrum, which shows the energies common to plutonium metal. Finally, we describe the application of the H2DPI in augmented reality, whereby real-time neutron and gamma-ray imaging data can be streamed to the Microsoft HoloLens to allow the user to visualize the source locations.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Local structure analysis of low-temperature neutron pair distribution function coupled with molecular dynamics simulations of CH 4 and CO 2 hydrates from 2 to 210 K

CH 4 hydrates occur naturally and are an abundant potential fuel source with a corresponding risk of potent greenhouse gas release, due to their stability conditions at low temperature and high pressure. Byproduct CO 2 can be exchanged with CH 4 in these natural deposits and can potentially stabilize the hydrates at higher temperatures. CH 4 , CO 2 , and mixed CH 4 -CO 2 hydrates are studied with in situ neutron pair distribution function experiments from 2 to 210 K to investigate the impact of varying the CH 4 -CO 2 guest composition in the gas hydrate structure. These experiments combined with Reverse Monte Carlo analysis allow for the characterization of intermolecular CH 4 and CO 2 interactions with the water molecules which form the hydrate lattice and how they impact the local structure of the lattice itself. Results indicate that when CH 4 and CO 2 co-occupy the hydrate, the host is more strongly distorted than in the pure CH 4 and CO 2 hydrates, but this becomes less defined with increasing temperature. The presence of CO 2 in mixed hydrate increases the stability range and creates a barrier for CH 4 to completely leave the structure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Constraints on light QCD and CP-violating axions from the death line of rotation-powered pulsars

For axions that couple to nucleons, the presence of dense nuclear matter can displace the axion from its vacuum minimum, sourcing large field gradients around neutron stars (and, more generally, compact objects). These gradients, which we refer to as axion hair, couple to the local background magnetic field, inducing a large voltage drop near the surface of the star; here, we demonstrate that the presence of axion hair decouples local near-field particle acceleration in the open magnetic field line bundle from the rotational frequency of the pulsar itself. This is significant as the non-observation of old slowly-rotating pulsars is attributed to the fact the rotationally-induced electric fields are not strong enough to sustain $e^\pm$ pair production. In this work, we review the evidence for the existence for `pulsar death', i.e. the threshold at which $e^\pm$ pair production (and thus, by association, coherent radio emission) ceases, and demonstrate using both semi-analytics and particle-in-cell simulations that the existence of axion hair can dramatically extend pulsar lifetimes. We show that the non-observation of extremely old, slowly rotating, pulsars allows for a new probe of light QCD and CP-violating axions. We also demonstrate how the observation of emission from both poles of pulsars with nearly orthogonal rotational and magnetic axes, as seen e.g. in PSR J1906+0746, can be used to set competitive limits on CP-violating axion-nucleon interactions.

Witte, Samuel J. [Oxford U., Theor. Phys.; DESY; H↗

First-collision source treatment for ray effect mitigation in discrete-ordinate radiation transport solutions

Deterministic transport codes play a fundamental role in the modelling and simulation of neutron transport. One of the most common deterministic methods is the method of discrete ordinates, also known as the S method. While offering significant advantages over other deterministic methods or stochastic methods like Monte Carlo, the method of discrete ordinates suffers from non-physical artifacts in its local solution due to its discretization of angle. These artifacts, referred to as ray effects because of their ray-like appearance, tend to be worse in problems with small sources in areas with little scattering. Significant effort has gone into developing methods to mitigate ray effects, such as the first-collision source treatment, which separates the angular flux into the uncollided and collided fluxes and solving them using non-traditional techniques such as ray tracing. One such code capable of doing this is Lawrence Livermore National Laboratory's deterministic transport code ARDRA. Current ray tracing methods typically trace to a set of points inside a zone to compute an overall flux. However, this approach has significant drawbacks, such as a low order of convergence and not being conservative. Therefore, a new method has been developed that traces instead to a set of points on each of a zone's surfaces and computing the currents, before using these to obtain the flux. A comparison between these two ray tracing methods showed significant advantages to the new surface method, including inherent conservation, a higher convergence rate, and an increase in calculable information like leakage. This work performed under the auspices of the U.S. Department of Energy by Lawrence Liver- more National Laboratory under Contract DE-AC52-07NA27344. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Embedding Neural Thermal Scattering (NeTS) Modules in SERPENT for Higher Fidelity Advanced Reactor Analysis

When a neutron born in fission thermalizes to the order of $k$ $B$ $T$, it’s de-Broglie wavelength and energy approach the order of inter-atomic spacing and elementary lattice oscillations, respectively. $S$($a,β,t$) or the scattering law, uuantify these temperature-dependent crystallographic contributions to total cross section (or reaction rate). In a Monte Carlo analysis, cumulative distribution functions (CDFs) of $S$($a,β,t$) are loaded to memory from “A Compact ENDF” (ACE) files for stochastically selecting thermal scattered neutron trajectories. In this work, novel neural thermal scattering (NeTS) modules for $S$($a,β,t$) CDFs are designed, trained, serialized and embedded within SERPENT using Python’s limited C-API for on-the-fly deployment of crystalline graphite $S$($a,β,t$) sampling. Torchscript tracing and Numba just-in-time (JIT) compilation streamline neural inference on NVIDIA GPUs with CUDA libraries. Demonstrations of bare sphere thermalization of fast and thermal sources show excellent agreement between embedded NeTS in SERPENT and MCNP. With an explicit model of the reactor, NeTS can predict on-the-fly changes in TREAT neutron spectra as a function of local temperature, which can serve to improve transient and accident predictions in a multiphysics analysis framework. This framework can be further extended to account on-the-fly for changes in local graphitic microstructure to scattering cross sections, and outlines a novel coupling of modern machine learning with state-of-the-art reactor physics methods.

97 MATHEMATICS AND COMPUTING↗

El-CID: a filter for gravitational-wave electromagnetic counterpart identification

ABSTRACT As gravitational-wave (GW) interferometers become more sensitive and probe ever more distant reaches, the number of detected binary neutron star mergers will increase. However, detecting more events farther away with GWs does not guarantee corresponding increase in the number of electromagnetic counterparts of these events. Current and upcoming wide-field surveys that participate in GW follow-up operations will have to contend with distinguishing the kilonova (KN) from the ever increasing number of transients they detect, many of which will be consistent with the GW sky-localization. We have developed a novel tool based on a temporal convolutional neural network architecture, trained on sparse early-time photometry and contextual information for Electromagnetic Counterpart Identification (El-CID). The overarching goal for El-CID is to slice through list of new transient candidates that are consistent with the GW sky localization, and determine which sources are consistent with KNe, allowing limited target-of-opportunity resources to be used judiciously. In addition to verifying the performance of our algorithm on an extensive testing sample, we validate it on AT2017gfo – the only EM counterpart of a binary neutron star merger discovered to date – and AT2019npv – a supernova that was initially suspected as a counterpart of the GW event, GW190814, but was later ruled out after further analysis.

79 ASTRONOMY AND ASTROPHYSICS↗

Time series anomaly detection in power electronics signals with recurrent and ConvLSTM autoencoders

The anomalies in the high voltage converter modulator (HVCM) remain a major down time for the spallation neutron source facility, that delivers the most intense neutron beam in the world for scientific materials research. In this work, we propose neural network architectures based on Recurrent AutoEncoders (RAE) to detect anomalies ahead of time in the power signals coming from the HVCM. Bi-directional gated recurrent unit, bi-directional long-short term memory (LSTM), and convolutional LSTM (ConvLSTM) are developed, trained, and tested using real experimental signals from the HVCM module. The results show a good performance of the proposed RAE models, achieving precision up to 91%, recall up to 88%, false omission rate as low as 20% (i.e. 80% of the anomalies were detected), and area under the ROC curve up to 0.9. The three RAE models provide very comparable performance, with LSTM showing slightly better performance than GRU and ConvLSTM. The RAE models are benchmarked against other anomaly detection methods, including isolation forest, support vector machine, local outlier factor, feedforward and convolutional autoencoders, and others; showing a better performance. Here, the results of this study demonstrate the promising potential of RAE in anomaly detection for real-world power systems, and for increasing the reliability of the HVCM modules in the spallation neutron source.

42 ENGINEERING↗

Constrained variational optimization of counting-time allocation in sequential scattering measurements: Application to Bonse–Hart USANS

Sequential scattering measurements are often performed under a fixed experimental-time budget, even though the expected count rate varies strongly across the measured coordinate. When the dwell time at each measurement position can be controlled independently, this variation creates a general resource-allocation problem: how should the available time be distributed to minimize the uncertainty of the reconstructed profile? We formulate this problem as a constrained variational optimization for measurements governed by Poisson counting statistics. When each measurement is treated independently, minimizing the averaged squared relative uncertainty yields an inverse-square-root intensity allocation. The formulation is then generalized to include correlations between neighboring measurements and an instrumental resolution operator, leading to an allocation criterion that equalizes the marginal reduction in posterior uncertainty per unit measurement time. Bonse–Hart ultra-small-angle neutron scattering (USANS), in which reciprocal space is sampled sequentially through analyzer-angle stepping, provides an experimentally grounded application. Computational benchmarking shows that the optimized allocation outperforms uniform-time and constant-relative-error strategies, while application to an experimentally measured graphite USANS profile from the Spallation Neutron Source, using Poisson resampling under alternative schedules, demonstrates how counting time should be redistributed toward weak-intensity regions under an identical total duration. The resulting framework applies to sequential scattering and related scanning measurements whenever local dwell times are adjustable and directly determine the measurement uncertainties, and when the relevant correlation and instrumental-response models are available.

Tung, Chi-Huan [ORNL] (ORCID:0000000221972074)↗

R-process-rich Stellar Streams in the Milky Way*

We present high-resolution Magellan/MIKE spectra of 22 bright (9 < V < 13.5) metal-poor stars (-3.18 < [Fe/H] < -1.37) in three different stellar streams, the Helmi debris stream, the Helmi trail stream, and the ω Centauri progenitor stream. We augment our Helmi debris sample with results for 10 stars by Roederer et al. for a total of 32 stars. Detailed chemical abundances of light elements as well as heavy neutron-capture elements have been determined for our 22 stars. All three streams contain carbon-enhanced stars. For 13 stars, neutron-capture element lines were detectable, and they all show signatures in agreement with the scaled solar r -process pattern, albeit with a large spread of -0.5 < [Eu/Fe] < +1.3. Eight of these stars show an additional small s -process contribution superposed onto their r -process pattern. This could be discerned because of the relatively high signal-to-noise ratio of the spectra given that the stars are close by in the halo. Our results suggest that the progenitors of these streams experienced one or more r -process events early on, such as a neutron star merger or another prolific r -process source. This widely enriched these host systems before their accretion by the Milky Way. The small s -process contribution suggests the presence of asymptotic giant branch stars and associated local (inhomogeneous) enrichment as part of the ongoing chemical evolution by low-mass stars. Stars in stellar streams may thus be a promising avenue for studying the detailed history of large dwarf galaxies and their role in halo assembly with easily accessible targets for high-quality spectra of many stars.

Astronomy & Astrophysics↗

Origin, parameters, and underlying deformation mechanisms of propagating deformation bands in irradiated 316L stainless steel

Lüders-type propagating deformation bands were observed in specimens of irradiated 316L stainless steel samples removed from Spallation Neutron Source target vessels after service. Mechanical testing with digital image correlation (DIC) and in-situ tensile testing with scanning electron microscopy electron backscatter diffraction showed that the observed Lüders-type behavior was not related to the known transformation-induced plasticity or twinning-induced plasticity behavior. Instead, the phenomenon occurs at small local strain values before a significant amount of martensite or deformation twins appear in the microstructure. Microstructural analysis and in-situ mechanical test results suggest Lüders-type band formation and propagation were related to the appearance and evolution of defect-free channels—analogous to slip bands. A modified Swift equation with a Ludwigson-like component was offered to rationalize the phenomenon and model the strain-softening processes at small strain values. Finally, the results indicate complex microstructural processes and deformation mechanisms were active at small strain values and underline the benefits of advanced mechanical test approaches such as DIC.

36 MATERIALS SCIENCE↗

Local Structure Analysis and Modelling of Lignin‐Based Carbon Composites through the Hierarchical Decomposition of the Radial Distribution Function

Abstract Carbonized lignin has been proposed as a sustainable and domestic source of activated, amorphous, graphitic, and nanostructured carbon for many industrial applications as the structure can be tuned through processing conditions. However, the inherent variability of lignin and its complex physicochemical structure resulting from feedstock and pulping selection make the Process‐Structure‐Property‐Performance (PSPP) relationships hard to define. In this work, radial distribution functions (RDFs) from synchrotron X‐ray and neutron scattering of lignin‐based carbon composites (LBCCs) are investigated using the Hierarchical Decomposition of the Radial Distribution Function (HDRDF) modelling method to characterize the local atomic environment and develop quantitative PSPP relationships. PSPP relationships for LBCCs defined by this work include crystallite size dependence on lignin feedstock as well as increasing crystalline volume fraction, nanoscale composite density, and crystallite size with increasing reduction temperature.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Search for Ultra-high-energy Photons from Gravitational Wave Sources with the Pierre Auger Observatory

A search for time-directional coincidences of ultra-high-energy (UHE) photons above 10 EeV with gravitational wave (GW) events from the LIGO/Virgo runs O1 to O3 is conducted with the Pierre Auger Observatory. Due to the distinctive properties of photon interactions and to the background expected from hadronic showers, a subset of the most interesting GW events is selected based on their localization quality and distance. Time periods of 1000 s around and 1 day after the GW events are analyzed. No coincidences are observed. Upper limits on the UHE photon fluence from a GW event are derived that are typically at ~7 MeV cm -2 (time period 1000 s) and ~35 MeV cm -2 (time period 1 day). Due to the proximity of the binary neutron star merger GW170817, the energy of the source transferred into UHE photons above 40 EeV is constrained to be less than 20% of its total GW energy. These are the first limits on UHE photons from GW sources.

79 ASTRONOMY AND ASTROPHYSICS↗

Local structure and distortions of mixed methane-carbon dioxide hydrates

A vast source of methane is found in gas hydrate deposits, which form naturally dispersed throughout ocean sediments and arctic permafrost. Methane may be obtained from hydrates by exchange with hydrocarbon byproduct carbon dioxide. It is imperative for the development of safe methane extraction and carbon dioxide sequestration to understand how methane and carbon dioxide co-occupy the same hydrate structure. Pair distribution functions (PDFs) provide atomic-scale structural insight into intermolecular interactions in methane and carbon dioxide hydrates. We present experimental neutron PDFs of methane, carbon dioxide and mixed methane-carbon dioxide hydrates at 10 K analyzed with complementing classical molecular dynamics simulations and Reverse Monte Carlo fitting. Mixed hydrate, which forms during the exchange process, is more locally disordered than methane or carbon dioxide hydrates. The behavior of mixed gas species cannot be interpolated from properties of pure compounds, and PDF measurements provide important understanding of how the guest composition impacts overall order in the hydrate structure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Axisymmetric hydrodynamics in numerical relativity using a multipatch method

In this paper, we describe a method of implementing the axisymmetric evolution of general-relativistic hydrodynamics and magnetohydrodynamics through modification of a multipatch grid scheme. In order to ease the computational requirements required to evolve the post-merger phase of systems involving binary compact massive objects in numerical relativity, it is often beneficial to take advantage of these system's tendency to rapidly settle into states that are nearly axisymmetric, allowing for 2D evolution of secular timescales. We implement this scheme in the spectral Einstein code and show the results of application of this method to four test systems including viscosity, magnetic fields, and neutrino radiation transport. Our results show that this method can be used to quickly allow already existing 3D infrastructure that makes use of local coordinate system transformations to be made to run in axisymmetric 2D with the flexible grid creation capabilities of multipatch methods. Our code tests include a simple model of a binary neutron star postmerger remnant, for which we confirm the formation of a massive torus which is a promising source of post-merger ejecta.

79 ASTRONOMY AND ASTROPHYSICS↗

GW190814: Gravitational Waves from the Coalescence of a 23 Solar Mass Black Hole with a 2.6 Solar Mass Compact Object

We report the observation of a compact binary coalescence involving a 22.2-24.3 M ⊙ black hole and a compact object with a mass of 2.50-2.67 M ⊙ (all measurements quoted at the 90% credible level). The gravitational-wave signal, GW190814, was observed during LIGO's and Virgo's third observing run on 2019 August 14 at 21:10:39 UTC and has a signal-to-noise ratio of 25 in the three-detector network. The source was localized to 18.5 deg 2 at a distance of 241 +41-41 Mpc; no electromagnetic counterpart has been confirmed to date. The source has the most unequal mass ratio yet measured with gravitational waves , 0.112+0.0090.008, , and its secondary component is either the lightest black hole or the heaviest neutron star ever discovered in a double compact-object system. The dimensionless spin of the primary black hole is tightly constrained to ≤0.07. Tests of general relativity reveal no measurable deviations from the theory, and its prediction of higher-multipole emission is confirmed at high confidence. We estimate a merger rate density of 1-23 Gpc -3 yr -1 for the new class of binary coalescence sources that GW190814 represents. Astrophysical models predict that binaries with mass ratios similar to this event can form through several channels, but are unlikely to have formed in globular clusters. However, the combination of mass ratio, component masses, and the inferred merger rate for this event challenges all current models of the formation and mass distribution of compact-object binaries.

79 ASTRONOMY AND ASTROPHYSICS↗

Confined magnons

Magnetic structures are known to possess magnon excitations confined to their surfaces and interfaces, but these spatially localized modes are often not resolved in spectroscopy experiments. In this work, we develop a theory to calculate the confined magnon spectra and its associated spin scattering function, which is the physical observable in neutron and electron scattering, and a proxy for photon spectroscopy based on x-ray, Raman, and terahertz (THz) sources. We show that extra anisotropy at the surface or interface plays a key role in magnon confinement. We obtain analytical expressions for the confinement length scale and show that it is qualitatively similar for ferromagnets and antiferromagnets in dimension d≥2. For d=1, we find remarkable differences between ferromagnetic and antiferromagnetic models. The theory indicates the presence of several confined magnon resonances in addition to the usual magnons thought to explain the excitations of magnetic nanostructures. Detecting these modes may elucidate the impact of the interface on spin anisotropy and magnetic order.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗