Search NASA⌕ Search

SEARCH · Search NASA

Results for “error 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 397 records · Page 22

Validation Practices for Satellite Soil Moisture retrievals: What are (the) errors?

This paper presents a community effort to develop good practice guidelines for the validation of global coarse-scale satellite soil moisture products. We provide theoretical background, a review of state-of-the-art methodologies for estimating errors in soil moisture data sets, practical recommendations on data pre-processing and presentation of statistical results, and a recommended validation protocol that is supplemented with an example validation exercise focused on microwave-based surface soil moisture products. We conclude by identifying research gaps that should be addressed in the near future.

A. Gruber↗

Measurement updating using the U-D factorization

A new mechanization of the Kalman updating algorithm based on a U-D factorization of the estimate error covariance is introduced. Efficient and stable updating recursions are developed for the unit upper triangular factor U and the diagonal factor D, treating only the parameter estimation problem. Properties of the factorization update performed here include efficient one point at a time processing that requires little more computation than the optimal but numerically unstable conventional Kalman measurement update algorithm, and stability that compares with the square root filter.

Bierman, G. J.↗

Application of the operator spline technique to nonlinear estimation and control of moving elastic systems

A bilinear model of the vibrational dynamics of a deformable maneuvering body is described. Estimates of the deformation state are generated through a low dimensional operator spline interpolator of bilinear systems combined with a feedback linearized based observer. Upper bounds on error estimates are also generated through the operator spline, and potential application to shaping control purposes is highlighted.

Karray, Fakhreddine↗

An Astronomical Test of CCD Photometric Precision

This article considers a posteriori error estimation of specified functionals for first-order systems of conservation laws discretized using the discontinuous Galerkin (DG) finite element method. Using duality techniques. we derive exact error representation formulas for both linear and nonlinear functionals given an associated bilinear or nonlinear variational form. Weighted residual approximations of the exact error representation formula are then proposed and numerically evaluated for Ringleb flow, an exact solution of the 2-D Euler equations.

Koch, David↗

Anisotropic Goal-Based Mesh Adaptation Metric Clarification and Development

Adaptive unstructured mesh techniques have a limited, but growing impact on production analysis workflows to control discretization error for reliable simulation results. Multiple independent implementations of flow solvers, anisotropic metric construction methods, and anisotropic mesh adaptation mechanics have matured. Goal-based metrics target estimated error in output functions, such as lift and drag, through the guidance of an adjoint solution. A unification of goal-based anisotropic metrics is presented for steady viscous flows, which is an active area of research. These goal-based metrics drive robust and efficient anisotropic mesh adaptation for the calculation of output functions. The super-convergent functional output error behavior of stabilized finite-element methods is exploited without a formal proof, and evidence of super-convergence is shown in numerical experiments. Mesh adapted drag and lift outputs for two simple bodies in compressible viscous flow show convergence of error to less than a single drag count. Asymptotic behavior established for relatively coarse meshes shows the efficiency of this goal-based metric when compared to solution interpolation error control and expert-guided meshing. Anisotropic mesh adaptation techniques are applied to a transport aircraft in a high-lift configuration where variation between approaches decreases with mesh refinement, but asymptotic behavior is not observed with available resources.

goal-based↗

Determination of intercontinental baselines and Earth orientation using VLBI

A series of experiments was conducted during the last decade to explore the capability of very long baseline interferometry (VLBI) to measure the crustal and rotational motions of the Earth with accuracies at the centimeter level. The observing stations are those of NASA's Deep Space Network in California, Spain and Australia. A multiparameter fit to the observed values of delay and delay rate yields radio source positions, polar motion, universal time, the precession constant, baseline vectors, and solid Earth tides. Source positions are obtained with formal errors of the order of 0''.01. UT1-UTC and polar motion are determined at 49 epochs, with formal error estimates for the more recent data of 0.5 msec for UT1-UTC and 2 to 6 mas for polar motion. Intercontinental baseline lengths are determined with formal errors of 5 to 10 cm. The Love numbers and Earth tide phase lag agree with the commonly accepted values.

Sovers, O. J.↗

Asteroid orbital error analysis: Theory and application

We present a rigorous Bayesian theory for asteroid orbital error estimation in which the probability density of the orbital elements is derived from the noise statistics of the observations. For Gaussian noise in a linearized approximation the probability density is also Gaussian, and the errors of the orbital elements at a given epoch are fully described by the covariance matrix. The law of error propagation can then be applied to calculate past and future positional uncertainty ellipsoids (Cappellari et al. 1976, Yeomans et al. 1987, Whipple et al. 1991). To our knowledge, this is the first time a Bayesian approach has been formulated for orbital element estimation. In contrast to the classical Fisherian school of statistics, the Bayesian school allows a priori information to be formally present in the final estimation. However, Bayesian estimation does give the same results as Fisherian estimation when no priori information is assumed (Lehtinen 1988, and reference therein).

Muinonen, K.↗

Approximate effect of parameter pseudonoise intensity on rate of convergence for EKF parameter estimators

When using parameter estimation methods based on extended Kalman filter (EKF) theory, it is common practice to assume that the unknown parameter values behave like a random process, such as a random walk, in order to guarantee their identifiability by the filter. The present work is the result of an ongoing effort to quantitatively describe the effect that the assumption of a fictitious noise (called pseudonoise) driving the unknown parameter values has on the parameter estimate convergence rate in filter-based parameter estimators. The initial approach is to examine a first-order system described by one state variable with one parameter to be estimated. The intent is to derive analytical results for this simple system that might offer insight into the effect of the pseudonoise assumption for more complex systems. Such results would make it possible to predict the estimator error convergence behavior as a function of the assumed pseudonoise intensity, and this leads to the natural application of the results to the design of filter-based parameter estimators. The results obtained show that the analytical description of the convergence behavior is very difficult.

Hill, Bryon K.↗

Altitude Registration of Limb-Scattered Radiation

One of the largest constraints to the retrieval of accurate ozone profiles from UV backscatter limb sounding sensors is altitude registration. Two methods, the Rayleigh scattering attitude sensing (RSAS) and absolute radiance residual method (ARRM), are able to determine altitude registration to the accuracy necessary for long-term ozone monitoring. The methods compare model calculations of radiances to measured radiances and are independent of onboard tracking devices. RSAS determines absolute altitude errors, but, because the method is susceptible to aerosol interference, it is limited to latitudes and time periods with minimal aerosol contamination. ARRM, a new technique introduced in this paper, can be applied across all seasons and altitudes. However, it is only appropriate for relative altitude error estimates. The application of RSAS to Limb Profiler (LP) measurements from the Ozone Mapping and Profiler Suite (OMPS) on board the Suomi NPP (SNPP) satellite indicates tangent height (TH) errors greater than 1 km with an absolute accuracy of +/-200 m. Results using ARRM indicate a approx. 300 to 400m intra-orbital TH change varying seasonally +/-100 m, likely due to either errors in the spacecraft pointing or in the geopotential height (GPH) data that we use in our analysis. ARRM shows a change of approx. 200m over 5 years with a relative accuracy (a long-term accuracy) of 100m outside the polar regions.

Ozone profiles↗

Convergence Analysis of Turbulent Flow Solutions

Data from the "Turbulence Modeling Resource" website for turbulent flow over an NACA-0012 airfoil is analyzed to determine the convergence behavior of three second-order CFD (Computational Fluid Dynamics) codes: CFL3D (Computational Fluids Lab 3 Dimensional flow solver), FUN3D (Fully Unstructured Navier-stokes flow solver), and TAU (German Aerospace Center (DLR) 2 dimensional code for unstructured hybrid grids solving the Reynolds-Averaged Navier-Stokes equations or the Euler equations). The convergence of both integrated properties and pointwise data are examined. Several different methods for estimating errors and computing convergence rates are compared. A high-order extension to the Richardson extrapolation is developed that improves the accuracy of the mesh limit values and provides a quantitative estimate of the threshold of the asymptotic regime. The coefficient of total drag exhibits second-order convergence for all three codes, and convergence is monotone over a sequence of 7 grids. Other force coefficients are not so well behaved. The convergence rates of the viscous component of drag on the three nest grids ranges from 3:0 for CFL3D to 1:0 for FUN3D. The three codes are converging to similar but not identical solutions. The largest differences between the codes are in the coefficient of lift for which the difference between CFL3D and FUN3D is greater than 10 (sup minus 4). The best agreement occurs in the viscous component of drag, which is the only force component for which all three codes are converging towards each other at a rate of second-order. The agreement between the two unstructured grid codes is good with all properties except lift converging towards common values at a rate of second-order. No one code was universally better than the other. The TAU code has the lowest error in total drag, FUN3D has the lowest error in lift, and CFL3D has the lowest error in the viscous component of drag.

Atkins, Harold L.↗

Covariance operator estimation via adaptive thresholding

This paper studies sparse covariance operator estimation for nonstationary processes with sharply varying marginal variance and small correlation lengthscale. We introduce a covariance operator estimator that adaptively thresholds the sample covariance function using an estimate of the variance component. Building on recent results from empirical process theory, we derive an operator norm bound on the estimation error in terms of the sparsity level of the covariance and the expected supremum of a normalized process. Furthermore, our theory and numerical simulations demonstrate the advantage of adaptive threshold estimators over universal threshold and sample covariance estimators in nonstationary settings.

Al-Ghattas, Omar [University of Chicago, IL (Unite↗

TRMM Field Campaigns: Objectives and Status Report

The Tropical Rainfall Measuring Mission (TRMM) satellite has been sending valuable data since launch in November 1997. Some of the key goals of the joint NASA (US) and NASDA (Japan) mission are: (1) to estimate the four-dimensional diabatic heating in the tropical and subtropical atmosphere, (2) understand the role of latent heating in driving tropical and extratropical circulations, (3) obtain monthly area-averaged estimates of rainfall over the data-sparse oceans, and (4) estimate the relative contribution of convective and stratiform precipitation over different regions during different seasons. The overarching scientific objective is to understand and improve estimates of rainfall and latent heating profiles throughout the global tropics. This requires observations for fundamental understanding of cloud dynamics and microphysics, as well as for validation, testing assumptions and error estimates of cloud-resolving models, forward radiative transfer models, algorithms used to estimate rainfall statistics and vertical structure of precipitation from surface-based radar, and from satellites. Field experiments designed to contribute to this understanding have been conducted in Texas and the South China Sea in spring of 1998, Florida in summer of 1998, and interior Brazil in (boreal) winter 1999. In summer 1999, a major oceanic campaign will be based at Kwajalein Atoll. Some early results will be highlighted, noting some significant contrasts between oceanic and continental convective systems.

Zipser, Edward I.↗

Hubble Space Telescope Magnetometer and Two-Gyro Control Law Design, Implementation, and On-Orbit Performance

For f i h years, the science mission of the Hubble Space Telescope (HST) required using at least three of six rate gyros for attitude control. In the past, HST has mitigated gyro hardware failures by replacement of the failed units through Space Shuttle Servicing Missions. Following the tragic loss of Space Shuttle Columbia on STS-107, the desire to have a safe haven for astronauts during missions has resulted in the cancellation of all planned maxu14 missions to HST. While a robotic servicing mission is being currently being planned, controlling with alternate sensors to replace failed gyros can extend the HST Science mission until the robotic mission can be performed and extend science at HST s end of life. A two-gym control law has been designed and implemented using magnetometers (Magnetic Sensing System - MSS), fixed head star trackers (FHSTs), and Fine Guidance Sensors (FGSs) to control vehicle rate about the missing gyro axis. The three aforementioned sensors are used in succession to reduce HST boresight jitter to less than 7 milli-arcseconds rms and attitude error to less than 10 milli-arcseconds prior to science imaging. The MSS and 2-Gyro (M2G) control law is used for large angle maneuvers and attitude control during earth occultation of FHSTs and FGSs. The Tracker and 2-Gyro (T2G) control law dampens M2G rates and corrects the majority of attitude error in preparation for guide star acquisition with the FGSs. The Fine Guidance Sensor and 2-Gyro (F2G) control law d a m p T2G rates and controls HST attitude during science imaging. This paper describes the M2G control law. Details of M2G algorithms are presented, including computation of the HST 3-axis attitude error estimate, design of the M2G control law, SISO hear stability analyses, and restrictions on operations to maintain the h d t h and safety requirement of a 10degree maximum attitude error. Results of simulations performed in HSTSIM, a high-fidelity non-linear time domain simulation, are presented to predict HST on-orbit performance in attitude hold and maneuver modes. Simulation results are compared to on-orbit data from M2G flight tests performed in November and December 2004 and February 2005. Flight telemetry, using a currently available third gyro, shows that HST attitude error with the new M2G control law is maintained below the 10-degree requirement, and attitude errors are under 2 degrees for 95% of the time.

Wirzburger, John H.↗

Measurement of D-T Neutron Source Flux by SelfNeutron Activation Analysis with a LaBr3 Detector

Conclusions: • Decay of 78Br produced from the 79Br(n,2n) reaction in LaBr3 is present in the spectrum plotted in the time domain. • The neutron production rate reconstructed from preliminary LaBr3 estimate is smaller than the generator output. • Generator characterization may become more routine if a simple and robust method is made available for applications. Ongoing work: • Implementing bootstrapping to improve error estimates. • Examining threshold and edge effects of LaBr3

SMITH, MURRY↗

Crack Turning and Arrest Mechanisms for Integral Structure

In the course of several years of research efforts to predict crack turning and flapping in aircraft fuselage structures and other problems related to crack turning, the 2nd order maximum tangential stress theory has been identified as the theory most capable of predicting the observed test results. This theory requires knowledge of a material specific characteristic length, and also a computation of the stress intensity factors and the T-stress, or second order term in the asymptotic stress field in the vicinity of the crack tip. A characteristic length, r(sub c), is proposed for ductile materials pertaining to the onset of plastic instability, as opposed to the void spacing theories espoused by previous investigators. For the plane stress case, an approximate estimate of r(sub c), is obtained from the asymptotic field for strain hardening materials given by Hutchinson, Rice and Rosengren (HRR). A previous study using of high order finite element methods to calculate T-stresses by contour integrals resulted in extremely high accuracy values obtained for selected test specimen geometries, and a theoretical error estimation parameter was defined. In the present study, it is shown that a large portion of the error in finite element computations of both K and T are systematic, and can be corrected after the initial solution if the finite element implementation utilizes a similar crack tip discretization scheme for all problems. This scheme is applied for two-dimensional problems to a both a p-version finite element code, showing that sufficiently accurate values of both K(sub I) and T can be obtained with fairly low order elements if correction is used. T-stress correction coefficients are also developed for the singular crack tip rosette utilized in the adaptive mesh finite element code FRANC2D, and shown to reduce the error in the computed T-stress significantly. Stress intensity factor correction was not attempted for FRANC2D because it employs a highly accurate quarter-point scheme to obtain stress intensity factors.

Pettit, Richard↗

Range validation using Kalman filter techniques

Range pseudo-residuals may be improved to the level required for data validation by a measurement updating process which utilizes Bierman's adaptation of the Kalman filter measurement updating algorithms together with process noise compensation to account for model errors. This algorithm involves combining the currently available range predictions and measurements to produce an updated range residual measurement whose accuracy is constrained by the range data quality and by the estimated error in the prediction. The algorithm is compact and fast, and is thus suitable for on-line applications in network control or at the station.

Madrid, G. A.↗

Five-Year Wilkinson Microwave Anisotropy Probe (WMAP1) Observations: Galactic Foreground Emission

We present a new estimate of foreground emission in the WMAP data, using a Markov chain Monte Carlo (MCMC) method. The new technique delivers maps of each foreground component for a variety of foreground models, error estimates of the uncertainty of each foreground component, and provides an overall goodness-of-fit measurement. The resulting foreground maps are in broad agreement with those from previous techniques used both within the collaboration and by other authors. We find that for WMAP data, a simple model with power-law synchrotron, free-free, and thermal dust components fits 90% of the sky with a reduced X(sup 2) (sub v) of 1.14. However, the model does not work well inside the Galactic plane. The addition of either synchrotron steepening or a modified spinning dust model improves the fit. This component may account for up to 14% of the total flux at Ka-band (33 GHz). We find no evidence for foreground contamination of the CMB temperature map in the 85% of the sky used for cosmological analysis.

Gold, B.↗

A passive microwave technique for estimating rainfall and vertical structure information from space. Part 2: Applications to SSM/I data

A multi channel physical approach for retrieving rainfall and its vertical structure from Special Sensor Microwave/Imager (SSM/I) observations is examined. While a companion paper was devoted exclusively to the description of the algorithm, its strengths, and its limitations, the main focus of this paper is to report on the results, applicability, and expected accuraciesfrom this algorithm. Some examples are given that compare retrieved results with ground-based radar data from different geographical regions to illustrate the performance and utility of the algorithm under distinct rainfall conditions. More quantitative validation is accomplished using two months of radar data from Darwin, Australia, and the radar network over Japan. Instantaneous comparisons at Darwin indicate that root-mean-square errors for 1.25 deg areas over water are 0.09 mm/h compared to the mean rainfall value of 0.224 mm/h while the correlation exceeds 0.9. Similar results are obtained over the Japanese validation site with rms errors of 0.615 mm/h compared to the mean of 0.0880 mm/h and a correlation of 0.9. Results are less encouraging over land with root-mean-square errors somewhat larger than the mean rain rates and correlations of only 0.71 and 0.62 for Darwin and Japan, respectively. These validation studies are further used in combination with the theoretical treatment of expected accuracies developed in the companion paper to define error estimates on a broader scale than individual radar sites from which the errors may be analyzed. Comparisons with simpler techniques that are based on either emission or scattering measurements are used to illustrate the fact that the current algorithm, while better correlated with the emission methods over water, cannot be reduced to either of these simpler methods.

Kummerow, Christian↗