Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “random 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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Preferred Void Orientation in Uniaxially Pressed PBX 9502

Ultra- and small-angle neutron scattering techniques were used to study the morphology of voids in directions parallel and perpendicular to the compaction plane of a uniaxially pressed cylinder of PBX 9502. Two-dimensional scattering patterns from samples cut parallel to the compaction plane had circular equal-intensity contours, implying that along this view the voids are circular or randomly oriented. Identical measurements of samples cut perpendicular to the compaction plane produced elliptical equal-intensity contours, indicating preferred orientation of voids. Analysis of the two dimensional scattering patterns confirmed that the voids are non-spherical with an average aspect ratio of 1.20, flattened so that the larger dimensions of the voids lie within the compaction plane. Model analysis of the combined ultra- and small-angle neutron scattering curves revealed that the PBX 9502 porosity is interconnected, in the form of aggregates comprised of non-spherical, primary voids. The volume dimension was found to be the same in both directions, indicating that the aggregate is randomly oriented throughout the sample, and that it is the primary voids that are preferentially oriented. Furthermore, void interfaces with the surrounding matrix were found to be rough or polydisperse, with a common roughness dimension or polydispersity in all directions.

42 ENGINEERING↗
DOE OSTI2023

Unbalanced Nested Random Effects Estimation of Variance Components

To investigate the contributing factors of variance in the measurement of an iodine 127 (127I) sample, we implement a nested random effects analysis of variance (ANOVA). Historically, the reported uncertainty on a measurement of 127I has been obtained by methods of forward uncertainty propagation because there is typically only one replicate of a given sample for which to estimate the uncertainty. When samples are processed there are several types of quality control (QC) standards analyzed along-side the unknown samples with two to five replicates for each. Assuming the variance observed in these replicate QC standards is representative of that of the unknown samples, we use these data in a nested random effects ANOVA to estimate the total uncertainty of a measurement. We demonstrate this approach with two sets of measurements from Idaho National Laboratory and compare the results with the forward uncertainty propagation approach. Variance component estimates for the coarsest level of the nesting structure were most imprecise because replicates were most limited at these levels. We find that the results largely agree between forward propagation and ANOVA, and the greatest contributors of variance are due to instrument variation and chemical processing, with human processing being among the smallest contributors. This analysis provides reassurance that the reported uncertainties using forward propagation are reasonable and the process is well controlled. We propose a future designed experiment to increase replicates at the coarsest level of the hierarchy to improve estimates of these variance components.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

UQpy: A general purpose Python package and development environment for uncertainty quantification

In this paper, we present the UQpy software toolbox, an open-source Python package for general uncertainty quantification (UQ) in mathematical and physical systems. The software serves as both a user-ready toolbox that includes many of the latest methods for UQ in computational modeling and a convenient development environment for Python programmers advancing the field of UQ. The paper presents an introduction to the software's architecture and existing capabilities, divided in the code in a set of modules centered around different UQ tasks such as sampling methods, generation of random processes and random fields, probabilistic inverse modeling, reliability analysis, surrogate modeling, and active learning. The paper also highlights the importance of the RunModel module, which is used to drive simulations in the uncertainty analyses performed in UQpy. This module conveniently allows the user to define computational models directly in Python, or to run simulations from a third-party software in serial or in parallel. To illustrate the various capabilities, two examples are tracked throughout the paper and analyzed repeatedly for various UQ tasks. The first is a Python model solving a nonlinear structural dynamics problem, used to illustrate UQpy's capabilities in sampling and forward propagation of high dimensional random vectors (stochastic processes), and probabilistic inference. The second model is a third-party Abaqus finite element model solving the thermomechanical response of a beam structure. This example is used to illustrate UQpy's capabilities in variance reduction sampling techniques, reliability analysis, surrogate modeling and active learning techniques.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Fast estimation of the look-elsewhere effect using Gaussian random fields

Abstract We discuss the use of Gaussian random fields to estimate the look-elsewhere effect correction. We show that Gaussian random fields can be used to model the null-hypothesis significance maps from a large set of statistical problems commonly encountered in physics, such as template matching and likelihood ratio tests. Some specific examples are searches for dark matter using pixel arrays, searches for astronomical transients, and searches for fast-radio bursts. Gaussian random fields can be sampled efficiently in the frequency domain, and the excursion probability can be fitted with these samples to extend any estimation of the look-elsewhere effect to lower p values. In addition, in cases where the Gaussian random field is stationary and the parameter space is Euclidean, the look-elsewhere effect correction can be computed analytically. We demonstrate these methods using two example template matching problems. Finally, we apply these methods to estimate the trial factor of a $$4^3$$ 4 3 accelerometer array for the detection of dark matter tracks in the Windchime project. When a global significance of $$3\sigma $$ 3 σ is required, the estimated trial factor for such an accelerometer array is $$10^{14}$$ 10 14 for a one-second search, and $$10^{22}$$ 10 22 for a 1-year search.

Qin, Juehang (ORCID:0000000182288949)↗
DOE OSTI2025

Intelligent Sampling of Extreme-Scale Turbulence Datasets for Accurate and Efficient Spatiotemporal Model Training

With the end of Moore’s law and Dennard scaling, efficient training increasingly requires rethinking data volume. Can we train better models with significantly less data via intelligent subsampling? To explore this, we develop SICKLE, a sparse intelligent curation framework for efficient learning, featuring a novel maximum entropy (MaxEnt) sampling approach, scalable training, and energy benchmarking. We compare MaxEnt with random and phase-space sampling on large direct numerical simulation (DNS) datasets of turbulence. Evaluating SICKLE at scale on Frontier, we show that subsampling as a preprocessing step can, in many cases, improve model accuracy and substantially lower energy consumption, with observed reductions of up to 38×.

Brewer, Wes [ORNL] (ORCID:0000000236393956)↗
NASA NTRS2020

Designing an Optimal Ensemble Strategy for GMAO S2S Forecast System

The NASA Global Modeling and Assimilation Office (GMAO) Sub-seasonal to Seasonal (S2S) prediction system is being readied for a major upgrade. An important factor in successful extended range forecasting is the definition of the ensemble. Our overall strategy is to run a relatively large ensemble of about 40 members up to 3 months (focusing on the sub-seasonal forecast problem), after which we sub-sample the ensemble, and continue the forecast with about 10 members (up to 12 months). Here we present the results of our testing of various ways to generate the initial perturbations and the validation of a stratified sampling approach for choosing the members of the smaller ensemble. For the initialization of the ensemble we propose a combination of lagged and burst initial conditions. To generate perturbations for the burst ensemble members we used scaled differences of pairs of analysis states (chosen randomly from the corresponding season) separated by 1-10 days. We consider perturbing separately the atmosphere and the ocean, or both. By varying the separation times between the analysis states, we are able to produce perturbations that resemble well-known modes of variability. Focusing on the ENSO SST indices, we found that all types of perturbations are important for the ensemble spread with, however, considerable differences in the timing of the impacts on spread for the atmospheric and oceanic perturbations.Our initial (larger) ensemble size was determined so as to maximize the skill of predicting some of the leading modes of boreal winter atmospheric modes (namely the NAO, PNA and AO). Since it is not feasible for us to run with the larger ensemble beyond about 3 months, we employ a stratified sampling procedure that identifies the emerging directions of error growth to subset the ensemble. By comparing the results from the stratified ensemble with that of the randomly sampled ensemble of the same size, we find that the former provides substantially better estimates the mean of the original large ensemble.

Borovikov, Anna↗
DOE OSTI2025

Topography restoration in white-light interferometry using an instrument transfer function evaluated with binary pseudo-random arrays

Accurate knowledge of the instrument transfer function (ITF) is vital for topography measurements using white-light interferometry (WLI). To this end, we derive a complete set of analytical expressions for the power spectral density (PSD) of a discretely sampled binary pseudo-random array (BPRA) as a theoretical benchmark. We then determine the ITF by comparing this theoretical PSD with the measured PSD of the BPRA. For the Zygo ZeGage Pro HR with a objective, the determined ITF closely matches the nominal modulation transfer function (MTF). Accordingly, we integrate the nominal MTF into Fourier-domain restoration filters and apply them to practical WLI topography data, yielding restored topographies that show improved agreement with atomic force microscopy measurements. Overall, our integrated BPRA-based ITF measurement and filtering methodology offers a robust yet practical means of characterizing WLI performance and enhancing the accuracy of surface measurements.

Hirose, Shigenobu [Japan Agency for Marine Earth S↗
DOE OSTI2026

Reducing Randomized Quantum Algorithm Cost [Slides]

We derive the optimal sampling strategy for minimizing total resource cost in randomized quantum algorithms. Our framework is completely general, allowing for resources as diverse as gate counts circuit depth, runtime, or even dissipated energy.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1988

Sampling inhomogeneous turbulent fields

The reconstruction of an inhomogeneous random process from a finite number of discrete samples can be performed in terms of the Karhunen-Loeve (KL) expansion for that process. The n(th) eigenfunction has n - 1 zero crossings which are the sampling points for the inhomogeneous process. The rapid variation of the KL eigenfunctions makes it unnecessary to have a high density of sampling (or grid points) near the wall. However, this result should not be construed to indicate that with spectral simulations significantly fewer grid points are required with the KL expansion as compared to other orthogonal expansions. Moin and Moser (1989) have shown that the advantage of the KL expansion over Chebychev expansion rapidly diminishes when high percentage (say 90 percent) energy recovery is demanded.

Adrian, R. J.↗
NASA NTRS1983

Spatial inventory integrating raster databases and point sample data

A timber inventory of the Eldorado National Forest, located in east-central California, provides an example of the use of a Geographic Information System (GIS) to stratify large areas of land for sampling and the collection of statistical data. The raster-based GIS format of the VICAR/IBIS software system allows simple and rapid tabulation of areas, and facilitates the selection of random locations for ground sampling. Algorithms that simplify the complex spatial pattern of raster-based information, and convert raster format data to strings of coordinate vectors, provide a link to conventional vector-based geographic information systems.

Strahler, A. H.↗
NASA NTRS1997

Two Extremely Red Galaxies

This report concerns one of the major observational studies in the ISO Central Programme, the ISO Normal Galaxy Survey. This is a survey of an unbiased sample of spiral and lenticular galaxies selected from the Revised Shapley-Ames Catalog. It is therefore optically-selected, with a brightness limit of blue magnitude = 12, and otherwise randomly chosen. The original sample included 150 galaxies, but this was reduced to 74 when the allocated observing time was expended because the ISO overheads encountered in flight were much larger than predicted.

Joseph, Robert D.↗
DOE OSTI2025

National serosurvey and risk mapping reveal widespread distribution of Coxiella burnetii in Kenya

Coxiella burnetii, the causative agent of Q fever, is an emerging pathogen that has the potential to cause severe chronic infections in animals and humans worldwide. The detrimental impact on public health is projected to be higher in the low- and middle-income countries given their lower capacity to sustain effective surveillance and response measures. We implemented a national serosurvey of cattle in Kenya to map the spatial distribution of the pathogen. The study used serum samples that were collected from randomly selected cattle in different ago-ecological zones across the country. These samples were screened for the pathogen using PrioCHECK Ruminant Q Fever AB Plate ELISA kit. The laboratory findings were analyzed using INLA package to identify risk factors for C. burnetii exposure from herd- and animal-level factors, area, and bioclimatic datasets accessed from online databases. A total of 6,593 cattle were recruited for the study; of these, 7.9% (95% CI; 7.2–8.5) were seropositive. Outputs from the multivariable analysis revealed that the animal age and some of the geographical variables including wind speed, area under shrubs and “petric calcisols” type of soil were significantly associated with C. burnetii seropositivity. Being a calf, weaner or subadult was associated with lower odds of exposure compared to being an adult by 0.24 (credibility interval: 2.5% and 97.5%), 0.41 (0.30–0.55) and 0.51 (0.38–0.69), respectively. In addition, a unit increase in the wind speed increased the odds of C. burnetii seropositivity by 1.27 (1.05–1.52) while an increase on the land area under shrubs was associated with lower odds of exposure (0.67 [0.47–0.69]). The effect of petric calcisols was non-linear; an increase of the land area with this soil type was associated with an exponential increase in C. burnetii seropositivity. This study provides new data on C. burnetii seroprevalence, information of its risk factors and a prevalence map that can be used for C. burnetii risk surveillance and control. The identification of environmental risk factors for C. burnetii exposure, and the increasing awareness of the zoonotic potential of the pathogen, calls for the need to enhance the existing collaborations for the surveillance and control of C. burnetii in line with the One Health framework. The evidence generated on the potential role of environmental factors can also be used to design nature-based interventions, such as replacement of vegetation in denuded areas, to reduce potential for the aerosolization of the pathogen. Livestock vaccination in the hotspots would also reduce animal infections and hence the contamination of the environment.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Circulating levels of micronutrients and risk of osteomyelitis: a Mendelian randomization study

Background Few observational studies have investigated the effect of micronutrients on osteomyelitis, and these findings are limited by confounding and conflicting results. Therefore, we conducted Mendelian randomization (MR) analyses to evaluate the association between blood levels of eight micronutrients (copper, selenium, zinc, vitamin B12, vitamin C, and vitamin D, vitamin B6, vitamin E) and the risk of osteomyelitis. Methods We performed the two-sample and multivariable Mendelian randomization (MVMR) to investigate causation, where instrument variables for the predictor (micronutrients) were derived from the summary data of micronutrients from independent cohorts of European ancestry. The outcome instrumental variables were used from the summary data of European-ancestry individuals ( n = 486,484). The threshold of statistical significance was set at p < 0.00625. Results We found a significant causal association that elevated zinc heightens the risk of developing osteomyelitis in European ancestry individuals OR = 1.23 [95% confidence interval (CI) [1.07, 1.43]; p = 4.26E-03]. Similarly, vitamin B6 showed a similar significant causal effect on osteomyelitis as a risk factor OR = 2.78 (95% CI [1.34, 5.76]; p = 6.04E-03; in the secondary analysis). Post-hoc analysis suggested this result (vitamin B6). However, the multivariable Mendelian randomization (MVMR) provides evidence against the causal association between zinc and osteomyelitis OR = 0.98(95% CI [−0.11, 0.07]; p = 7.20E-1). After searching in PhenoScanner, no SNP with confounding factors was found in the analysis of vitamin B6. There was no evidence of a reverse causal impact of osteomyelitis on zinc and vitamin B6. Conclusion This study supported a strong causal association between vitamin B6 and osteomyelitis while reporting a dubious causal association between zinc and osteomyelitis.

Zhang, Xu↗
DOE OSTI2020

UQ Toolkit v 2.0

The Uncertainty Quantification (UQ) Toolkit is a software library for the characterizaton and propagation of uncertainties in computational models. For the characterization of uncertainties, Bayesian inference tools are provided to infer uncertain model parameters, as well as Bayesian compressive sensing methods for discovering sparse representations of high-dimensional input-output response surfaces, and also Karhunen-Loève expansions for representing stochastic processes. Uncertain parameters are treated as random variables and represented with Polynomial Chaos expansions (PCEs). The library implements several spectral basis function types (e.g. Hermite basis functions in terms of Gaussian random variables or Legendre basis functions in terms of uniform random variables) that can be used to represent random variables with PCEs. For propagation of uncertainty, tools are provided to propagate PCEs that describe the input uncertainty through the computational model using either intrusive methods (Galerkin projection of equations onto basis functions) or non-intrusive methods (perform deterministic operation at sampled values of the random values and project the obtained results onto basis functions).

Safta, Cosmin↗
NASA NTRS1994

Two-sample discrimination of Poisson means

This paper presents a statistical test for detecting significant differences between two random count accumulations. The null hypothesis is that the two samples share a common random arrival process with a mean count proportional to each sample's exposure. The model represents the partition of N total events into two counts, A and B, as a sequence of N independent Bernoulli trials whose partition fraction, f, is determined by the ratio of the exposures of A and B. The detection of a significant difference is claimed when the background (null) hypothesis is rejected, which occurs when the observed sample falls in a critical region of (A, B) space. The critical region depends on f and the desired significance level, alpha. The model correctly takes into account the fluctuations in both the signals and the background data, including the important case of small numbers of counts in the signal, the background, or both. The significance can be exactly determined from the cumulative binomial distribution, which in turn can be inverted to determine the critical A(B) or B(A) contour. This paper gives efficient implementations of these tests, based on lookup tables. Applications include the detection of clustering of astronomical objects, the detection of faint emission or absorption lines in photon-limited spectroscopy, the detection of faint emitters or absorbers in photon-limited imaging, and dosimetry.

Lampton, M.↗
DOE OSTI2024

Computationally efficient and error aware surrogate construction for numerical solutions of subsurface flow through porous media

Limiting the injection rate to restrict the pressure below a threshold at a critical location can be an important goal of simulations that model the subsurface pressure between injection and extraction wells. The pressure is approximated by the solution of Darcy’s partial differential equation for a given permeability field. The subsurface permeability is modeled as a random field since it is known only up to statistical properties. This induces uncertainty in the computed pressure. Solving the partial differential equation for an ensemble of random permeability simulations enables estimating a probability distribution for the pressure at the critical location. These simulations are computationally expensive, and practitioners often need rapid online guidance for real-time pressure management. An ensemble of numerical partial differential equation solutions is used to construct a Gaussian process regression model that can quickly predict the pressure at the critical location as a function of the extraction rate and permeability realization. The Gaussian process surrogate analyzes the ensemble of numerical pressure solutions at the critical location as noisy observations of the true pressure solution, enabling robust inference using the conditional Gaussian process distribution. Our first novel contribution is to identify a sampling methodology for the random environment and matching kernel technology for which fitting the Gaussian process regression model scales as O ( n log n ) instead of the typical O ( n 3 ) rate in the number of samples n used to fit the surrogate. The surrogate model allows almost instantaneous predictions for the pressure at the critical location as a function of the extraction rate and permeability realization. Our second contribution is a novel algorithm to calibrate the uncertainty in the surrogate model to the discrepancy between the true pressure solution of Darcy’s equation and the numerical solution. Finally, although our method is derived for building a surrogate for the solution of Darcy’s equation with a random permeability field, the framework broadly applies to solutions of other partial differential equations with random coefficients.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Elemental profiling and genomewide association studies reveal genomic variants modulating ionomic composition in Populus trichocarpa leaves

Samples were collected from a population of 1,089 black cottonwood genotypes (P. trichocarpa) assembled from native stands to encompass the central portion of the natural range of the species, stretching from 38.8° to 54.3° N Q13 latitude from California, USA, to British Columbia, Canada. Establishment of the common garden, growth conditions, and site maintenance have been described by Muchero et al (2015). In this study, leaf samples for ionomic profiling were collected from 4-year-old trees, during the growing season, in July 2012, from a field located in Clatskanie, Oregon, USA (46°6′11″N 123°12′13″W). The field site was located in a protected alluvial floodplain containing a uniform Wauna-Locoda silt loam soil area characterized by an acidic pH, in Columbia County, Oregon. A subset of 584 out of the 1,089 P. trichocarpa genotypes were represented in this sampling. These genotypes were randomly selected to represent the geographical distribution of the population. A single fully mature (LPI 7-9) leaf on the south side of the tree exposed to full sunlight conditions was removed from the tree within a 6-hour window centering on solar noon and immediately frozen under dry ice before processing. Leaf samples of 584 P. trichocarpa genotypes were finely ground to 40 mm particle size using a mortar and pestle, and ionomic composition was analyzed using ICP-MS. In total, 20 elements were profiled, including aluminum (Al27), arsenic (As75), boron (B11), cadmium (Cd111), calcium (Ca43), cobalt (Co), copper (Cu), iron (Fe57), magnesium (Mg25), manganese (Mn55), molybdenum (Mo), nickel (Ni60), phosphorus (P31), potassium (K39), rubidium (Rb85), selenium (Se82), sodium (Na23), strontium (Sr88), sulfur (S34), and zinc (Zn66), following a protocol established by Ziegler et al. (2013). For each sample, 75mg of powder was digested overnight in 2.5 mL HNO3 containing 20 parts per billion (ppb) indium as an internal standard, following the protocol described in Ziegler et al. (2013). Following a dilution, concentration of the 20 elements was measured using an Elan 6000 DRC-e mass spectrometer (Perkin-Elmer SCIEX) connected to a PFA microflow nebulizer (Elemental Scientific) and Apex HF desolvator (Elemental Scientific). One measurement per sample per genotype was done. For subsequent analyses, the quantifications were converted to total element concentration.

CBI ionomics, GWAS, plasma-mass spectrometry, neut↗
NASA NTRS1985

Long term monitoring of methane in the atmosphere by multiplex gas chromatography

Methane is of interest in the study of the Earth's atmosphere because of its implication in the future global warming of the surface. This warming is produced by the absorption of infrared energy by trace gases. It has been estimated that in the next 40 to 50 years, methane could contribute 20 to 25 pct. as much atmospheric warming as that expected from carbon dioxide increases. Studies to examine sources, sinks, and cycles of methane will require analytical methods capable of continuous unattended measurement with temporal resolution of an hour or less for weeks at a time. Gas chromatography (GC) is one of the most practical methods available to conduct the analysis of air, but limitations in this technique still exist which can be alleviated with multiplex GC (MGC). MGC is a technique where many samples are pseudo-randomly introduced to the chromatograph without regard to the length of time required for an analysis. The resulting data must then be reduced using computational methods such as cross correlation. In the technique reported, a tube packed with silver oxide was used at the inlet of the GC column to create concentration pulses of methane in a sample stream of air. By using only one carrier, i.e., ambient air, an effective and accurate method to monitor the variations in concentration of methane in the atmosphere over long periods of time was developed. Methane in ambient air was monitored for an eight day period and an interesting temporal variability was found. This work has shown the utility of a relatively simple MGC for the analysis of a real environmental sample.

Valentin, Jose R.↗