Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “interpretability”

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 361 records · Page 20

Universal and interpretable classification of atomistic structural transitions via unsupervised graph learning

Materials processing often occurs under extreme dynamic conditions leading to a multitude of unique structural environments. These structural environments generally occur at high temperatures and/or high pressures, often under non-equilibrium conditions, which results in drastic changes in the material's structure over time. Computational techniques, such as molecular dynamics simulations, can probe the atomic regime under these extreme conditions. However, characterizing the resulting diverse atomistic structures as a material undergoes extreme changes in its structure has proved challenging due to the inherently non-linear relationship between structures as large-scale changes occur. Here, we introduce SODAS++, a universal graph neural network framework, that can accurately and intuitively quantify the atomistic structural evolution corresponding to the transition between any two arbitrary phases. We showcase SODAS++ for both solid–solid and solid–liquid transitions for systems of increasing geometric and chemical complexity, such as colloidal systems, elemental Al, rutile and amorphous TiO 2 , and the non-stoichiometric ternary alloy Ag 26 Au 5 Cu 19 . Finally, we show that SODAS++ can accurately quantify all transitions in a physically interpretable manner, showcasing the power of unsupervised graph neural network encodings for capturing the complex and non-linear pathway, a material's structure takes as it evolves.

36 MATERIALS SCIENCE↗

Interpretation of autoencoder-learned collective variables using Morse–Smale complex and sublevelset persistent homology: An application on molecular trajectories

Dimensionality reduction often serves as the first step toward a minimalist understanding of physical systems as well as the accelerated simulations of them. In particular, neural network-based nonlinear dimensionality reduction methods, such as autoencoders, have shown promising outcomes in uncovering collective variables (CVs). However, the physical meaning of these CVs remains largely elusive. In this work, we constructed a framework that (1) determines the optimal number of CVs needed to capture the essential molecular motions using an ensemble of hierarchical autoencoders and (2) provides topology-based interpretations to the autoencoder-learned CVs with Morse–Smale complex and sublevelset persistent homology. Furthermore, this approach was exemplified using a series of n-alkanes and can be regarded as a general, explainable nonlinear dimensionality reduction method.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Novel data interpretation method for DIII-D divertor retarding field energy analyzer with 3-D particle-in-cell simulations

A novel data interpretation process that utilizes comprehensive particle-in-cell (PIC) simulations is developed for the new retarding field energy analyzer (RFEA) currently being constructed at DIII-D for the lower divertor using the Divertor Material Evaluation System. Furthermore, this probe is expected to survive a heat load of up to 100 MW/m 2 for up to 5 s and reliably measure the main ion temperature (T i ) on the divertor target ranging from 10 to 200 eV. These extreme conditions posed significant engineering limitations on the probe geometry, thus extensive validation work has been performed. The conventional fitting method for the RFEA I–V characteristics is based on a simplified 1-D model without considering the ion space charge inside the probe cavity and may not be sufficient for probes designed for the DIII-D divertor environment. In this article, a more realistic description of the particle propagation process within the RFEA cavity is achieved by including both 3-D geometric effects and ion space charge in the PIC simulations, and the capability to reconstruct the ion energy distribution functions is demonstrated with reasonable consistency.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

New interpretation of ion cyclotron emission from a tokamak

The tempting interpretation of ion cyclotron emission in terms of compressional Alfvén eigenmodes involving energetic ions is inconsistent with recent TCV experimental observations in some important aspects, such as (i) the perturbed poloidal field is exceeding the parallel perturbed magnetic field significantly, and (ii) the modes are near cyclotron harmonic and exhibit Alfvèn scaling of their frequency. We show that these characteristics can be explained by considering finite Larmor radius effects of thermal ions in shear Alfvén waves that allow such waves to exist well above the ion cyclotron frequency in the form of wave-packets bouncing within the plasma volume.

Plasma waves↗

Stability-driven nonnegative matrix factorization to interpret spatial gene expression and build local gene networks

Spatial gene expression patterns enable the detection of local covariability and are extremely useful for identifying local gene interactions during normal development. The abundance of spatial expression data in recent years has led to the modeling and analysis of regulatory networks. The inherent complexity of such data makes it a challenge to extract biological information. We developed staNMF, a method that combines a scalable implementation of nonnegative matrix factorization (NMF) with a new stability-driven model selection criterion. When applied to a set of Drosophila early embryonic spatial gene expression images, one of the largest datasets of its kind, staNMF identified 21 principal patterns (PP). Providing a compact yet biologically interpretable representation of Drosophila expression patterns, PP are comparable to a fate map generated experimentally by laser ablation and show exceptional promise as a data-driven alternative to manual annotations. Our analysis mapped genes to cell-fate programs and assigned putative biological roles to uncharacterized genes. Finally, we used the PP to generate local transcription factor regulatory networks. Spatially local correlation networks were constructed for six PP that span along the embryonic anterior-posterior axis. Using a two-tail 5% cutoff on correlation, we reproduced 10 of the 11 links in the well-studied gap gene network. In conclusion, the performance of PP with the Drosophila data suggests that staNMF provides informative decompositions and constitutes a useful computational lens through which to extract biological insight from complex and often noisy gene expression data.

59 BASIC BIOLOGICAL SCIENCES↗

An interpretable machine learning framework to understand bikeshare demand before and during the COVID-19 pandemic in New York City

In recent years, bikesharing systems have become increasingly popular as affordable and sustainable micromobility solutions. Advanced mathematical models such as machine learning are required to generate good forecasts for bikeshare demand. Here, this study proposes a machine learning modeling framework to estimate hourly demand in a large-scale bikesharing system. Two Extreme Gradient Boosting models were developed: one using data from before the COVID-19 pandemic (March 2019 to February 2020) and the other using data from during the pandemic (March 2020 to February 2021). Furthermore, a model interpretation framework based on SHapley Additive exPlanations was implemented. Based on the relative importance of the explanatory variables considered in this study, share of female users and hour of day were the two most important explanatory variables in both models. However, the month variable had higher importance in the pandemic model than in the pre-pandemic model.

99 GENERAL AND MISCELLANEOUS↗

Progress Toward Interpretable Machine Learning–Based Disruption Predictors Across Tokamaks

Here in this paper we lay the groundwork for a robust cross-device comparison of data-driven disruption prediction algorithms on DIII-D and JET tokamaks. In order to consistently carry on a comparative analysis, we define physics-based indicators of disruption precursors based on temperature, density, and radiation profiles that are currently not used in many other machine learning predictors for DIII-D data. These profile-based indicators are shown to well-describe impurity accumulation events in both DIII-D and JET discharges that eventually disrupt. The univariate analysis of the features used as input signals in the data-driven algorithms applied on the data of both tokamaks statistically highlights the differences in the dominant disruption precursors. JET with its ITER-like wall is more prone to impurity accumulation events, while DIII-D is more subject to edge-cooling mechanisms that destabilize dangerous magnetohydrodynamic modes. Even though the analyzed data sets are characterized by such intrinsic differences, we show through a few examples that the inclusion of physics-based disruption markers in data-driven algorithms is a promising path toward the realization of a uniform framework to predict and interpret disruptive scenarios across different tokamaks. As long as the destabilizing precursors are diagnosed in a device-independent way, the knowledge that data-driven algorithms learn on one device can be re-used to explain a disruptive behavior on another device.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

THE EFFECT OF POPULATION-WIDE MASS-TO-RADIUS RELATIONSHIPS ON THE INTERPRETATION OF KEPLER AND HARPS SUPER-EARTH OCCURRENCE RATES

The ongoing High Accuracy Radial velocity Planet Search (HARPS) has found that 30%-50% of G and K dwarfs in the solar neighborhood host planets with M{sub pl} {approx}< M{sub Nep} in orbits of P {<=} 50 days. At first glance, this high overall occurrence rate seems at best to be marginally consistent with the planet frequency measured during Q0-Q2 of the Kepler mission, whose 1235 detected planetary candidates naively imply that {approx}15% of main-sequence dwarfs harbor a short-period planet with R{sub pl} < 4 R{sub Circled-Plus }. A rigorous comparison between the two surveys is difficult, however, as they observe different stellar populations, measure different planetary physical properties, and are subject to radically different sampling plans. In this article, we report the results of a Monte Carlo study which seeks to partially overcome this apparent discrepancy by identifying plausible planetary population distributions which can jointly conform to the results of the two surveys. We find that, given the HARPS occurrence rate, either a population subject to a mass-density relationship extrapolated from our solar system or a population concurrently consisting of dense silicate-iron planets and low-density gaseous worlds can produce total numbers of planet candidates consistent with those actually detected by Kepler. However, these two mass-to-radius relationships (M-Rs) resolve the apparent occurrence rate discrepancy with different mass and period distributions, enabling future observations to rule out M-Rs that do not fully describe the observed planet population in a global sense. Extracting information of this nature from the transit-radial-velocity comparison has significant implications for the interpretation of planet occurrence rates: if a multi-valued M-R, which allows planets of similar mass to have significantly different radii, emerges from observational data, then multiple formation mechanisms are likely at work. In this event, HARPS may be detecting a large population of dense low-mass planets, while Kepler detects a large population of gaseous sub-Neptunes.

79 ASTRONOMY AND ASTROPHYSICS↗

LORENTZ-FACTOR-ISOTROPIC-LUMINOSITY/ENERGY CORRELATIONS OF GAMMA-RAY BURSTS AND THEIR INTERPRETATION

The bulk Lorentz factor of the gamma-ray burst (GRB) ejecta ({Gamma}{sub 0}) is a key parameter to understanding GRB physics. Liang et al. have discovered a correlation between {Gamma}{sub 0} and isotropic {gamma}-ray energy: {Gamma}{sub 0}{proportional_to}E{sup 0.25}{sub {gamma},iso,52}. By including more GRBs with updated data and more methods to derive {Gamma}{sub 0}, we confirm this correlation and obtain {Gamma}{sub 0} {approx_equal} 91E{sup 0.29}{sub {gamma},iso,52}. Evaluating the mean isotropic {gamma}-ray luminosities L{sub {gamma},iso} of the GRBs in the same sample, we discover an even tighter correlation {Gamma}{sub 0} {approx_equal} 249L{sup 0.30}{sub {gamma},iso,52}. We propose an interpretation to this later correlation. Invoking a neutrino-cooled hyperaccretion disk around a stellar mass black hole as the central engine of GRBs, we derive jet luminosity powered by neutrino annihilation and baryon loading from a neutrino-driven wind. Applying beaming correction, we finally derive {Gamma}{sub 0}{proportional_to}L{sup 0.22}{sub {gamma},iso}, which is consistent with the data. This suggests that the central engine of long GRBs is likely a stellar mass black hole surrounded by a hyper-accreting disk.

79 ASTRONOMY AND ASTROPHYSICS↗

INTERPRETATION OF INFRARED VIBRATION-ROTATION SPECTRA OF INTERSTELLAR AND CIRCUMSTELLAR MOLECULES

Infrared vibration-rotation lines can be valuable probes of interstellar and circumstellar molecules, especially symmetric molecules, which have no pure rotational transitions. But most such observations have been interpreted with an isothermal absorbing slab model, which leaves out important radiative transfer and molecular excitation effects. A more realistic non-LTE and non-isothermal radiative transfer model has been constructed. The results of this model are in much better agreement with the observations, including cases where lines in one branch of a vibration-rotation band are in absorption and another in emission. In general, conclusions based on the isothermal absorbing slab model can be very misleading, but the assumption of LTE may not lead to such large errors, particularly if the radiation field temperature is close to the gas temperature.

74 ATOMIC AND MOLECULAR PHYSICS↗

INTERPRETING ERUPTIVE BEHAVIOR IN NOAA AR 11158 VIA THE REGION'S MAGNETIC ENERGY AND RELATIVE-HELICITY BUDGETS

In previous works, we introduced a nonlinear force-free method that self-consistently calculates the instantaneous budgets of free magnetic energy and relative magnetic helicity in solar active regions (ARs). Calculation is expedient and practical, using only a single vector magnetogram per computation. We apply this method to a time series of 600 high-cadence vector magnetograms of the eruptive NOAA AR 11158 acquired by the Helioseismic and Magnetic Imager on board the Solar Dynamics Observatory over a five-day observing interval. Besides testing our method extensively, we use it to interpret the dynamical evolution in the AR, including eruptions. We find that the AR builds large budgets of both free magnetic energy and relative magnetic helicity, sufficient to power many more eruptions than the ones it gave within the interval of interest. For each of these major eruptions, we find eruption-related decreases and subsequent free-energy and helicity budgets that are consistent with the observed eruption (flare and coronal mass ejection (CME)) sizes. In addition, we find that (1) evolution in the AR is consistent with the recently proposed (free) energy-(relative) helicity diagram of solar ARs, (2) eruption-related decreases occur before the flare and the projected CME-launch times, suggesting that CME progenitors precede flares, and (3) self terms of free energy and relative helicity most likely originate from respective mutual terms, following a progressive mutual-to-self conversion pattern that most likely stems from magnetic reconnection. This results in the non-ideal formation of increasingly helical pre-eruption structures and instigates further research on the triggering of solar eruptions with magnetic helicity firmly placed in the eruption cadre.

79 ASTRONOMY AND ASTROPHYSICS↗

A PROPELLER-EFFECT INTERPRETATION OF MAXI/GSC LIGHT CURVES OF 4U 1608-52 AND Aql X-1 AND APPLICATION TO XTE J1701-462

We present the luminosity dwell-time distributions during the hard states of two low-mass X-ray binaries containing a neutron star (NS), 4U 1608-52 and Aql X-1, observed with MAXI/GSC. The luminosity distributions show a steep cutoff on the low-luminosity side at {approx}1.0 Multiplication-Sign 10{sup 36} erg s{sup -1} in both sources. The cutoff implies a rapid luminosity decrease in their outburst decay phases and this decrease can be interpreted as being due to the propeller effect. We estimate the surface magnetic field of 4U 1608-52 to be (0.5-1.6) Multiplication-Sign 10{sup 8} G and Aql X-1 to be (0.6-1.9) Multiplication-Sign 10{sup 8} G from the cutoff luminosity and apply the same propeller mechanism to the similar rapid luminosity decrease observed in the transient Z source, XTE J1701-462, with RXTE/ASM. Assuming that the spin period of the NS is on the order of milliseconds, the observed cutoff luminosity implies a surface magnetic field on the order of 10{sup 9} G.

79 ASTRONOMY AND ASTROPHYSICS↗

INTERPRETING THE GLOBAL 21 cm SIGNAL FROM HIGH REDSHIFTS. I. MODEL-INDEPENDENT CONSTRAINTS

The sky-averaged (global) 21 cm signal is a powerful probe of the intergalactic medium (IGM) prior to the completion of reionization. However, so far it has been unclear whether it will provide more than crude estimates of when the universe's first stars and black holes formed, even in the best case scenario in which the signal is accurately extracted from the foregrounds. In contrast to previous work, which has focused on predicting the 21 cm signatures of the first luminous objects, we investigate an arbitrary realization of the signal and attempt to translate its features to the physical properties of the IGM. Within a simplified global framework, the 21 cm signal yields quantitative constraints on the Lyα background intensity, net heat deposition, ionized fraction, and their time derivatives without invoking models for the astrophysical sources themselves. The 21 cm absorption signal is most easily interpreted, setting strong limits on the heating rate density of the universe with a measurement of its redshift alone, independent of the ionization history or details of the Lyα background evolution. In a companion paper, we extend these results, focusing on the confidence with which one can infer source emissivities from IGM properties.

79 ASTRONOMY AND ASTROPHYSICS↗

An interpretation of flare-induced and decayless coronal-loop oscillations as interference patterns

We present an alternative model of coronal-loop oscillations, which considers that the waves are trapped in a two-dimensional waveguide formed by the entire arcade of field lines. This differs from the standard one-dimensional model which treats the waves as the resonant oscillations of just the visible bundle of field lines. Within the framework of our two-dimensional model, the two types of oscillations that have been observationally identified, flare-induced waves and 'decayless' oscillations, can both be attributed to MHD fast waves. The two components of the signal differ only because of the duration and spatial extent of the source that creates them. The flare-induced waves are generated by strong localized sources of short duration, while the decayless background can be excited by a continuous, stochastic source. Further, the oscillatory signal arising from a localized, short-duration source can be interpreted as a pattern of interference fringes produced by waves that have traveled diverse routes of various pathlengths through the waveguide. The resulting amplitude of the fringes slowly decays in time with an inverse square root dependence. The details of the interference pattern depend on the shape of the arcade and the spatial variation of the Alfvén speed. The rapid decay of this wave component, which has previously been attributed to physical damping mechanisms that remove energy from resonant oscillations, occurs as a natural consequence of the interference process without the need for local dissipation.

79 ASTRONOMY AND ASTROPHYSICS↗

OGLE-2012-BLG-0455/MOA-2012-BLG-206: Microlensing event with ambiguity in planetary interpretations caused by incomplete coverage of planetary signal

Characterizing a microlensing planet is done by modeling an observed lensing light curve. In this process, it is often confronted that solutions of different lensing parameters result in similar light curves, causing difficulties in uniquely interpreting the lens system, and thus understanding the causes of different types of degeneracy is important. In this work, we show that incomplete coverage of a planetary perturbation can result in degenerate solutions even for events where the planetary signal is detected with a high level of statistical significance. We demonstrate the degeneracy for an actually observed event OGLE-2012-BLG-0455/MOA-2012-BLG-206. The peak of this high-magnification event (A {sub max} ∼ 400) exhibits very strong deviation from a point-lens model with Δχ{sup 2} ≳ 4000 for data sets with a total of 6963 measurements. From detailed modeling of the light curve, we find that the deviation can be explained by four distinct solutions, i.e., two very different sets of solutions, each with a twofold degeneracy. While the twofold (so-called close/wide) degeneracy is well understood, the degeneracy between the radically different solutions is not previously known. The model light curves of this degeneracy differ substantially in the parts that were not covered by observation, indicating that the degeneracy is caused by the incomplete coverage of the perturbation. It is expected that the frequency of the degeneracy introduced in this work will be greatly reduced with the improvement of the current lensing survey and follow-up experiments and the advent of new surveys.

79 ASTRONOMY AND ASTROPHYSICS↗

RED NOISE VERSUS PLANETARY INTERPRETATIONS IN THE MICROLENSING EVENT OGLE-2013-BLG-446

For all exoplanet candidates, the reliability of a claimed detection needs to be assessed through a careful study of systematic errors in the data to minimize the false positives rate. We present a method to investigate such systematics in microlensing data sets using the microlensing event OGLE-2013-BLG-0446 as a case study. The event was observed from multiple sites around the world and its high magnification (A{sub max} ∼ 3000) allowed us to investigate the effects of terrestrial and annual parallax. Real-time modeling of the event while it was still ongoing suggested the presence of an extremely low-mass companion (∼3M{sub ⨁}) to the lensing star, leading to substantial follow-up coverage of the light curve. We test and compare different models for the light curve and conclude that the data do not favor the planetary interpretation when systematic errors are taken into account.

79 ASTRONOMY AND ASTROPHYSICS↗

Interpretation of rotation and momentum transport in the DIII-D edge plasma and comparison with neoclassical theory

Here, a low-confinement mode discharge which optimizes the capability of the new main-ion chargeexchange-recombination spectroscopy system on DIII-D to measure deuterium toroidal velocity is interpretted in comparison with the predictions of neoclassical theory, with an emphasis on the plasma edge region. The observed peaking in the deuterium toroidal velocity near the separatrix is shown to be consistent with intrinsic co-rotation due to ion orbit loss. In general, the standard neoclassical toroidal and poloidal momentum transport rates are found to be smaller than those inferred from experiment, but a comparison has not yet been made with the more recent extended neoclassical theory that calculates the effects of poloidal asymmetries using an elongated flux surface representation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

AN INTERPRETATION OF GLE71 CONCURRENT CME-DRIVEN SHOCK WAVE

Particle accelerations in solar flares and CME-driven shocks can sometimes result in very high-energy particle events (≥1 GeV) that are known as ground level enhancements (GLEs). Recent studies on the first GLE event (GLE71 2012 May 17 01:50 UT) of solar cycle 24 suggested that CME-driven shock played a leading role in causing the event. To verify this claim, we have made an effort to interpret the GLE71 concurrent shock wave. For this, we have deduced the possible speed and height of the shock wave in terms of the frequency (MHz) of the solar radio type II burst and its drift rate (MHz min{sup –1}), and studied the temporal evolution of the particle intensity profiles at different heights of the solar corona. For a better perception of the particle acceleration in the shock, we have studied the solar radio type II burst with concurrent solar radio and electron fluxes. When the particle intensity profiles are necessarily shifted in time at ∼1 AU, it is found that the growth phases of the electron and cosmic ray intensity fluxes are strongly correlated (>0.91; ≥0.87) with the frequency drift rate of the type II burst, which is also consistent with the intensive particle accelerations at upper coronal heights (∼≥0.80 R {sub S} < 1.10 R {sub S}). Thus, we conclude that the CME-driven shock was possibly capable of producing the high-energy particle event. However, since the peaks of some flare components are found to be strongly associated with the fundamental phase of the type II burst, the preceding flare is supposed to contribute to the shock acceleration process.

79 ASTRONOMY AND ASTROPHYSICS↗