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 1,027 records · Page 57
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

A high temperature instrument for consecutive measurements of thermal conductivity, electrical conductivity, and Seebeck coefficient

A device for measuring a plurality of material properties is designed to include accurate sensors configured to consecutively obtain thermal conductivity, electrical conductivity, and Seebeck coefficient of a single sample while maintaining a vacuum or inert gas environment. Four major design factors are confirmed as sample-heat spreader mismatch, radiation losses, parasitic losses, and sample surface temperature variance. The design is analyzed using finite element methods for high temperature ranges up to 1000°C as well as ultra-high temperatures up to 2500°C. A temperature uncertainty of 0.46% was estimated for a sample with cold and hot sides at 905.1 and 908.5°C, respectively. The uncertainty at 1000°C was calculated to be 0.7% for a ΔT of 5°C between the hot and cold sides. The thermal conductivity uncertainty was calculated to be -8.6% at ~900°C for a case with radiative gains, and +8.2% at ~1000°C for a case with radiative losses, indicating the sensitivity of the measurement to the temperature of the thermal guard in relation to the heat spreader and sample temperature. Lower limits of -17 and -13% error in thermal conductivity measurements were estimated at the ultra-high temperature of ~2500°C for a single-stage and double-stage radiation shield, respectively. Lastly, it is noted that this design is not limited to electro-thermal characterization and will enable measurement of ionic conductivity and surface temperatures of energy materials under realistic operating conditions in extreme temperature environments.

47 OTHER INSTRUMENTATION↗
NASA NTRS2006

Estimating Random Errors Due to Shot Noise in Backscatter Lidar Observations

In this paper, we discuss the estimation of random errors due to shot noise in backscatter lidar observations that use either photomultiplier tube (PMT) or avalanche photodiode (APD) detectors. The statistical characteristics of photodetection are reviewed, and photon count distributions of solar background signals and laser backscatter signals are examined using airborne lidar observations at 532 nm using a photon-counting mode APD. Both distributions appear to be Poisson, indicating that the arrival at the photodetector of photons for these signals is a Poisson stochastic process. For Poisson-distributed signals, a proportional, one-to-one relationship is known to exist between the mean of a distribution and its variance. Although the multiplied photocurrent no longer follows a strict Poisson distribution in analog-mode APD and PMT detectors, the proportionality still exists between the mean and the variance of the multiplied photocurrent. We make use of this relationship by introducing the noise scale factor (NSF), which quantifies the constant of proportionality that exists between the root-mean-square of the random noise in a measurement and the square root of the mean signal. Using the NSF to estimate random errors in lidar measurements due to shot noise provides a significant advantage over the conventional error estimation techniques, in that with the NSF uncertainties can be reliably calculated from/for a single data sample. Methods for evaluating the NSF are presented. Algorithms to compute the NSF are developed for the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) lidar and tested using data from the Lidar In-space Technology Experiment (LITE). OCIS Codes:

Liu, Zhaoyan↗
DOE OSTI2025

Signal-preserving CMB component separation with machine learning

Analysis of microwave sky signals, such as the cosmic microwave background, often requires component separation using multifrequency methods, whereby different signals are isolated according to their different frequency behaviors. Many so-called blind methods, such as the internal linear combination (ILC), make minimal assumptions about the spatial distribution of the signal or contaminants, and only assume knowledge of the frequency dependence of the signal. The ILC produces a minimum-variance linear combination of the measured frequency maps. In the case of Gaussian, statistically isotropic fields, this is the optimal linear combination, as the variance is the only statistic of interest. However, in many cases the signal we wish to isolate, or the foregrounds we wish to remove, are non-Gaussian and/or statistically anisotropic (in particular for the case of Galactic foregrounds). In such cases, it is possible that machine learning (ML) techniques can be used to exploit the non-Gaussian features of the foregrounds and thereby improve component separation. However, many ML techniques require the use of complex, difficult-to-interpret operations on the data. We propose a hybrid method whereby we train an ML model using only combinations of the data that , and combine the resulting ML-predicted foreground estimate with the ILC solution to reduce the error from the ILC. We demonstrate our methods on simulations of extragalactic temperature and Galactic polarization foregrounds and show that our ML model can exploit non-Gaussian features, such as point sources and spatially varying spectral indices, to produce lower-variance maps than ILC—e.g., reducing the variance of the B-mode residual by factors of up to 5—while preserving the signal of interest in an unbiased manner. Moreover, we often find improved performance even when applying our ML technique to foreground models on which it was not trained. Published by the American Physical Society 2025

McCarthy, Fiona (ORCID:0000000253893565)↗
DOE OSTI2015

A model and variance reduction method for computing statistical outputs of stochastic elliptic partial differential equations

We present a model and variance reduction method for the fast and reliable computation of statistical outputs of stochastic elliptic partial differential equations. Our method consists of three main ingredients: (1) the hybridizable discontinuous Galerkin (HDG) discretization of elliptic partial differential equations (PDEs), which allows us to obtain high-order accurate solutions of the governing PDE; (2) the reduced basis method for a new HDG discretization of the underlying PDE to enable real-time solution of the parameterized PDE in the presence of stochastic parameters; and (3) a multilevel variance reduction method that exploits the statistical correlation among the different reduced basis approximations and the high-fidelity HDG discretization to accelerate the convergence of the Monte Carlo simulations. The multilevel variance reduction method provides efficient computation of the statistical outputs by shifting most of the computational burden from the high-fidelity HDG approximation to the reduced basis approximations. Furthermore, we develop a posteriori error estimates for our approximations of the statistical outputs. Based on these error estimates, we propose an algorithm for optimally choosing both the dimensions of the reduced basis approximations and the sizes of Monte Carlo samples to achieve a given error tolerance. We provide numerical examples to demonstrate the performance of the proposed method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

UV luminosity functions at Redshifts z ∼ 4 TO z ∼ 10: 10,000 galaxies from HST legacy fields

The remarkable Hubble Space Telescope (HST) data sets from the CANDELS, HUDF09, HUDF12, ERS, and BoRG/HIPPIES programs have allowed us to map the evolution of the rest-frame UV luminosity function (LF) from z∼10 to z∼4. We develop new color criteria that more optimally utilize the full wavelength coverage from the optical, near-IR, and mid-IR observations over our search fields, while simultaneously minimizing the incompleteness and eliminating redshift gaps. We have identified 5859, 3001, 857, 481, 217, and 6 galaxy candidates at z∼4, z∼5, z∼6, z∼7, z∼8, and z∼10, respectively, from the ∼1000 arcmin{sup 2} area covered by these data sets. This sample of >10,000 galaxy candidates at z⩾4 is by far the largest assembled to date with HST. The selection of z ∼ 4–8 candidates over the five CANDELS fields allows us to assess the cosmic variance; the largest variations are at z⩾7. Our new LF determinations at z∼4 and z∼5 span a 6 mag baseline and reach to –16 AB mag. These determinations agree well with previous estimates, but the larger samples and volumes probed here result in a more reliable sampling of >L{sup ∗} galaxies and allow us to reassess the form of the UV LFs. Our new LF results strengthen our earlier findings to 3.4σ significance for a steeper faint-end slope of the UV LF at z>4, with α evolving from α=−1.64±0.04 at z∼4 to α=−2.06±0.13 at z∼7 (and α=−2.02±0.23 at z∼8), consistent with that expected from the evolution of the halo mass function. We find less evolution in the characteristic magnitude M{sup *} from z∼7 to z∼4; the observed evolution in the LF is now largely represented by changes in ϕ{sup ∗}. No evidence for a non-Schechter-like form to the z ∼ 4–8 LFs is found. A simple conditional LF model based on halo growth and evolution in the M/L ratio (∝(1+z){sup −1.5}) of halos provides a good representation of the observed evolution.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1989

Physically-based parameterization of spatially variable soil and vegetation using satellite multispectral data

A stochastic-geometric landsurface reflectance model is formulated and tested for the parameterization of spatially variable vegetation and soil at subpixel scales using satellite multispectral images without ground truth. Landscapes are conceptualized as 3-D Lambertian reflecting surfaces consisting of plant canopies, represented by solid geometric figures, superposed on a flat soil background. A computer simulation program is developed to investigate image characteristics at various spatial aggregations representative of satellite observational scales, or pixels. The evolution of the shape and structure of the red-infrared space, or scattergram, of typical semivegetated scenes is investigated by sequentially introducing model variables into the simulation. The analytical moments of the total pixel reflectance, including the mean, variance, spatial covariance, and cross-spectral covariance, are derived in terms of the moments of the individual fractional cover and reflectance components. The moments are applied to the solution of the inverse problem: The estimation of subpixel landscape properties on a pixel-by-pixel basis, given only one multispectral image and limited assumptions on the structure of the landscape. The landsurface reflectance model and inversion technique are tested using actual aerial radiometric data collected over regularly spaced pecan trees, and using both aerial and LANDSAT Thematic Mapper data obtained over discontinuous, randomly spaced conifer canopies in a natural forested watershed. Different amounts of solar backscattered diffuse radiation are assumed and the sensitivity of the estimated landsurface parameters to those amounts is examined.

Jasinski, Michael F.↗
NASA NTRS1989

Pyroxene thermobarometry in LL-group chondrites and implications for parent body metamorphism

Geothermometry based on the compositions of clinopyroxenes in type 6 and 7 LL chondrites gives coherent results, but the estimated temperatures from coexisting orthopyroxenes are consistently lower than for clinopyroxenes. Orthopyroxene thermometry is suspect because of compositional effects of polymorphic inversions and/or unknown kinetic factors. Lack of clinopyroxene equilibration precludes accurate estimation of peak metamorphic temperatures for type 4 and 5 chondrites. There is no apparent correlation between Al content (a pressure-dependent variable) and equilibration temperature in chondritic pyroxenes. This finding, which is at variance with a previously published conclusion that temperature and pressure are correlated in metamorphosed chondrites, may have important implications for asteroid thermal models.

Mcsween, Harry Y., Jr.↗
NASA NTRS2001

The AutoBayes Program Synthesis System: System Description

AUTOBAYES is a fully automatic program synthesis system for the statistical data analysis domain. Its input is a concise description of a data analysis problem in the form of a statistical model; its output is optimized and fully documented C/C++ code which can be linked dynamically into the Matlab and Octave environments. AUTOBAYES synthesizes code by a schema-guided deductive process. Schemas (i.e., code templates with associated semantic constraints) are applied to the original problem and recursively to emerging subproblems. AUTOBAYES complements this approach by symbolic computation to derive closed-form solutions whenever possible. In this paper, we concentrate on the interaction between the symbolic computations and the deductive synthesis process. A statistical model specifies for each problem variable (i.e., data or parameter) its properties and dependencies in the form of a probability distribution, A typical data analysis task is to estimate the best possible parameter values from the given observations or measurements. The following example models normal-distributed data but takes prior information (e.g., from previous experiments) on the data's mean value and variance into account.

Fischer, Bernd↗
NASA NTRS2006

Tropical Rainfall Variability on Interannual-to-Interdecadal/Longer-Time Scales Derived from the GPCP Monthly Product

Global and large regional rainfall variations and possible long-term changes are examined using the 26-year (1979-2004) GPCP monthly dataset (Adler et al., 2003). Our emphasis is to discriminate among variations due to ENSO, volcanic events, and possible long-term climate changes in the tropics. Although the global linear change of precipitation in the data set is near zero during the time period, an increase in tropical rainfall is noted, with a weaker decrease over northern hemisphere middle latitudes. Focusing on the tropics (25degS-25degN), the data set indicates an upward trend (0.06 mm/day/decade) and a downward trend (-0.02 mm/day/decade) over tropical ocean and land, respectively. This corresponds to an about 4.9% increase (ocean) and 1.6% decrease (land) during the entire 26-year time period. Techniques are applied to isolate and quantify variations due to ENSO and two major volcanic eruptions (El Chichon, March 1982; Pinatubo, June 1991) in order to examine longer time-scale changes. The ENSO events generally do not impact the tropical total rainfall, but, of course, induce significant anomalies with opposite signs over tropical land and ocean. The impact of the two volcanic eruptions is estimated to be about a 5% reduction in tropical rainfall over both land and ocean. A modified data set (with ENSO and volcano effects removed) retains the same approximate linear change slopes, but with reduced variance, thereby increasing the confidence levels associated with the long-term rainfall changes in the tropics 2

Gu, Guojun↗
NASA NTRS2021

The role of analysis error in the convergence of reanalysis production streams in MERRA-2

Due to production time constraints, most reanalyses are produced in multiple parallel streams instead of a single continuous one. These streams cover separate segments of the reanalysis time period with short overlaps to allow reconstruction of the official record. A fundamental assumption justifying this approach is that the streams will be assimilating the same observations during the periods where they overlap, and so will eventually converge to a similar atmospheric state, making discontinuities at stream junctions negligible. This assumption is revisited in this work by examining the impact of analysis error on the differences between MERRA-2 overlapping streams in three historical periods. Comparison results are shown in terms of standard deviations of stream differences as well as the spectral decomposition of the variance of their differences. Residual differences were found at the end of each year of overlap, with larger values observed in the earlier segments of the presatellite era. By drawing parallels with analysis error statistics estimated from the GMAO OSSE system, these differences are shown to reflect the varying constraint of data with the varying observing network, and to further carry the imprint of errors that the data assimilation process is not able to mitigate. As such, they are unlikely to be reduced by longer spinup periods. The ability of data assimilation to ensure continuity in the parallel streams is put into question when the observing system coverage is inadequate or simply when the data assimilation system as a whole is suboptimal.

Amal El Akkraoui↗
NASA NTRS2019

Assessing Uncertainties of a Geophysical Approach to Estimate Surface Fine Particulate Matter Distributions from Satellite-Observed Aerosol Optical Depth

Health impact analyses are increasingly tapping the broad spatial coverage of satellite aerosol optical depth (AOD) products to estimate human exposure to fine particulate matter (PM2.5). We use a forward geophysical approach to derive ground-level PM2.5 distributions from satellite AOD at 1 km2(exp) resolution for 2011 over the northeastern US by applying relationships between surface PM2.5 and column AOD (calculated offline from speciated mass distributions) from a regional air quality model (CMAQ; 12×12 km2(exp) horizontal resolution). Seasonal average satellite-derived PM2.5 reveals more spatial detail and best captures observed surface PM2.5 levels during summer. At the daily scale, however, satellite-derived PM2.5 is not only subject to measurement uncertainties from satellite instruments, but more importantly to uncertainties in the relationship between surface PM2.5 and column AOD. Using 11 ground-based AOD measurements within 10 km of surface PM2.5 monitors, we show that uncertainties in modeled PM2.5∕AOD can explain more than 70 % of the spatial and temporal variance in the total uncertainty in daily satellite-derived PM2.5 evaluated at PM2.5 monitors. This finding implies that a successful geophysical approach to deriving daily PM2.5 from satellite AOD requires model skill at capturing day-to-day variations in PM2.5∕AOD relationships. Overall, we estimate that uncertainties in the modeled PM2.5∕AOD lead to an error of 11 µg m−3(exp) in daily satellite-derived PM2.5, and uncertainties in satellite AOD lead to an error of 8 µg m−3(exp). Using multi-platform ground, airborne, and radiosonde measurements, we show that uncertainties of modeled PM2.5∕AOD are mainly driven by model uncertainties in aerosol column mass and speciation, while model representation of relative humidity and aerosol vertical profile shape contributes some systematic biases. The parameterization of aerosol optical properties, which determines the mass extinction efficiency, also contributes to random uncertainty, with the size distribution being the largest source of uncertainty and hygroscopicity of inorganic salt the second largest. Future efforts to reduce uncertainty in geophysical approaches to derive surface PM2.5 from satellite AOD would thus benefit from improving model representation of aerosol vertical distribution and aerosol optical properties, to narrow uncertainty in satellite-derived PM2.5.

aerosol optical depth↗
DOE OSTI2014

Processing of X-ray snapshots from crystals in random orientations

A new method for the treatment of partial reflections from X-ray snapshots is implemented in the program package nXDS, which yields intensity data of almost the same quality as those obtained by the classical rotation method. A functional expression is introduced that relates scattered X-ray intensities from a still or a rotation snapshot to the corresponding structure-factor amplitudes. The new approach was implemented in the program nXDS for processing monochromatic diffraction images recorded by a multi-segment detector where each exposure could come from a different crystal. For images containing indexable spots, the intensities of the expected reflections and their variances are obtained by profile fitting after mapping the contributing pixel contents to the Ewald sphere. The varying intensity decline owing to the angular distance of the reflection from the surface of the Ewald sphere is estimated using a Gaussian rocking curve. This decline is dubbed ‘Ewald offset correction’, which is well defined even for still images. Together with an image-scaling factor and other corrections, an explicit expression is defined that predicts each recorded intensity from its corresponding structure-factor amplitude. All diffraction parameters, scaling and correction factors are improved by post-refinement. The ambiguous case of a lower point group than the lattice symmetry is resolved by a method reminiscent of the technique of ‘selective breeding’. It selects the indexing alternative for each image that yields, on average, the highest correlation with intensities from all other images. Processing a test set of rotation images by XDS and treating the same images by nXDS as snapshots of crystals in random orientations yields data of comparable quality, clearly indicating an anomalous signal from Se atoms.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1979

Microprocessor realizations of range rate filters

The performance of five digital range rate filters is evaluated. A range rate filter receives an input of range data from a radar unit and produces an output of smoothed range data and its estimated derivative range rate. The filters are compared through simulation on an IBM 370. Two of the filter designs are implemented on a 6800 microprocessor-based system. Comparisons are made on the bases of noise variance reduction ratios and convergence times of the filters in response to simulated range signals.

Source record↗
NASA NTRS1979

Geomagnetopause surface fluctuations observed by Voyager 1

The paper discusses the geomagnetopause surface fluctuations observed by Voyager 1. Normals to the magnetopause were determined for the crossings of Voyager 1 by minimum variance analysis of the internal magnetic field; the oscillating nature of the ecliptic plane component of these normals indicates that the multiple crossings were due to a wavelike surface disturbance moving tailward along the magnetopause. The estimated amplitude of these waves was small compared to their wavelengths; this conclusion is independent of any bulk normal motion of the magnetopause.

Lepping, R. P.↗
NASA NTRS1988

An improved error assessment for the GEM-T1 gravitational model

Several tests were designed to determine the correct error variances for the Goddard Earth Model (GEM)-T1 gravitational solution which was derived exclusively from satellite tracking data. The basic method employs both wholly independent and dependent subset data solutions and produces a full field coefficient estimate of the model uncertainties. The GEM-T1 errors were further analyzed using a method based upon eigenvalue-eigenvector analysis which calibrates the entire covariance matrix. Dependent satellite and independent altimetric and surface gravity data sets, as well as independent satellite deep resonance information, confirm essentially the same error assessment. These calibrations (utilizing each of the major data subsets within the solution) yield very stable calibration factors which vary by approximately 10 percent over the range of tests employed. Measurements of gravity anomalies obtained from altimetry were also used directly as observations to show that GEM-T1 is calibrated. The mathematical representation of the covariance error in the presence of unmodeled systematic error effects in the data is analyzed and an optimum weighting technique is developed for these conditions. This technique yields an internal self-calibration of the error model, a process which GEM-T1 is shown to approximate.

Lerch, F. J.↗
NASA NTRS1989

An improved error assessment for the GEM-T1 gravitational model

Several tests were designed to determine the correct error variances for the GEM-T1 gravitational solution which was derived exclusively from satellite tracking data. The basic method employs both wholly independent and dependent subset data solutions and produces a full field coefficient by coefficient estimate of the model uncertainties. The GEM-T1 errors were further analyzed using a method based upon eigenvalue-eigenvector analysis which calibrates the entire covariance matrix. Dependent satellite and independent altimetric and surface gravity data sets, as well as independent satellite deep resonance information, confirm essentially the same error assessment.

Lerch, F. J.↗
DOE OSTI2019

Uncovering the Forced Climate Response from a Single Ensemble Member Using Statistical Learning

Internal atmospheric variability fundamentally limits predictability of climate and obscures evidence of anthropogenic climate change regionally and on time scales of up to a few decades. Dynamical adjustment techniques estimate and subsequently remove the influence of atmospheric circulation variability on temperature or precipitation. The residual component is expected to contain the thermodynamical signal of the externally forced response but with less circulation-induced noise. Existing techniques have led to important insights into recent trends in regional (hydro-) climate and their drivers, but the variance explained by circulation is often low. Here, we develop a novel dynamical adjustment technique by implementing principles from statistical learning. We demonstrate in an ensemble of Community Earth System Model (CESM) simulations that statistical learning methods, such as regularized linear models, establish a clearer relationship between circulation variability and atmospheric target variables, and need relatively short periods of record for training (around 30 years). The method accounts for, on average, 83% and 78% of European monthly winter temperature and precipitation variability at gridcell level, and around 80% of global mean temperature and hemispheric precipitation variability. We show that the residuals retain forced thermodynamical contributions to temperature and precipitation variability. Accurate estimates of the total forced response can thus be recovered assuming that forced circulation changes are gradual over time. Overall, forced climate response estimates can be extracted at regional or global scales from approximately 3–5 times fewer ensemble members, or even a single realization, using statistical learning techniques. We anticipate the technique will contribute to reducing uncertainties around internal variability and facilitating climate change detection and attribution.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

A variance-based decomposition and global sensitivity index method for uncertainty quantification: Application to retrieved ice cloud properties

This study develops a novel uncertainty quantification (UQ) method for cloud microphysical property retrievals using variance-based decomposition and global sensitivity index. In this UQ framework, empirical orthogonal function (EOF) analysis is applied to the U.S. Department of Energy Atmospheric Radiation Measurement (ARM) ground-based observations, which are the inputs for the cloud retrieval studied here. The principal components (PCs) in the EOF expansion are parameterized as random input variables, and hence, the input dimension is greatly reduced (up to a factor of 50), allowing large ensemble of random samplings. The EOF expansion improves the accuracy of the uncertainty estimation by taking into account the cross correlations in the input data profiles. Here, this method enables a probabilistic representation of a retrieval process by adding normally distributed perturbations into PCs of sample means of input data profiles within a time window. Therefore, it effectively facilitates objective validation of climate models against cloud retrievals under a probabilistic framework for rigorous statistical inferences. Moreover, the variance-based global sensitivity index analysis, part of this method, attributes the output uncertainties to each individual source, thus providing directions for improving retrieval algorithms and observation strategies. For demonstration, we apply this method to quantify the uncertainties of the ARM program's baseline cloud retrieval algorithm for an ice cloud case observed at the Southern Great Plains site on 9 March 2000.

54 ENVIRONMENTAL SCIENCES↗