Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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,405 records · Page 78
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Stochastic modeling and statistical calibration with model error and scarce data

This paper introduces a procedure to assess the predictive accuracy of stochastic models subject to model error and sparse data. Model error is introduced as uncertainty on the coefficients of appropriate polynomial chaos expansions (PCE). The error associated with finite sample size allows us to conceive of these coefficients as statistics of the data that we describe as random variables whose influence on output quantities of interest is evaluated through the extended polynomial chaos expansion (EPCE). A Bayesian data assimilation scheme is introduced to update these expansions by considering the resulting nested chaos expansion as a hierarchical probabilistic model. Stochastic models of quantities of interest (QoI) are thus constructed and efficiently evaluated. Here, the Metropolis–Hastings Markov chain Monte Carlo procedure is used to sample the posterior. Two illustrative analytical and numerical problems are used to demonstrate the proposed approach.

Bayesian inference↗
DOE OSTI2024

Stochastic symplectic reduced-order modeling for model-form uncertainty quantification in molecular dynamics simulations in various statistical ensembles

Here, this work focuses on the representation of model-form uncertainties in molecular dynamics simulations in various statistical ensembles. In prior contributions, the modeling of such uncertainties was formalized and applied to quantify the impact of, and the error generated by, pair-potential selection in the microcanonical ensemble (NVE). In this work, we extend this formulation and present a linear-subspace reduced-order model for the canonical (NVT) and isobaric (NPT) ensembles. The symplectic reduced-order basis is randomized on the tangent space of the Stiefel manifold to provide topological relationships and capture model-form uncertainty. Using the Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS), we assess the relevance of these stochastic reduced-order atomistic models on canonical problems involving a Lennard-Jones fluid and an argon crystal melt.

42 ENGINEERING↗
DOE OSTI2021

Kinetic, metabolic, and statistical analytics: addressing metabolic transport limitations among organelles and microbial communities

Microbial organisms engage in a variety of metabolic interactions. A crucial part of these interactions is the exchange of molecules between different organelles, cells, and the environment. The main forces mediating this metabolic exchange are transporters. This transport can be difficult to measure experimentally because several transport mechanisms remain opaque. However, theoretical calculations about the inputs and outputs of cells via metabolic exchanges have enabled the successful inference of the workings of intra-organismal and inter-organismal systems. Kinetic, metabolic, and statistical modeling approaches in combination with omics data are enhancing our knowledge and understanding about metabolic exchange and mass resource allocation. Furthermore, this model-driven analytics approach can guide effective experimental design and yield new insights into biological function and control.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2016

GPU acceleration of the Locally Selfconsistent Multiple Scattering code for first principles calculation of the ground state and statistical physics of materials

The Locally Self-consistent Multiple Scattering (LSMS) code solves the first principles Density Functional theory Kohn–Sham equation for a wide range of materials with a special focus on metals, alloys and metallic nano-structures. It has traditionally exhibited near perfect scalability on massively parallel high performance computer architectures. In this paper, we present our efforts to exploit GPUs to accelerate the LSMS code to enable first principles calculations of O(100,000) atoms and statistical physics sampling of finite temperature properties. We reimplement the scattering matrix calculation for GPUs with a block matrix inversion algorithm that only uses accelerator memory. Finally, using the Cray XK7 system Titan at the Oak Ridge Leadership Computing Facility we achieve a sustained performance of 14.5PFlop/s and a speedup of 8.6 compared to the CPU only code.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Deconvolving distribution of relaxation times, resistances and inductance from electrochemical impedance spectroscopy via statistical model selection: Exploiting structural-sparsity regularization and data-driven parameter tuning

The distribution of relaxation times (DRT) has drawn increasing attention for interpreting electrochemical impedance spectroscopy (EIS). Deconvolution of DRT from EIS is a challenging ill-posed problem that requires regularization methods. In this work, we formulate DRT reconstruction task as a statistical model selection problem with structural-sparsity penalties. We utilize the Elastic net regularization that simultaneously benefits from Ridge and Lasso regularizations with optimal tuning parameter automatically determined by the information criteria. We benchmark our approach on four synthetic experiments (a ZARC element, ZARC mixtures, a RC circuit and a Fractal element) and two real EIS datasets of a Lithium ion battery and an organic-inorganic halide class of perovskites in oxygen environment at different gas pressures. We demonstrate the superiority of proposed model selection procedure, that is capable of eliminating pseudo peaks and representing asymmetries in DRT as well as precisely estimating resistances. Here, we highlight our approach is robust to reducing and subsampling EIS frequency range, making it a promising tool for timing-resolved, localized and large scale EIS data analysis. For the Lithium ion battery data analysis, we extend the classical DRT model to incorporate the inductive effect and illustrate DRT as a guidance for equivalent circuit modeling to refine impedance reconstruction at low risks of overfitting. Furthermore, the structural-sparsity regularization could be extended for multidimensional and Bayesian EIS data analysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Understanding and improving anode performance in an alkaline membrane electrolyzer using statistical design of experiments

The behavior of the oxygen-evolving positive electrode (i.e. anode) in the anion exchange membrane water electrolyzer (AEMEL) is complex and influenced by several factors. Very few studies have been performed to understand oxygen evolution reaction (OER) behavior by optimizing the individual factors that influence performance. Here, this study highlights the effects of catalyst loading, catalyst selection, porous transport layer (PTL) type and conductive additive content. The influence of each factor is elucidated through a design of experiments (DoE) approach with a full statistical analysis. Electrochemical data, alongside Pareto charts, parametric trends and their mutual interactions are discussed. This DoE approach is also helpful in making useful predictions and discovering new combinations to be tested. The end result was a high-performance AEMEL able to operate at a current density of 1.0 A/cm 2 at 1.80 V with IrOx OER and PtNi hydrogen evolution reaction (HER) catalysts using 0.3 M KOH fed to the anode. Even lower operating voltage was observed with PbRuOx catalyst at the anode, 1.64 V @ 1.0 A/cm 2 , though the cell decay rate was higher. Lastly, a IrOx/PtNi cell was stably operated continuously for 30 days (720 h) at 1.0 A/cm 2 . This study can serve as a guide for optimal electrode design with insights into component-performance compromises, which can aid in making design choices and performing techno-economic analyses.

08 HYDROGEN↗
DOE OSTI2014

Development and application of a statistical methodology to evaluate the predictive accuracy of building energy baseline models

This paper documents the development and application of a general statistical methodology to assess the accuracy of baseline energy models, focusing on its application to Measurement and Verification (M&V) of whole-­building energy savings. The methodology complements the principles addressed in resources such as ASHRAE Guideline 14 and the International Performance Measurement and Verification Protocol. It requires fitting a baseline model to data from a ``training period’’ and using the model to predict total electricity consumption during a subsequent ``prediction period.’’ We illustrate the methodology by evaluating five baseline models using data from 29 buildings. The training period and prediction period were varied, and model predictions of daily, weekly, and monthly energy consumption were compared to meter data to determine model accuracy. Several metrics were used to characterize the accuracy of the predictions, and in some cases the best-­performing model as judged by one metric was not the best performer when judged by another metric.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Nanotomographic observation and statistical analysis of overcharging induced cracks in LiCoO 2 single crystalline particles

Particle fracture represents an important mode of battery electrode degradation. Extensive cracking occurs in the layered oxide cathodes of lithium-ion rechargeable batteries when they are charged beyond safe operational voltages. In this paper, we analyze the morphology and distribution of cracks in over 100 overcharged LiCoO 2 (LCO) single-crystalline particles imaged via X-ray nanotomography. Herein, we observe a strong particle size dependence in both nucleation and growth of cracks, with percentage of cracked particles, crack size, and crack surface area per unit volume all increasing significantly with particle volume. Increasing the cutoff voltage above 4.6V is found to trigger much more pronounced crack growth, which we attribute to the large misfit stress caused by the H1-3 to O1 phase transition. This study offers valuable statistical insights into the fracture behavior of layered oxide cathodes and reveals the importance of particle size control on electromechanical degradation of battery electrodes.

25 ENERGY STORAGE↗
DOE OSTI2023

Statistical chronometry of Meteorites: II. Initial abundances and homogeneity of short-lived radionuclides

Astrophysical models of planet formation require accurate radiometric dating of meteoritic components by short-lived (Al-Mg, Mn-Cr, Hf-W) and long-lived (U-Pb) chronometers, to develop a timeline of such events in the solar nebula as formation of Ca-rich, Al-rich Inclusions (CAIs), chondrules, planetesimals, etc. CAIs formed mostly around a time (“t=0”) when the short-lived radionuclide 26 Al (t 1/2 = 0.72 Myr) was present and presumably homogeneously distributed at a known level we define as ( 26 Al/ 27 Al) SS = 5.23 × 10 -5 . The time of formation after t=0 of another object can be found by determining its initial ( 26 Al/ 27 Al) 0 ratio and comparing it to ( 26 Al/ 27 Al) SS . Dating of meteoritic objects using the Mn-Cr or Hf-W systems is hindered because the abundances ( 53 Mn/ 55 Mn) SS and ( 182 Hf/ 180 Hf) SS at t=0 are not known precisely. To constrain these quantities, we compile literature Al-Mg, Mn-Cr, Hf-W and Pb-Pb data for 13 achondrites and use novel statistical techniques to minimize the discrepancies between their times of formation across these systems. We find that for ( 53 Mn/ 55 Mn) SS = (7.80±0.36)×10 -6 , ( 182 Hf/ 180 Hf) SS = (10.41 ± 0.12) × 10-5, t SS = 4568.65 ± 0.10 Myr, and a 53 Mn half-life of 3.98 ± 0.22 Myr, these four free parameters make concordant 18 formation times recorded by the different systems in all six known volcanic achondrites (the quenched angrites D’Orbigny, Sahara 99555, Northwest Africa 1670; the psuedo-eucrites Asuka 881394 and Ibitira; and the basaltic achondrite Northwest Africa 7325). These parameters also make concordant the ages derived for chondrules from CB/CH achondrites, formed simultaneously in an impact. The other seven achondrites are not quite concordant, but are plutonic angrites or ‘carbonaceous achondrites’ for which simultaneous closure of the isotopic systems might not be expected. Our findings provide very strong support for homogeneity of 26 Al, 53 Mn, and 182 Hf in the solar nebula, and our approach offers a path for more precise chronometry

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Statistical Modeling Approach to Quantitative Analysis of Interobserver Variability in Breast Contouring

Purpose: To develop a new approach for interobserver variability analysis. Methods and Materials: Eight radiation oncologists specializing in breast cancer radiation therapy delineated a patient's left breast “from scratch” and from a template that was generated using deformable image registration. Three of the radiation oncologists had previously received training in Radiation Therapy Oncology Group consensus contouring for breast cancer atlas. The simultaneous truth and performance level estimation algorithm was applied to the 8 contours delineated “from scratch” to produce a group consensus contour. Individual Jaccard scores were fitted to a beta distribution model. We also applied this analysis to 2 or more patients, which were contoured by 9 breast radiation oncologists from 8 institutions. Results: The beta distribution model had a mean of 86.2%, standard deviation (SD) of ±5.9%, a skewness of −0.7, and excess kurtosis of 0.55, exemplifying broad interobserver variability. The 3 RTOG-trained physicians had higher agreement scores than average, indicating that their contours were close to the group consensus contour. One physician had high sensitivity but lower specificity than the others, which implies that this physician tended to contour a structure larger than those of the others. Two other physicians had low sensitivity but specificity similar to the others, which implies that they tended to contour a structure smaller than the others. With this information, they could adjust their contouring practice to be more consistent with others if desired. When contouring from the template, the beta distribution model had a mean of 92.3%, SD ± 3.4%, skewness of −0.79, and excess kurtosis of 0.83, which indicated a much better consistency among individual contours. Similar results were obtained for the analysis of 2 additional patients. Conclusions: The proposed statistical approach was able to measure interobserver variability quantitatively and to identify individuals who tended to contour differently from the others. The information could be useful as feedback to improve contouring consistency.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2018

Enhancing statistical moment calculations for stochastic Galerkin solutions with Monte Carlo techniques

Here, we provide a method for enhancing stochastic Galerkin moment calculations to the linear elliptic equation with random diffusivity using an ensemble of Monte Carlo solutions. This hybrid approach combines the accuracy of low-order stochastic Galerkin and the computational efficiency of Monte Carlo methods to provide statistical moment estimates which are significantly more accurate than performing each method individually. The hybrid approach involves computing a low-order stochastic Galerkin solution, after which Monte Carlo techniques are used to estimate the residual. We show that the combined stochastic Galerkin solution and residual is superior in both time and accuracy for a one-dimensional test problem and a more computational intensive two-dimensional linear elliptic problem for both the mean and variance quantities.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

A probabilistic drought forecasting framework: A combined dynamical and statistical approach

In order to improve drought forecasting skill, this study develops a probabilistic drought forecasting framework comprised of dynamical and statistical modeling components. The novelty of this study is to seek the use of data assimilation to quantify initial condition uncertainty with the Monte Carlo ensemble members, rather than relying entirely on the hydrologic model or land surface model to generate a single deterministic initial condition, as currently implemented in the operational drought forecasting systems. Next, the initial condition uncertainty is quantified through data assimilation and coupled with a newly developed probabilistic drought forecasting model using a copula function. The initial condition at each forecast start date are sampled from the data assimilation ensembles for forecast initialization. Finally, seasonal drought forecasting products are generated with the updated initial conditions. This study introduces the theory behind the proposed drought forecasting system, with an application in Columbia River Basin, Pacific Northwest, United States. Results from both synthetic and real case studies suggest that the proposed drought forecasting system significantly improves the seasonal drought forecasting skills and can facilitate the state drought preparation and declaration, at least three months before the official state drought declaration.

data assimilation↗
DOE OSTI2026

Capturing the fractocohesive length scale in elastomers through a statistical mechanics-based gradient enhanced damage model

Here, this study aims to examine modeling flaw sensitivity in elastomers. The direct incorporation of polymer chain statistical mechanics considerations into a continuum stretch-based gradient-enhanced damage formulation, in turn, allows a representation of diffuse chain damage and fracture events that align with known micromechanical mechanisms. Through a series of numerical experiments, we simulate crack propagation and extract the fracture energy as an output of the model, while keeping track of the micromechanical signatures of diffuse chain damage that accommodate fracture propagation and eventually influence flaw sensitivity. Finally, by combining the fracture toughness and the work to rupture, we identify a fractocohesive length of the material, corresponding to the full width of the damage process zone. As the damage-to-fracture cascade in the proposed GED model is influenced by the introduction of a length scale associated with network imperfection and long-range load transfer, the emerging relationship of the two length scales is discussed, providing a potential link between microscopic damage mechanisms and the observed macroscopic fracture response.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Application of the immobilized low-activity waste glass corrosion model to the static dissolution of 24 statistically-designed alkali-borosilicate waste glasses

Glass corrosion models that capture the complex mechanisms of the glass-water reaction enable the prediction of nuclear waste glass durability in disposal scenarios. Parameterization of such models is challenging because of the need to capture changes in corrosion behavior with time, reaction conditions, and glass composition. In this report we describe and employ the immobilized low activity waste (ILAW) glass corrosion model (IGCM) in geochemical simulations of static dissolution tests, at two temperatures (40 °C and 90 °C), for a matrix of 24 enhanced low activity waste (eLAW) glasses statistically designed to cover a processable composition space defined by 8 major glass components (Al 2 O 3 , B 2 O 3 , CaO, Na 2 O, SiO 2 , SnO 2 , ZrO 2 , and Others as defined in the text). The IGCM includes a first-order chemical affinity term and an ion-exchange term that represents the net exchange of Na+ ions in the pristine glass with protons in aqueous solution. Constant, time-dependent, and time- and pH-dependent functional forms of the ion-exchange term are evaluated to reproduce the change in corrosion behavior with time in saturated, static dissolution tests. The agreement with measured aqueous concentrations of the main glass components (B, Na, and Si) improved significantly upon addition of a time-dependent term, which therefore constitutes a simple representation of the glass-water reaction progress. Due to the limited changes in pH in the static dissolution tests, past a short initial period of rapid increase, addition of a pH dependent term did not appreciably improve the fits, indicating that comprehensive model parameterization requires more than one type of glass corrosion test to capture a wide range of solution chemistries. IGCM parameters were found to be dependent on glass composition, and the parameter sets generated in this work will enable the development of composition–parameter correlation models that offer the promise of predicting IGCM parameters, and thus glass corrosion behavior, solely based on glass composition.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2017

Failure statistics for commercial lithium ion batteries: A study of 24 pouch cells

There are relatively few publications that assess capacity decline in enough commercial cells to quantify cell-to-cell variation, but those that do show a surprisingly wide variability. Capacity curves cross each other often, a challenge for efforts to measure the state of health and predict the remaining useful life (RUL) of individual cells. We analyze capacity fade statistics for 24 commercial pouch cells, providing an estimate for the time to 5% failure. Our data indicate that RUL predictions based on remaining capacity or internal resistance are accurate only once the cells have already sorted themselves into “better” and “worse” ones. Analysis of our failure data, using maximum likelihood techniques, provide uniformly good fits for a variety of definitions of failure with normal and with 2- and 3-parameter Weibull probability density functions, but we argue against using a 3-parameter Weibull function for our data. pdf fitting parameters appear to converge after about 15 failures, although business objectives should ultimately determine whether data from a given number of batteries provides sufficient confidence to end lifecycle testing. Increased efforts to make batteries with more consistent lifetimes should lead to improvements in battery cost and safety.

25 ENERGY STORAGE↗
DOE OSTI2019

A Bayesian statistical approach to inferring particle dynamics from in-situ pump POC and chloropigment data from the Mediterranean Sea

Concentrations of chloropigments and particulate organic carbon (POC) in large-volume in-situ pump samples from the Mediterranean Sea were used to estimate rate constants of processes that control the fate of particles, and specifically chloropigments, in the water column. In this paper we introduce a Bayesian statistical inversion method that combines the data with a new box model and has the capacity to infer rate constants for POC respiration/dissolution, chlorophyll and pheopigment degradation, and particle aggregation and disaggregation. We use first-order kinetics to model disaggregation, and use both first-order and second-order kinetics to model aggregation. Using these methods, the estimated small-particle (1–70 μm) POC respiration rate constant was 2.44 -1.00 +1.69 yr -1 (0.41 yr). The estimated disaggregation and second-order aggregation rate constants were 85.6 -36.4 +63.4 yr -1 (1.17×10 -2 yr) and 2.78 -1.17 +2.01 μM -1 yr -1 , respectively. Using the optimal rate constants and the corresponding particle concentrations, disaggregation is ~4.2 times faster than the small-size POC dissolution rate, which indicates that disaggregation is a dominant processes at the time of sampling. More importantly, by comparing our results with those of previous studies, we conclude that sampling methods have less influence than tracers themselves on inferring particle dynamic rate constants. We previously introduced a somewhat similar approach to modeling SV sediment trap data, but large volume pumps are a much more common sample collection method in oceanographic surveys than SV sediment traps, and thus our new model should have a wider applicability.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2015

Segmentation, statistical analysis, and modelling of the wall system in ceramic foams

Closed walls in otherwise open foam structures may have a great impact on macroscopic properties of the materials. In this paper, we present two algorithms for the segmentation of such closed walls from micro-computed tomography images of the foam structure. The techniques are compared on simulated data and applied to tomographic images of ceramic filters. This allows for a detailed statistical analysis of the normal directions and sizes of the walls. Finally, we explain how the information derived from the segmented wall system can be included in a stochastic microstructure model for the foam.

36 MATERIALS SCIENCE↗