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 379 records · Page 21

Towards physics-informed explainable machine learning and causal models for materials research

From emergent material descriptions to estimation of properties stemming from structures to optimization of process parameters for achieving best performance – all key facets of materials science and related fields have experienced tremendous growth with the introduction of data-driven models. This gradual progression goes at par with developments of machine learning workflows, from purely data-driven shallow models to those that are well-capable in encoding more complex graphs, symbolic representations, invariances, and positional embeddings. Furthermore, this perspective aims at summarizing strategic aspects of such transitions while providing insights into the requirements of bringing in explainable, interpretable predictive models, and causal learning to aid in materials design and discovery. Although the focus remains on a variety of functional materials by providing a handful of case studies, the applications of such integrated methodologies are universal to facilitate fundamental understandings of materials physics while enabling autonomous experiments.

36 MATERIALS SCIENCE↗

Parameter Estimation for Electrode Degradation: Learning in the Face of Model-Experiment Discrepancies

Use of physics-based models to interpret battery degradation data over the course of cycling can provide deeper physical insight into the internal states of the system and how they evolve. We present a neural network trained on simulations generated by a previously published physics-based model for a lithium trivanadate (LVO) cathode to estimate parameters that evolve over the course of cycling. We focus on the robustness of the neural network through two case studies that probe different kinds of discrepancies between model and experiment: nonideal data and imperfect model. In the former, the experimental protocols do not meet the assumption made in the training data generated by the physics-based model, while in the latter, the physics-based model fails to describe all of the measured cathode behavior even under ideal conditions. When there is total model-experiment agreement, a neural network estimates parameters with improved accuracy compared to a maximum likelihood analysis using the same set of simulations. However, in both types of model-experiment discrepancy, the neural network returned biased parameter estimates. We introduce a data augmentation procedure into the neural network training to mitigate these effects and improve robustness, and employ it to estimate parameters for a cycling LVO cathode.

Mayilvahanan, Karthik S. (ORCID:0000000316282332)↗

High-throughput phase-field simulations and machine learning of resistive switching in resistive random-access memory

Metal oxide-based Resistive Random-Access Memory (RRAM) exhibits multiple resistance states, arising from the activation/deactivation of a conductive filament (CF) inside a switching layer. Understanding CF formation kinetics is critical to achieving optimal functionality of RRAM. Here a phase-field model is developed, based on materials properties determined by ab initio calculations, to investigate the role of electrical bias, heat transport and defect-induced Vegard strain in the resistive switching behavior, using MO 2- x systems such as HfO 2- x as a prototypical model system. It successfully captures the CF formation and resultant bipolar resistive switching characteristics. High-throughput simulations are performed for RRAMs with different material parameters to establish a dataset, based on which a compressed-sensing machine learning is conducted to derive interpretable analytical models for device performance (current on/off ratio and switching time) metrics in terms of key material parameters (electrical and thermal conductivities, Vegard strain coefficients). These analytical models reveal that optimal performance (i.e., high current on/off ratio and low switching time) can be achieved in materials with a low Lorenz number, a fundamental material constant. This work provides a fundamental understanding to the resistive switching in RRAM and demonstrates a computational data-driven methodology of materials selection for improved RRAM performance, which can also be applied to other electro-thermo-mechanical systems.

36 MATERIALS SCIENCE↗

On the Importance of Laboratory Astrophysics and Astrochemistry and Interdisciplinary Research: Two Success Stories

Here we present two examples that demonstrate how cross-disciplinary research projects, where experimentalists, modelers and observers work together to answer science questions, allow expertise to be shared and misconceptions or missing key elements to be tackled by looking at the problem from different perspectives. They also show how, by working in unison, the group can accomplish more than the sum of its parts by combining results into a higher-level understanding of the chemical processes taking place. The first project is a collaborative study between experimentalists, modelers and observers to1) produce laboratory analogs of cosmic grains and planetary aerosols (Titan, Pluto...) from different gas mixtures in cold astrophysically relevant conditions (<200 K); and 2) characterize them with scanning electron microscopy and visible-to-far-infrared spectroscopy in order to assess the impact of the precursors on their growth structure and optical properties. We can then produce and study analogs that are representative of different formation stages or environments, and provide their complex refractive indices to the scientific community. We will show how, by using these experimental optical constants of more representative analogs in radiative transfer and reflectance spectra models, better interpretations of (exo)planetary atmosphere- and surface observations are possible. The second project is an interdisciplinary study of the formation of benzene clouds in the atmosphere of Saturn’s largest moon, Titan. We will show how combining Earth and Planetary Science laboratory expertise, modeling and observations has led to providing to the scientific community, for the first time, experimental vapor pressures for benzene at cold temperatures(<200K) relevant to Titan’s atmospheric conditions. These have been used in microphysical models instead of the higher temperature extrapolations used previously, allowing a better match to observations.

Ella Sciamma-o'Brien↗

Heat Flow and Magnetization in the Oceanic Lithosphere

Two aspects of the processing and interpretation of satellite measurements of the geomagnetic field are described. One deals with the extraction of the part of the geomagnetic field that originates from sources in the earth's atmosphere. The other investigates the possibility of using the thermal state of the oceanic lithosphere to further constrain modelling and interpretation of magnetic anomalies. It is shown that some of the magnetic signal in crustal anomaly maps can be an artifact of the mathematical algorithms that have been used to separate the crustal field from the observed data. Strong magnetic anomalies can be distorted but are probably real, but weak magnetic anomalies can arise from leakage of power from short wavelengths, and will also appear in anomaly maps as repetitions of the strong crustal anomaly. The distortion and the ghost anomalies follow the magnetic dip lines in a way that is similar to actual MAGSAT anomaly fields. This phenomenon will also affect the lower degree spherical harmonic terms in the power spectrum of the crustal field. A model of the magnetic properties of the oceanic crust that has been derived from direct measurements of the rock magnetic properties of oceanic rocks is presented. The average intensity of magnetization in the oceanic crust is not strong enough to explain magnetic anomalies observed over oceanic areas. This is the case for both near surface observations (ship and aeromagnetic data) and satellite altitude observations. It is shown that magnetic sources in the part of the upper mantle that is situated above the Curie isotherm, if sufficiently strong, can produce satellite magnetic anomalies that are comparable to MAGSAT data. The method developed for the study of depth to the Curie isotherm and magnetic anomalies can also be used in inverse modelling of satellite magnetic anomalies when the model is to be adjusted with an annihilator.

Hayling, Kjell Lennart↗

Measurements and Model Improvement: Insight into NWP Model Error Using Doppler Lidar and Other WFIP2 Measurement Systems

Abstract Doppler-lidar wind-profile measurements at three sites were used to evaluate NWP model errors from two versions of NOAA’s 3-km-grid HRRR model, to see whether updates in the latest version 4 reduced errors when compared against the original version 1. Nested (750-m grid) versions of each were also tested to see how grid spacing affected forecast skill. The measurements were part of the field phase of the Second Wind Forecasting Improvement Project (WFIP2), an 18-month deployment into central Oregon–Washington, a major wind-energy-producing region. This study focuses on errors in simulating marine intrusions, a summertime, 600–800-m-deep, regional sea-breeze flow found to generate large errors. HRRR errors proved to be complex and site dependent. The most prominent error resulted from a premature drop in modeled marine-intrusion wind speeds after local midnight, when lidar-measured winds of greater than 8 m s −1 persisted through the next morning. These large negative errors were offset at low levels by positive errors due to excessive mixing, complicating the interpretation of model “improvement,” such that the updates to the full-scale versions produced mixed results, sometimes enhancing but sometimes degrading model skill. Nesting consistently improved model performance, with version 1’s nest producing the smallest errors overall. HRRR’s ability to represent the stages of sea-breeze forcing was evaluated using radiation budget, surface-energy balance, and near-surface temperature measurements available during WFIP2. The significant site-to-site differences in model error and the complex nature of these errors mean that field-measurement campaigns having dense arrays of profiling sensors are necessary to properly diagnose and characterize model errors, as part of a systematic approach to NWP model improvement. Significance Statement Dramatic increases in NWP model skill will be required over the coming decades. This paper describes the role of major deployments of accurate profiling sensors in achieving that goal and presents an example from the Second Wind Forecast Improvement Program (WFIP2). Wind-profile data from scanning Doppler lidars were used to evaluate two versions of HRRR, the original and an updated version, and nested versions of each. This study focuses on the ability of updated HRRR versions to improve upon predicting a regional sea-breeze flow, which was found to generate large errors by the original HRRR. Updates to the full-scale HRRR versions produced mixed results, but the finer-mesh versions consistently reduced model errors.

Meteorology & Atmospheric Sciences↗

Investigation of passive atmospheric sounding using millimeter and submillimeter wavelength channels

Progress by the Georgia Institute of Technology's Laboratory for Radio-science and Remote Sensing in developing techniques for passive microwave retrieval of water vapor profiles and cloud and precipitation parameters using millimeter and submillimeter wavelength channels is reviewed. Channels of particular interest are in the tropospheric transmission windows at 90, 166, 220, 340, and 410 GHz and centered around the water vapor lines at 183 and 325 GHz. Collectively, these channels have potential application in high-resolution precipitation mapping (e.g., from geosynchronous orbit), remote sensing of cloud and precipitation parameters, including cirrus ice mass, and improved retrieval of water vapor profiles. During the period from January 1, 1994 through June 30, 1994 research activities focussed on calibrating and interpreting data from the Millimeter-Wave Imaging Radiometer (MIR). The MIR was deployed on the NASA ER-2 during the Convective Atmospheric Moisture Experiment (CAMEX, September-October 1993) to obtain the first submillimeter-wave tropospheric imagery of convective precipitations. A 325-GHz radiometer consisted of a submillimeter-wave DSB receiver with three IF channels at +/- 1, 3, and 8.5 GHz, and approximately 14 dB DSB noise figure was successfully operated during these experiments. Activities supported under this grant include a study of the impact of local oscillator reflections from the MIR calibration loads, the development of optimal gain and offset filters for radiometric calibration, and the modeling and interpretation of the MIR 325-GHz data over both clear and cloudy atmospheres. In addition, polarimetric radiometer measurements and modeling for ocean surface and atmospheric cloud-ice studies_were supported.

Gasiewski, Albin J.↗

A new method for predicting hurricane rapid intensification based on co-occurring environmental parameters

Abstract Tropical cyclones (TCs) that undergo Rapid Intensification (RI) can pose serious socioeconomic threats and can potentially result in major damaging impacts along coastal areas. Considering the complexity of various physical mechanisms that play a role in RI and its relatively low probability of occurrence, predicting RI remains a major operational challenge. In this study, we propose a simple deterministic binary classification model based on the co-occurrence of environmental parameters (MCE) to predict an RI event. More specifically, the model determines the possibility of RI based on a simple count of the number of environmental predictors deemed favorable and unfavorable. We compare our model results to logistic regression (LR) and decision tree (DT) models, well-trained using the same set of environmental predictors. Results reveal that at an RI threshold of 30 kt, the MCE exhibits a critical success index score of 0.233 which is 14% higher than DT and LR model performances. When tested at multiple RI thresholds, the MCE displays relatively higher skill scores across multiple metrics. By simultaneously evaluating the favorability of predictors, the MCE is able to comparatively reduce the number of false alarms predicted when certain predictors are unfavorable toward RI. Interpreting these model results to gain a physical understanding of how co-occurring environmental parameters can affect RI, we highlight future directions for using models based on the MCE approach to understand and predict TC RI as well as other meteorological extremes.

54 ENVIRONMENTAL SCIENCES↗

A Crustal Structure Study of South America

Variations in the average crustal thickness of the South American platform were determined by measuring fundamental mode Rayleigh wave dispersion between nine pairs of South American WWSSN seismograph stations and inverting the measurements to obtain average shear wave velocity-depth models across the platform. Additional models were derived from the Rayleigh wave dispersion data previously measured over eastern South America. Both sets of models were interpreted in terms of average crustal thickness and average crustal shear wave velocity. The results obtained were depicted in a contour map of crustal thickness. Average values for seismic velocities in the crust and upper mantle were calculated and compared with maps of regional gravity models, MAGSAT magnetic anomalies, and major tectonic features. The correlations and relationships discovered are discussed.

Renbarger, K. S.↗

Electron inelastic mean free path in water

Liquid phase transmission electron microscopy (LPTEM) is rapidly developing as a powerful tool for probing processes in liquid environments with close to atomic resolution. Knowledge of the water thickness is needed for reliable interpretation and modelling of analytical studies in LPTEM, and is particularly essential when using thin liquid layers, required for achieving the highest spatial resolutions. Furthermore, the log-ratio method in electron energy-loss spectroscopy (EELS) is often applied in TEM to quantify the sample thickness, which is measured relative to the inelastic mean free path ( λ IMFP ). However, λ IMFP itself is dependent on sample material, the electron energy, and the convergence and divergence angles of the microscope electronoptics. Here, we present a detailed quantitative analysis of the λ IMFP of water as functions of the EELS collection angle ( β ) at 120 keV and 300 keV in a novel nanochannel liquid cell. We observe good agreement with earlier studies conducted on ice, but find that the most widely used theoretical models significantly underestimate λ IMFP of water. We determine an adjusted average energy-loss term E m, water , and characteristic scattering angle θ E, water that improve the accuracy. The results provide a comprehensive knowledge of the λ IMFP of water (or ice) for reliable interpretation and quantification of observations in LPTEM and cryo-TEM studies.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

SOHIP Abel Transform and Onion Peeling Model Module

This software provides tools for analyzing and modeling physical systems using mathematical transforms and layered models. It includes (1) functions for performing the Abel transform, which is used to relate measurements of bending angles to properties such as refractive index and radius in a medium. The code can compute bending angles from input profiles and also reconstruct these profiles from observed data; (2) the functions for modeling systems with multiple layers using an onion-peeling approach, allowing users to simulate and analyze the behavior of layered materials or structures. These capabilities are useful for researchers and engineers working in fields such as optics, atmospheric science, and materials analysis, enabling them to interpret and model data from experiments or simulations relates to refraction in spherical symmetric medium.

Xu, Shuang [Lawrence Livermore National Laboratory↗

Stable carbon isotope values of syndepositional carbonate spherules and micrite record spatial and temporal changes in photosynthesis intensity

Marine and lacustrine carbonate minerals preserve carbon cycle information, and their stable carbon isotope values (δ 13 C) are frequently used to infer and reconstruct paleoenvironmental changes. However, multiple processes can influence the δ 13 C values of bulk carbonates, confounding the interpretation of these values in terms of conditions at the time of mineral precipitation. Co-existing carbonate forms may represent different environmental conditions, yet few studies have analyzed δ 13 C values of syndepositional carbonate grains of varying morphologies to investigate their origins. Here, we combine stable isotope analyses, metagenomics, and geochemical modeling to interpret δ 13 C values of syndepositional carbonate spherules (>500 μm) and fine-grained micrite (<63 μm) from a ~1600-year-long sediment record of a hypersaline lake located on the coral atoll of Kiritimati, Republic of Kiribati (1.9°N, 157.4°W). Petrographic, mineralogic, and stable isotope results suggest that both carbonate fractions precipitate in situ with minor diagenetic alterations. The δ 13 C values of spherules are high compared to the syndepositional micrite and cannot be explained by mineral differences or external perturbations, suggesting a role for local biological processes. We use geochemical modeling to test the hypothesis that the spherules form in the surface microbial mat during peak diurnal photosynthesis when the δ 13 C value of dissolved inorganic carbon is elevated. In contrast, we hypothesize that the micrite may precipitate more continuously in the water as well as in sub-surface, heterotrophic layers of the microbial mat. Both metagenome and geochemical model results support a critical role for photosynthesis in influencing carbonate δ 13 C values. The down-core spherule–micrite offset in δ 13 C values also aligns with total organic carbon values, suggesting that the difference in the δ13C values of spherules and micrite may be a more robust, inorganic indicator of variability in productivity and local biological processes through time than the δ 13 C values of individual carbonate forms.

59 BASIC BIOLOGICAL SCIENCES↗

​Understanding a Stratigraphic Hydrothermal Resource – Geophysical Imaging at Steptoe Valley, Nevada

Sandia National Laboratories partnered with a multi-disciplinary group of subject matter experts to evaluate a stratigraphic geothermal resource in Steptoe Valley, Nevada using both established and novel geophysical imaging techniques. The stratigraphic reservoir in northern Steptoe Valley was previously discovered during oil and gas exploration. Subsequent studies, such as the Nevada Play Fairway Analysis, included data which further highlighted potential resource targets in the basin. Geophysical surveys, complemented with refined geologic mapping and geochemical sampling, were deployed to further characterize the resource. The resulting 3D geologic interpretation, conceptual model refinements, and reservoir simulations suggest that a >100MWe power-capable reservoir is likely economically accessible using conventional well placement and stimulation techniques in the Paleozoic carbonates of the deep/central basin of northern Steptoe Valley. Additional geophysical characterization and exploration drilling efforts are recommended to calibrate interpretation and determine where/how to potentially develop the northern Steptoe resource. The geophysical tools, interpretations, lessons learned, and public data generated by this study establish an exploration methodology to inform decisions for characterization and development of northern Steptoe Valley and other stratigraphic geothermal reservoirs in the western U.S.

15 GEOTHERMAL ENERGY↗

Carbon isotope effects associated with aceticlastic methanogenesis

The carbon isotope effects associated with synthesis of methane from acetate have been determined for Methanosarcina barkeri 227 and for methanogenic archaea in sediments of Wintergreen Lake, Michigan. At 37 degrees C, the 13C isotope effect for the reaction acetate (methyl carbon) --> methane, as measured in replicate experiments with M. barkeri, was - 21.3% +/- 0.3%. The isotope effect at the carboxyl portion of acetate was essentially equal, indicating participation of both positions in the rate-determining step, as expected for reactions catalyzed by carbon monoxide dehydrogenase. A similar isotope effect, - 19.2% +/- 0.3% was found for this reaction in the natural community (temperature = 20 degrees C). Given these observations, it has been possible to model the flow of carbon to methane within lake sediment communities and to account for carbon isotope compositions of evolving methane. Extension of the model allows interpretation of seasonal fluctuations in 13C contents of methane in other systems.

NASA Discipline Exobiology↗

K/TH in Achondrites and Interpretation of Grand Data for the Dawn Mission

The Dawn mission will explore 4 Vesta [1], a highly differentiated asteroid believed to be the parent body of the howardite, eucrite and diogenite (HED) meteorite suite [e.g. 2]. The Dawn spacecraft is equipped with a gamma-ray and neutron detector (GRaND), which will enable measurement and mapping of elemental abundances on Vesta s surface [3]. Drawing on HED geochemistry, Usui and McSween [4] proposed a linear mixing model for interpretation of GRaND data. However, the HED suite is not the only achondrite suite representing asteroidal basaltic crusts; others include the mesosiderites, angrites, NWA 011, and possibly Ibitira, each of which is thought to have a distinct parental asteroid [5]. Here we critically examine the variability of GRaND-analyzed elements, K and Th, in HED meteorites, and propose a method based on the K-Th systematics to distinguish between HED and the other differentiated achondrites. Maps of these elements might also recognize incompatible element enriched areas such as mapped locally on the Moon (KREEP) [6], and variations in K/Th ratios might indicate impact volatilization of K. We also propose a new mixing model using elements that will be most reliably measured by GRaND, including K.

Usui, T.↗

AK112: Full Waveform Inversion Tomography of Alaska Improves Waveform Fits While Imaging Crustal, Mantle, and Slab Structure

We report a full waveform inversion tomography model of Alaska and the surrounding regions, inferring radially anisotropic shear and isotropic compressional wavespeeds by fitting complete waveforms from 120 regional earthquakes. Our multiscale approach inverted time–frequency phase misfits (maximum period of 100 s), starting with a minimum period of 40 s and ending at 20 s in 7 stages and 112 total iterations. The model (AK112) was evaluated by computing the misfits for 36 independent validation events. We find that misfit reductions were large and equal (∼55%) for both the inversion and validation data sets, providing confidence in the model. AK112 also provides much better waveform fits compared to other reported models for the region, including an isotropic version of itself, highlighting the importance of anisotropy. The model resolves known crustal, upper mantle, and slab structure to depths of 100 km with new detail: sedimentary basins in the Alaskan Shelf, Cook Inlet, and Colville basins, among others; discontinuous lithospheric structure across major terrane boundaries; and subducting slab geometry and back‐arc volcanic sources. In addition to tectonic interpretations, the model enables full waveform simulations for long‐period earthquake ground motions and source characterization (e.g., moment tensor and finite‐fault inversion).

Rodgers, Arthur [Lawrence Livermore National Labor↗

Building kinetic models for metabolic engineering

Kinetic formalisms of metabolism link metabolic fluxes to enzyme levels, metabolite concentrations and their allosteric regulatory interactions. Though they require the identification of physiologically relevant values for numerous parameters, kinetic formalisms uniquely establish a mechanistic link across heterogeneous omics datasets and provide an overarching vantage point to effectively inform metabolic engineering strategies. Advances in computational power, gene annotation coverage, and formalism standardization have led to significant progress over the past few years. However, careful interpretation of model predictions, limited metabolic flux datasets, and assessment of parameter sensitivity remain as challenges. In this study we highlight fundamental considerations which influence model quality and prediction, advances in methodologies, and success stories of deploying kinetic models to guide metabolic engineering.

59 BASIC BIOLOGICAL SCIENCES↗

A Deep Gravity model for mobility flows generation

The movements of individuals within and among cities influence critical aspects of our society, such as well-being, the spreading of epidemics, and the quality of the environment. When information about mobility flows is not available for a particular region of interest, we must rely on mathematical models to generate them. In this work, we propose Deep Gravity, an effective model to generate flow probabilities that exploits many features (e.g., land use, road network, transport, food, health facilities) extracted from voluntary geographic data, and uses deep neural networks to discover non-linear relationships between those features and mobility flows. Our experiments, conducted on mobility flows in England, Italy, and New York State, show that Deep Gravity achieves a significant increase in performance, especially in densely populated regions of interest, with respect to the classic gravity model and models that do not use deep neural networks or geographic data. Deep Gravity has good generalization capability, generating realistic flows also for geographic areas for which there is no data availability for training. Finally, we show how flows generated by Deep Gravity may be explained in terms of the geographic features and highlight crucial differences among the three considered countries interpreting the model’s prediction with explainable AI techniques.

97 MATHEMATICS AND COMPUTING↗