Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Kalman filter (KF)”

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 19 records

On-Orbit Multi-Field Wavefront Control with a Kalman Filter

A document describes a multi-field wavefront control (WFC) procedure for the James Webb Space Telescope (JWST) on-orbit optical telescope element (OTE) fine-phasing using wavefront measurements at the NIRCam pupil. The control is applied to JWST primary mirror (PM) segments and secondary mirror (SM) simultaneously with a carefully selected ordering. Through computer simulations, the multi-field WFC procedure shows that it can reduce the initial system wavefront error (WFE), as caused by random initial system misalignments within the JWST fine-phasing error budget, from a few dozen micrometers to below 50 nm across the entire NIRCam Field of View, and the WFC procedure is also computationally stable as the Monte-Carlo simulations indicate. With the incorporation of a Kalman Filter (KF) as an optical state estimator into the WFC process, the robustness of the JWST OTE alignment process can be further improved. In the presence of some large optical misalignments, the Kalman state estimator can provide a reasonable estimate of the optical state, especially for those degrees of freedom that have a significant impact on the system WFE. The state estimate allows for a few corrections to the optical state to push the system towards its nominal state, and the result is that a large part of the WFE can be eliminated in this step. When the multi-field WFC procedure is applied after Kalman state estimate and correction, the stability of fine-phasing control is much more certain. Kalman Filter has been successfully applied to diverse applications as a robust and optimal state estimator. In the context of space-based optical system alignment based on wavefront measurements, a KF state estimator can combine all available wavefront measurements, past and present, as well as measurement and actuation error statistics to generate a Maximum-Likelihood optimal state estimator. The strength and flexibility of the KF algorithm make it attractive for use in real-time optical system alignment when WFC alone cannot effectively align the system.

Lou, John↗

The application of LQR synthesis techniques to the turboshaft engine control problem

A power turbine governor was designed for a recent-technology turboshaft engine coupled to a modern, articulated rotor system using Linear Quadratic Regulator (LQR) and Kalman Filter (KF) techniques. A linear, state-space model of the engine and rotor system was derived for six engine power settings from flight idle to maximum continuous. An integrator was appended to the fuel flow input to reduce the steady-state governor error to zero. Feedback gains were calculated for the system states at each power setting using the LQR technique. The main rotor tip speed state is not measurable, so a Kalman Filter of the rotor was used to estimate this state. The crossover of the system was increased to 10 rad/s compared to 2 rad/sec for a current governor. Initial computer simulations with a nonlinear engine model indicate a significant decrease in power turbine speed variation with the LQR governor compared to a conventional governor.

Pfeil, W. H.↗

The application of LQR synthesis techniques to the turboshaft engine control problem

A power turbine governor was designed for a recent-technology turboshaft engine coupled to a modern, articulated rotor system using Linear Quadratic Regulator (LQR) and Kalman Filter (KF) techniques. A linear, state-space model of the engine and rotor system was derived for six engine power settings from flight idle to maximum continuous. An integrator was appended to the fuel flow input to reduce the steady-state governor error to zero. Feedback gains were calculated for the system states at each power setting using the LQR technique. The main rotor tip speed state is not measurable, so a Kalman Filter of the rotor was used to estimate this state. The crossover of the system was increased to 10 rad/s compared to 2 rad/sec for a current governor. Initial computer simulations with a nonlinear engine model indicate a significant decrease in power turbine speed variation with the LQR governor compared to a conventional governor.

Pfeil, W. H.↗

Filtering Methods for Error Reduction in Spacecraft Attitude Estimation Using Quaternion Star Trackers

Precision attitude determination for recent and planned space missions typically includes quaternion star trackers (ST) and a three-axis inertial reference unit (IRU). Sensor selection is based on estimates of knowledge accuracy attainable from a Kalman filter (KF), which provides the optimal solution for the case of linear dynamics with measurement and process errors characterized by random Gaussian noise with white spectrum. Non-Gaussian systematic errors in quaternion STs are often quite large and have an unpredictable time-varying nature, particularly when used in non-inertial pointing applications. Two filtering methods are proposed to reduce the attitude estimation error resulting from ST systematic errors, 1) extended Kalman filter (EKF) augmented with Markov states, 2) Unscented Kalman filter (UKF) with a periodic measurement model. Realistic assessments of the attitude estimation performance gains are demonstrated with both simulation and flight telemetry data from the Lunar Reconnaissance Orbiter.

Calhoun, Philip C.↗

An Extension to the Kalman Filter for an Improved Detection of Unknown Behavior

The use of Kalman filter (KF) interferes with fault detection algorithms based on the residual between estimated and measured variables, since the measured values are used to update the estimates. This feedback results in the estimates being pulled closer to the measured values, influencing the residuals in the process. Here we present a fault detection scheme for systems that are being tracked by a KF. Our approach combines an open-loop prediction over an adaptive window and an information-based measure of the deviation of the Kalman estimate from the prediction to improve fault detection.

Benazera, Emmanuel↗

Estimating and Forecasting Time-Varying Groundwater Recharge in Fractured Rock: A State-Space Formulation with Preferential and Diffuse Flow to the Water Table

Rapid infiltration following precipitation may result in groundwater contamination from surface contaminants or pathogens. In fractured rock, contaminants can migrate rapidly to points of groundwater withdrawals. In contrast to the temporal availability of groundwater quality chemical indicators, meteorological and groundwater level observations are available in real-time to estimate time-varying recharge, which can act as a surrogate to identify periods of rapid infiltration that may indicate contamination susceptibility. Estimating recharge using methods, such as base-flow recession, unsaturated infiltration models, or Water-Table Fluctuations (WTF), cannot capitalize on currently available technologies and telecommunication infrastructure to conduct real-time recharge estimation at scales relevant to characterizing rapid infiltration. We present a linear, physics-based State-Space (SS) model of one-dimensional infiltration to estimate recharge, which includes preferential and diffuse-flow to the water table. The model can take advantage of real-time data for water-table altitude, precipitation, and evapotranspiration. Model parameters are calibrated over an observation period, and the Kalman Filter (KF) is subsequently applied to continuously update the observed (water-table altitude) and unobserved (groundwater recharge) system states and predict future states as new data become available. The SS/KF algorithm is demonstrated at the Masser Groundwater Recharge Site in Pennsylvania, USA and comparisons are made with recharge estimates from WTF methods. Model results indicate that the frequency of observations (daily versus sub-daily) dictates the allocation between preferential and diffuse flow. Additionally, because infiltration processes encompass many nonlinearities, model parameters estimated from observation periods need to be updated at least seasonally to account for changing recharge conditions.

groundwater, Recharge, infiltration, recursive est↗

Investigating the application of Kalman Filters for real-time accountancy in fusion fuel cycles

Here tritium accountancy in the fusion fuel cycle is a significant concern for the operation of commercial devices. It is expected that a limit on the maximum amount of tritium inventory in the system will be implemented, meaning that accountancy of the tritium inventory in the fuel cycle will need to be as accurate as possible. This is difficult since not all locations along the fuel cycle can benefit from the implementation of a tritium accountancy sensor and measurements will inherently contain error. Commercial fusion plants will also operate continuously, challenging current accountancy techniques that rely on static processes in well-controlled environments. In this paper a simple fusion fuel cycle concept is defined and the tritium inventory over time of each component is modelled using the Euler Approximation of a series of differential Equations in Python. This information is then used to simulate sensor measurements at specific points on the fuel cycle and then passed through a Kalman Filter (KF) to improve the accuracy of the true measurement of the sensors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A robust adaptive flightpath reconstruction technique

Computational schemes are presented that allow accurate reconstruction of an aircraft's flightpath in real-time. The reconstruction of the flightpath is formulated as a linear state reconstruction problem, which can be solved via Kalman filtering (KF) techniques. This imposes some conditions upon the flight-test equipment. A reliable square root covariance KF (SRCF) implementation is chosen and further developed into a fully adaptive flightpath reconstruction scheme. Therefore, the basic SRCF is modified in order to cope with several practical problems such as: the automatic control of the convergence of the recursive KF calculations, time varying zero-bias errors on the input signal of the system model used in the KF, and the changing aircraft dynamics owing to a change in reference flight condition. The developed solutions for these problems are all implemented in a numerically stable way, which guarantees the overall flightpath reconstruction scheme to be robust. Furthermore, some special features of the used system model are exploited to make the algorithmic implementation very efficient. An experimental simulation study using simulated flight test data demonstrated these different capabilities.

Verhaegen, M. H.↗

Geodesy by radio interferometry - Corrections to the IAU 1980 nutation series

The 1980 data on VLBI were analyzed using two techniques: weighted least squares (WLS) and Kalman filtering (KF), estimating corrections to the 14 terms in the nutation series that have the largest coefficients and that could be separated spectrally using the VLBI data available. The estimates of the coefficients from the two analyses are in good agreement, with the rms difference between the two sets being 0.07 milliarcsec. The largest corrections to the nutation amplitudes found were -(/1.89 + or - 0.17/ + i/0.49 + or 0.17/) milliarcsec (WLS) and -(/2.03 + or - 0.12/ + i/0.38 + or - 0.12/) milliarcsec (KF) for the retrograde annual nutation; for the prograde semiannual nutation, the largest corrections were (/0.45 + or - 0.13/ + i/0.31 + or - 0.13/) milliarcsec (WLS) and (/0.43 + or - 0.10/ + i/0.34 + or 0.10/) milliarcsec (KF). The corrections of the amplitudes for all other terms (both in-phase and out-of-phase) were less than 0.3 milliarcsec.

Herring, T. A.↗

Suboptimal schemes for atmospheric data assimilation based on the Kalman filter

This work is directed toward approximating the evolution of forecast error covariances for data assimilation. The performance of different algorithms based on simplification of the standard Kalman filter (KF) is studied. These are suboptimal schemes (SOSs) when compared to the KF, which is optimal for linear problems with known statistics. The SOSs considered here are several versions of optimal interpolation (OI), a scheme for height error variance advection, and a simplified KF in which the full height error covariance is advected. To employ a methodology for exact comparison among these schemes, a linear environment is maintained, in which a beta-plane shallow-water model linearized about a constant zonal flow is chosen for the test-bed dynamics. The results show that constructing dynamically balanced forecast error covariances rather than using conventional geostrophically balanced ones is essential for successful performance of any SOS. A posteriori initialization of SOSs to compensate for model - data imbalance sometimes results in poor performance. Instead, properly constructed dynamically balanced forecast error covariances eliminate the need for initialization. When the SOSs studied here make use of dynamically balanced forecast error covariances, the difference among their performances progresses naturally from conventional OI to the KF. In fact, the results suggest that even modest enhancements of OI, such as including an approximate dynamical equation for height error variances while leaving height error correlation structure homogeneous, go a long way toward achieving the performance of the KF, provided that dynamically balanced cross-covariances are constructed and that model errors are accounted for properly. The results indicate that such enhancements are necessary if unconventional data are to have a positive impact.

Todling, Ricardo↗

Calibration of a Spacecraft Gyro Quadruplet

This work presents a new approach to gyro calibration where, in addition to being used for computing attitude that is needed in the calibration process, the gyro outputs are also used as measurements in a Kalman filter. Gyro calibration as well as calibration of other instruments occurs in two steps. In the first step, the instrument error parameters are estimated. During the second stage, those errors are continuously removed from the gyro readings. In the classical approach to gyro calibration, the gyro outputs are used to maintain or compute body orientation rather than being used as measurements in the context of filtering. In inertial navigation, for example, gyro errors cause erroneous computation of velocity and position, and then when the latter are compared to measured velocity and position, a great portion of the computed velocity and position errors can be determined. The latter errors are then fed into a Kalman filter (KF) that uses the INS error model to infer on the gyro errors. Similarly, when applying the classical approach to spacecraft (SC) altitude determination, the gyro outputs are used to compute the attitude and then attitude measurements are used to determine the attitude errors, which again using a KF, indicates what the gyro errors are. In the approach adopted in this work the gyro outputs are used as angular rate measurements and are compared to estimated angular rate measurements. However, this approach requires the knowledge of the angular rate. In the past, the estimated angular rate was computed in a rather simplistic way assuming that the rate was constant. In the present work, the estimated angular rate is derived using a KF whose input can be any kind of attitude measurement, therefore the angular rate experienced by the SC can be continuously changing, and yet a good estimate of the rate, necessary for calibration, can be obtained.

Harman, Richard R.↗

Investigating On-Orbit Attitude Determination Anomalies for the Solar Dynamics Observatory Mission

The Solar Dynamics Observatory (SDO) was launched on February 11, 2010 from Kennedy Space Center on an Atlas V launch vehicle into a geosynchronous transfer orbit. SDO carries a suite of three scientific instruments, whose observations are intended to promote a more complete understanding of the Sun and its effects on the Earth's environment. After a successful launch, separation, and initial Sun acquisition, the launch and flight operations teams dove into a commissioning campaign that included, among other things, checkout and calibration of the fine attitude sensors and checkout of the Kalman filter (KF) and the spacecraft s inertial pointing and science control modes. In addition, initial calibration of the science instruments was also accomplished. During that process of KF and controller checkout, several interesting observations were noticed and investigated. The SDO fine attitude sensors consist of one Adcole Digital Sun Sensor (DSS), two Galileo Avionica (GA) quaternion-output Star Trackers (STs), and three Kearfott Two-Axis Rate Assemblies (hereafter called inertial reference units, or IRUs). Initial checkout of the fine attitude sensors indicated that all sensors appeared to be functioning properly. Initial calibration maneuvers were planned and executed to update scale factors, drift rate biases, and alignments of the IRUs. After updating the IRU parameters, the KF was initialized and quickly reached convergence. Over the next few hours, it became apparent that there was an oscillation in the sensor residuals and the KF estimation of the IRU bias. A concentrated investigation ensued to determine the cause of the oscillations, their effect on mission requirements, and how to mitigate them. The ensuing analysis determined that the oscillations seen were, in fact, due to an oscillation in the IRU biases. The low frequencies of the oscillations passed through the KF, were well within the controller bandwidth, and therefore the spacecraft was actually following the oscillating biases, resulting in movement of the spacecraft on the order of plus or minus 20 arcsec. Though this level of error met the ACS attitude knowledge requirement of [35, 70, 70] arcsec, 3 sigma, the desire of the ACS and instrument teams was to remove as much of the oscillation as possible. The Kearfott IRUs have an internal temperature controller, designed to maintain the IRU temperature at a constant temperature of approximately 70 C, thus minimizing the change in the bias drift and scale factors of the mechanical gyros. During ground testing of the observatory, it was discovered that the 83-Hz control cycle of the IRU heaters put a tremendous amount of stress on the spacecraft battery. Analysis by the power systems team indicated that the constant charge/discharge on the battery due to the IRU thermal control cycle could potentially limit the life of the battery. After much analysis, the decision was made not to run the internal IRU heaters. Analysis of on orbit data revealed that the oscillations in the IRU bias had a connection to the temperature of the IRU; changes in IRU temperature resulted in changes in the amplitude and period of the IRU biases. Several mitigating solutions were investigated, the result of which was to tune the KF with larger IRU noise assumptions which allows the KF to follow and correct for the time-varying IRU biases.

Vess, Melissa F.↗

Improving Surface PM2.5 Forecasts in the United States Using an Ensemble of Chemical Transport Model Outputs: 2. Bias Correction with Satellite Data for Rural Areas

This work serves as the second of a two-part study to improve surface PM2.5 forecasts in the continental U.S. through the integrated use of multi satellite aerosol optical depth (AOD) products (MODIS Terra/Aqua and VIIRS DT/DB), multichemical transport model (CTM) (GEOS-Chem, WRF-Chem, and CMAQ) outputs, and ground observations. In Part I of the study, an ensemble Kalman filter (KF) technique using three CTM outputs and ground observations was developed to correct forecast bias and generate a single best forecast of PM2.5 for next day over non rural areas that have surface PM2.5 measurements in the proximity of 125 km. Here, with AOD data, we extended the bias correction into rural areas where the closest air quality monitoring station is at least 125–300 km away. First, we ensembled all of satellite AOD products to yield the single best AOD. Second, we corrected daily PM2.5 in rural areas from multiple models through the AOD spatial pattern between these areas and non rural areas, referred to as “extended ground truth” or EGT, for the present day. Lastly, we applied the KF technique to reduce the forecast bias for next day using the EGT. Our results find that the ensemble of bias-corrected daily PM2.5 from three CTMs for both today and next day show the best performance. Together, the two-part study develops a multimodel and multi-AOD bias-correction technique that has the potential to improve PM2.5 forecasts in both rural and non rural areas in near real time, and be readily implemented at state levels.

Huanxin Zhang↗

Attitude Control Approach for Solar Cruiser, a Large, Deep Space Solar Sail

Solar Cruiser is a small satellite Technology Demonstration Mission (TDM) to mature solar sail propulsion technology using a solar sail larger than 1600 square meters, demonstrating performance as a propulsion system and a stable pointing platform for science observations in an artificial halo orbit sunward of the Sun-Earth Lagrange Point 1 (sub-L1). For the “sailcraft” to meet mission ob- jectives, there are several unique attitude control challenges that the Attitude De- termination and Control System (ADCS) must overcome. Large disturbance tor- ques, primarily due to sail deformations coupled with off-sun pointing angles, make it more difficult to maintain adequate controls performance and manage accumulated momentum on the control actuators. The large-amplitude, low- frequency flexible body modes of the sail, in concert with noisy sensors and ac- tuators, make it challenging to mitigate control-structure interactions and main- tain fine pointing capabilities. Stability and control performance during sail de- ployment is complicated by rapidly and widely varying inertias. The Solar Cruiser ADCS effectively addresses these challenges using a simple, traditional control system – including momentum management actuators that use bang- bang control to constrain internal accumulated momentum, a low-pass controls filter and an attitude determination Kalman Filter (KF) that blends multiple star tracker solutions, and a reaction wheel assembly (RWA) with controller gains tuned to specific configurations or inertias. Solar Cruiser’s attitude control ap- proach, and lessons learned from its development, establishes a state of the art of high value to future solar sail missions and other small spacecraft with large de- ployable structures operating in deep space.

solar sail↗

Attitude Control Approach for Solar Cruiser, a Large, Deep Space Solar Sail

Solar Cruiser is a small satellite Technology Demonstration Mission (TDM) to mature solar sail propulsion technology using a solar sail larger than 1600 square meters, demonstrating performance as a propulsion system and a stable pointing platform for science observations in an artificial halo orbit sunward of the Sun-Earth Lagrange Point 1 (sub-L1). For the “sailcraft” to meet mission objectives, there are several unique attitude control challenges that the Attitude Determination and Control System (ADCS) must overcome. Large disturbance torques, primarily due to sail deformations coupled with off-sun pointing angles, make it more difficult to maintain adequate controls performance and manage accumulated momentum on the control actuators. The large-amplitude, low-frequency flexible body modes of the sail, in concert with noisy sensors and actuators, make it challenging to mitigate control-structure interactions and maintain fine pointing capabilities. Stability and control performance during sail deployment is complicated by rapidly and widely varying inertias. The Solar Cruiser ADCS effectively addresses these challenges using a simple, traditional control system – including momentum management actuators that use bang-bang control to constrain internal accumulated momentum, a low-pass controls filter and an attitude determination Kalman Filter (KF) that blends multiple star tracker solutions, and a reaction wheel assembly (RWA) with controller gains tuned to specific configurations or inertias. Solar Cruiser’s attitude control approach, and lessons learned from its development, establishes a state of the art of high value to future solar sail missions and other small spacecraft with large deployable structures operating in deep space.

solar sail↗

Applying the Kalman filter particle method to strange and open charm hadron reconstruction in the STAR experiment

Here, we apply KF Particle, a Kalman Filter package for secondary vertex finding and fitting, to strange and open charm hadron reconstruction in heavy-ion collisions in the STAR experiment. Compared to the conventional helix swimming method used in STAR, the KF Particle method improves the reconstructed Λ, Ω and $D$ 0 significance considerably. At the same time, we demonstrate that Monte Carlo simulation with the STAR detector responses can well reproduce the topological variable distributions reconstructed in real data using the KF Particle method, therefore retaining good control on the reconstruction efficiency uncertainties for strange and open charm hadrons measurements in heavy-ion collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Heading Toward Launch with the Integrated Multi-Satellite Retrievals for GPM (IMERG)

The Day-l algorithm for computing combined precipitation estimates in GPM is the Integrated Multi-satellitE Retrievals for GPM (IMERG). We plan for the period of record to encompass both the TRMM and GPM eras, and the coverage to extend to fully global as experience is gained in the difficult high-latitude environment. IMERG is being developed as a unified U.S. algorithm that takes advantage of strengths in the three groups that are contributing expertise: 1) the TRMM Multi-satellite Precipitation Analysis (TMPA), which addresses inter-satellite calibration of precipitation estimates and monthly scale combination of satellite and gauge analyses; 2) the CPC Morphing algorithm with Kalman Filtering (KF-CMORPH), which provides quality-weighted time interpolation of precipitation patterns following cloud motion; and 3) the Precipitation Estimation from Remotely Sensed Information using Artificial Neural Networks using a Cloud Classification System (PERSIANN-CCS), which provides a neural-network-based scheme for generating microwave-calibrated precipitation estimates from geosynchronous infrared brightness temperatures. In this talk we summarize the major building blocks and important design issues driven by user needs and practical data issues. One concept being pioneered by the IMERG team is that the code system should produce estimates for the same time period but at different latencies to support the requirements of different groups of users. Another user requirement is that all these runs must be reprocessed as new IMERG versions are introduced. IMERG's status at meeting time will be summarized, and the processing scenario in the transition from TRMM to GPM will be laid out. Initially, IMERG will be run with TRMM-based calibration, and then a conversion to a GPM-based calibration will be employed after the GPM sensor products are validated. A complete reprocessing will be computed, which will complete the transition from TMPA.

Huffman, George J.↗

Dynamic Matrix Completion Based State Estimation in Distribution Grids

The power distribution network is undergoing tremendous transformation due to an increase in the penetration of renewable energy resources and electric vehicles. These changes have resulted in greater uncertainty and dynamics in the distribution grid states. Therefore, the ability to track and monitor system states has become a critical need for accurate and timely control actions. In this paper, we propose two dynamic sparsity-based state estimation approaches for distribution systems: (1) locally weighted matrix completion (LW-MC) and (2) Bayesian matrix completion with Kalman filter prediction (BMC-KF). The performance of the proposed dynamic state estimation strategies is compared with the classic/static matrix completion (static-MC) approach using the IEEE 37 and IEEE 123 bus test systems. Finally, results indicate that BMC-KF approach outperforms both LW-MC as well as static-MC even when 30% of the measurement data is available. Computational complexity associated with both approaches is quantified.

42 ENGINEERING↗