Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “sampling 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 1,261 records · Page 70

Statistical Treatise on Critical Biodiesel (B100) Quality Properties in the United States from 2004-2022

The quality of neat biodiesel (B100) is critical for ensuring biodiesel blends used in diesel-powered vehicles do not adversely impact engine performance. In the United States, B100 is required to meet ASTM International's purity and fuel property requirements in D6751, "Standard Specification for Biodiesel Fuel Blend Stock (B100) for Middle Distillate Fuels." Here we review the development of this standard for the different grades of B100. The BQ-9000 program, which currently covers over 90% of U.S. and Canadian production volumes, is also described. Engine and original equipment manufacturers have expressed a desire for credible, third-party data on values of various ASTM B100 properties in the commercial market to inform their efforts to address future emissions and durability requirements. To address this need, the National Renewable Energy Laboratory published data from analysis and testing of monthly samples from BQ-9000 producers from calendar years 2017-2022, resulting in 300-500 data points (i.e., individual fuel samples) per year. These data and a statistical evaluation of properties including the average, mean, standard deviation, and 95th or 5th percentile are presented. The results show very consistent and high-quality B100 produced in this time frame. Additionally, we review the results of earlier quality surveys conducted in the United States between 2004 and 2011, revealing how quality has changed over time in response to ASTM standard revisions, as well as the significant overall improvement in quality for modern biodiesel production.

ASTM standard↗

A COMPREHENSIVE STATISTICAL DESCRIPTION OF RADIO-THROUGH- γ -RAY SPECTRAL ENERGY DISTRIBUTIONS OF ALL KNOWN BLAZARS

We combined multi-wavelength data for blazars from the Roma-BZCAT catalog and analyzed hundreds of X-ray spectra. We present the fluxes and spectral energy distributions (SEDs), in 12 frequency bands from radio to γ -rays, for a final sample of 2214 blazars. Using a model-independent statistical approach, we looked for systematic trends in the SEDs; the most significant trends involved the radio luminosities and X-ray spectral indices of the blazars. We used a principal component analysis (PCA) to determine the basis vectors of the blazar SEDs and, in order to maximize the size of the sample, imputed missing fluxes using the K -nearest neighbors method. Using more than an order of magnitude more data than was available when Fossati et al. first reported trends of SED shape with blazar luminosity, we confirmed the anti-correlation between radio luminosity and synchrotron peak frequency, although with greater scatter than was seen in the smaller sample. The same trend can be seen between bolometric luminosity and synchrotron peak frequency. Finally, we used all of the available blazar data to determine an empirical SED description that depends only on the radio luminosity at 1.4 GHz and the redshift. We verified that this statistically significant relation was not a result of the luminosity–luminosity correlations that are natural in flux-limited samples (i.e., where the correlation is actually caused by the redshift rather than the luminosity).

79 ASTRONOMY AND ASTROPHYSICS↗

Perturbation theory models for LSST-era galaxy clustering: Tests with subpercent mock catalog measurements in Fourier and configuration space

Here, we analyze the clustering of galaxies using the z = 1.006 snapshot of the CosmoDC2 simulation, a high-fidelity synthetic galaxy catalog designed to validate analysis methods for the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST). We present sub-percent measurements of the galaxy auto-correlation and galaxy-dark matter cross correlations in Fourier space and configuration space for a magnitude-limited galaxy sample. At these levels of precision, the statistical errors of the measurement are comparable to the systematic effects present in the simulation and measurement procedure; nevertheless, using a hybrid-PT model, we are able to model non-linear galaxy bias with 0.5% precision up to scales of k max = 0.5 h/Mpc and r min = 4 Mpc/h. While the linear bias parameter is measured with 0.01% precision, other bias parameters are determined with considerably weaker constraints and sometimes bimodal posterior distributions. We compare our fiducial model with lower dimensional models where the higher-order bias parameters are fixed at their co-evolution values and find that leaving these parameters free provides significant improvements in our ability to model small scale information. We also compare bias parameters for galaxy samples defined using different magnitude bands and find agreement between samples containing equal numbers of galaxies. Finally, we compare bias parameters between Fourier space and configuration space and find moderate to significant tension between the two approaches. Although our model is often unable to fit the CosmoDC2 galaxy samples within the 0.1% precision of our measurements, our results suggest that the hybrid-PT model used in this analysis is capable of modeling non-linear galaxy bias within the percent level precision needed for upcoming galaxy surveys.

79 ASTRONOMY AND ASTROPHYSICS↗

Statistical methods for efficient design of community surveys of response to noise: Random coefficients regression models

Research studies of residents' responses to noise consist of interviews with samples of individuals who are drawn from a number of different compact study areas. The statistical techniques developed provide a basis for those sample design decisions. These techniques are suitable for a wide range of sample survey applications. A sample may consist of a random sample of residents selected from a sample of compact study areas, or in a more complex design, of a sample of residents selected from a sample of larger areas (e.g., cities). The techniques may be applied to estimates of the effects on annoyance of noise level, numbers of noise events, the time-of-day of the events, ambient noise levels, or other factors. Methods are provided for determining, in advance, how accurately these effects can be estimated for different sample sizes and study designs. Using a simple cost function, they also provide for optimum allocation of the sample across the stages of the design for estimating these effects. These techniques are developed via a regression model in which the regression coefficients are assumed to be random, with components of variance associated with the various stages of a multi-stage sample design.

Tomberlin, T. J.↗

Small-Sample Estimator Decisions – Certainly Uncertain in Human Spaceflight at NASA

Within statistics when estimating means, we rely heavily on the Central Limit Theorem (CLT) to aid in our inference, however in cases where the CLT does not apply we begin to be strongly limited in our choices. If a violation to the CLT comes in the form of small samples and unknown population variance, yet we maintain the assumption of normality, the use of t-distributions is perfectly valid. We explore cases where population normality is not assumed, where sample sizes are small, and where some actionable estimate of the central tendency of a distribution is needed such as for spaceflight operational decisions. Here we identify different scenarios based on measures taken from spaceflight, characterizing performance of estimators of central tendency within our set of simulations. We examine bias and variability of standard estimators of central tendency as they apply to varying small sample sizes and varying population distributions. Further work in this area is needed to develop a framework or set of guidelines for individuals set in these situations. Even with limited information, decisions need to be made.

Central Tendency↗

Investigation of oscillations above the acoustic cut-off frequency for a sample of Kepler stars. II. Identifying temporal frequency shifts

ABSTRACT Stellar oscillations with frequencies greater than the acoustic cut-off frequency are not trapped within the star’s interior. Geometric interference between these high-frequency waves produces a peak-like structure in the power spectrum of the star, known as the pseudo-modes. In this study, we aim to analyse 51 stars previously identified (from a sample of 91 Kepler stars) to contain statistically significant evidence for pseudo-modes, to determine if the pseudo-mode frequencies vary in time. We obtained these frequency variations and the uncertainty by the resampled periodogram approach, where periodogram realisations were created from successive, overlapping 90-d time segments, and a cross-correlation function was generated between them. The computed shifts were compared against temporal p-mode frequency shifts and a stellar magnetic activity proxy, $S_{\mathrm{ph}}$. We found that for 16 per cent of the stars pseudo-mode frequency shifts were significantly anticorrelated with p-mode shifts, as is the case for the Sun. However, we also found pseudo-mode and p-mode shifts to be significantly correlated in-phase in 8 per cent of our sample. The remainder showed no significant correlation. We also searched for trends between the maximum variation of pseudo-mode frequency shifts and stellar parameters. Whilst we found no correlation between the maximum variation and stellar photospheric and chromospheric proxies, we found frequency variations to have a greater amplitude for cooler and slower rotating stars, which is in opposition to p-mode frequency shift variation. Temporal changes to pseudo-mode frequencies may be attributable to variations in stellar magnetic activity, and therefore can contribute to our understanding of stellar activity cycles.

Millson, L. J. (ORCID:0009000342542676)↗

The use of latin hypercube sampling for the efficient estimation of confidence intervals

Latin hypercube sampling (LHS) has long been used as a way of assuring adequate sampling of the tails of distributions in a Monte Carlo analysis and provided the framework for the uncertainty analysis performed in the NUREG-1150 risk assessment. However, this technique has not often been used in the performance of regulatory analyses due to the inability to establish confidence levels on the quantiles of the output distribution. Recent work has demonstrated a method that makes this possible. This method is compared to the procedure of crude Monte Carlo using order statistics, which is currently used to establish confidence levels. The results of several statistical examples demonstrate that the LHS confidence interval method can provide a more accurate and precise solution, but issues remain when applying the technique generally. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Estimating Highway Volumes Using Vehicle Probe Data - Proof of Concept: Preprint

This paper examines the feasibility of using sampled commercial probe data in combination with validated continuous counter data to accurately estimate vehicle volume across the entire roadway network, for any hour during the year. Currently either real time or archived volume data for roadways at specific times are extremely sparse. Most volume data are average annual daily traffic (AADT) measures derived from the Highway Performance Monitoring System (HPMS). Although methods to factor the AADT to hourly averages for typical day of week exist, actual volume data is limited to a sparse collection of locations in which volumes are continuously recorded. This paper explores the use of commercial probe data to generate accurate volume measures that span the highway network providing ubiquitous coverage in space, and specific point-in-time measures for a specific date and time. The paper examines the need for the data, fundamental accuracy limitations based on a basic statistical model that take into account the sampling nature of probe data, and early results from a proof of concept exercise revealing the potential of probe type data calibrated with public continuous count data to meet end user expectations in terms of accuracy of volume estimates.

25 ENERGY STORAGE↗

Mid-IR Properties of an Unbiased AGN Sample of the Local Universe: Emission-Line Diagnostics - 1

\Ve compare mid-IR emission-lines properties, from high-resolution Spitzer IRS spectra of a statistically-complete hard X-ray (14-195 keV) selected sample of nearby (z < 0.05) AGN detected by the Burst Alert Telescope (BAT) aboard Swift. The luminosity distribution for the mid-infrared emission-lines, [O IV] 25.89 microns, [Ne II] 12.81 microns, [Ne III] 15.56 microns and [Ne V] 14.32 microns, and hard X-ray continuum show no differences between Seyfert 1 and Seyfert 2 populations, although six newly discovered BAT AGNs are shown to be under-luminous in [O IV], most likely the result of dust extinction in the host galaxy. The overall tightness of the mid-infrared correlations and BAT luminosities suggests that the emission lines primarily arise in gas ionized by the AGN. We also compared the mid-IR emission-lines in the BAT AGNs with those from published studies of star-forming galaxies and LINERs. We found that the BAT AGN fall into a distinctive region when comparing the [Ne III]/[Ne II] and the [O IV]/[Ne III] quantities. From this we found that sources that have been previously classified in the mid-infrared/optical as AGN have smaller emission line ratios than those found for the BAT AGNs, suggesting that, in our X-ray selected sample, the AGN represents the main contribution to the observed line emission. Overall, we present a different set of emission line diagnostics to distinguish between AGN and star forming galaxies that can be used as a tool to find new AGN.

Weaver, K. A.↗

Far-infrared properties of Markarian galaxies with multiple nuclei - Warm dust emission in mergers

An investigation of coadded IRAS data is performed on 187 Markarian galaxies where distinguishing morphological characteristics or multiple optical nuclei are present. The far-IR properties of Markarian galaxies are compared to the IRAS Bright Galaxy Sample, and a much higher median dust temperature is found in the multiple nucleus galaxies, suggesting that more far-IR luminosity results from active star formation. Both optical/UV and far-IR selection techniques are necessary to extract complete samples of AGNs since the far-IR two-color plane can miss up to 50 percent of the galaxies. A systematic increase in the contribution of warm dust emission due to active star formation and AGNs is found in a statistical comparison of merger candidates and other galaxy samples. The assumed nature of precursor galaxies determines the assumed enhancement of far-IR luminosity caused by galaxy collisions. A model is presented which describes the properties of the Markarian galaxies in terms of enhanced OB star formation and different grain size distributions. The results of the investigation are shown to be consistent with a 'subdued' interpretation of merging galaxies with high luminosities.

Mazzarella, Joseph M.↗

DETERMINING THE LUMINOSITY FUNCTION OF SWIFT LONG GAMMA-RAY BURSTS WITH PSEUDO-REDSHIFTS

The determination of the luminosity function (LF) of gamma-ray bursts (GRBs) is an important role for the cosmological applications of the GRBs, which, however, is seriously hindered by some selection effects due to redshift measurements. In order to avoid these selection effects, we suggest calculating pseudo-redshifts for Swift GRBs according to the empirical L-E{sub p} relationship. Here, such a L-E{sub p} relationship is determined by reconciling the distributions of pseudo- and real redshifts of redshift-known GRBs. The values of E{sub p} taken from Butler's GRB catalog are estimated with Bayesian statistics rather than observed. Using the GRB sample with pseudo-redshifts of a relatively large number, we fit the redshift-resolved luminosity distributions of the GRBs with a broken-power-law LF. The fitting results suggest that the LF could evolve with redshift by a redshift-dependent break luminosity, e.g., L{sub b} = 1.2 Multiplication-Sign 10{sup 51}(1 + z){sup 2} erg s{sup -1}. The low- and high-luminosity indices are constrained to 0.8 and 2.0, respectively. It is found that the proportional coefficient between the GRB event rate and the star formation rate should correspondingly decrease with increasing redshifts.

79 ASTRONOMY AND ASTROPHYSICS↗

Designing for Small Volume Assembly of Advanced Electronics Packages

We describe a general methodology to Design for Producibility and Reliability (DFPAR) for very small volume production runs. In cases where the entire volume for fabrication is less than five products, traditional Statistical Process Control (SPC) is inadequate due to reliance on statistics of much larger volumes and the Central Limit Theorem. Data acquisition for process parameter estimation from such a small sample size is difficult; however, it is critical to producing high reliability product.

Central Limit Theorem statistics high reliability ↗

Comparison of experimental methods for estimating matrix diffusion coefficients for contaminant transport modeling

Here, diffusion cell and diffusion wafer experiments were conducted to compare methods for estimating effective matrix diffusion coefficients in rock core samples from Pahute Mesa at the Nevada Nuclear Security Site (NNSS). A diffusion wafer method, in which a solute diffuses out of a rock matrix that is pre-saturated with water containing the solute, is presented as a simpler alternative to the traditional through-diffusion (diffusion cell) method. Both methods yielded estimates of effective matrix diffusion coefficients that were within the range of values previously reported for NNSS volcanic rocks. The difference between the estimates of the two methods ranged from 14 to 30%, and there was no systematic high or low bias of one method relative to the other. From a transport modeling perspective, these differences are relatively minor when one considers that other variables (e.g., fracture apertures, fracture spacings) influence matrix diffusion to a greater degree and tend to have greater uncertainty than effective matrix diffusion coefficients. For the same relative random errors in concentration measurements, the diffusion cell method yields effective matrix diffusion coefficient estimates that have less uncertainty than the wafer method. However, the wafer method is easier and less costly to implement and yields estimates more quickly, thus allowing a greater number of samples to be analyzed for the same cost and time. Given the relatively good agreement between the methods, and the lack of any apparent bias between the methods, the diffusion wafer method appears to offer advantages over the diffusion cell method if better statistical representation of a given set of rock samples is desired.

54 ENVIRONMENTAL SCIENCES↗

Comparison of Experimental Methods for Estimating Matrix Diffusion Coefficients for Contaminant Transport Modeling

Diffusion cell and diffusion wafer experiments were conducted to compare methods for estimating matrix diffusion coefficients in rock core samples from Pahute Mesa at the Nevada Nuclear Security Site (NNSS). A diffusion wafer method, in which a solute diffuses out of a rock matrix that is pre-saturated with water containing the solute, is presented as a simpler alternative to the traditional through-diffusion (diffusion cell) method. Both methods yielded estimates of matrix diffusion coefficients that were within the range of values previously reported for NNSS volcanic rocks. The difference between the estimates of the two methods ranged from 14 to 30%, and there was no systematic high or low bias of one method relative to the other. From a transport modeling perspective, these differences are relatively minor when one considers that other variables (e.g., fracture apertures, fracture spacings) influence matrix diffusion to a greater degree and tend to have greater uncertainty than diffusion coefficients. For the same relative random errors in concentration measurements, the diffusion cell method yields diffusion coefficient estimates that have less uncertainty than the wafer method. However, the wafer method is easier and less costly to implement and yields estimates more quickly, thus allowing a greater number of samples to be analyzed for the same cost and time. Given the relatively good agreement between the methods, and the lack of any apparent bias between the methods, the diffusion wafer method appears to offer advantages over the diffusion cell method if better statistical representation of a given set of rock samples is desired.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

Identifying structural changes with unsupervised machine learning methods

In this work, unsupervised machine learning methods are used to identify structural changes using the melting point transition in classical molecular dynamics simulations as an example application of the approach. Dimensionality reduction and clustering methods are applied to instantaneous radial distributions of atomic configurations from classical molecular dynamics simulations of metallic systems over a large temperature range. Principal component analysis is used to dramatically reduce the dimensionality of the feature space across the samples using an orthogonal linear transformation that preserves the statistical variance of the data under the condition that the new feature space is linearly independent. From there, k-means clustering is used to partition the samples into solid and liquid phases through a criterion motivated by the geometry of the reduced feature space of the samples, allowing for an estimation of the melting point transition. This pattern criterion is conceptually similar to how humans interpret the data but with far greater throughput, as the shapes of the radial distributions are different for each phase and easily distinguishable by humans. The transition temperature estimates derived from this machine learning approach produce comparable results to other methods on similarly small system sizes. These results show that machine learning approaches can be applied to structural changes in physical systems.

42 ENGINEERING↗

Conservative Estimation of Tail Probabilities from Limited Sample Data

Several sparse-sample uncertainty quantification (UQ) methods are compared for conservative but not overly conservative estimation of small tail probabilities involving responses that lay beyond specified thresholds in the tails of probability distributions. Sixteen very differently shaped distributions (or probability density functions, PDFs) and tail probability magnitudes ranging from 10 -5 to 10 -1 are considered in order for the study to be relevant to a wide range of risk analysis and quantification of margins and uncertainty (QMU) problems. The emphasis of the study is on limited data regimes ranging from N = 2 to 20 samples, reflective of most experimental and some expensive computational situations. Relatively simple sparse-sample UQ methods tested for this regime involve statistical tolerance interval "Equivalent Normal and related "Ensemble of Normals" and "Superdistribution (SD) approaches. (The independently derived SD is effectively equivalent to the Bayesian posterior predictive distribution given the assumptions of the derivation.) The performance of the methods was generally improved for N ≥ 5 samples with a generalized Jackknife resampling technique, which determines a tail probability estimate by averaging estimates from smaller sub-samples. Several quantitative metrics for method conservatism and accuracy of tail probability estimation are used to assess and rank the methods' performance over many random trials for each test PDF and probability magnitude. A variant of Bootstrap resampling was also tried, but did not significantly improve tail probability estimates in most cases. Detailed results are presented from over 100-million tests over the above factors that provide useful granular information on which methods or combination of methods perform best in various areas of the factor space.

97 MATHEMATICS AND COMPUTING↗

Shuttle Hypervelocity Impact Database

With three missions outstanding, the Shuttle Hypervelocity Impact Database has nearly 3000 entries. The data is divided into tables for crew module windows, payload bay door radiators and thermal protection system regions, with window impacts compromising just over half the records. In general, the database provides dimensions of hypervelocity impact damage, a component level location (i.e., window number or radiator panel number) and the orbiter mission when the impact occurred. Additional detail on the type of particle that produced the damage site is provided when sampling data and definitive analysis results are available. Details and insights on the contents of the database including examples of descriptive statistics will be provided. Post flight impact damage inspection and sampling techniques that were employed during the different observation campaigns will also be discussed. Potential enhancements to the database structure and availability of the data for other researchers will be addressed in the Future Work section. A related database of returned surfaces from the International Space Station will also be introduced.

Hyde, James L.↗

Using elliptical galaxy kinematics to compare the strength of gravity in cosmological regions of differing gravitational potential – a first look

ABSTRACT Various models of modified gravity invoke ‘screening’ mechanisms that are sensitive to the value of the local gravitational potential. This could have observable consequences for galaxies. These consequences might be seen by comparing two proxies for galaxy mass – their luminosity and their internal kinematics – as a function of local galaxy density. Motivated by this prospect, we have compared the observed properties of luminous red galaxies (LRGs) inside and outside of voids in the cosmic large scale structure. We used archival measurements of line widths, luminosities, redshifts, colours, and positions of galaxies in conjunction with recent void catalogues to construct comparison LRG samples inside and outside of voids. We fitted these two samples to the well-established fundamental plane of elliptical galaxies to constrain any differences between the inferred value of the Newtonian gravitational constant G for the two samples. We obtained a null result, with an upper limit on any fractional difference in G within and outside of cosmological voids to be α = δG/G ∼ 40 per cent. This upper bound is dominated by the small-number statistics of our N ∼ 100 within-void LRG sample. With the caveat that environmental effects could influence various parameters such as star formation, we estimate that a 1 per cent statistical limit on α could be attained with data from 105 elliptical galaxies within voids. This is within the reach of future photometric and spectroscopic surveys, both of which are required to pursue this method.

79 ASTRONOMY AND ASTROPHYSICS↗