Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 541 records · Page 30

Interpretable Machine Learning for Molecular Biosignatures: a Novel Single-Sample Feature Importance Method That Is Sensitive To Statistical Interactions

Isotope ratio mass spectrometry (IRMS) of volatiles (e.g., CO 2 ) promises to be a powerful tool for potential biosignature detection for future missions to ocean worlds (OW) such as Europa and Enceladus. Machine learning (ML) methods for IRMS data could enable science autonomy by onboard prediction of seawater chemistry and biosignature presence. However, ML models are likely to be complex and involve statistical interactions between features (variables), which can make predictions seem opaque and enigmatic. For ML predictions as significant as extraterrestrial biosignatures, we must place extraordinary confidence in models. It is therefore essential that these models make interpretable predictions (i.e., human-understandable) and include false-prediction diagnostics. We achieve high accuracy and interpretability in ML biosignature and seawater chemistry models for OW through a nearest-neighbors feature selection tool that detects statistical interactions between predictors, constructs interaction networks for visualization of selected features working together to make a prediction, and reports single-sample feature importance scores for false-detection diagnostics. Here we develop a novel single-sample nearest-neighbors projected distance regression(ssNPDR) feature selection method that improves upon existing single-sample algorithms through the inclusion of statistical interactions while providing false-prediction diagnostics for ML models.

geochemistry↗

Statistical Characterization of Medium-Duty Electric Vehicle Drive Cycles: Preprint

In an effort to help commercialize technologies for electric vehicles (EVs) through deployment and demonstration projects, the U.S. Department of Energy’s (DOE's) American Recovery and Reinvestment Act (ARRA) provided funding to participating U.S. companies to cover part of the cost of purchasing new EVs. Within the medium- and heavy-duty commercial vehicle segment, both Smith Electric Newton and and Navistar eStar vehicles qualified for such funding opportunities. In an effort to evaluate the performance characteristics of the new technologies deployed in these vehicles operating under real world conditions, data from Smith Electric and Navistar medium-duty EVs were collected, compiled, and analyzed by the National Renewable Energy Laboratory's (NREL) Fleet Test and Evaluation team over a period of 3 years. More than 430 Smith Newton EVs have provided data representing more than 150,000 days of operation. Similarly, data have been collected from more than 100 Navistar eStar EVs, resulting in a comparative total of more than 16,000 operating days. Combined, NREL has analyzed more than 6 million kilometers of driving and 4 million hours of charging data collected from commercially operating medium-duty electric vehicles in various configurations. In this paper, extensive duty-cycle statistical analyses are performed to examine and characterize common vehicle dynamics trends and relationships based on in-use field data. The results of these analyses statistically define the vehicle dynamic and kinematic requirements for each vehicle, aiding in the selection of representative chassis dynamometer test cycles and the development of custom drive cycles that emulate daily operation. In this paper, the methodology and accompanying results of the duty-cycle statistical analysis are presented and discussed. Results are presented in both graphical and tabular formats illustrating a number of key relationships between parameters observed within the data set that relate to medium duty EVs.

33 ADVANCED PROPULSION SYSTEMS↗

Ideas for Effective Communication of Statistical Results

Effective presentation of statistical results to those with less statistical training, including managers and decision-makers requires planning, anticipation and thoughtful delivery. Here are several recommendations for effectively presenting statistical results.

97 MATHEMATICS AND COMPUTING↗

Statistical Characterization of Medium-Duty Electric Vehicle Drive Cycles

In an effort to help commercialize technologies for electric vehicles (EVs) through deployment and demonstration projects, the U.S. Department of Energy's (DOE's) American Recovery and Reinvestment Act (ARRA) provided funding to participating U.S. companies to cover part of the cost of purchasing new EVs. Within the medium- and heavy-duty commercial vehicle segment, both Smith Electric Newton and and Navistar eStar vehicles qualified for such funding opportunities. In an effort to evaluate the performance characteristics of the new technologies deployed in these vehicles operating under real world conditions, data from Smith Electric and Navistar medium-duty EVs were collected, compiled, and analyzed by the National Renewable Energy Laboratory's (NREL) Fleet Test and Evaluation team over a period of 3 years. More than 430 Smith Newton EVs have provided data representing more than 150,000 days of operation. Similarly, data have been collected from more than 100 Navistar eStar EVs, resulting in a comparative total of more than 16,000 operating days. Combined, NREL has analyzed more than 6 million kilometers of driving and 4 million hours of charging data collected from commercially operating medium-duty electric vehicles in various configurations. In this paper, extensive duty-cycle statistical analyses are performed to examine and characterize common vehicle dynamics trends and relationships based on in-use field data. The results of these analyses statistically define the vehicle dynamic and kinematic requirements for each vehicle, aiding in the selection of representative chassis dynamometer test cycles and the development of custom drive cycles that emulate daily operation. In this paper, the methodology and accompanying results of the duty-cycle statistical analysis are presented and discussed. Results are presented in both graphical and tabular formats illustrating a number of key relationships between parameters observed within the data set that relate to medium duty EVs.

33 ADVANCED PROPULSION SYSTEMS↗

Statistical Techniques For Real-time Anomaly Detection Using Spark Over Multi-source VMware Performance Data

Anomaly detection refers to the identi cation of an irregular or unusual pat- tern which deviates from what is standard, normal, or expected. Such deviated patterns typically correspond to samples of interest and are assigned different labels in different domains, such as outliers, anomalies, exceptions, or malware. Detecting anomalies in fast, voluminous streams of data is a formidable chal- lenge. This paper presents a novel, generic, real-time distributed anomaly detection framework for heterogeneous streaming data where anomalies appear as a group. We have developed a distributed statistical approach to build a model and later use it to detect anomaly. As a case study, we investigate group anomaly de- tection for a VMware-based cloud data center, which maintains a large number of virtual machines (VMs). We have built our framework using Apache Spark to get higher throughput and lower data processing time on streaming data. We have developed a window-based statistical anomaly detection technique to detect anomalies that appear sporadically. We then relaxed this constraint with higher accuracy by implementing a cluster-based technique to detect sporadic and continuous anomalies. We conclude that our cluster-based technique out- performs other statistical techniques with higher accuracy and lower processing time.

Apache Spark↗

Statistical Sampling Plan for AGR-5/6/7 Fuel Materials

Kernels, coated particles, and compacts will be manufactured to predefined specifications for irradiation in the AGR 5/6/7 experiments. Many tests will be performed for quality assurance and most of these tests correspond to criteria that have to be met with 95% confidence, based on random samples. This report describes the statistical acceptance testing methods and gives minimum sample sizes for kernels to achieve 95% confidence. Coated particles and compacts will be addressed in revisions to this report. The methods used for sample size determination are described in detail in INL/EXT 05 00349, Statistical Methods Handbook for Advanced Gas Reactor Fuel Materials (Einerson 2005). The statistical approach is to collect a representative sample, apply an acceptance test to the sample, and accept or reject the product based on this test for each criterion. The criteria or “product specifications” are contained in SPC-1352, “AGR 5/6/7 Fuel Specification,” Revision 1 and are duplicated in this plan.

A Statistical Testing Approach for Quantifying Software Reliability; Application to an Example System

The U.S. Nuclear Regulatory Commission (NRC) encourages the use of probabilistic risk assessment (PRA) technology in all regulatory matters, to the extent supported by the state-of-the-art in PRA methods and data. Although much has been accomplished in the area of risk-informed regulation, risk assessment for digital systems has not been fully developed. The NRC established a plan for research on digital systems to identify and develop methods, analytical tools, and regulatory guidance for (1) including models of digital systems in the PRA's of nuclear power plants (NPPs), and, (2) incorporating digital systems in the NRC's risk-informed licensing and oversight activities. Under NRC's sponsorship, Brookhaven National Laboratory (BNL) explored approaches for addressing the failures of digital instrumentation and control (I and C) systems in the current NPP PRA framework. Specific areas investigated included PRA modeling digital hardware , development of a philosophical basis for defining software failure , and identification of desirable attributes of quantitative software reliability methods . Based on the earlier research, statistical testing is considered a promising method for quantifying software reliability. It is widely recognized that software failures are due to the triggering of pre-existing defects by the software's operational environment. Software defects can arise from errors made in user requirements or coding errors introduced during the developmental process. The software's operational environment includes factors such as the time history of digital system inputs, communication interfaces, the internal state of the digital system, and external conditions. Thus, software reliability is a function of both the number of pre-existing defects and the presence of a triggering condition caused by the manner in which the software is used. In this paper, we describe a statistical software-testing approach for quantifying software reliability and applied it to the loop-operating control system (LOCS) of an experimental loop of the Advanced Test Reactor (ATR) at Idaho National Laboratory (INL). The work involved collaboration between BNL and INL. The objectives of the study include: (1) Development of a statistical testing approach for estimating software failure probability on demand, the results of which are suitable for including in a probabilistic risk assessment (PRA); and, (2) Application of the approach to the LOCS to estimate its failure probability, and obtain insights into the feasibility, practicality, and usefulness of the estimation in models of digital systems for inclusion in nuclear power plants' PRAs. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Effects of statistical uncertainty in activation calculations based on Monte Carlo radiation transport - 25353

This paper focuses on the analysis of the uncertainty in isotopic concentrations determined by Monte Carlo depletion calculations performed for the purposes of activation source calculations. In this work, we use repeated, statistically independent, Monte Carlo depletion calculations to examine the convergence behavior and statistical uncertainty in isotopic densities. Performing many, small, Monte Carlo depletion simulations and then recombining them is a useful approach and enables statistical analysis that can help ensure convergence. Convergence behavior is discussed in terms of stability of the mean, trends in the 95% confidence interval, and distributions of the isotopic densities. Just as in the more familiar Monte Carlo simulations of fluxes or reaction rates, convergence of isotopic densities may exhibit expected behavior such as stable means and 95% confidence intervals that decrease like N{sup -1/2} , or may exhibit poorly-converging behaviors. The most poorly-converging cases in this work were found to be for the products of high-energy threshold reactions. Finally, we examine the conditions under which uncertainty in the isotopic density can be determined from reaction rate uncertainty, and demonstrate the need for considering correlated variates in isotopic uncertainty propagation. Correlation between reaction rate tallies appears most important when the cross sections for the reactions and neutron energies considered are not strongly varying. (authors)

61 RADIATION PROTECTION AND DOSIMETRY↗

Dark Energy Survey Year 3 results: optimized $w$CDM simulation-based inference with weak lensing map-level hybrid statistics

We present cosmological constraints from the Dark Energy Survey Year 3 (DES Y3) weak lensing data using hierarchical hybrid statistics within a Bayesian simulation-based inference framework that is based on the Gower Street simulations. To maximize the precision of the inference, we have developed a new, information-theory based, data compression of the weak lensing maps to just seven highly informative summary statistics. The hybrid scheme exploits the high information content of the power spectrum, compressing both the power spectrum and neural-based summaries that are designed to extract further information. Our simulation-based approach enables principled forward modelling of all major sources of systematic uncertainty and survey properties into realistic mock observations, including the survey mask, photometric redshift uncertainties, intrinsic galaxy alignments, multiplicative shear calibration bias, source galaxy clustering, non-Gaussian shape noise, and non-linear structure formation. The summary statistics are then used in a Bayesian simulation-based inference pipeline. The inference is validated through coverage tests and checks for robustness against baryonic feedback. Assuming a $w$CDM cosmology, our analysis yields $S_8 = 0.808 \pm 0.017$, $Ω_{\rm m} = 0.325 \pm 0.024$, and $w < -0.766$ (marginalized posterior 68 per cent credible intervals). This rigorous combination of information theory, physics- and neural network-based extreme data compression, and principled Bayesian analysis improves the figure of merit for $(Ω_{\rm m}, S_8, w)$ by 60 per cent over the previous state-of-the-art, and by almost a factor of 3 over two-point analyses of the same data. They are the most precise joint constraints on $(Ω_{\rm m}, S_8, w)$ from weak gravitational lensing data alone of any survey to date. We intend to apply this analysis to the more recent DES Y6 data.

Williamson, J. [University Coll. London]↗

Statistics of Point Vortex Turbulence in Non-neutral Flows and in Flows with Translational and Rotational Symmetries

A theory (Esler and Ashbee in J Fluid Mech 779:275–308, 2015) describing the statistics of N freely-evolving point vortices in a bounded two-dimensional domain is extended. First, the case of a non-neutral vortex gas is addressed, and it is shown that the density of states function can be identified with the probability density function of an infinite sum of independent non-central chi-squared random variables, the details of which depend only on the shape of the domain. Equations for the equilibrium energy spectrum and other statistical quantities follow, the validity of which are verified against direct numerical simulations of the equations of motion. Second, domains with additional conserved quantities associated with a symmetry (e.g., circle, periodic channel) are investigated, and it is shown that the treatment of the non-neutral case can be modified to account for the additional constraint.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

An efficient direct solver for rarefied gas flows with arbitrary statistics

A new numerical methodology associated with a unified treatment is presented to solve the Boltzmann–BGK equation of gas dynamics for the classical and quantum gases described by the Bose–Einstein and Fermi–Dirac statistics. Utilizing a class of globally-stiffly-accurate implicit–explicit Runge–Kutta scheme for the temporal evolution, associated with the discrete ordinate method for the quadratures in the momentum space and the weighted essentially non-oscillatory method for the spatial discretization, the proposed scheme is asymptotic-preserving and imposes no non-linear solver or requires the knowledge of fugacity and temperature to capture the flow structures in the hydrodynamic (Euler) limit. The proposed treatment overcomes the limitations found in the work by Yang and Muljadi (2011) [33] due to the non-linear nature of quantum relations, and can be applied in studying the dynamics of a gas with internal degrees of freedom with correct values of the ratio of specific heat for the flow regimes for all Knudsen numbers and energy wave lengths. The present methodology is numerically validated with the unified treatment by the one-dimensional shock tube problem and the two-dimensional Riemann problems for gases of arbitrary statistics. Descriptions of ideal quantum gases including rotational degrees of freedom have been successfully achieved under the proposed methodology.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Computational and analytical comparison of flux discretizations for the semiconductor device equations beyond Boltzmann statistics

We compare three thermodynamically consistent numerical fluxes known in the literature, appearing in a Voronoï finite volume discretization of the van Roosbroeck system with general charge carrier statistics. Our discussion includes an extension of the Scharfetter–Gummel scheme to non-Boltzmann (e.g. Fermi–Dirac) statistics. It is based on the analytical solution of a two-point boundary value problem obtained by projecting the continuous differential equation onto the interval between neighboring collocation points. Hence, it serves as a reference flux. The exact solution of the boundary value problem can be approximated by computationally cheaper fluxes which modify certain physical quantities. One alternative scheme averages the nonlinear diffusion (caused by the non-Boltzmann nature of the problem), another one modifies the effective density of states. To study the differences between these three schemes, we analyze the Taylor expansions, derive an error estimate, visualize the flux error and show how the schemes perform for a carefully designed p-i-n benchmark simulation. We present strong evidence that the flux discretization based on averaging the nonlinear diffusion has an edge over the scheme based on modifying the effective density of states.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

On the statistical theory of self-gravitating collisionless dark matter flow: Scale and redshift variation of velocity and density distributions

The statistics of velocity and density fields are crucial for cosmic structure formation and evolution. Here, this paper extends our previous work on the two-point second-order statistics for the velocity field [Phys. Fluids 35, 077105 (2023)] to one-point probability distributions for both density and velocity fields. The scale and redshift variation of density and velocity distributions are studied by a halo-based non-projection approach. First, all particles are divided into halo and out-of-halo particles so that the redshift variation can be studied via generalized kurtosis of distributions for halo and out-of-halo particles, respectively. Second, without projecting particle fields onto a structured grid, the scale variation is analyzed by identifying all particle pairs on different scales $r$. We demonstrate that: (i) Delaunay tessellation can be used to reconstruct the density field. The density correlation, spectrum, and dispersion functions were obtained, modeled, and compared with the N-body simulation; (ii) the velocity distributions are symmetric on both small and large scales and are non-symmetric with a negative skewness on intermediate scales due to the inverse energy cascade on small scales with a constant rate $\varepsilon_u$; (iii) On small scales, the even order moments of pairwise velocity $\Delta u_L$ follow a two-thirds law $\propto{(-\varepsilon_ur)}^{2/3}$, while the odd order moments follow a linear scaling $\langle(\Delta u_L)^{2n+1}\rangle=(2n+1)\langle(\Delta u_L)^{2n}\rangle\langle\Delta u_L\rangle\propto{r}$; (iv) The scale variation of the velocity distributions was studied for longitudinal velocities $u_L$ or $u_L^{'}$, pairwise velocity (velocity difference) $\Delta u_L$=$u_L^{'}$-$u_L$ and velocity sum $\Sigma u_L$=$u^{'}_L$+$u_L$. Fully developed velocity fields are never Gaussian on any scale, despite that they can initially be Gaussian; (v) On small scales, $u_L$ and $\Sigma u_L$ can be modeled by a $X$ distribution to maximize the entropy of the system. The distribution of $\Delta u_L$ can be different; (vi) On large scales, $\Delta u_L$ and $\Sigma u_L$ can be modeled by a logistic or a $X$ distribution, while $u_L$ has a different distribution; (vii) the redshift variation of the velocity distributions follows the evolution of the $X$ distribution involving a shape parameter $\alpha(z)$ decreasing with time.

79 ASTRONOMY AND ASTROPHYSICS↗

A NEW STATISTICAL PERSPECTIVE TO THE COSMIC VOID DISTRIBUTION

In this study, we obtain the size distribution of voids as a three-parameter redshift-independent log-normal void probability function (VPF) directly from the Cosmic Void Catalog (CVC). Although many statistical models of void distributions are based on the counts in randomly placed cells, the log-normal VPF that we obtain here is independent of the shape of the voids due to the parameter-free void finder of the CVC. We use three void populations drawn from the CVC generated by the Halo Occupation Distribution (HOD) Mocks, which are tuned to three mock SDSS samples to investigate the void distribution statistically and to investigate the effects of the environments on the size distribution. As a result, it is shown that void size distributions obtained from the HOD Mock samples are satisfied by the three-parameter log-normal distribution. In addition, we find that there may be a relation between the hierarchical formation, skewness, and kurtosis of the log-normal distribution for each catalog. We also show that the shape of the three-parameter distribution from the samples is strikingly similar to the galaxy log-normal mass distribution obtained from numerical studies. This similarity between void size and galaxy mass distributions may possibly indicate evidence of nonlinear mechanisms affecting both voids and galaxies, such as large-scale accretion and tidal effects. Considering the fact that in this study, all voids are generated by galaxy mocks and show hierarchical structures in different levels, it may be possible that the same nonlinear mechanisms of mass distribution affect the void size distribution.

79 ASTRONOMY AND ASTROPHYSICS↗

The Global Statistical Response of the Outer Radiation Belt During Geomagnetic Storms

Using the total radiation belt electron content calculated from Van Allen Probe phase space density, the time-dependent and global response of the outer radiation belt during storms is statistically studied. Using phase space density reduces the impacts of adiabatic changes in the main phase, allowing a separation of adiabatic and nonadiabatic effects and revealing a clear modality and repeatable sequence of events in storm time radiation belt electron dynamics. This sequence exhibits an important first adiabatic invariant (μ)-dependent behavior in the seed (150 MeV/G), relativistic (1,000 MeV/G), and ultrarelativistic (4,000 MeV/G) populations. The outer radiation belt statistically shows an initial phase dominated by loss followed by a second phase of rapid acceleration, while the seed population shows little loss and immediate enhancement. Lastly, the time sequence of the transition to the acceleration is also strongly μ dependent and occurs at low μ first, appearing to be repeatable from storm to storm.

58 GEOSCIENCES↗

A hierarchical analysis of the impact of methodological decisions on statistical downscaling of daily precipitation and air temperatures

Despite the widespread application of statistical downscaling tools, uncertainty remains regarding the role of model formulation in determining model skill for daily maximum and minimum temperature (Tmax and Tmin), and precipitation occurrence and intensity. Impacts of several key aspects of statistical transfer function form on model skill are evaluated using a framework resistant to model overspecification. We focus on: (a) model structure: simple (generalized linear models, GLMs) versus complex (artificial neural networks, ANNs) models. (b) Predictor selection: Fixed number of predictors chosen a priori versus stepwise selection of predictors and inclusion of grid point values versus predictors derived from application of principal components analysis (PCA) to spatial fields. We also examine the influence of domain size on model performance. For precipitation downscaling, we consider the role of the threshold used to characterize a wet day and apply three approaches (Poisson and Gamma distributions in GLM and ANN) to downscale wet-day precipitation amounts. While no downscaling formulation is optimal for all predictands and at 10 locations representing diverse U.S. climates, and due to the exclusion of variance inflation all of the downscaling formulations fail to reproduce the range of observed variability, models with larger suites of prospective predictors generally have higher skill. For temperature downscaling, ANNs generally outperform GLM, with greater improvements for T min than T max . Furthermore, use of PCA-derived predictors does not systematically improve model skill, but does improve skill for temperature extremes. Model skill for precipitation occurrence generally increases as the wet-day threshold increases and models using PCA-derived predictors tend to outperform those based on grid cell predictors. Each model for wet-day precipitation intensity overestimates annual total precipitation and underestimates the proportion derived from extreme precipitation events, but ANN-based models and those with larger predictor suites tend to have the smallest bias.

54 ENVIRONMENTAL SCIENCES↗

Highly accelerated life testing (HALT): A review from a statistical perspective

Despite its use in one form or another for at least four decades, HALT and related techniques [e.g., highly accelerated-stress screening (HASS) and stress audits (HASA)] are not well understood within the statistical community and remain controversial. This largely reflects a conflict in motivation between engineers, testing under harsh conditions to discover and eliminate failure modes, and statisticians, taking a more cautious approach to develop quantitative estimates of parameters such as mean time between failures (MTBF). Here, this review article will clarify HALT concepts and methods and explain where it fits within the universe of methods that involve the application of accelerating factors to compress the time required to evaluate or enhance product reliability. A major distinction is between methods such as HALT, a high-stress test-analyze-fix-test iterative process directed at improving reliability by discovering and fixing weak points in a design, and quantitative accelerated life testing (QALT), whose goal is the estimation of product life for a fixed design. We discuss methods such as physics of failure that offer some hope of bridging the gap between the qualitative nature of HALT, and purely quantitative statistical methods. We present a variety of engineering applications of HALT including metal fatigue, piping and pressure vessels, structural damage, radiation damage, and rotating machinery. We also discuss potential synergies between HALT and QALT, such as rapid identification, through HALT, of failure modes requiring quantitative analysis. For further study, extensive references to the applicable literature are provided as well as an appendix that describes related methods.

97 MATHEMATICS AND COMPUTING↗