Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “covariance and correlation”

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

Development and application of two-step uncertainty propagation and sensitivity analysis methodology for fast reactor safety analysis

Uncertainty quantification (UQ) in nuclear reactors for transients is directly linked with safety assessment through the cross-sections uncertainties, provided as a covariance matrix, which are propagated through the reactor system to output of interest pertaining to reactor safety, such as peak temperatures in fuel/clad/coolant. Using a two-step approach, uncertainties are first quantified and propagated from basic input variables (such as reaction cross-sections) to intermediate quantities (such as reactivity feedback coefficients) through lattice level calculations. Uncertainties of intermediate quantities (from the first step) are then propagated through the system transient calculations, in the second step, to obtain uncertainties on reactor safety output parameters of interest. The scope of this work consists of Uncertainty Quantification & Propagation of nuclear data uncertainties that are highly correlated through unprotected transient overpower and unprotected loss of flow to assess their impact on core safety parameters. This two-step approach in the presence of covariance renders the sensitivity analysis very challenging. In fact, usually the sensitivity analysis is restricted to each step, which limits its application since the sensitivities between the system output quantities and the basic input variables are difficult to obtain. Here, in this work, we address this issue by proposing a simple, general methodology to combine the sensitivity indices obtained in each step by assuming the model behavior being linear. For the first step Generalized Perturbation theory based indices are used while in the second step the recently studied Johnson indices. The uncertainty quantification and sensitivity methodologies discussed here are demonstrated on a generic LFR design which is based on the 500 MWth demonstration Lead-cooled fast reactor (DLFR) using oxide fuel, developed by Westinghouse Electric Company (WEC).

42 - ENGINEERING↗

Statistical temperature profile retrievals in clear-air using passive 118-GHz O2 observations

The clean-air temperature profile accuracy yielded by a localized linear statistical retrieval operator applied to passive aircraft-based 118-GHz spectra is demonstrated. A comparison of the statistically and physically derived correlation coefficients of antenna temperature and kinetic temperature furnishes a physical justification of the statistical retrieval technique. The atmospheric temperature mean and covariance significantly depend on such geophysical parameters as latitude, longitude, local season, and time, as well as the prevailing meteorological state and orographic effects.

Gasiewski, A. J.↗

Observations of geographically correlated orbit errors for TOPEX/Poseidon using the global positioning system

We have compared Global Positioning System (GPS)-based dynamic and reduced-dynamic TOPEX/Poseidon orbits over three 10-day repeat cycles of the ground-track. The results suggest that the prelaunch joint gravity model (JGM-1) introduces geographically correlated errors (GCEs) which have a strong meridional dependence. The global distribution and magnitude of these GCEs are consistent with a prelaunch covariance analysis, with estimated and predicted global rms error statistics of 2.3 and 2.4 cm rms, respectively. Repeating the analysis with the post-launch joint gravity model (JGM-2) suggests that a portion of the meridional dependence observed in JGM-1 still remains, with global rms error of 1.2 cm.

Christensen, E. J.↗

A Framework for Propagation of Uncertainties in the Kepler Data Analysis Pipeline

The Kepler space telescope is designed to detect Earth-like planets around Sun-like stars using transit photometry by simultaneously observing 100,000 stellar targets nearly continuously over a three and a half year period. The 96-megapixel focal plane consists of 42 charge-coupled devices (CCD) each containing two 1024 x 1100 pixel arrays. Cross-correlations between calibrated pixels are introduced by common calibrations performed on each CCD requiring downstream data products access to the calibrated pixel covariance matrix in order to properly estimate uncertainties. The prohibitively large covariance matrices corresponding to the ~75,000 calibrated pixels per CCD preclude calculating and storing the covariance in standard lock-step fashion. We present a novel framework used to implement standard propagation of uncertainties (POU) in the Kepler Science Operations Center (SOC) data processing pipeline. The POU framework captures the variance of the raw pixel data and the kernel of each subsequent calibration transformation allowing the full covariance matrix of any subset of calibrated pixels to be recalled on-the-fly at any step in the calibration process. Singular value decomposition (SVD) is used to compress and low-pass filter the raw uncertainty data as well as any data dependent kernels. The combination of POU framework and SVD compression provide downstream consumers of the calibrated pixel data access to the full covariance matrix of any subset of the calibrated pixels traceable to pixel level measurement uncertainties without having to store, retrieve and operate on prohibitively large covariance matrices. We describe the POU Framework and SVD compression scheme and its implementation in the Kepler SOC pipeline.

Clarke, Bruce D.↗

Multi-Pass Sequential Mini-Batch Stochastic Gradient Descent Algorithms for Noise Covariance Estimation in Adaptive Kalman Filtering

Estimation of unknown noise covariances in a Kalman filter is a problem of significant practical interest in a wide array of applications. Although this problem has a long history, reliable algorithms for their estimation were scant, and necessary and sufficient conditions for identifiability of the covariances were in dispute until recently. Necessary and sufficient conditions for covariance estimation and a batch estimation algorithm were presented in our previous study. This paper presents stochastic gradient descent algorithms for noise covariance estimation in adaptive Kalman filters that are an order of magnitude faster than the batch method for similar or better root mean square error. More significantly, these algorithms are applicable to non-stationary systems where the noise covariances can occasionally jump up or down by an unknown magnitude. The computational efficiency of the new algorithms stems from adaptive thresholds for convergence, recursive fading memory estimation of the sample cross-correlations of the innovations, and accelerated stochastic gradient descent algorithms. The comparative evaluation of the proposed methods on a number of test cases demonstrates their computational efficiency and accuracy.

Adaptive Kalman filtering↗

Monitoring agroecosystem productivity and phenology at a national scale: A metric assessment framework

Effective measurement of seasonal variations in the timing and amount of production is critical to managing spatially heterogeneous agroecosystems in a changing climate. Although numerous technologies for such measurements are available, their relationships to one another at a continental extent are unknown. Using data collected from across the Long-Term Agroecosystem Research (LTAR) network and other networks, we investigated correlations among key metrics representing primary production, phenology, and carbon fluxes in croplands, grazing lands, and crop-grazing integrated systems across the continental U.S. Metrics we examined included gross primary productivity (GPP) estimated from eddy covariance (EC) towers and modelled from the Landsat satellite, Landsat NDVI, and vegetation greenness (Green Chromatic Coordinate, GCC) from tower-mounted PhenoCams for 2017 and 2018. Overall, our analysis compared production dynamics estimated from three independent ground and remote platforms using data for 34 agricultural sites constituting 51 site-years of co-located time series. Pairwise sensor comparisons across all four metrics revealed stronger correlation and lower root mean square error (RMSE) between end of season (EOS) dates (Pearson R ranged from 0.6 to 0.7 and RMSE from 32.5 to 67.8) than start of season (SOS) dates (0.46 to 0.69 and 40.4 to 66.2). Overall, moderate to high correlations between SOS and EOS metrics complemented one another except at some lower productivity grazing land sites where estimating SOS can be challenging. Growing season length estimates derived from 16-day satellite GPP (179.1 days) were significantly longer than those from PhenoCam G CC (70.4 days, p adj < 0.0001) and EC GPP (79.6 days, p adj < 0.0001). Landscape heterogeneity did not explain differences in SOS and EOS estimates. Annual integrated estimates of productivity from EC GPP and PhenoCam G CC diverged from those estimated by Landsat GPP and NDVI at sites where annual production exceeds 1000 gC/m –2 yr –1 . Based on our results, we developed a “metric assessment framework” that articulates where and how metrics from satellite, eddy covariance and PhenoCams complement, diverge from, or are redundant with one another. The framework was designed to optimize instrumentation selection for monitoring, modeling, and forecasting ecosystem functioning with the ultimate goal of informing decision-making by land managers, policy-makers, and industry leaders working at multiple scales.

54 ENVIRONMENTAL SCIENCES↗

Scaling Analysis of Two–Phase Flow in Fractal Permeability Fields

Fluid mixing in permeable media is essential in many practical applications. The mixing process is a consequence of velocity fluctuations owing to geological heterogeneities and mobility contrast of fluids. Heterogeneities in natural rocks are often spatially correlated, and their properties, such as permeability, may be described using fractal distributions. This work models the fractal characteristics of such permeability fields in which the covariance function is expressed as a power-law function. A generalized scaling relation is derived relating various fractal permeability fields using the magnitude of their fluctuations. Here, this relation reveals the self-similar behavior of two-phase flow in such permeable media. To that end, a recently developed, high-resolution numerical simulator is employed to validate the analytically derived scaling relations. Two flow problems are considered in which flow is governed by 1) a linear, and 2) a nonlinear transport equation. Due to the probabilistic representation of the fractal permeability fields, a sensitivity study is conducted for each flow scenario to determine the number of realizations required for statistical convergence. Scaling analysis is performed using ensemble averages of simulated saturation profiles and their mixing lengths. Results support the validity of the developed scaling relation across the range of investigated flow conditions at intermediate times. The dynamics of linear flow in the asymptotic regime is affected by the correlation structure of heterogeneity. In nonlinear flow, scaling behavior appears to be dominated by the degree of nonlinearity.

58 GEOSCIENCES↗

Geographically correlated orbit error for the Topex satellite using GPS tracking

This paper reports on ongoing investigations of geographically correlated orbit error performed in support of precise orbit determination for the Topex/Poseidon mission. Tracking data is assumed to be generated with the Global Positioning System. Geographically correlated orbit error due to errors in modeling the gravity field is described and demonstrated. A method for rational selection of the gravity model coefficients which produce significant orbit error is also presented. The effects of these coefficient errors and other error sources is evaluated using a consider covariance approach. Finally, the covariance analysis is used to evaluate the orbit accuracy for Topex using dynamic, reduced dynamic, and kinematic GPS tracking. The reduced dynamic and kinematic approaches to orbit determination are found to have significant advantages with respect to reducing the level of geographically correlated orbit errors.

Rosborough, George↗

Natural marine cloud brightening in the Southern Ocean

Abstract. The number of cloud droplets per unit volume (Nd) is a fundamentally important property of marine boundary layer (MBL) liquid clouds that, at constant liquid water path, exerts considerable controls on albedo. Past work has shown that regional Nd has a direct correlation to marine primary productivity (PP) because of the role of seasonally varying, biogenically derived precursor gases in modulating secondary aerosol properties. These linkages are thought to be observable over the high-latitude oceans, where strong seasonal variability in aerosol and meteorology covary in mostly pristine environments. Here, we examine Nd variability derived from 5 years of MODIS Level 2-derived cloud properties in a broad region of the summer eastern Southern Ocean and adjacent marginal seas. We demonstrate latitudinal, longitudinal and temporal gradients in Nd that are strongly correlated with the passage of air masses over high-PP waters that are mostly concentrated along the Antarctic Shelf poleward of 60∘ S. We find that the albedo of MBL clouds in the latitudes south of 60∘ S is significantly higher than similar liquid water path (LWP) clouds north of this latitude.

54 ENVIRONMENTAL SCIENCES↗

Simple model to investigate jet quenching and correlated errors for centrality-dependent nuclear modification factors in relativistic heavy-ion collisions

Here, we apply Bayesian techniques to compare a simple, empirical model for jet quenching in heavy-ion collisions to centrality-dependent jet R AA measured by ATLAS for Pb + Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV. We find that the R AA values for central collisions are adequately described with a model for the mean p T -dependent jet energy loss using only two parameters. This model is extended by incorporating two-dimensional initial geometry information from TRENTo and compared to centrality-dependent R AA values. We find that the results are sensitive to the value of the jet-quenching formation time, τ ƒ , and that the optimal value of τ ƒ varies with the assumed path-length dependence of the energy loss. We construct a covariance error matrix for the data from the p T -dependent contributions to the ATLAS systematic errors and perform Bayesian calibrations for several different assumptions for the systematic error correlations. We show that the most-probable functions and $χ^2_d$ values are sensitive to assumptions made when fitting to correlated errors. This work demonstrates the utility of a simple model that can quickly demonstrate the constraining power of jet-quenching observables with corresponding uncertainties and guide future studies using more sophisticated models.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Adjusting solar-induced fluorescence to nadir-viewing provides a better proxy for GPP

Solar-induced fluorescence (SIF) provides key information for inferring terrestrial gross primary productivity (GPP). However, there is little research on analyzing the performance of nadir-adjusted SIF (SIF nadir ) on GPP, compared to total SIF emitted by all leaves (SIF total ) for reducing the viewing angle effects and estimating GPP. Besides, there have been controversial opinions on whether SIF total is better correlated to GPP than raw SIF observations (SIF obs ). Here we systematically analyzed the relationship between raw/nadir/total SIF (i.e., SIF obs , SIF nadir and SIF total ) and GPP and investigated the underlying mechanism, using multi-angular field measurements and eddy covariance data in wheat and corn crops at sub-daily scale. We further compared the performance of SIF obs , SIF nadir and SIF total in estimating GPP using the TROPOspheric Monitoring Instrument (TROPOMI) SIF and concurrent AmeriFlux measurements at daily scale. Furthermore, the results indicate that diurnal SIF nadir has stronger correlations to GPP than SIF obs for field measurements, with the increase of mean coefficient of determination (R 2 ) by 0.05–0.07 for far-red band and 0.11–0.20 for red band. SIF nadir shows comparable performance with SIF total for both far-red and red bands. Although the viewing angle effects dominate the difference in estimating GPP between SIF nadir and SIF obs , the correlation between light use efficiency (LUE) and fesc further determines the different performance of SIF nadir and SIF total in estimating GPP. TROPOMI-based analysis further confirms that SIF nadir overall has higher correlations to AmeriFlux GPP than SIF obs for different plant functional types and shows similar performance with SIF total . Compared to SIF total , the estimation of SIF nadir independent of viewing angle effects does not require any canopy structure parameters, and thus offers promising potential for reliably estimating regional and global GPP.

54 ENVIRONMENTAL SCIENCES↗

Chiral-like doublet band structure and octupole correlations in 104 Ag

The nature of the yrast negative-parity band and its chiral-like partner band in 104 Ag is investigated experimentally and theoretically. Lifetimes of states in the negative-parity yrast band and positive-parity band based on the 4424-keV level are measured using Doppler shift attenuation technique. Lifetimes of three more states have been determined along with the upper limit for the lifetime of the highest observed yrast states. Further, lifetimes known from earlier studies are determined with better precision. The level scheme of 104 Ag has also been extended with the addition of new enhanced E1 transitions linking the positive-parity band based on the 4424-keV levels and the yrast negative-parity and its partner band. B(E1) and/or B(E1)/B(M1) values for the transitions from the positive-parity band to the yrast and its partner band have been determined for the first time; these suggest strong octupole correlation between the positive-parity and the negative-parity bands. Here, calculations based on the triaxial projected shell model (TPSM) and covariant density functional theory (CDFT) have been performed to unravel the intrinsic structures of the partner band and the excited positive-parity band. TPSM calculations predict that doublet bands have significant angular momentum contributions along the three principle axes, suggesting that bands could have chiral symmetry breaking origin. The CDFT calculations predict a π(g 9/2 ) -1 Ⓧν(h 11/2 )(g 7/2 ,d 5/2 ) 2 aligned quasiparticle configuration for the negative-parity doublet bands with deformation parameters β≈0.20 and γ≈5°. The partner band could be interpreted as a chiral vibration mode built on top of the yrast band. The excited positive-parity band is predicted to have aligned four quasiparticle configurations, namely, π(g 9/2 ) -1 Ⓧν(h 11/2 ) 2 (g 7/2 ,d 5/2 ) 1 . Further, these calculations predict significant octupole softness in 104 Ag which could be the reason for enhanced E1 transitions between the four quasiparticle positive-parity bands and the doublet negative-parity bands.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Method of statistical filtering

Minimal formula for bounding the cross correlation between a random forcing function and the state error when this correlation is unknown is used in optimal linear filter theory applications. Use of the bound results in overestimation of the estimation-error covariance.

Battin, R. H.↗

Monitoring global geodynamical parameter variations with GPS and low earth satellites

A low earth satellite enhances global geodynamical parameters determination with GPS in two ways. First, it improves the GPS orbits, which in turn improve the estimates of other parameters. Secondly, a low earth satellite completes an orbit cycle in far shorter time (90 to 120 minutes) than do GPS satellites (12 hours); it observes more GPS satellites than a ground receiver does in shorter time and increases the correlation between GPS orbit errors. This reduces the error in the determination of nonrotational coordinate parameters, i.e., geocentric offset of ground tracking sites. Covariance results with the global geodynamical parameters modeled as constants, and as random-walk parameters to closer reflect the actual variations, are compared. The effects of using different GPS data quality and different ground tracking network are studied. The case of using two earth satellites in orthogonal orbital planes is also investigated.

Wu, S. C.↗

The Effects of Noise on Polarimetric SAR Data

Polarimetric SAR data can provide a great deal of information about the scattering behavior of the surface under observation. Polarimetric SAR systems often measure the scattering matrices of the areas under observation in linear polarizations (H and V). From the scattering matrix commonly used forms such as the covariance matrix and the Stokes matrix can be easily derived. Other measures derived from polarimetric SAR data include the standard deviation of texture, correlation coefficients between scattering matrix terms, and the mode and variance of phase differences between scattering matrix terms. The effects of additive system noise on these measurements is not often considered in the literature on this subject.

Polarimetric SAR Data analysis↗

The role of blocking circulation and emerging open water feedbacks on Greenland cold‐season air temperature variability over the last century

Substantial marine, terrestrial, and atmospheric changes have occurred over the Greenland region during the last century. Several studies have documented record‐levels of Greenland Ice Sheet (GrIS) summer melt extent during the 2000s and 2010s, but relatively little work has been carried out to assess regional climatic changes in other seasons. Here, we focus on the less studied cold‐season (i.e., autumn and winter) climate, tracing the long‐term (1873–2013) variability of Greenland's air temperatures through analyses of coastal observations and model‐derived outlet glacier series and their linkages with North Atlantic sea ice, sea surface temperature (SST), and atmospheric circulation indices. Through a statistical framework, large amounts of west and south Greenland temperature variance (up to r^2~50%) can be explained by the seasonally‐contemporaneous combination of the Greenland Blocking Index (GBI) and the North Atlantic Oscillation (NAO; hereafter the combination of GBI and NAO is termed GBI). Lagged and concomitant regional sea‐ice concentration (SIC) and the Atlantic Multidecadal Oscillation (AMO) seasonal indices account for small amounts of residual air temperature variance (r^2<~10%) relative to the GBI. The correlations between GBI and cold‐season temperatures are predominantly positive and statistically‐significant through time, while regional SIC conditions emerge as a significant covariate from the mid‐20th century through the conclusion of the study period. The inclusion of the cold‐season Pacific Decadal Oscillation (PDO) in multivariate analyses bolsters the air temperature variance explained by the North Atlantic regional predictors, suggesting the remote, background climate state is important to long‐term Greenland temperature variability. These findings imply that large‐scale tropospheric circulation has a strong control on surface temperature over Greenland through dynamic and thermodynamic impacts and stress the importance of understanding the evolving two‐way linkages between the North Atlantic marine and atmospheric environment in order to more accurately predict Greenland seasonal climate variability and change through the 21st century.

air temperatures↗

On the Use of the Log-Normal Particle Size Distribution to Characterize Global Rain

Although most parameterizations of the drop size distributions (DSD) use the gamma function, there are several advantages to the log-normal form, particularly if we want to characterize the large scale space-time variability of the DSD and rain rate. The advantages of the distribution are twofold: the logarithm of any moment can be expressed as a linear combination of the individual parameters of the distribution; the parameters of the distribution are approximately normally distributed. Since all radar and rainfall-related parameters can be written approximately as a moment of the DSD, the first property allows us to express the logarithm of any radar/rainfall variable as a linear combination of the individual DSD parameters. Another consequence is that any power law relationship between rain rate, reflectivity factor, specific attenuation or water content can be expressed in terms of the covariance matrix of the DSD parameters. The joint-normal property of the DSD parameters has applications to the description of the space-time variation of rainfall in the sense that any radar-rainfall quantity can be specified by the covariance matrix associated with the DSD parameters at two arbitrary space-time points. As such, the parameterization provides a means by which we can use the spaceborne radar-derived DSD parameters to specify in part the covariance matrices globally. However, since satellite observations have coarse temporal sampling, the specification of the temporal covariance must be derived from ancillary measurements and models. Work is presently underway to determine whether the use of instantaneous rain rate data from the TRMM Precipitation Radar can provide good estimates of the spatial correlation in rain rate from data collected in 5(sup 0)x 5(sup 0) x 1 month space-time boxes. To characterize the temporal characteristics of the DSD parameters, disdrometer data are being used from the Wallops Flight Facility site where as many as 4 disdrometers have been used to acquire data over a 2 km path. These data should help quantify the temporal form of the covariance matrix at this site.

Meneghini, Robert↗

Cross-scale covariance for material property prediction

A simulation can stand its ground against an experiment only if its prediction uncertainty is known. The unknown accuracy of interatomic potentials (IPs) is a major source of prediction uncertainty, severely limiting the use of large-scale classical atomistic simulations in a wide range of scientific and engineering applications. Here we explore covariance between predictions of metal plasticity, from 178 large-scale (~10 8 atoms) molecular dynamics (MD) simulations, and a variety of indicator properties computed at small-scales (≤10 2 atoms). All simulations use the same 178 IPs. In a manner similar to statistical studies in public health, we analyze correlations of strength with indicators, identify the best predictor properties, and build a cross-scale “strength-on-predictors” regression model. This model is then used to estimate regression error over the statistical pool of IPs. Small-scale predictors found to be highly covariant with strength are computed using expensive quantum-accurate calculations and used to predict flow strength, within the statistical error bounds established in our study.

36 MATERIALS SCIENCE↗