Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Kalman filtering”

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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
NASA NTRS2009

GPS Navigation for the Magnetospheric Multi-Scale Mission

In 2014. NASA is scheduled to launch the Magnetospheric Multiscale Mission (MMS), a four-satellite formation designed to monitor fluctuations in the Earth's magnetosphere. This mission has two planned phases with different orbits (1? x 12Re and 1.2 x 25Re) to allow for varying science regions of interest. To minimize ground resources and to mitigate the probability of collisions between formation members, an on-board orbit determination system consisting of a Global Positioning System (GPS) receiver and crosslink transceiver was desired. Candidate sensors would be required to acquire GPS signals both below and above the constellation while spinning at three revolutions-per-minute (RPM) and exchanging state and science information among the constellation. The Intersatellite Ranging and Alarm System (IRAS), developed by Goddard Space Flight Center (GSFC) was selected to meet this challenge. IRAS leverages the eight years of development GSFC has invested in the Navigator GPS receiver and its spacecraft communication expertise, culminating in a sensor capable of absolute and relative navigation as well as intersatellite communication. The Navigator is a state-of-the-art receiver designed to acquire and track weak GPS signals down to -147dBm. This innovation allows the receiver to track both the main lobe and the much weaker side lobe signals. The Navigator's four antenna inputs and 24 tracking channels, together with customized hardware and software, allow it to seamlessly maintain visibility while rotating. Additionally, an extended Kalman filter provides autonomous, near real-time, absolute state and time estimates. The Navigator made its maiden voyage on the Space Shuttle during the Hubble Servicing Mission, and is scheduled to fly on MMS as well as the Global Precipitation Measurement Mission (GPM). Additionally, Navigator's acquisition engine will be featured in the receiver being developed for the Orion vehicle. The crosslink transceiver is a 1/4 Watt transmitter utilizing a TDMA schedule to distribute a science quality message to all constellation members every ten seconds. Additionally the system generates one-way range measurements between formation members which is used as input to the Kalman filter. In preparation for the MMS Preliminary Design Review (PDR), the Navigator was required to pass a series of Technology Readiness Level (TRL) tests to earn the necessary TRL-6 classification. The TRL-6 level is achieved by demonstrating a prototype unit in a relevant end-to-end environment. The IRAS unit was able to meet all requirements during the testing phase, and has thus been TRL-6 qualified

Bamford, William↗
NASA NTRS2003

Nonlinear Observers for Gyro Calibration

High precision estimation and control algorithms, to achieve unprecedented levels of pointing accuracy, will be required to support future formation flying missions such as interferometry missions. Achieving high pointing accuracy requires precise knowledge of the spacecraft rotation rate. Typically, the rotation rate is measured by a gyro. The measured rates can be corrupted by errors in alignment and scale factor, gyro biases, and noise. In this work, we present nonlinear observers for gyro calibration. Nonlinear observers are superior to extended or pseudo-linear Kalman filter type approaches for large errors and global stability. Three nonlinear gyro calibration observers are developed. The first observer estimates a constant gyro bias. The second observer estimates scale factor errors. The third observer estimates the gyro alignment for three orthogonal gyros. The convergence properties of all three observers are discussed. Additionally, all three observers are coupled with a nonlinear control algorithm. The stability of each of the resulting closed loop systems is analyzed. The observers are then combined, and the gyro calibration parameters are estimated simultaneously. The stability of the combined observers is addressed, as well as the stability of the resulting closed loop systems. Simulated test results are presented for each scenario. Finally, the nonlinear observers are compared to a pseudo-linear Kalman filter.

Thienel, Julie↗
NASA NTRS1994

The short-term prediction of universal time and length of day using atmospheric angular momentum

The ability to predict short-term variations in the Earth's rotation has gained importance in recent years owing to more precise spacecraft tracking requirements. Universal time (UT1), that component of the Earth's orientation corresponding to the rotation angle, can be measured by number of high-precision space geodetic techniques. A Kalman filter developed at the Jet Propulsion Laboratory (JPL) optimally combines these different data sets and generates a smoothed times series and a set of predictions for UT1, as well as for additional Earth orientation components. These UT1 predictions utilize an empirically derived random walk stochastic model for the length of the day (LOD) and require frequent and up-to-date measurements of either UT1 or LOD to keep errors from quickly accumulating. Recent studies have shown that LOD variations are correlated with changes in the Earth's axial atmospheric angular momentum (AAM) over timescales of several years down to as little as 8 days. AAM estimates and forecasts out to 10 days are routinely available from meteorological analysis centers; these data can supplement geodetic measurements to improve the short-term prediction of LOD and have therefore been incorporated as independent data types in the JPL Kalman filter. We find that AAM and, to a lesser extent, AAM forecast data are extremely helpful in generating accurate near-real-time estimates of UT1 and LOD and in improving short-term predictions of these quantities out to about 10 days.

Freedman, A. P.↗
NASA NTRS2009

Chemical Source Inversion using Assimilated Constituent Observations in an Idealized Two-dimensional System

We present a source inversion technique for chemical constituents that uses assimilated constituent observations rather than directly using the observations. The method is tested with a simple model problem, which is a two-dimensional Fourier-Galerkin transport model combined with a Kalman filter for data assimilation. Inversion is carried out using a Green's function method and observations are simulated from a true state with added Gaussian noise. The forecast state uses the same spectral spectral model, but differs by an unbiased Gaussian model error, and emissions models with constant errors. The numerical experiments employ both simulated in situ and satellite observation networks. Source inversion was carried out by either direct use of synthetically generated observations with added noise, or by first assimilating the observations and using the analyses to extract observations. We have conducted 20 identical twin experiments for each set of source and observation configurations, and find that in the limiting cases of a very few localized observations, or an extremely large observation network there is little advantage to carrying out assimilation first. However, in intermediate observation densities, there decreases in source inversion error standard deviation using the Kalman filter algorithm followed by Green's function inversion by 50% to 95%.

Tangborn, Andrew↗
NASA NTRS1999

Improved Earth Sensor Performance Using a Sequentially Correlated Noise Model

Spacecraft attitude estimation by means of an extended Kalman filter requires a reasonably true model of the inherent noise of each sensor. For some sensors, the largest uncorrected noise comes not from the sensor itself, but from errors in the model for the predicted observations. This is certainly the case for Earth horizon sensors. The Earth horizon as seen from low Earth orbit is nearly a circle whose radius depends primarily on altitude. A straightforward latitude-dependent correction is added to this to account for the oblateness of the Earth. There also are both seasonal and stochastic variations in the horizon height. The seasonal variations can be predicted to some limited degree based on models derived from historical data. The stochastic component characteristically shows variations that are correlated both in time and space but which are unpredictable over long time spans. This work investigates whether Earth horizon sensor performance can be improved by solving for its systematic error as an augmentation of an attitude Kalman filter. It is found that using only Earth and Sun sensors, the augmented state is not fully observable. Even when magnetometer data is included, only the pitch axis component of the error can be improved; the roll component is unobservable.

Sedlak, J.↗
NASA NTRS1999

Improved Earth Sensor Performance Using a Sequentially Correlated Noise Model

Spacecraft attitude estimation by means of an extended Kalman filter requires a reasonably true model of the inherent noise of each sensor. For some sensors, the largest uncorrected noise comes not from the sensor itself, but from errors in the model for the predicted observations. This is certainly the case for Earth horizon sensors. The Earth horizon as seen from low Earth orbit is nearly a circle whose radius depends primarily on altitude. A straightforward latitude-dependent correction is added to this to account for the oblateness of the Earth. There also are both seasonal and stochastic variations in the horizon height. The seasonal variations can be predicted to some limited degree based on models de(ived from historical data. The stochastic component characteristically shows variations that are correlated both in time and space but which are unpredictable over long time spans. This work investigates whether Earth horizon sensor performance can be improved by solving for its systematic error as an augmentation of an attitude Kalman filter. It is found that using only Earth and Sun sensors, the augmented state is not fully observable. Even when magnetometer data is included, only the pitch axis component of the error can be improved; the roll component is unobservable.

Sedlak, Joseph↗
NASA NTRS2012

Towards A Model-Based Prognostics Methodology for Electrolytic Capacitors: A Case Study Based on Electrical Overstress Accelerated Aging

A remaining useful life prediction methodology for electrolytic capacitors is presented. This methodology is based on the Kalman filter framework and an empirical degradation model. Electrolytic capacitors are used in several applications ranging from power supplies on critical avionics equipment to power drivers for electro-mechanical actuators. These devices are known for their comparatively low reliability and given their criticality in electronics subsystems they are a good candidate for component level prognostics and health management. Prognostics provides a way to assess remaining useful life of a capacitor based on its current state of health and its anticipated future usage and operational conditions. We present here also, experimental results of an accelerated aging test under electrical stresses. The data obtained in this test form the basis for a remaining life prediction algorithm where a model of the degradation process is suggested. This preliminary remaining life prediction algorithm serves as a demonstration of how prognostics methodologies could be used for electrolytic capacitors. In addition, the use degradation progression data from accelerated aging, provides an avenue for validation of applications of the Kalman filter based prognostics methods typically used for remaining useful life predictions in other applications.

Celaya, Jose R.↗
NASA NTRS2011

A Model-based Prognostics Methodology for Electrolytic Capacitors Based on Electrical Overstress Accelerated Aging

A remaining useful life prediction methodology for electrolytic capacitors is presented. This methodology is based on the Kalman filter framework and an empirical degradation model. Electrolytic capacitors are used in several applications ranging from power supplies on critical avionics equipment to power drivers for electro-mechanical actuators. These devices are known for their comparatively low reliability and given their criticality in electronics subsystems they are a good candidate for component level prognostics and health management. Prognostics provides a way to assess remaining useful life of a capacitor based on its current state of health and its anticipated future usage and operational conditions. We present here also, experimental results of an accelerated aging test under electrical stresses. The data obtained in this test form the basis for a remaining life prediction algorithm where a model of the degradation process is suggested. This preliminary remaining life prediction algorithm serves as a demonstration of how prognostics methodologies could be used for electrolytic capacitors. In addition, the use degradation progression data from accelerated aging, provides an avenue for validation of applications of the Kalman filter based prognostics methods typically used for remaining useful life predictions in other applications.

Electrolytic capacitors↗
NASA NTRS1989

Range Filtering For Navigation By Satellite

Less processing and storage of data needed. Mathematical basis developed for Kalman filtering of sequential measurements of range from single-channel receiver in Global Positioning System (GPS) to each of several navigation satellites. Range filtering involves separate applications of kinematic Kalman tracking filter to measurements of signal from each satellite in use. State vector consists of pseudorange and time derivatives. Main advantages: simplification of tuning and decreases in required amounts of storage and processing of data for navigation in ships, airplanes, and ground vehicles.

Paielli, Russell↗
NASA NTRS2021

Comparisons of IMERG Version 06 Precipitation at and between Passive Microwave Overpasses in the Tropics

The Integrated Multisatellite Retrievals for Global Precipitation Measurement Mission (IMERG) is a global precipitation product that uses precipitation retrievals from the virtual constellation of satellites with passive microwave (PMW) sensors, as available. In the absence of PMW observations, IMERG uses a Kalman filter scheme to morph precipitation from one PMW observation to the next. In this study, an analysis of convective systems observed during the Convective Process Experiment (CPEX) suggests that IMERG precipitation depends more strongly on the availability of PMW observations than previously suspected. Following this evidence, we explore systematic biases in IMERG through bulk statistics. In two CPEX case studies, cloud photographs, pilot’s radar, and infrared imagery suggest that IMERG represents the spatial extent of precipitation relatively well when there is a PMW observation but sometimes produces spurious precipitation areas in the absence of PMW observations. Also, considering an observed convective system as a precipitation object in IMERG, the maximum rain rate peaked during PMW overpasses, with lower values between them. Bulk statistics reveal that these biases occur throughout IMERG Version 06. We find that locations and times without PMW observations have a higher frequency of light precipitation rates and a lower frequency of heavy precipitation rates due to retrieval artifacts. These results reveal deficiencies in the IMERG Kalman filter scheme, which have led to the development of the Scheme for Histogram Adjustment with Ranked Precipitation Estimates in the Neighborhood (SHARPEN; described in a companion paper) that will be applied in the next version of IMERG.

Manikandan Rajagopal↗
NASA NTRS2019

GOES-R Series GEO Side-Lobe Capable GPSR Post-Launch Refinements and Operational Capabilities

This paper addresses three topics: 1) EOPP (EOP (Earth Orientation Prediction) Parameters) file modification, 2) Kalman filter parameter tuning regarding maneuvers and 3) off-pointing GPS (Global Positioning System) tracking capability. GOES-R (Geostationary Operational Environmental Satellite-R Series) is the first in a 4-part series of new weather satellites set to replace and upgrade the older GOES constellation. Two GOES-R series have been launched to date, GOES-S and GOES-R. GOES-R is operational over the Eastern United States and GOES-S over the West. The Global Positioning System Receiver (GPSR) on board this geostationary weather satellite is a mission critical enabling technology which has been both tested on the ground and evaluated on-orbit to verify its effectivity. Since becoming operational in November 2016, the GPSR onboard has performed extremely well under nominal circumstances. Further refinements regarding a variety of facets have taken place since the launch of GOES-R. One such refinement was the implementation of a modified EOP (Earth Orientation Prediction) parameter set to improve ECEF (Earth Centered Earth Fixed) to ECI (Earth Centered Inertial) transformation by restoring zonal tides removed from the EOP parameter fit per tech note 36. Another relevant refinement combined thermal consideration with Kalman filter tuning to improve orbit determination performance during maneuvers. Now with two years of data and two vehicles in orbit many capabilities of the GPSR have been identified and defined to a higher degree. For example, metrics on side-lobe tracking and off-Nadir tracking capabilities have been quantified to a high degree. This paper will seek to supplement the ESA (European Space Agency) GNC 2017 GOES-R GPSR performance paper as a deeper dive on specific tracking capabilities and performance improvements now implemented on the GOES-R and GOES-S vehicles.

GPS Receiver↗
NASA NTRSDate not supplied

Examination of L-Band Brightness Temperature Forecasts in the SMAP Level-4 Soil Moisture Analysis

The NASA Soil Moisture Active Passive (SMAP) mission [4] has been providing L-band (1.4 GHz) passive microwave brightness temperature (Tb) observations since April 2015. By assimilating the Tb observations into the NASA Catchment land surface model [5] using a spatially distributed ensemble Kalman filter, the NASA Global Modeling and Assimilation Office generates the SMAP Level-4 Soil Moisture (L4_SM) product, which provides global, 3-hourly, 9-km resolution estimates of surface (0-5 cm) and root-zone (0-100 cm) soil moisture with ~2.5-day latency for use in research and applications [6]. The L4_SM product also includes estimates of soil temperature, land surface fluxes, and assimilation diagnostics such as the model forecast and observed Tb values [7]. The output from the L4_SM system is routinely monitored by the L4_SM team. Such monitoring provides valuable information; instances of unusually large Tb observation-minus-forecast (O-F) residuals can indicate events during which soil moisture conditions are poorly described in the land modeling system [7,8]. For example, repeated occurrences of very large Tb O-F values in central Australia were traced back to deficiencies in the gauge-based precipitation product used in the land modeling system through L4_SM Version 5; this discovery prompted the use, beginning in Version 6, of satellite- and gauge-based precipitation observations outside of North America [9]. Recently, a statistical analysis of the Tb O-F residuals revealed systematic errors in the “tau-omega” L-band radiative transfer model that converts the land-model simulated soil moisture and temperature into the forecast Tb prior to the L4_SM analysis. The L4_SM analysis is built on the ensemble Kalman filter and assumes unbiased forecast errors. The presence of systematic Tb forecast errors could thus adversely impact the quality of the analyzed soil moisture. In this paper, we examine the Tb O-F residuals of the latest Version 7 L4_SM data (Science Version ID Vv7030 and Vv7032) [10,11,12]. In the L4_SM Version 7 algorithm, key parameters of the L-band radiative transfer model, including soil roughness, scattering albedo, and a (seasonally varying) climatology of vegetation opacity, are obtained from the SMAP Level-2 Radiometer retrieval product (Version 5) [1]. Additionally, we also use ground measurements of surface soil moisture and soil temperature from the SMAP core validation sites as in situ reference of soil conditions [2,3]. As part of the L4_SM system calibration, the seasonally-varying bias between the model forecast Tb and the observed values is removed prior to the assimilation of the SMAP Tb observations [7]. That is, the L4_SM system is designed to only correct errors in synoptic-scale and interannual variations from the long-term mean seasonal cycle while maintaining the model’s (potentially erroneous) climatology. The Tb O-F residuals examined here are thus computed using SMAP Tb observations after they are rescaled to the mean seasonal cycle of Tb from the modeling system without SMAP data assimilation. Consequently, the long-term average of the Tb O-F residuals has a global mean of only 0.13 K and locally small values, ranging from -1 to 3 K. Despite the small time-average values of the Tb O-F residuals, the model forecast Tb was nevertheless found to exhibit undesirable systematic errors. At some locations, the time-average Tb O-F values strongly depend on surface soil moisture (SM). At the Yanco SMAP core validation site, for example, the Tb O-F residuals typically range from 5 to 15 K under dry soil moisture conditions (SM < 0.15 m3 m-3) yet are predominantly negative under wet soil moisture conditions (SM > 0.25 m3 m-3), with values ranging from 0 to -40 K. This results in soil moisture analysis increments that persistently make the soil drier under dry SM conditions and persistently make the soil wetter under wet SM conditions, suggesting an error in the dynamic range of the simulated Tb, soil moisture or soil temperature. In this paper, we describe the higher-order systematic Tb forecast errors in more detail, examine their impact on the L4_SM product quality, and explore potential avenues to further improve the L4_SM algorithm.

Rolf Reichle↗
NASA NTRS2006

TRMM On-Orbit Performance Reassessed After Control Change

The Tropical Rainfall Measuring Mission (TRMM) spacecraft, a joint mission between the U.S. and Japan, launched onboard an H-I1 rocket on November 27, 1997, and transitioned in August, 2001, from an average operating altitude of 350 kilometers to 402.5 kilometers. Due to problems using the Earth Sensor Assembly (ESA) at the higher altitude, TRMM switched to a backup attitude control mode. Prior to the orbit boost TRMM controlled pitch and roll to the local vertical using ESA measurements while using gyro data to propagate yaw attitude between yaw updates from the Sun sensors. After the orbit boost, a Kalman filter used 3-axis gyro data with Sun sensor and magnetometers to estimate onboard attitude. While originally intended to meet a degraded attitude accuracy of 0.7 degrees, the new control mode met the original 0.2 degree attitude accuracy requirement after improving onboard ephemeris prediction and adjusting the magnetometer calibration onboard. Independent roll attitude checks using a science instrument, the Precipitation Radar (PR) which was built in Japan, provided a novel insight into the pointing performance. The PR data helped identify the pointing errors after the orbit boost, track the performance improvements, and show subtle effects from ephemeris errors and gyro bias errors. It also helped identify average bias trends throughout the mission. Roll errors tracked by the PR from sample orbits pre-boost and post-boost are shown in Figure 1. Prior to the orbit boost, the largest attitude errors were due to occasional interference in the ESA. These errors were sometime larger than 0.2 degrees in pitch and roll, but usually less, as estimated from a comprehensive review of the attitude excursions using gyro data. Sudden jumps in the onboard roll show up as spikes in the reported attitude since the control responds within tens of seconds to null the pointing error. The PR estimated roll tracks well with an estimate of the roll history propagated using gyro data. After the orbit boost, the attitude errors shown by the PR roll have a smooth sine-wave type signal because of the way that attitude errors propagate with the use of gyro data. Yaw errors couple at orbit period to roll with 1/4 orbit lag. By tracking the amplitude, phase, and bias of the sinusoidal PR roll error signal, it was shown that the average pitch rotation axis tends to be offset from orbit normal in a direction perpendicular to the Sun direction, as shown in Figure 2 for a 200 day period following the orbit boost. This is a result of the higher accuracy and stability of the Sun sensor measurements relative to the magnetometer measurements used in the Kalman filter. In November, 2001 a magnetometer calibration adjustment was uploaded which improved the pointing performance, keeping the roll and yaw amplitudes within about 0.1 degrees.

Bilanow, Stephen↗
DOE OSTI2022

Rock Physics-Based Data Assimilation of Integrated Continuous Active-Source Seismic and Pressure Monitoring Data during Geological Carbon Storage

Summary There has been substantial controversy concerning the role of geological carbon storage (GCS) in sequestering anthropogenic carbon emissions to mitigate climate change and global warming. Arguments center on the inability to monitor a geological storage site precisely and continuously, especially highlighting the associated costs and spatiotemporal trade-offs when using conventional subsurface monitoring techniques (well logs, core samples, chemical tracers, and 4D seismics). Active surveillance of GCS sites is essential for managing and mitigating potential leaks but is also required by regulation. With the goal of enhancing the monitoring capability at GCS sites, we present a rock physics-based joint data assimilation model to study a popular GCS site at Cranfield, Mississippi, USA. Synthetic continuous active-source seismic monitoring (CASSM) data (in the form of Vp and Qp measurements) and wellbore pressure monitoring data are assimilated with an ensemble of reservoir realizations to monitor gas saturation and reservoir pressure changes over a period of 100 years. Synthetic seismic attributes are generated using rock physics models (RPMs) and wellbore pressure monitoring data are extracted from the ground truth. Two assimilation methods, ensemble Kalman filter (EnKF) and ensemble Kalman smoother (EnKS), are tested in an observation system simulation experiment (OSSE) environment to assess the prediction accuracy of the individual and composite observation systems. The joint monitoring system achieves more accurate estimates of gas saturation and pressure, across the time span from start of injection to end of forecast, as compared to a single type of monitoring tool and irrespective of data assimilation algorithm choice. These results indicate that jointly assimilated data from two types of sensors (in this case, crosswell seismic and downhole pressure) may lead to a more risk-reducing monitoring design. One would expect that more data, vis-à-vis inclusion of a new sensor type, will improve the accuracy of any GCS monitoring system. However, from a practical standpoint, one important question is whether such a gain in accuracy is worth the additional cost associated with the new sensor. This paper focuses on quantifying the gain in accuracy, such that a practitioner can answer this question.

Engineering↗
NASA NTRS1972

Adaptive control of stochastic linear systems with unknown parameters

The problem of optimal control of linear discrete-time stochastic dynamical system with unknown and, possibly, stochastically varying parameters is considered on the basis of noisy measurements. It is desired to minimize the expected value of a quadratic cost functional. Since the simultaneous estimation of the state and plant parameters is a nonlinear filtering problem, the extended Kalman filter algorithm is used. Several qualitative and asymptotic properties of the open loop feedback optimal control and the enforced separation scheme are discussed. Simulation results via Monte Carlo method show that, in terms of the performance measure, for stable systems the open loop feedback optimal control system is slightly better than the enforced separation scheme, while for unstable systems the latter scheme is far better.

Ku, R. T.↗
NASA NTRS1977

Attitude determination techniques for Nimbus 6 and Nimbus G

Techniques for measuring attitude deviations in the Nimbus 6 and Nimbus G 3-axis stabilized, zero-momentum, earth observation satellites are presented. The on-orbit telemetered data include scanner roll and pitch error signals, an illuminated DSAS which can read the azimuth and elevation of the sunline when the spacecraft is not in the umbra, the signal from the single gyro in the RMP, and the momenta of the reaction wheels. Mathematical methods for interpreting the measured data are described with attention to filter equations, especially the Kalman filter because of its ability to eliminate high frequency noise.

Das, A.↗