Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “distributed inference”

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 775 records · Page 43

Neutrino Beam Monitoring

Accelerator facilities produce neutrino beams from meson decays in a decay volume. Experiments measure event rates that depend on flux, cross sections, and detector response, so the flux is predicted using hadron production and beamline modeling and constrained by beam instrumentation, since near detectors alone cannot separate flux from cross section. Proton, hadron, and muon monitors can track the parent particle distributions and beam conditions, providing the inputs needed for flux predictions in long-baseline experiments such as NOvA, T2K, and DUNE. This talk reviews how beam monitors are used in practice to understand neutrino flux. Proton beam monitors tell where the beam hits the target and how stable it is. Farther downstream, hadron and muon monitors sample particles produced in meson decays. Because those muons come from the same parents as the neutrinos, their profiles reveal focusing, alignment shifts, and other changes in the beam, and they are routinely used to detect problems and guide flux predictions. The muon information can be used more quantitatively; for example, to infer the parent meson phase space, and fast radiation-hard timing detectors can add sensitivity to the momentum dependence of the focusing. These developments show both how tightly beam measurements can constrain the flux and where the current limits still lie. These approaches complement monitored-beam concepts, in which the decay region is instrumented to detect charged leptons from meson decays and to measure the neutrino flux directly.

Ganguly, Sudeshna [Fermilab] (ORCID:00000003163482↗

Cosmology with persistent homology: a Fisher forecast

Abstract Persistent homology naturally addresses the multi-scale topological characteristics of the large-scale structure as a distribution of clusters, loops, and voids. We apply this tool to the dark matter halo catalogs from theQuijotesimulations, and build a summary statistic for comparison with the joint power spectrum and bispectrum statistic regarding their information content on cosmological parameters and primordial non-Gaussianity. Through a Fisher analysis, we find that constraints from persistent homology are tighter for 8 out of the 10 parameters by margins of 13–50%. The complementarity of the two statistics breaks parameter degeneracies, allowing for a further gain in constraining power when combined. We run a series of consistency checks to consolidate our results, and conclude that our findings motivate incorporating persistent homology into inference pipelines for cosmological survey data.

Astronomy & Astrophysics↗

Dark matter haloes of massive elliptical galaxies at z ∼ 0.2 are well described by the Navarro–Frenk–White profile

ABSTRACT We investigate the internal structure of elliptical galaxies at z ∼ 0.2 from a joint lensing–dynamics analysis. We model Hubble Space Telescope images of a sample of 23 galaxy–galaxy lenses selected from the Sloan Lens ACS (SLACS) survey. Whereas the original SLACS analysis estimated the logarithmic slopes by combining the kinematics with the imaging data, we estimate the logarithmic slopes only from the imaging data. We find that the distribution of the lensing-only logarithmic slopes has a median 2.08c ± 0.03 and intrinsic scatter 0.13 ± 0.02, consistent with the original SLACS analysis. We combine the lensing constraints with the stellar kinematics and weak lensing measurements, and constrain the amount of adiabatic contraction in the dark matter (DM) haloes. We find that the DM haloes are well described by a standard Navarro–Frenk–White halo with no contraction on average for both of a constant stellar mass-to-light ratio (M/L) model and a stellar M/L gradient model. For the M/L gradient model, we find that most galaxies are consistent with no M/L gradient. Comparison of our inferred stellar masses with those obtained from the stellar population synthesis method supports a heavy initial mass function (IMF) such as the Salpeter IMF. We discuss our results in the context of previous observations and simulations, and argue that our result is consistent with a scenario in which active galactic nucleus feedback counteracts the baryonic-cooling-driven contraction in the DM haloes.

Shajib, Anowar J.↗

Functional and Evolutionary Integration of a Fungal Gene With a Bacterial Operon

Abstract Siderophores are crucial for iron-scavenging in microorganisms. While many yeasts can uptake siderophores produced by other organisms, they are typically unable to synthesize siderophores themselves. In contrast, Wickerhamiella/Starmerella (W/S) clade yeasts gained the capacity to make the siderophore enterobactin following the remarkable horizontal acquisition of a bacterial operon enabling enterobactin synthesis. Yet, how these yeasts absorb the iron bound by enterobactin remains unresolved. Here, we demonstrate that Enb1 is the key enterobactin importer in the W/S-clade species Starmerella bombicola. Through phylogenomic analyses, we show that ENB1 is present in all W/S clade yeast species that retained the enterobactin biosynthetic genes. Conversely, it is absent in species that lost the ent genes, except for Starmerella stellata, making this species the only cheater in the W/S clade that can utilize enterobactin without producing it. Through phylogenetic analyses, we infer that ENB1 is a fungal gene that likely existed in the W/S clade prior to the acquisition of the ent genes and subsequently experienced multiple gene losses and duplications. Through phylogenetic topology tests, we show that ENB1 likely underwent horizontal gene transfer from an ancient W/S clade yeast to the order Saccharomycetales, which includes the model yeast Saccharomyces cerevisiae, followed by extensive secondary losses. Taken together, these results suggest that the fungal ENB1 and bacterial ent genes were cooperatively integrated into a functional unit within the W/S clade that enabled adaptation to iron-limited environments. This integrated fungal-bacterial circuit and its dynamic evolution determine the extant distribution of yeast enterobactin producers and cheaters.

59 BASIC BIOLOGICAL SCIENCES↗

Physics-informed neural networks for solving forward and inverse flow problems via the Boltzmann-BGK formulation

We report the Boltzmann equation with the Bhatnagar-Gross-Krook collision model (Boltzmann-BGK equation) has been employed to describe multiscale flows, i.e., from the hydrodynamic limit to free molecular flow. In this study we employ physics-informed neural networks (PINNs) to solve forward and inverse problems via the Boltzmann-BGK formulation (PINN-BGK), enabling PINNs to model flows in both the continuum and rarefied regimes. In particular, the PINN-BGK is composed of three sub-networks, i.e., the first for approximating the equilibrium distribution function, the second for approximating the non-equilibrium distribution function, and the third one for encoding the Boltzmann-BGK equation as well as the corresponding boundary/initial conditions. By minimizing the residuals of the governing equations and the mismatch between the predicted and provided boundary/initial conditions, we can approximate the Boltzmann-BGK equation for both continuous and rarefied flows. For forward problems, the PINN-BGK is utilized to solve various benchmark flows given boundary/initial conditions, e.g. Kovasznay flow, Taylor-Green flow, cavity flow, and micro Couette flow for Knudsen number up to 5. For inverse problems, we focus on rarefied flows in which accurate boundary conditions are difficult to obtain. We employ the PINN-BGK to infer the flow field in the entire computational domain given a limited number of interior scattered measurements on the velocity without using the (unknown) boundary conditions. Results for the two-dimensional micro Couette and micro cavity flows with Knudsen numbers ranging from 0.1 to 10 indicate that the PINN-BGK can infer the velocity field in the entire domain with good accuracy. Finally, we also present some results on using transfer learning to accelerate the training process. Specifically, we can obtain a three-fold speedup comparted to the standard training process (e.g., Adam plus L-BFG-B) for the two-dimensional flow problems considered in our work.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

An acoustic investigation of the near-surface turbulence on Mars

The Perseverance rover is carrying out an original acoustic experiment on Mars: the SuperCam microphone records the spherical acoustic waves generated by laser sparks at distances from 2 m to more than 8 m. These N-shaped acoustic waves scatter from the multiple local heterogeneities of the turbulent atmosphere. Therefore, large and random fluctuations of sound travel time and intensity develop as the waves cross the medium. The variances of the travel times and the scintillation index (normalized variance of the sound intensity) are studied within the mathematical formalism of the propagation of spherical acoustic waves through thermal turbulence to infer statistical properties of the Mars atmospheric temperature fluctuation field. The comparison with the theory is made by simplifying assumptions that do not include wind fluctuations and diffraction effects. Two Earth years (about one Martian year) of observations acquired during the maximum convective period (10:00–14:00 Mars local time) show a good agreement between the dataset and the formalism: the travel time variance diverges from the linear Chernov solution exactly where the density of occurrence of the first caustic reaches its maximum. Moreover, on average, waves travel faster than the mean speed of sound due to a fast path effect, which is also observed on Earth. To account for the distribution of turbulent eddies, several power spectra are tested and the best match to observation is obtained with a generalized von Karman spectrum with a shallower slope than the Kolmogorov cascade, φ(k)∝(1 + k 2 L 2 ) -4/3 . It is associated with an outer scale of turbulence, L, of 11 cm at 2 m above the surface and a standard deviation of 6 K over 9 s for the temperature. These near-surface atmospheric properties are consistent with a weak to moderate wave scattering regime around noon with little saturation. Overall, this study presents an innovative and promising methodology to probe the near-surface atmospheric turbulence on Mars.

79 ASTRONOMY AND ASTROPHYSICS↗

Estimating β eff from static, late-time Feynman neutron noise measurements

The effective delayed neutron fraction, β eff , is a kinematic parameter describing the contribution of delayed neutrons to the effective multiplication factor, $\kappa$ eff , and is of great interest to the nuclear criticality safety community and beyond. In fact, $\kappa$ eff presently must be inferred (as opposed to directly measured) since variables like β eff must be obtained from simulation, reference tables, or a dynamic or pulsed measurement where the assembly of fissile material is perturbed and not in steady state. Several recent works have proposed neutron noise techniques to estimate β eff from static measurements (such as benchmark experiments). This work experimentally investigates an approach based on looking at late-time (up to thousands of seconds) correlated neutrons in Feynman histograms. This work utilizes measurements of weapons grade plutonium and highly enriched uranium to construct late-time Feynman histograms and investigate a proposed methodology of estimating β eff from a fit of the distribution. Data is shown from both He-3 detectors and an organic scintillator array. This work is preliminary in nature, demonstrating some experimental feasibility and proposing improvements for future investigations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Halo effective field theory analysis of one-neutron knockout reactions of Be 11 and C 15

Background: One-nucleon knockout reactions provide insightful information on the single-particle structure of nuclei. When applied to one-neutron halo nuclei, they are purely peripheral, suggesting that they could be properly modeled by describing the projectile within a halo effective field theory (halo-EFT). Purpose: We reanalyze the one-neutron knockout measurements of 11 Be and 15 C —both one-neutron halo nuclei—on beryllium at about 60 MeV/nucleon. We consider halo-EFT descriptions of these nuclei which already provide excellent agreement with breakup and transfer data. Method: Here we include a halo-EFT description of the projectile within an eikonal-based model of the reaction and compare its outcome to existing data. Results: Excellent agreement with experiment is found for both nuclei. The asymptotic normalization coefficients inferred from this comparison confirm predictions from ab initio nuclear-structure calculations and values deduced from transfer data. Conclusions: Halo-EFT can be reliably used to analyze one-neutron knockout reactions measured for halo nuclei and test predictions from state-of-the-art nuclear structure models on these experimental data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Using neutrons and x rays to measure plasma conditions in a solid sphere of deuterated polyethylene compressed to densities of 35 g/cc at temperatures of 2 keV and pressures of 40 Gbar

This paper describes an experiment that shock compresses the center of a solid deuterated polyethylene sphere, CD 2 , to densities of 35 g/cc and temperatures of 2 keV with corresponding pressure of 40 Gbar. The design employs a strong spherically converging shock launched through a solid ball of material using a Hohlraum radiation drive. As the shock coalesces at the center it produces a hot spot that we characterize by measuring the x-ray self-emission and 2.45 MeV neutrons emitted. Two-dimensional images and time-resolved measurements of the x rays emitted determine the size and time duration of the hot spot, leading to an estimated 2k eV electron temperature. The neutron time of flight spectrometer measures an average ion temperature of 1.06 ± 0.15 keV and neutron yield of 7.0 (±0.5) × 10 9 DD neutrons. Our new distribution function tool enables us to create a forward model of the experimental data based on 1D radiation-hydrodynamic simulations, leading to a better understanding of the plasma conditions that produce the measured neutrons and x rays. Furthermore, our simulations indicate that the x rays are produced in a short-lived hot-dense core over tens of picoseconds, whereas the neutron emission continues for about 200 ps, as the hot core starts to expand, thereby leading to a lower mean temperature of the plasma during neutron production. This finding is in agreement with the experimental data, and we therefore conclude that the forward-modeling is a useful tool for inferring the conditions of the hot spot in a laser-driven implosion during burn.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Search for photons above 10 19 eV with the surface detector of the Pierre Auger Observatory

We use the surface detector of the Pierre Auger Observatory to search for air showers initiated by photons with an energy above 10 19 eV. Photons in the zenith angle range from 30° to 60° can be identified in the overwhelming background of showers initiated by charged cosmic rays through the broader time structure of the signals induced in the water-Cherenkov detectors of the array and the steeper lateral distribution of shower particles reaching ground. Applying the search method to data collected between January 2004 and June 2020, upper limits at 95% CL are set to an E -2 diffuse flux of ultra-high energy photons above 10 19 eV, 2 × 10 19 eV and 4 × 10 19 eV amounting to 2.11 × 10 -3 , 3.12 × 10 -4 and 1.72 × 10 -4 km -2 sr -1 yr -1 , respectively. While the sensitivity of the present search around 2 × 10 19 eV approaches expectations of cosmogenic photon fluxes in the case of a pure-proton composition, it is one order of magnitude above those from more realistic mixed-composition models. The inferred limits have also implications for the search of super-heavy dark matter that are discussed and illustrated.

ultra high energy photons and neutrinos↗

Accurate and Rapid Forecasts for Geologic Carbon Storage via Learning-Based Inversion-Free Prediction

Carbon capture and storage (CCS) is one approach being studied by the U.S. Department of Energy to help mitigate global warming. The process involves capturing CO 2 emissions from industrial sources and permanently storing them in deep geologic formations (storage reservoirs). However, CCS projects generally target “green field sites,” where there is often little characterization data and therefore large uncertainty about the petrophysical properties and other geologic attributes of the storage reservoir. Consequently, ensemble-based approaches are often used to forecast multiple realizations prior to CO 2 injection to visualize a range of potential outcomes. In addition, monitoring data during injection operations are used to update the pre-injection forecasts and thereby improve agreement between forecasted and observed behavior. Thus, a system for generating accurate, timely forecasts of pressure buildup and CO 2 movement and distribution within the storage reservoir and for updating those forecasts via monitoring measurements becomes crucial. This study proposes a learning-based prediction method that can accurately and rapidly forecast spatial distribution of CO 2 concentration and pressure with uncertainty quantification without relying on traditional inverse modeling. The machine learning techniques include dimension reduction, multivariate data analysis, and Bayesian learning. The outcome is expected to provide CO 2 storage site operators with an effective tool for timely and informative decision making based on limited simulation and monitoring data.

58 GEOSCIENCES↗

Bayesian uncertainty quantification for nuclear matter incompressibility

Within a Bayesian statistical framework using the standard Skyrme-Hartree-Fock model, the maximum a posteriori (MAP) values and uncertainties of nuclear matter incompressibility and isovector interaction parameters are inferred from the experimental data of giant resonances and neutron-skin thicknesses of typical heavy nuclei. Furthermore, with the uncertainties of the isovector interaction parameters constrained by the data of the isovector giant dipole resonance and the neutron-skin thickness, we have obtained K 0 = 223$^{+7}_{–8}$ MeV at 68% confidence level using the data of the isoscalar giant monopole resonance in 208 Pb measured at the Research Center for Nuclear Physics (RCNP), Japan, and at the Texas A&M University, USA. Although the corresponding 120 Sn data gives a MAP value for K 0 about 5 MeV smaller than the 208 Pb data, there are significant overlaps in their posterior probability distribution functions.

190 ≤ A ≤ 219↗

Exploring Li-Ion Transport Properties of Li 3 TiCl 6 : A Machine Learning Molecular Dynamics Study

We performed large-scale molecular dynamics simulations based on a machine-learning force field (MLFF) to investigate the Li-ion transport mechanism in cation-disordered Li 3 TiCl 6 cathode at six different temperatures, ranging from 25°C to 100°C. In this work, deep neural network method and data generated by ab − initio molecular dynamics (AIMD) simulations were deployed to build a high-fidelity MLFF. Radial distribution functions, Li-ion mean square displacements (MSD), diffusion coefficients, ionic conductivity, activation energy, and crystallographic direction-dependent migration barriers were calculated and compared with corresponding AIMD and experimental data to benchmark the accuracy of the MLFF. From MSD analysis, we captured both the self and distinct parts of Li-ion dynamics. The latter reveals that the Li-ions are involved in anti-correlation motion that was rarely reported for solid-state materials. Similarly, the self and distinct parts of Li-ion dynamics were used to determine Haven’s ratio to describe the Li-ion transport mechanism in Li 3 TiCl 6 . Obtained trajectory from molecular dynamics infers that the Li-ion transportation is mainly through interstitial hopping which was confirmed by intra- and inter-layer Li-ion displacement with respect to simulation time. Ionic conductivity (1.06 mS/cm) and activation energy (0.29eV) calculated by our simulation are highly comparable with that of experimental values. Overall, the combination of machine-learning methods and AIMD simulations explains the intricate electrochemical properties of the Li 3 TiCl 6 cathode with remarkably reduced computational time. Thus, our work strongly suggests that the deep neural network-based MLFF could be a promising method for large-scale complex materials.

Selvaraj, Selva Chandrasekaran (ORCID:000000029023↗

General Bayesian Inference over the Stiefel Manifold via the Givens Transform

We introduce here an approach based on the Givens representation for posterior inference in statistical models with orthogonal matrix parameters, such as factor models and probabilistic principal component analysis (PPCA). We show how the Givens representation can be used to develop practical methods for transforming densities over the Stiefel manifold into densities over subsets of Euclidean space. We demonstrate how to deal with issues arising from the topology of the Stiefel manifold and how to inexpensively compute the change-of-measure terms. We introduce an auxiliary parameter approach that limits the impact of topological issues. We provide both analysis of our methods and numerical examples demonstrating the effectiveness of the approach. We also discuss how our Givens representation can be used to define general classes of distributions over the space of orthogonal matrices. We then give demonstrations on several examples showing how the Givens approach performs in practice in comparison with other methods.

97 MATHEMATICS AND COMPUTING↗

$\mathrm{SageNet}$: Fast Neural Network Emulation of the Stiff-amplified Gravitational Waves from Inflation

Accurate modeling of the inflationary gravitational waves (GWs) requires time-consuming, iterative numerical integrations of differential equations to take into account their backreaction on the expansion history. To improve computational efficiency while preserving accuracy, we present the Stiff-amplified Gravitational-wave Emulator Network (SageNet), a deep learning framework designed to replace conventional numerical solvers (code available at https://github.com/YifangLuo/SageNet). SageNet employs a long short-term memory architecture to emulate the present-day energy density spectrum of the inflationary GWs with possible stiff amplification, Ω GW (f). Trained on a data set of 25,689 numerically generated solutions, SageNet allows accurate reconstructions of Ω GW (f) and generalizes well to a wide range of cosmological parameters; 90.9% of the test emulations with randomly distributed parameters exhibit errors of under 4%. In addition, SageNet demonstrates its ability to learn and reproduce the artificial, adaptive sampling patterns in numerical calculations, which implement denser sampling of frequencies around changes in spectral indices in Ω GW (f). The dual capability of learning both physical and artificial features of the numerical GW spectra establishes SageNet as a robust alternative to exact numerical methods. Finally, our benchmark tests show that SageNet reduces the computation time from tens of seconds to milliseconds, achieving a speedup of ∼10 4 times over standard CPU-based numerical solvers with the potential for further acceleration on GPU hardware. These capabilities make SageNet a powerful tool for accelerating Bayesian inference procedures for extended cosmological models. In a broad sense, the SageNet framework offers a fast, accurate, and generalizable solution to modeling cosmological observables whose theoretical predictions demand costly differential equation solvers.

Astronomy data modeling↗

Goal-oriented a-posteriori estimation of model error as an aid to parameter estimation

In this work, a Bayesian model calibration framework is presented that utilizes goal-oriented a-posterior error estimates in quantities of interest (QoIs) for classes of high-fidelity models characterized by PDEs. It is shown that for a large class of computational models, it is possible to develop a computationally inexpensive procedure for calibrating parameters of high-fidelity models of physical events when the parameters of low-fidelity (surrogate) models are known with acceptable accuracy. The main ingredients in the proposed model calibration scheme are goal-oriented a-posteriori estimates of error in QoIs computed using a so-called lower fidelity model compared to those of an uncalibrated higher fidelity model. The estimates of error in QoIs are used to define likelihood functions in Bayesian inversion analysis. A standard Bayesian approach is employed to compute the posterior distribution of model parameters of high-fidelity models. As applications, parameters in a quasi-linear second-order elliptic boundary-value problem (BVP) are calibrated using a second-order linear elliptic BVP. In a second application, parameters of a tumor growth model involving nonlinear time-dependent PDEs are calibrated using a lower fidelity linear tumor growth model with known parameter values.

A-posterior estimates↗

Responses of the AC/DC Global Electric Circuit to Volcanic Electrical Activity in the Hunga Tonga-Hunga Ha'apai Eruption on 15 January 2022

Responses of the AC and DC global electric circuits (GECs) to the large eruption of the Hunga Tonga-Hunga Ha'apai (HT-HH) volcano on 15 January 2022 are discussed. The AC-related investigation is based on Schumann resonance (SR) measurements from six stations on four continents. The DC-related investigation utilizes atmospheric electric field (potential gradient, PG) measurements from six recording stations in Europe and the USA. According to data from the GLD360 and WWLLN lightning detection networks, the peak lightning stroke rate, 83/s, was dominated by negative polarity lightning, but the distributions of positive and negative lightning discharges in latitude and longitude around the volcano differed. A global intensification of SR is apparent in connection with the enhanced lightning activity caused by the eruption. SR data-based results confirm that the lightning activity in the eruption dominated the naturally occurring global activity for a period of about 1 hr. The highly localized increase in lightning activity over HT-HH was a unique point source of SR excitation. PG measurements suggest that impulse-like charging of the DC GEC, by ~15%, via negative cloud-to-ground lightning strokes took place twice during the eruption. A time constant of 7 or 8 min has been inferred for near-surface PG changes due to these enhancements. This could be the first direct measurement of the time constant of the GEC near the Earth's surface, as well as the first observation of the direct charging of the DC GEC by a unique atmospheric electrified source.

58 GEOSCIENCES↗

Physics-Based Machine Learning Methods for U-235 Forensics Signatures

Signatures of low-intensity U-235 sources have been recently studied by utilizing a variety of machine learning (ML) classifiers using features derived from gamma spectral measurements collectedunder structured campaigns. Several ML classifiers, such as ensemble of tress and classification trees, revealed misleadingly-optimistic training error due to over-fitting, and furthermore,their performance is not directly relatable to the physical properties due to their data-driven, opaque designs. We present a regression-based ML method that first estimates the inverse distanceto the source and then utilizes a threshold to infer its presence, by representing the background as a source located at an infinite distance. For the inverse distance estimation, we study the ensembleof trees and Gaussian process regression methods, and a hyper parameter auto-tuning and selection method that employs five regression estimators. These methods avoid the over-fittingobserved in several ML classifiers, while providing the classification error nearly comparable to them based on independent test data. Their error is directly related to estimates of the inversephysical distance to source, and the precision of error determines the seperability property that determines the false alarm and missed detection rates. The property of monotonic decrease of thesource strength with increasing detector distance combined with Poisson distribution of measurements is utilized to analytically validate these methods by deriving the generalization equations ofunderlying regression methods.

Rao, Nageswara↗