Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “error estimation and analysis”

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 487 records · Page 27

Semimajor Axis Estimation Strategies

This paper extends previous analysis on the impact of sensing noise for the navigation and control aspects of formation flying spacecraft. We analyze the use of Carrier-phase Differential GPS (CDGPS) in relative navigation filters, with a particular focus on the filter correlation coefficient. This work was motivated by previous publications which suggested that a "good" navigation filter would have a strong correlation (i.e., coefficient near -1) to reduce the semimajor axis (SMA) error, and therefore, the overall fuel use. However, practical experience with CDGPS-based filters has shown this strong correlation seldom occurs (typical correlations approx. -0.1), even when the estimation accuracies are very good. We derive an analytic estimate of the filter correlation coefficient and demonstrate that, for the process and sensor noises levels expected with CDGPS, the expected value will be very low. It is also demonstrated that this correlation can be improved by increasing the time step of the discrete Kalman filter, but since the balance condition is not satisfied, the SMA error also increases. These observations are verified with several linear simulations. The combination of these simulations and analysis provide new insights on the crucial role of the process noise in determining the semimajor axis knowledge.

How, Jonathan P.↗

The dynamics of global positioning system orbits and the determination of precise ephemerides

The suggestion made on the basis of the analytical orbit perturbation theory that the errors in the ephemerides of the GPS satellites are due mostly to resonant effects that can be corrected by adjusting a few parameters in a empirical acceleration formula is tested using simulations and actual data analysis. Data from the Spring 1985 Experiment were used to calculate improved ephemerides, and these ephemerides were used in the estimation of the coordinates of GPS stations within the continental United States, previously positioned with VLBI. The results of this test support the idea that the errors are mostly of a resonant nature and can be corrected.

Colombo, Oscar L.↗

The GEOS Ozone Data Assimilation System: Specification of Error Statistics

A global three-dimensional ozone data assimilation system has been developed at the Data Assimilation Office of the NASA/Goddard Space Flight Center. The Total Ozone Mapping Spectrometer (TOMS) total ozone and the Solar Backscatter Ultraviolet (SBUV) or (SBUV/2) partial ozone profile observations are assimilated. The assimilation, into an off-line ozone transport model, is done using the global Physical-space Statistical Analysis Scheme (PSAS). This system became operational in December 1999. A detailed description of the statistical analysis scheme, and in particular, the forecast and observation error covariance models is given. A new global anisotropic horizontal forecast error correlation model accounts for a varying distribution of observations with latitude. Correlations are largest in the zonal direction in the tropics where data is sparse. Forecast error variance model is proportional to the ozone field. The forecast error covariance parameters were determined by maximum likelihood estimation. The error covariance models are validated using x squared statistics. The analyzed ozone fields in the winter 1992 are validated against independent observations from ozone sondes and HALOE. There is better than 10% agreement between mean Halogen Occultation Experiment (HALOE) and analysis fields between 70 and 0.2 hPa. The global root-mean-square (RMS) difference between TOMS observed and forecast values is less than 4%. The global RMS difference between SBUV observed and analyzed ozone between 50 and 3 hPa is less than 15%.

Stajner, Ivanka↗

An L-Band Brightness Temperature Disaggregation Method Using S-Band Radiometer Data for the Water Cycle Observation Mission (WCOM)

The Water Cycle Observation Mission (WCOM) will build upon previous L and C band passive microwave soil moisture satellite missions. WCOM will consist of a passive microwave synthetic aperture radiometer operating at L, S, and C bands. The WCOM requirements for passive soil moisture are to estimate soil moisture in the top 5 cm of soil layer with an error less than 0.04 cu m/cu m, at 15 km resolution and with a 3-day revisit. A new set of algorithms for these multi-frequency platforms will need to be developed for estimating the data products at the desired resolution. To accomplish this, a brightness temperature (TB) downscaling methodology is developed that uses passive S-band TB (30 km) to downscale L-band TB (50 km) and to estimate soil moisture at a 30 km resolution, based on the linear relationships between the passive signals of L-band and S-band. To test this downscaling method, analysis was performed using PALS data from the Soil Moisture Experiments in 2002 (SMEX02). For this study, 4 km L-band observations were downscaled to 800 m. The root mean square errors (RMSE) between the downscaled TBL at 800m with the observed TBL at 800m are 2.63 K and 1.60 K for H and V polarizations respectively. The results also showed that it was possible to use these disaggregated TB to estimate soil moisture to meet the mission requirement of 0.04 cu m/cu m. These results showed that we can obtain higher resolution soil moisture from L band passive TB with a high accuracy (<0.04 cu m/cu m) by using S-band information.

Disaggregation↗

A variable step incremental procedure

Description of a variable step incremental procedure for the solution of nonlinear equations in finite element structural analysis. The proposed procedure is effective in improving the accuracy of the basic incremental technique and in providing, in addition, an accurate estimate of the discretization error. The proposed approach is highly appropriate for solving problems for which the user has no a priori estimate of the step size to use.

Thomas, G. R.↗

Microwave radiometer measurement of water vapor path delay: Data reduction techniques

Data reduction techniques are developed to compensate for water vapor path delay, a limiting error source in geodetic measurements made with very long baseline interferometry and in radio ranging to spacecraft. It is shown that water vapor path delay is proportional to a linear combination of saturation-corrected sky brightness temperatures, measured on and off the water vapor line. The effects of emission from liquid water droplets in clouds as well as most of the oxygen emission are removed by the off-line channel. Sky brightness temperatures are saturation-corrected or 'linearized' using estimates of effective sky temperatures made from surface temperature. Tipping curves are used to remove instrumental error. Coefficients are found by two methods: from a regression analysis of measured brightness temperatures versus radiosonde measured delay, and from a regression analysis of theoretical brightness temperatures versus radiosonde measured delay. In each case the coefficients are adjusted for differing climatic conditions by measurements of surface temperature, pressure, and relative humidity. Regression solutions are constrained to remove liquid water contributions and to give the correct slope for radiometer versus radiosonde path delay.

Claflin, E. S.↗

A unified perspective on robot control - The energy Lyapunov function approach

A unified framework for the stability analysis of robot tracking control is presented. By using an energy-motivated Lyapunov function candidate, the closed-loop stability is shown for a large family of control laws sharing a common structure of proportional and derivative feedback and a model-based feedforward. The feedforward can be zero, partial or complete linearized dynamics, partial or complete nonlinear dynamics, or linearized or nonlinear dynamics with parameter adaptation. As result, the dichotomous approaches to the robot control problem based on the open-loop linearization and nonlinear Lyapunov analysis are both included in this treatment. Furthermore, quantitative estimates of the trade-offs between different schemes in terms of the tracking performance, steady state error, domain of convergence, realtime computation load and required a prior model information are derived.

Wen, John T.↗

Fragmented Land Cover Types and Estimation of Area with Course Spatial Resolution Imagery

Imagery of coarse resolution, such weather satellite imagery with 1 sq km pixels, is increasingly used to monitor dynamic and fragmented types of land surface types, such as scars from recent fires and ponds in wetlands. Accurate estimates of these land cover types at regional to global scales are required to assess the roles of fires and wetlands in global warming, yet difficult to compute when much of the area is accounted for by fragments about the same size as the pixels. In previous research, we found that size distribution of the fragments in several example scenes fit simple two-parameter models and related effects of coarse resolution to errors in area estimates based on pixel counts. We report on progress to develop accurate area estimates based on modelling the size distribution of the fragments, including analysis of size distributions on an expanded set of maps developed from digital imagery and a test of a procedure to correct for effects of coarse spatial resolution.

Hlavka, Chris↗

Regularization and Error Characterization of GRACE Mascons

We present a new global time-variable gravity mascon solution derived from Gravity Recovery and Climate Experiment (GRACE) Level 1B data. The new product from the NASA Goddard Space Flight Center (GSFC) results from a novel approach that combines an iterative solution strategy with geographical binning of inter-satellite range-acceleration residuals in the construction of time-dependent regularization matrices applied in the inversion of mascon parameters. This estimation strategy is intentionally conservative as it seeks to maximize the role of the GRACE measurements on the final solution while minimizing the influence of the regularization design process. We fully reprocess the Level 1B data in the presence of the final mascon solution to generate true post-fit inter-satellite residuals, which are utilized to confirm solution convergence and to validate the mascon noise uncertainties. We also present the mathematical case that regularized mascon solutions are biased, and that this bias, or leakage, must be combined with the estimated noise variance to accurately assess total mascon uncertainties. The estimated leakage errors are determined from the monthly resolution operators. We present a simple approach to compute the total uncertainty for both individual mascon and regional analysis of the GSFC mascon product, and validate the results in comparison with independent mascon solutions and calibrated Stokes uncertainties. Lastly, we present the new solution and uncertainties with global analyses of the mass trends and annual amplitudes, and compute updated trends for the global ocean, and the respective contributions of the Greenland Ice Sheet, Antarctic Ice Sheet, Gulf of Alaska, and terrestrial water storage. This analysis highlights the successful closure of the global mean sea level budget, that is, the sum of global ocean mass from the GSFC mascons and the steric component from Argo floats agrees well with the total determined from sea surface altimetry.

B. D. Loomis↗

Regularization and Error Characterization of GRACE Mascons

We present a new global time-variable gravity mascon solution derived from Gravity Recovery and Climate Experiment (GRACE) Level 1B data. The new product from the NASA Goddard Space Flight Center (GSFC) results from a novel approach that combines an iterative solution strategy with geographical binning of inter-satellite range-acceleration residuals in the construction of time-dependent regularization matrices applied in the inversion of mascon parameters. This estimation strategy is intentionally conservative as it seeks to maximize the role of the GRACE measurements on the final solution while minimizing the influence of the regularization design process. We fully reprocess the Level 1B data in the presence of the final mascon solution to generate true post-fit inter-satellite residuals, which are utilized to confirm solution convergence and to validate the mascon noise uncertainties. We also present the mathematical case that regularized mascon solutions are biased, and that this bias, or leakage, must be combined with the estimated noise variance to accurately assess total mascon uncertainties. The estimated leakage errors are determined from the monthly resolution operators. We present a simple approach to compute the total uncertainty for both individual mascon and regional analysis of the GSFC mascon product, and validate the results in comparison with independent mascon solutions and calibrated Stokes uncertainties. Lastly, we present the new solution and uncertainties with global analyses of the mass trends and annual amplitudes, and compute updated trends for the global ocean, and the respective contributions of the Greenland Ice Sheet, Antarctic Ice Sheet, Gulf of Alaska, and terrestrial water storage. This analysis highlights the successful closure of the global mean sea level budget, that is, the sum of global ocean mass from the GSFC mascons and the steric component from Argo floats agrees well with the total determined from sea surface altimetry.

Loomis, B. D.↗

Sketch-To-Solution: An Exploration of Viscous CFD with Automatic Grids

Numerical simulation of the Reynolds-averaged Navier–Stokes (RANS) equations has become a critical tool for the design of aerospace vehicles. However, the issues that affect the grid convergence of three dimensional RANS solutions are not completely understood, as documented in the AIAA Drag Prediction Workshop series. Grid adaption methods have the potential for increasing the automation and discretization error control of RANS solutions to impact the aerospace design and certification process. The realization of the CFD Vision 2030 Study includes automated management of errors and uncertainties of physics-based, predictive modeling that can set the stage for ensuring a vehicle is in compliance with a regulation or specification by using analysis without demonstration in flight test (i.e., certification or qualification by analysis). For example, the Cart3D inviscid analysis package has automated Cartesian cut-cell gridding with output-based error control. Fueled by recent advances in the fields of anisotropic grid adaptation, error estimation, and geometry modeling, a similar work flow is explored for viscous CFD simulations; where a CFD application engineer provides geometry, boundary conditions, and flow parameters, and the sketch-to-solution process yields a CFD simulation through automatic, error-based, grid adaptation.

Kleb, William L.↗

Determining Aircraft Moments of Inertia from Flight Test Data

Flight test maneuver design and dynamic modeling techniques were developed for determining aircraft moments of inertia from flight test data. Full nonlinear rigid-body rotational equations of motion were used in the analysis, with aerodynamic moment dependencies modeled by linear expansions in the aircraft states and controls. Aerodynamic parameters were estimated simultaneously with inertia parameters using equation-error modeling applied to flight test data from maneuvers designed specifically for this problem. The approach was demonstrated using a nonlinear F-16 simulation, then applied to a remotely-piloted subscale aircraft flight test. Errors in the aircraft moment of inertia parameters determined from simulated F-16 flight test maneuver data were less than 6 percent compared to the true values in the simulation. Flight test results for the subscale aircraft were within 6 percent of ground-test values obtained using the same aircraft.

moments of inertia↗

Determining Aircraft Moments of Inertia from Flight Test Data

Flight test maneuver design and dynamic modeling techniques were developed for determining aircraft moments of inertia from flight test data. Full nonlinear rigid-body rotational equations of motion were used in the analysis, with aerodynamic moment dependencies modeled by linear expansions in the aircraft states and controls. Aerodynamic parameters were estimated simultaneously with inertia parameters using equation-error modeling applied to flight test data from maneuvers designed specifically for this problem. The approach was demonstrated using a nonlinear F 16 simulation, then applied to a remotely piloted subscale aircraft flight test. Errors in the aircraft moment of inertia parameters determined from simulated F-16 flight test maneuver data were less than 6 percent compared to the true values in the simulation. Flight test results for the subscale aircraft were within 6 percent of ground-test values obtained using the same aircraft.

Aircraft Moments of Inertia↗

Determining Aircraft Moments of Inertia from Flight Test Data

Flight test maneuvers and dynamic modeling techniques were developed for determining aircraft moments of inertia from flight test data. Full nonlinear rigid-body rotational equations of motion were used in the analysis, with aerodynamic moment dependencies modeled by linear expansions in the aircraft states and controls. Aerodynamic parameters were estimated simultaneously with inertia parameters using equation-error modeling applied to flight test data from maneuvers designed specifically for this problem. The approach was demonstrated using a nonlinear F 16 simulation, then applied to a remotely-piloted subscale aircraft flight test. Errors in the aircraft moment of inertia parameters determined from simulated F-16 flight test data were less than 6 percent compared to the true values in the simulation. Flight test results for the subscale aircraft were within 6 percent of ground-test values obtained using the same aircraft.

Parameter estimation↗

Towards an Optimal Estimation Retrieval of Cirrus Cloud Optical and Microphysical Properties Using Hyperspectral Shortwave Instruments and A Fast Radiative Transfer Algorithm

Cirrus cloud retrieval products (here, cloud optical depth, effective particle size, and cloud top height) are important inputs into numerical weather and climate models. Uncertainties in such retrieval products, as a matter of course, propagate downstream, impacting model calculations. Improvements in high-quality global cirrus cloud optical and microphysical data products from satellite observations are needed to understand and reduce retrieval uncertainties. Hyperspectral shortwave instruments produce high-resolution and information-dense spectra, thus offering the opportunity to reduce uncertainties in retrieval products. We are in the process of developing a retrieval that uses the very fast Principal Component Radiative Transfer Model in the solar spectral region (PCRTM-Solar) in the forward model calculations. This retrieval will use measured reflectances from the NASA Earth Surface Mineral Dust Source Investigation (EMIT) and the forthcoming Climate Absolute Radiance and Refractivity Observatory Pathfinder (CLARREO-Pathfinder) instruments. In this manuscript we present progress towards a reference retrieval employing a widely used, verified, accurate, yet computationally slower radiative transfer modeling technique. The reference retrieval, while too slow for using the complete hyperspectral measurement, will allow us to study the behavior of retrieval products and help us verify results from our in-development fast retrieval. Both retrievals will use the optimal estimation retrieval framework. In this manuscript we present results from an uncertainty analysis considering three uncertainty sources for a cirrus cloud retrieval in the form of error covariance matrices: reflectance uncertainty due to water vapor, the reflectance uncertainty due to ice crystal scattering assumptions, and the instrument measurement uncertainty. Results show that the uncertainty due to habit selection is the largest, while that due to water vapor is at most 0.6% relative to channel reflectance. As a first step, the retrieval is being designed for single layer ice clouds over open ocean water.

Cirrus cloud↗

Dynamic Analyses of the Proposed Habitable Exoplanet Astrophysics Facility

The proposed Habitable Exoplanet (HabEx) astrophysics facility is one of four large such facilities being proposed to the 2020 decadal. It is a large telescope that is sensitive to ultraviolet, optical, and near-infrared photons. The proposed design’s overall length is on the order of 17.2 m and its maximum cross section is on the order of 5.25 X 5.25 m. The primary mirror is 4 m in diameter. A transient dynamic analysis was performed to estimate the order of magnitude of ring down time after moving the telescope and pointing at a new target for science planning purposes. Without uncertainty factors, results from a simple re-pointing maneuver indicate that primary to secondary mirror LOS errors are on the order of 10-4 pico-m after 5 minutes. Also, a frequency response analysis was performed to predict the impact of planned micro-thruster vibrations on required stability. Based on provided noise level associated with the micro-thrusters and loading assumptions and without uncertainty factors, the assessed vibrations do not impact predicted performance requirements.

Knight, J. Brent↗

Benefit of Modeling the Observation Error in a Data Assimilation Framework Using Vegetation Information Obtained From Passive Based Microwave Data

A primary operational goal of the United States Department of Agriculture (USDA) is to improve foreign market access for U.S. agricultural products. A large fraction of this crop condition assessment is based on satellite imagery and ground data analysis. The baseline soil moisture estimates that are currently used for this analysis are based on output from the modified Palmer two-layer soil moisture model, updated to assimilate near-real time observations derived from the Soil Moisture Ocean Salinity (SMOS) satellite. The current data assimilation system is based on a 1-D Ensemble Kalman Filter approach, where the observation error is modeled as a function of vegetation density. This allows for offsetting errors in the soil moisture retrievals. The observation error is currently adjusted using Normalized Difference Vegetation Index (NDVI) climatology. In this paper we explore the possibility of utilizing microwave-based vegetation optical depth instead.

Vegetation↗