SEARCH · Search NASA

Results for “interpretable models”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 199 records · Page 11

All sources NASA NTRS DOE OSTI

NASA NTRS1984

Finite-difference fluid dynamics computer mathematical models for the design and interpretation of experiments for space flight

Numerical methods are used to design a spherical baroclinic flow model experiment of the large scale atmosphere flow for Spacelab. The dielectric simulation of radial gravity is only dominant in a low gravity environment. Computer codes are developed to study the processes at work in crystal growing systems which are also candidates for space flight. Crystalline materials rarely achieve their potential properties because of imperfections and component concentration variations. Thermosolutal convection in the liquid melt can be the cause of these imperfections. Such convection is suppressed in a low gravity environment. Two and three dimensional finite difference codes are being used for this work. Nonuniform meshes and implicit iterative methods are used. The iterative method for steady solutions is based on time stepping but has the options of different time steps for velocity and temperature and of a time step varying smoothly with position according to specified powers of the mesh spacings. This allows for more rapid convergence. The code being developed for the crystal growth studies allows for growth of the crystal as the solid-liquid interface. The moving interface is followed using finite differences; shape variations are permitted. For convenience in applying finite differences in the solid and liquid, a time dependent coordinate transformation is used to make this interface a coordinate surface.

Roberts, G. O.↗

DOE OSTI2025

Explainable Machine Learning for Functional Data

Black-box machine learning models are recognized as useful tools for prediction applications, but the algorithmic complexity of some models causes interpretation challenges. Explainability methods have been proposed to provide insight into these models, but there is little research focused on supervised modeling with functional data inputs. We argue that, especially in applications of high consequence, it is important to explicitly model the functional dependence in a black-box analysis to not obscure or misrepresent patterns in explanations. As such, we propose the V ariable importance E xplainable E lastic S hape A nalysis (VEESA) pipeline for training supervised machine learning models with functional inputs. The pipeline is an analysis process that includes the data preprocessing, modeling, and post-hoc explanations. The preprocessing is done using elastic functional principal components analysis, which accounts for vertical and horizontal variability in functional data and, ultimately, allows for explanations in the original data space that identify the important functional variability without bias due to correlated variables. Here, we demonstrate the pipeline on two high-consequence applications: explosives classification for national security and inkjet printer identification in forensic science. The applications exhibit the VEESA pipeline’s ability to provide an understanding of the characteristics of the functional data useful for prediction. Code for implementing the pipeline is available in the veesa R package (and supplemental python code).

Elastic Shape Analysis↗

DOE OSTI2026

Overview of IMPACT Data Acquisition System and Data Reduction Process

This report documents the development of the data acquisition system (DAS) and data reduction methodologies for the Irradiated Material Property Accelerated Characterization Test (IMPACT) experiment at the Advanced Test Reactor (ATR). The IMPACT experiment is designed to enable in-pile measurement of thermal conductivity in metallic nuclear fuels, specifically U-10Zr, using an instrumented thermal conductivity probe. The DAS supports both passive temperature monitoring and active thermal interrogation of the probe through controlled AC and DC excitation. Significant modifications to laboratory-scale systems were required to accommodate the higher resistance paths associated with the in-pile application. Custom electronics and relay-controlled measurement sequencing were developed to enable the measurement and sufficient power delivery to the sensing region. A reduced-order, axisymmetric thermal model based on the thermal quadrupoles method is presented to support data interpretation. This model enables efficient evaluation of transient heat transfer behavior and facilitates solution of the inverse problem required to extract thermal properties from measured signals. Multiple boundary condition formulations are discussed to address varying experimental time scales and geometries. Additionally, machine learning techniques are introduced to support data reduction and improve confidence in inverse solutions. Convolutional neural networks are applied to identify the presence of gas gaps and other evolving geometric features that significantly impact thermal response during irradiation. These efforts contribute to the broader integration of digital twin frameworks and real-time modeling capabilities within the Advanced Fuels Campaign.

46 - INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AN↗

NASA NTRS1980

Wheat forecast economics effect study

A model to assess the value of improved information regarding the inventories, productions, exports, and imports of crop on a worldwide basis is discussed. A previously proposed model is interpreted in a stochastic control setting and the underlying assumptions of the model are revealed. In solving the stochastic optimization problem, the Markov programming approach is much more powerful and exact as compared to the dynamic programming-simulation approach of the original model. The convergence of a dual variable Markov programming algorithm is shown to be fast and efficient. A computer program for the general model of multicountry-multiperiod is developed. As an example, the case of one country-two periods is treated and the results are presented in detail. A comparison with the original model results reveals certain interesting aspects of the algorithms and the dependence of the value of information on the incremental cost function.

Mehra, R. K.↗

NASA NTRS1999

Remote Thermal IR Spectroscopy of our Solar System

Indirect methods to detect extrasolar planets have been successful in identifying a number of stars with companion planets. No direct detection of an extrasolar planet has yet been reported. Spectroscopy in the thermal infrared region provides a potentially powerful approach to detection and characterization of planets and planetary systems. We can use knowledge of our own solar system, its planets and their atmospheres to model spectral characteristics of planets around other stars. Spectra derived from modeling our own solar system seen from an extrasolar perspective can be used to constrain detection strategies, identification of planetary class (terrestrial vs. gaseous) and retrieval of chemical, thermal and dynamical information. Emission from planets in our solar system peaks in the thermal infrared region, approximately 10 - 30 microns, substantially displaced from the maximum of the much brighter solar emission in the visible near 0.5 microns. This fact provides a relatively good contrast ratio to discriminate between stellar (solar) and planetary emission and optimize the delectability of planetary spectra. Important molecular constituents in planetary atmospheres have rotational-vibrational spectra in the thermal infrared region. Spectra from these molecules have been well characterized in the laboratory and studied in the atmospheres of solar system planets from ground-based and space platforms. The best example of such measurements are the studies with Fourier transform spectrometers, the Infrared Interferometer Spectrometers (IRIS), from spacecraft: Earth observed from NIMBUS 8, Mars observed from Mariner 9, and the outer planets observed from Voyager spacecraft. An Earth-like planet is characterized by atmospheric spectra of ozone, carbon dioxide, and water. Terrestrial planets have oxidizing atmospheres which are easily distinguished from reducing atmospheres of gaseous giant planets which lack oxygen-bearing species and are characterized by spectra containing hydrocarbons such as methane and ethane. Spectroscopic information on extrasolar planets thus can permit their classification. Spectra and spectral lines contain information on the temperature structure of the atmosphere. Line and band spectra can be used to identify the molecular constituents and retrieve species abundances, thereby classifying and characterizing the planet. At high enough spectral resolution characteristic planetary atmospheric dynamics and unique phenomena such as failure of local thermodynamic equilibrium can be identified. Dynamically induced effects such as planetary rotation and orbital velocity shift and change the shape of spectral features and must be modeled in detailed spectral studies. We will use our knowledge of the compositional, thermal and dynamical characteristics of planetary atmospheres in our own solar system to model spectra observed remotely on similar planets in extrasolar planetary systems. We will use a detailed radiative transfer and beam integration program developed for the modeling and interpretation of thermal infrared spectra measured from nearby planet planets to generate models of an extra-solar "Earth" and "Jupiter". From these models we will show how key spectral features distinguish between terrestrial and gaseous planets, what information can be obtained with different spectral resolution, what spectral features can be used to search for conditions for biogenic activity, and how dynamics and distance modify the observed spectra. We also will look at unique planetary phenomena such as atmospheric lasing and discuss their utility as probes for detection and identification of planets. Results of such studies will provide information to constrain design for instrumentation needed to directly detect extrasolar planets.

Kostiuk, Theodor↗

DOE OSTI2025

pyTCR: A tropical cyclone rainfall model for python

pyTCR is a climatology software package developed in the Python programming language. It integrates the capabilities of several legacy physical models and increases computational efficiency to allow rapid estimation of tropical cyclone (TC) rainfall consistent with the large-scale environment. Specifically, pyTCR implements a horizontally distributed and vertically integrated model [Zhu et al., 2013] for simulating rainfall driven by TCs. Along storm tracks, rainfall is estimated by computing the cross-boundary-layer, upward water vapor transport caused by different mechanisms including frictional convergence, vortex stretching, large-scale baroclinic effect (i.e., wind shear), topographic forcing, and radiative cooling [Lu et al., 2018]. The package provides essential functionalities for modeling and interpreting spatio-temporal TC rainfall data. pyTCR requires a limited number of model input parameters, making it a convenient and useful tool for analyzing rainfall mechanisms driven by TCs. To sample rare (most intense) rainfall events that are often of great societal interest, pyTCR adapts and leverages outputs from a statistical-dynamical TC downscaling model [Lin et al., 2023] capable of rapidly generating a large number of synthetic TCs given a certain climate. As a result, pyTCR significantly reduces computational effort and improves the efficiency in capturing extreme TC rainfall events at the tail of the distributions from limited datasets. Furthermore, the TC downscaling model is forced entirely by large-scale environmental conditions from reanalysis data or coupled General Circulation Models (GCMs), simplifying the projection of TC-induced rainfall and wind speed under future climate using pyTCR. Finally, pyTCR can be coupled with hydrological and wind models to assess risks associated with independent and compound events (e.g., storm surges and freshwater flooding).

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1985

A model for the proton spectrum and cosmic ray anisotropy

The problem of the origin of the cosmic rays is still uncertain. As a theory, it should explain the support of particles and energy, the mechanism of acceleration and propagation as well as some important features obtained directly from cosmic ray experiments, such as the power spectrum and the knee. There are two kinds of models for interpreting the knee of the cosmic ray spectrum. One is the leaky box model. Another model suggests that the cut-off rigidity of the main sources causes the knee. The present paper studies the spectrum and the anisotropy of cosmic rays in an isotropic diffuse model with explosive discrete sources in an infinite galaxy.

Xu, C.↗

NASA NTRS2012

GRB 120422A: a Low-Luminosity Gamma-Ray Burst Driven by a Central Engine

GRB 120422A is a low-luminosity gamma-ray burst (GRB) associated with a bright supernova, which distinguishesitself by its relatively short T(sub 90) (approximately 5 s) and an energetic and steep-decaying X-ray tail. We analyze the Swift BurstAlert Telescope and X-ray Telescope data and discuss the physical implications. We show that the steep declineearly in the X-ray light curve can be interpreted as the curvature tail of a late emission episode around 58-86 s,with a curved instantaneous spectrum at the end of the emission episode. Together with the main activity in thefirst 20 s and the weak emission from 40 s to 60 s, the prompt emission is variable, which points to a centralengine origin in contrast to a shock-breakout origin, which is used to interpret some other nearby low-luminosity supernova GRBs. Both the curvature effect model and interpreting the early shallow decay as the coasting externalforward shock emission in a wind medium provide a constraint on the bulk Lorentz factor to be around several.Comparing the properties ofGRB 120422A and other supernova GRBs,we find that themain criterion to distinguish engine-driven GRBs from shock-breakout GRBs is the time-averaged -ray luminosity. Engine-driven GRBs likelyhave a luminosity above approximately 10(sup 48) erg s(sup -1).

engine↗

DOE OSTI2026

Deciphering the Solvation Structure of Aqueous ZnCl 2 Solutions from X-ray Absorption Spectra Using the Interpretable Graph Neural Network

Machine learning (ML) provides powerful pathways for predicting spectroscopic observables from atomic structures, but its broader impact depends on making model predictions interpretable in terms of physical and chemical principles. Here, we introduce a physics-guided graph neural network (GNN) model that predicts Zn K-edge X-ray spectroscopy (XAS) spectra of aqueous ZnCl 2 solutions. Training data are generated from ab initio XAS calculations on molecular dynamics snapshots obtained using a machine learning interatomic potential. The GNN reproduces experimental spectra across concentrations from dilute (<0.1 m) to highly concentrated (30 m, “water-in-salt”) regimes and scales efficiently to large, disordered liquid systems beyond the reach of conventional ab initio approaches. Gradient-based attribution analysis reveals that the model learns physically meaningful structure-spectrum relationships. Ligand-specific attributions reflect orbital hybridization patterns and the origin of the excitations derived from the density functional theory. Bond-length attributions recover spectral shifts consistent with multiple-scattering theory. Finally, this work bridges data-driven prediction with electronic-structure theory, establishing a general paradigm for interpretable ML that links atomic structure, electronic structure, and spectroscopic observables.

25 ENERGY STORAGE↗

DOE OSTI2026

A Unified Workflow for Sensitivity-Based Kinetic Analysis in Microkinetic Models

Degrees of rate control (DRC), apparent activation energies, and apparent reaction orders are established local sensitivity diagnostics for interpreting microkinetic models, but applying them routinely to large mechanisms often requires substantial reaction-specific bookkeeping, perturbation design, and postprocessing. Here, in this study, we present a unified derivative-based workflow that evaluates these quantities from a single compiled reaction-network model and target-rate definition. For any user-provided microkinetic model, the workflow compiles the mechanism into stoichiometrically consistent mass-action rate equations, solves the surface dynamics, and uses automatic differentiation to compute sensitivities with respect to rate constants, temperature, and gas partial pressures. By combining their calculations in the same framework, the workflow clearly demonstrates the relationships between different DRCs and the apparent activation energy. Using existing examples of propylene partial oxidation and methane oxidation on Pd(100), we verify expected transient redistribution of rate control, distinguish net Campbell DRCs from one-sided directional sensitivities, and show how apparent activation energy can be reconstructed either from one-sided DRCs or from state-based DRCs while critical mechanistic insights are obtained consistently. In the methane oxidation case, a pathway-subset test further illustrates how a simplified mechanism preserves key kinetic signatures of a full model, showing the potential of our user-friendly tool for model construction beyond kinetic analysis.

36 MATERIALS SCIENCE↗

NASA NTRS1980

Interpretation of postseismic deformation with a viscoelastic relaxation model

A viscoelastic relaxation model is used to interpret postseismic surface deformation for the large-magnitude Aleutian earthquakes of 1957 and 1965. The lithosphere and asthenosphere are modeled as elastic solids with an anomalous viscoelastic inclusion below the island arc volcanoes. It is found that the observed postseismic surface deformation is a corollary of the known island arc structure. A satisfactory fit to the uplift following the 1957 earthquake is found for a viscoelastic volume with 80-km width extending from a depth of 50 to 200 km.

Wahr, J.↗

DOE OSTI2026

A procedure for rule extraction from a Self-Organising plasma disruption predictor for JET

In a previous paper, a Self-Organizing Map had proven to be able to identify the regions of the plasma operative space characterizing the pre-disruptive phase at JET without relying on any a priori information. One of the strengths of this disruption predictor lies in its inherent self-organization capability. The Self-Organizing Map discovers non-trivial relationships and captures the complicated interplay of device diagnostics on the internal plasma states directly from the experimental data. Moreover, the provided model allows the visualization of high-dimensional plasma parameters and facilitates easy interrogation of the model to understand the reasons behind its correlations. In this paper, an additional step is taken towards the interpretability of models for predicting disruptions by training a Decision Tree to classify the plasma states according to the interpretation provided by the Self-Organizing Map (stable or at high risk of disruptions). The Decision tree provides a set of rules which describe the transition of the plasma towards the pre-disruptive phase as visualized in the Self-Organizing Map. The obtained rules for the database explored in the study identify four regions in the map, two of which are at risk of disruption. These regions correspond to partitions of a 3D space based on the peaking factors of the core and divertor radiation, as well as the Locked Mode. The agreement between the Self-Organizing Map answers and the rules supplied by the Decision Tree is confirmed by the comparison of the performance exhibited by the two models in the prediction of disruptions.

Setzu, Samuele [Univ. of Cagliari, Monserrato, Cag↗

NASA NTRS2011

The 5 Hour Pulse Period and Broadband Spectrum of the Symbiotic X-Ray Binary 3A 1954+319

We present an analysis of the highly variable accreting X-ray pulsar 3A 1954+319 using 2005-2009 monitoring data obtained with INTEGRAL and Swift. This considerably extends the pulse period history and covers flaring episodes in 2005 and 2008. In 2006 the source was identified as one of only a few known symbiotic X-ray binaries, Le" systems composed of a neutron star accreting from the inhomogeneous medium around an M-giant star. The extremely long pulse period of approximately 5.3 h is directly visible in the 2008 INTEGRAL-ISGRI outburst light curve. The pulse profile is double peaked and not significantly energy dependent. During the outburst a strong spin-up of -1.8 x 10(exp -4) h h(exp -1) occurred. Between 2005 and 2008 a long term spin-down trend of 2.1 x 10(exp -5) h h(exp -1) was observed for the first time for this source. The 3-80 keV pulse peak spectrum of 3A 1954+319 during the 2008 flare could be well described by a thermal Comptonization model. We interpret the results within the framework of a recently developed quasi-spherical accretion model for symbiotic X-ray binaries.

Marcu, Diana M.↗

NASA NTRS1991

Relationships of soil, grass, and bedrock over the Kaweah serpentine melange through spectral mixture analysis of AVIRIS data

A linear mixing model is used to model the spectral variability of an AVIRIS scene from the western foothills of the Sierra Nevada and calibrate these radiance data to reflectance. Five spectral endmembers from the AVIRIS data, plus an ideal 'shade' endmember were required to model the continuum reflectance of each pixel in the image. Three of the endmembers were interpreted to model the surface constituents green vegetation, dry grass, and illumination. These are the main transient surface constituents that are expected to change with shifts in land use or climatic influences and viewing conditions ('shade' only). The spectral distinction between the other three endmembers is very small, yet the spatial distributions are coherent and interpretable. These distributions cross anthropogenic and vegetation boundaries and are best interpreted as different soil types. Comparison of the fraction images to the bedrock geology maps indicates that substrate composition must be a factor contributing to the spectral properties of these endmembers. Detailed examination of the reflectance spectra of the three soil endmembers reveals that differences in the amount of ferric and ferrous iron and/or organic constituents in the soils is largely responsible for the differences in spectral properties of these endmembers.

Mustard, John F.↗

NASA NTRS2021

Integrating continuous atmospheric boundary layer and tower-based flux measurements to advance understanding of land-atmosphere interactions

The atmospheric boundary layer mediates the exchange of energy, matter, and momentum between the land surface and the free troposphere, integrating a range of physical, chemical, and biological processes and is defined as the lowest layer of the atmosphere (ranging from a few meters to 3 km). In this review, we investigate how continuous, automated observations of the atmospheric boundary layer can enhance the scientific value of co-located eddy covariance measurements of land-atmosphere fluxes of carbon, water, and energy, as are being made at FLUXNET sites worldwide. We highlight four key opportunities to integrate tower-based flux measurements with continuous, long-term atmospheric boundary layer measurements: (1) to interpret surface flux and atmospheric boundary layer exchange dynamics and feedbacks at flux tower sites, (2) to support flux footprint modelling, the interpretation of surface fluxes in heterogeneous and mountainous terrain, and quality control of eddy covariance flux measurements, (3) to support regional-scale modeling and upscaling of surface fluxes to continental scales, and (4) to quantify land-atmosphere coupling and validate its representation in Earth system models. Adding a suite of atmospheric boundary layer measurements to eddy covariance flux tower sites, and supporting the sharing of these data to tower networks, would allow the Earth science community to address new emerging research questions, better interpret ongoing flux tower measurements, and would present novel opportunities for collaborations between FLUXNET scientists and atmospheric and remote sensing scientists.

Manuel Helbig↗

DOE OSTI2025

2D reactive transport model of shale chemical weathering and biogeochemical fluxes along a mountainous hillslope, East River Watershed, Colorado: Input files and simulation results

This data package contains input files and simulation results for a two-dimensional (2D) reactive transport model used to quantitatively analyze the coupled hydrological and biogeochemical processes governing shale weathering and associated biogeochemical fluxes under realistic environmental conditions in the high-elevation East River Watershed. These data support the conclusions presented in Stolze et al. (Water Resources Research, under review), "Model-based interpretation of solute exports and carbon partitioning during shale weathering in a mountainous hillslope". The model simulates atmospheric-subsurface gas exchange, subsurface water flow, and shale weathering processes under dynamic, year-scale conditions along a shale-underlain hillslope located in the East River watershed. The simulations were performed using the PFLOTRAN flow and reactive transport code and executed on the Perlmutter supercomputer to leverage its large-scale parallel computing capabilities. The data package contains two zipped folders, "model_input_files" and "simulation_results", and one readme.txt file. "model_input_files" contains the necessary input files to run the calibrated base-base model presented in Stolze et al. (Water Resources Research, under review). "simulation_results" contains a single hdf5 file ("Output_2D_hillslope_model.h5") which includes the results of simulation performed using the base-case model. This file can be opened with HDFView 3.1.4, Python, or MATLAB. "readme.txt" contains relevant information about the base-case model and provides guidelines on how to run the associated input files provided in the folder "model_input_files". Furthermore, readme.txt provides information regarding the model results provided in "Output_2D_hillslope_model.h5" such as matrix dimensionality and output units. Field datasets used to evaluate model performance were collected at three monitoring wells located along a hillslope transect (PLM1, PLM2, and PLM3). Dissolved ion concentration data were collected from November 2016 to October 2021 for Ca, Mg, DIC, Na, K, SO4 (Dong et al., 2025 - dic_npoc_data_2014_2024.zip - DOI:10.15485/1660459; Williams et al., 2025 - anion_data_2014_2024.zip - DOI:10.15485/1668054; Dong et al., 2025 - cation_data_2014_2024.zip - DOI:10.15485/1668055). Note that we used the files named er_PLM1_xx_yy, er_PLM2_xx_yy, and er_PLM3_xx_yy where xx stands for the name of the aqueous species and yy stands for the depth where the measurements were performed. Soil water content ([0 - 1] m) and water table depth were collected from November 2016 to October 2021 (Wan et al., 2024 - Dynamic_water_table__depthsFig2b.csv and Soil_water_content_Fig4e.csv - DOI:10.15485/2322567). Gaseous CO2 concentration were collected from October 2020 to December 2021(Wan et al., 2024 - Soil_CO2_concentrations_Fig4h.csv - DOI:10.15485/2322567) Gaseous CO2 flux from the subsurface to the atmosphere were collected in the vicinity of PLM2 from October 2019 to May 2022 (Wu et al., 2025). Soil microbial biomass concentration was measured from August 2016 to June 2017 (Sorensen et al., 2019 - 2017_East_River_Pumphouse_Microbial_Biomass__1_.csv - DOI:10.15485/1577267) All field data are published as CSV files compatible with Microsoft Excel, MATLAB, and Python, or as text files. The coordinates of the monitoring wells and the CO2(g) flux sensor in the coordinate system WGS84 are: -PLM1: [38.9197710 ; -106.9492750] -PLM2: [38.9201580 ; -106.9487170] -PLM3: [38.9207843 ; -106.9483668] -PLM4: 38.9210060 ; -106.9479528] -CO2(g) flux sensor: [38.9199180 ; -106.9489906] ------------------------------------------------------------------------------------------- This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231. This research used resources of the National Energy Research Scientific Computing Center (NERSC), a Department of Energy User Facility using NERSC award BER-ERCAP 23980, BER-ERCAP 28550, and BER-ERCAP 33789.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1993

Relationships of soil, grass, and bedrock over the Kaweah Serpentinite Melange through spectral mixture analysis of AVIRIS data

Mustard, John F.↗

NASA NTRS1990

Predicted distribution of visible and near-infrared radiant flux above and below a transmittant leaf

The effects are studied analytically of leaf size, leaf height, and background reflectance on the upward and downward radiant flux (RF) of a leaf. The leaf is horizontal and isotropically scattering in the computer model which examines the light environment in three regions about the leaf. The spectral properties of the leaf are based on measurements of the big-leaf maple, and the model is interpreted in terms of relative RF which is defined as a percentage of the total light in the model. The results demonstrate the dependence of upward relative RF on the light's wavelength and background reflectance with large variations in the NIR. Brightness varies directly with distance from background with maximum brightness achieved at lower heights for smaller leaves. These and other results suggest that NIR canopy reflectance due to leaves is highly dependent on the background reflectance.

Roberts, Dar A.↗