Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Model interpretation”

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 325 records · Page 18

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.↗

Search for electroweak production of charginos and neutralinos in multilepton final states in proton-proton collisions at $\sqrt{s}=$ 13 TeV

Results are presented from a search for the direct electroweak production of charginos and neutralinos in signatures with either two or more leptons (electrons or muons) of the same electric charge, or with three or more leptons, which can include up to two hadronically decaying tau leptons. The results are based on a sample of proton-proton collision data collected at $ \sqrt{s}=13 $ TeV, recorded with the CMS detector at the LHC, corresponding to an integrated luminosity of 35.9 fb$^{−1}$. The observed event yields are consistent with the expectations based on the standard model. The results are interpreted in simplified models of supersymmetry describing various scenarios for the production and decay of charginos and neutralinos. Depending on the model parameters chosen, mass values between 180 GeV and 1150 GeV are excluded at 95% CL. These results significantly extend the parameter space probed for these particles in searches at the LHC. In addition, results are presented in a form suitable for alternative theoretical interpretations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

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↗

To use or not to use synthetic stellar spectra in population synthesis models?

ABSTRACT Stellar population synthesis (SPS) models are invaluable to study star clusters and galaxies. They provide means to extract stellar masses, stellar ages, star formation histories, chemical enrichment, and dust content of galaxies from their integrated spectral energy distributions, colours, or spectra. As most models, they contain uncertainties that can hamper our ability to model and interpret observed spectra. This work aims at studying a specific source of model uncertainty: the choice of an empirical versus a synthetic stellar spectral library. Empirical libraries suffer from limited coverage of parameter space, while synthetic libraries suffer from modelling inaccuracies. Given our current inability to have both ideal stellar-parameter coverage with ideal stellar spectra, what should one favour: better coverage of the parameters (synthetic library) or better spectra on a star-by-star basis (empirical library)? To study this question, we build a synthetic stellar library mimicking the coverage of an empirical library, and SPS models with different choices of stellar library tailored to these investigations. Through the comparison of model predictions and the spectral fitting of a sample of nearby galaxies, we learned that predicted colours are more affected by the coverage effect than the choice of a synthetic versus empirical library; the effects on predicted spectral indices are multiple and defy simple conclusions; derived galaxy ages are virtually unaffected by the choice of the library, but are underestimated when SPS models with limited parameter coverage are used; metallicities are robust against limited HRD coverage, but are underestimated when using synthetic libraries.

Coelho, Paula R. T.↗

Searches for new phenomena in events with two leptons, jets, and missing transverse momentum in 139 fb –1 of √s = 13 TeV $pp$ collisions with the ATLAS detector

Searches for new phenomena inspired by supersymmetry in final states containing an e + e – or μ + μ – pair, jets, and missing transverse momentum are presented. These searches make use of proton–proton collision data with an integrated luminosity of 139 fb –1 , collected during 2015–2018 at a centre-of-mass energy √s = 13 TeV by the ATLAS detector at the Large Hadron Collider. Two searches target the pair production of charginos and neutralinos. One uses the recursive-jigsaw reconstruction technique to follow up on excesses observed in 36.1 fb –1 of data, and the other uses conventional event variables. The third search targets pair production of coloured supersymmetric particles (squarks or gluinos) decaying through the next-to-lightest neutralino ($\tilde{χ}$$^{0}_{2}$) via a slepton ($\tilde{ℓ}$) or Z boson into ℓ + ℓ – $\tilde{χ}$$^{0}_{1}$ , resulting in a kinematic endpoint or peak in the dilepton invariant mass spectrum. The data are found to be consistent with the Standard Model expectations. Results are interpreted using simplified models and exclude masses up to 900 GeV for electroweakinos, 1550 GeV for squarks, and 2250 GeV for gluinos.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Searches for electroweak production of supersymmetric particles with compressed mass spectra in s = 13 TeV p p collisions with the ATLAS detector

This paper presents results of searches for electroweak production of supersymmetric particles in models with compressed mass spectra. The searches use 139/fb of sqrt(s) = 13 TeV proton-proton collision data collected by the ATLAS experiment at the Large Hadron Collider. Events with missing transverse momentum and two same-flavor, oppositely charged, low transverse momentum leptons are selected, and are further categorized by the presence of hadronic activity from initial-state radiation or a topology compatible with vector-boson fusion processes. The data are found to be consistent with predictions from the Standard Model. The results are interpreted using simplified models of R-parity-conserving supersymmetry in which the lightest supersymmetric partner is a neutralino with a mass similar to the lightest chargino, the second-to-lightest neutralino or the slepton. Lower limits on the masses of charginos in different simplified models range from 193 GeV to 240 GeV for moderate mass splittings, and extend down to mass splittings of 1.5 GeV to 2.4 GeV at the LEP chargino bounds (92.4 GeV). Similar lower limits on degenerate light-flavor sleptons extend up to masses of 251 GeV and down to mass splittings of 550 MeV. Constraints on vector-boson fusion production of electroweak SUSY states are also presented.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Report of the Peer Review Panel for the Rainier Mesa/Shoshone Mountain Flow and Transport Model, Nevada National Security Site, Nevada, Revision 0

The stated purpose of the peer review is to determine whether the model, results, interpretations, and documentation are consistent with the requirements of the FFACO Alternative Modeling Strategy and Decision Process for closure of the CAU. The Peer Review Panel was asked to explicitly answer six questions that provide an overview to its conclusions and recommendations. These six questions and the panel’s responses are presented in the next section of this review report. The Peer Review Panel was asked to present and summarize its primary findings by answering six key peer review questions. These questions are as follows: 1.Is the implementation of the AMS technically and regulatorily sound and defensible? 2.Have conceptual models and/or alternatives for groundwater flow and radionuclide transport in the RM/SM system, permissible based upon available data and observations, been identified to support analyses of radionuclide migration away from the testing areas toward potential receptors? 3.Based on available data that include “maximum transport” and “minimum transport,” do the forecasts represent a range of plausible scenarios and address assumptions for parameters that are most sensitive for saturated zone flow and transport? 4.Can radionuclide transport forecasts be used as an aid in evaluating potential exposure pathways, as intended by the AMS? 5.Has the extent of contamination been bounded such that a regulatory decision can be made regarding CAU closure that is protective of human health and the environment? 6.Has anything that influences the protection of human health and the environment been missed?

54 ENVIRONMENTAL SCIENCES↗

Programmatic Advantages of Linear Equivalent Seismic Models

Underground explosions nonlinearly deform the surrounding earth material and can interact with the free surface to produce spall. However, at typical seismological observation distances the seismic wavefield can be accurately modeled using linear approximations. Although nonlinear algorithms can accurately simulate very near field ground motions, they are computationally expensive and potentially unnecessary for far field wave simulations. Conversely, linearized seismic wave propagation codes are orders of magnitude faster computationally and can accurately simulate the wavefield out to typical observational distances. Thus, devising a means of approximating a nonlinear source in terms of a linear equivalent source would be advantageous both for scenario modeling and for interpretation of seismic source models that are based on linear, far-field approximations. This allows fast linear seismic modeling that still incorporates many features of the nonlinear source mechanics built into the simulation results so that one can have many of the advantages of both types of simulations without the computational cost of the nonlinear computation. In this report we first show the computational advantage of using linear equivalent models, and then discuss how the near-source (within the nonlinear wavefield regime) environment affects linear source equivalents and how well we can fit seismic wavefields derived from nonlinear sources.

58 GEOSCIENCES↗

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↗

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.↗

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↗

Search for supersymmetry in final states with two or three soft leptons and missing transverse momentum in proton-proton collisions at $\sqrt{s}$ = 13 TeV

A search for supersymmetry in events with two or three low-momentum leptons and missing transverse momentum is performed. The search uses proton-proton collisions at $\sqrt{s}$ = 13 TeV collected in the three-year period 2016–2018 by the CMS experiment at the LHC and corresponding to an integrated luminosity of up to 137 fb -1 . The data are found to be in agreement with expectations from standard model processes. The results are interpreted in terms of electroweakino and top squark pair production with a small mass difference between the produced supersymmetric particles and the lightest neutralino. For the electroweakino interpretation, two simplified models are used, a wino-bino model and a higgsino model. Exclusion limits at 95% confidence level are set on $^{\sim0}_{χ2}/^{\sim±}_{χ1}$ masses up to 275 GeV for a mass difference of 10 GeV in the wino-bino case, and up to 205(150) GeV for a mass difference of 7.5 (3) GeV in the higgsino case. The results for the higgsino are further interpreted using a phenomenological minimal supersymmetric standard model, excluding the higgsino mass parameter μ up to 180 GeV with the bino mass parameter M 1 at 800 GeV. In the top squark interpretation, exclusion limits are set at top squark masses up to 540 GeV for four-body top squark decays and up to 480 GeV for chargino-mediated decays with a mass difference of 30 GeV.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Isothermal evaporation of a-pinene ozonolysis SOA: volatility, phase state, and oligomeric composition

We present measurements of the isothermal evaporation of a-pinene ozonolysis secondary organic aerosol (SOA). Using a novel, filter-based method, we replicate previous findings of the volume fraction remaining, presented here as carbon number weighted signal fraction remaining (SFR), as a function of evaporation time. We apply two detailed physical models to interpret the evaporative behavior of both the bulk SOA and individual components. Both models find that a combination of effectively non-volatile products, together with reversible oligomers (or otherwise reactive monomers) having a decomposition lifetime of 9 to 28 hours, best explains the evolution of composition and volatility as particles age in the absence of both organic vapors and oxidants, even under an assumption of relatively viscous (soft wax-like) particles. We find that the residence time in the SOA formation chamber, and time spent undergoing isothermal evaporation, both indicative of the physical age of the aerosol, are the most important experimental parameter determining the evaporation rate. The evolution of volatility observed in these experiments is compared to field measurements in a boreal forest site. The ambient monoterpene-dominated SOA volatility is only reproduced in the laboratory after 24 hours of extended aging in a dilute, dark, oxidant-free environment.

D'Ambro, Emma↗

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↗

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.↗

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↗

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↗

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↗