Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “methods: 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,009 records · Page 56

Statistical Comparison and Improvement of Methods for Combining Random and Harmonic Loads

Structures in many environments experience both random and harmonic excitation. A variety of closed-form techniques has been used in the aerospace industry to combine the loads resulting from the two sources. The resulting combined loads are then used to design for both yield ultimate strength and high cycle fatigue capability. This paper examines the cumulative distribution function (CDF) percentiles obtained using each method by integrating the joint probability density function of the sine and random components. A new Microsoft Excel spreadsheet macro that links with the software program Mathematics is then used to calculate the combined value corresponding to any desired percentile along with a curve fit to this value. Another Excel macro is used to calculate the combination using a Monte Carlo simulation. Unlike the traditional techniques, these methods quantify the calculated load value with a Consistent percentile. Using either of the presented methods can be extremely valuable in probabilistic design, which requires a statistical characterization of the loading. Also, since the CDF at high probability levels is very flat, the design value is extremely sensitive to the predetermined percentile; therefore, applying the new techniques can lower the design loading substantially without losing any of the identified structural reliability.

Brown, Andrew M.↗

Statistical error of reactor calculations by the Monte Carlo method

Algorithms for calculating the statistical error with allowance for intergenerational correlations are described. The algorithms are constructed on the basis of statistical analysis of the results of computations by the Monte Carlo method. As a result, simple rules for choosing the parameters of the computational techniques, such as the number of simulated generations necessary for attaining the required accuracy and the number of first skipped generations, are elaborated.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Simple and statistically sound recommendations for analysing physical theories

Physical theories that depend on many parameters or are tested against data from many different experiments pose unique challenges to statistical inference. Many models in particle physics, astrophysics and cosmology fall into one or both of these categories. These issues are often sidestepped with statistically unsound ad hoc methods, involving intersection of parameter intervals estimated by multiple experiments, and random or grid sampling of model parameters. Whilst these methods are easy to apply, they exhibit pathologies even in low-dimensional parameter spaces, and quickly become problematic to use and interpret in higher dimensions. In this article we give clear guidance for going beyond these procedures, suggesting where possible simple methods for performing statistically sound inference, and recommendations of readily-available software tools and standards that can assist in doing so. Our aim is to provide any physicists lacking comprehensive statistical training with recommendations for reaching correct scientific conclusions, with only a modest increase in analysis burden. Furthermore, our examples can be reproduced with the code publicly available at Zenodo.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Thermonuclear fusion rates for tritium + deuterium using Bayesian methods

The 3 H(d, n) 4 He reaction has a large low-energy cross section and will likely be utilized in future commercial fusion reactors. This reaction also takes place during Big Bang nucleosynthesis. Studies of both scenarios require accurate and precise fusion rates. To this end, we implement a one-level, two-channel R-matrix approximation into a Bayesian model. Our main goals are to predict reliable astrophysical S-factors and to estimate R-matrix parameters using the Bayesian approach. All relevant parameters are sampled in our study, including the channel radii, boundary condition parameters, and data set normalization factors. In addition, we take uncertainties in both measured bombarding energies and S-factors rigorously into account. Thermonuclear rates and reactivities of the 3 H(d, n) 4 He reaction are derived by numerically integrating the Bayesian S-factor samples. The present reaction rate uncertainties at temperatures between 1.0 MK and 1.0 GK are in the range of 0.2% to 0.6%. Our reaction rates differ from previous results by 2.9% near 1.0 GK.Our reactivities are smaller than previous results, with a maximum deviation of 2.9% near a thermal energy of 4 keV. The present rate or reactivity uncertainties are more reliable compared to previous studies that did not include the channel radii, boundary condition parameters, and data set normalization factors in the fitting. Finally, we investigate previous claims of electron screening effects in the published 3 H(d, n) 4 He data. No such effects are evident and only an upper limit for the electron screening potential can be obtained.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Statistical design for microwave systems

This paper presents an introduction to statistical system design. Basic ideas needed to understand statistical design and a method for implementing statistical design are presented. The nonlinear characteristics of the system amplifiers and mixers are accounted for in the given examples. The specification of group delay, signal-to-noise ratio and output power are considered in these statistical designs.

Cooke, Roland↗

Statistical retrieval techniques for sounding the meteorological structure of the atmosphere

Statistical retrieval methods for remote sounding are reviewed. Methods are given for constraining an essentially incomplete problem by means of the known statistical behavior of the solution. Information content of the observations and the meteorological structure are discussed. Linear versions of maximum probability and minimum variance methods are given in some detail, and extensions to the nonlinear case are described.

Rodgers, C. D.↗

Review and statistical analysis of the ultrasonic velocity method for estimating the porosity fraction in polycrystalline materials

A review and statistical analysis of the ultrasonic velocity method for estimating the porosity fraction in polycrystalline materials is presented. Initially, a semi-empirical model is developed showing the origin of the linear relationship between ultrasonic velocity and porosity fraction. Then, from a compilation of data produced by many researchers, scatter plots of velocity versus percent porosity data are shown for Al2O3, MgO, porcelain-based ceramics, PZT, SiC, Si3N4, steel, tungsten, UO2,(U0.30Pu0.70)C, and YBa2Cu3O(7-x). Linear regression analysis produced predicted slope, intercept, correlation coefficient, level of significance, and confidence interval statistics for the data. Velocity values predicted from regression analysis for fully-dense materials are in good agreement with those calculated from elastic properties.

Roth, D. J.↗

Fully iterative scatter corrected digital breast tomosynthesis using GPU-based fast Monte Carlo simulation and composition ratio update

Purpose: In digital breast tomosynthesis (DBT), scatter correction is highly desirable, as it improves image quality at low doses. Because the DBT detector panel is typically stationary during the source rotation, antiscatter grids are not generally compatible with DBT; thus, a software-based scatter correction is required. This work proposes a fully iterative scatter correction method that uses a novel fast Monte Carlo simulation (MCS) with a tissue-composition ratio estimation technique for DBT imaging. Methods: To apply MCS to scatter estimation, the material composition in each voxel should be known. To overcome the lack of prior accurate knowledge of tissue composition for DBT, a tissue-composition ratio is estimated based on the observation that the breast tissues are principally composed of adipose and glandular tissues. Using this approximation, the composition ratio can be estimated from the reconstructed attenuation coefficients, and the scatter distribution can then be estimated by MCS using the composition ratio. The scatter estimation and image reconstruction procedures can be performed iteratively until an acceptable accuracy is achieved. For practical use, (i) the authors have implemented a fast MCS using a graphics processing unit (GPU), (ii) the MCS is simplified to transport only x-rays in the energy range of 10–50 keV, modeling Rayleigh and Compton scattering and the photoelectric effect using the tissue-composition ratio of adipose and glandular tissues, and (iii) downsampling is used because the scatter distribution varies rather smoothly. Results: The authors have demonstrated that the proposed method can accurately estimate the scatter distribution, and that the contrast-to-noise ratio of the final reconstructed image is significantly improved. The authors validated the performance of the MCS by changing the tissue thickness, composition ratio, and x-ray energy. The authors confirmed that the tissue-composition ratio estimation was quite accurate under a variety of conditions. Our GPU-based fast MCS implementation took approximately 3 s to generate each angular projection for a 6 cm thick breast, which is believed to make this process acceptable for clinical applications. In addition, the clinical preferences of three radiologists were evaluated; the preference for the proposed method compared to the preference for the convolution-based method was statistically meaningful (p < 0.05, McNemar test). Conclusions: The proposed fully iterative scatter correction method and the GPU-based fast MCS using tissue-composition ratio estimation successfully improved the image quality within a reasonable computational time, which may potentially increase the clinical utility of DBT.

60 APPLIED LIFE SCIENCES↗

SeaWiFS Technical Report Series. Volume 7: Cloud screening for polar orbiting visible and infrared (IR) satellite sensors

Methods for detecting and screening cloud contamination from satellite derived visible and infrared data are reviewed in this document. The methods are applicable to past, present, and future polar orbiting satellite radiometers. Such instruments include the Coastal Zone Color Scanner (CZCS), operational from 1978 through 1986; the Advanced Very High Resolution Radiometer (AVHRR); the Sea-viewing Wide Field-of-view Sensor (SeaWiFS), scheduled for launch in August 1993; and the Moderate Resolution Imaging Spectrometer (IMODIS). Constant threshold methods are the least demanding computationally, and often provide adequate results. An improvement to these methods are the least demanding computationally, and often provide adequate results. An improvement to these methods is to determine the thresholds dynamically by adjusting them according to the areal and temporal distributions of the surrounding pixels. Spatial coherence methods set thresholds based on the expected spatial variability of the data. Other statistically derived methods and various combinations of basic methods are also reviewed. The complexity of the methods is ultimately limited by the computing resources. Finally, some criteria for evaluating cloud screening methods are discussed.

Darzi, Michael↗

Reaction analogy based forcing for incompressible scalar turbulence

Here, we present a novel reaction analogy (RA) based forcing method for generating statistically stationary scalar fields in incompressible turbulence. The new method can produce more general scalar probability density functions (PDFs), for example, quasi-double- δ PDF, than current methods, while ensuring that scalar fields remain bounded, unlike existent forcing methodologies that can potentially violate naturally existing bounds. Such features are useful for generating initial fields in nonpremixed combustion, inlet conditions for spatially developing flows, or for studying non-Gaussian scalar turbulence. The RA method mathematically models hypothetical chemical reactions that convert reactants in a mixed state back into its pure unmixed components. Various types of chemical reactions are formulated and the corresponding mathematical expressions derived such that the reaction term is smooth in scalar space and is consistent with mass conservation. For large values of the scalar forcing rate, the method produces statistically stationary quasi-double- δ scalar PDFs. Quasiuniform, Gaussian, and stretched exponential scalar statistics are recovered for smaller values of the scalar forcing rate. The shape of the scalar PDF can be further controlled by changing the stoichiometric coefficients of the reaction. Finally, the ability of the new method to produce fully developed passive scalar fields with quasi-Gaussian PDFs is also investigated, by exploring the convergence of the scalar variance spectrum to the Obukhov-Corrsin scaling and of the third-order mixed structure function to the “four-thirds” Yaglom's law.

42 ENGINEERING↗

How galaxies populate halos in extreme density environments: An analysis of the halo occupation distribution in SDSS

Recent works have shown that the properties of galaxy populations in simulated dark matter halos vary with large-scale environments. These results suggest a variation in the halo occupation distribution (HOD) in extreme density environments since the dynamical and astrophysical conditions prevailing in these regions may significantly affect the formation and evolution of their halos and residing galaxies, influencing the mean number of satellite galaxies. To analyse these effects, we identify cosmic voids and future virialised structures (FVS) in the Sloan Digital Sky Server Data Release 12 (SDSS-DR12) and estimate the HOD within these super-structures using group catalogues as dark matter halo proxies. Our goal is to use observational galaxy data to characterise the HOD within cosmic voids and FVS, explore the different properties of these galaxies populations and compare them with the general results outside these super-structures. We use a publicly available observational galaxy catalogue with information on redshifts, positions, magnitudes and other astrophysical features to build a volume complete galaxy sample and identify cosmic voids and FVS. Using a publicly available galaxy group catalogue as a proxy to dark matter halos, we compute the HOD within both types of super-structures for different absolute magnitude thresholds. We also study the dependence on the results on the main void and FVS properties, density and volume. We also analysed the main characteristics of the stellar content of galaxies inside these extreme-density regions such as the mean stellar age and the galaxy light concentration index. In all cases, we compare the results with those derived from the Field sample, defined by objects outside both types of environments. Inside cosmic voids, we find a strong decrease in HOD concerning the Field results. In the most extreme cases, the mean number of satellites fall to ~ 50%. Inside FVS, the HOD shows a significant increase to the Field, with a ~ 40% excess in the mean number of satellites. These results are present for the different galaxy luminosity ranges explored. In both environments, the differences with respect to the Field increases for the extreme values of the density environments. However, we obtain no signs of variations related to intrinsic characteristics of the super-structures, indicating that the effects depend mainly on the density of the large-scale environment. In addition, we find that the cumulative distribution of the mean age of stars of the central galaxy also varies in the different regions, this suggests that the history of the formation of the dark matter halos may be different. Finally, we explore the HOD for the 25% youngest (oldest) galaxies, based on the mean age of their stars. We find that for the low-mass groups the youngest galaxies are only present inside voids and are generally central galaxies. On the other hand, for the high-mass groups the FVS environments show the same increase in the HOD concerning the Field as previously mentioned. We find that cosmic voids lack a significant fraction of galaxies with the oldest stellar population.

79 ASTRONOMY AND ASTROPHYSICS↗

rascalc : a jackknife approach to estimating single- and multitracer galaxy covariance matrices

ABSTRACT To make use of clustering statistics from large cosmological surveys, accurate and precise covariance matrices are needed. We present a new code to estimate large-scale galaxy two-point correlation function (2PCF) covariances in arbitrary survey geometries that, due to new sampling techniques, runs ∼104 times faster than previous codes, computing finely binned covariance matrices with negligible noise in less than 100 CPU-hours. As in previous works, non-Gaussianity is approximated via a small rescaling of shot noise in the theoretical model, calibrated by comparing jackknife survey covariances to an associated jackknife model. The flexible code, rascalc, has been publicly released, and automatically takes care of all necessary pre- and post-processing, requiring only a single input data set (without a prior 2PCF model). Deviations between large-scale model covariances from a mock survey and those from a large suite of mocks are found to be indistinguishable from noise. In addition, the choice of input mock is shown to be irrelevant for desired noise levels below ∼105 mocks. Coupled with its generalization to multitracer data sets, this shows the algorithm to be an excellent tool for analysis, reducing the need for large numbers of mock simulations to be computed.

79 ASTRONOMY AND ASTROPHYSICS↗

Statistical dependency in visual scanning

A method to identify statistical dependencies in the positions of eye fixations is developed and applied to eye movement data from subjects who viewed dynamic displays of air traffic and judged future relative position of aircraft. Analysis of approximately 23,000 fixations on points of interest on the display identified statistical dependencies in scanning that were independent of the physical placement of the points of interest. Identification of these dependencies is inconsistent with random-sampling-based theories used to model visual search and information seeking.

Ellis, Stephen R.↗

Advanced Gear Alloys for Ultra High Strength Applications

Single tooth bending fatigue (STBF) test data of UHS Ferrium C61 and C64 alloys are presented in comparison with historical test data of conventional gear steels (9310 and Pyrowear 53) with comparable statistical analysis methods. Pitting and scoring tests of C61 and C64 are works in progress. Boeing statistical analysis of STBF test data for the four gear steels (C61, C64, 9310 and Pyrowear 53) indicates that the UHS grades exhibit increases in fatigue strength in the low cycle fatigue (LCF) regime. In the high cycle fatigue (HCF) regime, the UHS steels exhibit better mean fatigue strength endurance limit behavior (particularly as compared to Pyrowear 53). However, due to considerable scatter in the UHS test data, the anticipated overall benefits of the UHS grades in bending fatigue have not been fully demonstrated. Based on all the test data and on Boeing s analysis, C61 has been selected by Boeing as the gear steel for the final ERDS demonstrator test gearboxes. In terms of potential follow-up work, detailed physics-based, micromechanical analysis and modeling of the fatigue data would allow for a better understanding of the causes of the experimental scatter, and of the transition from high-stress LCF (surface-dominated) to low-stress HCF (subsurface-dominated) fatigue failure. Additional STBF test data and failure analysis work, particularly in the HCF regime and around the endurance limit stress, could allow for better statistical confidence and could reduce the observed effects of experimental test scatter. Finally, the need for further optimization of the residual compressive stress profiles of the UHS steels (resulting from carburization and peening) is noted, particularly for the case of the higher hardness C64 material.

Shen, Tony↗