Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “statistical sampling techniques”

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 127 records · Page 7

Quantitative Particle Analysis of Neptunium-237 Oxides: Optimization of MAMA Analysis for Modified Direct Denitration Products

The production of plutonium-238 through irradiation of neptunium-237 ( 237 Np) target materials for the use in radioisotope thermoelectric generators is paramount for continued deep space exploration. This work employs scanning electron microscopy to analyze 237 Np materials coupled with a well-developed image analysis framework (Morphological Analysis for Material Attribution, or MAMA) to determine the degree of micron-scale homogeneity in the materials. This work demonstrated how the quantification of particle characteristics can validate production materials and affirm the qualitative similarities observed in micrographs. The 237 Np oxide particle analysis determined that the materials from five production runs were quantitatively homogenous (significant at α = 0.05) in particle area, circularity, equivalent circular diameter, and ellipse aspect ratio, with two of the sampling dates having statistically significant different means for one of the four characteristics. Furthermore, these metrics not only confirm general homogeneity of the material but also expand the application of MAMA workflows to 237 Np materials, demonstrating the utility of MAMA analysis for a wider breadth of nuclear materials than previously reported. In the open literature, this study is the first time that these microanalytical techniques were applied to 237 Np materials to this degree.

MAMA↗

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↗

Subgrain geometrically necessary dislocation density mapping in spalled Ta in three dimensions

Here, in this work, the intragranular distributions of the densities of geometrically necessary dislocations (GNDs) generated in a polycrystalline tantalum (Ta) sample loaded in shock compression to spall are investigated. The sample studied is a highly resolved three-dimensional (3D) map of the orientations, grain boundary morphology, and voids across a 0.276mm 3 (with sample dimensions of 688.5µm x 669µm x 600µm) region of the microstructure about the spall plane obtained by TriBeam tomography. We present a formulation for calculating the crystallographic screw and edge character GND densities at a material point within body centered cubic (BCC) crystals when provided any number of potential slip systems or modes. Combining the 3D characterization, GND formulation, and a sample containing ≈ 6000 grains, we sought statistically relevant correlations between GND density per grain with grain size, grain orientation, void formation, and proximity from the spall plane. We show that GND density per grain increases with decreasing distance from the spall plane. Most grains had heterogeneously distributed GND densities but were free of any detectable voids, at least within the ≈ 1.5µm resolution of the characterization technique. It is found that the GND density for most grains is insensitive to grain size and grain orientation. Only the fraction of grains with the highest GND densities tended to have a [001] parallel to shock-direction orientation. Intragranular 3D maps of GNDs reveals that grains that contain voids have high GND density concentrations in the intragranular region immediately surrounding the void. The local intragranular GND mapping also showed voids that are connected to grain boundaries preferred to grow into the grain with the larger size, regardless of orientation.

36 MATERIALS SCIENCE↗

Fission fragment decay simulations with the CGMF code

The CGMF code implements the Hauser-Feshbach statistical nuclear reaction model to follow the de-excitation of fission fragments by successive emissions of prompt neutrons and γ rays. The Monte Carlo technique is used to facilitate the analysis of complex distributions and correlations among the prompt fission observables. Starting from initial configurations for the fission fragments in mass, charge, kinetic energy, excitation energy, spin, and parity, Y(A,Z,KE,U,J,π), CGMF samples neutron and γ-ray probability distributions at each stage of the decay process, conserving energy, spin and parity. Nuclear structure and reaction input data from the RIPL3 library are used to describe fission fragment properties and decay probabilities. Characteristics of prompt fission neutrons, prompt fission γ rays, and independent fission yields can be studied consistently. Correlations in energy, angle and multiplicity among the emitted neutrons and γ rays can be easily analyzed as a function of the emitting fragments. Finally, in its current version, CGMF is limited to certain spontaneous fission reactions on Pu and Cf isotopes, as well as neutron-induced fission reactions up to 20 MeV on U, Np, and Pu isotopes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

How Many Elements Matter?

Some studies of stars' multielement abundance distributions suggest at least 5–7 significant dimensions, but others show that many elemental abundances can be predicted to high accuracy from [Fe/H] and [Mg/Fe] (or [Fe/H] and age) alone. We show that both propositions can be, and are, simultaneously true. We adopt a machine-learning technique known as normalizing flow to reconstruct the probability distribution of Milky Way disk stars in the space of 15 elemental abundances measured by APOGEE. Conditioning on T eff and $\mathrm{log}\,g$ minimizes the differential systematics. After further conditioning on [Fe/H] and [Mg/Fe], the residual scatter for most abundances is σ [X/H] ≲ 0.02 dex, consistent with APOGEE's reported statistical uncertainties of ~0.01–0.015 dex and intrinsic scatter of 0.01–0.02 dex. Despite the small scatter, residual abundances display clear correlations between elements, which we show are too large to be explained by measurement uncertainties or by the finite sampling noise. We must condition on at least seven elements to reduce the correlations to a level consistent with the observational uncertainties. Our results demonstrate that cross-element correlations are a much more sensitive probe of a hidden structure than dispersion, and they can be measured precisely in a large sample even if the star-by-star measurement noise is comparable to the intrinsic scatter. We conclude that many elements have an independent story to tell, even for the mundane disk stars and elements produced by the core-collapse and Type Ia supernovae. The only way to learn these lessons is to measure the abundances directly, and not merely infer them.

79 ASTRONOMY AND ASTROPHYSICS↗

Laser-initiated p – 11 B fusion reactions in petawatt high-repetition-rate laser facilities

Driving of the nuclear fusion reaction p + 11 B → 3α + 8.7 MeV under laboratory conditions by interaction between high-power laser pulses and matter has become a popular field of research, owing to its numerous potential applications: as an alternative to deuterium–tritium for fusion energy production, astrophysics studies, and alpha-particle generation for medical treatment. One possible scheme for laser-driven p– 11 B reactions is to direct a beam of laser-accelerated protons onto a boron (B) sample (the so-called “pitcher-catcher” scheme). This technique has been successfully implemented on large high-energy lasers, yielding hundreds of joules per shot at low repetition. We present here a complementary approach, exploiting the high repetition rate of the VEGA III petawatt laser at CLPU (Spain), aiming at accumulating results from many interactions at much lower energy, to provide better control of the parameters and the statistics of the measurements. Despite a moderate energy per pulse, our experiment allowed exploration of the laser-driven fusion process with tens (up to hundreds) of laser shots. The experiment provided a clear signature of the reactions involved and of the fusion products, accumulated over many shots, leading to an improved optimization of the diagnostics for experimental campaigns of this type. In this paper, we discuss the effectiveness of laser-driven p– 11 B fusion in the pitcher–catcher scheme, at a high repetition rate, addressing the challenges of this experimental scheme and highlighting its critical aspects. Our proposed methodology allows evaluation of the performance of this scheme for laser-driven alpha particle production and can be adapted to high-repetition-rate laser facilities with higher energy and intensity.

Alpha particles↗

The completed SDSS-IV extended Baryon Oscillation Spectroscopic Survey: growth rate of structure measurement from cosmic voids

ABSTRACT We present a void clustering analysis in configuration-space using the completed Sloan Digital Sky Survey IV (SDSS-IV) extended Baryon Oscillation Spectroscopic Survey (eBOSS) DR16 samples. These samples consist of Luminous Red Galaxies (LRGs) combined with the high-redshift tail of the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS) DR12 CMASS galaxies (called as LRG+CMASS sample), Emission Line Galaxies (ELGs), and quasars (QSOs). We build void catalogues from the three eBOSS DR16 samples using a ZOBOV-based algorithm, providing 2814 voids, 1801 voids, and 4347 voids in the LRG+CMASS, ELG, and QSO samples, respectively, spanning the redshift range 0.6 < z < 2.2. We measure the redshift space distortions around voids using the anisotropic void-galaxy cross-correlation function and we extract the distortion parameter β. We test the methodology on realistic simulations before applying it to the data, and we investigate all our systematic errors on these mocks. We find βLRG(z = 0.74) = 0.415 ± 0.087, βELG(z = 0.85) = 0.665 ± 0.125 and βQSO(z = 1.48) = 0.313 ± 0.134, for the LRG+CMASS, ELG, and QSO sample, respectively. The quoted errors include systematic and statistical contributions. In order to convert our measurements in terms of the growth rate fσ8, we use consensus values of linear bias from the eBOSS DR16 companion papers, resulting in the following constraints: fσ8(z = 0.74) = 0.50 ± 0.11, fσ8(z = 0.85) = 0.52 ± 0.10, and fσ8(z = 1.48) = 0.30 ± 0.13. Our measurements are consistent with other measurements from eBOSS DR16 using conventional clustering techniques.

79 ASTRONOMY AND ASTROPHYSICS↗

Testing and signal identification for two-sample high-dimensional covariances via multi-level thresholding

The paper considers testing and signal identification for covariance matrices from two populations of marginally sub-Gaussian distributed. A multi-level thresholding procedure is proposed for testing the equality of two high-dimensional covariance matrices, which is designed to detect sparse and faint differences between the covariances. A novel -statistic composition is developed to establish the asymptotic distribution of the thresholding statistics in conjunction with the matrix blocking and the coupling techniques. It is shown that the proposed test is more powerful than the existing tests in detecting sparse and weak signals in covariances. Multiple testing procedures are constructed to discover different covariances and the sub-groups of variables with different covariance structures between the two populations. The proposed procedures are based on the multi-level thresholding test, which are able to control the false discovery proportion () with high power. Simulation experiments and a case study on the returns of the S&P 500 stocks before and after the COVID-19 pandemic are conducted to demonstrate and compare the utilities of the proposed methods.

97 MATHEMATICS AND COMPUTING↗

Air, surface, and wastewater surveillance of SARS-CoV-2; a multimodal evaluation of COVID-19 detection in a built environment

Environmental surveillance of infectious organisms holds tremendous promise to reduce human-to-human transmission in indoor spaces through early detection. In this study we determined the applicability and limitations of wastewater, indoor high-touch surfaces, in-room air, and rooftop exhaust air sampling methods for detecting SARS-CoV-2 in a real world building occupied by residents recently diagnosed with COVID-19. We concurrently examined the results of three 24-hour environmental surveillance techniques, indoor surface sampling, exhaust air sampling and wastewater surveillance, to the known daily census fluctuations in a COVID-19 isolation dormitory. Additionally, we assessed the ability of aerosol samplers placed in the large volume lobby to detect SARS-CoV-2 multiple times per day. Our research reveals an increase in the number of individuals confirmed positive with COVID-19 as well as their estimated human viral load to be associated with statistically significant increases in viral loads detected in rooftop exhaust aerosol samples (p = 0.0413), wastewater samples (p = 0.0323,), and indoor high-touch surfaces (p < 0.001)). We also report that the viral load detected in lobby aerosol samples was statistically higher in samples collected during presence of occupants whose COVID-19 diagnostic tests were confirmed positive via qPCR compared to periods when the lobby was occupied by either contact-traced (suspected positive) individuals or during unoccupied periods (p = 0.0314 and <2e–16). We conclude that each daily (24h) surveillance method, rooftop exhaust air, indoor high-touch surfaces, and wastewater, provide useful detection signals for building owner/operator(s). Furthermore, we demonstrate that exhaust air sampling can provide spatially resolved signals based upon ventilation exhaust zones. Additionally, we find that indoor lobby air sampling can provide temporally resolved signals useful during short duration sampling periods (e.g., 2-4 hours) even with intermittent occupancy by occupants diagnosed with COVID-19.

60 APPLIED LIFE SCIENCES↗

Analysis of the Challenges in Developing Sample-Based Multi-fidelity Estimators for Non-deterministic Models

Multifidelity (MF) uncertainty quantification (UQ) seeks to leverage and fuse information from a collection of models to achieve greater statistical accuracy with respect to a single-fidelity counterpart, while maintaining an efficient use of computational resources. Despite many recent advancements in MF UQ, several challenges remain and these often limit its practical impact in certain application areas. In this manuscript, we focus on the challenges introduced by nondeterministic models to sampling MF UQ estimators. Nondeterministic models produce different responses for the same inputs, which means their outputs are effectively noisy. MF UQ is complicated by this noise since many state-of-the-art approaches rely on statistics, e.g., the correlation among models, to optimally fuse information and allocate computational resources. Here, we demonstrate how the statistics of the quantities of interest, which impact the design, effectiveness, and use of existing MF UQ techniques, change as functions of the noise. With this in hand, we extend the unifying approximate control variate framework to account for nondeterminism, providing for the first time a rigorous means of comparing the effect of nondeterminism on different multifidelity estimators and analyzing their performance with respect to one another. Numerical examples are presented throughout the manuscript to illustrate and discuss the consequences of the presented theoretical results.

97 MATHEMATICS AND COMPUTING↗

Rare Event Sampling Improves Mercury Instability Statistics

Due to the chaotic nature of planetary dynamics, there is a non-zero probability that Mercury's orbit will become unstable in the future. Previous efforts have estimated the probability of this happening between 3 and 5 billion years in the future using a large number of direct numerical simulations with an N-body code, but were not able to obtain accurate estimates before 3 billion years in the future because Mercury instability events are too rare. In this paper we use a new rare-event sampling technique, Quantile Diffusion Monte Carlo (QDMC), to estimate that the probability of a Mercury instability event in the next 2 billion years is approximately 10 –4 in the REBOUND N-body code. We show that QDMC provides unbiased probability estimates at a computational cost of up to 100 times less than direct numerical simulation. QDMC is easy to implement and could be applied to many problems in planetary dynamics in which it is necessary to estimate the probability of a rare event.

79 ASTRONOMY AND ASTROPHYSICS↗

Assessment of Potential Ergonomic Injury Risk in the Nuclear Material Processing Glovebox Environment [Capstone Project]

Gloveboxes are isolation barriers that are used within many different industries such as pharmaceuticals, electronic parts fabrication, nuclear, and biological. The research on glovebox ergonomics is currently limited with few ergonomic professionals that focusing exclusively on glovebox working environments. There is existing documentation on occupational injuries, such as musculoskeletal disorders (MSD), sustained as a direct result of working in gloveboxes. The main elements that pose ergonomic risks to glovebox workers are operational repetition, duration, force, vibration, lifting heavy (more than 15lbs with two hands) objects, and awkward postures. This paper examines a small sample of the potential causal or risk factors that lead to the ergonomic injuries. A meta-analysis utilizing a random effect model is used to examine data from several studies which focus on, dexterity and strength changes as result of glove thickness, and robotic assistive technology as a means to improve postural mechanics. With the meta-analysis technique, similar data sets from different research studies can be coalesced into a single weighted, statistically significant result for subsequent consideration. The results from the analysis show that increased glove thickness results in decreased dexterity for operators thus increasing ergonomic risk factors such as task duration. It is also shown that glove thickness decreases grip strength but a similar decrease in pinch strength is not definitively demonstrated. With decreased grip strength, operators will need to exert more force (a known ergonomic risk factor) on processing tools, etc. during operations thus increasing the risk of ergonomic injury. The use of robotic assistive technology as a means to improve operator posture (risk factor) was also examined. Although it may be intuitively assumed that human-robotic collaboration would be beneficial in reducing risk factors, the result from this study’s analysis was not statistically significant. It is inferred that with additional directed research on this topic another study/analysis could be statistically significant demonstrating the benefits of the technology.

99 GENERAL AND MISCELLANEOUS↗

Fast Grain Mapping with Sub-Nanometer Resolution Using 4D-STEM with Grain Classification by Principal Component Analysis and Non-Negative Matrix Factorization

High-throughput grain mapping with sub-nanometer spatial resolution is demonstrated using scanning nanobeam electron diffraction (also known as 4D scanning transmission electron microscopy, or 4D-STEM) combined with high-speed direct-electron detection. An electron probe size down to 0.5 nm in diameter is used and the sample investigated is a gold–palladium nanoparticle catalyst. Computational analysis of the 4D-STEM data sets is performed using a disk registration algorithm to identify the diffraction peaks followed by feature learning to map the individual grains. Two unsupervised feature learning techniques are compared: principal component analysis (PCA) and non-negative matrix factorization (NNMF). The characteristics of the PCA versus NNMF output are compared and the potential of the 4D-STEM approach for statistical analysis of grain orientations at high spatial resolution is discussed.

47 OTHER INSTRUMENTATION↗

On the uncertainty of long-period return values of extreme daily precipitation

Methods for calculating return values of extreme precipitation and their uncertainty are compared using daily precipitation rates over the Western U.S. and Southwestern Canada from a large ensemble of climate model simulations. The roles of return-value estimation procedures and sample size in uncertainty are evaluated for various return periods. We compare two different generalized extreme value (GEV) parameter estimation techniques, namely L-moments and maximum likelihood (MLE), as well as empirical techniques. Even for very large datasets, confidence intervals calculated using GEV techniques are narrower than those calculated using empirical methods. Furthermore, the more efficient L-moments parameter estimation techniques result in narrower confidence intervals than MLE parameter estimation techniques at small sample sizes, but similar best estimates. It should be noted that we do not claim that either parameter fitting technique is better calibrated than the other to estimate long period return values. While a non-stationary MLE methodology is readily available to estimate GEV parameters, it is not for the L-moments method. Comparison of uncertainty quantification methods are found to yield significantly different estimates for small sample sizes but converge to similar results as sample size increases. Finally, practical recommendations about the length and size of climate model ensemble simulations and the choice of statistical methods to robustly estimate long period return values of extreme daily precipitation statistics and quantify their uncertainty.

54 ENVIRONMENTAL SCIENCES↗

A probabilistic inverse prediction method for predicting plutonium processing conditions

In the past decade, nuclear chemists and physicists have been conducting studies to investigate the signatures associated with the production of special nuclear material (SNM). In particular, these studies aim to determine how various processing parameters impact the physical, chemical, and morphological properties of the resulting special nuclear material. By better understanding how these properties relate to the processing parameters, scientists can better contribute to nuclear forensics investigations by quantifying their results and ultimately shortening the forensic timeline. This paper aims to statistically analyze and quantify the relationships that exist between the processing conditions used in these experiments and the various properties of the nuclear end-product by invoking inverse methods. In particular, these methods make use of Bayesian Adaptive Spline Surface models in conjunction with Bayesian model calibration techniques to probabilistically determine processing conditions as an inverse function of morphological characteristics. Not only does the model presented in this paper allow for providing point estimates of a sample of special nuclear material, but it also incorporates uncertainty into these predictions. This model proves sufficient for predicting processing conditions within a standard deviation of the observed processing conditions, on average, provides a solid foundation for future work in predicting processing conditions of particles of special nuclear material using only their observed morphological characteristics, and is generalizable to the field of chemometrics for applicability across different materials.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Statistically elucidated responses from low-signal contrast mechanisms in ultrafast electron microscopy

The emergence of ultrafast electron microscopy (UEM) has enabled the discovery of strongly correlated dynamic mechanisms, including electron–phonon coupling, structural phase transitions, thermal transport, and electromagnetic deflection. Most UEM systems operate stroboscopically, meaning that the technique is susceptible to artifacts, mistakes, and misinterpretation of the data due to extensive experimental effort. In contrast to the ultrafast designation, data acquisition is extraordinarily slow because the electron beam has significantly reduced signal compared to traditional transmission electron microscopy due to pulsing the electron beam. Consequently, the sample may drift, tilt, or undergo irreversible structural changes that are independent of the time-resolved dynamics throughout the experimental time frame. Furthermore, these datasets require significant user interpretation that can be problematic when proper controls are not implemented thoroughly. Here, we demonstrate a new algorithm designed to separate ultrafast structural dynamics from long-term artifacts using a LiNbO 3 sample experiencing electrically driven surface acoustic wave propagation. Additionally, we provide examples of the impact of user bias when analyzing the data and provide a methodology, which enables the extraction of time-resolved responses when the image signal is extraordinarily low. Overall, the goal of this publication is to provide methods that validate the experimental results and reduce researcher biases during UEM data interpretation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Combining compositional data sets introduces error in covariance network reconstruction

Microbial communities are diverse biological systems that include taxa from across multiple kingdoms of life. Notably, interactions between bacteria and fungi play a significant role in determining community structure. However, these statistical associations across kingdoms are more difficult to infer than intra-kingdom associations due to the nature of the data involved using standard network inference techniques. We quantify the challenges of cross-kingdom network inference from both theoretical and practical points of view using synthetic and real-world microbiome data. We detail the theoretical issue presented by combining compositional data sets drawn from the same environment, e.g. 16S and ITS sequencing of a single set of samples, and we survey common network inference techniques for their ability to handle this error. We then test these techniques for the accuracy and usefulness of their intra- and inter-kingdom associations by inferring networks from a set of simulated samples for which a ground-truth set of associations is known. We show that while the two methods mitigate the error of cross-kingdom inference, there is little difference between techniques for key practical applications including identification of strong correlations and identification of possible keystone taxa (i.e. hub nodes in the network). Furthermore, we identify a signature of the error caused by transkingdom network inference and demonstrate that it appears in networks constructed using real-world environmental microbiome data.

59 BASIC BIOLOGICAL SCIENCES↗

Combining compositional data sets introduces error in covariance network reconstruction

Microbial communities are diverse biological systems that include taxa from across multiple kingdoms of life. Notably, interactions between bacteria and fungi play a significant role in determining community structure. However, these statistical associations across kingdoms are more difficult to infer than intra-kingdom associations due to the nature of the data involved using standard network inference techniques. We quantify the challenges of cross-kingdom network inference from both theoretical and practical points of view using synthetic and real-world microbiome data. We detail the theoretical issue presented by combining compositional data sets drawn from the same environment, e.g. 16S and ITS sequencing of a single set of samples, and we survey common network inference techniques for their ability to handle this error. We then test these techniques for the accuracy and usefulness of their intra- and interkingdom associations by inferring networks from a set of simulated samples for which a ground-truth set of associations is known. We show that while the two methods mitigate the error of cross-kingdom inference, there is little difference between techniques for key practical applications including identification of strong correlations and identification of possible keystone taxa (i.e. hub nodes in the network). Furthermore, we identify a signature of the error caused by transkingdom network inference and demonstrate that it appears in networks constructed using real-world environmental microbiome data.

59 BASIC BIOLOGICAL SCIENCES↗