Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Flux reconstruction”

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 235 records · Page 13

Artificial Neural Networks for Determining Magnetospheric Conditions

This chapter presents a neural-network-based technique that allows for the reconstruction of the global, time-varying distribution of some physical quantity Q, that has been sparsely sampled at various locations within the magnetosphere, and at different times. We begin with a general introduction to the problem of prediction and specification, and why it is important and difficult to achieve with existing methods. We then provide a basic introduction to neural networks, and describe our technique using the specific example of reconstructing the electron plasma density in the Earth's inner magnetosphere on the equatorial plane. We then show more advanced uses of the technique, including 3D reconstruction of the plasma density, specification of chorus and hiss waves, and energetic particle fluxes. We summarize and conclude with a general discussion of how machine learning techniques might be used to advance the state-of-the-art in space weather prediction, and insight discovery.

Bortnik, Jacob↗

Probing jet base emission of M87* with the 2021 Event Horizon Telescope observations

We investigate the presence and spatial characteristics of the jet base emission in M87* at 230 GHz, enabled by the significantly enhanced (u,v) coverage in the 2021 Event Horizon Telescope (EHT) observations. The integration of the 12−m Kitt Peak Telescope (USA) and NOEMA (France) stations into the array introduces two critical intermediate-length baselines to SMT (USA) and IRAM 30−m (Spain), providing sensitivity to emission structures at spatial scales of ∼250 μas and ∼2500 μas (∼ 0.02 pc and ∼ 0.02 pc). Without these new baselines, previous EHT observations of the source in 2017 and 2018 lacked the capability to constrain emission on large scales, where a “missing flux” of order ∼1 Jy is expected to reside. To probe these scales, we analyzed closure phases–robust against station-based gain calibration errors–and model the jet base emission using a simple Gaussian component offset from the compact ring emission at spatial separations > 100 μas. Our analysis revealed a Gaussian feature centered at (ΔRA ≈ 320 μ as, ΔDec. ≈ 60 μ as), projected separation of ≈ 5500 AU, with an estimated flux density of only ∼60 mJy, implying that most of the missing flux identified in previous EHT studies had to originate from different, larger scales. Brighter emission at the relevant spatial scales is firmly ruled out, and the data do not favor more complex models. This component aligns with the inferred position of the large-scale jet and is therefore physically consistent with the emission of the jet base. While our findings point to detectable jet base emission at 230 GHz, the limited coverage provided by only two intermediate baselines limits our ability to robustly reconstruct its morphology. Consequently, we treated the recovered Gaussian as an upper limit on the jet base flux density. Future EHT observations with expanded intermediate baseline coverage will be essential to constrain the structure and nature of this component with higher precision.

accretion, accretion disks↗

Neutrino energy reconstruction from semi-inclusive samples

We study neutrino-nucleus charged-current reactions on finite nuclei for the situation in which an outgoing muon and a proton are detected in coincidence, {\it i.e.,} we focus on semi-inclusive cross sections. We limit our attention to one-body current interactions (quasielastic scattering) and assess the impact of different nuclear effects in the determination of the neutrino energy. We identify the regions in phase space where the neutrino energy can be reconstructed relatively well, and study whether the cross section in those regions is significant. Our results indicate that it is possible to filter more than 50\% of all events according to the muon and proton kinematics, so that for the DUNE and T2K fluxes the neutrino energy can be determined with an uncertainty of less than 1\% and 3\%, respectively. Furthermore, we find that the reconstructed neutrino energy does not depend strongly on how one treats the final-state interactions and is not much affected by the description of the initial state. On the other hand, the estimations of the uncertainty on the neutrino energy show important sensitivity to the modeling of the initial state.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Important Role of Ion Flux Regulated by Separators in Lithium Metal Batteries

Abstract Polyolefin separators are the most common separators used in rechargeable lithium (Li)‐ion batteries. However, the influence of different polyolefin separators on the performance of Li metal batteries (LMBs) has not been well studied. By performing particle injection simulations on the reconstructed three‐dimensional pores of different polyethylene separators, it is revealed that the pore structure of the separator has a significant impact on the ion flux distribution, the Li deposition behavior, and consequently, the cycle life of LMBs. It is also discovered that the homogeneity factor of Li‐ion toward Li metal electrode is positively correlated to the longevity and reproducibility of LMBs. This work not only emphasizes the importance of the pore structure of polyolefin separators but also provides an economic and effective method to screen favorable separators for LMBs.

25 ENERGY STORAGE↗

Searches for neutrinos from cosmic-ray interactions in the Sun using seven years of IceCube data

Cosmic-ray interactions with the solar atmosphere are expected to produce particle showers which in turn produce neutrinos from weak decays of mesons. These solar atmospheric neutrinos (SAνs) have never been observed experimentally. A detection would be an important step in understanding cosmic-ray propagation in the inner solar system and the dynamics of solar magnetic fields. SAνs also represent an irreducible background to solar dark matter searches and a detection would allow precise characterization of this background. In this work, we present the first experimental search based on seven years of data collected from May 2010 to May 2017 in the austral winter with the IceCube Neutrino Observatory. An unbinned likelihood analysis is performed for events reconstructed within 5 degrees of the center of the Sun. No evidence for a SAν flux is observed. After inclusion of systematic uncertainties, we set a 90% upper limit of 1.02 +0.20 -0.18 ·10 -13 GeV-1cm -2 s -1 at 1 TeV.

79 ASTRONOMY AND ASTROPHYSICS↗

Light detection and cosmic rejection in the ICARUS LArTPC at Fermilab

Abstract The ICARUS-T600 detector is a 760-ton Liquid Argon Time Projection Chamber (LArTPC) currently operating at Fermilab as the Far Detector in the Short Baseline Neutrino (SBN) program. The SBN program is composed of three LArTPCs with a central goal of testing the sterile neutrino hypothesis. After operating for 3-years in the Gran Sasso Underground Laboratory, the ICARUS detector was shipped to CERN where it was outfitted with 360 8” Photomultiplier Tubes (PMTs) for a new optical detection system. The PMT system detects fast scintillation light from charged particles interacting in the Liquid Argon, generating the trigger signal for the full detector and allows 3D reconstruction of events. Now operating at shallow depth, the detector is exposed to a high flux of cosmic rays that can fake neutrino interactions. To mitigate this effect a Cosmic Ray Tagger (CRT) and a 3-meter-thick concrete were installed. Precise timing information from both the PMT and CRT subsystems can help to identify whether an interaction originated from inside or outside of the ICARUS cryostat. This paper reviews a method for cosmogenic background reduction and timing calibration of the CRT and PMT light detection systems in ICARUS.

Instruments & Instrumentation↗

High-Order Semi-Discrete Central-Upwind Schemes for Multi-Dimensional Hamilton-Jacobi Equations

We present the first fifth order, semi-discrete central upwind method for approximating solutions of multi-dimensional Hamilton-Jacobi equations. Unlike most of the commonly used high order upwind schemes, our scheme is formulated as a Godunov-type scheme. The scheme is based on the fluxes of Kurganov-Tadmor and Kurganov-Tadmor-Petrova, and is derived for an arbitrary number of space dimensions. A theorem establishing the monotonicity of these fluxes is provided. The spacial discretization is based on a weighted essentially non-oscillatory reconstruction of the derivative. The accuracy and stability properties of our scheme are demonstrated in a variety of examples. A comparison between our method and other fifth-order schemes for Hamilton-Jacobi equations shows that our method exhibits smaller errors without any increase in the complexity of the computations.

Bryson, Steve↗

Convolutional neural network based non-iterative reconstruction for accelerating neutron tomography *

Abstract Neutron computed tomography (NCT), a 3D non-destructive characterization technique, is carried out at nuclear reactor or spallation neutron source-based user facilities. Because neutrons are not severely attenuated by heavy elements and are sensitive to light elements like hydrogen, neutron radiography and computed tomography offer a complementary contrast to x-ray CT conducted at a synchrotron user facility. However, compared to synchrotron x-ray CT, the acquisition time for an NCT scan can be orders of magnitude higher due to lower source flux, low detector efficiency and the need to collect a large number of projection images for a high-quality reconstruction when using conventional algorithms. As a result of the long scan times for NCT, the number and type of experiments that can be conducted at a user facility is severely restricted. Recently, several deep convolutional neural network (DCNN) based algorithms have been introduced in the context of accelerating CT scans that can enable high quality reconstructions from sparse-view data. In this paper, we introduce DCNN algorithms to obtain high-quality reconstructions from sparse-view and low signal-to-noise ratio NCT data-sets thereby enabling accelerated scans. Our method is based on the supervised learning strategy of training a DCNN to map a low-quality reconstruction from sparse-view data to a higher quality reconstruction. Specifically, we evaluate the performance of two popular DCNN architectures—one based on using patches for training and the other on using the full images for training. We observe that both the DCNN architectures offer improvements in performance over classical multi-layer perceptron as well as conventional CT reconstruction algorithms. Our results illustrate that the DCNN can be a powerful tool to obtain high-quality NCT reconstructions from sparse-view data thereby enabling accelerated NCT scans for increasing user-facility throughput or enabling high-resolution time-resolved NCT scans.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Timing Analysis with INTEGRAL: Comparing Different Reconstruction Algorithms

INTEGRAL is one of the few instruments capable of detecting X-rays above 20keV. It is therefore in principle well suited for studying X-ray variability in this regime. Because INTEGRAL uses coded mask instruments for imaging, the reconstruction of light curves of X-ray sources is highly non-trivial. We present results from the comparison of two commonly employed algorithms, which primarily measure flux from mask deconvolution (ii-lc-extract) and from calculating the pixel illuminated fraction (ii-light). Both methods agree well for timescales above about 10 s, the highest time resolution for which image reconstruction is possible. For higher time resolution, ii-light produces meaningful results, although the overall variance of the lightcurves is not preserved.

Grinberg, V.↗

A substandard candle: the low-$$\nu $$ method at few-GeV neutrino energies

Abstract As accelerator-based neutrino oscillation experiments improve oscillation parameter constraints with more data, control over systematic uncertainties on the incoming neutrino flux and interaction models is increasingly important. The intense beams offered by modern experiments permit a variety of options to constrain the flux using in situ “standard candle” measurements. These standard candles must use very well understood interaction processes to avoid introducing additional interaction model dependence. One option often discussed in this context is the “low- $$\nu $$ ν ” method, which is designed to isolate neutrino interactions where there is low energy-transfer to the nucleus, such that the interaction cross section is expected to be approximately constant as a function of neutrino energy. The shape of the low-energy transfer event sample can then be used to extract the flux shape. Applications of the method at high neutrino energies (many tens of GeV) are well understood. However, the applicability of the method at the lower energies of current and future few-GeV accelerator neutrino experiments remains unclear due to the presence of nuclear and form-factor effects inherent in the interaction models.In this analysis we examine the prospects for improving constraints on the accelerator neutrino fluxes in situ with the low- $$\nu $$ ν method in an experiment-independent way, using (anti)neutrino interactions on argon and hydrocarbon targets from the GENIE, NEUT, NuWro and GiBUU event generators. We begin by investigating the extent to which deviations from the constant cross-section assumption are dependent on poorly understood aspects of the neutrino interaction model. We then assess whether a low energy-transfer event sample can be confidently identified using experimentally accessible observables. We finally consider how the practicalities of reconstructing the energy spectrum of interacting neutrinos in realistic detectors might further limit the utility of low- $$\nu $$ ν flux constraints. The results show that flux constraints from the low- $$\nu $$ ν method would be severely dependent on the interaction model assumptions used in an analysis of neutrinos with energies below 5 GeV, and anti-neutrinos below at least 15 GeV. The spread of model predictions show that a low- $$\nu $$ ν analysis is unlikely to offer much improvement on typical neutrino flux uncertainties, even with a perfect detector. Notably—running counter to the assumption inherent to the low- $$\nu $$ ν method—the model-dependence increases with decreasing energy transfer for experiments in the few-GeV region.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Performance of the Pair Spectrometer in Hall D at Jefferson Lab

This article describes the performance of the pair spectrometer installed in experimental Hall D at Jefferson Lab and its successfully operation in multiple experiments with the GlueX detector. The primary purpose of the pair spectrometer is the precise determination of the flux of beam photons incident on the GlueX target, a critical input for physics analyses such as absolute cross-section measurements. The photon energy spectrum is determined by reconstructing electron–positron pairs produced in a thin converter inserted into the photon beam. The spectrometer is integrated into the GlueX trigger system, whose electronics enable continuous real-time monitoring of the photon flux and the recording of e+e- pair candidates for offline analysis. In addition, the pair spectrometer provides a versatile test facility for evaluating calorimeter prototypes using leptons with well-defined energies produced via electromagnetic pair production.

Somov, Alexander [Jefferson Lab]↗

ANTARES Constrains a Blazar Origin of Two IceCube PeV Neutrino Events

Context. The source(s) of the neutrino excess reported by the IceCube Collaboration is unknown. The TANAMI Collaboration recently reported on the multiwavelength emission of six bright, variable blazars which are positionally coincident with two of the most energetic IceCube events. Objects like these are prime candidates to be the source of the highest-energy cosmic rays, and thus of associated neutrino emission. Aims. We present an analysis of neutrino emission from the six blazars using observations with the ANTARES neutrino telescope. Methods. The standard methods of the ANTARES candidate list search are applied to six years of data to search for an excess of muons - and hence their neutrino progenitors - from the directions of the six blazars described by the TANAMI Collaboration, and which are possibly associated with two IceCube events. Monte Carlo simulations of the detector response to both signal and background particle fluxes are used to estimate the sensitivity of this analysis for di erent possible source neutrino spectra. A maximum-likelihood approach, using the reconstructed energies and arrival directions of through-going muons, is used to identify events with properties consistent with a blazar origin. Results. Both blazars predicted to be the most neutrino-bright in the TANAMI sample (1653-329 and 1714-336) have a signal flux fitted by the likelihood analysis corresponding to approximately one event. This observation is consistent with the blazar-origin hypothesis of the IceCube event IC 14 for a broad range of blazar spectra, although an atmospheric origin cannot be excluded. No ANTARES events are observed from any of the other four blazars, including the three associated with IceCube event IC20. This excludes at a 90% confidence level the possibility that this event was produced by these blazars unless the neutrino spectrum is flatter than -2.4.

general↗

Cosmic Ray Background Rejection with Wire-Cell LArTPC Event Reconstruction in the MicroBooNE Detector

For a large liquid-argon time-projection chamber (LArTPC) operating on or near the Earth’s surface to detect neutrino interactions, the rejection of cosmogenic background is a critical and challenging task because of the large cosmic-ray flux and the long drift time of the time-projection chamber. We introduce a superior cosmic background rejection procedure based on the Wire-Cell three-dimensional (3D) event reconstruction for LArTPCs. From an initial 1:20 000 neutrino to cosmic-ray background ratio, we demonstrate these tools on data from the MicroBooNE experiment and create a high-performance generic neutrino event selection with a cosmic contamination of 14.9% (9.7%) for a visible energy region greater than O(200) MeV. The neutrino interaction selection efficiency is 80.4% and 87.6% for inclusive ν μ charged-current and ν e charged-current interactions, respectively. Here, this significantly improved performance compared with existing reconstruction algorithms marks a major milestone toward reaching the scientific goals of LArTPC neutrino oscillation experiments operating near the Earth’s surface.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Performance of high-order Godunov-type methods in simulations of astrophysical low Mach number flows

High-order Godunov methods for gas dynamics have become a standard tool for simulating different classes of astrophysical flows. Their accuracy is mostly determined by the spatial interpolant used to reconstruct the pair of Riemann states at cell interfaces and by the Riemann solver that computes the interface fluxes. In most Godunov-type methods, these two steps can be treated independently, so that many different schemes can in principle be built from the same numerical framework. Because astrophysical simulations often test out the limits of what is feasible with the computational resources available, it is essential to find the scheme that produces the numerical solution with the desired accuracy at the lowest computational cost. However, establishing the best combination of numerical options in a Godunov-type method to be used for simulating a complex hydrodynamic problem is a nontrivial task. In fact, formally more accurate schemes do not always outperform simpler and more diffusive methods, especially if sharp gradients are present in the flow. For this work, we used our fully compressible Seven-League Hydro (SLH) code to test the accuracy of six reconstruction methods and three approximate Riemann solvers on two- and three-dimensional (2D and 3D) problems involving subsonic flows only. We considered Mach numbers in the range from 10 −3 to 10 −1 , which are characteristic of many stellar and geophysical flows. In particular, we considered a well-posed, 2D, Kelvin–Helmholtz instability problem and a 3D turbulent convection zone that excites internal gravity waves in an overlying stable layer. Although the different combinations of numerical methods converge to the same solution with increasing grid resolution for most of the quantities analyzed here, we find that (i) there is a spread of almost four orders of magnitude in computational cost per fixed accuracy between the methods tested in this study, with the most performant method being a combination of a low-dissipation Riemann solver and a sextic reconstruction scheme; (ii) the low-dissipation solver always outperforms conventional Riemann solvers on a fixed grid when the reconstruction scheme is kept the same; (iii) in simulations of turbulent flows, increasing the order of spatial reconstruction reduces the characteristic dissipation length scale achieved on a given grid even if the overall scheme is only second order accurate; (iv) reconstruction methods based on slope-limiting techniques tend to generate artificial, high-frequency acoustic waves during the evolution of the flow; and (v) unlimited reconstruction methods introduce oscillations in the thermal stratification near the convective boundary, where the entropy gradient is steep.

79 ASTRONOMY AND ASTROPHYSICS↗

First experimental study of multiple orientation muon tomography, with image optimization in sparse data environments

Due to the high penetrating power of cosmic ray muons, they can be used to probe very thick and dense objects. As charged particles, they can be tracked by ionization detectors, determining the position and direction of the muons. With detectors on either side of an object, particle direction changes can be used to extract scattering information within an object. This can be used to produce a scattering intensity image within the object related to density and atomic number. Such imaging is typically performed with a single detector-object orientation, taking advantage of the more intense downward flux of muons, producing planar imaging with some depth-of-field information in the third dimension. Several simulation studies have been published with multi-orientation tomography, which can form a three-dimensional representation faster than a single orientation view. In this work we present the first experimental multiple orientation muon tomography study. Experimental muon-scatter based tomography was performed using a concrete filled steel drum with several different metal wedges inside, between detector planes. Data was collected from different detector-object orientations by rotating the steel drum. The data collected from each orientation were then combined using two different tomographic methods. Results showed that using a combination of multiple depth-of-field reconstructions, rather than a traditional inverse Radon transform approach used for CT, resulted in more useful images for sparser data. As cosmic ray muon flux imaging is rate limited, the imaging techniques were compared for sparse data. Using the combined depth-of-field reconstruction technique, fewer detector-object orientations were needed to reconstruct images that could be used to differentiate the metal wedge compositions.

Applied Physics (physics.app-ph)↗

Validating the simulation of beam-ion charge exchange in MAST Upgrade

Abstract Simulation of the impact of charge-exchange (CX) reactions on beam ions in the Mega Amp Spherical Tokamak (MAST) Upgrade was compared to measurements carried out with a fission chamber (neutron fluxes) and a fast ion deuterium-alpha (FIDA) diagnostic. A simple model was developed to reconstruct the outer-midplane neutral density based on measurements of deuterium-alpha emission from edge neutrals, and on Thomson scattering measurements of electron density and temperature. The main computational tools used were the ASCOT orbit-following code and the FIDASIM code for producing synthetic FIDA signals. The neutral density reconstruction agrees qualitatively with SOLPS-ITER modelling and yields a synthetic passive FIDA signal that is consistent with measurement. When CX losses of beam ions are accounted for, predicted neutron emission rates are quantitatively more consistent with measurement. It was necessary to account for CX losses of beam ions in simulations to reproduce the measured passive FIDA signal quantitatively and qualitatively. The results suggest that the neutral density reconstruction is a good approximation, that CX with edge neutrals causes significant beam-ion losses in MAST Upgrade, typically 20% of beam power, and that the ASCOT fast-ion CX model can be used to accurately predict the redistribution and loss of beam ions due to CX.

Physics↗

Constructing angular distributions of neutrinos in core-collapse supernovae from zeroth and first moments calibrated by full Boltzmann neutrino transport

Two-moment neutrino transport methods have been widely used for developing theoretical models of core-collapse supernovae (CCSN), since they substantially reduce the computational burden inherent in the multidimensional neutrino-radiation hydrodynamical simulations. The approximation, however, comes at a price; the detailed structure of angular distribution of neutrinos is sacrificed, that is the main drawback of this approach. In this paper, we develop a novel method by which to construct angular distributions of neutrinos from the zeroth and first angular moments. In our method, the angular distribution is expressed with two quadratic functions of the neutrino angle in a piecewise fashion. We determine the best parameters in the fitting function by comparing to the neutrino data in a spherically symmetric CCSN model with full Boltzmann neutrino transport. We demonstrate the capability of our method by using our recent 2D CCSN model. In this work, we find that the essential features of the angular distributions can be well reconstructed, whereas the angular distributions of incoming neutrinos tend to have large errors that increase with flux factor ($\textit{κ}$). This issue originates from the insensitiveness of incoming neutrinos to $\textit{κ}$, that is an intrinsic limitation in moment methods. Based on the results of the demonstration, we assess the reliability of electron- neutrinos lepton number (ELN)-crossing searches with the two-moment neutrino transport. This analysis is complementary to our [L. Johns and H. Nagakura arXiv:2104.04106] that scrutinizes the limitation of crossing searches with a few moments. We find that the systematic errors of angular distributions for incoming neutrinos lead to misjudgements of the crossing at $\textit{κ}$ ≳ 0.5. This casts doubt on the results of ELN-crossing searches based on two-moment methods in some previous studies.

79 ASTRONOMY AND ASTROPHYSICS↗

Chasing ghosts: characterization of artifact generation in coded aperture decoding due to experimental implementation

Coded aperture imaging is a form of lensless aperture imaging that projects multiple overlapping images of the source onto the detector, enhancing signal strength, which is advantageous for low-flux sources or high-resolution imaging. This technique requires decoding of the detector signal to reconstruct the original source, which involves convolving the detector data with the aperture pattern. When the signal is from a centered point source, the reconstructed source image is known as the point spread function (PSF). A clean PSF without artifacts is a Dirac delta function [Appl. Opt. 20, 1858 (1981)]. This paper examines the robustness of the decoding process against variations in experimental tolerances by analyzing artifact growth in the reconstructed PSF. We illustrate the effects of incorrect magnification, rotation, and detector size and find that aperture–detector rotational misalignment about the imaging axis is the most sensitive parameter, with significant artifact generation occurring with angular offsets of less than one degree. We discuss compensation methods for imperfect aperture placement, finding that small detector sizes produce uncompensatable artifact generation, and compare theoretical predictions with experimental PSF measurements of a rank , 6.8 mm thick (less than one mean free path) coded aperture with a 3.5 mm cell size, conducted at the MegaJOuLe Neutron Imaging Radiography dense plasma focus [IEEE Trans. Plasma Sci. 49, 3299 (2021)] using a 2.45 MeV neutron source. Based on our findings, we recommend using magnified coded apertures in the under-sampled regime, which allows for the inclusion of fiducial markers to characterize aperture–detector rotational offsets and the addition of mechanical coupling, where possible, to constrain rotational and magnification offsets.

Selwood, M. P. [Lawrence Livermore National Labora↗