Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “distributed estimation”

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 577 records · Page 32

On the use of maximum likelihood estimation for the assembly of Space Station Freedom

Distributed parameter models of the Solar Array Flight Experiment, the Mini-MAST truss, and Space Station Freedom assembly are discussed. The distributed parameter approach takes advantage of (1) the relatively small number of model parameters associated with partial differential equation models of structural dynamics, (2) maximum-likelihood estimation using both prelaunch and on-orbit test data, (3) the inclusion of control system dynamics in the same equations, and (4) the incremental growth of the structural configurations. Maximum-likelihood parameter estimates for distributed parameter models were based on static compliance test results and frequency response measurements. Because the Space Station Freedom does not yet exist, the NASA Mini-MAST truss was used to test the procedure of modeling and parameter estimation. The resulting distributed parameter model of the Mini-MAST truss successfully demonstrated the approach taken. The computer program PDEMOD enables any configuration that can be represented by a network of flexible beam elements and rigid bodies to be remodeled.

Taylor, Lawrence W., Jr.↗

Reinforcement Learning for Load-balanced Parallel Particle Tracing

We explore an online reinforcement learning (RL) paradigm to dynamically optimize parallel particle tracing performance in distributed-memory systems. Our method combines three novel components: (1) a work donation algorithm, (2) a high-order workload estimation model, and (3) a communication cost model. First, we design an RL-based work donation algorithm. Our algorithm monitors workloads of processes and creates RL agents to donate data blocks and particles from high-workload processes to low-workload processes to minimize program execution time. The agents learn the donation strategy on the fly based on reward and cost functions designed to consider processes' workload changes and data transfer costs of donation actions. Second, we propose a workload estimation model, helping RL agents estimate the workload distribution of processes in future computations. Third, we design a communication cost model that considers both block and particle data exchange costs, helping RL agents make effective decisions with minimized communication costs. We demonstrate that our algorithm adapts to different flow behaviors in large-scale fluid dynamics, ocean, and weather simulation data. Our algorithm improves parallel particle tracing performance in terms of parallel efficiency, load balance, and costs of I/O and communication for evaluations with up to 16,384 processors.

Distributed and parallel particle tracing↗

Survey of the Frequency Dependent Latitudinal Distribution of the Fast Magnetosonic Wave Mode from Van Allen Probes Electric and Magnetic Field Instrument and Integrated Science Waveform Receiver Plasma Wave Analysis

We present a statistical survey of the latitudinal structure of the fast magnetosonic wave mode detected by the Van Allen Probes spanning the time interval of 21 September 2012 to 1 August 2014. We show that statistically, the latitudinal occurrence of the wave frequency (f) normalized by the local proton cyclotron frequency (f(sub cP)) has a distinct funnel-shaped appearance in latitude about the magnetic equator similar to that found in case studies. By comparing the observed E/B ratios with the model E/B ratio, using the observed plasma density and background magnetic field magnitude as input to the model E/B ratio, we show that this mode is consistent with the extra-ordinary (whistler) mode at wave normal angles (theta(sub k)) near 90 deg. Performing polarization analysis on synthetic waveforms composed from a superposition of extra-ordinary mode plane waves with theta(sub k) randomly chosen between 87 and 90 deg, we show that the uncertainty in the derived wave normal is substantially broadened, with a tail extending down to theta(sub k) of 60 deg, suggesting that another approach is necessary to estimate the true distribution of theta(sub k). We find that the histograms of the synthetically derived ellipticities and theta(sub k) are consistent with the observations of ellipticities and theta(sub k) derived using polarization analysis.We make estimates of the median equatorial theta(sub k) by comparing observed and model ray tracing frequency-dependent probability occurrence with latitude and give preliminary frequency dependent estimates of the equatorial theta(sub k) distribution around noon and 4 R(sub E), with the median of approximately 4 to 7 deg from 90 deg at f/f(sub cP) = 2 and dropping to approximately 0.5 deg from 90 deg at f/f(sub cP) = 30. The occurrence of waves in this mode peaks around noon near the equator at all radial distances, and we find that the overall intensity of these waves increases with AE*, similar to findings of other studies.

Boardsen, Scott A.↗

Out-of-distribution detection with non-parametric density estimation for models predicting processing history of uranium ore concentrates

The rapid advancement in machine learning (ML) and computer vision (CV) coincides with the growth of interest in deploying these ML/CV models in numerous fields from medicine to social science. Similar to those areas, we have witnessed a great number of works in materials science employing ML/CV models – neural networks in particular – in their studies in recent years. These models have proven to obtain accurate performance in various tasks. However, these models struggle to attain a similar performance when encountering test samples coming from a distribution that is different from the training set. More importantly, they fail without providing any warning to the users. Therefore, we propose a framework for detecting out-of-distribution (OOD) samples to alert users when a human intervention might be necessary in this work. Specifically, we explore the use of a non-parametric density estimation method to detect OOD samples. Here, we assess OOD detection capability of the proposed framework on ML models developed for categorizing precipitation routes of U 3 O 8 when encountering OOD datasets that contain samples (1) undergone different imaging acquisition process, (2) undergone different material synthesis process, and (3) different materials than ID set. Through those experiments, we achieve an average area under the receiver operating characteristic (AUROC) of at least 91% on average in detecting OOD samples. With minimal overhead cost and superior performance, the proposed framework enables a reliable and safe system when deploying in real-world scenarios.

Convolutional neural networks↗

Estimating Watershed Subsurface Permeability From Stream Discharge Data Using Deep Neural Networks

Subsurface permeability is a key parameter in watershed models that controls the contribution from the subsurface flow to stream flows. Since the permeability is difficult and expensive to measure directly at the spatial extent and resolution required by fully distributed watershed models, estimation through inverse modeling has had a long history in subsurface hydrology. The wide availability of stream surface flow data, compared to groundwater monitoring data, provides a new data source to infer soil and geologic properties using integrated surface and subsurface hydrologic models. As most of the existing methods have shown difficulty in dealing with highly nonlinear inverse problems, we explore the use of deep neural networks for inversion owing to their successes in mapping complex, highly nonlinear relationships. We train various deep neural network (DNN) models with different architectures to predict subsurface permeability from stream discharge hydrograph at the watershed outlet. The training data are obtained from ensemble simulations of hydrographs corresponding to an permeability ensemble using a fully-distributed, integrated surface-subsurface hydrologic model. The trained model is then applied to estimate the permeability of the real watershed using its observed hydrograph at the outlet. Our study demonstrates that the permeabilities of the soil and geologic facies that make significant contributions to the outlet discharge can be more accurately estimated from the discharge data. Their estimations are also more robust with observation errors. Compared to the traditional ensemble smoother method, DNNs show stronger performance in capturing the nonlinear relationship between permeability and stream hydrograph to accurately estimate permeability. Our study sheds new light on the value of the emerging deep learning methods in assisting integrated watershed modeling by improving parameter estimation, which will eventually reduce the uncertainty in predictive watershed models.

54 ENVIRONMENTAL SCIENCES↗

Reconnaissance of marine resources

A test area along the NW Coast of Africa was used during the Skylab mission to study the distribution of temperature and plankton. The S190B Earth Terrain Camera with a spectral film response of 0.4-0.7 micrometers allowed qualitative estimates of the distribution patterns of suspended material. Differentiation between inorganic particles and phytoplankton could be made by comparing the green band and the red band of the S190A Camera System. The pictorial display of data obtained from the S191 scanning radiometer in the 10-11 micrometer atmospheric window allowed a detailed interpretation of the temperature distribution in the area where cold upwelled water reaches the euphotic zone. The comparison between infrared data and the imageries taken simultaneously indicated the origin of the cold water as well as the pathway within the Canary current. A fish survey carried out almost simultaneously in the area, by echosounding, showed high correlation between the position of good fishing grounds and the distribution of plankton as detected by remote sensing detectors on Skylab.

Szekielda, K.-H.↗

The Australian experiment with ETS-V

Land-mobile satellite propagation measurements were implemented at L Band (1.5 GHz) in South-Eastern Australia during an 11 day period in October 1988. Transmissions (CW) from both the Japanese ETS-5 and INMARSAT Pacific geostationary satellites were accessed. Previous measurements in this series were performed at both L Band (1.5 GHz) and UHF (870 MHz) in Central Maryland, North-Central Colorado, and the southern United States. The objectives of the Australian campaign were to expand the data base acquired in the U.S. to another continent, to validate a U.S. derived empirical model for estimating the fade distribution, to establish the effects of directive antennas, to assess the isolation between co- and cross-polarized transmissions, to derive estimates of fade as well as non-fade durations, and to evaluate diversity reception. All these objectives were met.

Vogel, Wolfhard J.↗

A new approach to importance sampling for the simulation of false alarms

In this paper a modified importance sampling technique for improving the convergence of Importance Sampling is given. By using this approach to estimate low false alarm rates in radar simulations, the number of Monte Carlo runs can be reduced significantly. For one-dimensional exponential, Weibull, and Rayleigh distributions, a uniformly minimum variance unbiased estimator is obtained. For Gaussian distribution the estimator in this approach is uniformly better than that of previously known Importance Sampling approach. For a cell averaging system, by combining this technique and group sampling, the reduction of Monte Carlo runs for a reference cell of 20 and false alarm rate of lE-6 is on the order of 170 as compared to the previously known Importance Sampling approach.

Lu, D.↗

Development and Evaluation of the Raindrop Size Distribution Parameters for the NASA Global Precipitation Measurement Mission Ground Validation Program

The National Aeronautics and Space Administration Global Precipitation Measurement (GPM) mission ground validation program uses dual-polarization radar moments to estimate raindrop size distribution (DSD) parameters, the mass-weighted mean drop diameter Dmass, and normalized intercept parameter NW, to validate the GPM Core Observatory–derived DSD parameters. The disdrometer-based Dmass and NW are derived through empirical relationships between Dmass and differential reflectivity ZDR, and between NW, reflectivity ZH, and Dmass. This study employs large datasets collected from two-dimensional video disdrometers (2DVD) during six different field studies to derive the requisite empirical relationships. The uncertainty of the derived Dmass(ZDR) relationship is evaluated through comparisons of 2DVD-calculated and ZDR-estimated Dmass, where ZDR is calculated directly from 2DVD observations. Similarly, the uncertainty of the NW(ZH, Dmass) relationship is evaluated through 2DVD-calculated and Dmass and ZH-estimated NW, where Dmass and ZH are directly calculated from 2DVD observations. This study also presents the sensitivity of Dmass(ZDR) relationships to climate regime and to disdrometer type after developing three additional Dmass(ZDR) relationships from second-generation Particle Size Velocity (PARSIVEL2) disdrometer (P2) observations collected in the Pacific Northwest, in Iowa, and at Kwajalein Atoll in the tropical Pacific Ocean. The application of P2-derived Dmass(ZDR) relationship based on precipitation in the northwestern United States to P2 observations collected over the tropical ocean resulted in the highest error among comparisons of the three datasets.

Ali Tokay↗

Magnetic flux in modeled magnetic clouds at 1 AU and some specific comparisons to associated photospheric flux

In order to better understand the solar origins of magnetic clouds, statistical distributions of the estimated axial magnetic flux of 30 magnetic clouds at 1 AU, separated according to their occurrence during the solar cycle, were obtained and a comparison was made of the magnetic flux of a magnetic cloud to the aggregate flux of apparently associated photospheric magnetic flux tubes, for some specific cases. The 30 magnetic clouds comprise 12 cases from WIND, and the remainder from IMP-8, earlier IMPs, the International Sun-Earth Explorer (ISEE) 3 and HELIOS. The total magnetic flux along the cloud axis was estimated using a constant alpha, cylindrical, force-free flux rope model to determine cloud diameter and axial magentic field strength. The distribution of magentic fluxes for the 30 clouds is shown to be in the form of a skewed Gaussian.

Lepping, R. P.↗

Adaptive Metropolis Sampling with Product Distributions

The Metropolis-Hastings (MH) algorithm is a way to sample a provided target distribution pi(z). It works by repeatedly sampling a separate proposal distribution T(x,x') to generate a random walk {x(t)}. We consider a modification of the MH algorithm in which T is dynamically updated during the walk. The update at time t uses the {x(t' less than t)} to estimate the product distribution that has the least Kullback-Leibler distance to pi. That estimate is the information-theoretically optimal mean-field approximation to pi. We demonstrate through computer experiments that our algorithm produces samples that are superior to those of the conventional MH algorithm.

Wolpert, David H.↗

A distributed fault-detection and diagnosis system using on-line parameter estimation

The development of a model-based fault-detection and diagnosis system (FDD) is reviewed. The system can be used as an integral part of an intelligent control system. It determines the faults of a system from comparison of the measurements of the system with a priori information represented by the model of the system. The method of modeling a complex system is described and a description of diagnosis models which include process faults is presented. There are three distinct classes of fault modes covered by the system performance model equation: actuator faults, sensor faults, and performance degradation. A system equation for a complete model that describes all three classes of faults is given. The strategy for detecting the fault and estimating the fault parameters using a distributed on-line parameter identification scheme is presented. A two-step approach is proposed. The first step is composed of a group of hypothesis testing modules, (HTM) in parallel processing to test each class of faults. The second step is the fault diagnosis module which checks all the information obtained from the HTM level, isolates the fault, and determines its magnitude. The proposed FDD system was demonstrated by applying it to detect actuator and sensor faults added to a simulation of the Space Shuttle Main Engine. The simulation results show that the proposed FDD system can adequately detect the faults and estimate their magnitudes.

Guo, T.-H.↗

An incremental database access method for autonomous interoperable databases

We investigated a number of design and performance issues of interoperable database management systems (DBMS's). The major results of our investigation were obtained in the areas of client-server database architectures for heterogeneous DBMS's, incremental computation models, buffer management techniques, and query optimization. We finished a prototype of an advanced client-server workstation-based DBMS which allows access to multiple heterogeneous commercial DBMS's. Experiments and simulations were then run to compare its performance with the standard client-server architectures. The focus of this research was on adaptive optimization methods of heterogeneous database systems. Adaptive buffer management accounts for the random and object-oriented access methods for which no known characterization of the access patterns exists. Adaptive query optimization means that value distributions and selectives, which play the most significant role in query plan evaluation, are continuously refined to reflect the actual values as opposed to static ones that are computed off-line. Query feedback is a concept that was first introduced to the literature by our group. We employed query feedback for both adaptive buffer management and for computing value distributions and selectivities. For adaptive buffer management, we use the page faults of prior executions to achieve more 'informed' management decisions. For the estimation of the distributions of the selectivities, we use curve-fitting techniques, such as least squares and splines, for regressing on these values.

Roussopoulos, Nicholas↗

Subcentimeter-size particle distribution functions in planetary rings from Voyager radio and photopolarimeter occultation data

Analysis of measurements of the scattered and direct components of Voyager 1 radio occultation signals at 3.5 and 13 cm wavelengths yield estimates of the distribution functions of supracentimeter-size particles and thickness of relatively broad regions in Saturn's rings. If mearurements of signal amplitude at a shorter wavelength are combined with the previously analyzed data, the shape of the distribution functions characterizing the smaller particles can be constrained. If size distributions of arbitrary form were considered, many solutions are found that are consistent with the three available observations of signal amplitude. The best-fit power law was calculated to the three observations at three wavelengths for several of the embedded Saturn ringlets. Mie scattering theory predicts that the measured phase of the radio occultation signal is highly sensitive to particles ranging from 0.1 to 1.0 wavelengths in size, thus additional constraints on the subcentimeter-size distribution functions for both the Saturn and Uranus rings can in principle be derived from radio phase measurements.

Zebker, Howard A.↗

Multi-Edge Graph Convolutional Networks for Power Systems

The exponential electrification of transportation has contributed to highly intermittent load variations in the distribution grid. This uncertainty has raised challenges for distribution system operation and control. Accurate nodal voltage estimation is highly essential for the safe and reliable operation of the grid. Graph convolutional networks have been used in machine-learning-based models for power grid applications like voltage estimation for their ability to capture the network topology of the grid. This paper presents a novel multi-edge graph convolutional layer that considers resistance and reactance as edge attributes. This layer is created by modifying the message-passing function within the graph convolutional network. The novel layer is then used to create a multi-edge graph convolutional network-based surrogate model for estimating voltage in the distribution network with highly uncertain electric vehicle loads. Results indicate improved performance of the multi-edge graph convolutional network model when compared to a standard graph convolutional network model.

Ravi, Abhijith↗

Geometrical evidence for dark matter: X-ray constraints on the mass of the elliptical galaxy NGC 720

We describe (1) a new test for dark matter and alternate theories of gravitation based on the relative geometries of the X-ray and optical surface brightness distributions and an assumed form for the potential, of the optical light, (2) a technique to measure the shapes of the total gravitating matter and dark matter of an ellipsoidal system which is insensitive to the precise value of the temperature of the gas and to modest temperature gradients, and (3) a new method to determine the ratio of dark mass to stellar mass that is dependent on the functional forms for the visible star, gas and dark matter distributions, but independent of the distance to the galaxy or the gas temperature. We apply these techniques to X-ray data from the ROSAT Position Sensitive Proportional Counter (PSPC) of the optically flattened elliptical galaxy NGC 720; the optical isophotes have ellipticity epsilon approximately 0.40 extending out to approximately 120 sec. The X-ray isophotes are significantly elongated, epsilon = 0.20-0.30 for semimajor axis a approximately 100 sec. The major axes of the optical and X-ray isophotes are misaligned by approximately 30 deg +/- 15 deg. Spectral analysis of the X-ray data reveals no evidence of temperature gradients or anisotropies and demonstrates that a single-temperature plasma (T approximately 0.6 keV) having subsolar heavy element abundances and a two-temperature model having solar abundances describe the spectrum equally well. Considering only the relative geometries of the X-ray and optical surface brightness distributions and an assumed functional form for the potential of the optical light, we conclude that matter distributed like the optical light cannot produce the observed ellipticities of the X-ray isophotes, independent of the gas pressure, the gas temperature, and the value of the stellar mass; this comparison assumes a state of quasi-hydrostatic equilibrium so that the three-dimensional surfaces of the gas emissivity trace the three-dimensional isopotential surfaces -- we discuss the viability of this assumption for NGC 720. Milgrom's Modification of Newtonian Dynamics (MOND) cannot dispel this manifestation of dark matter. Hence, geometrical considerations require, without mention of pressure or temperature, the presence of an extended, massive dark matter halo in NGC 720. Employing essentially the technique of Buote & Canizares (1992; Buote 1992) we use the shape of the X-ray surface brightness to constrain the shape of the total gravitating matter. The total matter is modeled as either an oblate or prolate spheriod of constant shape and orientation having either a Ferrers (rho approximately r(exp -n)) or Hernquist density. Assuming the X-ray gas is in hydrostatic equilibrium, we construct a model X-ray gas distribution for various temperature profiles. We determine the ellipticity of the total gravitating matter to be epsilon approximately 0.50-0.70. Using the single-temperature model we estimate a total mass approximately (0.41-1.4) x 10(exp 12) h(sub 80) solar mass interior to the ellipsoid of semimajor axis 43.6 h(sub 80) kpc. Ferrers densities as steep as r(exp -3) do not fit the data, but the r(exp -2) and Hernquist models yield excellent fits. We estimate the mass distributions of the stars and the gas and fit the dark matter directly. For a given gas equation of state and functional forms for the visible stars, gas, and dark matter, these models yield a distance-independent and temperature-independent measurement of the ratio of dark mass to stellar mass M(sub DM)/M(sub stars). We estimate a minimum M(sub DM)/M(sub stars) greater than or equal to 4 which corresponds to a total mass slightly greater than that derived from the single-temperature models for distance D = 20h(sub 80) Mpc.

Buote, David A.↗

Risk-Aware Measurement Synchronization and Recovery for DSSE With Heterogeneous Data Sources

Power distribution systems are increasingly integrating heterogeneous sensors with varying data reporting rates and types, which pose challenges to achieving observability at the desired temporal resolution of distribution system state estimation (DSSE). Multisensor failures caused by extreme events exacerbate these issues, introducing substantial uncertainties into DSSE. This article proposes a novel solution to these challenges by ensuring high-resolution system observability despite heterogeneous data sources and multisensor failures. First, a deep learning architecture combining long short-term memory (LSTM) and graph convolutional network (GCN) is employed to synchronize meters with different reporting rates, aiming to achieve system observability. A random-walk-model-based approach is introduced to generate pseudo-measurements while properly characterizing their uncertainties under multisensor failures. Finally, a disaster-risk-informed observability metric (RiOM) is defined to quantify the uncertainty associated with state estimation results. The proposed framework offers deeper insights into the system observability on the fly compared with conventional analysis. The effectiveness of the framework is demonstrated on an IEEE standard test case and a large-scale real-world distribution feeder in mid-Minnesota in the U.S.

97 MATHEMATICS AND COMPUTING↗

Finding the magnetopause location using soft X-ray observations and a statistical inverse method

Variability in the location and shape of the dayside magnetopause is attributed to magnetic reconnection, a fundamental process that enables the transfer of mass, energy, and momentum from the solar wind into the magnetosphere. The spatial and temporal properties of the magnetopause, under varying solar and magnetospheric conditions, remain largely unknown because empirical studies using in-situ observations are challenging to interpret. Global wide field-of-view (FOV) imaging is the only means to simultaneously observe the spatial distribution of the plasma properties over the vast dayside magnetospheric region and, subsequently, quantify the energy transport from the interplanetary medium into the terrestrial magnetosphere. Two upcoming missions, ESA/CAS SMILE and NASA’s LEXI will provide wide-field imagery of the dayside magnetosheath in soft X-rays, an emission generated by charge exchange interactions between high charge-state heavy ions of solar wind origin and exospheric neutral atoms. High-cadence two-dimensional observations of the magnetosheath will allow the estimation of dynamic properties of its inner boundary, the magnetopause, and enable studies of its response to changes in the solar wind dynamic pressure and interplanetary magnetic field orientation. This work introduces a statistically-based estimation approach based on inverse theory to estimate the spatial distribution of magnetosheath soft X-ray emissivities and, with this, identify the location of the magnetopause over the Sun−Earth line. To do so, we simulate the magnetosheath structure using the MHD-based OpenGGCM model and generate synthetic soft X-ray images using LEXI’s orbit and attitude information. Our results show that 3-D estimations using the described statistically-based technique are robust against Poisson-distributed shot noise inherent to soft X-ray images. Also, our proposed methodology shows that the accuracy of both three-dimensional (3-D) estimation and the magnetopause standoff distance calculation highly depends on the observational point.

Gonzalo Cucho-Padin↗