Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “uncertainty Quantification”

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,531 records · Page 85
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Quantifying Uncertainties In Neutron-Alpha Scattering With Chiral Nucleon-Nucleon And Three-Nucleon Forces

We report that modern ab initio theory combined with high-quality nucleon-nucleon (NN) and three-nucleon (3N) interactions from chiral effective field theory (EFT) can provide a predictive description of low-energy light-nuclei reactions relevant for astrophysics and fusion-energy applications. However, the high cost of computations has so far impeded a complete analysis of the uncertainty budget of such calculations. Starting from NN potentials up to fifth order (N 4 LO) combined with leading-order 3N forces, we study how the order-by-order convergence of the chiral expansion and confidence intervals for the 3N contact and contact-plus-one-pion-exchange low-energy constants (c E and c D ) contribute to the overall uncertainty budget of many-body calculations of neutron- 4 He (n-α) elastic scattering. We compute structure and reaction observables for three-, four-, and five-nucleon systems within the ab initio frameworks of the no-core shell model and no-core shell model with continuum. Using a small set of design runs, we construct a Gaussian process model (GPM) that acts as a statistical emulator for the theory. With this, we gain insight into how uncertainties in the 3N low-energy constants propagate throughout the calculation and determine the Bayesian posterior distribution of these parameters with Markov-Chain Monte Carlo. We find rapidly converging n-α phase shifts with respect to the chiral order. With the adopted leading-order 3N force, calculations based on the NN interaction at N 4 LO of Entem, Machleidt, and Nosyk are unable to reproduce the experimental phase shifts in the 3/2 - channel within the estimated chiral truncation errors. Closer agreement with empirical data is found when using an older parametrization of the NN interaction at order N 3 LO, and the position and width of the P-wave resonances can be used to reduce the uncertainty of the 3N low-energy constants. The present results point to a lack of spin-orbit strength when the newer parametrization of the chiral NN force up to fifth order is combined with the leading-order 3N force. The inclusion of higher-order 3N-force terms may be required to recover the missing strength. GPMs can act as fast and accurate emulators of ab initio many-body calculations of low-energy scattering and reactions of light nuclei, opening the way to a robust quantification of theoretical uncertainties grounded in the description of the underlying chiral Hamiltonian.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

High-Accuracy Heats of Formation for Alkane Oxidation: From Small to Large via the Automated CBH-ANL Method

It is generally challenging to obtain high-accuracy predictions for the heat of formation for species with more than a handful of heavy atoms, such as those of importance in standard combustion mechanisms. To this end, we construct the CBH-ANL approach and illustrate that, for a set of 194 alkane oxidation species, it can be used to produce ΔH f (o K) values with 2σ uncertainties of 0.2-0.5 kcal mol -1 . This set includes the alkanes, hydroperoxides, and alkyl, peroxy, and hydroperoxyalkyl radicals for 17 representative hydrocarbon fuels containing up to 10 heavy atoms with various degrees of branching in the alkane backbone. The CBH-ANL approach, automated in the QTC and AutoMech software suites, builds balanced chemical equations for the calculation of ΔH f (o K), in which the reference species may be up to five heavy atoms. The high-level ANLo and ANL1 reference ΔH f (o K) values are further refined for even the largest of these reference species with a novel laddering approach. We perform a comprehensive quantification of the uncertainties for both the individual reference species (the largest of which is 0.15 kcal mol -1 ) and the propagation of those uncertainties when used in the calculation of ΔH f (o K) for the 194 target species. Here we examine the sensitivity of the predicted ΔH f (o K) values to (i) electronic energies from various methods, including ωB97X-D/cc-pVTZ, B2PLYP-D3/cc-pVTZ, CCSD(T)-F12b/cc-pVDZ-F12//B2PLYP-D3/cc-pVTZ, and CCSD(T)-F12b/cc-pVTZ-F12//B2PLYP-D3/cc-pVTZ; (ii) the zero-point vibrational energies (ZPVEs), where we consider harmonic ZPVEs as well as two scaling-based estimates of the anharmonic ZPVEs, all implemented for both ωB97X-D/cc-pVTZ and B2PLYP-D3/cc-pVTZ calculations; (iii) the particular CBH-ANL scheme employed; and (iv) the procedure for choosing the reference conformer for the analyses. The discussion concludes with a summary of the estimated overall uncertainty in the predictions and a validation of the predictions for the alkane subset.

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

Quantification of oil and gas methane emissions in the Delaware and Marcellus basins using a network of continuous tower-based measurements

Abstract. According to the United States Environmental Protection Agency (US EPA), emissions from oil and gas infrastructure contribute 30 % of all anthropogenic methane (CH4) emissions in the US. Studies in the last decade have shown emissions from this sector to be substantially larger than bottom-up assessments, including the EPA inventory, highlighting both the increased importance of methane emissions from the oil and gas sector in terms of their overall climatological impact and the need for independent monitoring of these emissions. In this study we present continuous monitoring of regional methane emissions from two oil and gas basins using tower-based observing networks. Continuous methane measurements were taken at four tower sites in the northeastern Marcellus basin from May 2015 through December 2016 and five tower sites in the Delaware basin in the western Permian from March 2020 through April 2022. These measurements, an atmospheric transport model, and prior emission fields are combined using an atmospheric inversion to estimate monthly methane emissions in the two regions. This study finds the mean overall emission rate from the Delaware basin during the measurement period to be 146–210 Mg CH4 h−1 (energy-normalized loss rate of 1.1 %–1.5 %, gas-normalized rate of 2.5 %–3.5 %). Strong temporal variability in the emissions was present, with the lowest emission rates occurring during the onset of the COVID-19 pandemic. Additionally, a synthetic model–data experiment performed using the Delaware tower network shows that the presence of intermittent sources is not a significant source of uncertainty in monthly quantification of the mean emission rate. In the Marcellus, this study finds the overall mean emission rate to be 19–28 Mg CH4 h−1 (gas-normalized loss rate of 0.30 %–0.45 %), with relative consistency in the emission rate over time. These totals align with aircraft top-down estimates from the same time periods. In both basins, the tower network was able to constrain monthly flux estimates within ±20 % uncertainty in the Delaware and ±24 % uncertainty in the Marcellus. The results from this study demonstrate the ability to monitor emissions continuously and detect changes in the emissions field, even in a basin with relatively low emissions and complex background conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Aerosol Hygroscopic Growth, Mixing State, and Cloud Condensation Nuclei Activity during TRACER (Field Campaign Report)

Convective clouds play a critical role in the Earth’s climate system. Recent research has shown that a realistic representation of convective processes is critical to constraining climate sensitivity in global climate models. Theoretical and modeling studies showed that aerosols could have strong dynamic feedback to convection in warm and humid environments through enhancing ice-related processes and condensational growth. A few observation-based studies also suggested the influence of aerosols on convective cloud and precipitation properties. However, robust observational quantification of an aerosol effect on convective clouds isolated from other factors remains elusive. Understanding the impact of aerosol on convective clouds requires knowledge of the cloud condensation nuclei (CCN) spectrum, which represents the number of particles that uptake water and form cloud droplets as a function of supersaturation. The water uptake by aerosol is also of critical importance for the direct interaction of aerosol with radiation (i.e., aerosol direct effect) due to light scattering and absorption by aerosol. While both droplet activation under supersaturated conditions (i.e., relative humidity RH > 100%) and the hygroscopic growth under sub-saturated conditions (i.e., RH < 100%) are strongly influenced by particle hygroscopicity, the thermodynamic regimes and measurement methods are quite different. Aerosol particles, especially organic particles, can exhibit higher hygroscopicity for droplet activation than that for hygroscopic growth. In the subsaturated regime, the hygroscopicity of organic particles can also vary strongly with RH. However, global climate models usually treat organic species in aerosols with a constant hygroscopicity, potentially introducing substantial uncertainties in the quantification of aerosol radiative effects.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Redefining Net-Zero Water for Resilient Aquifers

Groundwater has long been a critical source of potable water, supplying one-third of all global water withdrawals and providing strategic reserves in times of drought (Konikow, 2012; Famiglietti, 2014). The risk in relying on groundwater is that it is often pumped at rates far greater than it can be naturally replenished. With increasing population and a changing climate, the depletion of potable water supplies has become an ever-growing concern (Srinivasan et al., 2012). Recently, the concept of “resilience” has been applied to groundwater, leaving many questions to be explored in terms of what resilience should mean for aquifers (Foster and MacDonald, 2014; Grönwall and Oduro-Kwarteng, 2017). The term is commonly used to indicate resilience to impacts of climate change specifically, but here we consider resilience in the context of more general water insecurity. Relatedly, the concept of “sustainable yield” has emerged in an attempt to determine a metric that can ensure the long-term resilience of groundwater systems, but the definition is fraught with inherent uncertainties and its quantification poses challenges (Rudestam and Langridge, 2014; Grönwall and Oduro-Kwarteng, 2017). At the heart of these concepts is a focus on sustainable management of aquifers. A sustainable approach to managing groundwater resources usually involves balancing long-term abstraction with the long-term recharge taking into account the needs for environmental flows and implementing drought-preparedness measures (Robins et al. 2006; Foster and MacDonald, 2014). At a high level, a sustainable balance of groundwater abstraction versus aquifer recharge may be considered “Net Zero” (NZ) or “Net Positive” (NP) in cases where recharge is equal to or in excess of abstraction, respectively. However, NZ/NP water management programs that address these concerns have varied definitions of compliance, which include considerations for the application at hand (ILFI, 2020; USGBC, 2020). Locations with high risk of low or poor water supply may be particularly interested in pursuing NZ/NP water. The AQUEDUCT Water Risk Atlas (World Resources Institute, 2019) indicates long-term availability, providing risk levels for quantity and quality of water at specific locations across the globe, for existing and predicted future conditions. This paper follows the U.S. Department of Energy (DOE) Federal Energy Management Program (FEMP) NZ water definition adopted by the U.S. Army for use at their installations, including Army Reserve Centers (ARCs). This definition also relates to the U.S. Army’s current focus on water resilience, which includes short-term and long-term availability of water resources.

net zero water, HUC12, Saipan, American Samoa, Mau↗
DOE OSTI2024

Rapid assessment of the creep rupture life of metals: A model enabling experimental design

Prediction of the creep rupture life of engineering metals is critical for qualification and design of new materials. The use of long-term creep tests and the need to quantify the performance variability in a priori similar systems hinder the rapid creep assessment of a given material. Therefore, it is essential to develop methods that can extrapolate the long-term performance of alloys and the associated variability from short-term experiments. To this end, this study introduces a new model which enables the estimation of the rupture life of a material for a given stress and temperature. This model relies on two components. First, a new relation for the minimum creep rate (MCR) of materials is introduced. It includes a stress dependent stress exponent allowing the model to capture the variation of MCR across a wide range of temperatures and stresses. Second, employing the Monkman-Grant (MG) law, we establish a relation between stress, temperature and creep rupture life. Together, these two elements yield a new closed-form mathematical expression for the Larson Miller parameter as a function of stress and temperature. This expression captures the creep rupture time for many metals (Gr91, Copper, Gr122 and 347H) and compares favorably with alternate empirical approaches. The model is then used to assess the minimum duration of creep rates necessary to qualify the material up to 100000h. Furthermore, it is found that depending on the material system, creep tests as few as five limited to 5000 h for steels (Gr91, Gr122, 347H) and 100 h for copper are sufficient to model creep lifetimes. Finally, using a Bayesian inference-based approach to calibrate the model, we demonstrate that variability in rupture life can be captured via the quantification of the uncertainty in the model parameters and extrapolated from a limited number of short to moderately short creep tests; thereby paving the way for accelerated creep testing.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Measurement of the 239 Pu(n,f)/ 235 U(n,f) Cross-Section Ratio with the NIFFTE fission Time Projection Chamber

In this work, the 239 Pu(n,f)/ 235 U(n,f) cross-section ratio has been measured with the fission Time Projection Chamber (fissionTPC) from 100 keV to 100 MeV. The fissionTPC provides three-dimensional reconstruction of fission-fragment ionization profiles, allowing for a precise quantification of measurement uncertainties. The measurement was performed at the Los Alamos Neutron Science Center which provides a pulsed white source of neutrons. The data are recommended to be used as a cross-section ratio shape. A discussion of the status of the absolute normalization and comparisons to ENDF evaluations and previous measurements is included.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Quantifying heavy quark transport coefficients with an improved transport model

We report that the heavy-flavor transport coefficients contain important information on the strong interaction at finite temperatures. The extraction of these numbers from experimental data requires dynamical modeling of heavy-flavor transport that is coupled to realistic medium evolution. Furthermore, meaningful extractions necessitate both a faithful implementation of the physical inputs to be tested and the quantification of model uncertainty. For these purposes, we have developed a partonic transport model LIDO. It has an improved treatment of in-medium parton bremsstrahlung, which has been calibrated to theoretical calculations in a simple medium to reduce modeling uncertainty. Regarding the interaction between heavy quark and the medium, few-body perturbative scatterings are applied to large-momentum transfer (q) processes, while a diffusion equation models the dynamics of small-q processes. Such a separation restricts the explicit use of medium quasi-particles to large-q processes only. Another advantage is that deviations from the leading-order probe-medium coupling can be parametrized as an additional contribution to the diffusion constant. The heavy quark transport coefficients are then extracted with uncertainty estimation from a Bayesian analysis including both the RHIC and the LHC data. The results are found to be consistent with earlier extraction of the light-quark transport coefficients at high momentum and be comparable with lattice calculations of the heavy-flavor diffusion constant in the static limit at low momentum.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Enhanced mid-latitude warming due to poleward shift of sea salt aerosols in a warmer future

Natural aerosols provide large uncertainties regarding the quantification of the climate effects of aerosols, especially those related to aerosol-cloud interactions. Advancing research on the interaction of natural aerosols and climate can provide a better assessment of the climate impacts of anthropogenic activities. Furthermore, as one of the most dominant natural aerosols, sea salt has the highest proportion of mass in natural aerosols. Sea salt can be regarded as a purely scattering aerosol, causing a negative radiative effect by scattering shortwave solar radiation and acting as cloud condensation nuclei (CCN) through aerosol-cloud interactions. Future climate change will influence sea salt distribution, thus modifying the radiative effects of sea salt, which will further affect the climate. However, few previous studies have explored the interactions between sea salt changes under warming scenarios and climate change.

Yang, Yang [Nanjing Univ. of Information Science a↗
DOE OSTI2020

Fast strength model characterization using Bayesian statistics

A variety of flow stress models exist with new models constantly being developed. These models aim to approximate the strength of materials in a variety of regimes from quasistatic loading through shock scenarios. All models contain an array of parameters which need to be tuned to the material under study. Some models perform well under limited conditions, requiring adjustment of the parameters when venturing outside of those predefined ranges. Other models perform well over a wide range of conditions with a set of parameters, but may be outperformed by other models optimized on a tighter range of conditions. Recent research by Los Alamos demonstrated the ability to optimize the Johnson Cook (JC) model using a set of 3 plate-impact experiments on Aluminum. They utilized Bayesian statistics and emulation to determine optimal parameters for the model with a quantification of parameter uncertainty. We present an extension of this capability to incorporate velocimetry from plate-impact tests, stress-strain data from split Hopkinson pressure bar and quasistatic compression tests, plus profiles from Taylor cylinders in a unified fashion. Statistically robust comparisons of the performance and uncertainty of different realizations of the JC flow stress model were carried out based on calibration to several possible combinations of these three different experiment types.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Assessing the numerical stability of physics models to equilibrium variation through database comparisons on DIII-D

High fidelity kinetic equilibria are crucial for tokamak modeling and analysis. Manual workflows for constructing kinetic equilibria are time consuming and subject to user error, motivating development of automated equilibrium reconstruction tools to provide accurate and consistent reconstructions for downstream physics analysis. These automated tools also provide access to kinetic equilibria at large database scales, which enables the quantification of general uncertainties arising from equilibrium reconstruction techniques. In this paper, we compare a large database of DIII-D kinetic equilibria generated manually by physics experts to equilibria from automated kinetic reconstruction tools, assessing the impact of reconstruction method on equilibrium parameters and resulting magnetohydrodynamic stability calculations. We find agreement among scalar parameters, whereas profile quantities, such as the bootstrap current, show larger disagreements. We analyze ideal kink and classical tearing stability with DCON and STRIDE respectively, finding that the kink stability calculation is generally more robust than the tearing index Δ' calculation. We find that in 90% of cases, both kink stability classifications are unchanged between the manual expert and automated kinetic equilibria.

CAKE↗
DOE OSTI2024

Assessing the design of integrated methane sensing networks

Abstract While methane is the second largest contributor to global warming after carbon dioxide, it has a larger warming effect over a much shorter lifetime. Despite accelerated technological efforts to radically reduce global carbon dioxide emissions, rapid reductions in methane emissions are needed to limit near-term warming. Being primarily emitted as a byproduct from agricultural activities and energy extraction, methane is currently monitored via bottom–up (i.e. activity level) or top–down (via airborne or satellite retrievals) approaches. However, significant methane leaks remain undetected and emission rates are challenging to characterize with current monitoring frameworks. In this paper, we study the design of a layered monitoring approach that combines bottom–up and top–down approaches as an integrated sensing network. By recognizing that varying meteorological conditions and emission rates impact the efficacy of bottom–up monitoring, we develop a probabilistic approach to optimal sensor placement in its bottom–up network. Subsequently, we derive an inverse Bayesian framework to quantify the improvement that a design-optimized integrated framework has on emission-rate quantifications and their uncertainties. We find that under realistic meteorological conditions, the overall error in estimating the true emission rates is approximately 1.3 times higher, with their uncertainties being approximately 2.4 times higher, when using a randomized network over an optimized network, highlighting the importance of optimizing the design of integrated methane sensing networks. Further, we find that optimized networks can improve scenario coverage fractions by more than a factor of 2 over experimentally-studied networks, and identify a budget threshold beyond which the rate of optimized-network coverage improvement exhibits diminishing returns, suggesting that strategic sensor placement is also crucial for maximizing network efficiency.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Implementation of a practical Markov chain Monte Carlo sampling algorithm in PyBioNetFit

Abstract Summary Bayesian inference in biological modeling commonly relies on Markov chain Monte Carlo (MCMC) sampling of a multidimensional and non-Gaussian posterior distribution that is not analytically tractable. Here, we present the implementation of a practical MCMC method in the open-source software package PyBioNetFit (PyBNF), which is designed to support parameterization of mathematical models for biological systems. The new MCMC method, am, incorporates an adaptive move proposal distribution. For warm starts, sampling can be initiated at a specified location in parameter space and with a multivariate Gaussian proposal distribution defined initially by a specified covariance matrix. Multiple chains can be generated in parallel using a computer cluster. We demonstrate that am can be used to successfully solve real-world Bayesian inference problems, including forecasting of new Coronavirus Disease 2019 case detection with Bayesian quantification of forecast uncertainty. Availability and implementation PyBNF version 1.1.9, the first stable release with am, is available at PyPI and can be installed using the pip package-management system on platforms that have a working installation of Python 3. PyBNF relies on libRoadRunner and BioNetGen for simulations (e.g. numerical integration of ordinary differential equations defined in SBML or BNGL files) and Dask.Distributed for task scheduling on Linux computer clusters. The Python source code can be freely downloaded/cloned from GitHub and used and modified under terms of the BSD-3 license (https://github.com/lanl/pybnf). Online documentation covering installation/usage is available (https://pybnf.readthedocs.io/en/latest/). A tutorial video is available on YouTube (https://www.youtube.com/watch?v=2aRqpqFOiS4&t=63s). Supplementary information Supplementary data are available at Bioinformatics online.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Electron scattering on A = 3 nuclei from quantum Monte Carlo based approaches

We perform first-principle calculations of electron-nucleus scattering on 3 He and 3 H using the Green's function Monte Carlo method and two approaches based on the factorization of the final hadronic state: the spectral-function formalism and the short-time approximation. These three methods are benchmarked among each other and compared to the experimental data for the longitudinal and transverse electromagnetic response functions of 3 He, and the inclusive cross sections of both 3 He and 3 H. Since these three approaches are based on the same description of nuclear dynamics of the initial target state, comparing their results enables a precise quantification of the uncertainties inherent to factorization schemes. At sufficiently large values of the momentum transfer, we find an excellent agreement of the Green's function Monte Carlo calculation with experimental data and with both the spectral-function formalism and the short-time approximation. Here, we also analyze the relevance of relativistic effects, whose inclusion becomes crucial to explain data at high momentum and energy transfer.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Predicting the neutrinoless double-β-decay matrix element of 136 Xe using a statistical approach

Calculation of the nuclear matrix elements (NMEs) for double-β decay is of paramount importance for guiding experiments and for analyzing and interpreting the experimental data, especially for the search of the neutrinoless double β decay mode (0v⁢ββ). However, there are currently still large differences between the NME values calculated by different methods, hence a quantification of their uncertainties is very much required. Here, in this paper, we propose a statistical analysis of 0⁢v⁢ββ NME for the 136 Xe isotope, based on the interacting shell model, but using three independent effective Hamiltonians, emphasizing the range of the NMEs' most probable values and its correlations with observables that can be obtained from the existing nuclear data. Consequently, we propose a common probability distribution function for the 0⁢vββ NME, which has a range of (1.55–2.65) at 90% confidence level, with a mean value of 1.99 and a standard deviation of 0.37.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Multimodal fission from self-consistent calculations

When multiple fission modes coexist in a given nucleus, distinct fragment yield distributions appear. Multimodal fission has been observed in a number of fissioning nuclei spanning the nuclear chart, and this phenomenon is expected to affect the nuclear abundances synthesized during the rapid neutron-capture process (𝑟-process). In this study, we generalize the previously proposed hybrid model for fission-fragment yield distributions to predict competing fission modes and estimate the resulting yield distributions. Here, our framework allows for a comprehensive large-scale calculation of fission-fragment yields suited for 𝑟-process nuclear network studies. Nuclear density functional theory is employed to obtain the potential energy and collective inertia tensor on a multidimensional collective space defined by mass multipole moments. Fission pathways and their relative probabilities are determined using the nudged elastic band method. Based on this information, mass and charge fission yields are predicted using the recently developed hybrid model. Fission properties of fermium isotopes are calculated in the axial quadrupole-octupole collective space for three energy density functionals (EDFs). Disagreement between the EDFs appears when multiple fission modes are present. Within our framework, the UNEDF⁢1 HFB EDF agrees best with experimental data. Calculations in the axial quadrupole-octupole-hexadecapole collective space improve the agreement with the experiment for SkM*. We also discuss the sensitivity of fission predictions on the choice of EDF for several superheavy nuclei. Fission-fragment yield predictions for nuclei with multiple fission modes are sensitive to the underlying EDF. For large-scale calculations in which a minimal number of collective coordinates is considered, UNEDF⁢1 HFB provides the best description of experimental data, though the sensitivity motivates robust quantification of the uncertainties of the theoretical model.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Data-driven reevaluation of f t values in superallowed β decays

We present a comprehensive reevaluation of the ft values in superallowed nuclear β decays crucial for the precise determination of V ud and low-energy tests of the electroweak standard model. It consists of the first, fully data-driven analysis of the nuclear β decay form factor, that utilizes isospin relations to connect the nuclear charged weak distribution to the measurable charge distributions. This prescription supersedes previous shell-model estimations, and allows for a rigorous quantification of theory uncertainties in f which is absent in the existing literature. Furthermore, our new evaluation shows an overall downward shift of the central values of f at the level of 0.01%.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

DEPRECATED AI-Batt-OS (Autonomous Identification of Battery Life Models - Open Source) [SWR 21-17]

DEPRECATED. This repository was archived by the owner on Jun 30, 2026. It is now read-only. Open source implementation of some of the methods utilized by AI-Batt, a battery lifetime modeling and analysis toolkit provided by the National Laboratory of the Rockies (NLR). This software demonstrates the use of bi-level optimization and symbolic regression techniques to semi-autonomously identify algebraic models predicting the capacity fade of lithium-ion batteries during calendar aging. Modeling the degradation of batteries is a complex task, due to the difficulty in separating the time-dependent and time-independent factors impacting cell level degradation, across multiple data series with different numbers of measurements and/or data quality. Bi-level optimization enables model parameters to be optimized to either the entire data set or to individual data series, allowing statistical disambiguation of global behaviors (data series independent) and local behaviors (data series dependent). Symbolic regression is used to automatically search for optimal low-dimesional models predicting the variation of locally optimized parameters versus time-independent experimental variables from millions of possible models, resulting in a more accurate and repeatable model identification process than is possible by a manual search. The provided tools also implement cross-validation and bootstrap resampling schemes, empowering statistical model comparison/selection and quantification of model uncertainties. An example script replicates the results from the manuscript "Challenging Practices of Algebraic Battery Life Models through Statistical Validation and Model Identification via Machine-Learning", submitted to ECS. All code is written in MATLAB. Requires the Statistics and Machine Learning Toolbox. Contact Dr. Paul Gasper at Paul.Gasper@nlr.gov for any questions.

Gasper, Paul↗