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 919 records · Page 51

Complementing a continuum thermodynamic approach to constitutive modeling with symbolic regression

Here, an interpretable machine learning method, physics-informed genetic programming-based symbolic regression (P-GPSR), is integrated into a continuum thermodynamic approach to developing constitutive models. The proposed strategy for combining a thermodynamic analysis with P-GPSR is demonstrated by generating a yield function for an idealized material with voids, i.e., the Gurson yield function. First, a thermodynamic-based analysis is used to derive model requirements that are exploited in a custom P-GPSR implementation as fitness criteria or are strongly enforced in the solution. The P-GPSR implementation improved accuracy, generalizability, and training time compared to the same GPSR code without physics-informed fitness criteria. The yield function generated through the P-GPSR framework is in the form of a composite function that describes a class of materials and is characteristically more interpretable than GPSR-derived equations. The physical significance of the input functions learned by P-GPSR within the composite function is acquired from the thermodynamic analysis. Fundamental explanations of why the implemented P-GPSR capabilities improve results over a conventional GPSR algorithm are provided.

42 ENGINEERING↗

Dynamic pressure model derived from an observation by Sakigake for Comet Halley on 31 December, 1985

An outstanding disconnection event (DE)-like knot was observed on 31 Dec. 1985 in P/Halley's tail. Analysis of the Sakigake/IMF data reveals that comet Halley did not encounter the heliospheric neutral sheet on the day, demanding a new explanation of the DE-like event. During this event, the comet encountered a high-speed solar wind from a coronal hole tongue of the Sun. The event can be explained by a dynamic pressure model, according to which the DE-like plasmoid was caused by a sudden increase in the dynamic pressure of the solar wind. A simulation result is found to support this interpretation. The dynamic pressure model for a comet can be compared with the mechanism of a possible geotail disturbance during a spacecraft triggered auroral substorm.

Saito, Takao↗

Application of theoretical models to active and passive remote sensing of saline ice

The random medium model is used to interpret the polarimetric active and passive measurements of saline ice. The ice layer is described as a host ice medium embedded with randomly distributed inhomogeneities, and the underlying sea water is considered as a homogeneous half-space. The scatterers in the ice layer are modeled with an ellipsoidal correlation function. The orientation of the scatterers is vertically aligned and azimuthally random. The strong permittivity fluctuation theory is employed to calculate the effective permittivity and the distorted Born approximation is used to obtain the polarimetric scattering coefficients. We also calculate the thermal emissions based on the reciprocity and energy conservation principles. The effects of the random roughness at the air-ice, and ice-water interfaces are accounted for by adding the surface scattering to the volume scattering return incoherently. The above theoretical model, which has been successfully applied to analyze the radar backscatter data of the first-year sea ice near Point Barrow, AK, is used to interpret the measurements performed in the CRRELEX program.

Han, H. C.↗

A Universal Machine Learning Model for Elemental Grain Boundary Energies

The grain boundary (GB) energy has a profound influence on the grain growth and properties of polycrystalline metals. Here, we show that the energy of a GB, normalized by the bulk cohesive energy, can be described purely by four geometric features. By machine learning on a large computed database of 361 small Σ (Σ<10) GBs of more than 50 metals, we develop a model that can predict the grain boundary energies to within a mean absolute error of 0.13 J m –2 . More importantly, this universal GB energy model can be extrapolated to the energies of high Σ GBs without loss in accuracy. These results highlight the importance of capturing fundamental scaling physics and domain knowledge in the design of interpretable, extrapolatable machine learning models for materials science.

36 MATERIALS SCIENCE↗

Measurable characteristics of lysozyme crystal growth

The behavior of protein crystal growth is estimated from measurements performed at both the microscopic and molecular levels. In the absence of solutal flow, it was determined that a model that balances the macromolecular flux toward the crystal surface with the flux of the crystal surface well characterizes crystal growth observed using microscopic methods. Namely, it was determined that the model provides accurate estimates for the crystal-growth velocities upon evaluation of crystal-growth measurements obtained in time. Growth velocities thus determined as a function of solution supersaturation were further interpreted using established deterministic models. From analyses of crystal-growth velocities, it was found that the mode of crystal growth varies with respect to increasing solution supersaturation, possibly owing to kinetic roughening. To verify further the hypothesis of kinetic roughening, crystal growth at the molecular level was examined using atomic force microscopy (AFM). From the AFM measurements, it was found that the magnitude of surface-height fluctuations, h(x), increases with increasing solution supersaturation. In contrast, the estimated characteristic length, xi, decreases rapidly upon increasing solution supersaturation. It was conjectured that the magnitude of both h(x) and xi could possibly determine the mode of crystal growth. Although the data precede any exact theory, the non-critical divergence of h(x) and xi with respect to increasing solution supersaturation was nevertheless preliminarily established. Moreover, approximate models to account for behavior of both h(x) and xi are also presented.

Muramidase/chemistry↗

Fingerprinting Interactions between Proteins and Ligands for Facilitating Machine Learning in Drug Discovery

Molecular recognition is fundamental in biology, underpinning intricate processes through specific protein–ligand interactions. This understanding is pivotal in drug discovery, yet traditional experimental methods face limitations in exploring the vast chemical space. Computational approaches, notably quantitative structure–activity/property relationship analysis, have gained prominence. Molecular fingerprints encode molecular structures and serve as property profiles, which are essential in drug discovery. While two-dimensional (2D) fingerprints are commonly used, three-dimensional (3D) structural interaction fingerprints offer enhanced structural features specific to target proteins. Machine learning models trained on interaction fingerprints enable precise binding prediction. Recent focus has shifted to structure-based predictive modeling, with machine-learning scoring functions excelling due to feature engineering guided by key interactions. Notably, 3D interaction fingerprints are gaining ground due to their robustness. Various structural interaction fingerprints have been developed and used in drug discovery, each with unique capabilities. This review recapitulates the developed structural interaction fingerprints and provides two case studies to illustrate the power of interaction fingerprint-driven machine learning. The first elucidates structure–activity relationships in β2 adrenoceptor ligands, demonstrating the ability to differentiate agonists and antagonists. The second employs a retrosynthesis-based pre-trained molecular representation to predict protein–ligand dissociation rates, offering insights into binding kinetics. Despite remarkable progress, challenges persist in interpreting complex machine learning models built on 3D fingerprints, emphasizing the need for strategies to make predictions interpretable. Binding site plasticity and induced fit effects pose additional complexities. Interaction fingerprints are promising but require continued research to harness their full potential.

3D structural interaction fingerprints↗

Co-optimized machine-learned manifold models for large eddy simulation of turbulent combustion

Many modeling approaches in large eddy simulation (LES) of turbulent combustion employ a projection of the thermochemical state onto a low-dimensional manifold within state space to reduce the number of transported variables and hence computational cost. Flamelet-generated manifolds (FGM) is an example of a well-established, physics-based approach, but increasingly, principal component analysis (PCA) is being used as a data-driven method for generating manifold models. For both approaches, the nonlinear relationship between the location on the predefined manifold and the outputs of interest, such as reaction rates, can be tabulated or encoded in a neural network. This work proposes a new approach for manifold modeling that extends these existing approaches. A modified neural network structure simultaneously encodes the definition of the manifold variables, the nonlinear mapping, and the subfilter closure for LES. This allows all three of these aspects of the model to be co-optimized, generating a model from any source of combustion thermochemical state data. The manifold parameterizing variables are constrained to be linear combinations of species, as in FGM and PCA-based models, to aid in interpretability and implementation. For LES, subfilter variances of the manifold variables are also included as inputs. Two types of a priori analysis are performed to evaluate the new approach. In the first, the model is trained on data from one-dimensional premixed flames. In this case, the approach recovers the behavior of flamelet-based manifold approaches, and in fact slightly improves performance by identifying an optimized progress variable. The approach is also applied to data from direct numerical simulations of spherical ignition kernels in isotropic turbulence. For any specified manifold dimensionality, the new approach provides substantially lower prediction errors than a PCA-based model developed from the same data set. Additionally, the LES formulation of the new approach can provide accurate predictions for filtered reaction rates across a variety of filter widths.

97 MATHEMATICS AND COMPUTING↗

Self-consistent models for triaxial galaxies with flat rotation curves - The disk case

We examine the possibility of constructing scale-free triaxial logarithmic potentials self-consistently, using Schwarzschild's linear programing method. In particular, we explore the limit of nonaxisymmetric disks. In this case it is possible to reduce the problem to the self-consistent reconstruction of the disk surface density on the unit circle, a considerably simpler problem than the usual 2D or 3D one. Models with surface densities of the form Sigma = (x exp n + (y/q) exp n) exp - 1/n with n = 2 or 4 are investigated. We show that the complicated shapes of the 'boxlet' orbit families (which replace the box orbit family found in potentials with smooth cores) limit the possibility of building self-consistent models, though elliptical disks of axis ratio above 0.7 and a restricted range of boxier models can be constructed. This result relies on using sufficiently fine bins, smaller than the 10 deg bins commonly used in 2D or 3D investigations. It also indicates the need for caution in interpreting N-body models of triaxial halos in which the core of the potential is numerically smoothed.

Kuijken, Konrad↗

TPSAS-NF1676L-10206-DND

NASA is at fore-front in developing active and passive remote sensing technologies and unique capabilities towards space-based observations for understanding the complexities and interactions among Earth system components. The world is facing significant environmental challenges and a robust, integrated, and flexible system of observations and models are needed for understanding the short-and long term impact on the Earth system. A fundamental challenge for the coming decade is to ensure that space-based observations, analyses, better interpretive understanding, enhanced predictive models, broadened international community participation, and improved means for information assimilation and disseminations are well coordinated to realize the full economic, societal, and security benefit of Earth science. This presentation will provide an overview of NASA?s future vision for Earth science missions for global observations and the challenges associated in applying them for societal benefit.

Upendra N Singh↗

NASA's Future Earth Science Missions for Global Observations

NASA is at fore-front in developing active and passive remote sensing technologies and unique capabilities towards space-based observations for understanding the complexities and interactions among Earth system components. The world is facing significant environmental challenges and a robust, integrated, and flexible system of observations and models are needed for understanding the short-and long term impact on the Earth system. A fundamental challenge for the coming decade is to ensure that space-based observations, analyses, better interpretive understanding, enhanced predictive models, broadened international community participation, and improved means for information assimilation and disseminations are well coordinated to realize the full economic, societal, and security benefit of Earth science. This presentation will provide an overview of enabling active and passive remote sensing technologies and techniques, NASA's future vision for Earth science missions for global observations, and the challenges associated in applying them for societal benefit.

Upendra N Singh↗

TPSAS-NF1676L-27491-DND

NASA is at fore-front in developing lidar technologies and unique active/passive remote sensing capabilities towards space-based observations for understanding the complexities and interactions among Earth system components. The world is facing significant environmental challenges and a robust, integrated, and flexible system of observations and models are needed for understanding the short-and long term impact on the Earth system. A fundamental challenge for the coming decade is to ensure that space-based observations, analyses, better interpretive understanding, enhanced predictive models, broadened international community participation, and improved means for information assimilation and disseminations are well coordinated to realize the full economic, societal, and security benefit of Earth science. This presentation will provide an overview of enabling lidar technologies and techniques from ground and space towards NASA's future vision for Earth science missions for global observations, and the challenges associated in applying them for societal benefit.

Upendra N Singh↗

Search for resonant pair production of Higgs bosons in the $\textrm{b}\overline{\textrm{b}}\textrm{b}\overline{\textrm{b}}$ final state using large-area jets in proton-proton collisions at $\sqrt{s}$ = 13 TeV

A search is presented for the resonant production of a pair of standard model-like Higgs bosons using data from proton-proton collisions at a centre-of-mass energy of 13 TeV, collected by the CMS experiment at the CERN LHC in 2016–2018, corresponding to an integrated luminosity of 138 fb −1 . The final state consists of two b quark-antiquark pairs. The search is conducted in the region of phase space where at least one of the pairs is highly Lorentz-boosted and is reconstructed as a single large-area jet. The other pair may be either similarly merged or resolved, the latter reconstructed using two b-tagged jets. The data are found to be consistent with standard model processes and are interpreted as 95% confidence level upper limits on the product of the cross sections and the branching fractions of the spin-0 radion and the spin-2 bulk graviton that arise in warped extradimensional models. The limits set are in the range 9.74–0.29 fb and 4.94–0.19 fb for a narrow radion and a graviton, respectively, with masses between 1 and 3 TeV. For a radion and for a bulk graviton with widths 10% of their masses, the limits are in the range 12.5–0.35 fb and 8.23–0.23 fb, respectively, for the same masses. These limits result in the exclusion of a narrow-width graviton with a mass below 1.2 TeV, and of narrow and 10%-width radions with masses below 2.6, and 2.9 TeV, respectively.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

White Paper: Scalable Digital Twin Capabilities for Aging and Surveillance of Engineered Systems

This white paper presents a multi-year initiative to develop practical, secure, and scalable digital twin capabilities for engineered systems in aging and surveillance contexts—an approach pioneered at the National Nuclear Security Administration (NNSA) Lawrence Livermore National Laboratory (LLNL) that maps directly onto the needs and ambitions of the Navy for ship- and fleet-level digital twins. LLNL’s work in building part- and process-level digital twins for advanced manufacturing, with a vision to scale up to entire factory floors and, ultimately, enterprise-wide digital twins, offers an adaptable pathway for the Navy as it seeks to modernize lifecycle management, readiness, and predictive maintenance across ships and fleets. For our application, we integrate physics-based modeling with automated data ingestion, processing, and AI-driven calibration, creating hybrid models that are both interpretable and data responsive. We modernized legacy workflows, established centralized data infrastructure, automated experimental pipelines, and demonstrated end-to-end coupling of accelerated aging data with finite element simulations via optimization and surrogate modeling. The result is a generalizable framework that supports part-level digital twins today and lays the groundwork for future system-level twins suitable for Navy applications.

36 MATERIALS SCIENCE↗

TPSAS-NF1676L-30074-DND

NASA is at forefront in developing unique active/passive remote sensing capabilities towards space-based observations for understanding the complexities and interactions among Earth system components. The world is facing significant environmental challenges and a robust, integrated, and flexible system of observations and models are needed for understanding the short-and long-term impact on the Earth system. A fundamental challenge for the coming decade is to ensure that space-based observations, analyses, better interpretive understanding, enhanced predictive models, broadened international community participation, and improved means for information assimilation and disseminations are well coordinated to realize the full economic, societal, and security benefit of Earth science. This presentation will provide an overview of active/passive remote sensing technologies and techniques from ground and space towards NASA’s future vision for Earth science missions for global observations, and the challenges associated in applying them for societal benefit.

Upendra N Singh↗

Theoretical interpretation of the HEAO-3 observations of Cygnus X-3 under the HEAO-3 Guest Investigator Program

A model of the galactic X-ray source Cygnus X-3 (Cyg X-3) is presented which deviates from previous models by positing that the X-rays originate in a jet rather than a binary system consiting of an ordinary star and a collapsed object. In the new model, the 4.8 hour period of Cyg X-3 is caused by variable absorption which occurs as the jet precesses. The primary role of the accretion disk corona (ADC) in modulating Cyg X-3 radiation is to make the observed intensity of a blob of material in a jet appear dimmer by absorption. The needed derivation of the positional dependence of the density of the ADC is freed of some complications by assuming that only the inner regions of the disk are precessing, with a period shorter than 4.8 hours. This model permits the secondary star to be a supergiant, as indicated by the luminosity of the system. The model is especially helpful in interpreting production of radio outbursts and very high energy gamma rays.

Marscher, Alan P.↗

VESIcal Part I: An open-source thermodynamic model engine for mixed volatile (H2O-CO2) solubility in silicate melts

Thermodynamics has been fundamental to the interpretation of geologic data and modeling of geologic systems for decades. However, more recent advancements in computational capabilities and a marked increase in researchers’ accessibility to computing tools has outpaced the functionality and extensibility of currently available modeling tools. Here we present VESIcal (Volatile Equilibria and Saturation Identification calculator): the first comprehensive modeling tool for H 2 O, CO 2 , and mixed (H 2 O-CO 2 ) solubility in silicate melts that: a) allows users access to seven commonly used models, plus easy inter-comparison between models; b) provides universal functionality for all models (e.g., functions for calculating saturation pressures, degassing paths, etc.); c) can process large datasets (1,000’s of samples) automatically; d) can output computed data into an excel spreadsheet for simple post-modeling analysis; e) integrates advanced plotting capabilities directly within the tool; and f) provides all of these within the framework of a python library, making the tool extensible by the user and allowing any of the model functions to be incorporated into any other code capable of calling python. The tool is presented within this manuscript, which is a Jupyter notebook containing worked examples accessible to python users with a range of skill levels. The basic functions of VESIcal can also be access via a web app (https://vesical.anvil.app). The VESIcal python library is open-source and available for download at https://github.com/kaylai/VESIcal.

K. Iacovino↗

Electronic Spectroscopy and Dissociation Dynamics of Vinyl-Substituted Criegee Intermediates: 2-Butenal Oxide and Comparison with Methyl Vinyl Ketone Oxide and Methacrolein Oxide Isomers

The 2-butenal oxide Criegee intermediate [(CH 3 CH=CH)CHOO], an isomer of the four-carbon unsaturated Criegee intermediates derived from isoprene ozonolysis, is characterized on its first π * ←π electronic transition and by the resultant dissociation dynamics to O ( 1 D) + 2-butenal [(CH 3 CH=CH)CHO] products. The electronic spectrum of 2-butenal oxide under jet-cooled conditions is observed to be broad and unstructured with peak absorption at 373 nm, spanning to half maxima at 320 and 420 nm, and in good accord with the computed vertical excitation energies and absorption spectra obtained for its lowest energy conformers. The distribution of total kinetic energy released to products is ascertained through velocity map imaging of the O ( 1 D) products. About half of the available energy, deduced from the theoretically computed asymptotic energy, is accommodated as internal excitation of the 2-butenal fragment. A reduced impulsive model is introduced to interpret the photodissociation dynamics, which accounts for the geometric changes between 2-butenal oxide and the 2-butenal fragment, and vibrational activation of associated modes in the 2-butenal product. Application of the reduced impulsive model to the photodissociation of isomeric methyl vinyl ketone oxide reveals greater internal activation of the methyl vinyl ketone product arising from methyl internal rotation and rock, which is distinctly different than the dissociation dynamics of 2-butenal oxide or methacrolein oxide.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Size-density relations in dark clouds: Non-LTE effects

One of the major goals of molecular astronomy has been to understand the physics and dynamics of dense interstellar clouds. Because the interpretation of observations of giant molecular clouds is complicated by their very complex structure and the dynamical effects of star formation, a number of studies have concentrated on dark clouds. Leung, Kutner and Mead (1982) (hereafter LKM) and Myers (1983), in studies of CO and NH3 emission, concluded that dark clouds exhibit significant correlations between linewidth and cloud radius of the form delta v varies as R(0.5) and between mean density and radius of the form n varies as R(-1), as originally suggested by Larson (1981). This result suggests that these objects are in virial equilibrium. However, the mean densities inferred from the CO data of LKM are based on an local thermodynamic equilibrium (LTE) analysis of their 13CO data. At the very low mean densities inferred by LKM for the larger clouds in their samples, the assumption of LTE becomes very questionable. As most of the range in R in the density-size correlation comes from the clouds observed in CO, it seems worthwhile to examine how non-LTE effects will influence the derived densities. One way to assess the validity of LTE-derived densities is to construct cloud models and then to interpret them in the same way as the observed data. Microturbulent models of inhomogeneous clouds of varying central concentration with the linewidth-size and mean density-size relations found by Myers show sub-thermal excitation of the 13CO line in the larger clouds, with the result that LTE analysis considerbly underestimates the actual column density. A more general approach which doesn't require detailed modeling of the clouds is to consider whether the observed T sub R*(13CO)/T sub R*(12CO) ratios in the clouds studied by LKM are in the range where the LTE-derived optical depths (and hence column densities) can be seriously in error due to sub-thermal excitation of the 13CO molecule.

Maloney, P.↗