Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model interpretability”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,585 records · Page 88

Virtual diffraction simulations using the quasi-coarse-grained dynamics method to understand and interpret plasticity contributions during in situ shock experiments

The dynamic deformation response of metallic materials has contributions from dislocations, deformation twinning, and plastic deformation. The current state-of-art techniques can detail the complex mechanistic history of deformation modes under shock loading in real-time using in situ X-ray diffraction (XRD). However, the capability of these experiments to unravel plasticity contributions is challenging due to limitations in interpreting results and the lack of validation from atomistic simulations. Molecular dynamics (MD) simulations can successfully capture various deformation modes in metals and complement experiments using simulated diffractograms at various stages of evolution. However, the difference in length and time scales of MD simulations and experiments is a substantial obstacle in completing and interpreting in situ diffractograms. Therefore, the existing modeling methods require various approximations to model defect evolution and interaction at the mesoscales. However, while using approximations to correlate the peak broadening behavior to the density of dislocations or the shifts/splitting due to the presence of twins, the interpretations of the plasticity contributions from diffractograms are non-trivial, especially when multiple modes of deformation may be operating. This viewpoint discusses combining a mesoscale modeling method called quasi-coarse-grained dynamics and virtual XRD to characterize the plasticity contributions in BCC metals from slip, twinning, and phase transformation behavior. Furthermore, the combined approach shows promise in bridging the mesoscale gap between the capabilities of atomic-scale simulations and in situ experiments to characterize the dynamic deformation of materials.

36 MATERIALS SCIENCE↗

Diagnosing the DARHT Electron Beam-Target Interaction and Hydrodynamic Expansion [Dissertation]

Quantitative electron beam-target interaction studies on the Dual-Axis Radiographic Hydrodynamic Test (DARHT) electron linear induction accelerators have only been performed recently. This includes characterization of the temperature, density, pressure, extent, and expansion velocity of the plasma plume. The results presented in this dissertation present a detailed and unified overview of the accelerator systems, target heating physics, beam transport, and diagnostic tools. Additional information includes calibration sources, radiation hydrodynamics, spectroscopic-quality radiation transport modeling, experimental measurements, and analyses of electron beam driven aluminum experiments. The first set of spatially and temporally resolved spectroscopic measurements of electron beam driven aluminum are presented. Contamination quantification analyses are used to understand the origin of the strong Na-I 3p-3s lines that are observed in absorption within the aluminum plasma continuum. These results inform the creation of the first spectroscopic-quality radiation transport model that links several atomic physics codes to interpret the conditions from which the Na-I lines originate. A good agreement is found between the surface analysis results and the model which confirms the concentration of the sodium present within the aluminum alloy foil material. It also demonstrates, for the first time on electron beam-driven target experiments, the ability to interpret plasma conditions from measured absorption lines. In a second experimental campaign which focuses on pure aluminum, the Al-I 3p-4s and 3p3d doublets are both measured in emission. A detailed analysis of the Al-I 3p-4s doublet reveals that the lines undergo moderate self-absorption. A simple model of the self-absorption effect is successfully used to match measured spectra at various temperature/density/plasma scale length combinations. These measurements led to the realization of the minimum density that can be resolved by the spectrometers for the Al-I 3p-4s lines due to the large slit width required to observe a signal on aluminum. These measurements also demonstrate the sensitivity of visible and long wave UV spectroscopy to minor changes in both temperature and density. The simple self-absorption model will be useful for analysis of other beam-target interaction experiments with spectra exhibiting either self-absorption or full self-reversal. Substantial headway has been made on the modeling front by linking together several codes needed to model both the energy deposition, hydrodynamic motion, and atomic kinetics to produce synthetic spectral calculations that are compared with experimental measurements. There exists ample space for improvement, especially with benchmarking the hydrodynamics codes and equation-of-state tables with experimental measurements. The X-ray diagnostics required to make these new measurements along with the simulation capabilities required to interpret the results are under development and will be the subject of future studies.

43 PARTICLE ACCELERATORS↗

The impact of clouds on the shortwave radiation budget of the surface-atmosphere system - Interfacing measurements and models

Two datasets have been combined to demonstrate how the availability of more comprehensive datasets could serve to elucidate the shortwave radiative impact of clouds on both the atmospheric column and the surface. These datasets consist of two measurements of net downward shortwave radiation: one of near-surface measurements made at the Boulder Atmospheric Observatory tower, and the other of collocated top-of-the-atmosphere measurements from the Earth Radiation Budget Experiment. Output from the European Centre for Medium Range Weather Forecasts General Circulation Model also has been used as an aid in interpreting the data, while the data have in turn been employed to validate the model's shortwave radiation code as it pertains to cloud radiation properties. Combined, the datasets and model demonstrate a strategy for determining under what conditions the shortwave radiative impact of clouds leads to a heating or cooling of the atmospheric column. The datasets also show, in terms of a linear slope-offset algorithm for retrieving the net downward shortwave radiation at the surface from satellite measurements, that the clouds present during this study produced a modest negative bias in the retrieved surface flux relative to that inferred from a clear-sky algorithm.

Cess, Robert D.↗

The Impact of Clouds on the Shortwave Radiation Budget of the Surface-Atmosphere System: Interfacing Measurements and Models

Two datasets have been combined to demonstrate how the availability of more comprehensive datasets could serve to elucidate the shortwave radiative impact of clouds on both the atmospheric column and the surface. These datasets consist of two measurements of net downward shortwave radiation: one of near-surface measurements made at the Boulder Atmospheric Observatory tower, and the other of collocated top-of-the-atmosphere measurements from the Earth Radiation Budget Experiment. Output from the European Centre for Medium-Range Weather Forecasts General Circulation Model also has been used as an aid in interpreting the data, while the data have in turn been employed to validate the model's shortwave radiation code as it pertains to cloud radiation properties. Combined, the datasets and model demonstrate a strategy for determining under what conditions the shortwave radiative impact of clouds leads to a heating or cooling of the atmospheric column. The datasets also show, in terms of a linear slope-offset algorithm for retrieving the net downward shortwave radiation at the surface from satellite measurements, that the clouds present during this study produced a modest negative bias in the retrieved surface flux relative to that inferred from a clear-sky algorithm.

Cess, Robert D.↗

Solar Wind Implantation Into the Lunar Regolith: Monte Carlo Simulations of H Retention in a Surface With Defects and the H2 Exosphere

The solar wind implants protons into the top 20–30 nm of lunar regolith grains, and the implanted hydrogen will diffuse out of the regolith but also interact with oxygen in the regolith oxides. We apply a statistical approach to estimate the diffusion of hydrogen in the regolith hindered by forming temporary bonds with regolith oxygen atoms. A Monte Carlo simulation was used to track the temporal evolution of bound OH surface content and the H2 exosphere. The model results are consistent with the interpretation of the Chandrayaan-1 M3 observations of infrared absorption spectra by surface hydroxyls as discussed in Li and Milliken (2017). Themodel reproduced the latitudinal concentration of OH by using a Gaussian energy distribution of f(U(o) = 0.5eV, U(w) = 0.078 – 0.1 eV) to characterize the activation energy barrier to the diffusion of hydrogen in space weathered regolith. In addition, the model results of the exospheric content of H2 are consistent with observations by the Lyman Alpha Mapping Project on the Lunar Reconnaissance Orbiter. Therefore, we provide support for hydroxyl formation by chemically trapped solar wind protons.

O. J. Tucker↗

Search for the decay of the Higgs boson to a $Z$ boson and a light pseudoscalar particle decaying to two photons

A search for the decay of the Higgs boson to a Z boson and a light, pseudoscalar particle, a, decaying respectively to two leptons and to two photons is reported. The search uses the full LHC Run 2 proton–proton collision data at $\sqrt{s}$ = 13 TeV, corresponding to 139 fb –1 collected by the ATLAS detector. This is one of the first searches for this specific decay mode of the Higgs boson, and it probes unexplored parameter space in models with axion-like particles (ALPs) and extended scalar sectors. The mass of the $a$ particle is assumed to be in the range 0.1–33 GeV. The data are analysed in two categories: a merged category where the photons from the $a$ decay are reconstructed in the ATLAS calorimeter as a single cluster, and a resolved category in which two separate photons are detected. The main background processes are from Standard Model Z boson production in association with photons or jets. The data are in agreement with the background predictions, and upper limits on the branching ratio of the Higgs boson decay to $Za$ times the branching ratio $a → γγ$ are derived at the 95% confidence level and they range from 0.08% to 2% depending on the mass of the $a$ particle. The results are also interpreted in the context of ALP models.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Evidence for the charge asymmetry in $pp$ → $t\bar{t}$ production at $\sqrt{s}$ = 13 TeV with the ATLAS detector

Inclusive and differential measurements of the top–antitop ($t\bar{t}$) charge asymmetry ${A}_{C}^{t\bar{t}}$ and the leptonic asymmetry ${A}_{C}^{\ell\bar{\ell}}$ are presented in proton–proton collisions at $\sqrt{s}$ = 13 TeV recorded by the ATLAS experiment at the CERN Large Hadron Collider. The measurement uses the complete Run 2 dataset, corresponding to an integrated luminosity of 139 fb –1 , combines data in the single-lepton and dilepton channels, and employs reconstruction techniques adapted to both the resolved and boosted topologies. A Bayesian unfolding procedure is performed to correct for detector resolution and acceptance effects. The combined inclusive $t\bar{t}$ charge asymmetry is measured to be ${A}_{C}^{t\bar{t}}$ = 0.0068 ± 0.0015, which differs from zero by 4.7 standard deviations. Differential measurements are performed as a function of the invariant mass, transverse momentum and longitudinal boost of the $t\bar{t}$ system. Both the inclusive and differential measurements are found to be compatible with the Standard Model predictions, at next-to-next-to-leading order in quantum chromodynamics perturbation theory with next-to-leading-order electroweak corrections. The measurements are interpreted in the framework of the Standard Model effective field theory, placing competitive bounds on several Wilson coefficients.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Topological heavy fermions in magnetic field

The recently introduced topological heavy fermion model (THFM) provides a means for interpreting the low-energy electronic degrees of freedom of the magic angle twisted bilayer graphene as hybridization amidst highly dispersing topological conduction and weakly dispersing localized heavy fermions. In order to understand the Landau quantization of the ensuing electronic spectrum, a generalization of THFM to include the magnetic field B is desired, but currently missing. Here we provide a systematic derivation of the THFM in B and solve the resulting model to obtain the interacting Hofstadter spectra for single particle charged excitations. While naive minimal substitution within THFM fails to correctly account for the total number of magnetic subbands within the narrow band i.e., its total Chern number, our method—based on projecting the light and heavy fermions onto the irreducible representations of the magnetic translation group— reproduces the correct total Chern number. Analytical results presented here offer an intuitive understanding of the nature of the (strongly interacting) Hofstadter bands.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Adaptive Recovery Model: Designing Systems for Testing Tracing and Vaccination to Support COVID-19 Recovery Planning.

This report documents a new approach to designing disease control policies that allocate scarce testing, contact tracing, and vaccination resources to better control community transmission of COVID19 or similar diseases. The Adaptive Recovery Model (ARM) combines a deterministic compartmental disease model with a stochastic network disease propagation model to enable us to simulate COVID-19 community spread through the lens of two complementary modeling motifs. ARM contact networks are derived from cell-phone location data that have been anonymized and interpreted as individual arrivals to specic public locations. Modeling disease spread over these networks allows us to identify locations within communities conducive to rapid disease spread. ARM applies this model- and data-derived abstractions of community transmission to evaluate the effectiveness of disease control measures including targeted social distancing, contact tracing, testing and vaccination. The architecture of ARM provides a unique capacity to help decision makers understand how best to deploy scarce testing, tracing and vaccination resources to minimize disease-spread potential in a community. This document details the novel mathematical formulations underlying ARM, presents a dynamical stability analysis of the deterministic model components, a sensitivity analysis of control parameters and network structure, and summarizes a process for deriving contact networks from cell-phone location data. An example use case steps through applying ARM to evaluate three targeted social distancing policies using Bernalillo County, New Mexico as an exemplar test locale. This step-by-step analysis demonstrates how ARM can be used to measure the relative performance of competing public health policies. Initial scenario tests of ARM shows that ARMs design focus on resource utilization rather than simple incidence prediction can provide decision makers with additional quantitative guidance for managing ongoing public health emergencies and planning future responses.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Model-based estimation and prediction of task-imposed mental workload

Mental workload has been an area of intensive research for better than a decade. One specific area of interest in aircrew related workload research is concerned with the development of quantitative indices of workload in aircraft piloting tasks. This paper presents a model-based approach for quantifying mental workload in operational terms. The suggested modeling framework is based on an interpreted Petri net characterization of a task in which 'places' are equated to specific task-related activities and 'transitions' are viewed as internal or external forcing events. It is shown that within this framework quantitative assessments can be made of both cumulative and instantaneous workload associated with the performance of a task and its individual component subtasks. It is suggested that insights gained from analyzing task-specific workload within this modeling paradigm can suggest plausible explanations for reconciling discrepancies between subjectively elicited workload estimates and behavioral/performance measures.

Madni, A. M.↗

Investigating Galactic Structure with COBE/DIRBE

In this work I applied the current version of the SKY model of the point source sky to the interpretation of the diffuse all-sky emission observed by COBE/DIRBE. The goal was to refine the SKY model using the all-sky DIRBE maps of the Galaxy, in order that a search could be made for an isotropic cosmic back-ground. Arendt et al. constructed their "Faint Source Model" (FSM) to remove Galactic foreground stars from the ZSMA products. The FSM mimics SKY version 1 but Arendt et al.'s concluded that it was inadequate to seek cosmic background emission because of the sizeable residual emission in the ZSMA products after this starlight subtraction. At this point I can only support Arendt et al.'s conclusion, namely that such models are currently inadequate to reveal a cosmic background. Even SKY5 yields the same disappointing result. On the positive side, our knowledge of the fundamental parameters of stars is increasing through the legacy of the Hipparcos mission. These parameters are critical inputs to any model of the stellar sky and it is my hope to revisit the problem once the new values are for visual-absolute magnitudes and stellar, space densities. I also hope to make progress through an independent analysis of Pioneer optical starlight with SKY, under a NASA Heliospheric Mission Guest Investigator proposal.

Cohen, Martin↗

Modeling the spatiotemporal organization of velocity storage in the vestibuloocular reflex by optokinetic studies

1. A generalized three-dimensional state space model of visual vestibular interaction was developed. Matrix and dynamical system operators associated with inputs from the semicircular canals, otolith velocity estimator, and the visual system have been incorporated into the model, which focus on their relationship to the velocity storage integrator. 2. A relationship was postulated between the eigenvalues and the direction of the eigenvectors of the system matrix and the orientation of the spatial vertical. It was assumed that the system matrix for a tilted position was a composition of two linear transformations of the system matrix for the upright position. One transformation modifies the eigenvalues of the system matrix, whereas another rotates the eigenvectors of the system matrix. The pitch and roll eigenvectors rotate with the head, whereas the yaw axis eigenvector remains approximately spatially invariant. 3. Based on the three-dimensional model, a computational procedure was formulated to identify the eigenvalues and eigenvectors of the system matrix with the use of a modification of the marquardt algorithm. With the use of data obtained from a monkey, it was shown that the three-dimensional behavior of velocity storage cannot be predicted solely in terms of its time constants, i.e., the inverse of its eigenvalues. With the use of the same eigenvalues the data could either be fit or not fit, depending on the eigenvector directions. Therefore, it is necessary to specify eigenvector directions when characterizing velocity storage in three dimensions. 4. Parameters found with the use of the Marquardt algorithm were incorporated into the model. Diagonal matrices in a head coordinate frame were introduced for coupling the visual system to the integrator and to the direct optokinetic pathway. Simulations of optokinetic nystagmus (OKN) and optokinetic after-nystagmus (OKAN) were run. The model predicted the behavior of yaw and pitch OKN and OKAN when the animal is upright. It also predicted the cross-coupling in the side down position. The trajectories in velocity space were also accurately simulated. 5. One of the predictions of the model is that when the stimulus direction is along an eigenvector, the trajectory in velocity space is a straight line. Using the "spectral width" of the residuals from a straight line sequence during OKAN, we developed a methodology to estimate how close the OKAN decay was to an eigenvector trajectory. 6. Thus we have developed a model-based approach for studying and interpreting the response characteristics of velocity storage in three dimensions.(ABSTRACT TRUNCATED AT 400 WORDS).

Non-NASA Center↗

The entropy of galaxy spectra: how much information is encoded?

Abstract The inverse problem of extracting the stellar population content of galaxy spectra is analysed here from a basic standpoint based on information theory. By interpreting spectra as probability distribution functions, we find that galaxy spectra have high entropy, thus leading to a rather low effective information content. The highest variation in entropy is unsurprisingly found in regions that have been well studied for decades with the conventional approach. We target a set of six spectral regions that show the highest variation in entropy – the 4000 Å break being the most informative one. As a test case with real data, we measure the entropy of a set of high-quality spectra from the Sloan Digital Sky Survey, and contrast entropy-based results with the traditional method based on line strengths. The data are classified into star-forming (SF), quiescent (Q), and active galactic nucleus (AGN) galaxies, and show – independently of any physical model – that AGN spectra can be interpreted as a transition between SF and Q galaxies, with SF galaxies featuring a more diverse variation in entropy. The high level of entanglement complicates the determination of population parameters in a robust, unbiased way, and affects traditional methods that compare models with observations, as well as machine learning (especially deep learning) algorithms that rely on the statistical properties of the data to assess the variations among spectra. Entropy provides a new avenue to improve population synthesis models so that they give a more faithful representation of real galaxy spectra.

Ferreras, Ignacio (ORCID:0000000345843127)↗

Neutrino Interaction Model and Uncertainties for MicroBooNE Analyses

The MicroBooNE experiment is searching for an excess of electron-like events at low neutrino energy, as well as making the first high-statistics measurements of neutrino interactions on argon. Interpretation of these results requires a reliable model of neutrino interaction physics with well-motivated uncertainties. This note describes the neutrino interaction model and uncertainties adopted for use in MicroBooNE’s cross section measurements and low-energy excess search.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Spin-coupling topology in the copper hexamer compounds A 2 Cu 3 O(SO 4 ) 3 ( A =Na, K)

The compounds A 2 Cu 3 O(SO 4 ) 3 (A=Na,K) are characterized by copper hexamers that are weakly coupled to realize antiferromagnetic order below TN≈3K. They constitute quantum spin systems with S=1 triplet ground states. In this work, we investigated the energy-level splittings of the copper hexamers by inelastic neutron scattering experiments covering the entire range of the magnetic excitation spectra. Additionally, the observed transitions are governed by very unusual selection rules that we ascribe to the underlying spin-coupling topology. This is rationalized by model calculations that allow an unambiguous interpretation of the magnetic excitations concerning both the peak assignments and the nature of the spin-coupling parameters.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Neural net modeling of equilibria in NSTX-U

Neural networks (NNs) offer a path towards synthesizing and interpreting data on faster timescales than traditional physics-informed computational models. In this work we develop two NNs relevant to equilibrium and shape control modeling, which are part of a suite of tools being developed for the National Spherical Torus Experiment-Upgrade for fast prediction, optimization, and visualization of plasma scenarios. The networks include Eqnet, a free-boundary equilibrium solver trained on the EFIT01 (Equilibrium FITtting 01) reconstruction algorithm, and Pertnet, which is trained on the Gspert code and predicts the non-rigid plasma response, a nonlinear term that arises in shape control modeling. The NNs are trained with different combinations of inputs and outputs in order to offer flexibility in use cases. In particular, Eqnet can use magnetic diagnostics as inputs and act as an EFIT-like reconstruction algorithm, or, by using pressure and current profile information the NN can act as a forward Grad–Shafranov equilibrium solver. This forward-mode version is envisioned to be implemented in the suite of tools for simulation of plasma scenarios. The reconstruction-mode version gives some performance improvements compared to the online reconstruction code real-time EFIT, especially when vessel eddy currents are significant. Here, we report strong performance for all NNs indicating that the models could reliably be used within closed-loop simulations or other applications. Some limitations are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Modeling the Role of Dislocation Substructure During Class M and Exponential Creep

The different substructures that form in the power-law and exponential creep regimes for single phase crystalline materials under various conditions of stress, temperature and strain are reviewed. The microstructure is correlated both qualitatively and quantitatively with power-law and exponential creep as well as with steady state and non-steady state deformation behavior. These observations suggest that creep is influenced by a complex interaction between several elements of the microstructure, such as dislocations, cells and subgrains. The stability of the creep substructure is examined in both of these creep regimes during stress and temperature change experiments. These observations are rationalized on the basis of a phenomenological model, where normal primary creep is interpreted as a series of constant structure exponential creep rate-stress relationships. The implications of this viewpoint on the magnitude of the stress exponent and steady state behavior are discussed. A theory is developed to predict the macroscopic creep behavior of a single phase material using quantitative microstructural data. In this technique the thermally activated deformation mechanisms proposed by dislocation physics are interlinked with a previously developed multiphase, three-dimensional. dislocation substructure creep model. This procedure leads to several coupled differential equations interrelating macroscopic creep plasticity with microstructural evolution.

Raj, S. V.↗

Resolving glacial isostatic adjustment (GIA) in response to modern and future ice loss at marine grounding lines in West Antarctica

Accurate glacial isostatic adjustment (GIA) modelling in the cryosphere is required for interpreting satellite, geophysical and geological records and for assessing the feedbacks of Earth deformation and sea-level change on marine ice-sheet grounding lines. GIA modelling in areas of active ice loss in West Antarctica is particularly challenging because the ice is underlain by laterally varying mantle viscosities that are up to several orders of magnitude lower than the global average, leading to a faster and more localised response of the solid Earth to ongoing and future ice-sheet retreat and necessitating GIA models that incorporate 3-D viscoelastic Earth structure. Improvements to GIA models allow for computation of the viscoelastic response of the Earth to surface ice loading at sub-kilometre resolution, and ice-sheet models and observational products now provide the inputs to GIA models at comparably unprecedented detail. However, the resolution required to accurately capture GIA in models remains poorly understood, and high-resolution calculations come at heavy computational expense. We adopt a 3-D GIA model with a range of Earth structure models based on recent seismic tomography and geodetic data to perform a comprehensive analysis of the influence of grid resolution on predictions of GIA in the Amundsen Sea Embayment (ASE) in West Antarctica. Through idealised sensitivity testing down to sub-kilometre resolution with spatially isolated ice loading changes, we find that a grid resolution of ~$\frac{1}{3}$ of the radius of the load or higher is required to accurately capture the elastic response of the Earth. However, when we consider more realistic, spatially coherent ice loss scenarios based on modern observational records and future ice-sheet model projections and adopt a viscoelastic Earth, we find that predicted deformation and sea-level change along the grounding line converge to within 5 % with grid resolutions of 7.5 km or higher, and to within 2% for grid resolutions of 3.75 km and higher, even when the input ice model is on a 1 km grid. Furthermore, we show that low mantle viscosities beneath the ASE lead to viscous deformation that contributes to the instrumental record on decadal timescales and equals or dominates over elastic effects by the end of the 21st century. Our findings suggest that for the range of resolutions of 1.9–15 km that we considered, the error due to adopting a coarser grid in this region is negligible compared to the effect of neglecting viscous effects and the uncertainty in the adopted mantle viscosity structure.

54 ENVIRONMENTAL SCIENCES↗