Modeling Fatigue Crack Growth Using Symbolic Regression
Explore the source record for details and available documents.
Keven Carlson↗
SEARCH · Engineering Papers
Results for “regression”
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.
Combining electrochemistry and data-sparse Gaussian process regression for lithium-ion battery hybrid modeling
Explore the source record for details and available documents.
Jackson Fogelquist↗
Overview of initial research into the effects of strong vortex flow on hybrid rocket combustion and performance
An examination of the effect of vortex flow on hybrid rocket combustion and performance is underway. Emphasis is on response of the fuel regression rate when subjected to vortex flow. Initial results show that there is a definite effect of the vortex on fuel regression rate. Future work will focus on quantitatively measuring this regression rate. This work is part of an overall program to develop an ultra low cost fuel system for hybrid rocket engines.
P Gloyer↗
An Automated Probabilistic Asteroid Prediscovery Pipeline
We present an automated and probabilistic method to make prediscovery detections of near-Earth asteroids (NEAs) in archival survey images, with the goal of reducing orbital uncertainty immediately after discovery. We refit the Minor Planet Center's astrometry and propagate the full six-parameter covariance to survey epochs to define search regions. We build low-threshold source catalogs for viable images and evaluate every detected source in a search region as a candidate prediscovery. We eliminate false positives by refitting a new orbit to each candidate and probabilistically linking detections across images using a likelihood ratio. Applied to the Zwicky Transient Facility's (ZTF) imaging, we identify approximately 3000 recently discovered NEAs with prediscovery potential, including a doubling of the observational arc for about 500. We use archival ZTF imaging to make prediscovery detections of the potentially hazardous asteroid 2021 DG1, extending its arc by 2.5 yr and reducing future apparition sky plane uncertainty from many degrees to arcseconds. We also recover 2025 FU24 nearly 7 yr before its first known observation, when its sky plane uncertainty covers hundreds of square degrees across thousands of ZTF images. The method is survey agnostic and scalable, enabling rapid orbit refinement for new discoveries from Rubin, NEO Surveyor, and NEOMIR.
79 ASTRONOMY AND ASTROPHYSICS↗
A 25,000 cycle life battery
Regression analysis of cycle life data for long life 20 ampere hour sealed nickel cadmium battery
Carr, E.↗
Spatially Refined Satellite Gravimetry Captures Human Signatures in Global Terrestrial Water Storage Trends
Human activities have directly altered the water cycle through water management, aquifer pumping, agricultural irrigation, and land use change. Although satellite gravimetry has transformed global hydrological research, its coarse resolution limits attribution of freshwater change to human activities at many management-relevant scales. Here we assessed global terrestrial water storage (TWS) trends from April 2002 to November 2025 using “stacked” regression of Level-1B intersatellite ranging data, which leverages temporal information and variability to dramatically improve effective spatial resolution relative to standard approaches. We combined this refined product with rigorous uncertainty analysis, autocorrelation-robust geostatistical methods, and literature assessment to evaluate TWS trend associations with land and water use, climate variability, and glacial mass loss. We identified TWS trend hotspots exhibiting significant spatial associations with anthropogenic 40 drivers, including groundwater and surface-water irrigation, rainfed agriculture, deforestation, and reservoir impoundment. Across these regions, cumulative TWS losses (3,122 Gt) substantially exceeded gains (2,432 Gt). Compared with traditional regression of monthly mascons, our approach yielded regional trend magnitudes that are on average 33% larger, revealing that global freshwater depletion, particularly from groundwater pumping, is considerably more acute than previously estimated. Multivariate regression models show that humans account for a significant share of the spatial variability in TWS trends on every non-polar continent except Australia. We detected localized TWS gains linked to rainfed agriculture, surface water irrigation, and reservoir filling that were unresolved in earlier gravimetric studies. The methodology provides a foundation for future gravity missions to independently track decadal freshwater change with unprecedented spatial fidelity.
groundwater↗
A Tutorial on Bayesian analysis of linear shock compression data
Gas gun and other shock compression experiments often produce shock wave velocity measurements that are linearly associated with particle velocity. Traditionally, this empirical relationship is quantified with a single Hugoniot curve that is estimated using least squares regression. However, for downstream modeling and simulation tasks, it is often more useful to have multiple Hugoniot curves in the pressure–volume plane that are consistent with the data. We employ Bayesian uncertainty quantification methods as a framework for propagating measurement uncertainty through to model parameters and predictions. Specifically, this Tutorial shows how to sample multiple Hugoniot curves in the pressure–volume plane that are consistent with the shock wave-particle velocity measurements in a two-step Bayesian approach. First, we obtain an analytical expression for the posterior distribution of the linear model parameters using Bayesian linear regression. Second, we propagate samples from the posterior distribution through the Rankine–Hugoniot equations to yield Hugoniot curves in the pressure–volume plane. The procedure is demonstrated with publicly available data on argon, copper, and nickel, and compared against bootstrapping and linear regression. The Bayesian procedure is shown to be interpretable, computationally inexpensive, and less sensitive than an alternative bootstrapping approach to the removal of the point in the copper dataset that has the largest particle velocity. As a Tutorial on Bayesian methodology for the shock compression community, we provide several derivations and explanations that make this paper self-contained, and make all code and data available at github.com/llnl/BALSCD.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Combustion Characteristics of a Solid Propellant With a Charring Binder
A brief investigation of the combustion characteristics of a solid propellant containing a binder which chars, as opposed to melting or volatizing, has been made. The burning rate of the propellant with the charring binder was significantly higher than similar propellants containing non-charring binders. High speed motion pictures of the burning propellant showed that the aluminum burned on the regressing surface, rather than a short distance from it as is typical with composite propellants.
D E Udlock↗
Estradiol associations with brain functional connectivity in postmenopausal women
Abstract Objective Previous studies have found that estrogens play a role in functional connectivity in the brain; however, little research has been done regarding how estradiol is associated with functional connectivity in postmenopausal women. The purpose of this study was to examine the relationship between estradiol and functional connectivity in postmenopausal women. Methods Structural and blood oxygenation level–dependent resting-state magnetic resonance imaging scans of 88 cognitively healthy postmenopausal individuals were obtained along with blood samples collected the same day as the magnetic resonance imaging to assess hormone levels. We generated connectivity values in CONN toolbox version 20.b, an SPM-based software. Results A regression analysis was run using estradiol level and regions of interest (ROI), including the hippocampus, parahippocampus, dorsolateral prefrontal cortex, and precuneus. Estradiol level was found to enhance parahippocampal gyrus anterior division left functional connectivity during ROI-to-ROI regression analysis. Estradiol enhanced functional connectivity between the parahippocampal gyrus anterior division left and the precuneus as well as the parahippocampal gyrus anterior division left and parahippocampal gyrus posterior division right. An exploratory analysis showed that years since the final menstrual period was related to enhanced connectivity between regions within the frontoparietal network. Conclusions These results illustrated the relationship between estradiol level and functional connectivity in postmenopausal women. They have implications for understanding how the functioning of the brain changes for individuals after menopause that may eventually lead to changes in cognition and behavior in older ages.
Obstetrics & Gynecology↗
A Data-Driven Method for Modeling Creep-Fatigue Stress- Strain Behavior Using Neural ODEs
In this paper, we introduce a data-driven machine learning approach for modeling one-dimensional stress–strain behavior under cyclic loading, utilizing experimental data from the nickel-based Alloy 617. The study employs uniaxial creep–fatigue test data acquired under various loading histories and compares two distinct neural network-based ODE models. The first model, known as the black-box model, comprehensively describes the strain–stress relationship using a Neural ODE equation. To interpret this black-box model, we apply the Sparse Identification of Nonlinear Dynamical Systems (SINDy) technique, transforming the black-box model into an equation-based model using symbolic regression. The second model, the Neural flow rule model, incorporates Hooke’s Law for the linear elastic component, with the nonlinear part characterized by a Neural ODE. Both models are trained with experimental data to accurately reflect the observed stress–strain behavior. We conduct a detailed comparison with the standard Chaboche model, which includes three back stresses. Our results demonstrate that the neural network-based ODE models precisely capture the experimental creep–fatigue mechanical behavior, exceeding the standard Chaboche model’s accuracy. Furthermore, an interpretable model derived from the black-box neural ODE model through symbolic regression achieves accuracy comparable to the Chaboche model, enhancing its interpretability. The results highlight the potential of neural network-based ODE models to depict complex creep–fatigue behavior, eliminating the necessity for experts to define a specific, material-focused model form.
creep-fatigue↗
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 [National Renewable Energy Lab. (NREL↗
The cluster decomposition of the configurational energy of multicomponent alloys
Abstract The cluster expansion method (CEM) is a widely used lattice-based technique in the study of multicomponent alloys. Despite its prevalent use, a clear understanding of expansion terms is lacking. We present a modern mathematical formalism of the CEM and introduce thecluster decomposition—a unique and basis-independent decomposition for functions of the atomic configuration in a crystal. We identify the cluster decomposition as an invariant ANOVA decomposition; and demonstrate how functional analysis of variance and sensitivity analysis can be used to interpret interactions among species. Furthermore, we show how the mathematical structure of the cluster decomposition enables numerical evaluation that scales with the number of clusters and is independent of the number of species. Overall, our work enables rigorous interpretations of interactions among species, provides opportunities to explore parameter estimation beyond linear regression, introduces a numerical efficient implementation, and enables analysis of cluster expansions based on established mathematical and statistical principles.
Chemistry↗
AI-Batt (Autonomous Identification of Battery Life Models) [SWR 21-36]
Autonomous Identification of Battery Life Models (AI-Batt) AI-Batt is a MATLAB code base for developing lifetime models for batteries from accelerated aging data. The code base provides many functions for processing, visualizing, and modeling battery aging data, making the data processing, exploration, and modeling workflow substantially faster. These tools are tailored for working with battery aging data sets, which usually consist of many separate time-series for each cell, with many test conditions and possible replicates at each condition, which makes it difficult to simply process or visualize the data set. Complex modeling tasks, such as cross-validation, sensitivity analysis, and uncertainty quantification have been implemented to enable thorough statistical investigation of model predictions. Additionally, several machine-learning algorithms are implemented to autonomously identify suitable models via symbolic regression. Data processing functions automatically cast data from the struct data type, which is commonly used to store experimental data, but is not an acceptable input for most algorithms, to the table data type, which can be easily used as input to any optimization algorithm. Also, the data can be separated into time-invariant and time-variant data tables, which is helpful for exploring the data set as well as developing separate models for time-variant and time-invariant aging mechanisms. For example, in aging tests with constant temperature, temperature is a time-invariant experimental condition. Visualization tools enable plotting of data, model fits, and model simulations possible with single-line function calls, empowering data exploration of complex data sets with both time-varying and time-invariant trends. Plots can be automatically generated for the whole data set, or separated by data group (groups of test replicates) or individual data series. Data points or data series can be automatically colored by the value of a variable with a variety of color maps, and model predictions can also be colored by the value of a fit statistic. Comparisons between data sets and the predictions/simulations of different models on the same data set can be easily plotted as well. Distributions of parameter values from bootstrap resampling can be plotted to visualize the reliability of parameter estimation, or determine any correlations between parameters. Modeling tools handle the complex task of creating and parsing symbolic equations for modeling battery lifetime. Equations are parsed to grab relevant data variables, parameter values, or specified sub-models for input into optimization, evaluation, or simulation functions. Models can be optimized locally (one set of parameters for each data series), bi-level (some parameters shared across the data set), or globally (single set of parameters for all data). Functions implementing symbolic regression algorithms help users to discover effective model equations, even in poorly sampled, high-dimensional data.
Smith, Kandler [National Renewable Energy Lab. (NR↗
Flash Heating of a Hollow Rod to Determine Thermal Properties
In a nuclear reactor, prolonged exposure to high neutron flux levels can change material properties of reactor components. Due to the danger posed to personnel by highly radioactive components such as spent fuel, noncontact methods must be used to determine these properties. The objective of the current research is to test the viability of using a flash-heating method to determine the thermal diffusivity of a silicon carbide (SiC) tube, which could be used as cladding for a nuclear fuel rod. The flash-heating experiment was conducted at Oak Ridge National Laboratory and temperature results were analyzed as a two-Dimensional cross section which assumed axial symmetry. The model was fitted iteratively to the temperature measurements using nonlinear regression which required values for thermal diffusivity, Biot number and a term to determine the magnitude of the heat absorbed by the flash. The method could be expanded in the future to estimate the thermal diffusivity of spent nuclear fuel rods in an effort to increase the certainty of new reactor designs using precise thermal properties throughout the core life. A statistical analysis of the results of this work is provided as part of the analysis.
McMasters, Robert L [Virginia Military Institute, ↗
A physically interpretable precursor framework for sub-seasonal prediction of Northern Hemisphere flash flourishing
Flash flourishing describes rapid vegetation increases that can quickly reshape land–atmosphere exchanges and impacts on ecosystem, yet its large-scale precursors, circulation context, and sub-seasonal predictability remain poorly understood. Here, we identified onset-stage circulation regimes across northern extratropical latitudes (NEL; >30°N) using 200 and 1000 hPa geopotential height, and examined their regional expressions over eastern Asia, western North America, and Europe. Flash flourishing onset in East Asian was associated with a baroclinic circulation regime and was preceded by a North Atlantic sea surface temperature (SST) precursor at a four-pentad lead. In contrast, onset in western North American and European preferentially occurred under barotropic regimes, preconditioned by Great Plains soil moisture at three-pentad lead and North Atlantic SST at a four-pentad lead, respectively. Ridge regression forecasts revealed regime-dependent sub-seasonal predictability, with mean out-of-sample R 2 exceeding 0.3 up to lead times of two pentads in East Asia, three pentads in western North America, and four pentads in Europe. Together, these findings established a mechanistic and regionally specific framework for anticipating rapid vegetation greening at sub-seasonal timescales.
Kong, Xiangxu [Nanjing Univ. of Information Scienc↗
Improving Adhesive Bondline Time of Flight Predictions During Autoclave Cure Utilizing Machine Learning
Composite materials are increasingly being used in aerospace applications due to their superior strength-to-weight ratio compared to commonly used metals. A current limitation to widespread adoption is the certification of adhesively bonded joints. One approach to improving adhesive bonding in composites is accurately measuring the thickness of adhesive bondlines in composite laminates. Precise bondline thickness control is essential for aerospace applications where adhesive layer thickness directly affects joint fracture properties and structural performance. This study focused on implementing machine learning techniques to determine the ultrasonic time of flight (directly correlated to thickness) in adhesive bondlines throughout autoclave cure cycles. A high-temperature (use up to 180°C) ultrasonic scanning system was deployed in an autoclave to provide time of flight data through composite panels. Three experiments were conducted on the curing of 305 mm × 305 mm unidirectional composite panels. In the first experiment, a piecewise function was fit for the temperature correction factor to account for changing autoclave temperatures. Due to deficiencies in the first calibration experiment, a second experiment was run, and the results were used to train a machine learning model. The revised experiment, in combination with the machine learning model, significantly increased the accuracy of the bondline time of flight predictions (~14% error reduced to <1%). Data was processed using the Regression Learner Application in MATLAB®, with a Support Vector Machine selected for the model. The result was a machine learning algorithm capable of reliably quantifying ultrasonic time of flight through adhesive bondlines. The third experiment provided independent test data for the machine learning model, demonstrating that the model produces accurate predictions from data beyond its training set.
Machine Learning↗
Optimal binning of correlated measurements
Experimental measurements are commonly represented on a discrete grid, requiring a balance between granularity and statistical noise. Two strategies have traditionally been used to improve such representations: selecting an appropriate bin width to control discretization error and applying kernel-based smoothing to suppress fluctuations. Despite their shared goal, these approaches have largely developed independently, without a unified statistical description of how discretization and correlation jointly determine measurement precision. Here, we extend the discussion of optimal interval averaging to a correlation-aware setting by Gaussian process regression, which explicitly accounts for correlations among neighboring bins. Starting from first principles, we derive the mean-squared error of discretized measurements and obtain closed-form asymptotic expressions for the optimal bin width and correlation length. When recast in reduced variables, the theory reveals distinct universal scaling laws governing the error in the correlation-free and correlation-controlled regimes. Characterized by intrinsically smooth intensity profiles and counting-based statistics, neutron scattering measurements are well suited for demonstrating the enhanced error contraction enabled by inter-bin correlations. We show that such improvement is achievable over the experimentally accessible Q-range and across multiple instruments and material systems. These results show that explicitly accounting for correlations systematically reshapes the limits of precision in discretized, noise-limited measurements. More broadly, the framework provides a transferable statistical foundation for optimizing data representation, inference, and experimental design across the physical and data sciences.
Tung, Chi-Huan [ORNL] (ORCID:0000000221972074)↗
Hybrid Data‐Driven Discovery of High‐Performance Silver Selenide‐Based Thermoelectric Composites
Optimizing material compositions often enhances thermoelectric performances. However, the large selection of possible base elements and dopants results in a vast composition design space that is too large to systematically search using solely domain knowledge. To address this challenge, a hybrid data-driven strategy that integrates Bayesian optimization (BO) and Gaussian process regression (GPR) is proposed to optimize the composition of five elements (Ag, Se, S, Cu, and Te) in AgSe-based thermoelectric materials. Data is collected from the literature to provide prior knowledge for the initial GPR model, which is updated by actively collected experimental data during the iteration between BO and experiments. Within seven iterations, the optimized AgSe-based materials prepared using a simple high-throughput ink mixing and blade coating method deliver a high power factor of 2100 µW m −1 K −2 , which is a 75% improvement from the baseline composite (nominal composition of Ag 2 Se 1 ). In conclusion, the success of this study provides opportunities to generalize the demonstrated active machine learning technique to accelerate the development and optimization of a wide range of material systems with reduced experimental trials.
36 MATERIALS SCIENCE↗