Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “applied 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,027 records · Page 57
All sourcesNASA NTRSDOE OSTI
NASA NTRS2019

Evaluation of the Impacts of Biomass Burning Aerosols on the NASA GEOS Sub-Seasonal Climate Forecasts

Emission sources of trace gases and aerosol particles in the South American (SA)and African (Af) continents have a strong seasonal and space variability associated with the extensive vegetation fires activities. In both continents, during the austral winter, the fires affect mainly tropical forest and savannah-type biomes and are mostly associated with deforestation and agricultural/pasture land management. Smoke aerosol particles, on average, contribute to at least 90% of the total aerosol optical depth (AOD) in the visible spectrum in the case of the South America regional smoke. Smoke aerosols also act as cloud condensation nuclei affecting cloud microphysics properties and therefore, changing the radiation budget, hydrological cycle and global circulation patterns over disturbed areas (Kaufman, 1995; Rosenfeld, 1999; Andreae,et al., 2004; Koren et al., 2004, Zhang, 2008; Ott et al., 2010; Randles et al., 2013). This study aims to evaluate and quantify the impact of including a comprehensive emission field of biomass burning aerosol on the performance of a seasonal climate forecast system, not only regarding the AOD itself but mainly on the meteorological state variable (e.g., precipitation and temperature). To address the questions put above, we designed two numerical experiments: 1- named"AERO_CTL" which applies the Quick Fire Emissions Dataset (QFED) emissions estimated with intra-diurnal variation (hereafter, BBE), and 2- named "AERO_CLM" where the sourcee mission is based on a climatology of the QFED emissions, with only monthly variation(hereafter, BBCLIM). Hindcast simulations were produced using the Goddard Earth ObservingSystem global circulation model, version 5, sub-seasonal to seasonal (GEOS5-S2S) system with a nominal spatial resolution of 56km (Rienecker et al., 2008). In both experiments, the aerosol feedbacks from cloud developments and radiation interactions were accounted. The two experiments consisted of 4 members each and ran from June to November spanning over the years 2000 to 2015. Model performance was evaluated by calculating statistical metrics on the mean area of SA and Af. Our results demonstrated that the skill model in predicting AOD is significantly improve when BBE source emission is applied over SA, but not over the Afcontinent. Over SA, the correlation between the AERO_CTL model configuration and MERRA-2 is 0.93 (R2= 0.86, RMS=0.02, BIAS=0.01), while the AERO_CLM model presents a value of0.81 (R2= 0.65, RMS=0.04, BIAS=0.06). However, the AERO_CTL experiment better represents the inter-annual variability of the AOS in both regions. The gain of the skill in predicting the AOD by the AERO_CTL experiment is also seen in some meteorological variables. We observed an increase in the model skill in predicting the 2-meter temperature and precipitation of up to 0.3 for the AERO_CTL experiment in comparison to the AERO_CLM. AERO_CLM. According to the analyzed hindcast, we inferred that representing the BBE more realistically implies in a significant gain of skills in the seasonal climate forecasting over SA and Af continents.

Aerosol↗
DOE OSTI2022

Freely jointed chain models with extensible links

We report analytical relations for the mechanical response of single polymer chains are valuable for modeling purposes, on both the molecular and the continuum scale. These relations can be obtained using statistical thermodynamics and an idealized single-chain model, such as the freely jointed chain model. To include bond stretching, the rigid links in the freely jointed chain model can be made extensible, but this almost always renders the model analytically intractable. Here, an asymptotically correct statistical thermodynamic theory is used to develop analytic approximations for the single-chain mechanical response of this model. The accuracy of these approximations is demonstrated using several link potential energy functions. This approach can be applied to other single-chain models, and to molecular stretching in general.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2020

Advancing Methodologies for Applying Machine Learning and Evaluating Spatiotemporal Models of Fine Particulate Matter (PM 2.5 ) Using Satellite Data Over Large Regions

Reconstructing the distribution of fine particulate matter (PM 2.5 ) in space and time, even far from ground monitoring sites, is an important exposure science contribution to epidemiologic analyses of PM 2.5 health impacts. Flexible statistical methods for prediction have demonstrated the integration of satellite observations with other predictors, yet these algorithms are susceptible to overfitting the spatiotemporal structure of the training datasets. We present a new approach for predicting PM 2.5 using machine-learning methods and evaluating prediction models for the goal of making predictions where they were not previously available. We apply extreme gradient boosting (XGBoost) modeling to predict daily PM 2.5 on a 1 x 1 km 2 resolution for a 13 state region in the Northeastern USA for the years 2000–2015 using satellite-derived aerosol optical depth and implement a recursive feature selection to develop a parsimonious model. We demonstrate excellent predictions of withheld observations but also contrast an RMSE of 3.11 μg/m 3 in our spatial cross-validation withholding nearby sites versus an overfit RMSE of 2.10 μg/m 3 using a more conventional random ten-fold splitting of the dataset. As the field of exposure science moves forward with the use of advanced machine-learning approaches for spatiotemporal modeling of air pollutants, our results show the importance of addressing data leakage in training, overfitting to spatiotemporal structure, and the impact of the predominance of ground monitoring sites in dense urban sub-networks on model evaluation. The strengths of our resultant modeling approach for exposure in epidemiologic studies of PM 2.5 include improved efficiency, parsimony, and interpretability with robust validation while still accommodating complex spatiotemporal relationships.

air pollution↗
DOE OSTI2026

Enhancing ZFP: A Statistical Approach to Understanding and Reducing Error Bias in a Lossy Floating-Point Compression Algorithm

The amount of data generated and gathered in scientific simulations and data collection applications is continuously growing, putting mounting pressure on storage and bandwidth concerns. A means of reducing such issues is data compression; but, lossless data compression is typically ineffective when applied to floating-point data. Thus, users tend to apply a lossy data compressor, which allows for small deviations from the original data. It is essential to understand how the error from lossy compression impacts the accuracy of the data analytics. Thus, we must analyze not only the compression properties but the error as well. In this paper, we provide a statistical analysis of the error caused by ZFP compression, a state-of-the-art, lossy compression algorithm explicitly designed for floating-point data. We show that the error is indeed biased and propose simple modifications to the algorithm to neutralize the bias and further reduce the resulting error.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2014

A Reliability Model for Ni-BaTiO3-Based (BME) Ceramic Capacitors

The evaluation of multilayer ceramic capacitors (MLCCs) with base-metal electrodes (BMEs) for potential NASA space project applications requires an in-depth understanding of their reliability. The reliability of an MLCC is defined as the ability of the dielectric material to retain its insulating properties under stated environmental and operational conditions for a specified period of time t. In this presentation, a general mathematic expression of a reliability model for a BME MLCC is developed and discussed. The reliability model consists of three parts: (1) a statistical distribution that describes the individual variation of properties in a test group of samples (Weibull, log normal, normal, etc.), (2) an acceleration function that describes how a capacitors reliability responds to external stresses such as applied voltage and temperature (All units in the test group should follow the same acceleration function if they share the same failure mode, independent of individual units), and (3) the effect and contribution of the structural and constructional characteristics of a multilayer capacitor device, such as the number of dielectric layers N, dielectric thickness d, average grain size r, and capacitor chip size S. In general, a two-parameter Weibull statistical distribution model is used in the description of a BME capacitors reliability as a function of time. The acceleration function that relates a capacitors reliability to external stresses is dependent on the failure mode. Two failure modes have been identified in BME MLCCs: catastrophic and slow degradation. A catastrophic failure is characterized by a time-accelerating increase in leakage current that is mainly due to existing processing defects (voids, cracks, delamination, etc.), or the extrinsic defects. A slow degradation failure is characterized by a near-linear increase in leakage current against the stress time; this is caused by the electromigration of oxygen vacancies (intrinsic defects). The two identified failure modes follow different acceleration functions. Catastrophic failures follow the traditional power-law relationship to the applied voltage. Slow degradation failures fit well to an exponential law relationship to the applied electrical field. Finally, the impact of capacitor structure on the reliability of BME capacitors is discussed with respect to the number of dielectric layers in an MLCC unit, the number of BaTiO3 grains per dielectric layer, and the chip size of the capacitor device.

reliability↗
DOE OSTI2025

Testing convolutional neural network based deep learning systems: a statistical metamorphic approach

Machine learning technology spans many areas and today plays a significant role in addressing a wide range of problems in critical domains,i.e., healthcare, autonomous driving, finance, manufacturing, cybersecurity,etc. Metamorphic testing (MT) is considered a simple but very powerful approach in testing such computationally complex systems for which either an oracle is not available or is available but difficult to apply. Conventional metamorphic testing techniques have certain limitations in verifying deep learning-based models (i.e., convolutional neural networks (CNNs)) that have a stochastic nature (because of randomly initializing the network weights) in their training. In this article, we attempt to address this problem by using a statistical metamorphic testing (SMT) technique that does not require software testers to worry about fixing the random seeds (to get deterministic results) to verify the metamorphic relations (MRs). We propose seven MRs combined with different statistical methods to statistically verify whether the program under test adheres to the relation(s) specified in the MR(s). We further use mutation testing techniques to show the usefulness of the proposed approach in the healthcare space and test two CNN-based deep learning models (used for pneumonia detection among patients). The empirical results show that our proposed approach uncovers 85.71% of the implementation faults in the classifiers under test (CUT). Furthermore, we also propose an MRs minimization algorithm for the CUT, thus saving computational costs and organizational testing resources.

Computer Science↗
NASA NTRS1989

Rain core structure statistics derived from radar and disdrometer measurements in the mid-Atlantic coast of the US

During a period spanning more than 5 years, low elevation radar measurements of rain were systematically obtained in the mid-Atlantic coast of the U.S. Drop size distribution measurements with a disdrometer were also acquired on the same rain days. The drop size data were utilized to convert the radar reflectivity factors to estimated rain rates for the respective rain days of operation. Applying high level algorithms to the rain data, core values of rain intensities were identified (peak rain rates), and families of rain rate isopleths analyzed. In particular, equicircle diameters of the family of isopleths enveloping peak rain intensities were statistically characterized. The presented results represents the analysis of two rain days, 12 radar scans, corresponding to 430 culled rain rate isopleths from an available data base of 22,000 contours, approximately 100 scans encompassing 17 rain days. The results presented show trends of the average rain rate vs. contour scale dimensions, and cumulative distributions of rain cell dimensions which belong to core families of precipitation.

Goldhirsh, Julius↗
DOE OSTI2022

Recurrent neural networks for short-term and long-term prediction of geothermal reservoirs

Accurate prediction of geothermal reservoir responses to alternative energy production scenarios is critical for optimizing the development of the underlying resources. While the conventional physics-based models offer a comprehensive prediction tool, data-driven models provide an efficient alternative to build fit-for-purpose predictive models by extracting and using the statistical patterns in the collected data to make predictions. The recurrent neural network (RNN) is a data-driven model that is commonly applied to predict time series sequences. This paper presents a variant of RNN that also utilizes the efficiency of convolutional neural networks (CNN) for the prediction of energy production from geothermal reservoirs. Specifically, a CNN–RNN architecture is developed that takes historical well controls as input (features) and their corresponding production response data as output (labels) to learn an input-output mapping that can predict the future well production responses/performance for any given future well control inputs. The model is paired with a labeling scheme to handle real field disturbances that create data gaps. In addition to the model structure, we introduce a thorough workflow for applying the model, which includes data pre-processing, feature selection, as well as different training strategies for short-term and long-term prediction. Finally, the performance and accuracy of the model are evaluated by applying it to multiple datasets, including a field reservoir model.

15 GEOTHERMAL ENERGY↗
DOE OSTI2025

Going Off Grid: A Comparative Study of the Lagrangian and Eulerian Perspectives of New Particle Formation Events

New particle formation and growth (NPF&G) is the process by which ultrafine particles are formed from gas-phase precursors. NPF&G is the dominant source of global aerosol number with important influences on climate. Most observations of NPF&G events are conducted at stationary sites; however, NPF&G observed from stationary sites is influenced by gradual or rapid changes in the air masses passing over the site, complicating NPF&G analysis. In this work, we use observations and a 3D aerosol model to compare aerosol size distributions at a stationary site (Southern Great Plains [SGP] observatory, Oklahoma, USA) and along Lagrangian trajectories crossing the site. The model simulates the NPF&G events reasonably well at SGP. Using the model to compare the Lagrangian and stationary perspectives, we can explain previously unanalyzable days with some evidence of NPF&G as either non-event or analyzable NPF&G days. We find most of the unanalyzable NPF&G days are due to isolated and inhomogeneous NPF&G occurring upwind of the stationary site, often in the outflow of urban regions. Finally, we compare formation rates of 3 nm particles, growth rates, and the survival probability of 3 nm particles growing to 25 nm between the stationary and Lagrangian perspectives. Because of the much larger number of analyzable days along the Lagrangian trajectories, this perspective potentially provides more robust statistics and better characterization of NPF&G event extremes. Our method for extracting chemical/physical properties along Lagrangian trajectories from 3D models can be applied to a wide range of science questions.

O’Donnell, Samuel E. [Colorado State Univ., Fort C↗
DOE OSTI2023

Self-testing of a single quantum system from theory to experiment

Abstract Self-testing allows one to characterise quantum systems under minimal assumptions. However, existing schemes rely on quantum nonlocality and cannot be applied to systems that are not entangled. Here, we introduce a robust method that achieves self-testing of individual systems by taking advantage of contextuality. The scheme is based on the simplest contextuality witness for the simplest contextual quantum system—the Klyachko-Can-Binicioğlu-Shumovsky inequality for the qutrit. We establish a lower bound on the fidelity of the state and the measurements as a function of the value of the witness under a pragmatic assumption on the measurements. We apply the method in an experiment on a single trapped 40 Ca + using randomly chosen measurements and perfect detection efficiency. Using the observed statistics, we obtain an experimental demonstration of self-testing of a single quantum system.

Physics↗
DOE OSTI2021

Statistical tools for a better optical model

Background: Modern statistical tools provide the ability to compare the information content of observables and provide a path to explore which experiments would be most useful to give insight into and constrain theoretical models. Purpose: Here we study three such tools in the context of nuclear reactions with the goal of constraining the optical potential. Method: The three statistical tools examined are (i) the principal component analysis, (ii) the sensitivity analysis based on derivatives, and (iii) the Bayesian evidence. We first apply these tools to a toy-model case, comparing the form of the imaginary part of the optical potential. Then we consider two different reaction observables, elastic angular distributions and polarization data for reactions on 48 Ca and 208 Pb at two different beam energies. Results: For the toy-model case, we find significant discrimination power in the sensitivities and the Bayesian evidence, showing clearly that the volume imaginary term is more useful to describe scattering at higher energies. When comparing between elastic cross sections and polarization data using realistic optical models, sensitivity studies indicate that both observables are roughly equally sensitive but the variability of the optical model parameters is strongly angle dependent. The Bayesian evidence shows some variability between the two observables, but the Bayes factor obtained is not sufficient to discriminate between angular distributions and polarization. Conclusions: From the cases considered, we conclude that, in general, elastic scattering angular distributions have similar impact in constraining the optical potential parameters compared with the polarization data. The angular ranges for the optimum experimental constraints can vary significantly with the observable considered.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Experimental Measurements of the Wake of a Sphere at Subcritical Reynolds Numbers

This work experimentally investigated the flow phenomena and vortex structures in the wake of a sphere located in a water loop at Reynolds numbers of Re = 850, 1,250, and 1,700. Velocity fields in the wake region were obtained by applying the time-resolved stereoscopic particle image velocimetry (TR-SPIV) technique. From the acquired TR-SPIV velocity vector fields, the statistical values of mean and fluctuating velocities were computed. Spectral analysis, two-point velocity–velocity cross-correlation, proper orthogonal decomposition (POD) and vortex identification analyses were also performed. The velocity fields show a recirculation region that decreases in length with an increase of Reynolds numbers. The power spectra from the spectral analysis had peaks corresponding to a Strouhal number of St = 0.2, which is a value commonly found in the literature studies of flow over a sphere. The two-point cross-correlation analysis revealed elliptical structures in the wake, with estimated integral length scales ranging between 12% and 63% of the sphere diameter. The POD analysis revealed the statistically dominant flow structures that captured the most flow kinetic energy. It is seen that the flow kinetic energy captured in the smaller scale flow structures increased as Reynolds number increased. The POD modes contained smaller structure as the Reynolds number increased and as mode order increased. In addition, spectral analysis performed on the POD temporal coefficients revealed peaks corresponding to St = 0.2, similar to the spectral analysis on the fluctuating velocity. The ability of POD to produce low-order reconstructions of the flow was also utilized to facilitate vortex identification analysis, which identified average vortex sizes of 0.41D for Re1, 0.33D for Re2, and 0.32D for Re3.

Engineering↗
DOE OSTI2020

Systematic investigation of the particle spectra in heavy-ion collisions at the Large Hadron Collider

We investigate the charged particle spectra produced in the heavy-ion collisions at nine centralities from different systems, i.e. Pb+Pb at √ s NN = 2.76 TeV and 5.02 TeV as well as Xe+Xe at √ s NN = 5.44 TeV, at Large Hadron Collider (LHC) using one empirical formula inspired by the stationary solution of the Fokker-Planck equation, dubbed as the generalized Fokker-Planck solution (GFPS). Our results show that the GFPS can reproduce the experimental particle spectrum up to transverse momentum pT about 45 GeV/c with the maximum discrepancy 30% covering 10 orders of magnitude. Here, the discrepancy between the data and the results from the GFPS decreases to 15% when the maximum of the charged particle transverse momentum is cut to 20 GeV/c. We confirmed that the Tsallis distribution derived from the non-extensive statistics, which can reproduce the particle spectra produced in small collision systems, such as p+p, up to few hundreds GeV/c, can only apply to systematically study the particle spectra up to 8 GeV/c in A+A collisions at LHC, as pointed out in the study of identified particle spectra in Pb+Pb collisions at √ s NN = 2.76 TeV. A brief discussion on GFPS is also given.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Comparison of All Solid Cancer Mortality and Incidence Dose-Response in the Life Span Study of Atomic Bomb Survivors, 1958–2009

Recent analysis of all solid cancer incidence (1958–2009) in the Life Span Study (LSS) revealed evidence of upward curvature in the radiation dose response among males but not females. Upward curvature in sex-averaged excess relative risk (ERR) for all solid cancer mortality (1950–2003) was also observed in the 0–2 Gy dose range. As reasons for non-linearity in the LSS are not completely understood, we conducted dose response analyses for all solid cancer mortality and incidence applying similar methods (1958–2009 follow-up, DS02R1 doses, including subjects notin-city (NIC) at the time of the bombing) and statistical models. Incident cancers were ascertained from Hiroshima and Nagasaki cancer registries, while cause of death was ascertained from death certificates over entire Japan. The study included 105,444 LSS subjects who were alive and not known to have cancer before Jan 1, 1958 (80,205 with dose estimates and 25,239 NIC subjects). Between 1958 and 2009, there were 3.1 million person-years (PY) and 22,538 solid cancers for incidence analysis and 3.8 million PY and 15,419 solid cancer deaths for mortality analysis. We fitted sex-specific ERR models adjusted for smoking to both types of data. Over the entire range of doses, solid cancer mortality dose response exhibited a borderline significant upward curvature among males (P=0.062) and significant upward curvature among females (P=0.010); for solid cancer incidence, as before, we found a significant upward curvature among males (P=0.001) but not among females (P=0.624). The sex difference in magnitude of dose response curvature was statistically significant for cancer incidence (P=0.017) but not for cancer mortality (P=0.781). The results of analyses in the 0–2 Gy range and restricted lower dose ranges generally supported inferences made about the sex-specific dose response shape over the entire range of doses for each outcome. Patterns of sex-specific curvature by calendar period (1958–1987 vs 1988–2009) and age at exposure (0–19 vs 20–83) varied between mortality and incidence data, particularly among females, although for each outcome there was an indication of curvature among 0–19 year old male survivors in both calendar periods and among 0–19 year old female survivors in the recent period. Collectively, our findings indicate that the upward curvature in all solid cancer dose response in the LSS is neither specific to males nor to incidence data; it appears to depend on composition of case series and age at exposure or time. Further follow-up and site-specific analyses of cancer mortality and incidence will be important to confirm the emerging trend in dose response curvature among young survivors and unveil the contributing factors and sites.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗
DOE OSTI2021

Estimating the Adequacy of a Multi-Objective Optimization

Multi-objective optimization methods can be criticized for lacking a statistically valid measure of the quality and representativeness of a solution. This stance is especially relevant to metaheuristic optimization approaches but can also apply to other methods that typically might only report a small representative subset of a Pareto frontier. Here we present a method to address this deficiency based on random sampling of a solution space to determine, with a specified level of confidence, the fraction of the solution space that is surpassed by an optimization. The Superiority of Multi-Objective Optimization to Random Sampling, or SMORS method, can evaluate quality and representativeness using dominance or other measures, e.g., a spacing measure for high-dimensional spaces. SMORS has been tested in a combinatorial optimization context using a genetic algorithm but could be useful for other optimization methods.

42 ENGINEERING↗
DOE OSTI2026

DESI Strong Lens Foundry. I. HST Observations and Modeling with GIGA-Lens

We present the Dark Energy Spectroscopic Instrument (DESI) Strong Lens Foundry. We discovered ∼3500 new strong gravitational lens candidates in the DESI Legacy Imaging Surveys using residual neural networks (ResNet). We observed a subset (51) of our candidates using the Hubble Space Telescope (HST). Except for one ambiguous case, we have confirmed 50 of the 51 candidates to be strong lenses. We also briefly describe spectroscopic follow-up observations by DESI and Keck NIRES programs. From this very rich data set, a number of studies will be carried out, including evaluating the quality of the ResNet search candidates and lens modeling. In this paper, we present our initial effort in these directions. In particular, as a demonstration, we present the lens model for DESI-165.4754−06.0423, with imaging data from HST, and lens and source redshifts from DESI and Keck NIRES, respectively. In this effort, we have applied a fully forward-modeling Bayesian approach (GIGA-Lens), using multiple GPUs, to a strong lens with HST data, and achieved statistical convergence.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1971

Statistical studies in stellar rotation 2: A method of analyzing rotational coupling in double stars and an introduction to its applications

The correlation between the equatorial velocities of the components of double stars is studied from a statistical standpoint. A theory of rotational correlation is developed and discussed with regard to its applicability to existing observations. The theory is then applied to a sample of visual binaries which are the least studied for rotational coupling. Consideration of eclipsing systems and spectroscopic binaries is limited to show how the degrees of freedom in the spin parallelism problem can be reduced. The analysis lends support to the existence of synchronism in closely spaced binaries.

Bernacca, P. L.↗
NASA NTRS1976

Nondeterministic data base for computerized visual perception

A description is given of the knowledge representation data base in the perception subsystem of the Mars robot vehicle prototype. Two types of information are stored. The first is generic information that represents general rules that are conformed to by structures in the expected environments. The second kind of information is a specific description of a structure, i.e., the properties and relations of objects in the specific case being analyzed. The generic knowledge is represented so that it can be applied to extract and infer the description of specific structures. The generic model of the rules is substantially a Bayesian representation of the statistics of the environment, which means it is geared to representation of nondeterministic rules relating properties of, and relations between, objects. The description of a specific structure is also nondeterministic in the sense that all properties and relations may take a range of values with an associated probability distribution.

Yakimovsky, Y.↗