Search NASA⌕ Search

SEARCH · Search NASA

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.

At least 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Online LIBS–ML Framework for Dynamic Characterization of Heterogeneous Waste-Derived Gasification Feedstocks

LIBS−ML framework for real time feedstock characterization during continuous conveyor transport Heterogeneous waste derived feedstocks (e.g., waste coal, biomass and blends) introduce rapid variability in heating value and ash chemistry that affect gasifier operation, yet conventional laboratory characterization techniques are too slow to support proactive control. To address this gap, this study reports on an online, in situ, dynamic characterization framework that couple’s laser-induced breakdown spectroscopy (LIBS) with leakage safe machine learning (ML) regression to deliver real time, decision quality predictions of gasifier relevant properties. A controlled sample matrix spanning two different waste coals, two different biomasses, and engineered blends under two particle size conditions were constructed and benchmarked using standardized laboratory analyses for proximate/ultimate properties and ash composition. LIBS spectra were acquired dynamically as material flowed on a conveyor belt, using high energy 1064 nm laser ablation and shot averaging to improve repeatability and precision. Supervised regression models (multi layer perceptron (MLP) /artificial neural network (ANN), random forest (RF), and support vector regression (SVR)) and an optimized weighted ensemble were trained on emission line feature sets using nested cross validation with Bayesian hyperparameter tuning and validated against an independent hold out set. The proposed LIBS−ML workflow achieves near laboratory predictive fidelity across parametric targets (including higher heating value (HHV), ash content, fixed carbon, sulfur, major ash forming oxides, and initial deformation temperature (IDT)), with the weighted ensemble providing a robust default predictor under dynamic measurement conditions. These results demonstrate a practical pathway for real time feedstock characterization that can enable feedforward adjustments and more resilient gasifier operation for variable quality waste derived fuels.

Biomass↗
DOE OSTI2026

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↗
DOE OSTI2024

Local primordial non-Gaussianity from the large-scale clustering of photometric DESI luminous red galaxies

ABSTRACT We use angular clustering of luminous red galaxies from the Dark Energy Spectroscopic Instrument (DESI) imaging surveys to constrain the local primordial non-Gaussianity parameter fNL. Our sample comprises over 12 million targets, covering 14 000 deg2 of the sky, with redshifts in the range 0.2 < z < 1.35. We identify Galactic extinction, survey depth, and astronomical seeing as the primary sources of systematic error, and employ linear regression and artificial neural networks to alleviate non-cosmological excess clustering on large scales. Our methods are tested against simulations with and without fNL and systematics, showing superior performance of the neural network treatment. The neural network with a set of nine imaging property maps passes our systematic null test criteria, and is chosen as the fiducial treatment. Assuming the universality relation, we find $f_{\rm NL} = 34^{+24(+50)}_{-44(-73)}$ at 68 per cent (95 per cent) confidence. We apply a series of robustness tests (e.g. cuts on imaging, declination, or scales used) that show consistency in the obtained constraints. We study how the regression method biases the measured angular power spectrum and degrades the fNL constraining power. The use of the nine maps more than doubles the uncertainty compared to using only the three primary maps in the regression. Our results thus motivate the development of more efficient methods that avoid overcorrection, protect large-scale clustering information, and preserve constraining power. Additionally, our results encourage further studies of fNL with DESI spectroscopic samples, where the inclusion of 3D clustering modes should help separate imaging systematics and lessen the degradation in the fNL uncertainty.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

A Machine Learning based Approach of Estimating Equivalent Circuit Model Parameters at Different SoCs of Li-ion Batteries from Voltage Relaxation

Abstract: In this study, an approach of estimating the equivalent circuit model (ECM) parameters for Li-ion batteries (LIBs) is proposed based on the voltage value at different intervals while relaxing the LIB after discharge. The typical approach for estimating ECM parameters of a LIB is to conduct electrochemical impedance spectroscopy (EIS) measurements at different frequencies and fit them to a predefined circuit model, which requires additional measuring arrangements and specialized devices. The proposed methodology utilizes four different voltages at 0s, 60s, 360s, and 1800s alongside the specific state of charge (SoC) value for a specific constant discharge current value of ~1C until the relaxation stage to train and evaluate three regression-based machine learning models— Support Vector Regression (SVR), Extreme Gradient Boosting (XGBoost), and Gaussian Process Regression (GPR)—for estimating the ECM parameters of the selected model. Bayesian optimization is employed for hyperparameter tuning to achieve optimal performance for all the regressor models, among which, the GPR provided the best performance with the root-mean-squared error (RMSE) of less than 4x10-4 on average for the resistive components and less than 0.27 for capacitive components with excellent R2 scores. The simplicity of the approach enables it to eliminate the need for sophisticated measuring equipment and computation power.

Sagar, Md. Samiul [The University of Alabama (UA)]↗
DOE OSTI2025

Physical Interpretation of Early Battery Life Prediction Models

Early battery life prediction models are most useful for R&D if they help us understand the early changes in battery electrochemical response that correspond with long-term degradation and failure. Linear regression models such as Fused lasso and Partial Least Squares can fit coefficients directly to high-dimensional electrochemical data like capacity-voltage and ΔV–state-of-charge, i.e., Q(V) and ΔV(SOC) curves, learning coefficients that can be physically interpreted. We leverage the ISU-ILCC battery aging data set to learn high-dimensional coefficients for early battery life prediction from traditional slow-rate capacity check data, demonstrating learning on Q(V), d Q· d V −1 , and ΔV(SOC) curves. A thorough study on the dependence of coefficient values on train/test size and data preprocessing methods is made, demonstrating the reliability of high-dimensional regression approaches unless very small amounts of data are used for model training. For this data set, coefficients from Q(V) and d Q· d V −1 models highlight changes in electrode stoichiometry due to lithium loss, while ΔV(SOC) coefficients highlight changes in positive electrode diffusivity due to particle cracking as well as electrode stoichiometry shifts. By directly interpreting the coefficients of a regression model, we make physical insights into battery degradation mechanisms without requiring the assumptions of traditional battery data analysis methods.

25 ENERGY STORAGE↗
DOE OSTI2024

SysCaps (Language Interfaces for Simulation Surrogates of Complex Systems) [SWR-24-97]

You've found the official code repository for the paper "SysCaps: Language Interfaces for Simulation Surrogates of Complex Systems," presented at the Foundation Models for Science: Progress, Opportunities, and Challenges workshop at NeurIPS 2024. Our paper conjectures that interfaces (both text templates as well as conversational) makes interacting with simulation surrogate models for complex systems more intuitive and accessible for both non-experts and experts. "System captions", or SysCaps, are text-based descriptions of systems based on information contained in simulation metadata. Our paper's goal is to train multimodal regression models that take text inputs (SysCaps) and timeseries inputs (exogenous system conditions such as hourly weather) and regress timeseries simulation outputs (e.g. hourly building energy consumption). The experiments in our paper with building and wind farm simulators, which can be reproduced using this codebase, aim to help us understand whether a) accurate regression in this setting is possible and b) if so, how well can we do it. Paper: https://arxiv.org/abs/2405.19653

Emami, Patrick↗
DOE OSTI2025

Decayheatml

This code is designed to predict and analyze the decay heat generated in molten salt reactors (MSRs) using a hybrid approach that combines machine learning and segmented polynomial fitting. The accurate prediction of decay heat is essential for reactor safety and the optimization of spent fuel storage. The code operates through several key components: 1) Data Architecture: It incorporates a modular data architecture that handles various MSR-specific operational parameters such as power density, humidity content, and air ingress. These parameters are sampled using Sobol sequences to ensure comprehensive coverage of operational uncertainties. 2) Machine Learning Framework: The code employs a diverse set of machine learning models, including polynomial regression, decision trees, random forests, gradient boosting, support vector regression, k-nearest neighbors, multi-layer perceptrons, and symbolic regression. These models are trained to predict decay heat over a wide temporal range, from immediate shutdown up to 10,000 years. 3) Region-Optimized Training: The temporal domain is divided into multiple regions, each modeled separately to capture distinct decay heat characteristics across different time scales. This approach significantly improves the accuracy and interpretability of predictions. 4) Segmented Polynomial Interpretation (SPI): The SPI method translates machine learning predictions into piecewise polynomial equations. These equations are physically interpretable and can be directly integrated into existing engineering workflows and safety analyses. 5) Front-End Interfaces: The code includes both a Jupyter notebook interface for research development and a Streamlit web application for operational deployment. These interfaces allow users to interactively explore decay heat predictions, adjust operational parameters, and visualize results in real-time. 6) Applications: The framework supports various applications, including safety system validation and spent fuel container optimization. It enables real-time evaluation of worst-case decay heat scenarios, informing the design of passive safety systems and optimizing container designs for long-term storage. Overall, this code provides a robust, accurate, and user-friendly tool for predicting decay heat in MSRs, enhancing reactor safety, and optimizing spent fuel management.

Retamales, Mauricio Eduardo Tano [Idaho National L↗
DOE OSTI2025

Field Validation of Thermoelectric Generation System at Holcim Cement Plant in Alpena, Michigan

Executive Summary Project Background The Industrial Technology Validation (ITV) program aims to identify and demonstrate the performance of new, emerging, and underutilized energy-saving technologies in the industrial sector to help inform decisions to help accelerate their commercialization and deployment, as well as to help make industries more competitive. This ITV demonstration evaluated a thermoelectric generation (TEG) technology at a cement plant, aiming to reduce energy demand in the cement industry. A median cement plant consumes 5.73 million British thermal units per ton of clinker production (resulting in 0.838 metric tons of carbon dioxide [CO₂] emissions per ton of clinker) (Boyd and Zhang 2011, EPA 2021), equivalent to approximately 6.9 trillion British thermal units (TBtu) per year in energy consumption at a cement plant producing 3,300 tons of clinker per day.¹ Collaborating with Holcim, Advanced Thermovoltaic Systems (ATS) developed and deployed a pilot-scale thermoelectric power system to efficiently capture and convert waste heat to electricity. The system leverages the Seebeck effect to convert temperature differences on two sides of semiconductor cartridges into electrical power (ScienceDirect, n.d.). This generation is realized with minimal moving parts compared to existing waste-heat-to-generation solutions and allows capture from heat sources with temperatures as low as 150°C. This project aimed to validate a scalable solution applicable for capturing medium-temperature waste heat, including ambient losses from other high-temperature processes, and high-temperature sources less suitable for other waste-heat-to-power solutions. By recovering this otherwise wasted heat, this project intends to validate improvements to overall process efficiency through reduction in purchased electricity, thereby reducing operational costs while enhancing resiliency and competitiveness. Description and Scope This study evaluated the performance of a TEG system from ATS as a solution to convert waste heat into useful power at a Holcim cement plant in Alpena, Michigan. This plant is a fully integrated cement plant that has been operating since 1907. The facility operates continuously (24/7/365) with approximately 250 employees and five long dry kilns, yielding a total production capacity of 7,852 tons of cement per day (EPA 2023). Currently, the Alpena plant uses waste heat boilers to convert waste heat from the exhaust of each kiln into steam, which drives steam turbine generators. The ATS TEG is being evaluated for its potential to supplement the steam turbines by capturing the remaining lower grade heat. This technology is also being considered for other Holcim plants where steam turbines are not a viable option. ATS installed a pilot-scale TEG unit with an array of 582 individual thermoelectric semiconductor cartridges, of which 573 were operational. The cartridges are sandwiched between 48 hot plates and 49 cold plates. Each cartridge is designed to generate 20 watts (W) of gross power at a hot-side temperature of 240°C and cold-side temperature of 20°C. As such, the total gross generation capacity of the installed system is 11.5 kilowatts (kW) at design conditions. The system configuration for the evaluation was designed to prioritize convenience of installation and minimize disruption to production at the site, while ensuring that the heat required can be obtained for evaluating the TEG system at various operational conditions. To accomplish this, a portion of the steam supplied to Alpena’s steam turbine generation system was diverted to be used as the heat source for the TEG system, while water was supplied to the cold side of the system from nearby Lake Huron. This configuration was designed for the evaluation of the pilot-scale system to assess the performance at different conditions. A commercial-scale system will likely vary from the pilot system depending on typical configurations, including both scale and application. Future commercial applications of the ATS system would involve integrating the system into the exhaust from kiln preheaters, clinker coolers, or radiant heat capture from kiln shells for the heat source. For the cold source, a range of cooling solutions can be considered, including a mechanical cooling system, depending on the location and the application. To increase the generation capacity for commercial applications, the technology provider is working toward developing a commercial-scale TEG system, which would combine multiple TEG units (each similar in design to the pilot system) together. The scope of this evaluation includes the pilot-scale TEG system and all impacted equipment including pumps, controllers, and power handling equipment. Study Objectives The evaluation's goal was to assess the potential of the ATS TEG system to generate useful electrical power by capturing waste heat from cement production kilns. The objectives of this study are to evaluate and verify the following claims made by ATS regarding the pilot-scale system installed at the Holcim Alpena plant. The following design parameters and claims are also outlined in Table ES- 1 and Table ES- 2: • Gross Power: The thermoelectric system converts heat into power to create gross power, the total measured power generated by the system. The 573 active cartridge pilot-scale system is expected to generate 11.5 kW of gross power at the designed hot-side temperature of 240°C and cold-side temperature of 20°C. Power production is dependent on the temperature difference between the heat source (ultimately from the waste heat) and cold temperature supply source. • Net Power: The net power is the total usable power provided to the site by the TEG system after deducting parasitic power loads from the gross generated power. Supplementary equipment is required to operate the TEG system including pumps, controllers, and, in certain anticipated applications, mechanical cooling, which introduce parasitic loads to system operation. After deducting the parasitic loads from the gross power generation, ATS anticipates achieving a net power generation of 7.5 kW from the pilot-scale system. • Thermal Efficiency: The thermal efficiency is the percent of the total heat transferred to the TEG system that is converted to gross power. Historically, TEGs have a thermal efficiency of 2%–5% (DOE 2008). Prior industrial-scale TEG systems, such as the E1 TEG offered by Alphabet Energy, operated at an efficiency of 2.5% (Lamonica, 2014). ATS anticipates achieving an average efficiency of 4.8% or higher in converting heat energy to usable electricity. • Cartridge Performance: The TEG system comprises 573 active individual semiconductor cartridges, each of which generates a portion of the total power. Cartridge optimization and selection is an important design consideration for potential future TEG system design performance. Therefore, understanding the distribution of gross power and efficiency within the pilot system is vital to understanding what is achievable. At a design hot-side temperature of 240°C and cold-side temperature of 20°C, ATS anticipates a cartridge performance of 20 W of gross power per cartridge at an efficiency of 4.8% per cartridge. In addition to evaluating the claimed performance of the TEG pilot-scale unit, the study estimated the potential annual impacts of a scaled-up commercial system used to capture kiln waste heat over annual operations. The evaluation estimated the gross and net annual electric generation achievable by capturing heat from the two proposed tap-in points: the kiln exhaust and the clinker cooler exhaust; see Section 2.1 for details. Two use cases were examined: • Holcim Alpena: The Holcim Alpena site consists of long dry kilns with superheater boilers, which differs from the rest of Holcim’s cement plant portfolio and results in lower waste heat temperatures. The study estimates gross and net annual generation using the superheater boiler exhaust and clinker cooler exhaust, based on 2023 operational data. • Typical Installation: Common cement plants have preheater kilns with higher exhaust temperatures than Holcim Alpena across a range of production rates. The study estimates gross and net annual generation using the preheater exhaust and clinker cooler exhaust, with a sensitivity analysis to account for the typical range of preheater exhaust temperatures, clinker cooler exhaust temperatures, and clinker production rates. Methodology The evaluation methodology followed a measurement and verification (M&V) strategy based on the International Performance Measurement and Verification Protocol Option B through comprehensive measurements and analyses of the affected systems. Evaluation data was collected from March 9 to March 11, 2024, the test period of the pilot TEG system. During the test period, in coordination with the ITV team, the ATS team adjusted system operations to capture the range of variability expected for each of the variables pertinent to performance of the system. The methodology consisted of two parts: evaluating the performance of the pilot unit's TEG system and estimating the annual TEG impact in terms of gross and net power based on a given waste heat profile. First, the evaluation of the thermoelectric generation performance of the pilot unit relative to the claims was performed by analyzing the collected test data. Gross power of the pilot TEG system was directly measured. Net power was determined by deducting the measured parasitic power from the gross power. The gross power generation was compared to heat transferred to the system by the working fluid (which was heated by steam generated from the kiln waste heat) to calculate the thermal efficiency achieved by the system. Performance of individual semiconductor cartridges within the pilot array was also assessed in terms of measured gross cartridge power and calculated cartridge thermal efficiency. The second part of the evaluation estimated the annual TEG impacts in terms of gross power and net power (calculated from the difference between gross power and parasitic power). This analysis comprised development of mathematical regression models for gross power and parasitic power, with assessment of each model’s goodness-of-fit characteristics to ensure satisfaction of statistical requirements. The models predicted the gross power generation, the parasitic load based on the temperature difference between the hot working fluid and the cold-side fluid (cold water from Lake Huron) entering the system, the volumetric flow rate of the cold-side fluid at the inlet, and the volumetric flow rate of the hot working fluid at the inlet. The annual impact analysis considered a theoretical commercial-scale system sized to capture the available waste heat at a cement plant, consisting of linked pilot-scale units that receive heat from a theoretical gas-to-working-fluid heat exchanger. To estimate annual impacts at the Alpena plant, the gross power and parasitic power regression models were applied to the arrays in the theoretical commercial-scale system. The heat supplied to the unit was calculated based on the kiln run time, annual production, kiln exhaust waste heat, and clinker cooler waste heat derived from 2023 Holcim Alpena kiln operational data. Net power impacts were calculated by deducting the resulting parasitic power from the estimated gross power. Inputs for the model were generated from a combination of hourly data, assumed design considerations for TEG system scale-up from the pilot-scale unit, and assumptions regarding TEG system operations. This analysis was then used as the basis for estimating annual impacts of typical TEG installation at cement plants, by applying sensitivity analyses to key kiln operational characteristics including kiln preheater exhaust temperatures, cooler clinker exhaust temperatures, and plant daily production rates across a range of expected values. Project Results/Findings Table ES- 2 and Table ES- 2 provide a summary of the operating conditions and evaluation results compared to the stated claims from the technology provider. Key takeaways include: • Gross Power: The peak gross power achieved during the testing period was 10.0 kW, compared to the 11.5 kW expected for 573 active cartridges. The claimed gross power was associated with a target hot side of 240°C; however, the system only received a maximum hot-side mean plate temperature of 212°C during the testing period. • Net Power: The pilot-scale unit exceeded the claims for net power, achieving a peak of 7.7 kW net compared to a claim of 7.5 kW. One factor contributing to the higher achieved net power is the relatively high water pressure available through Lake Huron. The pilot TEG system did not require cold-side pumps during the test, whereas most installations would. This reduced the parasitic loads on the system, ultimately contributing to higher net power relative to the gross power. • Thermal Efficiency: The pilot-scale unit outperformed the claimed efficiency, achieving a peak system efficiency of 5.0% thermal efficiency compared to the stated 4.8%. • Cartridge Performance: To compare cartridge performance against claims, the study focused on the third day of testing, which aimed for conditions closest to the design specifications, with a hot side of 240°C and cold-side exit temperature of 6.4°–30°C. On this day, the mean gross power observed in the cartridges within the TEG array was 18.1 W/cartridge, and the peak performance was 34.7 W/cartridge. The estimated mean cartridge efficiency was 5.2%, and the estimated efficiency at peak gross cartridge power was 10%. The regression models developed for gross power generation and parasitic loads were used to estimate the generation impact for given heat input to the TEG from the working fluid (captured from the waste heat) and from the cold loop (Lake Huron) on an hourly basis for a year of operation. Based on this analysis, installation of a commercial-scale TEG system at the Holcim cement plant in Alpena, Michigan, with a waste heat exchanger of 0.85 effectiveness, would generate up to 391 kW of net power, translating to between 920,000 and 1,800,000 kilowatt-hours (kWh) in net electricity per year. Based on typical grid emissions for Alpena, this would avoid estimated net emissions by 752 metric tons of CO₂ annually.² The sensitivity analysis estimated that typical TEG system installations at cement plants could generate an average of 56–1,040 kW of net power, or between 488,000 and 9,110,000 kWh of net energy. This generation potential is most significantly affected by plant production rates and also influenced by preheater and clinker cooler exhaust temperatures. Applying the national average emission rate, typical commercial-scale installations at Holcim plants are projected to avoid between 182 and 3,401 metric tons of CO₂ annually per site. Table ES- 3 shows a summary of the estimated annual impacts.³ While parasitic loads are significant and vary by application, this analysis assumed the use of heating loop pumps and access to Lake Huron as a cold sink. This setup assumed no need for cooling loop pumps due to the available water pressure at the test site. Applications that require cooling towers or additional equipment are likely to experience higher parasitic loads. Therefore, the study’s estimates are most applicable to scenarios with similar parasitic load configurations—namely, access to a high-pressure cold sink. Applicability to other locations may be limited, as differing conditions could necessitate additional pumps and cooling systems, potentially impacting performance significantly.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

Performance Comparison of Machine Learning Models for Ultrasonic Nondestructive Evaluation of Alkali-Silica Reaction in Concrete

Alkali-silica reaction (ASR) causes concrete degradation, leading to cracking, rebar corrosion, and reduced structural integrity, which raises safety concerns. Ultrasonic nondestructive evaluation (NDE) effectively assesses concrete properties and monitors ASR progression. However, its deployment and analysis require specialized expertise and subjective interpretation. As computational power increases, artificial intelligence (AI) and machine learning (ML) algorithms are increasingly being used to automate NDE data analysis across various industries for AI-assisted automation. Regulatory agencies are adapting to this technological shift, prompting a need to evaluate current ML technologies’ capabilities and limitations in assessing concrete material properties and damage. This report presents a comparative analysis of four ML regression models for predicting concrete material damage induced by ASR expansion using long-term ultrasonic data monitoring. The models investigated include linear regression (LR), support vector regression (SVR), shallow neural networks (NN), and deep neural networks (DNN). LR, SVR, and shallow NN models use features extracted from ultrasonic signals, whereas the DNN model processes time-domain ultrasonic signals and frequency spectra directly. The study systematically compared the models’ performance from various perspectives, including model input, prediction performance, and generalization ability. The findings indicate significant variability in model performance, with some ML algorithms achieving very high or very low prediction accuracy depending on the preprocessing and feature engineering (extraction and selection) applied. Key insights include the observation that shallow ML models (LR, SVR, and shallow NNs) require meticulous preprocessing and feature extraction to achieve high accuracy. In contrast, the DNN model, although it bypasses the need for feature engineering, necessitates extensive preprocessing to mitigate noise and computational demands. The SVR model emerged as the top performer among the shallow models, and the DNN model exhibited superior performance on specific datasets but struggled with generalization across specimens from different batches. Additionally, the SVR model is sensitive to temperature variations, whereas the DNN model is robust in this regard. Using recurrent neural networks is recommended for future ASR expansion prediction studies. Recurrent neural networks’ inherent ability to capture temporal dependencies and long-term patterns makes them well suited for analyzing sequential ultrasonic monitoring data. Overall, the results and conclusions of this study could provide insights into the capabilities and effectiveness of ML when applied to ultrasonic NDE data and help identify best practices for using ML for ultrasonic NDE of concrete material properties.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Characterization of Arsenic and Selenium in Coal Fly Ash to Improve Evaluations for Disposal and Reuse Potential (Final Technical Report)

Coal fly ash is a high volume waste material that is discarded in landfills and surface water impoundments across the U.S. and is also widely recycled for a variety of applications. The leaching of potential of contaminants of concern, such as arsenic (As) and selenium (Se), is often the driver of risk assessments for coal ash disposal and reuse. The extent of leachable As and Se depends on several factors related to environmental conditions and fly ash characteristics. Previous studies employed various methods to delineate the concentration, chemical form, and distribution of As and Se in fly ash materials. However, few studies have attempted to directly correlate these properties to mobilization parameters relevant to disposal and reuse. Instead, the coal residuals industries often rely upon standardized leaching protocols that can be laborious or involve hazardous chemicals. The goals of the project were to: 1) Develop and evaluate a characterization protocol that can be used to screen fly ash samples for leachability of As and Se; 2) Characterize As, Se, and associated constituents of fly ash particles at multiple length scales (nanometer to micrometer) to determine if elemental associations differ as a function of the resolution of characterization; and 3) Establish a predictive model for the chemical composition of coal ash produced annually at major U.S. coal fired power facilities on 50-year national coal supply records. For the first objective, we performed leaching experiments with 52 fly ash samples collected from 15 different U.S. power plants and representing coal feedstocks from the three major domestic coal regions. For this work, we assessed the mobilization potential of As and Se in fly ash based on standardized leaching protocols and performed multivariate and lasso regression analyses to explore correlations of leachable As and Se contents with characteristics such as major element contents, loss on ignition (LOI) and pH. The results of regression models indicated that major elements (Fe, Ca, Al) for a wide range of fly ashes can serve as predictor variables for the leaching potential of As, but not for Se. LOI and pH were not important predictive variables in the models. Both regression approaches resulted in relatively strong fits for leachable As (correlation coefficient R 2 = 0.78 for both models) compared to models for leachable Se (R 2 = 0.49). Overall, these results suggest that correlation models combined with on-site elemental analysis with portable analyzers may enable a screening method for leachable As in coal ash. For the second objective, we utilized nanoscale 2-D imaging (30-50 nm spot size) with the Hard X-ray Nanoprobe (HXN) in combination with microprobe X-ray capabilities (~5 µm resolution) to determine As and Se elemental associations in fly ash particles. Speciation of As and Se was also measured at the nano- to microscale with X-ray absorption spectroscopy. The enhanced resolution of HXN showed As and Se that were diffusely located around or comingled with Ca- and Fe-rich particles. The results also showed nanoparticles of Se attached to the surface of fly ash grains. Overall, a comparison of As and Se species across scales highlights the heterogeneity and complexity of chemical associations for these trace elements of concern in coal fly ash. For the final objective, we developed a predictive model for major element composition of coal ash in reserve at disposal sites of major U.S. coal fired power plants. This model was constructed from coal purchase records of 705 power stations from 1973-2022 and was trained on coal ash composition data showing that coal ash elemental composition is strongly associated with the source of feedstock coal. The model showed regional shifts in the major element contents of ash produced by power plants in the last 50 years, particularly for calcium and iron (expressed as %CaO and %Fe 2 O 3 ), as coal-fired power stations changed their source of coal over this time frame. Our approach enables an estimation of coal ash chemical composition that is stored in waste impoundments at individual power stations. Such information can help delineate the regional market potential for material applications that would utilize coal ash harvested from disposal sites across the U.S.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2025

Machine Learning Analysis of Temperature-Strain Relationships for Structural Health Monitoring of Pipes: Self-powered wireless sensor system for health monitoring of liquid-sodium cooled fast reactors

This report presents machine learning (ML) analysis of temperature-strain relationships for structural health monitoring of nuclear reactor stainless steel (SS) pipes with the strain gauge sensor directly printed on the pipe with a 3D conformal aerosol jet printer. We investigate correlations for two sensor pairs installed on the same SS304 pipe: commercial K-type thermocouple with a printed gold strain gauge (TC3-SG3), and commercial K-type thermocouple with commercial Kyowa strain gauge (TC0-SG0). The temperature ranges for the sensor pairs TC0-SG0 and TC3-SG3 are 20.00°C to 266.37°C and 39.95°C to 219.28°C respectively. ML algorithms in this study include Linear Regression (baseline method), Ridge Regression, Lasso Regression, and Gradient Boosting. Performance evaluation metrics include Root Mean Square Error (RMSE), Mean Square Error (MSE), Mean Absolute Error (MAE), R 2 Score, and Explained Variance. Using advanced feature engineering techniques, we extracted 27 temperature-based features and 30 strategic inclusion features. The best performance was obtained with the Gradient Boosting method, which achieves prediction accuracy of R 2 = 0.9999 and RMSE = 7.69 μStrain for TC0-SG0, and R 2 = 0.9998 and RMSE = 18.03 μStrain for TC3-SG3. While the temperature-strain correlations are weaker for the gauge directly printed on the pipe than for the commercial strain gauge, deployment-ready performance exceeding industry standards is achieved for both sensor pairs.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Leveraging design of experiments to build chemometric models for the quantification of uranium (VI) and HNO3 by Raman spectroscopy

Partial least squares regression (PLSR) and support vector regression (SVR) models were optimized for the quantification of U(VI) (10–320 g L −1 ) and HNO 3 (0.6–6 M) by Raman spectroscopy with optimized calibration sets chosen by optimal design of experiments. The designed approach effectively minimized the number of samples in the calibration set for PLSR and SVR by selecting sample concentrations with a quadratic process model, despite complex confounding and covarying spectral features in the spectra. The top PLS2 model resulted in percent root mean square errors of prediction for U(VI), HNO 3 , and NO 3 − of 3.7%, 3.6%, and 2.9%, respectively. PLS1 models performed similarly despite modeling an analyte with a majority linear response (i.e., uranyl symmetric stretch) and another with more covarying vibrational modes (i.e., HNO 3 ). Partial least squares (PLS) model loadings and regression coefficients were evaluated to better understand the relationship between weaker Raman bands and covarying spectral features. Support vector machine models outperformed PLS1 models, resulting in percent root mean square error of prediction values for U(VI) and HNO 3 of 1.5% and 3.1%, respectively. The optimal nonlinear SVR model was trained using a similar number of samples (11) compared with the PLSR model, even though PLS is a linear modeling approach. The generic D-optimal design presented in this work provides a robust statistical framework for selecting training set samples in disparate two-factor systems. This approach reinforces Raman spectroscopy for the quantification of species relevant to the nuclear fuel cycle and provides a robust chemometric modeling approach to bolster online monitoring in challenging process environments.

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

Next-Level Energy Management in Manufacturing: Facility-Level Energy Digital Twin Framework Based on Machine Learning and Automated Data Collection

This research introduces an energy prediction framework at the facility level supported by automated data collection and machine learning models. It investigates whether reducing the prediction time scale allows for applying more complex machine learning techniques and if those techniques improve the prediction accuracy. The primary advantages of this framework lie in its automation of the energy prediction process and its provision of real-time energy data suitable for use in energy dashboards or digital twins. A sitewide dataset was created by combining 15 min energy and daily production data of five shops—assembly, battery, body (electric), body (gas), and paint—from a globally recognized electric vehicle manufacturer. Various machine learning models were evaluated on daily, weekly, and monthly datasets, including, in increasingly complex order: naïve, simple linear regression, net regularized generalized linear regression, principal component regression, k-nearest neighbor, random forest, and Bayesian regularized neural network. Compared to the current state-of-the-art energy consumption prediction for the industrial facility level, this research investigates more complex models and smaller time intervals for higher accuracy. The findings revealed that the more complex monthly models require a minimum of a year and a half of data to operate, while weekly models demand a year of data to achieve improved accuracy. Daily models can operate with only six months of data but exhibit poor performance due to reduced prediction accuracy of production. Key challenges identified include access to reliable, high-quality energy and production data and the initial demand for human labor.

digital twin↗
DOE OSTI2026

Mechanistic Modeling of TEG Dehydrator Emissions in Oil and Gas Industry

This work presents a mechanistic modeling approach for simulating methane emissions from triethylene glycol (TEG) dehydrators used in oil & gas (O&G) operations. The model was developed as a modular component of the Mechanistic Air Emissions Simulator (MAES) tool, incorporating species-specific absorption and emission dynamics through two-level, second-order polynomial regression (PR) models trained on ProMax simulation data: (1) species-level regression models that track the transfer rates of individual gas species within the dehydrator unit streams, and (2) outlet flow stream regression models that predict the fraction of inlet gas distributed among the outlet streams of the dehydrator unit. These behaviors were characterized over a range of glycol circulation ratios, wet gas pressures, and temperatures. The model was validated using root mean square error (RMSE) analysis. The species-level PR achieved low root mean square error (RMSE) values (<0.03) for light hydrocarbon species across all dehydrator components, ranging from 0.0009 for methane to 0.029 for normal pentane. Similarly, the outlet-level PR yielded RMSE values below 0.002 for the dry gas fraction, 0.001 for the flash tank fraction, and 0.002 for the still vent fraction, demonstrating strong agreement between predicted and reference ProMax values. When deployed at field facilities, the model significantly improved MAES-simulated dehydrator emissions, revealing that gas-assisted glycol pump emissions are the dominant contributors to both dehydrator-level and site-level methane emissions under uncontrolled conditions. Further analysis of the 154 dehydrator units reported by operators under the AMI 2024 project showed that 54 units (31%) used gas-driven glycol pumps, of which 6 units (11%) operated with uncontrolled flash tanks, and 22 units (40.7%) were identified as potentially oversized. Of the six dehydrator units with uncontrolled gas-assisted pumps, pump emissions accounted for 90.25% of total dehydrator emissions and 63.10% of total site-level emissions. These findings highlight substantial opportunities for emissions mitigation through equipment upgrades.

MAES↗
DOE OSTI2025

Image-Driven Hybrid Structural Analysis Based on Continuum Point Cloud Method with Boundary Capturing Technique

Conventional approaches for the structural health monitoring of infrastructures often rely on physical sensors or targets attached to structural members, which require considerable preparation, maintenance, and operational effort, including continuous on-site adjustments. This paper presents an image-driven hybrid structural analysis technique that combines digital image processing (DIP) and regression analysis with a continuum point cloud method (CPCM) built on a particle-based strong formulation. Polynomial regressions capture the boundary shape change due to the structural loading and precisely identify the edge and corner coordinates of the deformed structure. The captured edge profiles are transformed into essential boundary conditions. This allows the construction of a strongly formulated boundary value problem (BVP), classified as the Dirichlet problem. Capturing boundary conditions from the digital image is novel, although a similar approach was applied to the point cloud data. It was shown that the CPCM is more efficient in this hybrid simulation framework than the weak-form-based numerical schemes. Unlike the finite element method (FEM), it can avoid aligning boundary nodes with regression points. A three-point bending test of a rubber beam was simulated to validate the developed technique. The simulation results were benchmarked against numerical results by ANSYS and various relevant numerical schemes. The technique can effectively solve the Dirichlet-type BVP, yielding accurate deformation, stress, and strain values across the entire problem domain when employing a linear strain model and increasing the number of CPCM nodes. In addition, comparative analysis with conventional displacement tracking techniques verifies the developed technique’s robustness. The proposed technique effectively circumvents the inherent limitations of traditional monitoring methods resulting from the reliance on physical gauges or target markers so that a robust and non-contact solution for remote structural health monitoring in real-scale infrastructures can be provided, even in unfavorable experimental environments.

Chemistry↗
DOE OSTI2026

Advancing the quantification of aerosol-cloud interactions with the CALIPSO-CloudSat-Aqua/MODIS record

Aerosol-cloud-precipitation interactions are assessed over the non-polar ocean using more than 11 years of combined Aqua-MODIS, CALIPSO-CALIOP, and CloudSat products. The analysis first shows the benefit of incorporating vertically resolved aerosol extinction coefficient (σext) in aerosol-cloud interactions (ACI) assessments, demonstrating that: σext vertically collocated with the cloud layer () correlates best with cloud droplet number concentration (Nd), column-integrated aerosol optical depth (AOD) cannot explain the Nd variability in the extratropics, and the S-shape of the AOD-Nd relationship reported in previous studies is not replicated when using instead of AOD, with a Nd- linearity more consistent with in-situ studies over the ocean. ACI metric, estimated as the log-scale regression between CALIOP and MODIS Nd reveals that the eastern Pacific is the region with the strongest ACI, followed by the Southern Ocean. The susceptibility of clouds to changes in their liquid water path (LWP) and frequency of precipitation followed a 2-step calculation by combining the Nd- regression (ACI) with the regression between these macrophysical variables and Nd. LWP susceptibility is negative (LWP decreases with aerosol loading) and statistically significant over the eastern Pacific, eastern Atlantic, and extratropics. In contrast, vast areas of the tropical and subtropical ocean feature negligible changes in LWP with aerosol. Precipitation frequency susceptibility is negative, but the values are only significant over the coastal eastern Pacific and Atlantic. The findings suggest that previous modeling assessments relying on AOD may need to be revisited by taking advantage of the synergy between passive and active sensors.

Li, Zhujun↗
DOE OSTI2024

A robust approach to Gaussian process implementation

Abstract. Gaussian process (GP) regression is a flexible modeling technique used to predict outputs and to capture uncertainty in the predictions. However, the GP regression process becomes computationally intensive when the training spatial dataset has a large number of observations. To address this challenge, we introduce a scalable GP algorithm, termed MuyGPs, which incorporates nearest-neighbor and leave-one-out cross-validation during training. This approach enables the evaluation of large spatial datasets with state-of-the-art accuracy and speed in certain spatial problems. Despite these advantages, conventional quadratic loss functions used in the MuyGPs optimization, such as root mean squared error (RMSE), are highly influenced by outliers. We explore the behavior of MuyGPs in cases involving outlying observations and, subsequently, develop a robust approach to handle and mitigate their impact. Specifically, we introduce a novel leave-one-out loss function based on the pseudo-Huber function (LOOPH) that effectively accounts for outliers in large spatial datasets within the MuyGPs framework. Our simulation study shows that the LOOPH loss method maintains accuracy despite outlying observations, establishing MuyGPs as a powerful tool for mitigating unusual observation impacts in the large data regime. In the analysis of US ozone data, MuyGPs provides accurate predictions and uncertainty quantification, demonstrating its utility in managing data anomalies. Through these efforts, we advance the understanding of GP regression in spatial contexts.

Mukangango, Juliette↗
DOE OSTI2025

Beyond pinball loss: Quantile methods for calibrated uncertainty quantification

Amongthemanywaysofquantifying uncertainty in a regression setting, specifying the full quantile function is attractive, as quantiles are amenable to interpretation and evaluation. A model that predicts the true conditional quantiles for each input, at all quantile levels, presents a correct and efficient representation of the underlying uncertainty. To achieve this, many current quantile-based methods focus on optimizing the pinball loss. However, this loss restricts the scope of applicable regression models, limits the ability to target many desirable properties (e.g. calibration, sharpness, centered intervals), and may produce poor conditional quantiles. In this work, we develop new quantile methods that address these shortcomings. In particular, we propose methods that can apply to any class of regression model, select an explicit balance between calibration and sharpness, optimize for calibration of centered intervals, and produce more accurate conditional quantiles. We provide a thorough experimental evaluation of our methods, which includes a high dimensional uncertainty quantification task in nuclear fusion.

97 MATHEMATICS AND COMPUTING↗