Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “statistical”

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 1,027 records · Page 57

Determination of a most representative cycle from cylinder pressure ensembles via statistical method using distribution skewness

In internal combustion engine research, cylinder pressure measurements provide valuable information about the underlying thermodynamic and combustion processes, and are typically collected in ensembles of several 100 traces. Although in some particular fields of combustion research all traces are analyzed, in most cases only one trace is studied because analyzing all the traces is impractical due to the large number of collected samples. Instead, an ensemble-averaged pressure trace is commonly calculated and used for analysis. However, this pressure trace is highly smoothed and dynamic information is lost during the averaging process. With the average trace, pressure rise rates are lower and pressure oscillations such as the ones resulting from combustion knock are lost. In this work, a statistical method was developed to determine the “most representative cycle,” which is the cycle from the ensemble that has the pressure trace most representative of the engine operating condition. Eleven characteristic parameters are computed from each pressure trace and probabilistic distributions are obtained for each of the parameters using all the traces in the ensemble. Finally, the most representative cycle is selected by means of a cost function minimization. The benefits of this method are illustrated using experimental data from four very different engine platforms, under four different combustion modes and over a range of operating conditions.

42 ENGINEERING↗

Improving multiwell petrophysical interpretation from well logs via machine learning and statistical models

Well-log interpretation estimates in situ rock properties along well trajectory, such as porosity, water saturation, and permeability, to support reserve-volume estimation, production forecasts, and decision making in reservoir development. However, due to measurement errors, variability of well logs caused by multiple measurement vendors, different borehole tools, and nonuniform drilling/borehole conditions, estimations of rock properties with original well logs without proper preprocessing may not be accurate, especially in the context of multiwell estimation. Well-log normalization techniques such as two-point scaling and mean-variance normalization are commonly used to improve the robustness of multiwell rock-property estimation. However, these techniques do not consider the correlation between well logs and require subjective knowledge for their effective implementation. To reduce uncertainties and processing time associated with multiwell rock-property estimation from well logs, we develop discriminative adversarial (DA) and linear constraint models for well-log normalization and rock-property estimation. The DA neural network model developed for well-log normalization and interpretation can perform linear and nonlinear well-log normalization while considering the joint distribution of each well log and rock properties. However, the linear constraint model uses an ensemble of predictions from linear models to constrain well-log normalization and rock-property estimation. We also develop a divergence-based type well identification method to select type (training) wells for a test well based on the statistical similarity of associated well-log distributions instead of the interwell distance. We apply the DA model to perform well-log normalization and prediction of permeability for the Seminole San Andres Unit carbonate reservoir. Compared with the permeability predicted with the classical machine learning model without well-log normalization and models with two-point scaling normalization, the DA model yields the most accurate permeability prediction by decreasing the mean-squared error of permeability prediction by 20%–50%.

Geochemistry & Geophysics↗

A Generalizable Evaluated Approach, Applying Advanced Geospatial Statistical Methods, to Identify High Lead Exposure Locations at Census Tract Scale: Michigan Case Study

BACKGROUND: Despite great progress in reducing environmental lead (Pb) levels, many children in the United States are still being exposed. OBJECTIVE: Our aim was to develop a generalizable approach for systematically identifying, verifying, and analyzing locations with high prevalence of children’s elevated blood Pb levels (EBLLs) and to assess available Pb models/indices as surrogates, using a Michigan case study. METHODS: We obtained ~1:9 million BLL test results of children <6 years of age in Michigan from 2006–2016; we then evaluated them for data representativeness by comparing two percentage EBLL (%EBLL) rates (number of children tested with EBLL divided by both number of children tested and total population). We analyzed %EBLLs across census tracts over three time periods and between two EBLL reference values (≥5 vs. ≥10 μg/dL) to evaluate consistency. Locations with high %EBLLs were identified by a top 20 percentile method and a Getis-Ord Gi* geospatial cluster “hotspot” analysis. For the locations identified, we analyzed convergences with three available Pb exposure models/indices based on old housing and sociodemographics. RESULTS: Analyses of 2014–2016 %EBLL data identified 11 Michigan locations via cluster analysis and 80 additional locations via the top 20 percentile method and their associated census tracts. Data representativeness and consistency were supported by a 0.93 correlation coefficient between the two EBLL rates over 11 y, and a Kappa score of ~0:8 of %EBLL hotspots across the time periods (2014–2016) and reference values. Many EBLL hotspot locations converge with current Pb exposure models/indices; others diverge, suggesting additional Pb sources for targeted interventions. DISCUSSION: This analysis confirmed known Pb hotspot locations and revealed new ones at a finer geographic resolution than previously available, using advanced geospatial statistical methods and mapping/visualization. It also assessed the utility of surrogates in the absence of blood Pb data. This approach could be applied to other states to inform Pb mitigation and prevention efforts. https://doi.org/10.1289/EHP9705

54 ENVIRONMENTAL SCIENCES↗

National Summary Statistics for Depot-Based Medium- and Heavy-Duty Vehicle Operations

This file provides nationally aggregated summary statistics to characterize the daily operations of medium- and heavy-duty vehicles. It provides typical daily driving distances, domicile dwell durations, and one version of potential normalized electric vehicle depot charging load curves. The methodology to create these data was designed with national representativeness in mind, and the data are suitable for national analysis.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Manipulation and statistical analysis of the fluid flow of polymer semiconductor solutions during meniscus-guided coating

Recent work in structure–processing relationships of polymer semiconductors have demonstrated the versatility and control of thin-film microstructure offered by meniscusguided coating (MGC) techniques. Here, we analyze the qualitative and quantitative aspects of solution shearing, a model MGC method, using coating blades augmented with arrays of pillars. The pillars induce local regions of high strain rates—both shear and extensional—not otherwise possible with unmodified blades, and we use fluid mechanical simulations to model and study a variety of pillar spacings and densities. We then perform a statistical analysis of 130 simulation variables to find correlations with three dependent variables of interest: thin-film degree of crystallinity and transistor field-effect mobilities for charge-transport parallel (Μ para ) and perpendicular (Μ perp ) to the coating direction. Our study suggests that simple fluid mechanical models can reproduce substantive correlations between the induced fluid flow and important performance metrics, providing a methodology for optimizing blade design.

36 MATERIALS SCIENCE↗

A Robust Statistical Algorithm for Boundary Detection in Liquid Sprays

This work presents a statistical algorithm to detect phase boundaries in liquid sprays. Experimental data were obtained in a nonevaporating acetone spray at atmospheric conditions to minimize the refractive index gradients and associated beam steering. Mie scattering and diffuse background illumination (DBI) diagnostics are implemented to demonstrate the robustness of the approach. The analysis shows that the present algorithm yields similar results for the two different diagnostics, unlike existing threshold-based methodologies providing inconsistent results. The algorithm can be used to identify partial liquid boundaries between single-phase liquid and vapor boundaries. Thus, for each case, liquid boundaries were defined in terms of liquid probability of 100% or less than 100%. Sprays from two different injectors were analyzed to demonstrate the versatility of the methodology. Finally, a parametric study was conducted to demonstrate that motion blur can affect an accurate detection of the phase boundaries, and hence, the exposure time in the experiment must be sufficiently small.

42 ENGINEERING↗

Deterministic and Statistical Modeling of a New Thermal Breakout Technology for Measuring the Maximum Horizontal In-Situ Stress

The current state-of-the-art technology for in-situ stress measurements involves an integrated approach that combines borehole breakout observations, drilling-induced tensile fractures, and hydraulic fracturing tests (i.e., “mini-fracs”). This methodology has achieved wide application in the oil and gas industry but has several limitations that often prevent successful in-situ stress measurements. One major limitation is that breakouts do not appear in all boreholes and are generally only a natural occurrence that cannot easily be controlled. Because borehole breakouts are used to directly measure the maximum horizontal in-situ stress magnitude, the absence of borehole breakouts presents a major data gap for in-situ stress measurements. In response to this data gap, a new US Department of Energy (US DOE)-sponsored thermal breakout technology that will provide a method for thermally inducing borehole breakouts and allow the consistent measurement of the maximum horizontal stress magnitude is currently in development. This thermal breakout technology involves heating the borehole and increasing the thermoelastic compressive stress in the rock until a breakout develops, which can be directly correlated to the maximum horizontal stress magnitude. The first step in this project was an analytical modeling study of the thermal breakout process. Based on the Kirsch solution (Kirsch 1898), a deterministic and statistical analysis was performed on the pertinent parameters that influence the maximum horizontal stress calculation. As a result of the analysis, the findings indicate that the thermal breakout technology is feasible and provides improved accuracy and/or an enhanced ability to measure the maximum horizontal stress. Future work as part of this US DOE-sponsored project includes additional validation through more detailed numerical modeling, laboratory testing, and field testing of the thermal breakout technology.

58 GEOSCIENCES↗

Evaluation of Component Reliability in Photovoltaic Systems using Field Failure Statistics

Ongoing operations and maintenance (O&M) are needed to ensure photovoltaic (PV) systems continue to operate and meet production targets over the lifecycle of the system. Although average costs to operate and maintain PV systems have been decreasing over time, reported costs can vary significantly at the plant level. Estimating O&M costs accurately is important for informing financial planning and tracking activities, and subsequently lowering the levelized cost of electricity (LCOE) of PV systems. This report describes a methodology for improving O&M planning estimates by using empirically-derived failure statistics to capture component reliability in the field. The report also summarizes failure patterns observed for specific PV components and local environmental conditions observed in Sandia's PV Reliability, Operations & Maintenance (PVROM) database, a collection of field records across 800+ systems in the U.S. Where system-specific or fleet-specific data are lacking, PVROM-derived failure distribution values can be used to inform cost modeling and other reliability analyses to evaluate opportunities for performance improvements.

14 SOLAR ENERGY↗

Measurement of Spontaneous Undulator Radiation Statistics Generated by a Single Electron (URSSE)

As previously stated in the experiment proposal [1], the goal of the URSSE experiment was to study statistical properties of the undulator radiation generated by a single electron, circulating inthe IOTA storage ring. Namely, the idea was to study the distribution of the number of photons detected in a certain time window, and the distribution of the inter-arrival times of the photons. We wanted to see if the photostatistics was exactly Poissonian, as predicted or if there was some deviation towards Sub-Poissonian light, similar to what was observed in another study.

43 PARTICLE ACCELERATORS↗

The Idiot's Guide to the Statistical Theory of Fission Chains

What follows are my personal notes developed as I’ve struggled to understand the Statistical Theory of Fission Chains. The theoretical work detailed here is principally the work of R. P. Feynman, G. I. Bell, K. Boehnel, W. Hage, D. M. Cifarelli, M. Prasad, and N. Snyderman, plus a little of my own dablings around the edges.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Statistical Hauser-Feshbach model calculation for (α, n) reaction cross section and discrete level population

We calculate the α-particle induced reactions on 17,18 O, 19 F, and 23 Na, in the energy range 0 ≤ E α ≤ 10 MeV, with a particular attention to the branching ratios of populated discrete levels in the (α,n) channels. Since there are too many open channels to employ the R-matrix theory, we apply the statistical Hauser-Feshbach model to calculate these reaction cross sections and branching ratios. The branching ratio is defined as b n = $\frac{σ_n}{Σ_iσ_i + σ_c}$ , where σ n is the n-th level production cross section after the neutron emission, and σ c is the production of the continuum state. In the cases of our target nuclei and the energy range of interest, the residual nuclei of the neutron emission channel are always in their discrete states, so that σ c can be negligible.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Statistical Distributions for Mesh Independent Solutions in ALEGRA

The representation of material heterogeneity (also referred to as "spatial variation") plays a key role in the material failure simulation method used in ALEGRA. ALEGRA is an arbitrary Lagrangian-Eulerian shock and multiphysics code developed at Sandia National Laboratories and contains several methods for incorporating spatial variation into simulations. A desirable property of a spatial variation method is that it should produce consistent stochastic behavior regardless of the mesh used (a property referred to as "mesh independence"). However, mesh dependence has been reported using the Weibull distribution with ALEGRA's spatial variation method. This report describes efforts towards providing additional insight into both the theory and numerical experiments investigating such mesh dependence. In particular, we have implemented a discrete minimum order statistic model with properties that are theoretically mesh independent.

36 MATERIALS SCIENCE↗

Studying Rate of Fire Spread Using Landscape Statistics and FIRETEC

Wildfires are multi-system and multi-scale processes that are influenced by the land and atmosphere. Here, we study the influence of land, in terms of topography and vegetation, on uphill and downhill fire spread using data of the state of California as a testbed. The California landscape is first quantified by reclassifying its dominant vegetation types into grass, trees, and shrubs. Next, statistics of slope, vegetation height, and live fuel moisture percentage for each dominant vegetation type are gathered to create ensembles of idealized computational fluid dynamics (CFD) model simulations representative of the landscape. Rate of spread values for each CFD simulation are recorded to observe the differences in perturbing slope and fuel moisture for each vegetation type. We have found that all vegetation types prefer slopes of less than 30 degrees. Trees in the California data are 10 to 40 meters high, whereas shrub heights are less than 3 meters and grass heights are less than 1 meter. The live fuel moisture percentage of both trees and shrubs fluctuates around approximately 100% over one year (2019-2020). Using these parameters as inputs into the CFD model, fire spread rate is found to be fastest in the grass ensemble (maximum 12 m s -1 ), followed by shrubs (maximum 10 m s -1 ) and lastly by trees (maximum 7 m s -1 ). Increasing slope increases the spread rate of the uphill fire but the downhill spread decreases for slopes greater than 15 degrees. Increasing fuel moisture decreases the spread rate both up and down the hill. This study highlights the impacts that topography and vegetation have on fire spread.

54 ENVIRONMENTAL SCIENCES↗

Quantifying Uncertainty in All-to-All Estimates of Space Object Conjunction Probabilities using U-Statistics

Predicting space object conjunctions is inherently probabilistic due to initial state and orbit model uncertainty. A commonly considered Monte Carlo estimator of the conjunction probability is the ’all-to-all’ estimator. Given independent random samples of the trajectories of both objects, the estimator is the percentage of all pairs of trajectories that result in a conjunction. Intuitively, the all-to-all estimator is the best possible estimator of the conjunction probability since it considers all pairs of Monte Carlo samples. However, its distribution is not available in closed-form, which limits its use in practice and makes this intuition difficult to make rigorous. In this paper, the all-to-all estimator is identified as a U-statistic, which implies that it has several favorable properties. Specifically, the estimator is the minimum variance unbiased estimator of the conjunction probability and is asymptotically Gaussian distributed. An approximate confidence interval for the conjunction probability is obtained from an estimate of the asymptotic Gaussian distribution. We show how to efficiently compute the confidence interval and demonstrate that the interval has the nominal coverage level. The confidence intervals are also seen to be narrower than those based on the commonly-used each-to-each estimator. Furthermore, the all-to-all estimator is shown to allow different Monte Carlo sample sizes, whereas the each-to-each estimator requires equal sample sizes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The STAT7 Code for Statistical Propagation of Uncertainties In Steady-State Thermal Hydraulics Analysis of Plate-Fueled Reactors

STAT7 was written to automate many of the steady-state thermal hydraulic safety calculations for the MIT research reactor, both for conversion of the reactor from highly enriched uranium fuel to low-enriched uranium fuel and for future fuel re-loads after the conversion. A Monte-Carlo statistical propagation approach is used to treat uncertainties in important parameters in the analysis. These safety calculations are ultimately intended to protect against high fuel plate temperatures due to critical heat flux or departure from nucleate boiling or onset of flow instability; but additional margin is obtained by basing the limiting safety settings on avoiding onset of nucleate boiling. STAT7 can simultaneously analyze all of the axial nodes of all of the fuel plates and all of the coolant channels for one stripe of a fuel element. The stripes run the length of the fuel, from the bottom to the top. Power splits are calculated for each axial node of each plate to determine how much of the power goes out each face of the plate. By running STAT7 multiple times, full core analysis can be performed by analyzing the margin to ONB for each axial node of each stripe of each plate of each element in the core.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Synthesis of observed and simulated rain microphysics to inform a new Bayesian statistical framework for microphysical parameterization in climate models (Final Scientific Report)

This project aimed to investigate a new approach to bulk microphysics schemes. We identified the need to move beyond the fixed structural assumptions and approximation of existing schemes. For example, most schemes assume some functional form for the rain drop size distribution (e.g. a gamma or exponential distribution). Most schemes then evolve some number of statistical moments of that distribution via various processes, such as evaporation, sedimentation, collision-coalescence, and collisional breakup. These microphysical processes, in turn, are typically some fixed functional form derived from either some other (more detailed) model, or using some single-particle rates that are then integrated over the size distribution. While these process rate functions may have some free parameters to adjust, in other cases doing so is impossible (e.g. when the process rates are analytical or piecewise solutions to some target function). Our goal was to build and test a scheme that assumed no size distribution form, and used a series of power laws as the basis for the process rates. Moments of the size distribution would be predicted, but no underlying size distribution would be specified. Any number or choice of size distribution moments could be used, and any number of power laws could be employed to model the process rates. Thus, our approach is flexible, and can seamlessly scale across levels of complexity, as demanded by the data. Bayesian inference then would provide the formalism to estimate the model parameters and structure, allowing for robust uncertainty quantification. We proposed testing this framework in idealized simulations using bin microphysical schemes as a data source. We also proposed using real data to inform our microphysics scheme, and also that we would integrate our scheme into WRF.

58 GEOSCIENCES↗

Solving Statistical Flowgraph Models: A Brief Review

Statistical flowgraph models (SFGMs) have proven useful for analysis and modeling of multistate processes in various fields. In this paper we briefly review flowgraph models, describe how to set up a model, and provide a fully worked example from biostatistics. References are provided for more details.

97 MATHEMATICS AND COMPUTING↗