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 883 records · Page 49

The inverse problem of constructing a gravimetric geoid

Computation of a single geoidal height from gravity acceleration data formally requires that the latter be known everywhere on the earth. A computational procedure based on linear inverse theory for estimating geoidal heights from incomplete sets of data is presented. The same scheme can be used to estimate gravity accelerations from altimetry-derived geoids. The systematic error owing to lack of data and the choice of a particular inverse operator is described by using resolution functions and their spherical harmonic expansions. An rms value of this error is also estimated by assuming a spectrum for the unknown geoid. The influence of the size of the data region, the spacing between data, the filtering applied to the data, and the model weighting function chosen are all quantified in a spherical geometry. The examples presented show that when low degree spherical harmonic coefficients are available - from satellite orbit analysis - a band-passed version of the geoid can be constructed from local gravity data, even with a relatively restricted data set.

Zlotnicki, V.↗

Future possibilities in objective analysis and data assimilation for atmospheric dynamics

The role that dynamics plays in estimating the state of the atmosphere from incomplete and noisy data is reviewed. Objective analysis represents an attempt at relying mostly on the data and minimizing the role of dynamics in the estimation. Data assimilation tries to balance properly the roles of dynamical and observational information. Sequential estimation is presented as the proper framework for understanding this balance, and the Kalman filter as the ideal, optimal procedure for data assimilation. The optimal filter computes forecast error covariances of a given atmospheric model exactly, and hence data assimilation should be closely connected with predictability studies. This connection is described, and consequences drawn for currently active areas of the atmospheric and related sciences, namely, mesoscale meteorology, long range forecasting, and upper ocean dynamics. Possibilities offered by judicious data assimilation in understanding barotropic adjustment, a phenomenon that appears to play a crucial role in atmospheric behavior on the scale of weeks to months, and hence in long range forecasting are addressed.

Ghil, M.↗

Color rendering of color camera data

Conditions under which a computational procedure can be applied to arbitrary camera sensors to permit estimation of the human photoreceptor response are considered. The adopted procedures recover the effective surface reflectance at the time of measurement, and the reflectance estimates depend not only on the surface, but upon the viewing geometry. The present method for color rendering assumes that the observer's state of adaptation at the time of viewing the original and the rendered images is the same. The analysis aids in specifying which classes of surfaces are required to be accurately rendered, and for which surfaces some error can be tolerated.

Wandell, Brian A.↗

Kalman estimation for SEDS measurements

The first on-orbit experiment of the Small Expendable Deployer System (SEDS) for tethered satellites will collect telemetry data for tethered length, rate of deployment, and tether tension. The post-flight analysis will use this data to reconstruct the deployment history and determine payload position and tether shape. Two Kalman estimator algorithms were written, and output using simulated measurement data was compared. Both estimators exhibited the same estimated state histories, indicating that numerical instability in the traditional algorithm was not the cause of filter divergence. Estimation of acceleration biases was added, which reduced the error but did not correct the divergence. An add-a-bead estimator that adds lumped masses as the tether is deployed was written, which provides a state model that matches the BEADSIM simulation providing the true measurements and states. This twenty-one bead estimator produced state histories similar to those of the two-bead estimator, indicating that the filter divergence was not caused by a reduced-order model. The noise models used to date are relatively simple and may be the source of estimator divergence. The investigation of colored noise models, cross-correlated measurement and process covariances, and noise-adaptive filter techniques is recommended.

Carrington, Connie K.↗

Elevation Change of the Southern Greenland Ice Sheet from Satellite Radar Altimeter Data

Long-term changes in the thickness of the polar ice sheets are important indicators of climate change. Understanding the contributions to the global water mass balance from the accumulation or ablation of grounded ice in Greenland and Antarctica is considered crucial for determining the source of the about 2 mm/yr sea-level rise in the last century. Though the Antarctic ice sheet is much larger than its northern counterpart, the Greenland ice sheet is more likely to undergo dramatic changes in response to a warming trend. This can be attributed to the warmer Greenland climate, as well as a potential for amplification of a global warming trend in the polar regions of the Northern Hemisphere. In collaboration with Drs. Curt Davis and Craig Kluever of the University of Missouri, we are using data from satellite radar altimeters to measure changes in the elevation of the Southern Greenland ice sheet from 1978 to the present. Difficulties with systematic altimeter measurement errors, particularly in intersatellite comparisons, beset earlier studies of the Greenland ice sheet thickness. We use altimeter data collected contemporaneously over the global ocean to establish a reference for correcting ice-sheet data. In addition, the waveform data from the ice-sheet radar returns are reprocessed to better determine the range from the satellite to the ice surface. At JPL, we are focusing our efforts principally on the reduction of orbit errors and range biases in the measurement systems on the various altimeter missions. Our approach emphasizes global characterization and reduction of the long-period orbit errors and range biases using altimeter data from NASA's Ocean Pathfinder program. Along-track sea-height residuals are sequentially filtered and backwards smoothed, and the radial orbit errors are modeled as sinusoids with a wavelength equal to one revolution of the satellite. The amplitudes of the sinusoids are treated as exponentially-correlated noise processes with a time-constant of six days. Measurement errors (e.g., altimeter range bias) are simultaneously recovered as constant parameters. The corrections derived from the global ocean analysis are then applied over the Greenland ice sheet. The orbit error and measurement bias corrections for different missions are developed in a single framework to enable robust linkage of ice-sheet measurements from 1978 to the present. In 1998, we completed our re-evaluation of the 1978 Seasat and 1985-1989 Geosat Exact Repeat Mission data. The estimates of ice thickness over Southern Greenland (south of 72N and above 2000 m) from 1978 to 1988 show large regional variations (+/-18 cm/yr), but yield an overall rate of +1.5 +/- 0.5 cm/yr (one standard error). Accounting for systematic errors, the estimate may not be significantly different from the null growth rate. The average elevation change from 1978 to 1988 is too small to assess whether the Greenland ice sheet is undergoing a long-term change.

Haines, Bruce J.↗

An Airborne Scanning Laser Altimetry Survey of Long Valley, California

Between 28 September and 7 October 1995, we conducted an airborne laser altimetry experiment over the Long Valley caldera, California, in which each of two scanning laser altimeters (dubbed SLICER and RASCAL) were flown in a NASA T-39 jet aircraft. Operating concurrently were a Global Positioning System (GPS) guidance system and dual frequency receivers for precise navigation and post-flight calculation or the airplane trajectory relative to a ground station, and an inertial navigation system (INS) for attitude determination. Reduction of raw laser ranges requires merging the differential kinematic GPS aircraft trajectory and the INS data with the laser data, and determination of the atmospheric delay. Data geolocation consists of obtaining the centre location and the mean elevation within each footprint in a geodetic coordinate system. The elevation of Crowley Lake is recovered to an accuracy of approximately 3 cm or better from 3 km above ground level and crossover analysis indicates that the elevation estimates are consistent from pass to pass. We test our geolocation procedures by comparing laser-derived elevations with those determined in situ for recognizable ground features. A comparison of laser and GPS-derived positions shows that the horizontal accuracy is better than the diameter of the footprint and vertical accuracy is within the error inherent in the range measurement. A comparison of SLICER elevation data with digital elevation models (DEMs) of the region shows that the DEM data provides surface topography to within stated accuracy limits. Although research continues to utilize the full potential of laser altimetry data, our results constitute a successful demonstration that the technique may be used to perform geodetic monitoring of surface topographic changes.

Hofton, M. A.↗

Precise Clock Solutions Using Carrier Phase from GPS Receivers in the International GPS Service

As one of its activities as an Analysis Center in the International GPS Service (IGS), the Jet Propulsion Laboratory (JPL) uses data from a globally distributed network of geodetic-quality GPS receivers to estimate precise clock solutions, relative to a chosen reference, for both the GPS satellites and GPS receiver internal clocks, every day. The GPS constellation and ground network provide geometrical strength resulting in formal errors of about 100 p sec for these estimates. Some of the receivers in the global IGS network contain high quality frequency references, such as hydrogen masers. The clock solutions for such receivers are smooth at the 20-p sec level on time scales of a few minutes. There are occasional (daily to weekly) shifts at the microsec level, symptomatic of receiver resets, and 200-p sec-level discontinuities at midnight due to 1-day processing boundaries. Relative clock solutions among 22 IGS sites proposed as "fiducial" in the IGS/BIPM pilot project have been examined over a recent 4-week period. This allows a quantitative measure of receiver reset frequency as a function of site. For days and-sites without resets, the Allan deviation of the relative clock solutions is also computed for subdaily values of tau..

Zumberge, J. F.↗

RFI-SIM: RFI Simulation Package

RFI-SIM simulates the RFI environment to estimate the interference from terrestrial emitters into spacecraft, or vice versa. A high-fidelity simulation of the RFI environment has been developed by employing all antenna- related and radar system-related parameters of multiple emitters, as well as that of the desired spacecraft. In the simulation, the real-time analysis of the interference and its effects on error budgets of a desired radar system is taken into account. This provides a reliable tool for radar system design to deal with RFI issues and to evaluate the sensitivity of various parts of a radar system including antenna pattern, RF front-end and digital processing to RFI signals.

Ghaemi, Hirad↗

Studies of the Effects of Actuator Errors on the HCIT/PIAA Contrast Performance

The High Contrast Imaging Testbed Phase Induced Amplitude Apodization (HCIT/PIAA) coronagraph system at JPL relies on an Electric-Field Conjugation (EFC) wavefront correction algorithm to obtain the required high contrast. This algorithm works with one deformable mirror (DM) to estimate the electric-field to be controlled, and with one or multiple DM's to create a "dark-hole" in the image plane. We have investigated the effects of DM actuator errors on the efficiency of the EFC algorithm. The structural design of the optical system as well as the parameters of various optical elements used in the analysis are drawn from those of the HCIT/PIAA system that have been and will be implemented with one or two DM's. The simulation takes into account the surface errors of various optical elements. In this paper, we report our findings in the case of narrowband wavelength light.

PIAA↗

Surrogate Modeling and Optimization of a Combustor with an Interdigitated Flushwall Injector

Design and Analysis of Computer Experiments (DACE) methods are applied and used to perform surrogate modeling and optimization of a simplified combustor flowpath with an interdigitated flushwall injector. The objectives of the optimization are the thrust potential and combustion efficiency, which are evaluated across a range of flight Mach numbers, duct heights, spanwise spacings, and injection angles. The focus of this work is to highlight the application of a sequential learning approach, in order to learn about the responses of the objective functions over the design space and to identify local regions of interest for further analysis. This approach is contrasted to a previous effort where only a single sampling set was used to fit surrogate models and perform optimization. The optimal solutions resulting from the previous and present approaches are different, due to surrogate model-guided local refinement of the design space allowed by the sequential learning method. The values of the global error estimates between the previous and present approaches are comparable, but the sequential method proved more computationally cost-effective. Further refinement in the optimal regions might be needed to improve predictive capability of surrogate models and to obtain the optimal solutions sets.

Rajiv R Shenoy↗

Investigation of spectral analysis techniques for randomly sampled velocimetry data

It is well known that velocimetry (LV) generates individual realization velocity data that are randomly or unevenly sampled in time. Spectral analysis of such data to obtain the turbulence spectra, and hence turbulence scales information, requires special techniques. The 'slotting' technique of Mayo et al, also described by Roberts and Ajmani, and the 'Direct Transform' method of Gaster and Roberts are well known in the LV community. The slotting technique is faster than the direct transform method in computation. There are practical limitations, however, as to how a high frequency and accurate estimate can be made for a given mean sampling rate. These high frequency estimates are important in obtaining the microscale information of turbulence structure. It was found from previous studies that reliable spectral estimates can be made up to about the mean sampling frequency (mean data rate) or less. If the data were evenly samples, the frequency range would be half the sampling frequency (i.e. up to Nyquist frequency); otherwise, aliasing problem would occur. The mean data rate and the sample size (total number of points) basically limit the frequency range. Also, there are large variabilities or errors associated with the high frequency estimates from randomly sampled signals. Roberts and Ajmani proposed certain pre-filtering techniques to reduce these variabilities, but at the cost of low frequency estimates. The prefiltering acts as a high-pass filter. Further, Shapiro and Silverman showed theoretically that, for Poisson sampled signals, it is possible to obtain alias-free spectral estimates far beyond the mean sampling frequency. But the question is, how far? During his tenure under 1993 NASA-ASEE Summer Faculty Fellowship Program, the author investigated from his studies on the spectral analysis techniques for randomly sampled signals that the spectral estimates can be enhanced or improved up to about 4-5 times the mean sampling frequency by using a suitable prefiltering technique. But, this increased bandwidth comes at the cost of the lower frequency estimates. The studies further showed that large data sets of the order of 100,000 points, or more, high data rates, and Poisson sampling are very crucial for obtaining reliable spectral estimates from randomly sampled data, such as LV data. Some of the results of the current study are presented.

Sree, Dave↗

Numerical analysis of discrete fractional integrodifferential structural dampers

This paper develops solution algorithms enabling the handling of the dynamic response of nonlinear structures contained discretely attached dampers modeled by fractional integrodifferential operators of the Grunwald-Liouville-Riemann type. The development consists of two levels of formulation, namely: (1) numerical approximations of fractional operators and, (2) the establishment of global level implicit schemes enabling the solution to nonlinear structural formulations. To generalize the overall results, error estimates are derived for the fractional operator approximation algorithm. These enable an ongoing optimization of solution efficiency for a given error tolerance. To benchmark the scheme, the results of several numerical experiments are presented. These illustrate the numerical characteristics of the overall formulation.

Padovan, Joseph↗

HIFiRE Flight 2 Unstart Reliability Analysis

The objective of this work was to assess the unstart reliability of the HIFiRE Flight 2 system. To do this, a quantification of margins and uncertainties framework was used for comparing model predictions of the predicted combustion induced shock location to the predicted last stable shock location within the isolator. Uncertainty sources included parametric uncertainty in the flight conditions, the heat release model, and turbulence modeling, as well as model verification errors. Additionally, an estimate of model-form uncertainty was established by comparing the model to measured ground test data. A computationally efficient non-intrusive polynomial chaos approach was used to propagate parametric uncertainty through the computational flu-ids dynamics models of both the ground test configuration and the flight vehicle. Compared to direct-connect ground test data, computational fluid dynamics predictions yielded about two duct heights of model-form uncertainty. This was applied to a prediction of flight vehicle unstart margin at the Mach 6.5 flight condition. Building up all of the computational model uncertainty, including parametric uncertainty, verification errors, and the determined model-form uncertainty, a 95% probability level based confidence ratio, or a ratio of total uncertainty to a statistical margin measure, was found to be 0.31 for the flight system.

Thomas K West↗

Wind speeds in lower atmosphere of Venus: Status report on possible measurement via differential VLBI tracking of entry probes

The potential of very-long-baseline interferometry (VLBI) is examined for use in the determination of wind speeds in Venus' lower atmosphere via the differential tracking of entry probes. A simplified mathematical model is presented in detail. An incomplete error analysis based on this model permits an educated guess to be made: an uncertainty in wind speed determination of no more than about 100/t m/sec, where t l is the corresponding time resolution in seconds, is an achievable goal, without the use of transponders on the miniprobes. If transponders are available on all probes, there should be little difficulty in estimating wind speeds with useful precision.

Shapiro, I. I.↗

A study of oceanic surface heat fluxes in the Greenland, Norwegian, and Barents Seas

This study examines oceanic surface heat fluxes in the Norwegian, Greenland, and Barents seas using the gridded Navy Fleet Numerical Oceanography Central surface analysis and the First GARP Global Experiment (FGGE) IIc cloudiness data bases. Monthly and annual means of net and turbulent heat fluxes are computed for the FGGE year 1979. The FGGE IIb data base consisting of individual observations provides particularly good data coverage in this region for a comparison with the gridded Navy winds and air temperatures. The standard errors of estimate between the Navy and FGGE IIb winds and air temperatures are 3.6 m/s and 2.5 C, respectively. The computations for the latent and sensible heat fluxes are based on bulk formulas with the same constant heat exchange coefficient of 0.0015. The results show extremely strong wintertime heat fluxes in the northern Greenland Sea and especially in the Barents Sea in contrast to previous studies.

Hakkinen, Sirpa↗

Estimation of the diurnal cycle of oceanic precipitation from SSM/I data

A study of differences between the morning and evening monthly rainfall for 5 x 5 deg cells over the oceans from the SSM/I data has been conducted. The monthly rainfalls are estimated from the technique given by Wilheit et al. The difference between the morning and evening monthly rainfall arises due to the various random errors involved in the retrieval process, the sampling error in the observations, and the diurnal component of oceanic rainfall. The diurnal component is weak but clearly visible when averaged over large areas and for long time periods. The analysis shows that morning rainfall is consistently greater than evening rainfall. The Northern Hemisphere seem to have a larger diurnal variation than does the Southern Hemisphere. The maximum ratio between the morning and evening monthly rainfall is 1.7 while 1.2 is the more typical value.

Sharma, Awdhesh K.↗

Variational Assimilation of Global Microwave Rainfall Retrievals: Physical and Dynamical Impact on GEOS Analyses and Forecasts

Global microwave rainfall retrievals from a 5-satellite constellation, including TMI from TRMM, SSWI from DMSP F13, F14 and F15, and AMSR-E from EOS-AQUA, are assimilated into the NASA Goddard Earth Observing System (GEOS) Data Assimilation System (DAS) using a 1-D variational continuous assimilation (VCA) algorithm. The physical and dynamical impact of rainfall assimilation on GEOS analyses and forecasts is examined at various temporal and spatial scales. This study demonstrates that the 1-D VCA algorithm, which was originally developed and evaluated for rainfall assimilations over tropical oceans, can effectively assimilate satellite microwave rainfall retrievals and improve GEOS analyses over both the Tropics and the extratropics where the atmospheric processes are dominated by different large-scale dynamics and moist physics, and also over the land, where rainfall estimates from passive microwave radiometers are believed to be less accurate. Results show that rainfall assimilation renders the GEOS analysis physically and dynamically more consistent with the observed precipitation at the monthly-mean and 6-hour time scales. Over regions where the model precipitation tends to misbehave in distinctly different rainy regimes, the 1-D VCA algorithm, by compensating for errors in the model s moist time-tendency in a 6-h analysis window, is able to bring the rainfall analysis closer to the observed. The radiation and cloud fields also tend to be in better agreement with independent satellite observations in the rainfall-assimilation m especially over regions where rainfall analyses indicate large improvements. Assimilation experiments with and without rainfall data for a midlatitude frontal system clearly indicates that the GEOS analysis is improved through changes in the thermodynamic and dynamic fields that respond to the rainfall assimilation. The synoptic structures of temperature, moisture, winds, divergence, and vertical motion, as well as vorticity are more realistically captured across the front. Short-term forecasts using initial conditions assimilated with rainfall data also show slight improvements. 1

Lin, Xin↗

A Pulsar-Based Time-Scale from the International Pulsar Timing Array

We have constructed a new time-scale, TT(IPTA16), based on observations of radio pulsars presented in the first data release from the International Pulsar Timing Array (IPTA). We used two analysis techniques with independent estimates of the noise models for the pulsar observations and different algorithms for obtaining the pulsar time-scale. The two analyses agree within the estimated uncertainties and both agree with TT(BIPM17), a post-corrected time-scale produced by the Bureau International des Poids et Mesures (BIPM). We show that both methods could detect significant errors in TT(BIPM17) if they were present. We estimate the stability of the atomic clocks from which TT(BIPM17) is derived using observations of four rubidium fountain clocks at the US Naval Observatory. Comparing the power spectrum of TT(IPTA16) with that of these fountain clocks suggests that pulsar-based time-scales are unlikely to contribute to the stability of the best time-scales over the next decade, but they will remain a valuable independent check on atomic time-scales. We also find that the stability of the pulsar-based time-scale is likely to be limited by our knowledge of solar-system dynamics, and that errors in TT(BIPM17) will not be a limiting factor for the primary goal of the IPTA, which is to search for the signatures of nano-Hertz gravitational waves.

G. Hobbs↗