Search NASA⌕ Search

SEARCH · Search NASA

Results for “randomized sampling”

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 811 records · Page 45

Deciphering the source of primary biological aerosol particles: a pollen case study

Primary biological aerosol particles (PBAPs) are microscopic solids suspended in the atmosphere emitted by biological systems and play critical roles in the atmosphere and at the atmosphere-biosphere interface, impacting human health, climate, and the ecosystem function. Understanding the sources of PBAPs is necessary to decipher the mechanistic interactions between aerosols, climate, and distinct ecosystem components. However, the detection of specific PBAPs in complex ambient aerosol samples is challenging. We performed metabolomics analyses of pollen from three pollinating tree species and ambient samples collected during the peak pollination period of each species. Random Forest and sPLS-DA machine learning methods were employed to evaluate whether metabolic signatures of ambient samples can reveal the source of the main pollen particles present in the atmosphere. Our results suggest that atmospheric eco-metabolomics techniques combined with sophisticated statistical methods can decipher the origin of abundant PBAPs from complex ambient samples. Developing complete libraries containing high-resolution metabolomic fingerprints of the major PBAPs present in the atmosphere would significantly advance future research to accurately understand the role of PBAPs in the atmosphere, ecosystems and human health.

Rivas-Ubach, Albert↗

Assessing overall building energy performance of a large population of residential single-family homes using limited field data

The U.S. Department of Energy (DOE), through the Building Energy Codes Program (BECP), supports building energy code development and implementation processes to achieve the maximum practicable, cost-effective improvements in energy efficiency while providing safe, healthy buildings for occupants. Currently, there is a lack of statistically representative data that shows how energy codes are implemented in the field and the energy savings opportunities that exist. To address this gap, DOE initiated a research study targeted at new residential single-family homes in 8 states starting in the fall of 2014. The study comprises of a baseline phase to determine current levels of compliance with the key energy provisions of the residential building energy code, identify issues, and calculate potential energy savings from improved compliance; a second phase implementing education, training and outreach activities designed to increase targeted code compliance; and a post-study phase to measure the impact of the training and outreach activities compared to the baseline phase. Initial analyses determined the code requirements with the largest energy impact and the minimum number of observations required to detect statistically significant differences in energy use. Data was collected in the field during a single site visit to avoid bias from repeat visits. However, the single site visit principle prevented energy simulation of individual homes due to the lack of all required model inputs. To address this, energy simulation was conducted using Pacific Northwest National Laboratory’s prototype building models and employing bootstrap random sampling on the data collected from all homes visited. The modeled energy consumption using the observed field data was evaluated relative to that associated with the minimum code requirements and potential savings associated with improved compliance of select items were assessed. This paper discusses the methodology developed for this field study, and presents the results and findings from the baseline phase.

Xie, YuLong↗

A Markov Random Field Model for Texture-based Segmentation of Small Cracks in Thin InconelTubes

Markov random fields have been used for image segmentation since their introduction in the1980s. This work applies a method from Principal Component Thermography to enhance thecontrast in damage regions in 3D images derived from X-ray computed tomography (CT)inspections. The developed method is applied to sizing of small cracks in thin Inconel tubes designed as probability of detection (POD) samples for radiographic inspection.Misclassification errors arising from artifacts due to beam-hardening are reduced by fitting the boundary of the segmented damage region with the arc of an ellipse. Results are comparedagainst those obtained through manual inspection.

Image segmentation↗

Density estimation via measure transport: Outlook for applications in the biological sciences

Abstract One among several advantages of measure transport methods is that they allow or a unified framework for processing and analysis of data distributed according to a wide class of probability measures. Within this context, we present results from computational studies aimed at assessing the potential of measure transport techniques, specifically, the use of triangular transport maps, as part of a workflow intended to support research in the biological sciences. Scenarios characterized by the availability of limited amount of sample data, which are common in domains such as radiation biology, are of particular interest. We find that when estimating a distribution density function given limited amount of sample data, adaptive transport maps are advantageous. In particular, statistics gathered from computing series of adaptive transport maps, trained on a series of randomly chosen subsets of the set of available data samples, leads to uncovering information hidden in the data. As a result, in the radiation biology application considered here, this approach provides a tool for generating hypotheses about gene relationships and their dynamics under radiation exposure.

gene expression data↗

Galaxy–galaxy lensing estimators and their covariance properties

Here, we study the covariance properties of real space correlation function estimators – primarily galaxy–shear correlations, or galaxy–galaxy lensing – using SDSS data for both shear catalogues and lenses (specifically the BOSS LOWZ sample). Using mock catalogues of lenses and sources, we disentangle the various contributions to the covariance matrix and compare them with a simple analytical model. We show that not subtracting the lensing measurement around random points from the measurement around the lens sample is equivalent to performing the measurement using the lens density field instead of the lens overdensity field. While the measurement using the lens density field is unbiased (in the absence of systematics), its error is significantly larger due to an additional term in the covariance. Therefore, this subtraction should be performed regardless of its beneficial effects on systematics. Comparing the error estimates from data and mocks for estimators that involve the overdensity, we find that the errors are dominated by the shape noise and lens clustering, which empirically estimated covariances (jackknife and standard deviation across mocks) that are consistent with theoretical estimates, and that both the connected parts of the four-point function and the supersample covariance can be neglected for the current levels of noise. While the trade-off between different terms in the covariance depends on the survey configuration (area, source number density), the diagnostics that we use in this work should be useful for future works to test their empirically determined covariances.

79 ASTRONOMY AND ASTROPHYSICS↗

Probabilistic effects of porosity and chemical kinetics on the shock initiation of an octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (HMX) based explosive

In this work, we investigate the effects of porosity and chemical reactivity on the shock initiation of an HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine)-based polymer bonded explosive due to thin flyer impact using two-dimensional mesoscale simulations. Microstructural samples are computationally generated, and the calculated shock response is compared to experiments. The sample dimensions are 0.1×1.05 mm 2 , and each sample accounts for grain scale heterogeneities, including randomized grain morphology and pore spacing. Simulations are performed using the Lawrence Livermore National Laboratory (LLNL) arbitrary Lagrangian–Eulerian multi-physics code ALE3D. Chemistry and equation of states are handled with LLNL's thermochemical code Cheetah. We employ a modified Arrhenius-based chemical kinetics model that scales the temperature exponent fitting parameter with flyer shock conditions. This new model is calibrated to the mean response of experimental shock initiation thresholds for three separate Kapton flyer thicknesses. A simplified shock initiation criterion is proposed, which classifies the macroscale initiation response based on cross-sectional temperature profiles of the shocked sample. This criterion is used to further investigate the probabilistic response of the sample to variations in the mean porosity and chemical kinetics rate, and we find an inverse correlation between the shock pressure and the confidence interval of the experimental threshold velocity. Furthermore, the methodology presented in this work for quantifying the probabilistic nature of shock initiation is material agnostic. These mesoscale studies may be used to inform continuum scale models and expedite the development of future, tailored energetic formulations.

36 MATERIALS SCIENCE↗

A Practical Methodology for Quantifying Random and Systematic Components of Unexplained Variance in a Wind Tunnel

This paper documents a check standard wind tunnel test conducted in the Langley 0.3-Meter Transonic Cryogenic Tunnel (0.3M TCT) that was designed and analyzed using the Modern Design of Experiments (MDOE). The test designed to partition the unexplained variance of typical wind tunnel data samples into two constituent components, one attributable to ordinary random error, and one attributable to systematic error induced by covariate effects. Covariate effects in wind tunnel testing are discussed, with examples. The impact of systematic (non-random) unexplained variance on the statistical independence of sequential measurements is reviewed. The corresponding correlation among experimental errors is discussed, as is the impact of such correlation on experimental results generally. The specific experiment documented herein was organized as a formal test for the presence of unexplained variance in representative samples of wind tunnel data, in order to quantify the frequency with which such systematic error was detected, and its magnitude relative to ordinary random error. Levels of systematic and random error reported here are representative of those quantified in other facilities, as cited in the references.

Deloach, Richard↗

Investigation of properties of CdSeTe/Sb and Sb/CdSeTe stacks

The effect of Sb inclusion in CdSeTe thin films was examined. CdSeTe alloys were grown by co-evaporation, either on 5 nm thick electron-beam evaporated Sb layers (Sb/CdSeTe) or coated with Sb post-deposition (CdSeTe/Sb). Sb inclusion improved CdSeTe grain structure, a critical factor for device performance, with average grain size increasing from 0.6 to 1.5 μm. All samples crystallized on annealing to the cubic structure, with a change in film orientation from (111) to random or (220). The presence of Sb in CdSeTe/Sb or Sb/CdSeTe samples, however, had a limited effect on the carrier concentration compared to Sb-free samples. Time resolved photoluminescence showed Sb inclusion increased decay lifetimes, with τ1 = 0.8 ns for and τ2 = 6.9 ns for the Sb/CdSeTe sample.

14 SOLAR ENERGY↗

Combined Uncertainty and A-Posteriori Error Bound Estimates for General CFD Calculations: Theory and Software Implementation

This workshop presentation discusses the design and implementation of numerical methods for the quantification of statistical uncertainty, including a-posteriori error bounds, for output quantities computed using CFD methods. Hydrodynamic realizations often contain numerical error arising from finite-dimensional approximation (e.g. numerical methods using grids, basis functions, particles) and statistical uncertainty arising from incomplete information and/or statistical characterization of model parameters and random fields. The first task at hand is to derive formal error bounds for statistics given realizations containing finite-dimensional numerical error [1]. The error in computed output statistics contains contributions from both realization error and the error resulting from the calculation of statistics integrals using a numerical method. A second task is to devise computable a-posteriori error bounds by numerically approximating all terms arising in the error bound estimates. For the same reason that CFD calculations including error bounds but omitting uncertainty modeling are only of limited value, CFD calculations including uncertainty modeling but omitting error bounds are only of limited value. To gain maximum value from CFD calculations, a general software package for uncertainty quantification with quantified error bounds has been developed at NASA. The package provides implementations for a suite of numerical methods used in uncertainty quantification: Dense tensorization basis methods [3] and a subscale recovery variant [1] for non-smooth data, Sparse tensorization methods[2] utilizing node-nested hierarchies, Sampling methods[4] for high-dimensional random variable spaces.

CFD↗

Derivative-free stochastic optimization via adaptive sampling strategies

In this paper, we present a novel derivative-free framework for solving unconstrained stochastic optimization problems. Many problems in fields ranging from simulation optimization to reinforcement learning to quantum computing involve settings where only stochastic function values are obtained via a zeroth-order oracle, which has no available gradient information and necessitates the usage of derivative-free optimization methodologies. Our approach includes estimating gradients using stochastic function evaluations and integrating adaptive sampling techniques to control the accuracy in these stochastic approximations. Our framework encapsulates several gradient estimation techniques, including standard finite-difference, Gaussian smoothing, sphere smoothing, randomized coordinate finite-difference, and randomized subspace finite-difference methods. We provide theoretical convergence guarantees for our framework and analyze the worst-case iteration and sample complexities associated with each gradient estimation method. Finally, we demonstrate the empirical performance of the methods on logistic regression and nonlinear least squares problems.

Adaptive sampling↗

Information density and efficiency of two-dimensional /2-D/ sampled imagery

Information density and efficiency (i.e., the ratio of information density to data density) are used as criteria for assessing the quality of 2-D sampled and quantized imagery as a function of the statistical properties of random radiance fields, the spatial response (PSF or MTF) and sensitivity of imaging systems, and the sampling and quantization intervals. Computational results are intuitively satisfying: they are consistent with experimental and theoretical results obtained by earlier investigators concerned with the performance of TV cameras, and they provide useful guidelines for optimizing the design of line-scan and sensor-array imaging systems, especially if these systems use a digital communication link for transmitting data.

Huck, F. O.↗

Stochastic Stability of Sampled Data Systems with a Jump Linear Controller

In this paper an equivalence between the stochastic stability of a sampled-data system and its associated discrete-time representation is established. The sampled-data system consists of a deterministic, linear, time-invariant, continuous-time plant and a stochastic, linear, time-invariant, discrete-time, jump linear controller. The jump linear controller models computer systems and communication networks that are subject to stochastic upsets or disruptions. This sampled-data model has been used in the analysis and design of fault-tolerant systems and computer-control systems with random communication delays without taking into account the inter-sample response. This paper shows that the known equivalence between the stability of a deterministic sampled-data system and the associated discrete-time representation holds even in a stochastic framework.

Gonzalez, Oscar R.↗

Stochastic Stability of Nonlinear Sampled Data Systems with a Jump Linear Controller

This paper analyzes the stability of a sampled- data system consisting of a deterministic, nonlinear, time- invariant, continuous-time plant and a stochastic, discrete- time, jump linear controller. The jump linear controller mod- els, for example, computer systems and communication net- works that are subject to stochastic upsets or disruptions. This sampled-data model has been used in the analysis and design of fault-tolerant systems and computer-control systems with random communication delays without taking into account the inter-sample response. To analyze stability, appropriate topologies are introduced for the signal spaces of the sampled- data system. With these topologies, the ideal sampling and zero-order-hold operators are shown to be measurable maps. This paper shows that the known equivalence between the stability of a deterministic, linear sampled-data system and its associated discrete-time representation as well as between a nonlinear sampled-data system and a linearized representation holds even in a stochastic framework.

Gonzalez, Oscar R.↗

Information transfer as a framework for optimized phase imaging

In order to efficiently image a non-absorbing sample (a phase object), dedicated phase contrast optics are required. Typically, these optics are designed with the assumption that the sample is weakly scattering, implying a linear relation between a sample’s phase and its transmission function. In the strongly scattering, nonlinear case, the standard optics are ineffective, and the transfer functions used to characterize them are uninformative. We use the Fisher information (FI) to assess the efficiency of various phase imaging schemes and to calculate an information transfer function (ITF). We show that a generalized version of Zernike phase contrast is efficient given sufficient prior knowledge of the sample. We show that with no prior knowledge, a random sensing measurement yields a significant fraction of the available information. Finally, we introduce a generalized approach to common path interferometry that can be optimized to prioritize sensitivity to particular sample features. Each of these measurements can be performed using Fourier lenses and phase masks.

Koppell, Stewart↗

X-ray thermal diffuse scattering as a texture-robust temperature diagnostic for dynamically compressed solids

We present a model of x-ray thermal diffuse scattering (TDS) from a cubic polycrystal with an arbitrary crystallographic texture, based on the classic approach of Warren [B. E. Warren, Acta Crystallogr. 6, 803 (1953)]. We compare the predictions of our model with femtosecond x-ray diffraction patterns gathered from ambient and dynamically compressed rolled copper foils obtained at the High Energy Density instrument of the European X-Ray Free-Electron Laser facility and find that the texture-aware TDS model yields more accurate results than does the conventional powder model owed to Warren. Nevertheless, we further show: with sufficient angular detector coverage, the TDS signal is largely unchanged by sample orientation and in all cases strongly resembles the signal from a perfectly random powder; shot-to-shot fluctuations in the TDS signal resulting from grain-sampling statistics are at the percent level, in stark contrast to the fluctuations in the Bragg-peak intensities (which are over an order of magnitude greater); and TDS is largely unchanged even following texture evolution caused by compression-induced plastic deformation. We conclude that TDS is robust against texture variation, making it a flexible temperature diagnostic applicable just as well to off-the-shelf commercial foils as to ideal powders.

Crystal lattices↗

Bayesian estimation of Karhunen–Loève expansions; A random subspace approach

One of the most widely-used statistical procedures for dimensionality reduction of high dimensional random fields is Principal Component Analysis (PCA), which is based on the Karhunen-Lo eve expansion (KLE) of a stochastic process with finite variance. The KLE is analogous to a Fourier series expansion for a random process, where the goal is to find an orthogonal transformation for the data such that the projection of the data onto this orthogonal subspace is optimal in the L 2 sense, i.e, which minimizes the mean square error. In practice, this orthogonal transformation is determined by performing an SVD (Singular Value Decomposition) on the sample covariance matrix or on the data matrix itself. Sampling error is typically ignored when quantifying the principal components, or, equivalently, basis functions of the KLE. Furthermore, it is exacerbated when the sample size is much smaller than the dimension of the random field. In this paper, we introduce a Bayesian KLE procedure, allowing one to obtain a probabilistic model on the principal components, which can account for inaccuracies due to limited sample size. The probabilistic model is built via Bayesian inference, from which the posterior becomes the matrix Bingham density over the space of orthonormal matrices. We use a modified Gibbs sampling procedure to sample on this space and then build a probabilistic Karhunen-Lo eve expansions over random subspaces to obtain a set of low-dimensional surrogates of the stochastic process. We illustrate this probabilistic procedure with a finite dimensional stochastic process inspired by Brownian motion.

97 MATHEMATICS AND COMPUTING↗

Hybrid Ultrasonic TOFD Imaging of Weld Flaws Using Wavelet Transforms and Image Registration

Ultrasonic time of flight diffraction (TOFD) is an effective weld crack inspection technique. Due to the intensity of diffraction wave is rather weak compared with the lateral wave and the bottom echo wave, thus the signal-to-ratio (SNR) of TOFD image is low. A new dichotomous method is comprised of two steps that contains wavelet shrinkage and image registration is proposed in this paper to reduce the noise and improve the resolution of TOFD images as well. In order to evaluate the reliability of our proposed method in this paper, we have established the experiment system, and sampled a number of TOFD data with random distribution of noise characteristics. We adopted one-dimension wavelet transform and two-dimension wavelet transform in the very beginning of the first step of the proposed algorithm respectively. The SNR of the result obtained in this step is improved significantly compared with the classic algorithms. Next, the image registration is applied. After the registered images have been added to form a new one, then it comes to the final result that shows not only the SNR but also the definition of the image is enhanced effectively.

42 ENGINEERING↗

A statistical approach to screening isotopic signatures in monitoring for underground nuclear explosions

The ability to differentiate between atmospheric radionuclide signatures from underground nuclear explosions (UNEs) and signals from other sources, such as medical isotope-production facilities and nuclear reactors, can be critical to the detection and monitoring of unannounced, low-yield nuclear events. Signatures having anomalously high amplitudes, compared to background levels, remain the best indicator in screening for a UNE. However, isotopic composition can further validate a suspected UNE signature, but separation from any atmospheric background composition is first necessary. To date, evaluating the challenges of performing this separation has typically involved comparing an observed background with a highly idealized deterministic model of radioxenon signature production by a UNE that does not consider the influence of post-detonation chemical/physical processes in the detonation cavity or the subsequent gas transport mechanisms that can also affect the isotopic composition of the detected gas signature. In addition, purely deterministic models, as previously employed, overlook the uncertainty inherent in estimating critical parameters characterizing the UNE and its detonation environment. In this paper, we create detailed, multi-parameter models of radionuclide evolution using the widely accepted England and Rider post-detonation radionuclide decay-chain network coupled to detailed models simulating physical production and transport processes affecting the gas signature. Because these models are governed by uncertain parameters including barometric fluctuations, realistic ranges of variation for each of the parameters influencing isotopic composition are then defined. A Latin-Hypercube sampling approach is used to obtain a random distribution of isotopic production and gas transport results associated with a given value of each parameter. We apply these results to background histories of two stations, one providing 4-isotope background measurements and the other providing two-isotope measurements associated with the 2013 DPRK announced UNE.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗