Search NASA⌕ Search

SEARCH · Search NASA

Results for “sampling statistics”

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 631 records · Page 35

Entropy and its Relationship with Statistics

The purpose of our report is to discuss the notion of entropy and its relationship with statistics. Our goal is to provide a manner in which you can think about entropy, its central role within information theory and relationship with statistics. We review various relationships between information theory and statistics—nearly all are well-known but unfortunately are often not recognized. Entropy quantities the "average amount of surprise" in a random variable and lies at the heart of information theory, which studies the transmission, processing, extraction, and utilization of information. For us, data is information. What is the distinction between information theory and statistics? Information theorists work with probability distributions. Instead, statisticians work with samples. In so many words, information theory using samples is the practice of statistics.

97 MATHEMATICS AND COMPUTING↗

Sampling Functions from Gaussian Processes and Structured Covariance Gaussian Networks

When learning aerodynamic models from data, it is critical to incorporate estimates of model uncertainty. This motivates the design of probabilistic aerodynamic databases which can be sampled to generate physically and statistically plausible aerodynamic models. In this talk we discuss how to sample deterministic functions from two different kinds of probabilistic models and demonstrate their use. First, Gaussian Process Regressors (GPRs) are a widely used probabilistic kernel-based model which can be thought of as Gaussian distributions over functions. GPRs are generally trained by maximizing the marginal likelihood of seeing the training data over the kernel parameter space. Sample functions are easily generated by drawing points from the Gaussian distribution at desired input points. However, when the points are not known ahead of time, the classical sampling approach is not possible since successive function samples will generate different function realizations. We present an approach for sampling consistent function evaluations from a GPR over multiple samples. Second, we describe a neural network architecture which learns a conditional Gaussian distribution by maximizing the marginal likelihood at each point in the input space. We then discuss and compare several options for generating sample functions which match this distribution. Finally, we demonstrate the use of these probabilistic aerodynamic models in an atmospheric reentry simulation.

Gaussian process regression↗

The Complete Light-curve Sample of Spectroscopically Confirmed SNe Ia from Pan-STARRS1 and Cosmological Constraints from the Combined Pantheon Sample

Here, we present optical light curves, redshifts, and classifications for 365 spectroscopically confirmed Type Ia supernovae (SNe Ia) discovered by the Pan-STARRS1 (PS1) Medium Deep Survey. We detail improvements to the PS1 SN photometry, astrometry, and calibration that reduce the systematic uncertainties in the PS1 SN Ia distances. We combine the subset of 279 PS1 SNe Ia (0.03 < z < 0.68) with useful distance estimates of SNe Ia from the Sloan Digital Sky Survey (SDSS), SNLS, and various low-z and Hubble Space Telescope samples to form the largest combined sample of SNe Ia, consisting of a total of 1048 SNe Ia in the range of 0.01 < z < 2.3, which we call the Pantheon Sample. When combining Planck 2015 cosmic microwave background (CMB) measurements with the Pantheon SN sample, we find Ω m = 0.307 ± 0.012 and w = -1.026 ± 0.041 for the wCDM model. When the SN and CMB constraints are combined with constraints from BAO and local H 0 measurements, the analysis yields the most precise measurement of dark energy to date: w 0 = -1.007 ± 0.089 and w α = -0.222 ± 0.407 for the w 0 w α CDM model. Tension with a cosmological constant previously seen in an analysis of PS1 and low-z SNe has diminished after an increase of 2× in the statistics of the PS1 sample, improved calibration and photometry, and stricter light-curve quality cuts. We find that the systematic uncertainties in our measurements of dark energy are almost as large as the statistical uncertainties, primarily due to limitations of modeling the low-redshift sample. This must be addressed for future progress in using SNe Ia to measure dark energy.

79 ASTRONOMY AND ASTROPHYSICS↗

Evaluating cosmological biases using photometric redshifts for Type Ia Supernova cosmology with the Dark Energy Survey Supernova Program

Cosmological analyses with Type Ia Supernovae (SNe Ia) have traditionally been reliant on spectroscopy for both classifying the type of supernova and obtaining reliable redshifts to measure the distance–redshift relation. While obtaining a host-galaxy spectroscopic redshift for most SNe is feasible for small-area transient surveys, it will be too resource intensive for upcoming large-area surveys such as the Vera Rubin Observatory Legacy Survey of Space and Time, which will observe on the order of millions of SNe. Here, we use data from the Dark Energy Survey (DES) to address this problem with photometric redshifts (photo-z) inferred directly from the SN light curve in combination with Gaussian and full p(z) priors from host-galaxy photo-z estimates. Using the DES 5-yr photometrically classified SN sample, we consider several photo-z algorithms as host-galaxy photo-z priors, including the Self-Organizing Map redshifts (SOMPZ), Bayesian Photometric Redshifts (BPZ), and Directional-Neighbourhood Fitting (DNF) redshift estimates employed in the DES 3 × 2 point analyses. With detailed catalogue-level simulations of the DES 5-yr sample, we find that the simulated w can be recovered within ±0.02 when using SN+SOMPZ or DNF prior photo-z, smaller than the average statistical uncertainty for these samples of 0.03. With data, we obtain biases in w consistent with simulations within ~1σ for three of the five photo-z variants. We further evaluate how photo-z systematics interplay with photometric classification and find classification introduces a subdominant systematic component. This work lays the foundation for next-generation fully photometric SNe Ia cosmological analyses.

(cosmology:) dark energy↗

AGN Activity in Nucleated Galaxies as Measured by Chandra

Motivated by theoretical expectations that nuclear star clusters (NSCs) in galactic centers may provide a favorable environment for supermassive black holes to form and/or efficiently grow, we set out to measure the fraction of nearby nucleated galaxies that also host an active galactic nucleus. We targeted a distance-limited sample of 98 objects with the Chandra X-ray Telescope, down to a uniform X-ray luminosity threshold of ∼10{sup 38} erg s{sup −1}. The sample is composed of 47 late-types and 51 early-types, enabling us to further investigate the active fraction as a function of galactic morphology. After correcting for contamination to the nuclear X-ray signal from bright X-ray binaries, we measure an active fraction f=11.2%{sub −4.9}{sup +7.4} (1σ C.L.) across the whole sample, in agreement with previous estimates based on a heterogeneous combination of optical, X-ray, and radio diagnostics, by Seth et al. After accounting for the different stellar mass distributions in our samples, we find no statistically significant difference in the active fraction of early- versus late-type nucleated galaxies, with f=10.6%{sub −4.9}{sup +11.9} and 10.8%{sub −6.3}{sup +11.3}, respectively. For the early-type nucleated galaxies, we are able to carry out a controlled comparison with a parent sample of non-nucleated galaxies covering the same stellar mass range, again finding no statistically significant difference in the active fraction. Taken at face value, our findings suggest that the presence of an NSC does not facilitate or enhance accretion-powered emission from a nuclear supermassive black hole. This is true even for late-type nucleated galaxies, home to bluer NSCs and arguably larger gas reservoirs.

79 ASTRONOMY AND ASTROPHYSICS↗

An Analysis of the Statistics and Systematics of Limb Anomaly Detections in HST/STIS Transit Images of Europa

Several recent studies derived the existence of plumes on Jupiter’s moon Europa. The only technique that provided multiple detections is the far-ultraviolet imaging observations of Europa in transit of Jupiter taken by the Space Telescope Imaging Spectrograph (STIS) on the Hubble Space Telescope (HST). In this study, we reanalyze the three HST/STIS transit images in which Sparks et al. identified limb anomalies as evidence for Europa’s plume activity. After reproducing the results of Sparks et al., we find that positive outliers are similarly present in the images as the negative outliers that were attributed to plume absorption. A physical explanation for the positive outliers is missing. We then investigate the systematic uncertainties and statistics in the images and identify two factors that are crucial when searching for anomalies around the limb. One factor is the alignment between the actual and assumed locations of Europa on the detector. A misalignment introduces distorted statistics, most strongly affecting the limb above the darker trailing hemisphere where the plumes were detected. The second factor is a discrepancy between the observation and the model used for comparison, adding uncertainty in the statistics. When accounting for these two factors, the limb minima (and maxima) are consistent with random statistical occurrence in a sample size given by the number of pixels in the analyzed limb region. The plume candidate features in the three analyzed images can be explained by purely statistical fluctuations and do not provide evidence for absorption by plumes.

79 ASTRONOMY AND ASTROPHYSICS↗

An IUE survey of interstellar H I Ly alpha absorption. 2: Interpretations

We present an analysis of interstellar neutral hydrogen column densities obtained from archival Ly alpha absorption line data toward the complete sample of B2 and hotter stars observed at high spectral resolution with the IUE satellite. The full sample includes 554 stars and more than doubles the number of lines of sight previously observed. Many of the B2 and B1.5 stars exhibit contamination from stellar Ly alpha absorption, and were excluded from the subsequent analysis. The final working sample includes 393 stars. We present statistical averages for a number of quantities derived for the sample. The stars range in distance from 0.12 to 11 kpc with an average distance of 2.1 kpc. Values for average sight-line density range from 0.017 to 8.62 atoms cm(exp -3), with an overall average of 0.23 atoms cm(exp -3). Neutral hydrogen and dust are well correlated, with N(H I)/E(B - V) = 4.93 x 10(exp 21) cm(exp -2) mag(exp -1), and N(H I)/E(Bump) = 2.17 x 10(exp 21) cm(exp -2) mag(exp -1). There is an increase in the H I to dust ratios for the densest sight lines. Evidently, the H I to dust ratios in dense clouds are higher than in the lower density medium between the clouds and in interarm directions. The effect is probably an indication of dust modification in the dense clouds. The H I to dust ratios are approx. 17% smaller for the sight lines to O stars compared to the ratios found toward B stars. We have identified sight lines with anomalously large and small values of N(H I)/E(B - V) and H(H I)/E(Bump). Some of these sight lines are well-known examples of peculiar extinction (i.e., HD 147933 and HD 37061 (NU Ori)). Others are less well known and may provide important new examples of regions with highly modified dust. The number of stars with anomalously low H I to dust ratios is disproportionately larger than the number of stars with anomalously high H I to dust ratios. The distribution of the gas away form the Galactic plane for the sample of objects is complex. For the subset of 375 stars whose lines of sight have a statistically low probablility of intersecting a large cloud, the distribution may be roughly approximated by an exponential density distribution, with a midplane density of 0.366 atoms cm(exp -3), a scale height of 195 pc, and a random logarithmic scattering parameter sigma (sub p) of 0.159 dex. For the same sample of stars, interstellar dust is found to be confined closer to the plane than the H I, with (E(B - V)/r)(sub 0) = 0.257 mag kpc(exp -1), scale height h = 152 pc, and sigma(sub p) = 0.263 dex. The dust distribution is more inhomogeneous than the gas distribution. The observed distribution of H I away from the Galactic plane can also be fitted with a two-component exponential model including the effects of patchiness. Unfortunately, the foreground contamination produced by the cloudy and more confined component of H I makes it difficult to estimate the scale height of the extended component without invoking an independent estimate of its midplane density. Assuming n(sub 2)(0) = 0.16 cm(exp -3) from the Copernicus satellite H I/H(sub 2) survey, we obtain n(sub 1)(0) = 0.247 cm(exp -3), h(sub 1) = 73 PC, and h(sub 2) + 357 pc.

Diplas, Athanassios↗

Statistics of associations among IR galaxies

In the course of expanding the search of Kleinmann et. al. (1988) for distant, infrared-luminous objects, the authors noticed (as is often remarked) that a large number of infrared-selected galaxies have close neighbors or show merger characteristics (e.g., tidal tails, distorted disks). Because the sample size is large (567 infrared galaxies and 2182 field galaxies), this sample is ideal for statistically examining the importance of interactions among infrared galaxies. In particular, the authors compare the nearest-neighbor distribution and the two-point correlation function of their sample with that of a control sample of field galaxies.

Gallimore, Jack F.↗

A comparison of Landsat point and rectangular field training sets for land-use classification

Rectangular training fields of homogeneous spectroreflectance are commonly used in supervised pattern recognition efforts. Trial image classification with manually selected training sets gives irregular and misleading results due to statistical bias. A self-verifying, grid-sampled training point approach is proposed as a more statistically valid feature extraction technique. A systematic pixel sampling network of every ninth row and ninth column efficiently replaced the full image scene with smaller statistical vectors which preserved the necessary characteristics for classification. The composite second- and third-order average classification accuracy of 50.1 percent for 331,776 pixels in the full image substantially agreed with the 51 percent value predicted by the grid-sampled, 4,100-point training set.

Tom, C. H.↗

Sampling errors in satellite estimates of tropical rain

The GATE rainfall data set is used in a statistical study to estimate the sampling errors that might be expected for the type of snapshot sampling that a low earth-orbiting satellite makes. For averages over the entire 400-km square and for the duration of several weeks, strong evidence is found that sampling errors less than 10 percent can be expected in contributions from each of four rain rate categories which individually account for about one quarter of the total rain.

Mcconnell, Alan↗

Image correlation and sampling study

The development of analytical approaches for solving image correlation and image sampling of multispectral data is discussed. Relevant multispectral image statistics which are applicable to image correlation and sampling are identified. The general image statistics include intensity mean, variance, amplitude histogram, power spectral density function, and autocorrelation function. The translation problem associated with digital image registration and the analytical means for comparing commonly used correlation techniques are considered. General expressions for determining the reconstruction error for specific image sampling strategies are developed.

Popp, D. J.↗

Uncertainty quantification tools for multiphase gas-solid flow simulations using MFIX

Computational fluid dynamics (CFD) has been widely studied and used in the scientific community and in the industry. Various models were proposed to solve problems in different areas. However, all models deviate from reality. Uncertainty quantification (UQ) process evaluates the overall uncertainties associated with the prediction of quantities of interest. In particular it studies the propagation of input uncertainties to the outputs of the models so that confidence intervals can be provided for the simulation results. In the present work, a non-intrusive quadrature-based uncertainty quantification (QBUQ) approach is proposed. The probability distribution function (PDF) of the system response can be then reconstructed using extended quadrature method of moments (EQMOM) and extended conditional quadrature method of moments (ECQMOM). The report first explains the theory of QBUQ approach, including methods to generate samples for problems with single or multiple uncertain input parameters, low order statistics, and required number of samples. Then methods for univariate PDF reconstruction (EQMOM) and multivariate PDF reconstruction (ECQMOM) are explained. The implementation of QBUQ approach into the open-source CFD code MFIX is discussed next. At last, QBUQ approach is demonstrated in several applications. The method is first applied to two examples: a developing flow in a channel with uncertain viscosity, and an oblique shock problem with uncertain upstream Mach number. The error in the prediction of the moment response is studied as a function of the number of samples, and the accuracy of the moments required to reconstruct the PDF of the system response is discussed. The QBUQ approach is then demonstrated by considering a bubbling fluidized bed as example application. The mean particle size is assumed to be the uncertain input parameter. The system is simulated with a standard two-fluid model with kinetic theory closures for the particulate phase implemented into MFIX. The effect of uncertainty on the disperse-phase volume fraction, on the phase velocities and on the pressure drop inside the fluidized bed are examined, and the reconstructed PDFs are provided for the three quantities studied. Then the approach is applied to a bubbling fluidized bed with two uncertain parameters, particle-particle and particle-wall restitution coefficients. Contour plots of the mean and standard deviation of solid volume fraction, solid phase velocities and gas pressure are provided. The PDFs of the response are reconstructed using EQMOM with appropriate kernel density functions. The simulation results are compared to experimental data provided by the 2013 NETL small-scale challenge problem. Lastly, the proposed procedure is demonstrated by considering a riser of a circulating fluidized bed as an example application. The mean particle size is considered to be the uncertain input parameter. Contour plots of the mean and standard deviation of solid volume fraction, solid phase velocities, and granular temperature are provided. Mean values and confidence intervals of the quantities of interest are compared to the experiment results. The univariate and bivariate PDF reconstructions of the system response are performed using EQMOM and ECQMOM.

97 MATHEMATICS AND COMPUTING↗

Testing of Microchannels and Lab-Grown Stress Corrosion Cracks for Quantification of Aerosol Transmission

The formation of a stress corrosion crack (SCC) in the canister wall of a dry cask storage system (DCSS) has been identified as a potential issue for the long-term storage of spent nuclear fuel. The presence of an SCC in a storage system could represent a through-wall flow path from the canister interior to the environment. Modern, vertical DCSSs are of particular interest due to the commercial practice of using higher backfill pressures in the canister, up to approximately 800 kPa, compared to their horizontal counterparts. This pressure differential offers a relatively high driving potential for blowdown of any particulates that might be present in the canister. In this study, the rates of gas flow and aerosol transmission of a spent fuel surrogate through an engineered microchannel with dimensions representative of an SCC were evaluated experimentally using coupled mass flow and aerosol analyzers. The microchannel was formed by mating two gage blocks with a linearly tapering slot orifice nominally 13 μm (0.005 in.) tall on the upstream side and 25 μm (0.0010 in.) tall on the downstream side. The orifice is 12.7 mm (0.500 in.) wide by 8.86 mm (0.349 in.) long (flow length). Surrogate aerosols of cerium oxide, CeO 2 , were seeded and mixed with either helium or air inside a pressurized tank. The aerosol characteristics were measured immediately upstream and downstream of the simulated SCC at elevated and ambient pressures, respectively. These data sets are intended to add to previous testing that characterized SCCs under well-controlled boundary conditions through the inclusion of testing improvements that establish initial conditions in a more consistent way. While the engineered microchannel has dimensions similar to actual SCCs, it does not reproduce the tortuous path the aerosol laden flow would have to traverse for eventual transmission. SCCs can be rapidly grown in a laboratory setting given the right conditions, and initial characterization and clean-flow testing has begun on lab grown crack samples provided to Sandia National Laboratories (SNL). Many such samples are required to produce statistically relevant transmission results, and SNL is developing a procedure to produce samples in welded steel plates. These ongoing testing efforts are focused on understanding the evolution in both size and quantity of a hypothetical release of aerosolized spent fuel particles from failed fuel to the canister interior and ultimately through an SCC.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Within-field variability of plant and soil parameters

The variability of ground truth data collected for vegetation experiments was investigated. Two fields of wheat and one field of corn were sampled on two different dates. The variability of crop and soil parameters within a field, between two fields of the same type, and within a field over time were compared statistically. The number of samples from each test site required in order to be able to determine with confidence the mean and standard deviations for a given variable was determined. Eight samples were found to be adequate for plant height determinations, while twenty samples were required for plant moisture and soil moisture characterization. Eighteen samples were necessary for detecting within field variability over time and for between field variability for the same crop. The necessary sample sites vary according to the physiological growth stage of the crop and recent weather events that affect the moisture and/or height characteristics of the field in question.

Ulaby, F. T.↗

Variability of cloudiness at airline cruise altitudes from GASP measurements

Additional statistics relating to the climatology of cloud cover at airline cruise altitudes are presented. The data were obtained between 1975 and 1979 from commercial airliners participating in the Global Atmospheric Sampling Program (GASP). The statistics describe the seasonal, latitudinal and altitudinal variation in cloudiness parameters as well as differences in the high-altitude cloud structure attributed to cyclone and convective-cloud generation processes. The latitudinal distribution of cloud cover derived form the GASP data was found to agree with high-altitude satellite observations. The relationships between three different measures of cloudiness and the relative vorticity at high altitudes is also discussed.

Jasperson, W. H.↗

H I 21 centimeter observations and I-band CCD surface photometry of spiral galaxies behind the Virgo Cluster and toward its antipode

Sample selection, radio and optical data acquisition and reduction, and observation results are presented for spiral galaxies behind the Virgo Cluster and toward its antipode. I-band CCD photometry was obtained for all the bright galaxies and part of the sample of faint galaxies in the two local volumes was studied. The statistical properties of the galaxy samples are discussed.

Lu, N. Y.↗

GIS-Based Modeling of Contaminated Soil Volumes at Multiple Sites in the Formerly Utilized Sites Remedial Action Program - 20149

The remediation of hazardous, toxic, and radioactive waste (HTRW) sites produces cost-related risks associated with the estimation of contaminated soil or debris volumes. Historical risk-management techniques include cost contingencies to cover volume uncertainties that affect project budgeting and decision-making. The Buffalo District teamed with project partners to lessen volume uncertainty and reduce project risks at multiple HTRW sites managed under the Formerly Utilized Sites Remedial Action Program (FUSRAP). Historical remedial investigations under FUSRAP commonly identified the presence of radiological material in site media, the associated human health risk, and then areas of remediation. To manage remedial execution and reduce risk, pre-design or remediation-phase sampling essentially 'chased' contamination, which was not conducive to efficient predictive budgeting derived from Feasibility Study (FS) cost analyses. The Buffalo District first optimized their approach to better understand volume uncertainty by utilizing the Argonne National Laboratory's Bayesian Approaches for Adaptive Spatial Sampling (BAASS) software [1]. BAASS processed soft data (e.g., gamma walk-over data) and spatial sampling data to estimate the lateral extent of contaminated soil irrespective of depth (i.e., gross contamination extent) and define areas of contaminant uncertainty. The software performed a binary transformation of contaminant concentrations at all sampling points based upon remedial action goals or a sum of ratios approach (i.e., clean, impacted, or range of impacts in soil). The model produced two-dimensional (horizontal) contaminant probability contours and statistical uncertainty in the sampling coverage and resulting contaminant extents. This method was translated vertically by partitioning the sampling data into depth brackets that produced a stacked representation of contaminant extents and uncertainty in the subsurface (i.e., similar to construction lifts). The results commonly led to a better understanding of project uncertainty and the need for sampling strategies that produce high-confidence soil volumes, which control costs. The BAASS-based delineations were eventually replaced by Empirical Bayesian Kriging (EBK) methods available in ArcGIS Spatial or 3D Analysts [2]. The EBK method calculates contaminant probability zones derived from user-controlled semivariograms of the spatial datasets. The resulting probability zones (e.g., 50% or 80% of contaminant probability) represent the two-dimensional surface delineation of the overall horizontal remedial area, similarly to BAASS. However, unlike BAASS, the vertical sampling data within these probability zones became vertical control points to contour a subterranean surface that connects subsurface points to the land-surface delineations of contamination. The resulting representation of horizontal and vertical impacts within an enclosed envelop (volume) of soil included uncertainty distributions that are used to plan uncertainty-reduction sampling. These data-driven and math-based models of three-dimensional sampling results produced well-bounded remedial volumes for project planning and better uncertainty predictions during project budgeting. The EBK method was applied to several FUSRAP sites managed by the Buffalo District and compared to less rigorously modeled sites previously remediated by the District. The comparison of modeled to actual remediated volumes provide a basis for validating the volume-estimation method. This comparison is important to ensure modeled volumes match physical boundaries of site remediation. FUSRAP sites with denser investigative sampling and lesser volume uncertainty proved useful in remedial planning and contracting. The Buffalo District noted that historical sites with sparser sampling arrays had greater disparity between estimated volumes and final remedial volumes. The benefit achieved over the cost of detailed soil sampling appears positive for FUSRAP projects, especially where impacts vary widely and appear unbounded by investigation-phase sampling. The subsequent Empirical Bayesian Kriging of contamination coupled with vertical contouring for soil estimations reduces uncertainty in soil volumes or indicates where sampling is required to reduce uncertainty, which together optimize remedial planning and budgeting. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Effects of finite sampling on fatigue damage estimation of wind turbine components: A statistical study

The variability of the wind turbine loads complicates fatigue assessment in the design phase, as performing simulations covering the entire lifetime is computationally expensive. The current work provides important information for assessing the uncertainty in fatigue damage estimation due to finite data. Here we study the sample size effect on mean, variance, and skewness of damage in each wind bin, identify the important wind bins, and study the uncertainty propagation from each wind bin to the lifetime damage using 3600 aeroelastic simulations and bootstrapping. To achieve less than 1% error in the damage estimation across all load channels in the current case study, at least 100 turbulence seeds are needed. Damage in different wind bins follows a lognormal distribution when using the conventional approach of six seeds. The provided insights and information allow the designer to achieve a specific level of accuracy for a given computational cost using strategic bin sampling.

17 WIND ENERGY↗