SEARCH · Search NASA
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.
This paper describes relationships between Kalman-Bucy filter design parameters and the asymptotic eigenstructures of the resulting error dynamics. Simple, explicit formulas are developed which select noise intensity matrices to achieve pre-specified eigenstructures. The results also add basic understanding of the asymptotic properties of filters, particularly with respect to recently proposed robustness recovery processes.
Martingale theory and appropriate smoothing properties of Loeve (1953) have been used to develop a modified Gaussian second-order filter. The performance of the filter is evaluated through numerical simulation of a Jupiter flyby mission. The observations used in the simulation are on-board measurements of the angle between Jupiter and a fixed star taken at discrete time intervals. In the numerical study, the influence of each of the second-order terms is evaluated. Five filter algorithms are used in the simulations. Four of the filters are the modified Gaussian second-order filter and three approximations derived by neglecting one or more of the second-order terms in the equations. The fifth filter is the extended Kalman-Bucy filter which is obtained by neglecting all of the second-order terms.
This paper explores a class of multiple-model-based fault detection and identification (FDI) methods for bias-type faults in actuators and sensors. These methods employ banks of Kalman-Bucy filters to detect the faults, determine the fault pattern, and estimate the fault values, wherein each Kalman-Bucy filter is tuned to a different failure pattern. Necessary and sufficient conditions are presented for identifiability of actuator faults, sensor faults, and simultaneous actuator and sensor faults. It is shown that FDI of simultaneous actuator and sensor faults is not possible using these methods when all sensors have biases.
In this paper, a semiparametric modeling approach is introduced as a paradigm for addressing model error arising from unresolved physical phenomena. Our approach compensates for model error by learning an auxiliary dynamical model for the unknown parameters. Practically, the proposed approach consists of the following steps. Given a physics-based model and a noisy data set of historical observations, a Bayesian filtering algorithm is used to extract a time-series of the parameter values. Subsequently, the diffusion forecast algorithm is applied to the retrieved time-series in order to construct the auxiliary model for the time evolving parameters. The semiparametric forecasting algorithm consists of integrating the existing physics-based model with an ensemble of parameters sampled from the probability density function of the diffusion forecast. To specify initial conditions for the diffusion forecast, a Bayesian semiparametric filtering method that extends the Kalman-based filtering framework is introduced. In difficult test examples, which introduce chaotically and stochastically evolving hidden parameters into the Lorenz-96 model, we show that our approach can effectively compensate for model error, with forecasting skill comparable to that of the perfect model.
This is the docker image for Pacific Northwest National Laboratory's DSSE (distribution system state estimator) used for the demo of OEDI-SI platform. To support the operation of modern distribution systems, operators require real-time visibility into system states. Due to a lack of measurements and unbalanced operation, the state estimation in distribution systems is challenging. This submission is related to an OEDI-SI use case which demonstrates the application of an extended Kalman filter-based state estimator on an IEEE-123 bus system. The state estimator uses the measurements to generate voltage estimates for the system.
The 25-m aperture Cornell Caltech Atacama Telescope (CCAT) will have a primary mirror that is divided into 162 individual segments, each of which is equipped with 3 positioning actuators. This paper presents a mathematical description of the telescope, its actuators and sensors, and uses it to derive control laws for figure maintenance. A Kalman Filter-based Optical State Estimator is used to continuously estimate the aberrations of the telescope; these are used in a state-feedback controller to maintain image quality. This approach provides the means to correct for the optical effects of errors that occur in unactuated degrees of freedom, such as lateral translations of the segments. The control laws are exercised in Monte Carlo and covariance analysis, to bound the closed-loop performance of the telescope and to conduct control design trades.
The NSTX-U Shorted Turn Protection (STP) system is a vital safety feature designed to protect the tokamak’s coils during operations. It employs a Kalman filter-based algorithm to monitor coil currents and detect faults in real-time. If an anomaly is detected, the system swiftly terminates the tokamak pulse to prevent damage. Developed with Matlab and Simulink, the STP system has been validated through real-time simulations, ensuring its effectiveness for the upcoming upgrades of the NSTX-U facility. In conclusion, this article presents the theory for implementation of the STP algorithm and the results from simulation and real-time implementation.
Land-atmosphere carbon exchange is known to be extremely heterogeneous in arctic ice-wedge polygonal tundra regions. In this study, a Kalman filter-based method was developed to estimate the spatio-temporal dynamics of daytime average net ecosystem exchange (NEEday) at 0.5-m resolution over a 550 m by 700 m study site. We integrated multi-scale, multi-type datasets, including normalized difference vegetation indices (NDVIs) obtained from a novel automated mobile sensor system (or tram system) and a greenness index map obtained from airborne imagery. We took advantage of the significant correlations between NDVI and NEEday identified based on flux chamber measurements. The weighted average of the estimated NEEday within the flux-tower footprint agreed with the flux tower data in term of its seasonal dynamics. We then evaluated the spatial variability of the growing season average NEEday, as a function of polygon geomorphic classes; i.e., the combination of polygon types—which are known to present different degradation stages associated with permafrost thaw—and microtopographic features (i.e., troughs, centers and rims). Our study suggests the importance of considering microtopographic features and their spatial coverage in computing spatially aggregated carbon exchange.
A Kalman-Schmidt filter is used to estimate atmospheric and trajectory parameters for entry into the Venusian atmosphere. A significantly improved version of the Landing Trajectory Reconstruction (LTR) computer program, used to obtain the estimates is described. Major improvements involve precision and linearity control of numerical differencing, corrected perturbation modeling, and incorporation of a refractivity model. Important results show that gyroscopic data are not usable with LTR, that atmospheric properties are generally well estimated for altitudes less than 100 km, and that the two LTR modes of operation are complementary in performance.
Application of Kalman statistical filter to nonlinear mathematical model
Kalman-Schmidt filter stability in orbit calculation, dependence on transition matrix, initial covariance matrix, observations, and covariance of noise in measurements
Optimum bounding filters are derived for a specific version (steady state time-invariant with scalar observations) of the Kalman-Bucy filtering problem with inexactly known system parameters and for the Wiener filtering problem with inexactly known spectral densities. The designed filter obtains a bound on the actual error covariance which is not known, and it also prevents apparent divergence. Conditions are derived for the design of the optimum bounding filter within a permissible class of solutions; this turns out to be the min-max mean-square error filter for an extended class of solutions. The bounding filter can be of lower order than the original system, and a technique is devised for reducing the order of the filtering system and concurrently obtaining a figure of merit for its performance.
In this paper, a preliminary shadowgraph-based analysis of dust particles re-suspension due to loss of vacuum accident (LOVA) in ITER-like nuclear fusion reactors has been presented. Dust particles are produced through different mechanisms in nuclear fusion devices, one of the main issues is that dust particles are capable of being re-suspended in case of events such as LOVA. Shadowgraph is based on an expanded collimated beam of light emitted by a laser or a lamp that emits light transversely compared to the flow field direction. In the STARDUST facility, the dust moves in the flow, and it causes variations of refractive index that can be detected by using a CCD camera. The STARDUST fast camera setup allows to detect and to track dust particles moving in the vessel and then to obtain information about the velocity field of dust mobilized. In particular, the acquired images are processed such that per each frame the moving dust particles are detected by applying a background subtraction technique based on the mixture of Gaussian algorithm. The obtained foreground masks are eventually filtered with morphological operations. Finally, a multi-object tracking algorithm is used to track the detected particles along the experiment. For each particle, a Kalman filter-based tracker is applied; the particles dynamic is described by taking into account position, velocity, and acceleration as state variable. The results demonstrate that it is possible to obtain dust particles’ velocity field during LOVA by automatically processing the data obtained with the shadowgraph approach.
A software module which designs and tests controllers and filters in Kalman Estimator form, based on a polynomial state-space model is discussed. The user-friendly program employs an interactive graphics approach to simplify the design process. A variety of input methods are provided to test the effectiveness of the estimator. Utilities are provided which address important issues in filter design such as graphical analysis, statistical analysis, and calculation time. The program also provides the user with the ability to save filter parameters, inputs, and outputs for future use.
The structure of the upper atmosphere can be indirectly probed by light in order to determine the global density structure of ozone, aerosols, and neutral atmosphere. Scattered and directly transmitted light is measured by a satellite and is shown to be a nonlinear function of the state which is defined to be a point-wise decomposition of the density profiles. Dynamics are imposed on the state vector and a structured estimation problem is developed. The estimation of these densities is then performed using a linearized Kalman-Bucy filter and a linearized Kushner-Stratonovich filter.
The Dawn spacecraft spent over a year in orbit around Vesta from July 2011 through August 2012. In order to maintain the designated science reference orbits and enable the transfers between those orbits, precise and timely orbit determination was required. Challenges included low-thrust ion propulsion modeling, estimation of relatively unknown Vesta gravity and rotation models, track-ing data limitations, incorporation of real-time telemetry into dynamics model updates, and rapid maneuver design cycles during transfers. This paper discusses the dynamics models, filter configuration, and data processing implemented to deliver a rapid orbit determination capability to the Dawn project.
Terrain-relative navigation (TRN) encompasses a wide variety of algorithms that perform localization with respect to the terrain below a flying vehicle. In traditional approaches, measurements of the terrain are matched to a map carried onboard. This work presents a terrain-relative navigation filter with a position measurement inspired by neural activity associated with positioning in nature. Here, the filter is shown to produce accurate position measurements that outperform popular optimization and template matching methods given poor prior knowledge of the position. The proposed method is also better-suited to distributed implementation than optimization-based methods.
Thermal power plants experience cycling duty leading to the fatigue of the boiler and heat exchanger tubes. As a result, tube failures occur frequently in coal fired fleets leading to forced outages. Furthermore, because the tube leaks have been the major source of unwanted shutdowns and the number of outages is increasing, present work focuses on the detection and isolation of the leak in a subcritical boiler based upon the process data from a commercial coalfired power plant. The mass balance equation around the steam drum was analyzed using timeseries data collected from a 300 MW power plant. The ratio of the feed water mass flow rate to the steam mass flow rate was defined as a key parameter for detecting leaks. The difference in slope between the feedwater and steam mass flow rate during the normal and faulty operations was established as the upper control limit for real time monitoring. To reduce false alarm rates that arise when raw signal is directly compared against the threshold due to common process fluctuations, an optimal filter was derived for smoothing. It was found that the optimal filter reacted much more quickly to process changes than an exponential moving average filter, around 8 h earlier on average. Occurrence of relatively high false alarm rates even in the filtered responses was related to the cycling of the boiler from the base load condition. Variable threshold was established to keep false alarm rates to the minimum while maintaining the leak detection rate. Finally, the leak was located at the economizer and this could readily be isolated by investigating the magnitude of the mass flow rates ratio and the temperature at the economizer outlet.