Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model 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 1,603 records · Page 89

Geological terrain models

The initial phase of a program to determine the best interpretation strategy and sensor configuration for a radar remote sensing system for geologic applications is discussed. In this phase, terrain modeling and radar image simulation were used to perform parametric sensitivity studies. A relatively simple computer-generated terrain model is presented, and the data base, backscatter file, and transfer function for digital image simulation are described. Sets of images are presented that simulate the results obtained with an X-band radar from an altitude of 800 km and at three different terrain-illumination angles. The simulations include power maps, slant-range images, ground-range images, and ground-range images with statistical noise incorporated. It is concluded that digital image simulation and computer modeling provide cost-effective methods for evaluating terrain variations and sensor parameter changes, for predicting results, and for defining optimum sensor parameters.

Kaupp, V. H.↗

Quantitative Power System Resilience Metrics and Evaluation Approach

Power system resilience is an emerging topic and plays an essential role in helping the power industry understand and respond to the increasing threats of extreme weather events. The first step of power system resilience analysis is to introduce metrics to quantify the resilience reasonably. Existing resilience metrics are typically restrained by the limited data for extreme event modeling and fall short in terms of physical interpretation and comparability. This paper develops novel quantitative metrics to evaluate power system resilience in pre- and post-event contexts. The developed metrics illustrate clear physical meanings and can be effectively used to compare resilience across different systems under different extreme events. Moreover, the developed metrics can be applied to both transmission and distribution systems. Simulation on a distribution system is employed to validate the effectiveness of the proposed resilience metrics and resilience evaluation approach.

power system resilience↗

Baryon number violation: from nuclear matrix elements to BSM physics

Processes that violate baryon number, most notably proton decay and $n\bar{n}$ transitions, are promising probes of physics beyond the Standard Model (BSM) needed to understand the lack of antimatter in the Universe. To interpret current and forthcoming experimental limits, theory input from nuclear matrix elements to UV complete models enters. Thus, an interplay of experiment, effective field theory, lattice QCD, and BSM model building is required to develop strategies to accurately extract information from current and future data and maximize the impact and sensitivity of next-generation experiments. Here, we briefly summarize the main results and discussions from the workshop ‘INT-25-91W: Baryon Number Violation: From Nuclear Matrix Elements to BSM Physics,’ held at the Institute for Nuclear Theory, University of Washington, Seattle, WA, 13–17 January 2025.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Quantification of the Impact of Acidified Brine on Fracture-Matrix Transport in a Naturally Fractured Shale Using in Situ Imaging and Modeling

Understanding flow, transport, chemical reactions, and hydromechanical processes in fractured geologic materials is key for optimizing a range of subsurface processes including carbon dioxide and hydrogen storage, unconventional energy resource extraction, and geothermal energy recovery. Flow and transport processes in naturally fractured shale rocks have been challenging to characterize due to experimental complexity and the multiscale nature of quantifying continuum scale descriptions of mass exchange between micrometer-scale fractures and nanometer-scale pores. In this study, we use positron emission tomography (PET) to image the transport of a conservative tracer in a naturally fractured Wolfcamp shale core before and after the core was exposed to low pH brine conditions. Image-based experimental observations are interpreted by fitting an analytical transport model to fracture-containing voxels in the core. Results of this analysis indicate subtle increases in matrix diffusivity and a slightly more uniform fracture velocity distribution following exposure to low pH conditions. These observations are compared with a multicomponent one-dimensional reactive transport model that indicates the capacity for a 10% increase in porosity at the fracture-matrix interface as a result of the low pH brine exposure. This porosity change is the result of the dissolution of carbonate minerals in the shale matrix to low pH conditions. Further, this image-based workflow represents a new approach for quantifying spatially resolved fracture-matrix transport processes and provides a foundation for future work to better understand the role of coupled transport, reaction, and mechanical processes in naturally fractured rocks.

58 GEOSCIENCES↗

Search for heavy Higgs bosons decaying to a top quark pair in proton-proton collisions at $ \sqrt{s} $ = 13 TeV

A search is presented for additional scalar (H) or pseudoscalar (A) Higgs bosons decaying to a top quark pair in proton-proton collisions at a center-of-mass energy of 13 TeV. The data set analyzed corresponds to an integrated luminosity of 35.9 fb$^{-1}$ collected by the CMS experiment at the LHC. Final states with one or two charged leptons are considered. The invariant mass of the reconstructed top quark pair system and variables that are sensitive to the spin of the particles decaying into the top quark pair are used to search for signatures of the H or A bosons. The interference with the standard model top quark pair background is taken into account. A moderate signal-like deviation compatible with an A boson with a mass of 400 GeV is observed with a global significance of 1.9 standard deviations. New stringent constraints are reported on the strength of the coupling of the hypothetical bosons to the top quark, with the mass of the bosons ranging from 400 to 750 GeV and their total relative width from 0.5 to 25%. The results of the search are also interpreted in a minimal supersymmetric standard model scenario. Values of $m_\mathrm{A}$ from 400 to 700 GeV are probed, and a region with values of $\tan\beta$ below 1.0 to 1.5, depending on $m_\mathrm{A}$, is excluded at 95% confidence level.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Interpreting the spatial distribution of soil properties with a physically-based distributed hydrological model

Digital soil maps are commonly data-driven as the development of physically-based models for soil mapping is difficult due to the complexity of soils. However, physically-based hydrologic models have been successful in simulating water dynamics. Since water movement is a major driver of pedogenesis, the physical rules that govern water movement might help explain and predict the spatial variation of soil properties. Here, we demonstrate the novel use of a physically-based, distributed hydrologic model to inform the spatial distribution of soil properties. The Distributed Hydrology Soil Vegetation Model (DHSVM) was utilized to simulate soil moisture content (SM) and water table depth (WTD) in two hillslope catchments under pasture and forest management wherein hydrologic model outputs were then compared with soil properties measured in situ. SM sensors and wells were installed in both catchments to validate simulations of soil water movement via Nash-Sutcliffe Efficiency (E). In-situ observations were made at 87 sites within both catchments to study the connection between simulated water movement (SM and WTD) and observed soil properties, namely the depth and thickness of the argillic (Bt), fragic (Btx), and C horizons, and the depth of redoximorphic features. The simulated time series of SM and WTD were also clustered per season using Dynamic Time Warping (DTW), which identified similarity among time series at varying timescales. Model validation suggested that simulations of surficial SM (0–20 cm) were reasonable (E = 0.45), however, simulated subsurface SM (45–60 cm) and WTD were not sufficiently accurate. The thickness of Btx horizons were spatially grouped into different populations by SM clusters from every season except spring. For the other properties, only SM dynamics of specific seasons grouped into significantly different populations, suggesting that the explanatory power of simulated water movement varies seasonally and was greater during winter. Here, we show clusters of simulated SM separated soil properties into statistically different populations, showing that hydrologic models could inform areas that followed different water dynamics related to pedogenic trajectories and related biogeochemical processes not necessarily simulated by the model. As such, physically-based modeling of water dynamics can, therefore, inform and advance digital soil mapping by linking water movement patterns stemming from hydrologic model outputs to spatial patterns of soil properties and pedogenesis.

54 ENVIRONMENTAL SCIENCES↗

A Unification Framework for Euclidean and Hyperbolic Graph Neural Networks

Hyperbolic neural networks have recently gained significant attention due to their promising results on several graph problems including node classification and link prediction. The primary reason for this success is the effectiveness of hyperbolic space in capturing the inherent hierarchy of graph datasets. However, they are limited in terms of generalization, scalability, and have inferior performance when applied to non-hierarchical datasets. In this paper, we take a completely different perspective for modeling hyperbolic networks and answer the following question: is an Euclidean model able to approximate a function or behavior in the hyperbolic space? Extending the universal approximation theory developed for Euclidean models, We draw an analogy from the hyperbolic components to the Euclidean counterparts and conclude that, in order to capture hierarchical features, it is possible to generalize hyperbolic models to be a special case of Euclidean models with the proposed Pseudo-Poincaré technique. We applied our non-linear hyperbolic normalization to the current state-of-the-art homogeneous and multi-relational graph networks and demonstrate significant improvements in performance compared to both Euclidean and hyperbolic counterparts. The primary impact of this work lies in its ability to capture hierarchical features in the Euclidean space, and thus, can replace hyperbolic networks without any loss in performance metrics while simultaneously leveraging the power of Euclidean networks such as interpretability and efficient execution of various model components.

Khatir, Mehrdad↗

An interpretation of epsilon aurigae.

Rotating opaque gaseous disk model to explain contradictory behavior of binary Epsilon Aurigae found in photometric and spectroscopic observations

ROTATING DISK↗

Magnetospheric substorms - A model.

An attempt is made to integrate and interpret various polar and magnetospheric substorm phenomena as consequences of the acceleration process of auroral particles which could result from the diversion of a part of the magnetotail current to the night side of the auroral oval. Thus, this paper is not intended to be a literature survey, but to present a way along which magnetospheric substorms might be studied in the future.

Akasofu, S.-I.↗

Magnetic models of solar flares.

A solar flare is interpreted as the explosive release of magnetic energy associated with a current sheet in the solar atmosphere. The release is due to field-line reconnection caused by the tearing-mode instability. Current sheets may be classified into closed and open, depending on whether the adjacent field lines are closed (contained in the sun's atmosphere) or open, extending into interplanetary space. A high-energy event requires an open current sheet. It is shown that certain types of photospheric motion may lead to the formation of current sheets. One of these is associated with high-energy events, another with filaments, and a third with surge flares.

Sturrock, P. A.↗

Beaming of Uranian kilometric radiation and implications for the planetary magnetic field

Radio waves from Uranus at 481.2 kHz detected by Voyager 2 show a strong beaming effect. When the data are analyzed according to the magnetic latitude of the observer in an offset tilted dipole field model, a sharply defined emission beam extending from 25 to 55 deg latitude is evident consistent with a source over the south magnetic pole. The lack of prominent emission from the north magnetic pole can be attributed to the large spatial offset of the magnetic dipole resulting in much lower source frequencies. For proper alignment of the beam, an advance of the dipole axis longitude by 23 deg from the standard model is required. This phase shift is interpreted as the result of higher order magnetic moments producing a skew of the field lines at low altitudes in the polar source region. This conclusion further substantiates models of the interior of Uranus which predict a complex intrinsic field with significant nondipolar contributions to the surface magnetic field.

Barbosa, D. D.↗

Model of hot-film sensor with substrate

A detailed model is constructed to investigate parameters influencing a hot-film anemometer system used to study cross-flow phenomena on swept wings. This anemometer is designed to detect vortices set up by cross-flow. The nature of the vortices is not well known, and a well-verified instrumentation model is considered essential to allow accurate interpretation of experimental results. Previous investigators have indicated that heat conduction to the film's substrate is significant. Two variations on the inclusion of substrate effects in the model are discussed. In the first, the substrate is assumed to have a constant temperature throughput: in the second more elaborate model, the temperature of the substrate is allowed to vary. This variation is related to the convective heat removal at the substrate/gas interface.

Mcree, Griffith J.↗

The role of particle diffusion in the lower transition region: Revised interpretation of emission measures

Our energy-balance models of the lower transit ion region were presented in a previous paper (New Models of the Chromosphere and Transition Region). Here we show the influence of particle diffusion on the calculated hydrogen and helium number densities for a given temperature-density model (model C in the preceding paper). We have also solved the statistical equilibrium and radiative transfer equations for a 13-level He I atom (22 radiative transitions) and a 6-level He II ion (15 radiative transitions) together with He III. The resulting He I and He II level-1 number densities and He III density are shown as functions of temperature. Diffusion substantially increases eta(sub HeI) for T greater than 35,000 K and decreases eta(sub HeI) (while increasing eta(sub(HeII)) between 9,000 and 25,000 K. Including the effects of diffusion also increases eta(sub HeIII) for T less than 60,000 K.

Avrett, Eugene H.↗

Stormtime ring current and radiation belt ion transport: Simulations and interpretations

We use a dynamical guiding-center model to investigate the stormtime transport of ring current and radiation-belt ions. We trace the motion of representative ions' guiding centers in response to model substorm-associated impulses in the convection electric field for a range of ion energies. Our simple magnetospheric model allows us to compare our numerical results quantitatively with analytical descriptions of particle transport, (e.g., with the quasilinear theory of radial diffusion). We find that 10-145-keV ions gain access to L approximately 3, where they can form the stormtime ring current, mainly from outside the (trapping) region in which particles execute closed drift paths. Conversely, the transport of higher-energy ions (approximately greater than 145 keV at L approximately 3) turns out to resemble radial diffusion. The quasilinear diffusion coefficient calculated for our model storm does not vary smoothly with particle energy, since our impulses occur at specific (although randomly determined) times. Despite the spectral irregularity, quasilinear theory provides a surprisingly accurate description of the transport process for approximately greater than 145-keV ions, even for the case of an individual storm. For 4 different realizations of our model storm, the geometric mean discrepancies between diffusion coefficients D(sup sim, sub LL) obtained from the simulations and the quasilinear diffusion coefficient D(sup ql, sub LL) amount to factors of 2.3, 2.3, 1.5, and 3.0, respectively. We have found that these discrepancies between D(sup sim, sub LL) and D(sup ql, sub LL) can be reduced slightly by invoking drift-resonance broadening to smooth out the sharp minima and maxima in D(sup ql, sub LL). The mean of the remaining discrepancies between D(sup sim, sub LL) and D(sup ql, sub LL) for the 4 different storms then amount to factors of 1.9, 2.1, 1.5, and 2.7, respectively. We find even better agreement when we reduce the impulse amplitudes systematically in a given model storm (e.g., reduction of all the impulse amplitudes by half reduces the discrepancy factor by at least its square root) and also when we average our results over an ensemble of 20 model storms (agreement is within a factor of 1.2 without impulse-amplitude reduction). We use our simulation results also to map phase-space densities f in accordance with Liouville's theorem. We find that the stormtime transport of approximately greater than 145-keV ions produces little change in f-bar the drift-averaged phase-space density on any drift shell of interest. However, the stormtime transport produces a major enhancement from the pre-storm phase-space density at energies approximately 30-145 keV, which are representative of the stormtime ring current.

Lyons, Larry R.↗

Revealing the hidden structure of disordered materials by parameterizing their local structural manifold

Abstract Durable interest in developing a framework for the detailed structure of glassy materials has produced numerous structural descriptors that trade off between general applicability and interpretability. However, none approach the combination of simplicity and wide-ranging predictive power of the lattice-grain-defect framework for crystalline materials. Working from the hypothesis that the local atomic environments of a glassy material are constrained by enthalpy minimization to a low-dimensional manifold in atomic coordinate space, we develop a generalized distance function, the Gaussian Integral Inner Product (GIIP) distance, in connection with agglomerative clustering and diffusion maps, to parameterize that manifold. Applying this approach to a two-dimensional model crystal and a three-dimensional binary model metallic glass results in parameters interpretable as coordination number, composition, volumetric strain, and local symmetry. In particular, we show that a more slowly quenched glass has a higher degree of local tetrahedral symmetry at the expense of cyclic symmetry. While these descriptors require post-hoc interpretation, they minimize bias rooted in crystalline materials science and illuminate a range of structural trends that might otherwise be missed.

36 MATERIALS SCIENCE↗

Advances in 3D transient plasma dynamics and control through MHD and hybrid fluid-kinetic simulations with JOREK

Transient phenomena and their control are of high relevance in magnetic confinement fusion plasmas to guarantee a stable and safe plasma operation. Interpretative simulations can maximize the insights gained from experiments on present machines and predictive simulations can help in the preparation of design, mitigation techniques and operational scenarios for future devices. In this article, we provide an overview of recent advances and novel scientific results obtained with the 3D non-linear hybrid fluid-kinetic code JOREK, covering physics of plasma transients from the core to the scrape-off layer (SOL) both for tokamak and stellarator devices. Substantial progress was made in the physics understanding, model validation with experiments and experiment interpretation, thus, giving confidence for predictions to devices like DTT, ITER and DEMO. The topics addressed comprise a wide range: the edge physics of new operation scenarios and edge localized mode suppression; major disruptions with a focus on runaway electrons and vertical displacement events as well as disruption mitigation by shattered pellet injection; the physics mechanisms and operational limits of the flux pumping regime for sawtooth control; MHD limits of stellarators and work towards incorporating advanced edge/SOL/exhaust dynamics; continuing improvements of the code for more efficient hybrid simulations on conventional and accelerated high performance computing architectures.

disruptions↗

Unraveling Hydrogen Induced Geochemical Reaction Mechanisms through Coupled Geochemical Modeling and Machine Learning

Underground hydrogen storage (UHS) provides a promising large-scale, long-term energy storage solution. A reasonable recovery of stored hydrogen is critical for a successful storage scheme. However, in subsurface reservoirs hydrogen is subject to active geochemical reactions that might result in hydrogen loss. In this study, we implemented a geochemical modeling approach coupled with an unsupervised machine learning technique called non-negative matrix factorization (NMF) to unravel the complex brine-rock-H 2 geochemical processes responsible for hydrogen losses, with particular focus on sulfate reduction reactions. NMF is applied to modeled mineral evolution and fluid component profiles to retrieve profiles that can be interpreted to more easily assess competing processes. NMF decouples simulated competing equilibrium reactions. This facilitates separation of overlapping reaction profiles from redox processes, dissolution fronts, and secondary precipitation while considering the effects of simulation parameters such as salinity, temperature, and total H 2 pressure. NMF successfully discriminates these competing effects in nonlinear ways, allowing robust interpretation. In addition, NMF reveals subtle coupled mineral associations and reaction fronts that are invisible to conventional model analysis. This integrated approach strengthens the conceptual understanding of complex nonlinear hydrogen-brine-rock interactions and advances geochemical research on UHS systems to resolve complexities in modeled geochemical systems without the need for direct experiments or prior knowledge. Furthermore, this study highlights the efficacy of combining geochemical modeling with machine learning techniques to enhance the interpretability of the intricate geochemical simulation output through deciphering the overlapping reaction path that cannot be achieved only using conventional analysis of geochemical models alone.

08 HYDROGEN↗