Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parallel state 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 37 records · Page 2

Discovering causal structure with reproducing-kernel Hilbert space ε -machines

We merge computational mechanics’ definition of causal states (predictively equivalent histories) with reproducing-kernel Hilbert space (RKHS) representation inference. The result is a widely applicable method that infers causal structure directly from observations of a system’s behaviors whether they are over discrete or continuous events or time. A structural representation—a finite- or infinite-state kernel ϵ-machine—is extracted by a reduced-dimension transform that gives an efficient representation of causal states and their topology. In this way, the system dynamics are represented by a stochastic (ordinary or partial) differential equation that acts on causal states. We introduce an algorithm to estimate the associated evolution operator. Paralleling the Fokker–Planck equation, it efficiently evolves causal-state distributions and makes predictions in the original data space via an RKHS functional mapping. We demonstrate these techniques, together with their predictive abilities, on discrete-time, discrete-value infinite Markov-order processes generated by finite-state hidden Markov models with (i) finite or (ii) uncountably infinite causal states and (iii) continuous-time, continuous-value processes generated by thermally driven chaotic flows. The method robustly estimates causal structure in the presence of varying external and measurement noise levels and for very high-dimensional data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

2020 International Conference on Smart Grids and Energy Systems (SGES)

This paper is concerned about phasor-based linear state estimation of three-phase power grids which contain non transposed lines and unsymmetrical branches. It builds upon our previous work where the linear measurement equations were decoupled into three independent sets via the symmetrical component transformation. Such a transformation may however not fully decouple the measurement equations if there are nontransposed transmission lines and unsymmetrical branches in the system. This paper addresses this limitation by deriving additive terms that compensate for the remaining coupling between sequence measurements after the application of the conventional transformation. It will be shown that general unbalanced threephase power grids can still be transformed into three decoupled subsystems and their state estimation solution can be obtained in parallel by the proposed approach.

Khalili, Ramtin↗

Quadratic optimal cooperative control synthesis with flight control application

An optimal control-law synthesis approach is presented that involves simultaneous solution for two cooperating controllers operating in parallel. One controller's structure includes stochastic state estimation and linear feedback of the state estimates, while the other controller involves direct linear feedback of selected system output measurements. This structure is shown to be optimal under the constraint of linear feedback of system outputs in one controller. Furthermore, it is appropriate for flight control synthesis where the full-state optimal stochastic controller can be adjusted to be representative of an optimal control model of the human pilot in a stochastic regulation task. The method is experimentally verified in the case of the selection of pitch-damper gain for optimum pitch tracking, where optimum implies the best subjective pilot rating in the task. Finally, results from application of the method to synthesize a controller for a multivariable fighter aircraft are presented, and implications of the results of this method regarding the optimal plant dynamics for tracking are discussed.

Schmidt, D. K.↗

Satellite Angular Rate Estimation From Vector Measurements

This paper presents an algorithm for estimating the angular rate vector of a satellite which is based on the time derivatives of vector measurements expressed in a reference and body coordinate. The computed derivatives are fed into a spacial Kalman filter which yields an estimate of the spacecraft angular velocity. The filter, named Extended Interlaced Kalman Filter (EIKF), is an extension of the Kalman filter which, although being linear, estimates the state of a nonlinear dynamic system. It consists of two or three parallel Kalman filters whose individual estimates are fed to one another and are considered as known inputs by the other parallel filter(s). The nonlinear dynamics stem from the nonlinear differential equation that describes the rotation of a three dimensional body. Initial results, using simulated data, and real Rossi X ray Timing Explorer (RXTE) data indicate that the algorithm is efficient and robust.

Azor, Ruth↗

Graph-Learning-Assisted State and Event Tracking for Solar-Penetrated Power Grids with Heterogeneous Data Sources

Unlike transmission systems, distribution systems do not typically contain sufficient metering to enable real-time state estimation. The lack of sufficient real-time measurements prohibits accurate and timely monitoring of the state of distribution systems. As a result, control and optimal operation of distribution systems, especially those containing large numbers of renewable generation units are not possible without proper data and information about the current state of the system. The main motivation of this project is to address this shortcoming by developing an approach which provides “predicted” real-time measurements so that they can be used to execute a distribution system state estimator. Thus, the objective of the project is to make the distribution systems fully observable, such that the hosting capacity for solar generation can be accurately estimated, and unnecessary solar curtailments can be avoided. In order to accomplish this goal, the project investigated the use of a grid-model-informed machine learning (ML) tool which integrates heterogeneous data streams obtained from AMI meters, SCADA as well as PMU measurements and created synchronous measurement snapshots for the state estimator (SE); and developed a hybrid robust SE which provides not only accurate state estimates but also real-time feedback for the ML model refinement.

14 SOLAR ENERGY↗

A Shift Selection Strategy for Parallel Shift-invert Spectrum Slicing in Symmetric Self-consistent Eigenvalue Computation

The central importance of large-scale eigenvalue problems in scientific computation necessitates the development of massively parallel algorithms for their solution. Recent advances in dense numerical linear algebra have enabled the routine treatment of eigenvalue problems with dimensions on the order of hundreds of thousands on the world’s largest supercomputers. In cases where dense treatments are not feasible, Krylov subspace methods offer an attractive alternative due to the fact that they do not require storage of the problem matrices. However, demonstration of scalability of either of these classes of eigenvalue algorithms on computing architectures capable of expressing massive parallelism is non-trivial due to communication requirements and serial bottlenecks, respectively. In this work, we introduce the SISLICE method: a parallel shift-invert algorithm for the solution of the symmetric self-consistent field (SCF) eigenvalue problem. The SISLICE method drastically reduces the communication requirement of current parallel shift-invert eigenvalue algorithms through various shift selection and migration techniques based on density of states estimation and k-means clustering, respectively. This work demonstrates the robustness and parallel performance of the SISLICE method on a representative set of SCF eigenvalue problems and outlines research directions that will be explored in future work.

97 MATHEMATICS AND COMPUTING↗

A Shift Selection Strategy for Parallel Shift-Invert Spectrum Slicing in Symmetric Self-Consistent Eigenvalue Computation

The central importance of large scale eigenvalue problems in scientific computation necessitates the development massively parallel algorithms for their solution. Recent advances in dense numerical linear algebra have enabled the routine treatment of eigenvalue problems with dimensions on the order of hundreds of thousands on the world's largest supercomputers. In cases where dense treatments are not feasible, Krylov subspace methods offer an attractive alternative due to the fact that they do not require storage of the problem matrices. However, demonstration of scalability of either of these classes of eigenvalue algorithms on computing architectures capable of expressing excessive parallelism is non-trivial due to communication requirements and serial bottlenecks, respectively. In this work, we introduce the SISLICE method: a parallel shift-invert algorithm for the solution of the symmetric self-consistent field (SCF) eigenvalue problem. The SISLICE method drastically reduces the communication requirement of current parallel shift-invert eigenvalue algorithms through various shift selection and migration techniques based on density of states estimation and k-means clustering, respectively. This work demonstrates the robustness and parallel performance of the SISLICE method on a representative set of SCF eigenvalue problems and outlines research directions which will be explored in future work.

97 MATHEMATICS AND COMPUTING↗

SELF-CONSISTENT ION CYCLOTRON ANISOTROPY-BETA RELATION FOR SOLAR WIND PROTONS

We derive a set of self-consistent marginally stable states for a system of ion-cyclotron waves propagating parallel to the large-scale magnetic field through a homogeneous proton-electron plasma. The proton distributions and the wave dispersions are related through the condition that no further ion-cyclotron resonant particle scattering or wave growth/damping may take place. The thermal anisotropy of the protons in these states therefore defines the threshold value for triggering the proton-cyclotron anisotropy instability. A number of recent papers have noted that the anisotropy of solar wind protons at 1 AU does not seem to be limited by the proton-cyclotron anisotropy threshold, even at low plasma beta. However, this puzzle seems to be due solely to the estimation of this anisotropy threshold under the assumption that the protons have a bi-Maxwellian distribution. We note that bi-Maxwellian distributions are never marginally stable to the resonant cyclotron interaction, so these estimates do not represent physically valid thresholds. The threshold anisotropies obtained from our marginally stable states are much larger, as a function of proton parallel beta, than the bi-Maxwellian estimates, and we show that the measured data remains below these more rigorous thresholds. Thus, the results of this paper resolve the apparent contradiction presented by the solar wind anisotropy observations at 1 AU: the bi-Maxwellian anisotropies are not rigorous thresholds, and so do not limit the proton distributions in the solar wind.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

On a useful functional representation of control system structure

An alternative structure for control systems is proposed. The structure is represented by a three-element block diagram and three functional definitions. It is argued that the three functional elements form a canonical set. The set includes the functions description, estimation and control. General overlay of the structure on parallel state and nested-state control systems is discussed. Breakdown of two real nested-state control systems into the proposed functional format is displayed. Application of the process to the mapping of complex control systems R and D efforts is explained with the Mars Rover Sample and Return mission as an example. A previous application of this basic functional structure to Space Station performance requirements organization is discussed.

Malchow, Harvey L.↗

Infrared emission from desorbed NO2(*) and NO(*)

Infrared photons from the radiative cascade accompany both the gas phase NO2 continuum chemiluminescence (which originates from its 2B2 and 2B1 states) and the NO beta bands. When these upper electronic states are excited by recombination/desorption at surfaces of low Earth orbiting spacecraft, similar IR emission spectrums will be observed. The principal NO2 features (other than the long wavelength tail of its electronic transitions) are the nu sub 3 fundamental sequence near 6.2 microns and nu sub 1 + nu sub 3 intercombination bands near 3.6 microns; NO would emit the delta v=1 and delta v=2 systems above 5.3 and 2.7 microns. Because of the long radiative lifetimes of the upper vibrational states, the infrared radiances in projections parallel to the vehicle surface (which we estimate) are substantially less than those of the visible and ultraviolet glows.

Kofsky, I. L.↗

Airborne radar technology for windshear detection

The objectives and accomplishments of the two-and-a-half year effort to describe how returns from on-board Doppler radar are to be used to detect the presence of a wind shear are reported. The problem is modeled as one of first passage in terms of state variables, the state estimates are generated by a bank of extended Kalman filters working in parallel, and the decision strategy involves the use of a voting algorithm for a series of likelihood ratio tests. The performance issue for filtering is addressed in terms of error-covariance reduction and filter divergence, and the performance issue for detection is addressed in terms of using a probability measure transformation to derive theoretical expressions for the error probabilities of a false alarm and a miss.

Hibey, Joseph L.↗

GridOPTICS/GridPACK

GridPACK is a software framework consisting of a set of modules designed to simplify the development of programs that model the power grid and run on parallel, high performance computing platforms. It also contains several fully developed applications, including powerflow, dynamic simulation, state estimation, Kalman filter analysis (dynamic state estimation), contingency analysis and real time path rating. These applications can be used either standalone or as components in more complicated workflows that combine several different types of application together. The framework modules are available as a combination of libraries and software templates and consist of components for setting up and distributing power grid networks, support for modeling the behavior of individual buses and branches in the network, converting the network models to the corresponding algebraic equations, and parallel routines for manipulating and solving large algebraic systems. The framework also contains a module for distributing tasks evenly amongst computing resources, even if individual tasks vary widely in their execution times. Additional modules support input and output, basic statistical analysis of contingency based calculations, distributed data structures, as well as basic profiling and error management.

Palmer, Bruce↗

Towards Resilient Autonomous Navigation of Drones

Robots and particularly drones are especially useful in exploring extreme environments that pose hazards to humans. To ensure safe operations in these situations, usually perceptually degraded and without good GNSS, it is critical to have a reliable and robust state estimation solution. The main body of literature in robot state estimation focuses on developing complex algorithms favoring accuracy. Typically, these approaches rely on a strong underlying assumption: the main estimation engine will not fail during operation. In contrast, we propose an architecture that pursues robustness in state estimation by considering redundancy and heterogeneity in both sensing and estimation algorithms. The architecture is designed to expect and detect failures and adapt the behavior of the system to ensure safety. To this end, we present HeRO (Heterogeneous Redundant Odometry): a stack of estimation algorithms running in parallel supervised by a resiliency logic. This logic carries out three main functions: a) perform confidence tests both in data quality and algorithm health; b) re-initialize those algorithms that might be malfunctioning; c) generate a smooth state estimate by multiplexing the inputs based on their quality. The state and quality estimates are used by the guidance and control modules to adapt the mobility behaviors of the system. The validation and utility of the approach are shown with real experiments on a ying robot for the use case of autonomous exploration of subterranean environments, with particular results from the STIX event of the DARPA Subterranean Challenge.

Agha-mohammadi, Ali-akbar↗

A meshless stochastic method for Poisson–Nernst–Planck equations

A plethora of biological, physical, and chemical phenomena involve transport of charged particles (ions). Its continuum-scale description relies on the Poisson–Nernst–Planck (PNP) system, which encapsulates the conservation of mass and charge. The numerical solution of these coupled partial differential equations is challenging and suffers from both the curse of dimensionality and difficulty in efficiently parallelizing. We present a novel particle-based framework to solve the full PNP system by simulating a drift–diffusion process with time- and space-varying drift. We leverage Green’s functions, kernel-independent fast multipole methods, and kernel density estimation to solve the PNP system in a meshless manner, capable of handling discontinuous initial states. The method is embarrassingly parallel, and the computational cost scales linearly with the number of particles and dimension. We use a series of numerical experiments to demonstrate both the method’s convergence with respect to the number of particles and computational cost vis-à-vis a traditional partial differential equation solver.

Chemistry↗

Technical Advances for Sustainable Remediation: In Situ Monitoring and Climate Resiliency - 18436

Recent advances in sustainable remediation techniques show promise as solutions for contaminant removal in soil and groundwater media. Although passive remediation (or natural attenuation) has gained momentum as an acceptable technique, there are still challenges in regulatory frameworks as well as in public acceptance. Leaving contaminants in subsurface requires increased burden of proof to show that plumes are stable and residual contaminants do not pose a significant health risk. Climate changes and frequent extreme climate events contribute another level of complexity and make it difficult to demonstrate system stability and predict future plume conditions. In this study, we show two coupled approaches to tackle this challenge: (1) in situ groundwater monitoring based on advanced sensor and data analytics technologies (i.e., Big Data), and (2) state-of-art modeling for optimizing monitoring and site closure strategies. In the new monitoring approach, in situ sensors are used to continuously measure the key variables that control contaminant plume mobility and its spatial and temporal distribution (such as pH, redox potential, electrical conductivity, and groundwater level). Coupled with a limited number of groundwater sampling, the data analytics methods - data mining and machine learning - allow us to estimate contaminant concentrations, and also to detect any significant changes associated with the plume mobility. In parallel, the state-of-art numerical flow and reactive transport model - ASCEM software under DOE Office of Environmental Management - provides a physical and mechanistic understanding of the system behavior, and predicts the long-term plume migration for optimizing the monitoring and closure strategies. In particular, we focus on the effects of extreme precipitation or drought on the site closure as well as monitoring strategies. We demonstrate this approach at the Savannah River Site (SRS) F-Area, where groundwater is contaminated by various radionuclides, including uranium, tritium and technetium. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Computing the partition function, ensemble averages, and density of states for lattice spin systems by sampling the mean

An algorithm to approximately calculate the partition function (and subsequently ensemble averages) and density of states of lattice spin systems through non-Monte-Carlo random sampling is developed. This algorithm (called the sampling-the-mean algorithm) can be applied to models where the up or down spins at lattice nodes interact to change the spin states of other lattice nodes, especially non-Ising-like models with long-range interactions such as the biological model considered here. Because it is based on the Central Limit Theorem of probability, the sampling-the-mean algorithm also gives estimates of the error in the partition function, ensemble averages, and density of states. Easily implemented parallelization strategies and error minimizing sampling strategies are discussed. The sampling-the-mean method works especially well for relatively small systems, systems with a density of energy states that contains sharp spikes or oscillations, or systems with little a priori knowledge of the density of states.

97 MATHEMATICS AND COMPUTING↗

Exploring the Connection Between Sampling Problems in Bayesian Inference and Statistical Mechanics

The Bayesian and statistical mechanical communities often share the same objective in their work - estimating and integrating probability distribution functions (pdfs) describing stochastic systems, models or processes. Frequently, these pdfs are complex functions of random variables exhibiting multiple, well separated local minima. Conventional strategies for sampling such pdfs are inefficient, sometimes leading to an apparent non-ergodic behavior. Several recently developed techniques for handling this problem have been successfully applied in statistical mechanics. In the multicanonical and Wang-Landau Monte Carlo (MC) methods, the correct pdfs are recovered from uniform sampling of the parameter space by iteratively establishing proper weighting factors connecting these distributions. Trivial generalizations allow for sampling from any chosen pdf. The closely related transition matrix method relies on estimating transition probabilities between different states. All these methods proved to generate estimates of pdfs with high statistical accuracy. In another MC technique, parallel tempering, several random walks, each corresponding to a different value of a parameter (e.g. "temperature"), are generated and occasionally exchanged using the Metropolis criterion. This method can be considered as a statistically correct version of simulated annealing. An alternative approach is to represent the set of independent variables as a Hamiltonian system. Considerab!e progress has been made in understanding how to ensure that the system obeys the equipartition theorem or, equivalently, that coupling between the variables is correctly described. Then a host of techniques developed for dynamical systems can be used. Among them, probably the most powerful is the Adaptive Biasing Force method, in which thermodynamic integration and biased sampling are combined to yield very efficient estimates of pdfs. The third class of methods deals with transitions between states described by rate constants. These problems are isomorphic with chemical kinetics problems. Recently, several efficient techniques for this purpose have been developed based on the approach originally proposed by Gillespie. Although the utility of the techniques mentioned above for Bayesian problems has not been determined, further research along these lines is warranted

Pohorille, Andrew↗

Thin-Film Thermal Conductivity Measurements Using Superconducting Nanowires

We present a simple experimental scheme for estimating the cryogenic thermal transport properties of thin films using superconducting nanowires. In a parallel array of nanowires, the heat from one nanowire in the normal state changes the local temperature around adjacent nanowires, reducing their switching current. Calibration of this change in switching current as a function of bath temperature provides an estimate of the temperature as a function of displacement from the heater. This provides a method of determining the contribution of substrate heat transport to the cooling time of superconducting nanowire single-photon detectors. Understanding this process is necessary for successful electrothermal modeling of superconducting nanowire systems.

77 NANOSCIENCE AND NANOTECHNOLOGY↗