Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “testing 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 109 records · Page 6

Grid-based minimization at scale: Feldman-Cousins corrections for light sterile neutrino search

High Energy Physics (HEP) experiments generally employ sophisticated statistical methods to present results in searches of new physics. In the problem of searching for sterile neutrinos, likelihood ratio tests are applied to short-baseline neutrino oscillation experiments to construct confidence intervals for the parameters of interest. The test statistics of the form Δχ2 is often used to form the confidence intervals, however, this approach can lead to statistical inaccuracies due to the small signal rate in the region-of-interest. In this paper, we present a computational model for the computationally expensive Feldman-Cousins corrections to construct a statistically accurate confidence interval for neutrino oscillation analysis. The program performs a grid-based minimization over oscillation parameters and is written in C++. Our algorithms make use of vectorization through Eigen3, yielding a single-core speed-up of 350 compared to the original implementation, and achieve MPI data parallelism by employing DIY. We demonstrate the strong scaling of the application at High-Performance Computing (HPC) sites. We utilize HDF5 along with HighFive to write the results of the calculation to file.

Wospakrik, Marianette↗

First Searches for Axions and Axionlike Particles with the LUX Experiment

The first searches for axions and axionlike particles with the Large Underground Xenon experiment are presented. Under the assumption of an axioelectric interaction in xenon, the coupling constant between axions and electrons g Ae is tested using data collected in 2013 with an exposure totaling 95 live days ×118 kg. A double-sided, profile likelihood ratio statistic test excludes g Ae larger than 3.5 × 10 –12 (90% C.L.) for solar axions. Assuming the Dine-Fischler-Srednicki-Zhitnitsky theoretical description, the upper limit in coupling corresponds to an upper limit on axion mass of 0.12 eV/c 2 , while for the Kim-Shifman-Vainshtein-Zhakharov description masses above 36.6 eV/c 2 are excluded. For galactic axionlike particles, values of g Ae larger than 4.2 × 10 –13 are excluded for particle masses in the range 1–16 keV/c 2 . As a result, these are the most stringent constraints to date for these interactions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

SU-F-I-80: Correction for Bias in a Channelized Hotelling Model Observer Caused by Temporally Variable Non-Stationary Noise

Purpose: Application of a channelized Hotelling model observer (CHO) over a wide range of x-ray angiography detector target dose (DTD) levels demonstrated substantial bias for conditions yielding low detectability indices (d’), including low DTD and small test objects. The purpose of this work was to develop theory and methods to correct this bias. Methods: A hypothesis was developed wherein the measured detectability index (d’b) for a known test object is positively biased by temporally variable non-stationary noise in the images. Hotelling’s T2 test statistic provided the foundation for a mathematical theory which accounts for independent contributions to the measured d’b value from both the test object (d’o) and non-stationary noise (d’ns). Experimental methods were developed to directly estimate d’o by determining d’ns and subtracting it from d’b, in accordance with the theory. Specifically, d’ns was determined from two sets of images from which the traditional test object was withheld. This method was applied to angiography images with DTD levels in the range 0 to 240 nGy and for disk-shaped iodine-based contrast targets with diameters 0.5 to 4.0 mm. Results: Bias in d’ was evidenced by d’b values which exceeded values expected from a quantum limited imaging system and decreasing object size and DTD. d’ns increased with decreasing DTD, reaching a maximum of 2.6 for DTD = 0. Bias-corrected d’o estimates demonstrated sub-quantum limited performance of the x-ray angiography for low DTD. Findings demonstrated that the source of non-stationary noise was detector electronic readout noise. Conclusion: Theory and methods to estimate and correct bias in CHO measurements from temporally variable non-stationary noise were presented. The temporal non-stationary noise was shown to be due to electronic readout noise. This method facilitates accurate estimates of d’ values over a large range of object size and detector target dose.

60 APPLIED LIFE SCIENCES↗

CORRELATING THE INTERSTELLAR MAGNETIC FIELD WITH PROTOSTELLAR JETS AND ITS SOURCES

This paper combines new CCD polarimetric data with previous information about protostellar objects in a search for correlations involving the interstellar magnetic field (ISMF). Specifically, we carried out an optical polarimetric study of a sample of 28 fields of 10' Multiplication-Sign 10' located in the neighborhood of protostellar jets and randomly spread over the Galaxy. The polarimetry of a large number of field stars is used to estimate both the average and dispersion of the ISMF direction in each region. The results of the applied statistical tests are as follows. Concerning the alignment between the jet direction and the ISMF, the whole sample does not show alignment. There is, however, a statistically significant alignment for objects of Classes 0 and I. Regarding the ISMF dispersion, our sample presents values slightly larger for regions containing T Tauri objects than for those harboring younger protostars. Moreover, the ISMF dispersion in regions containing high-mass objects tends to be larger than in those including only low-mass protostars. In our sample, the mean interstellar polarization as a function of the average interstellar extinction in a region reaches a maximum value around 3% for A(V) = 5, after which it decreases. Our data also show a clear correlation of the mean value of the interstellar polarization with the dispersion of the ISMF: the larger the dispersion, the smaller the polarization. Based on a comparison of our and previous results, we suggest that the dispersion in regions forming stars is larger than in quiescent regions.

79 ASTRONOMY AND ASTROPHYSICS↗

DefectTrack: a deep learning-based multi-object tracking algorithm for quantitative defect analysis of in-situ TEM videos in real-time

Abstract In-situ irradiation transmission electron microscopy (TEM) offers unique insights into the millisecond-timescale post-cascade process, such as the lifetime and thermal stability of defect clusters, vital to the mechanistic understanding of irradiation damage in nuclear materials. Converting in-situ irradiation TEM video data into meaningful information on defect cluster dynamic properties (e.g., lifetime) has become the major technical bottleneck. Here, we present a solution called the DefectTrack , the first dedicated deep learning-based one-shot multi-object tracking (MOT) model capable of tracking cascade-induced defect clusters in in-situ TEM videos in real-time. DefectTrack has achieved a Multi-Object Tracking Accuracy (MOTA) of 66.43% and a Mostly Tracked (MT) of 67.81% on the test set, which are comparable to state-of-the-art MOT algorithms. We discuss the MOT framework, model selection, training, and evaluation strategies for in-situ TEM applications. Further, we compare the DefectTrack with four human experts in quantifying defect cluster lifetime distributions using statistical tests and discuss the relationship between the material science domain metrics and MOT metrics. Our statistical evaluations on the defect lifetime distribution suggest that the DefectTrack outperforms human experts in accuracy and speed.

42 ENGINEERING↗

Emissions and Flame Structure Assessment of Aeroengine Micromixing Injectors for Lean Direct Injection of Hydrogen and Hydrogen/Natural Gas Blends

Abstract The present study focuses on whether adapting low emissions liquid fuel injection concepts can be applied to gaseous mixtures. The goal is to develop an injector with specific geometric characteristics that will achieve the lowest emissions. This is conducted by evaluating emissions and flame diagnostics on a set of test hardware that has been developed for combustion of pure hydrogen and hydrogen/methane mixtures. The current test hardware adapts lean direct injection concepts from aeroengine injectors to a ground-based engine based on target flows associated with a Solar Turbines Industrial Engine gas turbine. A statistical test plan was developed to assess these geometric characteristics using a three level Box Behnken design. The test plan investigates the effect of injector flow splits, air swirl, fuel swirl, pressure drop, preheat temperature, fuel composition, and the adiabatic flame temperature on the flame structure and emissions performance. A total of 16 injectors were considered for a total of 688 test points for emission and imaging data. Emissions data were reported on a mass basis in ng/J and chemical reactor network analysis were conducted to further interpret emissions behavior. Images were also captured to analyze the flame structure features of interest including reaction volume, flame length, flame spread, and flame width. With these results, a predictive model for emissions was made to conclude the most optimal injector configuration is to have the highest level of air split (more volume of air through the injector plate versus the injector) and air swirl angle for both methane and hydrogen gas mixtures. The predictive emissions level for all of these cases were below 10 ng/J.

08 HYDROGEN↗

Cervical Precancers and Cancers Attributed to HPV Types by Race and Ethnicity: Implications for Vaccination, Screening, and Management

Abstract Background Racial and ethnic variations in attribution of cervical precancer and cancer to human papillomavirus (HPV) types may result in different HPV vaccine protection, screening test coverage, and clinical management. Methods Pooling data from 7 US studies, we calculated the proportional attribution of precancers and cancers to HPV types using HPV DNA typing from diagnosis. All statistical tests were 2-sided. Results For all racial and ethnic groups, most cases of cervical intraepithelial neoplasia grade 3 (CIN3) (84.2%-90.8% of 5526) and squamous cell carcinoma (SCC) (90.4%-93.8% of 1138) were attributed to types targeted by the 9-valent vaccine. A higher proportion of CIN3s were attributed to nonvaccine HPV types among non-Hispanic Black women (15.8%) compared with non-Hispanic Asian or Pacific Islander (9.7%; P = .002), non-Hispanic White (9.2%; P < .001), and Hispanic (11.3%; P = .004) women. The proportion of SCCs attributed to 9-valent types was similar by race and ethnicity (P = .80). A higher proportion of CIN3s were attributed to nonvaccine HPV35 among non-Hispanic Black (9.0%) compared with non-Hispanic Asian or Pacific Islander (2.2%), non-Hispanic White (2.5%), and Hispanic (3.0%; all P < .001) women. Compared with CIN3, the proportion of SCCs attributed to HPV35 among non-Hispanic Black women (3.2%) was lower and closer to other groups (0.3%-2.1%; P = .70). Conclusion The 9-valent HPV vaccine will prevent nearly all cervical precancers and invasive cancers among major racial and ethnic groups in the United States. Adding HPV35 to vaccines could prevent a small percentage of CIN3s and SCCs, with greater potential impact for CIN3s among Black women. HPV screening tests target high-risk HPV types, including HPV35. Future genotyping triage strategies could consider the importance of HPV35- and other HPV16-related types.

Mix, Jacqueline (ORCID:0000000164187706)↗

Sandia Fracture Challenge 3: detailing the Sandia Team Q failure prediction strategy

The third Sandia Fracture Challenge illustrated the geometric and material uncertainties presented by modern additive manufacturing techniques. Tasked with the challenge of predicting failure of a complex additively-manufactured geometry made of 316L stainless steel, we combined a rigorous material calibration scheme with a number of statistical assessments of problem uncertainties. Particularly, we used optimization techniques to calibrate a rate-dependent and anisotropic Hill plasticity model to represent material deformation coupled with a damage model driven by void growth and nucleation. Through targeted simulation studies we assessed the influence of internal voids and surface flaws on the specimens of interest in the challenge which guided our material modeling choices. Employing the Kolmogorov–Smirnov test statistic, we developed a representative suite of simulations to account for the geometric variability of test specimens and the variability introduced by material parameter uncertainty. This method allowed the team to successfully predict the failure mode of the experimental test population as well as the global response with a high degree of accuracy.

316L stainless steel↗

The clustering of the SDSS-IV extended Baryon Oscillation Spectroscopic Survey quasar sample: testing observational systematics on the Baryon Acoustic Oscillation measurement

ABSTRACT Baryon Acoustic Oscillations (BAO) are considered to be a very robust standard ruler against various systematics. This premise has been tested against observational systematics, but not to the level required for the next generation of galaxy surveys such as the Dark Energy Spectroscopic Instrument (DESI) and Euclid. In this paper, we investigate the effect of observational systematics on the BAO measurement of the final sample of quasars from the extended Baryon Oscillation Spectroscopic Survey Data Release 16 in order to prepare and hone a similar analysis for upcoming surveys. We employ catalogues with various treatments of imaging systematic effects using linear and neural network-based non-linear approaches and consider how the BAO measurement changes. We also test how the variations to the BAO fitting model respond to the observational systematics. As expected, we confirm that the BAO measurements obtained from the DR16 quasar sample are robust against imaging systematics well within the statistical error, while reporting slightly modified constraints that shift the line-of-sight BAO signal by less than 1.1 per cent. We use realistic simulations with similar redshift and angular distributions as the DR16 sample to conduct statistical tests for validating the pipeline, quantifying the significance of differences, and estimating the expected bias on the BAO scale in future high-precision data sets. Although we find a marginal impact for the eBOSS QSO data, the work presented here is of vital importance for constraining the nature of dark energy with the BAO feature in the new era of big data cosmology.

79 ASTRONOMY AND ASTROPHYSICS↗

Filtering a statistically exactly solvable test model for turbulent tracers from partial observations

A statistically exactly solvable model for passive tracers is introduced as a test model for the authors' Nonlinear Extended Kalman Filter (NEKF) as well as other filtering algorithms. The model involves a Gaussian velocity field and a passive tracer governed by the advection-diffusion equation with an imposed mean gradient. The model has direct relevance to engineering problems such as the spread of pollutants in the air or contaminants in the water as well as climate change problems concerning the transport of greenhouse gases such as carbon dioxide with strongly intermittent probability distributions consistent with the actual observations of the atmosphere. One of the attractive properties of the model is the existence of the exact statistical solution. In particular, this unique feature of the model provides an opportunity to design and test fast and efficient algorithms for real-time data assimilation based on rigorous mathematical theory for a turbulence model problem with many active spatiotemporal scales. Here, we extensively study the performance of the NEKF which uses the exact first and second order nonlinear statistics without any approximations due to linearization. The role of partial and sparse observations, the frequency of observations and the observation noise strength in recovering the true signal, its spectrum, and fat tail probability distribution are the central issues discussed here. The results of our study provide useful guidelines for filtering realistic turbulent systems with passive tracers through partial observations.

54 ENVIRONMENTAL SCIENCES↗

AFIP6-MkII and RERTR-12 Porosity Data Collection and Analysis for Modeling and Simulation

Gathering data for the improvement of nuclear fuel modeling and simulation efforts is the primary driver for this work. Mechanistic models allow for a better understanding of the material on a micro- and macrostructural level while saving time and money over traditional experiment efforts. Historically, summarized data and correlations are the inputs for empirical material models and model validation. When improving these models for nuclear fuels with experimental results, there is a lack of reliable data readily available. Experiments - RERTR-12 and AFIP6-MkII - were conducted to understand the irradiation behavior of metallic U-10Mo monolithic fuels for use in extreme reactor environments such as research reactors like the Advanced Test Reactor (ATR) or the High Flux Isotope Reactor (HFIR). Microstructural characteristics of fission gas pores (FGP) in each experiment are collected using an automated image analysis technique developed at the University of Florida and presented here. A series of statistical tests are performed to explore the reliability of the results, as well as understand where the data is lacking and what future data collection is necessary to provide sufficient information to assist modeling efforts. The focus is on the porosity, pore size, and eccentricity of FGPs formed during irradiation in three AFIP6-MkII samples and one RERTR-12 sample. From the analysis, it is clear there are substantial impacts of fission density on the pore structure, but there also exist also underlying connections between each sample and the behavior observed in the pores. Further analyses of the pre- and post-irradiation microstructure are needed to improve the understanding of these connections. An early method for microstructural data analysis is presented within and is currently being expanded to include other microstructure data.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Bayesian Analysis of Congruence of Core Genes in Prochlorococcus and Synechococcus and Implications on Horizontal Gene Transfer

It is often suggested that horizontal gene transfer is so ubiquitous in microbes that the concept of a phylogenetic tree representing the pattern of vertical inheritance is oversimplified or even positively misleading. “Universal proteins” have been used to infer the organismal phylogeny, but have been criticized as being only the “tree of one percent.” Currently, few options exist for those wishing to rigorously assess how well a universal protein phylogeny, based on a relative handful of well-conserved genes, represents the phylogenetic histories of hundreds of genes. Here, we address this problem by proposing a visualization method and a statistical test within a Bayesian framework. We use the genomes of marine cyanobacteria, a group thought to exhibit substantial amounts of HGT, as a test case. We take 379 orthologous gene families from 28 cyanobacteria genomes and estimate the Bayesian posterior distributions of trees – a “treecloud” – for each, as well as for a concatenated dataset based on putative “universal proteins.” We then calculate the average distance between trees within and between all treeclouds on various metrics and visualize this high-dimensional space with non-metric multidimensional scaling (NMMDS). We show that the tree space is strongly clustered and that the universal protein treecloud is statistically significantly closer to the center of this tree space than any individual gene treecloud. We apply several commonly-used tests for incongruence/HGT and show that they agree HGT is rare in this dataset, but make different choices about which genes were subject to HGT. Our results show that the question of the representativeness of the “tree of one percent” is a quantitative empirical question, and that the phylogenetic central tendency is a meaningful observation even if many individual genes disagree due to the various sources of incongruence.

59 BASIC BIOLOGICAL SCIENCES↗

Detectability of neutrino-signal fluctuations induced by the hadron-quark phase transition in failing core-collapse supernovae

Here we introduce a systematic and quantitative methodology for establishing the presence of neutrino oscillatory signals due to the hadron-quark phase transition (PT) in failing core-collapse supernovae from the observed neutrino event rate in water- or ice-based neutrino detectors. The methodology uses a likelihood ratio in the frequency domain as a test-statistic; it is employed for quantitative analysis of neutrino signals without assuming the frequency, amplitude, starting time, and duration of the PT-induced oscillations present in the neutrino events and thus it is suitable for analyzing neutrino signals from a wide variety of numerical simulations. We test the validity of this method by using a core-collapse simulation of a 17 solar-mass star by Zha et al. [Astrophys. J. 911, 74 (2021) ]. Based on this model, we further report the presence of a PT-induced oscillations quantitatively for a core-collapse supernovae out to a distance of ~10 kpc, ~5 kpc for IceCube and to a distance of ~10 kpc, ~5 kpc, and ~1 kpc for a 0.4 Mt mass water Cherenkov detector. This methodology will aid the investigation of a future galactic supernova and the study of hadron-quark phase in the core of core-collapse supernovae.

79 ASTRONOMY AND ASTROPHYSICS↗

Image Based Characterization of Formal and Informal Neighborhoods in an Urban Landscape

The high rate of global urbanization has resulted in a rapid increase in informal settlements, which can be de ned as unplanned, unauthorized, and/or unstructured housing. Techniques for ef ciently mapping these settlement boundaries can bene t various decision making bodies. From a remote sensing perspective, informal settlements share unique spatial characteristics that distinguish them from other types of structures (e.g., industrial, commercial, and formal residential). These spatial characteristics are often captured in high spatial resolution satellite imagery. We analyzed the role of spatial, structural, and contextual features (e.g., GLCM, Histogram of Oriented Gradients, Line Support Regions, Lacunarity) for urban neighborhood mapping, and computed several low-level image features at multiple scales to characterize local neighborhoods. The decision parameters to classify formal-, informal-, and non-settlement classes were learned under Decision Trees and a supervised classi cation framework. Experiments were conducted on high-resolution satellite imagery from the CitySphere collection, and four different cities (i.e., Caracas, Kabul, Kandahar, and La Paz) with varying spatial characteristics were represented. Overall accuracy ranged from 85% in La Paz, Bolivia, to 92% in Kandahar, Afghanistan. While the disparities between formal and informal neighborhoods varied greatly, many of the image statistics tested proved robust.

58 GEOSCIENCES↗

Observational tests of the black hole area increase law

The black hole area theorem implies that when two black holes merge, the area of the final black hole should be greater than the sum of the areas of the two original black holes. In this study, we examine how this prediction can be tested with gravitational-wave observations of binary black holes. By separately fitting the early inspiral and final ringdown stages, we calculate the posterior distributions for the masses and spins of the two initial and the final black holes. This yields posterior distributions for the change in the area and thus a statistical test of the validity of the area increase law. We illustrate this method with a GW150914-like binary black hole waveform calculated using numerical relativity, and detector sensitivities representative of both the first observing run and the design configuration of Advanced LIGO. We obtain a ~74.6 % probability that the simulated signal is consistent with the area theorem with current sensitivity, improving to ~99.9 % when Advanced LIGO reaches design sensitivity. Lastly, an important ingredient in our test is a method of estimating when the postmerger signal is well fit by a damped sinusoid ringdown waveform.

79 ASTRONOMY AND ASTROPHYSICS↗

Multiscale spatial analysis of fracture arrangement and pattern reconstruction using Ripley's K-function

This work presents novel multiscale spatial data analytics using Ripley's K-function, as a measure of spatial interaction, to study one-dimensional arrangement of fractures. Fracture spatial arrangements are classified into clustered, anticlustered, or indistinguishable from random by testing statistical significance of the calculated Ripley's K-function. Characterizations of fracture arrangements are performed as a function of length scale and position. Analysis of the K-function along the study interval identifies where fracture clustering and anticlustering occur. Here, a simulation technique is also introduced here to statistically reconstruct spatial arrangements and to generate fracture realizations that are spatially similar to the fractures observed in the field. With this simulation technique, one can also fill spatial gaps in fracture measurements where data are absent, unreliable, or unused. Synthetic as well as field-measured 1D fracture datasets are used for testing and demonstration. Methods introduced in this work can be readily applied to fracture datasets observed in outcrops, borehole image logs, and cores.

58 GEOSCIENCES↗

PNNL-Predictive-Phenomics/ProteoMeter

ProteoMeter is a Python package that assists in the statistical analysis of global proteomics, protein post-translation modification (PTM), and limited proteolysis (LiP) data. It contains batch correction, normalization, and statistical testing methods, as well as functions that "roll up" peptide-level data to the single-site level. It has a robust user configuration system, allowing it to flexibly integrate different types of experiment designs. For basic usage, a simple configuration file provides the essential functionality. Advanced users have access to the entire statistical pipeline for fine-tuning analyses. Processed data is easily exported to many common spreadsheet and data-frame formats.

Rozum, Jordan [Pacific Northwest National Lab]↗

Integrating Symbolic and Statistical Methods for Testing Intelligent Systems Applications to Machine Learning and Computer Vision

Embedded intelligent systems ranging from tiny im- plantable biomedical devices to large swarms of autonomous un- manned aerial systems are becoming pervasive in our daily lives. While we depend on the flawless functioning of such intelligent systems, and often take their behavioral correctness and safety for granted, it is notoriously difficult to generate test cases that expose subtle errors in the implementations of machine learning algorithms. Hence, the validation of intelligent systems is usually achieved by studying their behavior on representative data sets, using methods such as cross-validation and bootstrapping.In this paper, we present a new testing methodology for studying the correctness of intelligent systems. Our approach uses symbolic decision procedures coupled with statistical hypothesis testing to. We also use our algorithm to analyze the robustness of a human detection algorithm built using the OpenCV open-source computer vision library. We show that the human detection implementation can fail to detect humans in perturbed video frames even when the perturbations are so small that the corresponding frames look identical to the naked eye.

Jha, Sumit Kumar↗