Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “predictive 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 91 records · Page 5

On-Orbit Model Refinement for Controller Redesign

High performance control design for a flexible space structure is challenging since high fidelity plant models are difficult to obtain a priori. Uncertainty in the control design models typically require a very robust, low performance control design which must be tuned on-orbit to achieve the required performance. A new procedure for refining a multivariable open loop plant model based on closed-loop response data is presented. Using a minimal representation of the state space dynamics, a least squares prediction error method is employed to estimate the plant parameters. This control-relevant system identification procedure stresses the joint nature of the system identification and control design problem by seeking to obtain a model that minimizes the difference between the predicted and actual closed-loop performance. This paper presents an algorithm for iterative closed-loop system identification and controller redesign along with illustrative examples.

Whorton, Mark S.↗

InSight orbit determination

The InSight mission relied on accurate deep-space navigation for a successful Mars landing on November 26, 2018. In this paper, we discuss the role of the cruise Orbit Determination team, whose responsibilities included determining the spacecraft state, predicting the future trajectory, and quantifying the uncertainty associated with those estimates. In particular, we will focus on spacecraft dynamic modeling, small forces due to attitude control, radiometric tracking data, filter strategies, uncertainty quantification, and responses to unexpected flight situations. We will also provide analysis of reconstructed maneuvers, small forces, and delivery accuracy at Mars arrival.

Seubert, Jill↗

An efficient method for estimating neighboring steady-state numerical solutions to the Euler equations

The paper concentrates on a linear approximation method for predicting the changes occurring in steady-state numerical solutions of the Euler equations as a consequence of small changes in the independent variables which control the problem. The importance of proper boundary-condition treatment and other issues concerning the problem are covered along with the importance of proper algorithm selection for a fully supersonic inviscid flow. The method is applied to a subsonic nozzle involving variation of the pressure on the outflow boundary and to a supersonic inlet involving variation of the inflow Mach number. In the subsonic test case, the comparisons between the predicted and conventional numerical solutions are shown to be good, while in the supersonic test case, the agreement between the approximation method and conventional numerical solution starts out well but rapidly degenerates at some point in the flowfield as the perturbation of the boundary conditions is increased.

Taylor, Arthur C., III↗

Creep of a Silicon Nitride Under Various Specimen/Loading Configurations

Extensive creep testing of a hot-pressed silicon nitride (NC132) was performed at 1300 C in air using five different specimen/loading configurations, including pure tension, pure compression, four-point uniaxial flexure, ball-on-ring biaxial flexure, and ring-on-ring biaxial flexure. Nominal creep strain and its rate for a given nominal applied stress were greatest in tension, least in compression, and intermediate in uniaxial and biaxial flexure. Except for the case of compressive loading, nominal creep strain generally decreased with time, resulting in less-defined steady-state condition. Of the four different creep formulations - power-law, hyperbolic sine, step, redistribution models - the conventional power-law model still provides the most convenient and reasonable means to estimate simple, quantitative creep parameters of the material. Predictions of creep deformation for the case of multiaxial stress state (biaxial flexure) were made based on pure tension and compression creep data by using the design code CARES/Creep.

Choi, Sung R.↗

Analysis of surface pressure distributions on two elliptic missile configurations

The state-of-the-art methods for predicting missile aerodynamic characteristics do not accurately predict the loads of missile configurations with bodies of elliptic cross section. An investigation of this problem found significant nonlinear flow disturbance on the windward surface of a 3:1 elliptic body at Mach 2.50 in addition to the nonlinear vortical flows which develop on the leeside. A nonlinear full-potential flow method (NCOREL) was found to provide extremely accurate pressure estimates for attached-flow conditions and the vortex prediction method contained in the state-of-the-art method (NOSEVTX) was shown to accurately calculate body vortices and leeside pressures.

Allen, J. M.↗

The Principle of Energetic Consistency: Application to the Shallow-Water Equations

If the complete state of the earth's atmosphere (e.g., pressure, temperature, winds and humidity, everywhere throughout the atmosphere) were known at any particular initial time, then solving the equations that govern the dynamical behavior of the atmosphere would give the complete state at all subsequent times. Part of the difficulty of weather prediction is that the governing equations can only be solved approximately, which is what weather prediction models do. But weather forecasts would still be far from perfect even if the equations could be solved exactly, because the atmospheric state is not and cannot be known completely at any initial forecast time. Rather, the initial state for a weather forecast can only be estimated from incomplete observations taken near the initial time, through a process known as data assimilation. Weather prediction models carry out their computations on a grid of points covering the earth's atmosphere. The formulation of these models is guided by a mathematical convergence theory which guarantees that, given the exact initial state, the model solution approaches the exact solution of the governing equations as the computational grid is made more fine. For the data assimilation process, however, there does not yet exist a convergence theory. This book chapter represents an effort to begin establishing a convergence theory for data assimilation methods. The main result, which is called the principle of energetic consistency, provides a necessary condition that a convergent method must satisfy. Current methods violate this principle, as shown in earlier work of the author, and therefore are not convergent. The principle is illustrated by showing how to apply it as a simple test of convergence for proposed methods.

Cohn, Stephen E.↗

Accumulation of the planets

In modeling the accumulation of planetesimals into planets, it is appropriate to distinguish between two stages: an early stage, during which approximately 10 km diameter planetesimals accumulate locally to form bodies approximate 10 to the 25th g in mass; and a later stage in which the approximately 10 to the 25th g planetesimals accumulate into the final planets. In the terrestrial planet region, an initial planetesimal swarm corresponding to the critical mass of dust layer gravitational instabilities is considered. In order to better understand the accumulation history of Mercury-sized bodies, 19 Monte-Carlo simulations of terrestrial planet growth were calculated. A Monte Carlo technique was used to investigate the orbital evolution of asteroidal collision debris produced interior to 2.6 AU. It was found that there are two regions primarily responsible for production of Earth-crossing meteoritic material and Apollo objects. The same techniques were extended to include the origin of Earth-approaching asteroidal bodies. It is found that these same two resonant mechanisms predict a steady-state number of Apollo-Amor about 1/2 that estimated based on astronomical observations.

Wetherill, G. W.↗

Control of Thermo-Acoustics Instabilities: The Multi-Scale Extended Kalman Approach

"Multi-Scale Extended Kalman" (MSEK) is a novel model-based control approach recently found to be effective for suppressing combustion instabilities in gas turbines. A control law formulated in this approach for fuel modulation demonstrated steady suppression of a high-frequency combustion instability (less than 500Hz) in a liquid-fuel combustion test rig under engine-realistic conditions. To make-up for severe transport-delays on control effect, the MSEK controller combines a wavelet -like Multi-Scale analysis and an Extended Kalman Observer to predict the thermo-acoustic states of combustion pressure perturbations. The commanded fuel modulation is composed of a damper action based on the predicted states, and a tones suppression action based on the Multi-Scale estimation of thermal excitations and other transient disturbances. The controller performs automatic adjustments of the gain and phase of these actions to minimize the Time-Scale Averaged Variances of the pressures inside the combustion zone and upstream of the injector. The successful demonstration of Active Combustion Control with this MSEK controller completed an important NASA milestone for the current research in advanced combustion technologies.

Le, Dzu K.↗

Multiple Damage Progression Paths in Model-Based Prognostics

Model-based prognostics approaches employ domain knowledge about a system, its components, and how they fail through the use of physics-based models. Component wear is driven by several different degradation phenomena, each resulting in their own damage progression path, overlapping to contribute to the overall degradation of the component. We develop a model-based prognostics methodology using particle filters, in which the problem of characterizing multiple damage progression paths is cast as a joint state-parameter estimation problem. The estimate is represented as a probability distribution, allowing the prediction of end of life and remaining useful life within a probabilistic framework that supports uncertainty management. We also develop a novel variance control mechanism that maintains an uncertainty bound around the hidden parameters to limit the amount of estimation uncertainty and, consequently, reduce prediction uncertainty. We construct a detailed physics-based model of a centrifugal pump, to which we apply our model-based prognostics algorithms. We illustrate the operation of the prognostic solution with a number of simulation-based experiments and demonstrate the performance of the chosen approach when multiple damage mechanisms are active

Daigle, Matthew↗

Prognostics for Ground Support Systems: Case Study on Pneumatic Valves

Prognostics technologies determine the health (or damage) state of a component or sub-system, and make end of life (EOL) and remaining useful life (RUL) predictions. Such information enables system operators to make informed maintenance decisions and streamline operational and mission-level activities. We develop a model-based prognostics methodology for pneumatic valves used in ground support equipment for cryogenic propellant loading operations. These valves are used to control the flow of propellant, so failures may have a significant impact on launch availability. Therefore, correctly predicting when valves will fail enables timely maintenance that avoids launch delays and aborts. The approach utilizes mathematical models describing the underlying physics of valve degradation, and, employing the particle filtering algorithm for joint state-parameter estimation, determines the health state of the valve and the rate of damage progression, from which EOL and RUL predictions are made. We develop a prototype user interface for valve prognostics, and demonstrate the prognostics approach using historical pneumatic valve data from the Space Shuttle refueling system.

Daigle, Matthew↗

An Uncertainty Quantification Framework for Prognostics and Condition-Based Monitoring

This paper presents a computational framework for uncertainty quantification in prognostics in the context of condition-based monitoring of aerospace systems. The different sources of uncertainty and the various uncertainty quantification activities in condition-based prognostics are outlined in detail, and it is demonstrated that the Bayesian subjective approach is suitable for interpreting uncertainty in online monitoring. A state-space model-based framework for prognostics, that can rigorously account for the various sources of uncertainty, is presented. Prognostics consists of two important steps. First, the state of the system is estimated using Bayesian tracking, and then, the future states of the system are predicted until failure, thereby computing the remaining useful life of the system. The proposed framework is illustrated using the power system of a planetary rover test-bed, which is being developed and studied at NASA Ames Research Center.

Health Monitoring↗

A Residuals Approach to Filtering, Smoothing and Identification for Static Distributed Systems

An approach for state estimation and identification of spatially distributed parameters embedded in static distributed (elliptic) system models is advanced. The method of maximum likelihood is used to find parameter values that maximize a likelihood functional for the system model, or equivalently, that minimize the negative logarithm of this functional. To find the minimum, a Newton-Raphson search is conducted that from an initial estimate generates a convergent sequence of parameter estimates. For simplicity, a Gauss-Markov approach is used to approximate the Hessian in terms of products of first derivatives. The gradient and approximate Hessian are computed by first arranging the negative log likelihood functional into a form based on the square root factorization of the predicted covariance of the measurement process. The resulting data processing approach, referred to here by the new term of predicted data covariance square root filtering, makes the gradient and approximate Hessian calculations very simple. A closely related set of state estimates is also produced by the maximum likelihood method: smoothed estimates that are optimal in a conditional mean sense and filtered estimates that emerge from the predicted data covariance square root filter.

Rodriguez, G.↗

Orbit Determination Accuracy Analysis of the Magnetospheric Multiscale Mission During Perigee Raise

The Goddard Space Flight Center (GSFC) Flight Dynamics Facility (FDF) will provide orbit determination and prediction support for the Magnetospheric Multiscale (MMS) mission during the missions commissioning period. The spacecraft will launch into a highly elliptical Earth orbit in 2015. Starting approximately four days after launch, a series of five large perigee-raising maneuvers will be executed near apogee on a nearly every-other-orbit cadence. This perigee-raise operations concept requires a high-accuracy estimate of the orbital state within one orbit following the maneuver for performance evaluation and a high-accuracy orbit prediction to correctly plan and execute the next maneuver in the sequence. During early mission design, a linear covariance analysis method was used to study orbit determination and prediction accuracy for this perigee-raising campaign. This paper provides a higher fidelity Monte Carlo analysis using the operational COTS extended Kalman filter implementation that was performed to validate the linear covariance analysis estimates and to better characterize orbit determination performance for actively maneuvering spacecraft in a highly elliptical orbit. The study finds that the COTS extended Kalman filter tool converges on accurate definitive orbit solutions quickly, but prediction accuracy through orbits with very low altitude perigees is degraded by the unpredictability of atmospheric density variation.

Determination↗

Orbit Determination Accuracy Analysis of the Magnetospheric Multiscale Mission During Perigee Raise

The Goddard Space Flight Center (GSFC) Flight Dynamics Facility (FDF) will provide orbit determination and prediction support for the Magnetospheric Multiscale (MMS) mission during the mission's commissioning period. The spacecraft will launch into a highly elliptical Earth orbit in 2015. Starting approximately four days after launch, a series of five large perigee-raising maneuvers will be executed near apogee on a nearly every-other-orbit cadence. This perigee-raise operations concept requires a high-accuracy estimate of the orbital state within one orbit following the maneuver for performance evaluation and a high-accuracy orbit prediction to correctly plan and execute the next maneuver in the sequence. During early mission design, a linear covariance analysis method was used to study orbit determination and prediction accuracy for this perigee-raising campaign. This paper provides a higher fidelity Monte Carlo analysis using the operational COTS extended Kalman filter implementation that was performed to validate the linear covariance analysis estimates and to better characterize orbit determination performance for actively maneuvering spacecraft in a highly elliptical orbit. The study finds that the COTS extended Kalman filter tool converges on accurate definitive orbit solutions quickly, but prediction accuracy through orbits with very low altitude perigees is degraded by the unpredictability of atmospheric density variation.

Multiscale↗

Climatological Structures of the GRIPS Models: Mean States and Forcing

The GCM-Reality Intercomparison Project for SPARC (GRIPS) is assessing and monitoring the performance of state-of-the-art general circulation models (GCMs). A wide variety of tasks have been initiated. These are designed to: (1) assess the ability of the GCMs to represent the current climatological structure of the troposphere and middle atmosphere,(2) to compare their response to imposed forcing anomalies, and (3) to estimate the certainty with which future climate perturbations can be predicted. This paper is concerned with assessments of the climatological states in the GCM simulations. Comparing the simulations with observational datasets reveals considerable discrepancies in the modelled fields. While it might be anticipated that certain types of biases in the model simulations might be related to the formulation of different aspects of the numerical package (dynamical schemes, cloud schemes, radiation transfer, inclusion of gravity wave drag), there is no clear relationship between these features. This paper attempts to draw a more comprehensive picture of the GCMs'performance than has previously been shown, by comparing the dominant forcing mechanisms in the models with observational estimates, and relating model deficiencies to the differences in the physical mechanisms in the GCMS.

Pawson, Steven↗

Mars Science Laboratory Navigation Results

The Mars Science Laboratory (MSL), carrying the Curiosity rover to Mars, was launched on November 26, 2011, from Cape Canaveral, Florida. The target for MSL was selected to be Gale Crater, near the equator of Mars, with an arrival date in early August 2012. The two main interplanetary navigation tasks for the mission were to deliver the spacecraft to an entry interface point that would allow the rover to safely reach the landing area, and to tell the spacecraft where it entered the atmosphere of Mars, so it could guide itself accurately to close proximity of the landing target. MSL used entry guidance as it slowed down from the entry speed to a speed low enough to allow for a successful parachute deployment, and this guidance allowed shrinking the landing ellipse to a 99% conservative estimate of 7 by 20 kilometers. Since there is no global positioning system in Mars, achieving this accuracy was predicated on flying a trajectory that closely matched the reference trajectory used to design the guidance algorithm, and on initializing the guidance system with an accurate Mars-relative entry state that could be used as the starting point to integrate the inertial measurement unit data during entry and descent. The pre-launch entry flight path angle (EFPA) delivery requirement was +/- 0.20 deg, but after launch a smaller threshold of +/- 0.05 deg was used as the criteria for late trajectory correction maneuver (TCM) decisions. The pre-launch requirement for entry state knowledge was 2.8 kilometers in position error and 2 meters per second in velocity error, but also smaller thresholds were defined after launch to evaluate entry state update opportunities. The biggest challenge for the navigation team was to accurately predict the trajectory of the spacecraft, so the estimates of the entry conditions could be stable, and late trajectory correction maneuvers or entry parameter updates could be waved off. As a matter of fact, the prediction accuracy was such that the last TCM performed was a small burn executed eight days before landing, and the entry state that was calculated just 36 hours after that TCM, and that was uploaded to the spacecraft the same day, did not need to be updated. The final EFPA was 0.013 deg shallower than the -15.5 deg target, and the on-board entry state was just 200 meters in position and 0.11 meters per second in velocity from the post-landing reconstructed entry state. Overall the entry delivery and knowledge requirements were fulfilled with a margin of more than 90% with respect to the pre-launch thresholds. This excellent accuracy contributed to a very successful and accurate entry, descent, and landing, and surface mission.

navigation↗

Equilibrium properties of a one-dimensional kinetic system.

One-dimensional systems of N = 500 and 250 particles in equilibrium are numerically simulated utilizing the method of molecular dynamics. Periodic boundary conditions are imposed. The classical two-body interaction potential is short range, repulsive and has a corresponding finite force. The equations of state are determined for densities both less and greater than one. Corresponding theoretical isochores are determined from models based on nearest-neighbor interactions and on a truncated virial expansion, and a comparison is made with the experimental isochores. Time independent radial distributions are constructed numerically and discussed. A change of state from a solidlike state to a fluid-gas state based on the penetrability of the particles is predicted. The transition temperatures are estimated from the radial distribution functions and the nearest-neighbor model. Self-diffusion is observed and the corresponding constants are determined from the velocity autocorrelation functions.

Williams, J. H.↗

Some New Mathematical Methods for Variational Objective Analysis

New and/or improved variational methods for simultaneously combining forecast, heterogeneous observational data, a priori climatology, and physics to obtain improved estimates of the initial state of the atmosphere for the purpose of numerical weather prediction are developed. Cross validated spline methods are applied to atmospheric data for the purpose of improved description and analysis of atmospheric phenomena such as the tropopause and frontal boundary surfaces.

Wahba, G.↗