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 361 records · Page 20

Structure classification and melting temperature prediction in octet AB solids via machine learning

Machine learning methods are being increasingly used in condensed matter physics and materials science to classify crystals structures and predict material properties. However, the reliability of these methods for a given problem, especially when large data sets are unavailable, has not been well studied. By addressing the tasks of classifying crystal structure and predicting melting temperatures of the octet subset of AB solids, we performed such a study and found potential problems with using machine learning methods on relatively small data sets. At the same time, however, we can reaffirm the potential power of such methods for these tasks. In particular, we uncovered an important new material feature, the excess Born effective charge, that significantly increased the accuracy of the predictions for the classification problem we defined. This discovery leads us to propose a new scale for the degree of ionicity and covalency in these solids. More specifically, we partitioned the crystal structures of a set of 75 octet solids into those that are ionic and covalent bonded and thus performed a binary classification task. We found that using the standard indices (r σ ,r π ), suggested by St. John and Bloch several decades ago, enabled an average success in classification of 92%. Using just r σ and the excess Born effective charge ΔZ A of the A atom enabled an average success of 97%, but we also found relatively large variations about these averages that were dependent on how certain machine learning methods were used and for which a standard deviation was not a proper measure of the degree of confidence we can place in either average. Instead, we calculated and report with 95 % confidence that the traditional classification pair predicts an accuracy in the interval [ 89%, 95%] and the accuracy of the new pair lies in the interval [96 %, 99%]. For melting temperature predictions, the size of our data set was 46. We estimate the root-mean-squared error of our resulting model to be 11% of the mean melting temperature of the data, but we note that if the accuracy of this predicted error is itself measured, our estimated fitting error itself has a root-mean-square error of 50%. In short, what we illustrate is that classification and regression predictions can vary significantly, depending on the details of how machine learning methods are applied to small data sets. This variation makes it important, if not essential, to average the predictions and compute confidence intervals about these averages to report results meaningfully. However, when properly used, these statistical methods can advance our understanding and improve predictions of material properties even for small data sets.

36 MATERIALS SCIENCE↗

Seeing through noise in power laws

Despite widespread claims of power laws across the natural and social sciences, evidence in data is often equivocal. Modern data and statistical methods reject even classic power laws such as Pareto’s law of wealth and the Gutenberg–Richter law for earthquake magnitudes. We show that the maximum-likelihood estimators and Kolmogorov–Smirnov (K-S) statistics in widespread use are unexpectedly sensitive to ubiquitous errors in data such as measurement noise, quantization noise, heaping and censorship of small values. This sensitivity causes spurious rejection of power laws and biases parameter estimates even in arbitrarily large samples, which explains inconsistencies between theory and data. We show that logarithmic binning by powers of λ > 1 attenuates these errors in a manner analogous to noise averaging in normal statistics and that λ thereby tunes a trade-off between accuracy and precision in estimation. Binning also removes potentially misleading within-scale information while preserving information about the shape of a distribution over powers of λ, and we show that some amount of binning can improve sensitivity and specificity of K-S tests without any cost, while more extreme binning tunes a trade-off between sensitivity and specificity. We therefore advocate logarithmic binning as a simple essential step in power-law inference.

97 MATHEMATICS AND COMPUTING↗

Probabilistic-learning-based stochastic surrogate model from small incomplete datasets for nonlinear dynamical systems

We consider a high-dimensional nonlinear computational model of a dynamical system, parameterized by a vector-valued control parameter, in the presence of uncertainties represented by an uncontrolled parameter modeled by a vector-valued random variable, and possibly with stochastic excitation. The objective is to construct a statistical surrogate model where the input is any deterministic value of the control parameter, and the output is a vector-valued observation of the computational model, which is a random vector whose probability measure is updated using a target dataset. To construct this statistical surrogate model, the stochastic response of the computational model must be built, which is a vector-valued time-discretized stochastic process in high dimension, depending on the control parameter. It is assumed that the computational cost of a single evaluation of the deterministic model is high. For the probabilistic updating, we consider a subset of the components of the observation of the computational model, defined as the “identification observation” of the computational model, for which a small target dataset is available. Therefore, the target dataset is associated with partial observability, corresponding to an incomplete data case. Given a prior probability model of the random control and uncontrolled parameters, a training dataset is constructed, consisting of realizations of the random triplet composed of the stochastic response, the random identification observation, and the random control parameter. Since the computational cost of a single evaluation of the deterministic model is assumed to be large, the training dataset is also of small size. The main challenges in this problem are the high dimensionality, partial observability leading to incomplete data in the target dataset for the identification observation of the computational model (which is not sufficient to identify the computational stochastic responses), and the availability of a small training dataset. To address these challenges, we propose a methodology based on statistical methods for constructing necessary reduced representations, direct probabilistic learning under constraints using probabilistic learning on manifolds (PLoM) constrained by the target dataset, and the use of a weak formulation of the Fourier transform of probability measures. Statistical conditioning is also employed to explore the learned dataset. The constructed predictive statistical surrogate model can be implemented in the context of online computation. Here, we apply this approach to a problem of nonlinear stochastic dynamics in high dimensions within the framework of deformable solids mechanics.

Engineering↗

Deciphering the source of primary biological aerosol particles: a pollen case study

Primary biological aerosol particles (PBAPs) are microscopic solids suspended in the atmosphere emitted by biological systems and play critical roles in the atmosphere and at the atmosphere-biosphere interface, impacting human health, climate, and the ecosystem function. Understanding the sources of PBAPs is necessary to decipher the mechanistic interactions between aerosols, climate, and distinct ecosystem components. However, the detection of specific PBAPs in complex ambient aerosol samples is challenging. We performed metabolomics analyses of pollen from three pollinating tree species and ambient samples collected during the peak pollination period of each species. Random Forest and sPLS-DA machine learning methods were employed to evaluate whether metabolic signatures of ambient samples can reveal the source of the main pollen particles present in the atmosphere. Our results suggest that atmospheric eco-metabolomics techniques combined with sophisticated statistical methods can decipher the origin of abundant PBAPs from complex ambient samples. Developing complete libraries containing high-resolution metabolomic fingerprints of the major PBAPs present in the atmosphere would significantly advance future research to accurately understand the role of PBAPs in the atmosphere, ecosystems and human health.

Rivas-Ubach, Albert↗

Source Data Applicability Impacts on Epistemic Uncertainty for Launch Vehicle Fault Tree Models

Launch vehicle systems are designed and developed using both heritage and new hardware. Design modifications to the heritage hardware to fit new functional system requirements can impact the applicability of heritage reliability data. Risk estimates for newly designed systems must be developed from generic data sources such as commercially available reliability databases using reliability prediction methodologies, such as those addressed in MIL-HDBK-217F. Failure estimates must be converted from the generic environment to the specific operating environment of the system where it is used. In addition, some qualification of applicability for the data source to the current system should be made. Characterizing data applicability under these circumstances is crucial to developing model estimations that support confident decisions on design changes and trade studies. This paper will demonstrate a data-source applicability classification method for assigning uncertainty to a target vehicle based on the source and operating environment of the originating data. The source applicability is determined using heuristic guidelines while translation of operating environments is accomplished by applying statistical methods to MIL-HDK-217F tables. The paper will provide a case study example by translating Ground Benign (GB) and Ground Mobile (GM) to the Airborne Uninhabited Fighter (AUF) environment for three electronic components often found in space launch vehicle control systems. The classification method will be followed by uncertainty-importance routines to assess the need to for more applicable data to reduce uncertainty.

Al Hassan, Mohammad↗

Phase Space Dissimilarity Measures for Structural Health Monitoring

A novel method for structural health monitoring (SHM), known as the Phase Space Dissimilarity Measures (PSDM) approach, is proposed and developed. The patented PSDM approach has already been developed and demonstrated for a variety of equipment and biomedical applications. Here, we investigate SHM of bridges via analysis of time serial accelerometer measurements. This work has four aspects. The first is algorithm scalability, which was found to scale linearly from one processing core to four cores. Second, the same data are analyzed to determine how the use of the PSDM approach affects sensor placement. We found that a relatively low-density placement sufficiently captures the dynamics of the structure. Third, the same data are analyzed by unique combinations of accelerometer axes (vertical, longitudinal, and lateral with respect to the bridge) to determine how the choice of axes affects the analysis. The vertical axis is found to provide satisfactory SHM data. Fourth, statistical methods were investigated to validate the PSDM approach for this application, yielding statistically significant results.

47 OTHER INSTRUMENTATION↗

Statistical Analysis of Factors Riving Surface Ozone Variability over Continental South Africa

Statistical relationships between surface ozone (O3) concentration, precursor species and meteorological conditions in continental South Africa were examined from data obtained from measurement stations in north-eastern South Africa. Three multivariate statistical methods were applied in the investigation, i.e. multiple linear regression (MLR), principal component analysis (PCA) and –regression (PCR), and generalised additive model (GAM) analysis. The daily maximum 8-h moving average O3 concentrations were considered in these statistical models (dependent variable). MLR models indicated that meteorology and precursor species concentrations are able to explain ~50% of the variability in daily maximum O3 levels. MLR analysis revealed that atmospheric carbon monoxide (CO), temperature and relative humidity were the strongest factors affecting the daily O3 variability. In summer, daily O3 variances were mostly associated with relative humidity, while winter O3 levels were mostly linked to temperature and CO. PCA indicated that CO, temperature and relative humidity were not strongly collinear. GAM also identified CO, temperature and relative humidity as the strongest factors affecting the daily variation of O3. Partial residual plots found that temperature, radiation and nitrogen oxides most likely have a non-linear relationship with O3,while the relationship with relative humidity and CO is probably linear. An inter-comparison between O3 levels modelled with the three statistical models compared to measured O3 concentrations showed that the GAM model offered a slight improvement over the MLR model. These findings emphasise the critical role of regional-scale O3 precursors coupled with meteorological conditions in daily variances of O3 levels in continental South Africa.

multiple linear regression (MLR)↗

Uncertainty Management in Remote Sensing of Climate Data. Summary of A Workshop

Great advances have been made in our understanding of the climate system over the past few decades, and remotely sensed data have played a key role in supporting many of these advances. Improvements in satellites and in computational and data-handling techniques have yielded high quality, readily accessible data. However, rapid increases in data volume have also led to large and complex datasets that pose significant challenges in data analysis (NRC, 2007). Uncertainty characterization is needed for every satellite mission and scientists continue to be challenged by the need to reduce the uncertainty in remotely sensed climate records and projections. The approaches currently used to quantify the uncertainty in remotely sensed data, including statistical methods used to calibrate and validate satellite instruments, lack an overall mathematically based framework.

Remote sensing↗

Modeling Ka-band low elevation angle propagation statistics

The statistical variability of the secondary atmospheric propagation effects on satellite communications cannot be ignored at frequencies of 20 GHz or higher, particularly if the propagation margin allocation is such that link availability falls below 99 percent. The secondary effects considered in this paper are gaseous absorption, cloud absorption, and tropospheric scintillation; rain attenuation is the primary effect. Techniques and example results are presented for estimation of the overall combined impact of the atmosphere on satellite communications reliability. Statistical methods are employed throughout and the most widely accepted models for the individual effects are used wherever possible. The degree of correlation between the effects is addressed and some bounds on the expected variability in the combined effects statistics are derived from the expected variability in correlation. Example estimates are presented of combined effects statistics in the Washington D.C. area of 20 GHz and 5 deg elevation angle. The statistics of water vapor are shown to be sufficient for estimation of the statistics of gaseous absorption at 20 GHz. A computer model based on monthly surface weather is described and tested. Significant improvement in prediction of absorption extremes is demonstrated with the use of path weather data instead of surface data.

Russell, Thomas A.↗

Source Data Impacts on Epistemic Uncertainty for Launch Vehicle Fault Tree Models

Launch vehicle systems are designed and developed using both heritage and new hardware. Design modifications to the heritage hardware to fit new functional system requirements can impact the applicability of heritage reliability data. Risk estimates for newly designed systems must be developed from generic data sources such as commercially available reliability databases using reliability prediction methodologies, such as those addressed in MIL-HDBK-217F. Failure estimates must be converted from the generic environment to the specific operating environment of the system in which it is used. In addition, some qualification of applicability for the data source to the current system should be made. Characterizing data applicability under these circumstances is crucial to developing model estimations that support confident decisions on design changes and trade studies. This paper will demonstrate a data-source applicability classification method for suggesting epistemic component uncertainty to a target vehicle based on the source and operating environment of the originating data. The source applicability is determined using heuristic guidelines while translation of operating environments is accomplished by applying statistical methods to MIL-HDK-217F tables. The paper will provide one example for assigning environmental factors uncertainty when translating between operating environments for the microelectronic part-type components. The heuristic guidelines will be followed by uncertainty-importance routines to assess the need for more applicable data to reduce model uncertainty.

Al Hassan, Mohammad↗

Periodic Verification of the Scaling Factor for Radwastes in Korean NPPs - 13294

According to the acceptance criteria for a low and intermediate level radioactive waste (LILW) listed in Notice No. 2012-53 of the Nuclear Safety and Security Commission (NSSC), specific concentrations of radionuclides inside a drum has to be identified and quantified. In 5 years of effort, scaling factors were derived through destructive radiochemical analysis, and the dry active waste, spent resin, concentration bottom, spent filter, and sludge drums generated during 2004 ∼ 2008 were evaluated to identify radionuclide inventories. Eventually, only dry active waste among LILWs generated from Korean NPPs were first shipped to a permanent disposal facility on December 2010. For the LILWs generated after 2009, the radionuclides are being radiochemically quantified because the Notice clarifies that the certifications of the scaling factors should be verified biennially. During the operation of NPP, the radionuclides designated in the Notice are formed by neutron activation of primary coolant, reactor structural materials, corrosion products, and fission products released into primary coolant through defects or failures in fuel cladding. Eventually, since the radionuclides released into primary coolant are transported into the numerous auxiliary and support systems connected to primary system, the LILWs can be contaminated, and the radionuclides can have various concentration distributions. Thus, radioactive wastes, such as spent resin and dry active waste generated at various Korean NPP sites, were sampled at each site, and the activities of the regulated radionuclides present in the sample were determined using radiochemical methods. The scaling factors were driven on the basis of the activity ratios between a or β-emitting nuclides and γ-emitting nuclides. The resulting concentrations were directly compared with the established scaling factors' data using statistical methods. In conclusions, the established scaling factors were verified with a reliability of within 2σ, and the scaling factors will be applied for newly analyzed LILWs to evaluate the radionuclide inventories. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Stereotactic Body Radiation Therapy and the Influence of Chemotherapy on Overall Survival for Large (≥5 Centimeter) Non-Small Cell Lung Cancer

Purpose: Stereotactic body radiation therapy (SBRT) for ≥5 cm lesions is poorly defined, largely owing to the low sample sizes in existing studies. The present analysis examined the SBRT outcomes and assessed the effect of chemotherapy in this population. Methods and Materials: The National Cancer Data Base was queried for primary non-small cell lung cancer ≥5 cm treated with SBRT (≤10 fractions). Patient, tumor, and treatment parameters were extracted. The primary outcome was overall survival (OS). Statistical methods involved Kaplan-Meier analysis and multivariable Cox proportional hazards modeling. Results: From 2004 to 2012, data from 201 patients were analyzed. The median follow-up was 41.1 months. The median tumor size was 5.5 cm (interquartile range 5.0-6.0), with cT2a, cT2b, and cT3 disease in 24.9%, 53.2%, and 21.9%, respectively. The median total SBRT dose and fractionation was 50 Gy in 4 fractions, and 92.5% of the patients underwent SBRT with ≤5 fractions. The median OS was 25.1 months. Of the 201 patients, 15% received chemotherapy. The receipt of chemotherapy was associated with longer OS (median 30.6 vs 23.4 months; P=.027). On multivariable analysis, worse OS was seen with increasing age (hazard ratio [HR] 1.03; P=.012), poorly differentiated tumors (HR 2.06; P=.049), and T3 classification (HR 2.13; P=.005). On multivariable analysis, chemotherapy remained independently associated with improved OS (HR 0.57; P=.039). Conclusions: SBRT has utility in the setting of tumors ≥5 cm, with chemotherapy associated with improved OS in this subset. These hypothesis-generating data now raise the necessity of performing prospective analyses to determine whether chemotherapy confers outcome benefits after SBRT.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Statistical Approach To Extraction Of Texture In SAR

Improved statistical method of extraction of textural features in synthetic-aperture-radar (SAR) images takes account of effects of scheme used to sample raw SAR data, system noise, resolution of radar equipment, and speckle. Treatment of speckle incorporated into overall statistical treatment of speckle, system noise, and natural variations in texture. One computes speckle auto-correlation function from system transfer function that expresses effect of radar aperature and incorporates range and azimuth resolutions.

Rignot, Eric J.↗

Multimodel Predictive System for Carbon Dioxide Solubility in Saline Formation Waters

The prediction of carbon dioxide solubility in brine at conditions relevant to carbon sequestration (i.e., high temperature, pressure, and salt concentration (T-P-X)) is crucial when this technology is applied. Eleven mathematical models for predicting CO{sub 2} solubility in brine are compared and considered for inclusion in a multimodel predictive system. Model goodness of fit is evaluated over the temperature range 304–433 K, pressure range 74–500 bar, and salt concentration range 0–7 m (NaCl equivalent), using 173 published CO{sub 2} solubility measurements, particularly selected for those conditions. The performance of each model is assessed using various statistical methods, including the Akaike Information Criterion (AIC) and the Bayesian Information Criterion (BIC). Different models emerge as best fits for different subranges of the input conditions. A classification tree is generated using machine learning methods to predict the best-performing model under different T-P-X subranges, allowing development of a multimodel predictive system (MMoPS) that selects and applies the model expected to yield the most accurate CO{sub 2} solubility prediction. Statistical analysis of the MMoPS predictions, including a stratified 5-fold cross validation, shows that MMoPS outperforms each individual model and increases the overall accuracy of CO{sub 2} solubility prediction across the range of T-P-X conditions likely to be encountered in carbon sequestration applications.

54 ENVIRONMENTAL SCIENCES↗

Grid-based minimization at scale: Feldman-Cousins corrections for light sterile neutrino search

High Energy Physics (HEP) experiments generally employ sophisticated statistical methods to present results in searches of new physics. In the problem of searching for sterile neutrinos, likelihood ratio tests are applied to short-baseline neutrino oscillation experiments to construct confidence intervals for the parameters of interest. The test statistics of the form Δχ2 is often used to form the confidence intervals, however, this approach can lead to statistical inaccuracies due to the small signal rate in the region-of-interest. In this paper, we present a computational model for the computationally expensive Feldman-Cousins corrections to construct a statistically accurate confidence interval for neutrino oscillation analysis. The program performs a grid-based minimization over oscillation parameters and is written in C++. Our algorithms make use of vectorization through Eigen3, yielding a single-core speed-up of 350 compared to the original implementation, and achieve MPI data parallelism by employing DIY. We demonstrate the strong scaling of the application at High-Performance Computing (HPC) sites. We utilize HDF5 along with HighFive to write the results of the calculation to file.

Wospakrik, Marianette↗

Creating a Face from a DNA Sample

The idea of creating a recognizable face from a sample of DNA is moving closer to reality as Oak Ridge National Laboratory (ORNL) researchers employ a novel statistical representation of the face and a genome-wide association study (GWAS) to relate facial characteristics to genetic mutations. The reduction in cost of DNA sequencing and genotyping technology is leading to the availability of massive datasets of human genome sequences and human genotypes representing mutational variation among people. These large-scale data provide for the first time the potential to correlate an individual’s specific genetic code with their physical appearance. Using a specialized statistical method for representing the entire 3D shape of faces, ORNL researchers performed a GWAS to discover associations between genetic mutations and aspects of face shape. Data from 3,556 subjects resulted in the discovery of a very large set of statistically significant associations, which is strong evidence for the viability of predicting faces from DNA.

59 BASIC BIOLOGICAL SCIENCES↗

Active galactic nucleus X-ray variability in the XMM-COSMOS survey

We used the observations carried out by XMM in the COSMOS field over 3.5 yr to study the long term variability of a large sample of active galactic nuclei (AGNs) (638 sources) in a wide range of redshifts (0.1 < z < 3.5) and X-ray luminosities (10{sup 41} < L {sub 0.5-10} <10{sup 45.5}). Both a simple statistical method to assess the significance of variability and the Normalized Excess Variance (σ{sub rms}{sup 2}) parameter were used to obtain a quantitative measurement of the variability. Variability is found to be prevalent in most AGNs, whenever we have good statistics to measure it, and no significant differences between type 1 and type 2 AGNs were found. A flat (slope –0.23 ± 0.03) anti-correlation between σ{sub rms}{sup 2} and X-ray luminosity is found when all significantly variable sources are considered together. When divided into three redshift bins, the anti-correlation becomes stronger and evolving with z, with higher redshift AGNs being more variable. We prove, however, that this effect is due to the pre-selection of variable sources: when considering all of the sources with an available σ{sub rms}{sup 2} measurement, the evolution in redshift disappears. For the first time, we were also able to study long term X-ray variability as a function of M {sub BH} and Eddington ratio for a large sample of AGNs spanning a wide range of redshifts. An anti-correlation between σ{sub rms}{sup 2} and M {sub BH} is found, with the same slope of anti-correlation between σ{sub rms}{sup 2} and X-ray luminosity, suggesting that the latter may be a by-product of the former. No clear correlation is found between σ{sub rms}{sup 2} and the Eddington ratio in our sample. Finally, no correlation is found between the X-ray σ{sub rms}{sup 2} and optical variability.

79 ASTRONOMY AND ASTROPHYSICS↗

Quantitative methods and modeling to assess COVID–19–interrupted in vivo pharmacokinetic bioequivalence studies with two reference batches

The coronavirus disease 2019 (COVID-19) has presented unprecedented challenges to the generic drug development, including interruptions in bioequivalence (BE) studies. Per guidance published by the US Food and Drug Administration (FDA) during the COVID-19 public health emergency, any protocol changes or alternative statistical analysis plan for COVID-19-interrupted BE study should be accompanied with adequate justifications and not lead to biased equivalence determination. In this study, we used a modeling and simulation approach to assess the potential impact of study outcomes when two different batches of a Reference Standard (RS) were to be used in an in vivo pharmacokinetic BE study due to the RS expiration during the COVID-19 pandemic. Simulations were performed with hypothetical drugs under two scenarios: (1) uninterrupted study using a single batch of an RS, and (2) interrupted study using two batches of an RS. The acceptability of BE outcomes was evaluated by comparing the results obtained from interrupted studies with those from uninterrupted studies. The simulation results demonstrated that using a conventional statistical approach to evaluate BE for COVID-19-interrupted studies may be acceptable based on the pooled data from two batches. An alternative statistical method which includes a “batch” effect to the mixed effects model may be used when a significant “batch” effect was found in interrupted four-way crossover studies. However, such alternative method is not applicable for interrupted two-way crossover studies. Overall, the simulated scenarios are only for demonstration purpose, the acceptability of BE outcomes for the COVID19-interrupted studies could be case-specific.

60 APPLIED LIFE SCIENCES↗