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 793 records · Page 44

Cross-correlation of point series using a new method

Traditional methods of cross-correlation of two time series do not apply to point time series. Here, a new method, devised specifically for point series, utilizes a correlation measure that is based in the rms difference (or, alternatively, the median absolute difference) between nearest neightbors in overlapped segments of the two series. Error estimates for the observed locations of the points, as well as a systematic shift of one series with respect to the other to accommodate a constant, but unknown, lead or lag, are easily incorporated into the analysis using Monte Carlo techniques. A methodological restriction adopted here is that one series be treated as a template series against which the other, called the target series, is cross-correlated. To estimate a significance level for the correlation measure, the adopted alternative (null) hypothesis is that the target series arises from a homogeneous Poisson process. The new method is applied to cross-correlating the times of the greatest geomagnetic storms with the times of maximum in the undecennial solar activity cycle.

Strothers, Richard B.↗

Assimilation of Quality Controlled AIRS Temperature Profiles using the NCEP GFS

We have previously conducted a number of data assimilation experiments using AIRS Version-5 quality controlled temperature profiles as a step toward finding an optimum balance of spatial coverage and sounding accuracy with regard to improving forecast skill. The data assimilation and forecast system we used was the Goddard Earth Observing System Model , Version-5 (GEOS-5) Data Assimilation System (DAS), which represents a combination of the NASA GEOS-5 forecast model with the National Centers for Environmental Prediction (NCEP) operational Grid Point Statistical Interpolation (GSI) global analysis scheme. All analyses and forecasts were run at a 0.5deg x 0.625deg spatial resolution. Data assimilation experiments were conducted in four different seasons, each in a different year. Three different sets of data assimilation experiments were run during each time period: Control; AIRS T(p); and AIRS Radiance. In the "Control" analysis, all the data used operationally by NCEP was assimilated, but no AIRS data was assimilated. Radiances from the Aqua AMSU-A instrument were also assimilated operationally by NCEP and are included in the "Control". The AIRS Radiance assimilation adds AIRS observed radiance observations for a select set of channels to the data set being assimilated, as done operationally by NCEP. In the AIRS T(p) assimilation, all information used in the Control was assimilated as well as Quality Controlled AIRS Version-5 temperature profiles, i.e., AIRS T(p) information was substituted for AIRS radiance information. The AIRS Version-5 temperature profiles were presented to the GSI analysis as rawinsonde profiles, assimilated down to a case-by-case appropriate pressure level p(sub best) determined using the Quality Control procedure. Version-5 also determines case-by-case, level-by-level error estimates of the temperature profiles, which were used as the uncertainty of each temperature measurement. These experiments using GEOS-5 have shown that forecasts resulting from analyses using the AIRS T(p) assimilation system were superior to those from the Radiance assimilation system, both with regard to global 7 day forecast skill and also the ability to predict storm tracks and intensity.

Susskind, Joel↗

Camera modeling, centroiding performance, and geometric camera calibration on ASTERIA

The Arcsecond Space Telescope Enabling Research in Astrophysics (ASTERIA) is a 10-kg, 6U CubeSat in low-Earth orbit that was able to achieve subarcsecond pointing stability and repeatability. To date, this is the best pointing on a spacecraft of its size. This paper will analyze various aspects of the performance of its key piece of hardware—the payload. First,a model of the optics and imager, which is used to simulate stellar images, will be presented. The imager parameters used in this model were derived from simple ground measurements. Next, a centroiding algorithm is provided and used on the simulated images to predict centroiding performance. These results will be shown to match on-orbit telemetry of centroiding performance, validating the modeling approach. This paper will then describe an approach for and results of a geometric camera calibration algorithm to estimate the focal length,distortion,and alignment parameters. The modeling, analyses, and results presented in this paper provide key information that can be used in a time domain pointing simulation or a frequency-domain pointing error analysis.

Smith, Matthew W.↗

Stratospheric Assimilation of Chemical Tracer Observations Using a Kalman Filter: Chi-Square Validated Results and Analysis of Variance and Correlation Dynamics - Pt. 2

A Kalman filter system designed for the assimilation of limb-sounding observations of stratospheric chemical tracers, which has four tunable covariance parameters, was developed in Part I (Menard et al. 1998) The assimilation results of CH4 observations from the Cryogenic Limb Array Etalon Sounder instrument (CLAES) and the Halogen Observation Experiment instrument (HALOE) on board of the Upper Atmosphere Research Satellite are described in this paper. A robust (chi)(sup 2) criterion, which provides a statistical validation of the forecast and observational error covariances, was used to estimate the tunable variance parameters of the system. In particular, an estimate of the model error variance was obtained. The effect of model error on the forecast error variance became critical after only three days of assimilation of CLAES observations, although it took 14 days of forecast to double the initial error variance. We further found that the model error due to numerical discretization as arising in the standard Kalman filter algorithm, is comparable in size to the physical model error due to wind and transport modeling errors together. Separate assimilations of CLAES and HALOE observations were compared to validate the state estimate away from the observed locations. A wave-breaking event that took place several thousands of kilometers away from the HALOE observation locations was well captured by the Kalman filter due to highly anisotropic forecast error correlations. The forecast error correlation in the assimilation of the CLAES observations was found to have a structure similar to that in pure forecast mode except for smaller length scales. Finally, we have conducted an analysis of the variance and correlation dynamics to determine their relative importance in chemical tracer assimilation problems. Results show that the optimality of a tracer assimilation system depends, for the most part, on having flow-dependent error correlation rather than on evolving the error variance.

Menard, Richard↗

Towards robust autonomous impedance spectroscopy analysis: A calibrated hierarchical Bayesian approach for electrochemical impedance spectroscopy (EIS) inversion

Distribution-based analyses, such as the distribution of relaxation times (DRT) and the distribution of diffusion times (DDT), present model-free alternatives to equivalent circuit modeling for analysis of electrochemical impedance spectroscopy (EIS) data. However, reconstructing such distributions from noisy impedance data is an ill-posed problem that must be solved with specialized inversion algorithms, requiring careful control and tuning. Furthermore, most inversion algorithms developed to date can only solve problems of limited complexity. Herein, we present a new hierarchical Bayesian method for EIS inversion, leveraging efficient algorithms for optimization and Hamiltonian Monte Carlo (HMC) sampling to solve models of arbitrary complexity. We overcome the challenge of ad-hoc parameter tuning by encoding intrinsic characteristics of the DRT and DDT into flexible prior distributions and “pre-calibrating” the model to simulated data. This approach is versatile, highly robust to noise, and provides quantitative estimates of both the error structure of the data and the uncertainty in the recovered distributions. The model is validated with simulated data to demonstrate accurate recovery of the DRT and the DDT. The method also shows promise for simultaneous recovery of multiple distributions, raising the intriguing possibility of semi-autonomous EIS analysis and ad-hoc model construction. Finally, the practical utility of the method is illustrated with experimental data. Throughout, we draw comparisons to several recently published EIS inversion methodologies.

36 MATERIALS SCIENCE↗

Magnetospheric Multiscale (MMS) Mission Commissioning Phase Orbit Determination Error Analysis

The Magnetospheric MultiScale (MMS) mission commissioning phase starts in a 185 km altitude x 12 Earth radii (RE) injection orbit and lasts until the Phase 1 mission orbits and orientation to the Earth-Sun li ne are achieved. During a limited time period in the early part of co mmissioning, five maneuvers are performed to raise the perigee radius to 1.2 R E, with a maneuver every other apogee. The current baseline is for the Goddard Space Flight Center Flight Dynamics Facility to p rovide MMS orbit determination support during the early commissioning phase using all available two-way range and Doppler tracking from bo th the Deep Space Network and Space Network. This paper summarizes th e results from a linear covariance analysis to determine the type and amount of tracking data required to accurately estimate the spacecraf t state, plan each perigee raising maneuver, and support thruster cal ibration during this phase. The primary focus of this study is the na vigation accuracy required to plan the first and the final perigee ra ising maneuvers. Absolute and relative position and velocity error hi stories are generated for all cases and summarized in terms of the ma ximum root-sum-square consider and measurement noise error contributi ons over the definitive and predictive arcs and at discrete times inc luding the maneuver planning and execution times. Details of the meth odology, orbital characteristics, maneuver timeline, error models, and error sensitivities are provided.

Chung, Lauren R.↗

Estimating Contrail Climate Effects from Satellite Data

An automated contrail detection algorithm (CDA) is developed to exploit six of the infrared channels on the 1-km MODerate-resolution Imaging Spectroradiometer (MODIS) on the Terra and Aqua satellites. The CDA is refined and balanced using visual error analysis. It is applied to MODIS data taken by Terra and Aqua over the United States during 2006 and 2008. The results are consistent with flight track data, but differ markedly from earlier analyses. Contrail coverage is a factor of 4 less than other retrievals and the retrieved contrail optical depths and radiative forcing are smaller by approx.30%. The discrepancies appear to be due to the inability to detect wider, older contrails that comprise a significant amount of the contrail coverage. An example of applying the algorithm to MODIS data over the entire Northern Hemisphere is also presented. Overestimates of contrail coverage are apparent in some tropical regions. Methods for improving the algorithm are discussed and are to be implemented before analyzing large amounts of Northern Hemisphere data. The results should be valuable for guiding and validating climate models seeking to account for aviation effects on climate.

Minnis, Patrick↗

Common View Time Transfer Using Worldwide GPS and DMA Monitor Stations

Analysis of the on-orbit Navstar clocks and the Global Positioning System (GPS) monitor station reference clocks is performed by the Naval Research Laboratory using both broadcast and postprocessed precise ephemerides. The precise ephemerides are produced by the Defense Mapping Agency (DMA) for each of the GPS space vehicles from pseudo-range measurements collected at five GPS and at five DMA monitor stations spaced around the world. Recently, DMA established an additional site co-located with the US Naval Observatory precise time site. The time reference for the new DMA site is the DoD Master Clock. Now, for the first time, it is possible to transfer time every 15 minutes via common view from the DoD Master Clock to the 11 GPS and DMA monitor stations. The estimated precision of a single common-view time transfer measurement taken over a 15-minute interval was between 1.4 and 2.7 nanoseconds. Using the measurements from all Navstar space vehicles in common view during the 15-minute interval, typically 3-7 space vehicles, improved the estimate of the precision to between 0.65 and 1.13 nanoseconds. The mean phase error obtained from closure of the time transfer around the world using the 11 monitor stations and the 25 space vehicle clocks over a period of 4 months had a magnitude of 31 picoseconds. Analysis of the low noise time transfer from the DoD Master Clock to each of the monitor stations yields not only the bias in the time of the reference clock, but also focuses attention on structure in the behaviour of the reference clock not previously seen. Furthermore, the time transfer provides a a uniformly sampled database of 15-minute measurements that make possible, for the first time, the direct and exhaustive computation of the frequency stability of the monitor station reference clocks. To lend perspective to the analysis, a summary is given of the discontinuities in phase and frequency that occurred in the reference clock at the Master Control Station during the period covered by the analysis.

Reid, Wilson G.↗

Estimating Soil Thermal Inertia Profiles From the Passive Equilibration of a Temperature Probe

Knowledge of the distribution of soil thermal properties is important for understanding subsurface hydrological and biogeochemical processes. This study describes and evaluates quick thermal profiling (QTP), a new measurement technique aimed at providing rapid, depth-resolved measurements of soil thermal inertia at numerous locations across the landscape. A cylindrical probe with temperature sensors at multiple depths is quickly inserted into the ground, and soil thermal inertia is estimated from how quickly the probe temperature equilibrates with the soil. To this end, a finite volume heat transfer model is used to generate temperature equilibration time series across combinations of controlling factors, and a gridded search inversion approach is applied to infer soil thermal inertia. Field tests in the Arctic indicate that QTP measurements have a minimum uncertainty of 0.14 J m −2 K −1 s −1/2 and covary with dual-probe heat pulse thermal analyzer measurements (concordance correlation coefficient = 0.56) with a root-mean-square error of 0.40 J m −2 K −1 s −1/2 . Besides demonstrating the value of QTP for estimating thermal inertia, this study identifies various sources of measurement uncertainty, particularly probe-soil contact resistance and frictional heating. Further, analysis of soil samples indicates that thermal inertia can be used to estimate thermal conductivity and dry bulk density in the studied area, although such inferences are highly site-specific. Overall, the QTP method holds promise to generate thermal inertia data products and to complement other characterization approaches for advancing understanding of soil properties across far more locations than is currently possible.

Lamb, J. R. [Lawrence Berkeley National Laboratory↗

Enhanced shear stabilization of turbulence in NSTX

In studying a particular non-stationary NSTX L-mode plasma, we observed unexpectedly high levels of flux—first with quasilinear (TGLF) modeling and subsequently with nonlinear gyrokinetic simulations. Upon more detailed analysis, a novel confinement regime was discovered in which a modest increase in E x B shear (beyond baseline experimental estimates) rapidly reduced turbulent transport to levels consistent with power balance. This modest increase is plausible given the errors inherent to the estimation of shearing rates, and the added complexity of the non-stationary (time-dependent) power balance. Remarkably, an additional small increase in shear yields the familiar ion-neoclassical transport level with what appears to be the onset of high-k electron transport only. Although analyses using the TGLF-SAT2 model successfully capture numerous parametric dependencies of this plasma, TGLF does not reproduce the rapid E x B stabilization seen in CGYRO. We believe the results presented should help to better characterize the nonlinear physics of spherical tokamak confinement regimes, provide useful ST datasets for reduced model development, and motivate more accurate experimental diagnosis of E x B shearing rates.

Atomic and molecular collisions↗

Long and short arc altitude determination for GEOS-C

The accuracy with which the GEOS-C altitude may be estimated over long (7 day) and short (40 minute) orbital arcs is investigated. Over the long are excellent agreement was attained between a simulation of the orbit determination process and a covariance analysis. Both approaches yielded RMS altitude errors of about 1.5 meters over the Caribbean calibration area and approximately 7.5 meters overall. The geopotential was identified as the largest error source. For the short arc, the covariance analysis revealed that the propagated altitude error is linearly dependent upon station survey component errors which are also the largest source of altitude errors. An Appendix contains the mathematics of covariance analysis as applied to orbit determination.

Koch, D.↗

Stochastic process approximation for recursive estimation with guaranteed bound on the error covariance

An approach, is proposed for the design of approximate, fixed order, discrete time realizations of stochastic processes from the output covariance over a finite time interval, was proposed. No restrictive assumptions are imposed on the process; it can be nonstationary and lead to a high dimension realization. Classes of fixed order models are defined, having the joint covariance matrix of the combined vector of the outputs in the interval of definition greater or equal than the process covariance; (the difference matrix is nonnegative definite). The design is achieved by minimizing, in one of those classes, a measure of the approximation between the model and the process evaluated by the trace of the difference of the respective covariance matrices. Models belonging to these classes have the notable property that, under the same measurement system and estimator structure, the output estimation error covariance matrix computed on the model is an upper bound of the corresponding covariance on the real process. An application of the approach is illustrated by the modeling of random meteorological wind profiles from the statistical analysis of historical data.

Menga, G.↗

Application of in-flight thrust determination uncertainty

A numerical example is given of a previously proposed methodology for the evaluation of in-flight thrust measurement uncertainty, using data extracted from a performance report comparing two different missile prototypes under a variety of flight conditions. Attention is given to the data for the AGM-68B Air Launched Cruise Missile, which is powered by the F107 dual-spool, mixed flow turbofan engine. Assessments are made of the definition of the measurement process, instrumentation error estimation, the propagation of errors to thrust calculation, mathematical model errors, the in-flight thrust error component, and correction to standard conditions. It is concluded that in-flight thrust measurement uncertainty limits can be evaluated from measurement system error analysis results and test data for the missile evaluation process presently described.

Adams, G. R.↗

Statistical magnetic anomalies from satellite measurements for geologic analysis

The errors of numerically averaging satellite magnetic anomaly data for geologic analysis are investigated using orbital anomaly simulations of crustal magnetic sources by Gauss-Legendre quadrature integration. These simulations suggest that numerical averaging errors constitute small and relatively minor contributions to the total error-budget of higher orbital estimates (approx. 400 km), whereas for lower orbital estimates the error of averaging may increase substantially. Least-squares collocation is also investigated as an alternative to numerical averaging and found to produce substantially more accurate anomaly estimates as the elevation of prediction is decreased towards the crustal sources.

Goyal, H. K.↗

Modeling non-gravitational forces acting on TOPEX/POSEIDON: The early days

TOPEX/POSEIDON is a satellite mission that will use altimetry to make precise measurements of sea-level. The principal goal is to measure sea-level with unprecedented accuracy such that small-amplitude, basin wide sea-level changes caused by large-scale ocean circulation can be detected. To reach this goal, the sensor system and orbit must measure sea-level with decimeter accuracy. This requires that the radial component of the orbit be known to the decimeter level. Orbital errors are dominated by mismodelled gravitational and non-gravitational forces. This paper presents our analysis of non-gravitational forces acting on the satellite during the early days of the mission. Studies were conducted by comparing direct estimates of these forces with observed perturbations in the mean orbital elements. The results show that the satellite is experiencing an unexpected along-track acceleration. Hypotheses range from out-gassing to thruster leaks and drag forces to radiative force. Currently, these issues have not been resolved; however, the evidence suggests that out-gassing was dominant during the first weeks of the mission and that thermal imbalances persist.

Christensen, E. J.↗

Short arc orbit determination for altimeter calibration and validation on TOPEX/POSEIDON

TOPEX/POSEIDON (T/P) is a joint mission of United States' National Aeronautics and Space Administration (NASA) and French Centre National d'Etudes Spatiales (CNES) design launched August 10, 1992. It carries two radar altimeters which alternately share a common antenna. There are two project designated verification sites, a NASA site off the coast at Pt. Conception, CA and a CNES site near Lampedusa Island in the Mediterranean Sea. Altimeter calibration and validation for T/P is performed over these highly instrumented sites by comparing the spacecraft's altimeter radar range to computed range based on in situ measurements which include the estimated orbit position. This paper presents selected results of orbit determination over each of these sites to support altimeter verification. A short arc orbit determination technique is used to estimate a locally accurate position determination of T/P from less than one revolution of satellite laser ranging (SLR) data. This technique is relatively insensitive to gravitational and non-gravitational force modeling errors and is demonstrated by covariance analysis and by comparison to orbits determined from longer arcs of data and other tracking data types, such as Doppler Orbitography and Radiopositioning Integrated by Satellite (DORIS) and Global Positioning System Demonstration Receiver (GPSDR) data.

Williams, B. G.↗

The Global Precipitation Climatology Project (GPCP) Combined Precipitation Dataset

The Global Precipitation Climatology Project (GPCP) has released the GPCP Version 1 Combined Precipitation Data Set, a global, monthly precipitation dataset covering the period July 1987 through December 1995. The primary product in the dataset is a merged analysis incorporating precipitation estimates from low-orbit-satellite microwave data, geosynchronous-orbit -satellite infrared data, and rain gauge observations. The dataset also contains the individual input fields, a combination of the microwave and infrared satellite estimates, and error estimates for each field. The data are provided on 2.5 deg x 2.5 deg latitude-longitude global grids. Preliminary analyses show general agreement with prior studies of global precipitation and extends prior studies of El Nino-Southern Oscillation precipitation patterns. At the regional scale there are systematic differences with standard climatologies.

Huffman, George J.↗

Improved Flight Test Procedures for Flutter Clearance

Flight flutter testing is an integral part of flight envelope clearance. This paper discusses advancements in several areas that are being investigated to improve efficiency and safety of flight test programs. Results are presented from recent flight testing of the F/A-18 Systems Research Aircraft. A wingtip excitation system was used to generate aeroelastic response data. This system worked well for many flight conditions but still displayed some anomalies. Wavelet processing is used to analyze the flight data. Filtered transfer functions are generated that greatly improve system identification. A flutter margin is formulated that accounts for errors between a model and flight data. Worst-case flutter margins are computed to demonstrate the flutter boundary may lie closer to the flight envelope than previously estimated. This paper concludes with developments for a distributed flight analysis environment and on-line health monitoring.

Lind, Rick C.↗