Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Uncertainty Propagation”

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 271 records · Page 15

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.↗

Solar and Heliospheric Data Requirements: Going Further Than L1

Current operational space weather forecasting relies on solar wind observations made by the ACE spacecraft located at the L1 point providing 30-40 minutes warning time. Some use is also made of SOHO and STEREO solar imaging that potentially can give multiple days of warning time. However, our understanding of the propagation and evolution of solar wind transients is still limited resulting in a typical timing uncertainty of approximately 10 hours. In order to improve this critical understanding, a number of NASA missions are being planned. Specifically the Solar Probe Plus and Solar Orbiter missions will investigate the inner Heliospheric evolution of coronal mass ejections and the acceleration and propagation of solar energetic particles. In addition, a number of multi-spacecraft concepts have been studied that have the potential to significantly improve the accuracy of long-term space weather forecasts.

Szabo, A.↗

Interval prediction in structural dynamic analysis

Methods for assessing the predictive accuracy of structural dynamic models are examined with attention given to the effects of modal mass, stiffness, and damping uncertainties. The methods are based on a nondeterministic analysis called 'interval prediction' in which interval variables are used to describe parameters and responses that are unknown. Statistical databases for generic modeling uncertainties are derived from experimental data and incorporated analytically to evaluate responses. Covariance matrices of modal mass, stiffness, and damping parameters are propagated numerically in models of large space structures by means of three methods. The test data tend to fall within the predicted intervals of uncertainty determined by the statistical databases. The present findings demonstrate the suitability of using data from previously analyzed and tested space structures for assessing the predictive accuracy of an analytical model.

Hasselman, Timothy K.↗

Seismic Study of the Dynamics of the Solar Subsurface from SoHO Observations

In collaboration with Dr. Baudin, we have developed and refined the new observational methodology for local helioseismology known as time-distance analysis. Global helioseismology study the solar oscillations as a superposition of resonant modes, whose properties (mode frequencies) reflect the global structure of the sun (sound speed stratification, rotation rate, etc). In contrast, local helioseismology look at the solar oscillations as wave packets whose propagation will be affected by perturbations of the media sampled. These local perturbations (sound speed or velocity flows) will modify the propagation time, that in turn can be used as a diagnostic tool for a given region. From a data reduction perspective, the processing of solar dopplergrams that result in time-distance maps, i.e. propagation times as a function of distance between bounces at the surface, is radically different from the methodology used for global mode analysis. We have, in a first step, further develop the programs needed to carry out such analysis. We have then applied them to NMI data set, and explore the trade-off between various averaging and filtering approaches - steps required to improve the signal-to-noise ratio of correlation maps - and the resulting stability and precision of the fitted propagation times. While excessive averaging (whether over space, propagation distance, or time) will reduce the diagnostic potential of the method, insufficient averaging lead to unstable fits, or uncertainties so large as to hide the information we seek. In a second phase, we have developed the analysis methodology required to infer local properties from perturbation in time propagation. Namely, we have developed time-distance inversion techniques, with an emphasis on inferences of velocity flows from time anomalies. Note also that during the period covered by this grant, all the investigators on this proposal (i.e., Drs. Baudin, Eff-Darwich, Korzennik, and Noyes) took part in the organization of the SOHO 6 /GONG 99 Workshop: Structure and Dynamics of the Interior of the Sun and Sunlike Stars, held on June 1-4 1999 at the Boston Park Plaza Hotel in Boston, Massachusetts, USA. it was very well attended by more than 160 participants from 26 countries from all over the world. The proceedings were published in two volumes as ESA SP-418, with Sessions I-III in Volume 1, and Sessions IV-VI in Volume 2 (1,000 pages in total). The complete contents are also included in digital form on a CD-ROM included with Volume 1. This CD-ROM also contains additional multi-media material that complements some of the contributions.

Korzennik, Sylvain G.↗

The Impact of Model and Rainfall Forcing Errors on Characterizing Soil Moisture Uncertainty in Land Surface Modeling

The contribution of rainfall forcing errors relative to model (structural and parameter) uncertainty in the prediction of soil moisture is investigated by integrating the NASA Catchment Land Surface Model (CLSM), forced with hydro-meteorological data, in the Oklahoma region. Rainfall-forcing uncertainty is introduced using a stochastic error model that generates ensemble rainfall fields from satellite rainfall products. The ensemble satellite rain fields are propagated through CLSM to produce soil moisture ensembles. Errors in CLSM are modeled with two different approaches: either by perturbing model parameters (representing model parameter uncertainty) or by adding randomly generated noise (representing model structure and parameter uncertainty) to the model prognostic variables. Our findings highlight that the method currently used in the NASA GEOS-5 Land Data Assimilation System to perturb CLSM variables poorly describes the uncertainty in the predicted soil moisture, even when combined with rainfall model perturbations. On the other hand, by adding model parameter perturbations to rainfall forcing perturbations, a better characterization of uncertainty in soil moisture simulations is observed. Specifically, an analysis of the rank histograms shows that the most consistent ensemble of soil moisture is obtained by combining rainfall and model parameter perturbations. When rainfall forcing and model prognostic perturbations are added, the rank histogram shows a U-shape at the domain average scale, which corresponds to a lack of variability in the forecast ensemble. The more accurate estimation of the soil moisture prediction uncertainty obtained by combining rainfall and parameter perturbations is encouraging for the application of this approach in ensemble data assimilation systems.

Maggioni, V.↗

External control of the Saturn kilometric radiation by the solar wind - Comparison between Voyager 1 and 2 observations

The long-term modulation of Saturn's nonthermal radio emission in the kilometric wavelength range has been studied based upon data obtained by Voyagers 1 and 2. A comparison of the ballistic and hydrodynamic propagation of solar wind features from the spacecraft to Saturn allows the uncertainty inherent in the projection to be determined. The results confirm the previous suggestion that momentum, ram pressure, and kinetic energy flux are the primary solar wind parameters that drive the nonthermal radio emission. It is suggested that, under certain conditions and for limited periods of time, the magnetic properties and time derivatives of the solar wind have increased importance.

Rucker, H. O.↗

Embedding Temporal Constraints For Coordinated Execution in Habitat Automation

Future NASA plans call for long-duration deep space missions with human crews. Because of light-time delay and other considerations, increased autonomy will be needed. This will necessitate integration of tools in such areas as anomaly detection, diagnosis, planning, and execution. In this paper we investigate an approach that integrates planning and execution by embedding planner-derived temporal constraints in an execution procedure. To avoid the need for propagation, we convert the temporal constraints to dispatchable form. We handle some uncertainty in the durations without it affecting the execution; larger variations may cause activities to be skipped.

Morris, Paul↗

Methodology for error analysis and simulation of lidar aerosol measurements

A methodology is presented for objective and automated determination of the uncertainty in lidar aerosol measurements. This methodology is based on standard error-propagation procedures, a large data base on atmospheric behavior, and long experience in lidar data processing. Algebraic expressions for probable error are derived as a function of the relevant parameters. The validity of these expressions is then tested by making simulated measurements and analyses in which random errors of appropriate size are injected at proper steps of the measurement and analysis process. An illustrative example is given where the methodology is applied to a new lidar system now being used for airborne measurements of the stratospheric aerosol.

Russell, P. B.↗

On Combining Thermal-Infrared and Radio-Occultation Data of Saturn's Atmosphere

Radio-occultation and thermal-infrared measurements are complementary investigations for sounding planetary atmospheres. The vertical resolution afforded by radio occultations is typically approximately 1 km or better, whereas that from infrared sounding is often comparable to a scale height. On the other hand, an instrument like CIRS can easily generate global maps of temperature and composition, whereas occultation soundings are usually distributed more sparsely. The starting point for radio-occultation inversions is determining the residual Doppler-shifted frequency, that is the shift in frequency from what it would be in the absence of the atmosphere. Hence the positions and relative velocities of the spacecraft, target atmosphere, and DSN receiving station must be known to high accuracy. It is not surprising that the inversions can be susceptible to sources of systematic errors. Stratospheric temperature profiles on Titan retrieved from Cassini radio occultations were found to be very susceptible to errors in the reconstructed spacecraft velocities (approximately equal to 1 mm/s). Here the ability to adjust the spacecraft ephemeris so that the profiles matched those retrieved from CIRS limb sounding proved to be critical in mitigating this error. A similar procedure can be used for Saturn, although the sensitivity of its retrieved profiles to this type of error seems to be smaller. One issue that has appeared in inverting the Cassini occultations by Saturn is the uncertainty in its equatorial bulge, that is, the shape in its iso-density surfaces at low latitudes. Typically one approximates that surface as a geopotential surface by assuming a barotropic atmosphere. However, the recent controversy in the equatorial winds, i.e., whether they changed between the Voyager (1981) era and later (after 1996) epochs of Cassini and some Hubble observations, has made it difficult to know the exact shape of the surface, and it leads to uncertainties in the retrieved temperature profiles of one to a few kelvins. This propagates into errors in the retrieved helium abundance, which makes use of thermal-infrared spectra and synthetic spectra computed with retrieved radio-occultation temperature profiles. The highest abundances are retrieved with the faster Voyager-era winds, but even these abundances are somewhat smaller than those retrieved from the thermal-infrared data alone (albeit with larger formal errors). The helium abundance determination is most sensitive to temperatures in the upper troposphere. Further progress may include matching the radio-occultation profiles with those from CIRS limb sounding in the upper stratosphere.

Flasar, F. M.↗

Estimating Sonic Boom Metrics Across a Community Using a Kalman Filter

As part of the Quesst mission, NASA will fly the X‑59 aircraft over selected communities to evaluate community responses to low-intensity sonic booms. The purpose of these community tests is to determine the dose-response relationship between the noise exposure metrics and the community response. The independent variables for the dose-response relationship are the noise exposure metrics experienced by survey respondents within the community. Two sources of noise exposure metrics are available in each community: measurements at sparse locations throughout the community, and calculations from propagation models across the community. Both the measurements and calculations are subject to uncertainty. A Kalman filter is proposed to combine the measured and calculated noise exposure metrics to obtain the best estimate of the true noise exposure metrics across the community. The noise exposure metrics estimated by the Kalman filter have lower uncertainty than either the measured or calculated noise exposure metrics alone. Simulations demonstrate that the Kalman filter produces a more accurate estimate of the true noise exposure metrics than other noise estimation methods.

Sonic boom↗

Quantitative estimates of the effects of cross section uncertainties on the derivation of GCR source composition

Cosmic-ray source abundance uncertainties, resulting from the cross-section uncertainties, are calculated for the elements carbon through nickel. A significant dominance of cross-section errors in the source abundance uncertainties is found for most of the elements considered, even with uncorrelated partial cross-section errors; much larger uncertainties than those reported previously are found for Ca, N, Na, Al, and S. Propagation errors are noted to preclude a significant determination of source abundances for F, Cl, and Mn. A need for an assessment of the existing cosmic-ray propagation models - especially the propagation cross-sections which they employ - is expressed, motivated by recent improvements in the resolution and statistical accuracy of observations.

Hinshaw, G. F.↗

Uncertainty Analysis for the Miniaturized Laser Heterodyne Radiometer (mini-LHR)

Presented here is a sensitivity analysis for the miniaturized laser heterodyne radiometer (mini-LHR). This passive, ground-based instrument measures carbon dioxide (CO2) in the atmospheric column and has been under development at NASA/GSFC since 2009. The goal of this development is to produce a low-cost, easily-deployable instrument that can extend current ground measurement networks in order to (1) validate column satellite observations, (2) provide coverage in regions of limited satellite observations, (3) target regions of interest such as thawing permafrost, and (4) support the continuity of a long-term climate record. In this paper an uncertainty analysis of the instrument performance is presented and compared with results from three sets of field measurements. The signal-to-noise ratio (SNR) and corresponding uncertainty for a single scan are calculated to be 329.4+/-1.3 by deploying error propagation through the equation governing the SNR. Reported is an absorbance noise of 0.0024 for 6 averaged scans of field data, for an instrument precision of approximately 0.2 ppmv for CO2.

carbon dioxide (CO2)↗

Rao-Blackwellization for Adaptive Gaussian Sum Nonlinear Model Propagation

When dealing with imperfect data and general models of dynamic systems, the best estimate is always sought in the presence of uncertainty or unknown parameters. In many cases, as the first attempt, the Extended Kalman filter (EKF) provides sufficient solutions to handling issues arising from nonlinear and non-Gaussian estimation problems. But these issues may lead unacceptable performance and even divergence. In order to accurately capture the nonlinearities of most real-world dynamic systems, advanced filtering methods have been created to reduce filter divergence while enhancing performance. Approaches, such as Gaussian sum filtering, grid based Bayesian methods and particle filters are well-known examples of advanced methods used to represent and recursively reproduce an approximation to the state probability density function (pdf). Some of these filtering methods were conceptually developed years before their widespread uses were realized. Advanced nonlinear filtering methods currently benefit from the computing advancements in computational speeds, memory, and parallel processing. Grid based methods, multiple-model approaches and Gaussian sum filtering are numerical solutions that take advantage of different state coordinates or multiple-model methods that reduced the amount of approximations used. Choosing an efficient grid is very difficult for multi-dimensional state spaces, and oftentimes expensive computations must be done at each point. For the original Gaussian sum filter, a weighted sum of Gaussian density functions approximates the pdf but suffers at the update step for the individual component weight selections. In order to improve upon the original Gaussian sum filter, Ref. [2] introduces a weight update approach at the filter propagation stage instead of the measurement update stage. This weight update is performed by minimizing the integral square difference between the true forecast pdf and its Gaussian sum approximation. By adaptively updating each component weight during the nonlinear propagation stage an approximation of the true pdf can be successfully reconstructed. Particle filtering (PF) methods have gained popularity recently for solving nonlinear estimation problems due to their straightforward approach and the processing capabilities mentioned above. The basic concept behind PF is to represent any pdf as a set of random samples. As the number of samples increases, they will theoretically converge to the exact, equivalent representation of the desired pdf. When the estimated qth moment is needed, the samples are used for its construction allowing further analysis of the pdf characteristics. However, filter performance deteriorates as the dimension of the state vector increases. To overcome this problem Ref. [5] applies a marginalization technique for PF methods, decreasing complexity of the system to one linear and another nonlinear state estimation problem. The marginalization theory was originally developed by Rao and Blackwell independently. According to Ref. [6] it improves any given estimator under every convex loss function. The improvement comes from calculating a conditional expected value, often involving integrating out a supportive statistic. In other words, Rao-Blackwellization allows for smaller but separate computations to be carried out while reaching the main objective of the estimator. In the case of improving an estimator's variance, any supporting statistic can be removed and its variance determined. Next, any other information that dependents on the supporting statistic is found along with its respective variance. A new approach is developed here by utilizing the strengths of the adaptive Gaussian sum propagation in Ref. [2] and a marginalization approach used for PF methods found in Ref. [7]. In the following sections a modified filtering approach is presented based on a special state-space model within nonlinear systems to reduce the dimensionality of the optimization problem in Ref. [2]. First, the adaptive Gaussian sum propagation is explained and then the new marginalized adaptive Gaussian sum propagation is derived. Finally, an example simulation is presented.

state estimation↗

Charge-exchange plasma environment for an ion drive spacecraft

The charge exchange plasma environment around a spacecraft that uses mercury ion thrusters for propulsion is described. The interactions between the plasma environment and the spacecraft are determined and a model which describes the propagation of the mercury charge exchange plasma is discussed. The model is extended to describe the flow of the molybdenum component of the charge exchange plasma. The uncertainties in the models for various conditions are discussed and current drain to the solar array, charge exchange plasma material deposition, and the effects of space plasma on the charge exchange plasma propagation are addressed.

Kaufman, H. R.↗

Charge-exchange plasma environment for an ion drive spacecraft

A model was reviewed which describes the propagation of the mercury charge-exchange plasma and extended to describe the flow of the molybdenum component of the charge-exchange plasma. The uncertainties in the models for various conditions are discussed. Such topics as current drain to the solar array, charge-exchange plasma material deposition, and the effects of space plasma on the charge-exchange plasma propagation are addressed.

Kaufman, H. R.↗

Use of In Situ and Airborne Multiangle Data to Assess MODIS- and Landsat-based Estimates of Surface Albedo

The quantification of uncertainty of global surface albedo data and products is a critical part of producing complete, physically consistent, and decadal land property data records for studying ecosystem change. A current challenge in validating satellite retrievals of surface albedo is the ability to overcome the spatial scaling errors that can contribute on the order of 20% disagreement between satellite and field-measured values. Here, we present the results from an uncertain ty analysis of MODerate Resolution Imaging Spectroradiometer (MODIS) and Landsat albedo retrievals, based on collocated comparisons with tower and airborne multi-angular measurements collected at the Atmospheric Radiation Measurement Program s (ARM) Cloud and Radiation Testbed (CART) site during the 2007 Cloud and Land Surface Interaction Campaign (CLAS33 IC 07). Using standard error propagation techniques, airborne measurements obtained by NASA s Cloud Absorption Radiometer (CAR) were used to quantify the uncertainties associated with MODIS and Landsat albedos across a broad range of mixed vegetation and structural types. Initial focus was on evaluating inter-sensor consistency through assessments of temporal stability, as well as examining the overall performance of satellite-derived albedos obtained at all diurnal solar zenith angles. In general, the accuracy of the MODIS and Landsat albedos remained under a 10% margin of error in the SW(0.3 - 5.0 m) domain. However, results reveal a high degree of variability in the RMSE (root mean square error) and bias of albedos in both the visible (0.3 - 0.7 m) and near-infrared (0.3 - 5.0 m) broadband channels; where, in some cases, retrieval uncertainties were found to be in excess of 20%. For the period of CLASIC 07, the primary factors that contributed to uncertainties in the satellite-derived albedo values include: (1) the assumption of temporal stability in the retrieval of 500 m MODIS BRDF values over extended periods of cloud-contaminated observations; and (2) the assumption of spatial 45 and structural uniformity at the Landsat (30 m) pixel scale.

Roman, Miguel O.↗

A secondary tracer approach to the derivation of galactic cosmic-ray source isotopic abundances

A formalism has been developed for deriving cosmic-ray source isotopic abundances from observed local abundances using a purely secondary nuclide as a tracer of spallation production during propagation. Although the formalism is based on the leaky-box model of cosmic-ray propagation, it is shown that source abundances derived by the tracer technique are reasonably independent of detailed propagation models. The tracer formalism also permits a quantitative evaluation of the effects of observational uncertainties on deduced source abundances. It is shown that statistical errors in the observed abundances and uncertainties in the spallation cross sections are at present the dominant sources of uncertainty. The latter error can be reduced with increased detector size or exposure time, while the former can be minimized by measurements of the relative production cross sections. As a specific example, the tracer technique is applied to the isotopes of sulfur and calcium, and the level of uncertainties which must be achieved to distinguish evolutionary differences between solar-system material and cosmic ray-source material are established.

Stone, E. C.↗