Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “data inference”

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 613 records · Page 34

Aseismic crustal deformation in the Transverse Ranges of southern California

Using geodetic data at a plate boundary, the subsurface displacements that accompany plate motion were modeled as the sum of rigid block motion and the elastic effects of frictional drag between blocks in a model consisting of 12 blocks and 27 rectangular fault patches. The block motion inferred from the geodetic data has the same order of magnitude as the geologic estimates, and for many faults the agreement was excellent. The model implies about 20 mm/y of crustal shortening normal to the trend of the San Andreas fault.

Cheng, Abe↗

Quantum-Inspired Bayesian Sampling for Uncertainty Quantification and Machine Learning (Final Technical Report)

With increasing simulation and measurement data, machine learning and artificial intelligence have been widely used in computational decision-making of complex engineering systems. The resulting tools, such as uncertainty quantification solvers, reinforcement learning, and physics-informed machine learning, have achieved great success in critical DOE tasks such as material discovery and design, energy system modeling and control, and numerical weather and climate prediction. A core topic in scientific machine learning and artificial intelligence is Bayesian inference: given an observed data set, people want to estimate the posterior distribution of a (possibly large) number of hidden parameters. Due to the flexibility and weak assumptions, Bayesian sampling has been the mainstream Bayesian inference solvers despite the rapid progress of approximate Bayesian inference. Classical Bayesian sampling methods such as Markov-chain Monte Carlo suffer from a low-acceptance rate due to the random walk nature, therefore state-of-the-art techniques use Hamiltonian Monte Carlo and its variants to efficiently draw posterior samples in a high dimension. The key idea of Hamiltonian Monte Carlo and its variants is to simulate the Hamiltonian dynamics of a classical particle with a fixed mass, and their performance significantly degrades when the posterior distribution is highly spiky or has multiple modes. Leveraging the idea of quantum physics, this project has investigated new theory, algorithms and applications of Bayesian inference (especially Bayesian sampling). The main results include: (1) novel quantum-inspired Bayesian sampling methods that can lead to better accuracy for challenging multi-modal or spiky distributions, (2) more scalable machine learning framework leveraging tensor-compressed Bayesian inference, and (3) Bayesian and sampling approaches for verifying the robustness of continuous and binary neural networks.

97 MATHEMATICS AND COMPUTING↗

Integration of a DER Management System in Riverside. Final report

The tasks in this project covered various aspects, including algorithm development, algorithm integration into a commercial Active Network Management (ANM) platform, hardware-in-the-loop (HIL) testing in an industry-standard testing platform, pilot demonstration in Riverside, California, and also cost and benefit analysis. The DERMS platform in this project can host different algorithms developed on different platforms (e.g., MATLAB and Python) and it can interact with different hardware devices (e.g., different PV inverters, battery inverters, and different sensors). The DER control solution are based on an advanced model-free, layered, and clustered DER control paradigm. At the core of the DER control algorithms was the concept of Extremum Seeking (ES), which is a model-free probing-based control technique. The ES-based control algorithms were tested on major real-world inverters; both individually and in a cluster. It was shown that even legacy equipment (or when paired with a few additional advanced equipment) can support such advanced control. The monitoring algorithms utilize a heterogeneous set of legacy and advanced sensor measurements, such as behind-the-meter DER sensors, distribution-level Phase Measurement Units, distribution-substation Supervisory Control and Data Acquisition (SCADA), and line current sensors, with their limited availability; in order to infer practical network conditions. Sensor data are utilized to achieve resource forecasting, phase identification, and distribution system state estimation. The technology that was developed and demonstrated in this project could be transformational to utilities, including the smaller municipal utilities such as in Riverside, which may not have the resources to deploy advanced distribution system and DERMS solutions in order to support high penetration of solar power integration. This project created a real-world prototype to provide utilities with an assessment of smart grid monitoring and control technologies.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Measuring cosmological distances using cluster edges as a standard ruler

ABSTRACT The line-of-sight velocity dispersion profile of galaxy clusters exhibits a ‘kink’ corresponding to the spatial extent of orbiting galaxies. Because the spatial extent of a cluster is correlated with the amplitude of the velocity dispersion profile, we can utilize this feature as a gravity-calibrated standard ruler. Specifically, the amplitude of the velocity dispersion data allows us to infer the physical cluster size. Consequently, observations of the angular scale of the ‘kink’ in the profile can be translated into a distance measurement to the cluster. Assuming the relation between cluster radius and cluster velocity dispersion can be calibrated from simulations, we forecast that with existing data from the Sloan Digital Sky Survey we will be able to measure the Hubble constant with 3.0 per cent precision. Implementing our method with data from the Dark Energy Spectroscopic Instrument (DESI) will result in a 1.3 per cent measurement of the Hubble constant. Adding cosmological supernova data improves the uncertainty of the DESI measurement to 0.7 per cent. While these error estimates are statistical only, they provide strong motivation for pursuing the necessary simulation program required to characterize and calibrate the systematic uncertainties impacting our proposed measurement. Whether or not our proposed measurement can in fact result in competitive H0 constraints will depend on what the eventual systematics floor for this method is.

79 ASTRONOMY AND ASTROPHYSICS↗

Shedding light on shadow generalized parton distributions

The feasibility of extracting generalized parton distributions (GPDs) from deeply virtual Compton scattering (DVCS) data has recently been questioned because of the existence of an infinite set of so-called “shadow GPDs” (SGPDs). These SGPDs depend on the process and manifest as multiple solutions (at a fixed scale 𝑄 2 ) to the inverse problem that needs to be solved to infer GPDs from DVCS data. SGPDs therefore pose a significant challenge for extracting GPDs from DVCS data. With this motivation we study the extent to which QCD evolution can provide constraints on SGPDs. This is possible because the known classes of SGPDs begin to contribute to observables after evolution, and can then be constrained (at the input scale 𝑄$^{2}_{0}$) by data that has a finite 𝑄 2 range. The impact that SGPDs could have on determining the total angular momentum, pressure and sheer force distributions, and tomography is also discussed. We find that 𝑄 2 evolution, given empirical data over a wide range of skewness (𝜉 ) and 𝑄 2 scales, constrains the SGPDs studied at low 𝑥 and low 𝜉. If this finding generalizes to all SGPDs—and we have not been able to find any SGPDs with qualitatively different behavior—then this potentially makes extraction of GPDs from DVCS possible over a limited range of 𝑥 and 𝜉.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Mixture Model Framework for Traumatic Brain Injury Prognosis Using Heterogeneous Clinical and Outcome Data

Prognoses of Traumatic Brain Injury (TBI) outcomes are neither easily nor accurately determined from clinical indicators. This is due in part to the heterogeneity of damage inflicted to the brain, ultimately resulting in diverse and complex outcomes. Using a data-driven approach on many distinct data elements may be necessary to describe this large set of outcomes and thereby robustly depict the nuanced differences among TBI patients’ recovery. In this work, we develop a method for modeling large heterogeneous data types relevant to TBI. Our approach is geared toward the probabilistic representation of mixed continuous and discrete variables with missing values. The model is trained on a dataset encompassing a variety of data types, including demographics, blood-based biomarkers, and imaging findings. In addition, it includes a set of clinical outcome assessments at 3, 6, and 12 months post-injury. The model is used to stratify patients into distinct groups in an unsupervised learning setting. We use the model to infer outcomes using input data, and show that the collection of input data reduces uncertainty of outcomes over a baseline approach. In addition, we quantify the performance of a likelihood scoring technique that can be used to self-evaluate the extrapolation risk of prognosis on unseen patients.

97 MATHEMATICS AND COMPUTING↗

How well is Martian ozone inferred with reflectance spectroscopy?

The Mariner 9 UV spectrometer scanned from 2100 to 3500 Angstroms in one of its two spectral channels every 3 seconds with a spectral resolution of 15 Angstroms and an effective field-of-view of approximately 300 km(exp 2). The only gaseous absorption in the 2000 to 3000 Angstrom region was assumed to come from the Hartley band system of ozone, and therefore the amount of ozone was inferred by fitting this absorption feature with laboratory data of ozone absorption. Mars O3 as inferred from these spectra shows strong seasonal and latitudinal variation, with column abundances ranging from 0.2 microns at equatorial latitudes to 60 microns over the northern winter polar latitudes. The detectability limit of the spectrometer was approximately 3 microns. I use a radiative transfer model based on the discrete ordinate method to calculate synthetic radiance spectra. When typical amounts of dust and cloud are present, significant underestimation of O3 occurs. A factor of 3 times as much O3 is needed to generate the same spectrum for cloudy, dusty atmospheres as for a clear atmosphere. If the scattering properties of Martian clouds and dust were well known, then their appearance would not be a problem, as a model would be capable of retrieving the O3 abundance. However, these properties are not well known, which raises doubts about the effectiveness of the current UV spectroscopy technique used to measure O3. Spatial and temporal variability in temperature and water vapor have been claimed to account for the scatter of the data points. However, water vapor is a small source of odd hydrogen in the winter polar atmosphere, and may not account for most of the variability. Masking by clouds and dust may also account for some of the observed O3 variability, because the nature and opacity of the clouds and dust in the polar hood change dramatically in latitude and even on a day-to-day basis. As the maximum O3 abundance resides near the surface, spacecraft must be able to observe through the entire cloud and dust abundance in order to actually see the total O3 column abundance. If reflectance spectroscopy is used, as on Mariner 9, then the cloud and the airborne dust must be traversed twice; first by the incoming solar flux down to the surface, and then once again upon reflection from the surface out to the spacecraft.

Lindner, Bernhard Lee↗

Libpanda: A High Performance Library for Vehicle Data Collection

Cyber-Physical Systems (CPS) generally involve time-critical components due to physical dynamics, therefore necessitating high-performance subsystems. This is also true in data collection scenarios to infer physical phenomena. This paper covers Libpanda as an example of a component that has been designed to address performance issues in CPS implementations. Libpanda is a C++ library that interfaces software with a Comma.ai Panda device. Pandas are used for installation in modern vehicles to read the vehicle CAN bus, providing rich sensor data and limited vehicle control through message injection. The motivation to design lib-panda stems from the lack of performance in Python-based code that runs on inexpensive hardware like a Raspberry Pi. In such situations, Python code would result in utilizing 92% CPU while also dropping around 40% of the CAN packet due to bottlenecks. Without using different tools, inconsistent data collection means a loss of time-based vehicle state interpretation. Libpanda addresses these issues through implementation in a different language and implementation of different design paradigms involving asynchronous calls and multithreading. The Panda also features a GPS module that allows multiple instances to synchronize clocks for large-scale data collection scenarios. Libpanda has been designed with time-synchronization in mind to aid in the measurement of inter-vehicle dynamics. The performance improvements of libpanda have resulted in it becoming an important component in automotive dynamics research that requires a higher technical performance in large-scale experiments.

Bunting, Matthew↗

International Ultraviolet Explorer observations of amorphous hot galaxies

In order to better understand the nature of star formation processes in amorphous galaxies, short wavelength International Ultraviolet Explorer (IUE) spectra of galaxies NGC 1705 and NGC 1800 were obtained. The IUE data for NGC 1705 were of excellent quality while the low signal-to-noise NGC 1800 observation was useful only as a rough guide to the ultraviolet energy distribution. It was found that NGC 1705 contains a normal mix of OB stars, which is consistent with the nearly constant recent star formation rate inferred from new optical data. The NGC 1800 is likely to have similar properties, and blue galaxies with amorphous structures thus do not show evidence for anomalies in stellar mass distributions. The UV spectra of amorphous galaxies and a variety of other hot extragalactic stellar systems have similar characteristics, which suggests OB stellar populations often are homogeneous in their properties.

Lamb, S. A.↗

ERPS to Monitor Non-conscious Mentation

Event Related Brain Potentials (or ERPs) are extracted from the EEG that can be recorded between a pair of electrodes placed on a person's scalp. The EEG is recorded as a continual fluctuation in voltage. It is the results of the integration of the potential fields generated by a multitude of neuronal ensembles that are active as the brain goes about its business. Within this ongoing signal it is possible to distinguish voltage fluctuations that are triggered in neural structures by the occurrence of specific events. This activity, evoked as it is by an external event, is known as the Evoked, or Event Related, Potential. The ERPs provide a unique opportunity to monitor non-conscious mentation. The inferences that can be based on ERP data are described and the limits of these inferences are emphasized. This, however, will not be an exhaustive review of the use of ERPs in Engineering Psychology. The application, its scope, and its limitations will be illustrated by means of one example. This example is preceded by a brief technical introduction to the methodology used in the study of ERPs. The manner in which ERPs are used to study cognition is described.

Donchin, E.↗

System and method for creating expert systems

A system and method provides for the creation of a highly graphical expert system without the need for programming in code. An expert system is created by initially building a data interface, defining appropriate Mission, User-Defined, Inferred, and externally-generated GenSAA (EGG) data variables whose data values will be updated and input into the expert system. Next, rules of the expert system are created by building appropriate conditions of the rules which must be satisfied and then by building appropriate actions of rules which are to be executed upon corresponding conditions being satisfied. Finally, an appropriate user interface is built which can be highly graphical in nature and which can include appropriate message display and/or modification of display characteristics of a graphical display object, to visually alert a user of the expert system of varying data values, upon conditions of a created rule being satisfied. The data interface building, rule building, and user interface building are done in an efficient manner and can be created without the need for programming in code.

Hughes, Peter M.↗

Molybdenum Valence in Basaltic Silicate Melts: Effects of Temperature and Pressure

The metal-silicate partitioning behavior of molybdenum has been used as a test for equilibrium core formation hypotheses [for example, 1-6]. However, current models that apply experimental data to equilibrium core-mantle differentiation infer the oxidation state of molybdenum from solubility data or from multivariable coefficients from metal-silicate partitioning data [1,3,7]. Molybdenum, a multi-valent element with a valence transition near the fO2 of interest for core formation (approx.IW-2) will be sensitive to changes in fO2 of the system and silicate melt structure. In a silicate melt, Mo can occur in either 4+ or 6+ valence state, and Mo(6+) can be either octahedrally or tetrahedrally coordinated. Here we present X-ray absorption near edge structure (XANES) measurements of Mo valence in basaltic run products at a range of P, T, and fO2 and further quantify the valence transition of Mo.

Danielson, L. R.↗

Sources of field-aligned currents in the auroral plasma

Data from the Dynamics Explorer 1 High Altitude Plasma Instrument (HAPI) and magnetometer are used to investigate the sources of field-aligned currents in the nightside auroral zone. It is found that the formula developed by S. Knight predicts the field-aligned current density fairly accurately in regions where a significant potential drop can be inferred from the HAPI data; there are, however, regions in which the proportionality between potential drop and field-aligned current does not hold. In particular, occurrences of strong upward field-aligned current associated not with inverted-V events but instead with suprathermal bursts are noted. In addition, upward field-aligned currents are often observed to peak near the edges of inverted-V events, rather than in the center as would be predicted by Knight.

Marshall, J. A.↗

Impact of Scattering Model on Disdrometer Derived Attenuation Scaling

NASA Glenn Research Center (GRC), the Air Force Research Laboratory (AFRL), and the Politecnico di Milano (POLIMI) are currently entering the third year of a joint propagation study in Milan, Italy utilizing the 20 and 40 GHz beacons of the Alphasat TDP#5 Aldo Paraboni scientific payload. The Ka- and Q-band beacon receivers were installed at the POLIMI campus in June of 2014 and provide direct measurements of signal attenuation at each frequency. Collocated weather instrumentation provides concurrent measurement of atmospheric conditions at the receiver; included among these weather instruments is a Thies Clima Laser Precipitation Monitor (optical disdrometer) which records droplet size distributions (DSD) and droplet velocity distributions (DVD) during precipitation events. This information can be used to derive the specific attenuation at frequencies of interest and thereby scale measured attenuation data from one frequency to another. Given the ability to both predict the 40 gigahertz attenuation from the disdrometer and the 20 gigahertz time-series as well as to directly measure the 40 gigahertz attenuation with the beacon receiver, the Milan terminal is uniquely able to assess these scaling techniques and refine the methods used to infer attenuation from disdrometer data. In order to derive specific attenuation from the DSD, the forward scattering coefficient must be computed. In previous work, this has been done using the Mie scattering model, however, this assumes a spherical droplet shape. The primary goal of this analysis is to assess the impact of the scattering model and droplet shape on disdrometer-derived attenuation predictions by comparing the use of the Mie scattering model to the use of the T-matrix method, which does not assume a spherical droplet. In particular, this paper will investigate the impact of these two scattering approaches on the error of the resulting predictions as well as on the relationship between prediction error and rain rate.

Radio Transmission↗

The International Satellite Cloud Climatology Project (ISCCP) - The first project of the World Climate Research Programme

The first project of the World Climate Research Program is the International Satellite Cloud Climatology Project, (ISCCP) whose objective is the collection and analysis of satellite radiance data in order to infer the global distribution of cloud radiative properties and improve the modeling of cloud effects on climate. The operational component of ISCCP takes advantage of the global coverage provided by the current and planned international array of geostationary and polar-orbiting meteorological satellites in the 1980s. It will produce a five-year global radiance and cloud data set. The research component of ISCCP will coordinate studies to validate climatology, improve cloud analysis algorithms, improve cloud effects modelling, and investigate the role of clouds in the atmospheric radiation budget and hydrologic cycle.

Schiffer, R. A.↗

Robust Deep Gaussian Process-Based Probabilistic Electrical Load Forecasting Against Anomalous Events

The abnormal events, such as the unprecedented COVID-19 pandemic, can significantly change the load behaviors, leading to huge challenges for traditional short-term forecasting methods. This article proposes a robust deep Gaussian processes (DGP)-based probabilistic load forecasting method using a limited number of data. Since the proposed method only requires a limited number of training samples for load forecasting, it allows us to deal with extreme scenarios that cause short-term load behavior changes. In particular, the load forecasting at the beginning of abnormal event is cast as a regression problem with limited training samples and solved by double stochastic variational inference DGP. The mobility data are also utilized to deal with the uncertainties and pattern changes and enhance the flexibility of the forecasting model. The proposed method can quantify the uncertainties of load forecasting outcomes, which would be essential under uncertain inputs. Extensive comparison results with other state-of-the-art point and probabilistic forecasting methods show that our proposed approach can achieve high forecasting accuracies with only a limited number of data while maintaining the excellent performance of capturing the forecasting uncertainties.

anomalous events↗

Mapping Isoprene Emissions over North America using Formaldehyde Column Observations from Space

I] We present a methodology for deriving emissions of volatile organic compounds (VOC) using space-based column observations of formaldehyde (HCHO) and apply it to data from the Global Ozone Monitoring Experiment (GOME) satellite instrument over North America during July 1996. The HCHO column is related to local VOC emissions, with a spatial smearing that increases with the VOC lifetime. lsoprene is the dominant HCHO precursor over North America in summer, and its lifetime (approx. = 1 hour) is sufficiently short that the smearing can be neglected. We use the Goddard Earth Observing System global 3-D model of tropospheric chemistry (GEOS-CHEM) to derive the relationship between isoprene emissions and HCHO columns over North America and use these relationships to convert the GOME HCHO columns to isoprene emissions. We also use the GEOS-CHEM model as an intermediary to validate the GOME HCHO column measurements by comparison with in situ observations. The GEOS-CHEM model including the Global Emissions Inventory Activity (GEIA) isoprene emission inventory provides a good simulation of both the GOME data (r(sup 2) = 0.69, n = 756, bias = +l1 %) and the in situ summertime HCHO measurements over North America (r(sup 2) = 0.47, n = 10, bias = -3%). The GOME observations show high values over regions of known high isoprene emissions and a day-to-day variability that is consistent with the temperature dependence of isoprene emission. Isoprene emissions inferred from the GOME data are 20% less than GEIA on average over North America and twice those from the U S . EPA Biogenic Emissions Inventory System (BEIS2) inventory. The GOME isoprene inventory when implemented in the GEOS-CHEM model provides a better simulation of the HCHO in situ measurements thaneitherGEIAorBEIS2 (r(sup 2) = 0.71,n= 10, bias = -10 %).

Palmer, Paul I.↗

Improvement of two-phase closure models in CTF using Bayesian inference

Under the Consortium for Advanced Simulation of Light Water Reactors (CASL) program, extensive capabilities have been developed in CTF to analyze light-water reactors (LWRs) for normal operating conditions, departure from nucleate boiling (DNB), and system transients. However, further improvements are required in the modeling and simulation of boiling water reactors (BWRs), which is a focus of the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. In this work, CTF validation results were used to optimize selected modeling coefficients by calibrating to experimental data using a Bayesian inference approach. Here, calibration studies were conducted to improve (vapor) void fraction prediction without worsening the two-phase pressure drop prediction, as well as to improve the two-phase pressure drop prediction. Calibration was performed for interfacial drag and wall shear models. Surrogates were developed to alleviate the computational expense required for sampling the parameter space using Markov chain Monte Carlo (MCMC). An assessment performed with calibrated models demonstrated an improvement of CTF in its prediction of key parameters such as void fraction and two-phase pressure drop.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗