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,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Comment on 'The Global Resurfacing of Venus' by R. G. Strom, G.G. Schaber, and D.D. Dawson

The distribution of impact craters on Venus has been the subject of a great deal of analysis since the return of Magellan data. Phillips el al. (1992) performed Monte Carlo two-dimensional (2-D) modeling of the areal distribution of craters, and the results of that exercise allowed a restricted, but still quite large, range of possible planetary resurfacing histories, including the possibility that the crater, were emplaced on a geologically inactive planet. However, the nonrandom distribution of embayed and deformed craters (Phillips el al., 1992), the hypsometric distribution of craters (Herrick and Phillips, 1994), the varied degradation states of craters (Izenberg et al., 1994), their nonrandom distribution with different geologic terrain types (Namiki and Solomon, 1994; Price et al, 1994), and three-dimensional resurfacing modeling (Bullock el al., 1993) all seem to argue against that particular possibility. In contrast, Strom el al. (1994) have collected a refined and more comprehensive data set of impact features, and they input these data into more sophisticated 2-D Monte Carlo modeling and statistical analyses of the areal distribution of craters, the hypsometric distribution of craters, and the number of embayed craters. They concluded that 'Venus experienced a global resurfacing event about 300 m.y. ago followed by a dramatic reduction of volcanism and tectonism. This global resurfacing event ended abruptly (less than 10 m.y.). The present crater population has accumulated since then and remains largely intact . . . only about 4%-6% of the planet has been volcanically resurfaced since the global event . . .' If these conclusions are well-founded, this work certainly represents a significant advancement in restricting tile number of plausible resurfacing histories for the planet. If Strom et al. (1994) are correct, it would also mean that all of the other aforementioned works are in error to various degrees, or at least represent overzealous interpretation of the data. However, we have identified apparent flaws in the observations, modeling, and interpretations presented by Strom el al. (1994) that lead us to question whether their conclusions are warranted. We limit our comments to three areas of their analysis: (1) observations pertaining to the number and area of disrupted and pristine craters and crater-related features, (2) modeling of the areal and elevation distribution of craters, and (3) interpretations of resurfacing models.

Herrick, Robert R.↗
DOE OSTI2025

eDNAjoint: An R package for interpreting paired or semi‐paired environmental DNA and traditional survey data in a Bayesian framework

Abstract Environmental DNA (eDNA) sampling is increasingly used in surveys of species distribution as a potentially sensitive and efficient monitoring method. Yet access to modelling tools designed specifically for interpreting this new data type lags behind its ubiquity. While occupancy modelling software has dominated the analytical landscape for eDNA data analysis of single species, this type of model may not always be the most appropriate. The rate of eDNA detection often corresponds to species density, rather than just occupancy, and researchers often have access to observations from non‐genetic sampling methods at the same sites. To provide users access to a modelling framework designed to maximize the use of all available data, we developed an R package, eDNAjoint . The package provides an easy‐to‐use interface for fitting a ‘joint’ model that integrates data from paired or semi‐paired eDNA and traditional surveys in a Bayesian framework. The model can be used to estimate parameters like the probability of a false positive eDNA detection and mean catch rate at a site, and the package allows access to multiple model variations and Bayesian prior customization. Additional functionality can be used for model selection, summarising posteriors and comparing the relative sensitivities of the two survey methods. We demonstrate the use of eDNAjoint by fitting a variation of the model with site‐level covariates that scale the sensitivity of eDNA sampling relative to traditional sampling. The example workflow uses binary eDNA and seine count data for the endangered tidewater goby ( Eucyclogobius newberryi ) from a study by Schmelzle and Kinziger (2016). This use case includes a prior sensitivity analysis and an evaluation of the relationship between detection rates and environmental variables. eDNAjoint has the potential to greatly increase the range of users who will be able to rigorously analyse eDNA and traditional survey data in a Bayesian framework, understand if and how eDNA can improve monitoring practices, and gain confidence in the interpretability of eDNA data.

Keller, Abigail G. [Department of Environment Scie↗
NASA NTRS1994

Tropical intraseasonal oscillation, super cloud clusters, and cumulus convection schemes

A new framework for interpreting the origin of the tropical intraseasonal oscillation (TISO), which avoids the speed and scale selection problems in the previous theories, is proposed in this study. In this interpretation TISO is viewed as an oscillation driven by an eastward moving convective region. This convective region consists of one or more super cloud clusters originating in the Indian Ocean and terminating in mid-Pacific, and is then followed by another convective region arising in the Indian Ocean in a period of 40-50 days. Additionally, a formal analogy is pointed out between super cloud clusters and the middle-latitude baroclinic wave packets. This study includes a simulation of TISO in a 2D model to support our interpretation. Experiments were conducted with four different convection schemes. The authors advocate that the successful simulation of TISO depends on the successful simulation of super clouds cluster, which in turn depends on the successful simulation of the life cycle of cloud clusters, which further in turn depends on the choice of cumulus convection scheme. What makes a cumulus convection scheme successful in simulating TISO is discussed.

Chao, Winston C.↗
DOE OSTI2024

Investigating the effects of precise mass measurements of Ru and Pd isotopes on machine learning mass modeling

Atomic masses are a foundational quantity in our understanding of nuclear structure, astrophysics, and fundamental symmetries. The longstanding goal of creating a predictive global model for the binding energy of a nucleus remains a significant challenge, however, and prompts the need for precise measurements of atomic masses to serve as anchor points for model developments. We present precise mass measurements of neutron-rich Ru and Pd isotopes performed at the Californium Rare Isotope Breeder Upgrade facility at Argonne National Laboratory using the Canadian Penning Trap mass spectrometer. The masses of 108 Ru, 110 Ru, and 116 Pd were measured to a relative mass precision $\delta$⁢$m/m$ ≈ 10 -8 via the phase-imaging ion-cyclotron-resonance technique, and represent an improvement of approximately an order of magnitude over previous measurements. Further, these mass data were used in conjunction with the physically interpretable machine learning (PIML) model, which uses a mixture density neural network to model mass excesses via a mixture of Gaussian distributions. The effects of our new mass data on a Bayesian-updating of a PIML model are presented.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Multi-jet gas cooling of in-beam foils or specimens: CFD predictions of the convective heat-transfer coefficient

An experiment was designed to investigate the possible enhancement of the convective heat-transfer coefficient by utilizing multiple, parallel jets in the cooling of a small heated surface, such as typically induced by an accelerated ion beam on a thin foil or specimen. The hot spot was provided using a small electrically heated plate. It was found that heat-transfer calculations by means of simple empirical methods based on dimensional analysis are not useful in this case and that advanced computational fluid dynamics (CFD) modelling is essential to interpret the results. It is shown that enhanced convective cooling can indeed be obtained with a multi-jet configuration as compared to a single-jet configuration but only under very selective conditions. An improperly designed multi-jet configuration can also provide significantly reduced cooling relative to the single-jet case and the estimation of the behavior of any particular jet geometry is eminently non-intuitive. CFD provides acceptable quantitative results and seem to be the only tool available to gain an understanding of these complex flows where simple models and “rules of thumb” cannot be relied upon.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS1978

Transport of contaminants in the planetary boundary layer

A planetary boundary layer model is described and used to simulate PBL phenomena including cloud formation and pollution transport in the San Francisco Bay Area. The effect of events in the PBL on air pollution is considered, and governing equations for the average momentum, potential temperature, water vapor mixing ratio, and air contaminants are presented. These equations are derived by integrating the basic equations vertically through the mixed layer. Characteristics of the day selected for simulation are reported, and the results suggest that the diurnally cyclic features of the mesoscale motion, including clouds and air pollution, can be simulated in a readily interpretable way with the model.

Lee, I. Y.↗
NASA NTRS1987

Theoretical studies of Resonance Enhanced Stimulated Raman Scattering (RESRS) of frequency doubled Alexandrite laser wavelength in cesium vapor

The third phase of research will focus on the propagation and energy extraction of the pump and SERS beams in a variety of configurations including oscillator structures. In order to address these questions a numerical code capable of allowing for saturation and full transverse beam evolution is required. The method proposed is based on a discretized propagation energy extraction model which uses a Kirchoff integral propagator coupled to the three level Raman model already developed. The model will have the resolution required by diffraction limits and will use the previous density matrix results in the adiabatic following limit. Owing to its large computational requirements, such a code must be implemented on a vector array processor. One code on the Cyber is being tested by using previously understood two-level laser models as guidelines for interpreting the results. Two tests were implemented: the evolution of modes in a passive resonator and the evolution of a stable state of the adiabatically eliminated laser equations. These results show mode shapes and diffraction losses for the first case and relaxation oscillations for the second one. Finally, in order to clarify the computing methodology used to exploit the speed of the Cyber's computational speed, the time it takes to perform both of the computations previously mentioned to run on the Cyber and VAX 730 must be measured. Also included is a short description of the current laser model (CAVITY.FOR) and a flow chart of the test computations.

Lawandy, Nabil M.↗
DOE OSTI2021

Machine-Learned Manifold-Based Models for Large Eddy Simulation of Turbulent Combustion

Reduced-order manifold combustion models are commonly used to lower the cost of reacting Large Eddy Simulations (LES) and may be derived either from physical principles as in flamelet models or from data-driven methods like Principal Component Analysis (PCA). In either case, neural networks are increasingly used as part of these models to provide a nonlinear mapping between a small set of pre-defined variables that parameterize the manifold and outputs of interest, such as reaction rates. In this work, we propose a new manifold-based modeling approach that combines the definition of the manifold-parameterizing variables (linear combinations of species), the nonlinear mapping to the outputs, and closure of filtered quantities for LES into the structure of a single neural network. This allows the process used to train the neural network to simultaneously optimize both the functional form of the model and the identities of the inputs to the model. The new approach can flexibly incorporate thermochemical data from any combustion system; if trained on data from 1D flames it can be interpreted as an optimized flamelet model, but it can also be used to learn models from data from more complex configurations. This work presents a priori evaluations of the new approach in both contexts. Evaluation using data from 1D premixed flames demonstrates the physical interpretability of the manifold variables generated by the new approach. Evaluation using data from direct numerical simulations of turbulent flames shows improved predictions relative to either flamelet or PCA-based models in a more complex configuration.

47 OTHER INSTRUMENTATION↗
NASA NTRS1994

Assessment of impact damage of composite rocket motor cases

This contract reviewed the available literature on mechanisms of low velocity impact damage in filament wound rocket motor cases, MDE methods to quantify damage, critical coupon level test methods, manufacturing and material process variables and empirical and analytical modeling off impact damage. The critical design properties for rocket motor cases are biaxial hoop and axial tensile strength. Low velocity impact damage is insidious because it can create serious nonvisible damage at very low impact velocities. In thick rocket motor cases the prevalent low velocity impact damage is fiber fracture and matrix cracking adjacent to the front face. In contrast, low velocity loading of thin wall cylinders induces flexure, depending on span length and the flexure induces delamination and tensile cracking on the back face wall opposed to impact occurs due to flexural stresses imposed by impact loading. Important NDE methods for rocket motor cases are non-contacting methods that allow inspection from one side. Among these are vibrothermography, and pulse-echo methods based on acoustic-ultrasonic methods. High resolution techniques such as x-ray computed tomography appear to have merit for accurate geometrical characterization of local damage to support development of analytical models of micromechanics. The challenge of coupon level testing is to reproduce the biaxial stress state that the full scale article experiences, and to determine how to scale the composite structure to model full sized behavior. Biaxial tensile testing has been performed by uniaxially tensile loading internally pressurized cylinders. This is experimentally difficult due to gripping problems and pressure containment. Much prior work focused on uniaxial tensile testing of model filament wound cylinders. Interpretation of the results of some studies is complicated by the fact that the fabrication process did not duplicate full scale manufacturing. It is difficult to scale results from testing subscale cylinders since there are significant differences in out time of the resins relative to full scale cylinder fabrication, differences in hoop fiber tensioning and unsatisfactory coupon configurations. It appears that development of a new test method for subscale cylinders is merited. Damage tolerance may be improved by material optimization that uses fiber treatments and matrix modifications to control the fiber matrix interface bonding. It is difficult to develop process optimization in subscale cylinders without also modeling the longer out times resins experience in full scale testing. A major breakthrough in characterizing the effect of impact damage on residual strength, and understanding how to scale results of subscale evaluations, will be a sound micromechanical model that described progressive failure of the composite. Such models will utilize a three dimensional stress analysis due to the complex nature of low velocity impact stresses in thick composites. When these models are coupled with non-contact NDE methods that geometrically characterize the damage and acoustic methods that characterize the effective local elastic properties, accurate assessment of residual strength from impact damage may be possible. Directions for further development are suggested.

Paris, Henry G.↗
NASA NTRS1971

Rocket-based measurement of Birkeland currents related to an auroral arc and electrojet.

A rocket-borne experiment performed to study currents associated with a quiet auroral arc is discussed. The magnetic field in the vicinity of the arc was measured with a vector magnetometer, while the orientation of the payload was determined with a lunar-aspect sensor. Possible current configurations were inferred by constructing model current systems that reproduced the magnetic field variations measured along the flight path. The data are interpreted in terms of a model current system consisting of a northwestward electrojet and two oppositely directed Birkeland sheet currents, all lying in planes approximately parallel to the auroral arc. The current density integrated through the 16-km thickness of each Birkeland current sheet was found to be roughly 0.16 A/m. The electrojet current was roughly 6000 A.

Park, R. J.↗
NASA NTRS1976

Measurements of field-aligned currents in a multiple auroral arc system

A rocket-borne experiment to study the currents associated with a system of multiple auroral arcs was conducted at Poker Flat, Alaska, at 1122 UT on Feb. 2, 1972. The magnetic field in the vicinity of the auroral system was measured with a cesium vector magnetometer. Possible configurations were inferred by constructing model current systems that reproduced the magnetic field variations measured along the flight path. The data are interpreted in terms of a model current system consisting of two eastward electrojets and one westward electrojet and three pairs of oppositely directed Birkeland sheet currents, all lying in a plane approximately parallel to the auroral arcs. Sheet thicknesses ranged from 20 to 60 km and current densities from 10 to 45 microampere/sq m; the electrojet currents ranged from 1000 to 2000 A. A possible alternate model consisted of four pairs of sheets whose thicknesses range from 10 to 40 km with current densities from 10 to 90 microampere/sq m. There was quite good agreement between the locations of the visual arcs and the upward current sheets. The overall current configuration is discussed in view of the theoretical models constructed by Atkinson and Sato and Holzer and of other observations.

Sesiano, J.↗
NASA NTRS1976

A comparison of one-dimensional theoretical models of stratospheric minor constituents

The density profiles of stratospheric minor constituents are compiled from several one-dimensional theoretical models that have appeared in the literature. The models are compared with each other and also with observations if data are available. The similarity and disparity of the models are clarified and discussed. If a large difference exists among the models, attempts are made to interpret it as much as possible in terms of the physical and chemical data employed in each model. The inadequacy of one-dimensional models in calculating the realistic ozone density distributions and the problems related to it are briefly discussed.

Shimazaki, T.↗
NASA NTRS1990

Proceedings of the MEVTV Workshop on The Evolution of Magma Bodies on Mars

The workshop focused on many of the diverse approaches related to the evolution of magma bodies on Mars that have been pursued during the course of the Mars Evolution of Volcanism, Tectonism, and Volatiles (MEVTV) Program. Approximately 35 scientists from the Mars volcanology, petrology, geochemistry, and modeling communities attended. Segments of the meeting concentrated of laboratory analyses and investigations of SNC meteorites, the interpretation of Viking Orbiter and Lander datasets, and the interpretation of computer codes that model volcanic and tectonic processes on Mars. Abstracts of these reports are presented.

Mouginis-Mark, P.↗
NASA NTRS1994

The Sun as a Variable Star: Solar and Stellar Irradiance Variations; Colloquium of the International Astronomical Union, 143rd, Boulder, CO, Jun. 20-25, 1993

Variations in solar and stellar irradiances have long been of interest. An International Astronomical Union (IAU) colloquium reviewed such relevant subjects as observations, theoretical interpretations, and empirical and physical models, with a special emphasis on climatic impact of solar irradiance variability. Specific topics discussed included: (1) General Reviews on Observations of Solar and Stellar Irradiance Variability; (2) Observational Programs for Solar and Stellar Irradiance Variability; (3) Variability of Solar and Stellar Irradiance Related to the Network, Active Regions (Sunspots and Plages), and Large-Scale Magnetic Structures; (4) Empirical Models of Solar Total and Spectral Irradiance Variability; (5) Solar and Stellar Oscillations, Irradiance Variations and their Interpretations; and (6) The Response of the Earth's Atmosphere to Solar Irradiance Variations and Sun-Climate Connections.

Pap, Judit M.↗
NASA NTRS2010

Simulation and 'TWINS Observations of the 22 July 2009 Storm

TWINS is the first mission to perform stereo imaging of the Earth's ring current. The magnetic storm on 22 July 2009 is the largest storm observed since TWINS began routine stereo imaging in June 2008. On 22 July 2009, the Dst dropped to nearly -80nT at 7:00 and 10:00 UT. During the main phase and at the peak of the storm, TWINS 1 and 2 were near apogee and moving from pre-dawn to post-dawn local time. The energetic neutral atom (ENA) imagers on the 2 spacecraft captured the storm intensification and the formation of the partial ring current. The peak of the ENA emissions was seen in the midnight-to-dawn local-time sector. The development of this storm has been simulated using the Comprehensive Ring Current Model (CRCM) to understand and interpret the observed signatures. We perform CRCM runs with constant and time-varying magnetic field. The model calculations are validated by comparing the simulated ENA and ion flux intensities with TWINS ENA images and in-situ ion data from THEMIS satellites. Simulation with static magnetic field produces a strong shielding electric field that skews the ion drift trajectories toward dawn. The model's corresponding peak ENA emissions are always eastward than those in the observed TWINS images. On the other hand, simulation with a dynamic magnetic field gives better spatial agreements with both ENA and insitu particle data, suggesting that temporal variations of the geomagnetic field exert a significant influence upon global ring current ion dynamics.

Fok, Mei-Ching↗
NASA NTRS2011

Simulation and Twins Observations of the 22 July 2009 Storm

TWINS is the first mission to perform stereo imaging of the Earth's ring current. The magnetic storm on 22 July 2009 was at the time the largest storm observed since TWINS began routine stereo imaging in June 2008. On 22 July 2009, the Dst dropped to nearly .80 nT at 0700 and 1000 UT. During the main phase, and at the peak of the storm, TWINS 1 and 2 were near apogee and moving between predawn and postdawn local time. The energetic neutral atom (ENA) imagers on the two spacecraft captured the storm intensification and the formation of the partial ring current. The peak of the high-altitude ENA emissions was seen in the midnight-to-dawn local time sector. The development of this storm has been simulated using the comprehensive ring current model (CRCM) to understand and interpret the observed signatures. We perform CRCM runs with constant and time-varying magnetic field. The model calculations are validated by comparing the simulated ENA and ion flux intensities with TWINS ENA images and in situ ion data from a THEMIS satellite. Simulation with a static magnetic field produces a strong shielding electric field that skews the ion drift trajectories toward dawn. The model's corresponding peak ENA emissions are always more eastward than those in the observed TWINS images. On the other hand, the simulation with a dynamic magnetic field gives better spatial agreement with both ENA and in situ particle data, suggesting that temporal variations of the geomagnetic field exert a significant influence upon global ring current ion dynamics.

Fok, M.-C.↗
DOE OSTI2024

Towards interpretable Cryo-EM: disentangling latent spaces of molecular conformations

Molecules are essential building blocks of life and their different conformations (i.e., shapes) crucially determine the functional role that they play in living organisms. Cryogenic Electron Microscopy (cryo-EM) allows for acquisition of large image datasets of individual molecules. Recent advances in computational cryo-EM have made it possible to learn latent variable models of conformation landscapes. However, interpreting these latent spaces remains a challenge as their individual dimensions are often arbitrary. The key message of our work is that this interpretation challenge can be viewed as an Independent Component Analysis (ICA) problem where we seek models that have the property of identifiability. That means, they have an essentially unique solution, representing a conformational latent space that separates the different degrees of freedom a molecule is equipped with in nature. Thus, we aim to advance the computational field of cryo-EM beyond visualizations as we connect it with the theoretical framework of (nonlinear) ICA and discuss the need for identifiable models, improved metrics, and benchmarks. Moving forward, we propose future directions for enhancing the disentanglement of latent spaces in cryo-EM, refining evaluation metrics and exploring techniques that leverage physics-based decoders of biomolecular systems. Moreover, we discuss how future technological developments in time-resolved single particle imaging may enable the application of nonlinear ICA models that can discover the true conformation changes of molecules in nature. The pursuit of interpretable conformational latent spaces will empower researchers to unravel complex biological processes and facilitate targeted interventions. This has significant implications for drug discovery and structural biology more broadly. More generally, latent variable models are deployed widely across many scientific disciplines. Thus, the argument we present in this work has much broader applications in AI for science if we want to move from impressive nonlinear neural network models to mathematically grounded methods that can help us learn something new about nature.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Search for electroweak production of vector-like leptons in $\tau$-lepton and b-jet final states in pp collisions at $\sqrt{s}$ = 13 TeV with the ATLAS detector

A search for pair-production of vector-like leptons is presented, considering their decays into a third-generation Standard Model (SM) quark and a vector leptoquark ( U 1 ) as predicted by an ultraviolet-complete extension of the SM, referred to as the ‘4321’ model. Given the assumed decay of U 1 into third-generation SM fermions, the final state can contain multiple τ-leptons and b-quarks. This search is based on a dataset of pp collisions at $\sqrt{s}=13$ TeV recorded with the ATLAS detector during Run 2 of the Large Hadron Collider, corresponding to an integrated luminosity of up to $140~\textrm{fb}^{-1}$. No significant excess above the SM background prediction is observed, and 95% confidence level limits on the cross-section times branching ratio are derived as a function of the vector-like lepton mass. A lower observed (expected) limit of 910 GeV (970 GeV) is set on the vector-like lepton mass. Additionally, the results are interpreted for a supersymmetric model with an R-parity violating coupling to the third-generation quarks and leptons. Lower observed (expected) limits are obtained on the higgsino mass at 880 GeV (940 GeV) and on the wino mass at 1170 GeV (1170 GeV).

Aad, G. [CPPM, Aix-Marseille Université, CNRS/IN2P↗