Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 1,099 records · Page 61

Estimation of rate constants of elementary processes - A review of the state of the art.

'Thermochemical kinetics,' the codification and extrapolation of empirical observations, as applied to certain elementary reactions of importance to combustion studies, is described. This approach allows the critical scrutiny of experimental data in areas where sufficient previous data exist, while, at the same time, illuminating those key areas where more experimentation is crucial. It is shown that combination of transition-state theory with an understanding of the molecular basis of entropy puts fairly rigid constraints on the values of the Arrhenius A-factor for most reactions. This, in turn, means that the activation energy is often the key datum that is missing, and that such data can be obtained with some degree of confidence, even from measurements of rate constants at only one temperature. In complex mechanisms, it is often possible to distinguish among alternate pathways and pinpoint key processes.

Golden, D. M.↗

An equation of state for oxygen and nitrogen

Preliminary equations of state are presented for oxygen and nitrogen which provide accurate representations of the available P-density-T data for both fluids. The equation for nitrogen is applicable for temperatures from 70 K to 1300 K at pressures to 10,000 atmospheres, and the equation for oxygen for temperatures from 70 K to 323 K at pressures to 350 atmospheres. Deviations of calculated densities from representative experimental data are included. A volume-explicit equation of state for oxygen to be used in estimating density values in the range of applicability of the equation of state is also presented.

Stewart, R. B.↗

Orbit determination accuracies using satellite-to-satellite tracking

The uncertainty in relay satellite sate is a significant error source which cannot be ignored in the reduction of satellite-to-satellite tracking data. Based on simulations and real data reductions, it is numerically impractical to use simultaneous unconstrained solutions to determine both relay and user satellite epoch states. A Bayesian or least squares estimation technique with an a priori procedure is presented which permits the adjustment of relay satellite epoch state in the reduction of satellite-to-satellite tracking data without the numerical difficulties introduced by an ill-conditioned normal matrix.

Vonbun, F. O.↗

Mars Phoenix EDL Trajectory and Atmosphere Reconstruction Using NewSTEP

This document describes the trajectory and atmosphere reconstruction of the Mars Phoenix Entry, Descent, and Landing using the New Statistical Trajectory Estimation Program. The approach utilizes a Kalman filter to blend inertial measurement unit data with initial conditions and radar altimetry to obtain the inertial trajectory of the entry vehicle. The nominal aerodynamic database is then used in combination with the sensed accelerations to obtain estimates of the atmosphere-relative state. The reconstructed atmosphere pro le is then blended with pre-flight models to construct an estimate of the as-flown atmosphere.

Karlgaard, Christopher D.↗

Estimating Ground-Level PM(sub 2.5) Concentrations in the Southeastern United States Using MAIAC AOD Retrievals and a Two-Stage Model

Previous studies showed that fine particulate matter (PM(sub 2.5), particles smaller than 2.5 micrometers in aerodynamic diameter) is associated with various health outcomes. Ground in situ measurements of PM(sub 2.5) concentrations are considered to be the gold standard, but are time-consuming and costly. Satellite-retrieved aerosol optical depth (AOD) products have the potential to supplement the ground monitoring networks to provide spatiotemporally-resolved PM(sub 2.5) exposure estimates. However, the coarse resolutions (e.g., 10 km) of the satellite AOD products used in previous studies make it very difficult to estimate urban-scale PM(sub 2.5) characteristics that are crucial to population-based PM(sub 2.5) health effects research. In this paper, a new aerosol product with 1 km spatial resolution derived by the Multi-Angle Implementation of Atmospheric Correction (MAIAC) algorithm was examined using a two-stage spatial statistical model with meteorological fields (e.g., wind speed) and land use parameters (e.g., forest cover, road length, elevation, and point emissions) as ancillary variables to estimate daily mean PM(sub 2.5) concentrations. The study area is the southeastern U.S., and data for 2003 were collected from various sources. A cross validation approach was implemented for model validation. We obtained R(sup 2) of 0.83, mean prediction error (MPE) of 1.89 micrograms/cu m, and square root of the mean squared prediction errors (RMSPE) of 2.73 micrograms/cu m in model fitting, and R(sup 2) of 0.67, MPE of 2.54 micrograms/cu m, and RMSPE of 3.88 micrograms/cu m in cross validation. Both model fitting and cross validation indicate a good fit between the dependent variable and predictor variables. The results showed that 1 km spatial resolution MAIAC AOD can be used to estimate PM(sub 2.5) concentrations.

aerosol optical depth↗

Study of a homotopy continuation method for early orbit determination with the Tracking and Data Relay Satellite System (TDRSS)

A recent mathematical technique for solving systems of equations is applied in a very general way to the orbit determination problem. The study of this technique, the homotopy continuation method, was motivated by the possible need to perform early orbit determination with the Tracking and Data Relay Satellite System (TDRSS), using range and Doppler tracking alone. Basically, a set of six tracking observations is continuously transformed from a set with known solution to the given set of observations with unknown solutions, and the corresponding orbit state vector is followed from the a priori estimate to the solutions. A numerical algorithm for following the state vector is developed and described in detail. Numerical examples using both real and simulated TDRSS tracking are given. A prototype early orbit determination algorithm for possible use in TDRSS orbit operations was extensively tested, and the results are described. Preliminary studies of two extensions of the method are discussed: generalization to a least-squares formulation and generalization to an exhaustive global method.

Smith, R. L.↗

A steady-state Kalman filter for assimilating data from a single polar orbiting satellite

A steady-state scheme for data assimilation in the context of a single, short period (relative to a day), sun-synchronous, polar-orbiting satellite is examined. If the satellite takes observations continuously, the gains, which are the weights for blending observations and predictions together, are steady in time. For a linear system forced by random noise, the optimal steady-state gains (Wiener gains) are equivalent to those of a Kalman filter. Computing the Kalman gains increases the computational cost of the model by a large factor, but computing the Wiener gains does not. The latter are computed by iteration using prior estimates of the gains to assimilate simulated observations of one run of the model, termed 'truth' into another run termed 'prediction'. At each stage, the prediction errors form the basis for the next estimate of the gains. Steady state is achieved after three or four iterations. Further simplification is achieved by making the gains depend on longitudinal distance from the observation point, not on absolute longitude. For a single-layer primitive equation model, the scheme works well even if only the mass field is observed but not the velocity field. Although the scheme was developed for Mars Observer, it should be applicable to data retrieved from Earth atmosphere satellites, for example, UARS.

Banfield, Don↗

An Approach for the Long-Term 30-m Land Surface Snow-Free Albedo Retrieval from Historic Landsat Surface Reflectance and MODIS-based A Priori Anisotropy Knowledge

Land surface albedo has been recognized by the Global Terrestrial Observing System (GTOS) as an essential climate variable crucial for accurate modeling and monitoring of the Earth's radiative budget. While global climate studies can leverage albedo datasets from MODIS, VIIRS, and other coarse-resolution sensors, many applications in heterogeneous environments can benefit from higher-resolution albedo products derived from Landsat. We previously developed a "MODIS-concurrent" approach for the 30-meter albedo estimation which relied on combining post-2000 Landsat data with MODIS Bidirectional Reflectance Distribution Function (BRDF) information. Here we present a "pre-MODIS era" approach to extend 30-m surface albedo generation in time back to the 1980s, through an a priori anisotropy Look-Up Table (LUT) built up from the high quality MCD43A BRDF estimates over representative homogenous regions. Each entry in the LUT reflects a unique combination of land cover, seasonality, terrain information, disturbance age and type, and Landsat optical spectral bands. An initial conceptual LUT was created for the Pacific Northwest (PNW) of the United States and provides BRDF shapes estimated from MODIS observations for undisturbed and disturbed surface types (including recovery trajectories of burned areas and non-fire disturbances). By accepting the assumption of a generally invariant BRDF shape for similar land surface structures as a priori information, spectral white-sky and black-sky albedos are derived through albedo-to-nadir reflectance ratios as a bridge between the Landsat and MODIS scale. A further narrow-to-broadband conversion based on radiative transfer simulations is adopted to produce broadband albedos at visible, near infrared, and shortwave regimes.We evaluate the accuracy of resultant Landsat albedo using available field measurements at forested AmeriFlux stations in the PNW region, and examine the consistency of the surface albedo generated by this approach respectively with that from the "concurrent" approach and the coincident MODIS operational surface albedo products. Using the tower measurements as reference, the derived Landsat 30-m snow-free shortwave broadband albedo yields an absolute accuracy of 0.02 with a root mean square error less than 0.016 and a bias of no more than 0.007. A further cross-comparison over individual scenes shows that the retrieved white sky shortwave albedo from the "pre-MODIS era" LUT approach is highly consistent (R(exp 2) = 0.988, the scene-averaged low RMSE = 0.009 and bias = −0.005) with that generated by the earlier "concurrent" approach. The Landsat albedo also exhibits more detailed landscape texture and a wider dynamic range of albedo values than the coincident 500-m MODIS operational products (MCD43A3), especially in the heterogeneous regions. Collectively, the "pre-MODIS" LUT and "concurrent" approaches provide a practical way to retrieve long-term Landsat albedo from the historic Landsat archives as far back as the 1980s, as well as the current Landsat-8 mission, and thus support investigations into the evolution of the albedo of terrestrial biomes at fine resolution.

Albedo algorithm↗

The estimation of the constituent densities of the upper atmosphere by means of a recursive filtering algorithm.

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.

Mcgarty, T. P.↗

Assessment of ICESat-2 Sea Ice Surface Classification with Sentinel-2 Imagery: Implications for Freeboard and New Estimates of Lead and Floe Geometry

NASA's Ice, Cloud, and Land Elevation Satellite-2 (ICESat-2) mission launched in September 2018 and is now providing high-resolution surface elevation profiling across the entire globe, including the sea ice cover of the Arctic and Southern Oceans. For sea ice applications, successfully discriminating returns between sea ice and open water is key for accurately determining freeboard (the extension of sea ice above local sea level) and new information regarding the geometry of sea ice floes and leads. We take advantage of near-coincident optical imagery obtained from the European Space Agency (ESA) Sentinel-2 (S-2) satellite over the Western Weddell Sea of the Southern Ocean in March 2019 and the Lincoln Sea of the Arctic Ocean in May 2019 to evaluate the surface classification scheme in the ICESat-2 ATL07 and ATL10 sea ice products. We find a high level of agreement between the ATL07 (specular) lead classification and visible leads in the S-2 imagery in these two coincident images across all six ICESat-2 beams, increasing our confidence in the freeboard products and deriving new estimates of the sea ice state. The S-2 overlays provide additional, albeit limited, evidence of the misclassification of dark leads due to clouds. Dark leads are no longer used to derive sea surface and thus freeboard as of the third release (r003) of the ICESat-2 sea ice products. We show estimates of lead fraction and more preliminary estimates of chord length (a proxy for floe size) using two metrics for classifying sea surface (lead) segments across both the Arctic and Southern Ocean for the first winter season of data collection.

A. A. Petty↗

On the sea-state bias of the Geosat altimeter

The sea-state bias in a satellite altimeter's range measurement is caused by the influence of ocean waves on the radar return pulse; it results in an estimate of sea level that is too low according to some function of the wave height. This bias is here estimated for Geosat by correlating collinear differences of altimetric sea-surface heights with collinear differences of significant wave heights (H1/3). Corrections for satellite orbit error are estimated simultaneously with the sea-state bias. Based on twenty 17-day repeat cycles of the Geosat Exact Repeat Mission, the solution for the sea-state bias is 2.6 + or - 0.2 percent of H1/3. The least-squares residuals, however, show a correlation with wind speed U, so the traditional model of the bias has been supplemented with a second term: H1/3 + alpha-2H1/3U. This second term produces a small, but statistically significant, reduction in variance of the residuals. Both systematic and random errors in H1/3 and U tend to bias the estimates of alpha-1 and alpha-2, which complicates comparisons of the results with ground-based measurements of the sea-state bias.

Ray, Richard D.↗

Ensemble Forecasting/Assimilation/Emissions Estimation with WRF-Chem/DART: Accomplishments, Lessons Learned, and Future Plans

Over the past ten years, we have been conducting research on regional ensemble atmospheric composition forecasting/data assimilation/emissions estimation with WRF-Chem/DART. WRF- Chem/DART integrates the Weather Research and Forecasting model (WRF) with online chemistry (WRF-Chem) into the Data Assimilation Research Testbed (DART). DART is an ensemble data assimilation system based on the ensemble adjustment Kalman filter (EAKF) with adaptive inflation, localization (physical and state space), and an optional non-Gaussian formulation of the EAKF. DART includes assimilation of meteorological and limited chemical observations. WRF-Chem/DART extends DART to include assimilation of: MOPITT CO; IASI CO and O3; MODIS AOD; OMI O3, NO2, and SO2; TROPOMI CO, O3, NO2, and SO2, TES CO, CO2 (research mode), O3, NH3, and CH4 (research mode); CrIS CO, O3, NH3, CH4 (research mode), and PAN; SCIAMACHY NO2; GOME2a NO2; MLS O3 and HNO3; and proxy TEMPO O3, and NO2 satellite retrievals as raw retrievals or as ‘compact phase space retrievals’ (CPSRs) for profile retrievals. WRF-Chem/DART also assimilates in situ atmospheric composition measurements and uses the ‘state augmentation method’ for emissions estimation. In our presentation, we will provide an overview of WRF-Chem/DART: • applications and results; • lessons learned from: (i) independent versus joint assimilation; (ii) total/partial column versus profile retrieval assimilation; (iii) joint in situ and retrieval assimilation; (iv) assimilation at grid resolutions ranging from 100 km to 4 km; (v) dynamic emissions estimation; and • future work related to wildfire emissions estimation and intercomparison of CMAQ (online)/JEDI and CMAQ (online and offline)/DART.

WRF-Chem/DART↗

Extraction of Lateral-Directional Stability and Control Derivatives for the Basic F-18 Aircraft at High Angles of Attack

The results of parameter identification to determine the lateral-directional stability and control derivatives of an F-18 research aircraft in its basic hardware and software configuration are presented. The derivatives are estimated from dynamic flight data using a specialized identification program developed at NASA Dryden Flight Research Center. The formulation uses the linearized aircraft equations of motions in their continuous/discrete form and a maximum likelihood estimator that accounts for both state and measurement noise. State noise is used to model the uncommanded forcing function caused by unsteady aerodynamics, such as separated and vortical flows, over the aircraft. The derivatives are plotted as functions of angle of attack between 3 deg and 47 deg and compared with wind-tunnel predictions. The quality of the derivative estimates obtained by parameter identification is somewhat degraded because the maneuvers were flown with the aircraft's control augmentation system engaged, which introduced relatively high correlations between the control variables and response variables as a result of control motions from the feedback control system.

Iliff, Kenneth W.↗

Development, test and evaluation of a computerized procedure for using Landsat data to estimate spring small grains acreage

A number of methods which can provide information concerning crop acreages on the basis of a utilization of multispectral scanner (MSS) data require for their implementation a comparatively large amount of labor. The present investigation is concerned with a project designed to improve the efficiency of analysis through increased automation. The Caesar technique was developed to realize this objective. The processability rates of the Caesar procedure versus the historical state-of-the-art proportion estimation procedures were determined in an experiment. Attention is given to the study site, the aggregation technology, the results of the aggregation test, and questions of error characterization. It is found that the Caesar procedure, which has been developed for the spring small grains region of North America, is highly efficient and provides accurate results.

Mohler, R. R. J.↗

Prognostics Applied to Electric Propulsion UAV

Health management plays an important role in operations of UAV. If there is equipment malfunction on critical components, safe operation of the UAV might possibly be compromised. A technology with particular promise in this arena is equipment prognostics. This technology provides a state assessment of the health of components of interest and, if a degraded state has been found, it estimates how long it will take before the equipment will reach a failure threshold, conditional on assumptions about future operating conditions and future environmental conditions. This chapter explores the technical underpinnings of how to perform prognostics and shows an implementation on the propulsion of an electric UAV. A particle filter is shown as the method of choice in performing state assessment and predicting future degradation. The method is then applied to the batteries that provide power to the propeller motors. An accurate run-time battery life prediction algorithm is of critical importance to ensure the safe operation of the vehicle if one wants to maximize in-air time. Current reliability based techniques turn out to be insufficient to manage the use of such batteries where loads vary frequently in uncertain environments.

Batteries↗

Comparison of Kalman filter and optimal smoother estimates of spacecraft attitude

Given a valid system model and adequate observability, a Kalman filter will converge toward the true system state with error statistics given by the estimated error covariance matrix. The errors generally do not continue to decrease. Rather, a balance is reached between the gain of information from new measurements and the loss of information during propagation. The errors can be further reduced, however, by a second pass through the data with an optimal smoother. This algorithm obtains the optimally weighted average of forward and backward propagating Kalman filters. It roughly halves the error covariance by including future as well as past measurements in each estimate. This paper investigates whether such benefits actually accrue in the application of an optimal smoother to spacecraft attitude determination. Tests are performed both with actual spacecraft data from the Extreme Ultraviolet Explorer (EUVE) and with simulated data for which the true state vector and noise statistics are exactly known.

Sedlak, J.↗

The Shock State of Itokawa Sample

One of the fundamental aspects of any astromaterial is its shock history, since this factor elucidates critical historical events, and also because shock metamorphism can alter primary mineralogical and petrographic features, and reset chronologies [1]. Failure to take shock history into proper account during characterization can result in seriously incorrect conclusions being drawn. Thus the Hayabusa Preliminary Examination Team (HASPET) made shock stage determination of the Itokawa samples a primary goal [2]. However, we faced several difficulties in this particular research. The shock state of ordinary chondrite materials is generally determined by simple optical petrographic observation of standard thin sections. The Itokawa samples available to the analysis team were mounted into plastic blocks, were polished on only one side, and were of non-standard and greatly varying thickness, all of which significantly complicated petrographic analysis but did not prevent it. We made an additional estimation of the sample shock state by a new technique for this analysis - electron back-scattered diffraction (EBSD) in addition to standard petrographic techniques. We are also investigating the crystallinity of Itokawa olivine by Synchrotron X-ray diffraction (SXRD).

Zolensky, Michael↗

Picosecond excite-and-probe absorption measurement of the intra-2E(g)E(3/2)-state vibrational relaxation time in Ti(3+):Al2O3

The Ti(3+)-doped Al2O3 has been recently demonstrated to be a tunable solid-state laser system with Ti(3+) as the laser-active ion. In this paper, the kinetics of vibrational transitions in the 2E(g)E(3/2) electronic state of Ti(3+):Al2O3a (crucial for characterizing new host materials for the Ti ion) was investigated. A 527-nm 5-ps pulse was used to excite a band of higher vibrational levels of the 2E(g)E(3/2) state, and the subsequent growth of population in the zero vibrational level and lower vibrational levels was monitored by a 3.9-micron picosecond probe pulse. The time evolution curve in the excited 2E(g)E(3/2) state at room temperature was found to be characterized by a sharp rise followed by a long decay, the long lifetime decay reflecting the depopulation of the zero and the lower vibrational levels of the 2E(g)E(3/2) state via radiative transitions. An upper limit of 3.5 ps was estimated for intra-2E(g)E(3/2)-state vibrational relaxation time.

Gayen, S. K.↗