Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “variance 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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Frequency domain analysis of errors in cross-correlations of ambient seismic noise

We analyse random errors (variances) in cross-correlations of ambient seismic noise in the frequency domain, which differ from previous time domain methods. Extending previous theoretical results on ensemble averaged cross-spectrum, we estimate confidence interval of stacked cross-spectrum of finite amount of data at each frequency using non-overlapping windows with fixed length. The extended theory also connects amplitude and phase variances with the variance of each complex spectrum value. Analysis of synthetic stationary ambient noise is used to estimate the confidence interval of stacked cross-spectrum obtained with different length of noise data corresponding to different number of evenly spaced windows of the same duration. This method allows estimating Signal/Noise Ratio (SNR) of noise cross-correlation in the frequency domain, without specifying filter bandwidth or signal/noise windows that are needed for time domain SNR estimations. Based on synthetic ambient noise data, we also compare the probability distributions, causal part amplitude and SNR of stacked cross-spectrum function using one-bit normalization or pre-whitening with those obtained without these pre-processing steps. Natural continuous noise records contain both ambient noise and small earthquakes that are inseparable from the noise with the existing pre-processing steps. Using probability distributions of random cross-spectrum values based on the theoretical results provides an effective way to exclude such small earthquakes, and additional data segments (outliers) contaminated by signals of different statistics (e.g. rain, cultural noise), from continuous noise waveforms. This technique is applied to constrain values and uncertainties of amplitude and phase velocity of stacked noise cross-spectrum at different frequencies, using data from southern California at both regional scale (~35 km) and dense linear array (~20 m) across the plate-boundary faults. A block bootstrap resampling method is used to account for temporal correlation of noise cross-spectrum at low frequencies (0.05–0.2 Hz) near the ocean microseismic peaks.

58 GEOSCIENCES↗
DOE OSTI2018

Evapotranspiration simulations in ISIMIP2a—Evaluation of spatio-temporal characteristics with a comprehensive ensemble of independent datasets

Actual land evapotranspiration (ET) is a key component of the global hydrological cycle and an essential variable determining the evolution of hydrological extreme events under different climate change scenarios. However, recently available ET products show persistent uncertainties that are impeding a precise attribution of human-induced climate change. Here, we aim at comparing a range of independent global monthly land ET estimates with historical model simulations from the global water, agriculture, and biomes sectors participating in the second phase of the Inter-Sectoral Impact Model Intercomparison Project (ISIMIP2a). Among the independent estimates, we use the EartH2Observe Tier-1 dataset (E2O), two commonly used reanalyses, a pre-compiled ensemble product (LandFlux-EVAL), and an updated collection of recently published datasets that algorithmically derive ET from observations or observations-based estimates (diagnostic datasets). A cluster analysis is applied in order to identify spatiotemporal differences among all datasets and to thus identify factors that dominate overall uncertainties. The clustering is controlled by several factors including the model choice, the meteorological forcing used to drive the assessed models, the data category (models participating in the different sectors of ISIMIP2a, E2O models, diagnostic estimates, reanalysis-based estimates or composite products), the ET scheme, and the number of soil layers in the models. By using these factors to explain spatial and spatio-temporal variabilities in ET, we and that the model choice mostly dominates (24%{40% of variance explained), except for spatio-temporal patterns of total ET, where the forcing explains the largest fraction of the variance (29%). The most dominant clusters of datasets are further compared with individual diagnostic and reanalysis-based estimates to assess their representation of selected heat waves and droughts in the Great Plains, Central Europe and Western Russia. Although most of the ET estimates capture these extreme events, the generally large spread among the entire ensemble indicates substantial uncertainties.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2022

Deep Interacting Multiple Model Filtering

In this paper, a deep learning-based multiple model estimation framework is presented for the state estimation of hybrid dynamical systems from high dimensional observations such as camera images. A low dimensional vector which represents the measurement of the latent dynamical system and its corresponding variance are learned using a deep encoder neural network. An Interacting Multiple Model (IMM) filter is used to generate the latent state estimates and covariances using multiple dynamical models, which can be learned using backpropagation through time. The state estimates of the dynamical system and the corresponding covariance matrix are generated from the latent state estimates and covariance using a deep decoder neural network. The whole network is trained in an end-to-end manner using a loss function which minimizes the negative log-likelihood of the neural network parameters. Simulation results are presented using a 2D bouncing ball example and estimation error statistics are computed which demonstrates the accuracy and consistency of the estimation.

Ghananeel Rotithor↗
DOE OSTI2020

Singular value decomposition of adjoint flux distributions for Monte Carlo variance reduction

Monte Carlo (MC) shielding calculations often use weight windows (WWs) and biased sources formed from a deterministic estimate of the adjoint flux to improve the convergence rate of tallies. This requires a significant amount of computer memory, which can limit the memory available for high-resolution tally output. A new method is proposed for reducing these memory requirements by using singular value decomposition (SVD) in linear or logarithmic space to approximate the adjoint flux. This method’s performance is evaluated using the Shift and Denovo codes for streaming and diffusion base case problems, followed by problems using the Westinghouse AP1000 and the Joint European Torus. The log SVD reduced WW memory requirements by an order of magnitude in all cases without a significant performance penalty. Additionally, the linear SVD reduced biased source memory requirements by an order of magnitude, but further investigation is needed to account for observed limitations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Utilizing the maximum likelihood estimator for flow analysis

We explore the possibility of evaluating flow harmonics by employing the maximum likelihood estimator (MLE). For a given finite multiplicity, the MLE simultaneously furnishes estimations for all the parameters of the underlying distribution function while efficiently suppressing the variance of measures. Also, the method provides a means to assess a specific class of mixed harmonics, which is not straightforwardly feasible by the approaches primarily based on particle correlations. The results are analyzed using the Wald, likelihood ratio, and score tests of hypotheses. Besides, the resultant flow harmonics obtained using MLE are compared with those derived using particle correlations and event plane methods. Here, the dependencies of extracted flow harmonics on the multiplicity of individual events and the total number of events are analyzed. It is shown that the proposed approach works efficiently to deal with the deficiency in detector acceptability. Moreover, we elaborate on a fictitious scenario where the event plane is not a well-defined quantity in the distribution function. For the latter case, the MLE is shown to largely perform better than the two-particle correlation estimator. In this regard, one concludes that the MLE furnishes a meaningful alternative to the existing approaches for flow analysis.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2015

The Pattern Across the Continental United States of Evapotranspiration Variability Associated with Water Availability

The spatial pattern across the continental United States of the interannual variance of warm season water-dependent evapotranspiration, a pattern of relevance to land-atmosphere feedback, cannot be measured directly. Alternative and indirect approaches to estimating the pattern, however, do exist, and given the uncertainty of each, we use several such approaches here. We first quantify the water dependent evapotranspiration variance pattern inherent in two derived evapotranspiration datasets available from the literature. We then search for the pattern in proxy geophysical variables (air temperature, stream flow, and NDVI) known to have strong ties to evapotranspiration. The variances inherent in all of the different (and mostly independent) data sources show some differences but are generally strongly consistent they all show a large variance signal down the center of the U.S., with lower variances toward the east and (for the most part) toward the west. The robustness of the pattern across the datasets suggests that it indeed represents the pattern operating in nature. Using Budykos hydroclimatic framework, we show that the pattern can largely be explained by the relative strength of water and energy controls on evapotranspiration across the continent.

Air Temperature↗
DOE OSTI2018

Tree cover shows strong sensitivity to precipitation variability across the global tropics

Abstract Aim Vegetation is sensitive to mean annual precipitation (MAP), but the sensitivity of vegetation to precipitation variability (PV) is less clear. Tropical ecosystems are likely to experience increased PV in the future. Here we assessed the importance, magnitude and mechanism of PV effects on tree cover in the context of covarying environmental drivers such as fire, temperature and soil properties. Location Tropical land. Time period 2000–2010. Major taxa studied Trees. Methods We compiled climate, soil and remotely‐sensed tree cover data over tropical land. We then comprehensively assessed the contribution of PV at different time‐scales to tropical tree cover variations and estimated the sensitivity of tree cover to PV changes by conducting rolling‐window regression and variance decomposition analyses. We further adopted a mechanistic modelling approach to test whether water competition between trees and grasses can explain the observed effect of PV. Results We find that PV contributes 33–56% to the total explained spatial variation (65–79%) in tree cover. The contribution of PV depends on MAP and is highest under intermediate MAP (500–1,500 mm). Tree cover generally increases with rainy day frequency and wet season length but shows mixed responses to inter‐annual PV. Based on the estimated sensitivity, tropical tree cover can decrease by 3–5% overall and by up to 20% in Amazonia under a 20% decrease in rainy days. Mechanistic modelling analysis reproduced the continental differences in tree cover along an MAP gradient. Main conclusions Under intermediate rainfall regimes (500–1,500 mm), PV can be a more important determinant of tropical tree cover than conventionally proposed drivers such as MAP and fire. The effect of PV likely results from the sensitivity of tree–grass competition to the temporal distribution of water resources. These results show that climate variability can strongly shape the biosphere.

Xu, Xiangtao↗
NASA NTRS1990

Estimates of daily rainfall over the Amazon basin

Five geostationary satellite rain estimation techniques were tested over Amazonia. Individually, the techniques explained 1/4 to 1/3 of the variance of daily gage rainfall. Based in large part on cost, one technique, which involves a nonlinear relation in temperature, was selected to provide a mapping of daily Amazonia rainfall between May 6 and 12, 1987. Accumulated over the 7 days, rainfall by this technique averaged 40 mm. It varied from zero in the southeast to more than 150 mm in the northwest. To the southwest the predominantly convective pattern of the rain image was overlaid by a streakines, implying some baroclinic influence. In maps combining gage observations with satellite estimates, rainfall varied significantly from day to day. Only over the largest scale did a trend emerge: a tendency for rain to withdraw from south to north.

Martin, David W.↗
NASA NTRS1976

Linear minimum variance filters applied to carrier tracking

A new approach is taken to the problem of tracking a fixed amplitude signal with a Brownian-motion phase process. Classically, a first-order phase-lock loop (PLL) is used; here, the problem is treated via estimation of the quadrature signal components. In this space, the state dynamics are linear with white multiplicative noise. Therefore, linear minimum-variance filters, which have a particularly simple mechanization, are suggested. The resulting error dynamics are linear at any signal/noise ratio, unlike the classical PLL. During synchronization, and above threshold, this filter with constant gains degrades by 3 per cent in output rms phase error with respect to the classical loop. However, up to 80 per cent of the maximum possible noise improvement is obtained below threshold, where the classical loop is nonoptimum, as demonstrated by a Monte Carlo analysis. Filter mechanizations are presented for both carrier and baseband operation.

Gustafson, D. E.↗
NASA NTRSDate not supplied

Evaluating Combined Lidar and Polarimeter Measurements of Cloud Top Parameters (Extinction, Scattering Cross Sections, and Droplet Number Density, Liquid Water Content)

We present a new method to derive profiles of extinction from a lidar which is then combined with polarimeter measurements to derive cloud droplet number density (CDNC) in the tops of warm clouds. The method employs polarization-sensitive elastic backscatter lidar measurements to estimate attenuation of the lidar signal within the cloud and polarimeter estimates of cloud droplet size distributions. The measurements used for this demonstration are from NASA Langley Research Center’s High Spectral Resolution Lidar – Generation 2 (HSRL-2) and NASA GISS’s Research Scanning Polarimeter (RSP). The measurements were acquired on NASA’s ACTIVATE mission, during which the instruments were deployed in a down-looking mode from a high-altitude aircraft, which flew in coordination with a low-flying aircraft acquiring coincident in situ measurements of cloud droplet size and number. The high vertical (1.25 m) and horizontal (~50 m) sampling resolution of the HSRL-2 data enabled retrievals of single-scattering extinction profiles to within ~2.5 optical depths of cloud top. Another key feature of the method was the well-calibrated measurement of backscatter at cloud top due to using the HSRL technique. The RSP retrievals of cloud droplet effective radius and variance were accomplished using the “cloud-bow” technique. Overall, the technique provides extinction profile estimates for the top 2.5 optical depths of water clouds that can be used to estimate the cloud droplet number density in various types of warm clouds. The mean extinction values at cloud top (0-1 optical depth) are compared against cloud drop size and number concentration acquired from wing-mounted probes (e.g., DMT Cloud Droplet Probe - CDP, SPEC Fast Cloud Droplet Probe FCDP, and SPEC 2D Stereo Probe - 2DS) that were deployed on the low-flying aircraft. Comparisons between the effective radius and variance from the polarimeter, lidar ratio (extinction to backscatter) from the lidar and polarimeter, and LWC are also presented. All measurements were acquired over the Western North Atlantic over 3 years from 2020 to 2022.

lidar↗
DOE OSTI2016

Angular dependence of primordial trispectra and CMB spectral distortions

Under the presence of anisotropic sources in the inflationary era, the trispectrum of the primordial curvature perturbation has a very specific angular dependence between each wavevector that is distinguishable from the one encountered when only scalar fields are present, characterized by an angular dependence described by Legendre polynomials. We examine the imprints left by curvature trispectra on the TT μ bispectrum, generated by the correlation between temperature anisotropies (T) and chemical potential spectral distortions (μ) of the Cosmic Microwave Background (CMB). Due to the angular dependence of the primordial signal, the corresponding TT μ bispectrum strongly differs in shape from TT μ sourced by the usual g {sub NL} or τ{sub NL} local trispectra, enabling us to obtain an unbiased estimation. From a Fisher matrix analysis, we find that, in a cosmic-variance-limited (CVL) survey of TT μ, a minimum detectable value of the quadrupolar Legendre coefficient is d {sub 2} ∼ 0.01, which is 4 orders of magnitude better than the best value attainable from the TTTT CMB trispectrum. In the case of an anisotropic inflationary model with a f (φ) F {sup 2} interaction (coupling the inflaton field φ with a vector kinetic term F {sup 2}), the size of the curvature trispectrum is related to that of quadrupolar power spectrum asymmetry, g {sub *}. In this case, a CVL measurement of TT μ makes it possible to measure g {sub *} down to 10{sup −3}.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Analytic marginalization of N(z) uncertainties in tomographic galaxy surveys

In this paper, we present a new method to marginalize over uncertainties in redshift distributions, N(z), within tomographic cosmological analyses applicable to current and upcoming photometric galaxy surveys. We allow for arbitrary deviations from the best-guess N(z) governed by a general covariance matrix describing the uncertainty in our knowledge of redshift distributions. In principle, this is marginalization over hundreds or thousands of new parameters describing potential deviations as a function of redshift and tomographic bin. However, by linearly expanding the theory predictions around a fiducial model, this marginalization can be performed analytically, resulting in a modified data covariance matrix that effectively downweights the modes of the data vector that are more sensitive to redshift distribution variations. We showcase this method by applying it to the galaxy clustering measurements from the Hyper Suprime-Cam first data release. We illustrate how to marginalize over sample variance of the calibration sample and a large general systematic uncertainty in photometric estimation methods, and explore the impact of priors imposing smoothness in the redshift distributions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Ground-state properties of light 4 n self-conjugate nuclei in ab initio no-core Monte Carlo shell model calculations with nonlocal N N interactions

We report J π = 0 + ground-state energies and point-proton radii of 4 He, 8 Be, 12 C, 16 O, and 20 Ne nuclei calculated by the ab initio no-core Monte Carlo shell model with the JISP16 and Daejeon16 nonlocal NN interactions. Ground-state energies are obtained in the basis spaces up to seven oscillator shells (N shell =7) with several oscillator energies (ℏω) around the optimal oscillator energy for the convergence of ground-state energies. These energy eigenvalues are extrapolated to obtain estimates of converged ground-state energies in each basis space using energy variances of computed energy eigenvalues. We further extrapolate these energy-variance-extrapolated energies obtained in the finite basis spaces to infinite basis-space results with an empirical exponential form. This form features a dependence on the basis-space size but is independent of the value of ℏω used for the harmonic-oscillator basis functions. Point-proton radii for these states of atomic nuclei are also calculated following techniques employed for the energies. From these results, it is found that the Daejeon16 NN interaction provides good agreement with experimental data up to approximately 16 O, while the JISP16 NN interaction provides good agreement with experimental data up to approximately 12 C. Beyond these nuclei, the interactions produce overbinding accompanied by radii that are too small. These findings suggest and encourage further revisions of nonlocal NN interactions towards the investigation of nuclear structure in heavier-mass regions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2015

Accelerated Sensitivity Analysis in High-Dimensional Stochastic Reaction Networks

Existing sensitivity analysis approaches are not able to handle efficiently stochastic reaction networks with a large number of parameters and species, which are typical in the modeling and simulation of complex biochemical phenomena. In this paper, a two-step strategy for parametric sensitivity analysis for such systems is proposed, exploiting advantages and synergies between two recently proposed sensitivity analysis methodologies for stochastic dynamics. The first method performs sensitivity analysis of the stochastic dynamics by means of the Fisher Information Matrix on the underlying distribution of the trajectories; the second method is a reduced-variance, finite-difference, gradient-type sensitivity approach relying on stochastic coupling techniques for variance reduction. Here in this paper we demonstrate that these two methods can be combined and deployed together by means of a new sensitivity bound which incorporates the variance of the quantity of interest as well as the Fisher Information Matrix estimated from the first method. The first step of the proposed strategy labels sensitivities using the bound and screens out the insensitive parameters in a controlled manner. In the second step of the proposed strategy, a finite-difference method is applied only for the sensitivity estimation of the (potentially) sensitive parameters that have not been screened out in the first step. Results on an epidermal growth factor network with fifty parameters and on a protein homeostasis with eighty parameters demonstrate that the proposed strategy is able to quickly discover and discard the insensitive parameters and in the remaining potentially sensitive parameters it accurately estimates the sensitivities. The new sensitivity strategy can be several times faster than current state-of-the-art approaches that test all parameters, especially in “sloppy” systems. In particular, the computational acceleration is quantified by the ratio between the total number of parameters over the number of the sensitive parameters.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1984

An analysis of the AVE-SESAME I period using statistical structure and correlation functions

Structure and correlation functions are used to describe atmospheric variability during the 10-11 April day of AVE-SESAME 1979 that coincided with the Red River Valley tornado outbreak. The special mesoscale rawinsonde data are employed in calculations involving temperature, geopotential height, horizontal wind speed and mixing ratio. Functional analyses are performed in both the lower and upper troposphere for the composite 24 h experiment period and at individual 3 h observation times. Results show that mesoscale features are prominent during the composite period. Fields of mixing ratio and horizontal wind speed exhibit the greatest amounts of small-scale variance, whereas temperature and geopotential height contain the least. Results for the nine individual times show that small-scale variance is greatest during the convective outbreak. The functions also are used to estimate random errors in the rawinsonde data. Finally, sensitivity analyses are presented to quantify confidence limits of the structure functions.

Fuelberg, H. E.↗
NASA NTRSDate not supplied

Redundancy: How Many Unreliable Spares are Needed for High Reliability and Confidence on a Time Limited Mission?

This paper investigates the number of redundant units needed to achieve high reliability with high confidence. The approach applies to the case where the unit failure rate is too high for a single unit to provide the required reliability over the mission duration. To achieve high reliability, the design then uses N redundant units, one operating unit and N – 1 spares. If the unit failure rate is f, the mission length is L, and f * L is small (not the case assumed here), the unit failure probability over the mission duration is F1 = f * L << 1. In this case, the probability that all N units will fail is FN = F1N, and the needed N = LN(FN)/LN(F1). For the case of large f * L assumed here, F1 = f * L > 1, and F1 is the expected number of failures during the mission. The needed redundancy, N, to achieve the specified N unit reliability, FN, can be computed using the cumulative Poisson distribution with mean equal to F1. The number of spares, N - 1, is increased until the probability - that the total number of failures will be less than N -1 - achieves the required reliability. The confidence that this reliability can be achieved can be computed using the cumulative Poisson distribution or the chi-square distribution. Since the measured unit failure rate, f, has some uncertainty, the confidence that the rate is not lower than the actual failure rate and the required reliability is not overestimated is about 50%. Adding more redundant units increases the confidence that the required reliability, FN, will be achieved. For a fixed number of redundant units, the expected reliability and confidence can be traded off, since lower reliability goals have higher confidence in being achieved. Both the required reliability and confidence can be specified initially and the needed number of redundant units computed using the measured failure rate. The unit failure rate is determined by initial reliability growth testing to remove design errors and to better estimate the final constant failure rate. Reducing the failure rate and reducing its variance both reduce the number of redundant units needed for the required reliability and confidence. Since the total cost is the sum of the costs of the units and of the testing, there is an optimum test time that produces minimum cost.

Harry W. Jones↗
NASA NTRS2014

Generalized Linear Covariance Analysis

This talk presents a comprehensive approach to filter modeling for generalized covariance analysis of both batch least-squares and sequential estimators. We review and extend in two directions the results of prior work that allowed for partitioning of the state space into solve-for'' and consider'' parameters, accounted for differences between the formal values and the true values of the measurement noise, process noise, and textita priori solve-for and consider covariances, and explicitly partitioned the errors into subspaces containing only the influence of the measurement noise, process noise, and solve-for and consider covariances. In this work, we explicitly add sensitivity analysis to this prior work, and relax an implicit assumption that the batch estimator's epoch time occurs prior to the definitive span. We also apply the method to an integrated orbit and attitude problem, in which gyro and accelerometer errors, though not estimated, influence the orbit determination performance. We illustrate our results using two graphical presentations, which we call the variance sandpile'' and the sensitivity mosaic,'' and we compare the linear covariance results to confidence intervals associated with ensemble statistics from a Monte Carlo analysis.

Navagation↗