Search NASA⌕ Search

SEARCH · Search NASA

Results for “error estimation”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

High-Spectral Resolution Lidar (HSRL) Instrument Handbook

High-spectral-resolution lidar (HSRL) systems provide vertical profiles of optical depth, backscatter cross-section, depolarization, and backscatter phase function. All HSRL measurements are calibrated by reference to molecular scattering, which is measured at each point in the lidar profile. Like the Raman lidar (RL), but unlike simple backscatter lidars such as the micropulse lidar (MPL), this enables the HSRL to measure backscatter cross-sections and optical depths without prior assumptions about the scattering properties of the atmosphere. The depolarization observations allow robust discrimination between ice and water clouds. Rigorous error estimates can be computed for all measurements. A very narrow angular field of view reduces multiple scattering contributions. The small field of view, coupled with a narrow optical bandwidth, nearly eliminates noise due to scattered sunlight. The laser transmitter is a diode-pumped, frequency-doubled Nd:YAG laser. Narrow-band, single-frequency operation is provided by injection seeding with a single-frequency, cw-diode-pumped diode laser. The main laser cavity is maintained in resonance with the seed laser by adjusting the cavity length to minimize the time between the Q-switch trigger and the emission of the laser pulse. The emission wavelength is tuned via temperature control of the seed laser crystal and is locked to line #1109 of the iodine absorption spectra. Locking is accomplished by minimizing the transmission through a 2-cm-long iodine absorption cell. Use of a high-repetition-rate laser and expansion of the transmitted beam through a 400-mm telescope reduces the transmitted energy density to eye-safe levels. It is possible to look directly into the output beam without hazard. The receiver and transmitter use the same afocal telescope, simplifying the maintenance of stable alignment of the transmitter and receiver although the angular FOV is only 100 μrad. The small FOV and the 4-kHz repetition rate also limit the near-field signal strength, making it possible to record continuous profiles that start at an altitude of ~100 m and extend to 30 km using photon counting detectors. The small FOV also suppresses multiple scattering contributions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Measuring Orbit Responses with Oscillating Trajectories in the Fermilab Linac

Recording changes in beam transverse positions re-ported by Beam Position Monitors (BPMs) in response to a beam deflection by an upstream dipole corrector (orbit response) is a powerful tool for analysis of accelerator optics and assisting with machine tuning. In the Fermilab Linac, orbit responses are recorded by oscillating the currents of up to 19 correctors, providing faster, drift-resistant measurements through frequency-domain analysis. This report describes the technique, including error estimations and consistency checks and shows an example of the measurements.

Shemyakin, Alexander V. [Fermilab] (ORCID:00000001↗
DOE OSTI2026

Development and implementation of high-throughput proteomic and metabolomics assays by using advanced chromatographic and mass spectrometric systems (CRADA Final Report)

The mission of this CRADA with Agilent was to couple powerful MS platforms (QQQ, IM-QTOFMS) with Agilent’s novel Ultra-High-Performance Liquid Chromatography (UHPLC) fast metabolomic workflows and perform ABF Machine Learning (ML) to generated datasets. Agilent transferred UHPLC methods to PNNL and LBNL and methods were implemented and demonstrated in both labs, achieving total acquisition times of < 10 min. Metabolites analyzed using Agilent’s shared methods included metabolites from central carbon metabolism, common across hosts, and metabolites unique to engineered strains. Standards were acquired in an UHPLC-Drift Tube Ion Mobility Mass Spectrometer (DTIMS) system for the first time within the context of ABF and methods were optimized based on Agilent’s protocols. Samples from ABF hosts Pseudomonas putida, Aspergillus pseudoterreus, Aspergillus niger and Rhodosporidium toruloides were analyzed using the UHPLC-DTIMS platform for a total of 276 runs. A data analysis workflow compatible with the Experimental Data Depot (EDD) and completely shareable was developed for the acquired UHPLC-DTIMS data. Samples were analyzed using a Data Independent Acquisition Approach (DIA), which for most of the standards provided more transitions therefore increasing detection confidence. Using the data acquired by PNNL, LBNL, and Agilent’s specifications from previous ML projects, SNL applied an ensemble ML strategy to pick the best performing model for automated LC-method selection. Finally, with the contribution of the participant labs and Agilent, SNL developed an Automated Method Selection (AMS) software tool to predict the best liquid chromatography method for analysis of any new molecules of interest. Samples with novel pathways and new metabolite targets of interest are generated at a high pace in the ABF. Overall, the project advanced rapid metabolomics by combining liquid chromatography, ion mobility spectrometry, and data-independent mass spectrometry with machine learning. This multidimensional approach uses retention time, collision cross-section, precursor mass, and fragment-ion information to distinguish chemically similar metabolites that can be difficult to resolve using conventional liquid- or gas-chromatography methods. The resulting workflow also provided automated metabolite-identification error estimates, addressing a recognized need for statistical confidence measures in metabolomics.

Petzold, Christopher [Lawrence Berkeley National L↗
DOE OSTI2025

Increased Ly α Opacity in White Dwarf Photospheres from Transient H − Resonances

The structure and emergent flux of hydrogen atmosphere white dwarfs depend on the opacity of the Lyα and Lyβ spectral lines. The opacity here is set by the strength and broadening of these lines; the latter is dictated by the far line wing, which is in the “quasi-static” limit of electron broadening, placing it in the incomplete collision regime, and describes the transient parts of electron and ion collisions. These transient stages of the collision form resonances: In the case of ions, they manifest as molecular resonances, while for electrons they are H − resonances, both of which can only be captured quantum-mechanically. Quantum-mechanical calculations have historically preserved only a handful of broadening terms that are most important near the center of the line. However, in the wings of the line, the previously neglected terms that describe the transient stages of the collision need to be included. This requirement arises because, in the line wings, the broadening from the 1s ground state, which is generally assumed to be extremely small compared to the broadening of the upper state, is no longer negligible within a quantum-mechanical model that takes into account exchange interactions. The inclusion of all the transient terms results in asymmetries and extra broadening. The increased broadening of Lyα increases the opacity at the energy where most of the flux leaves the star. The broader Lyα lines also impact the visible flux, raising it by an amount that exceeds previously estimated errors.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Efficient screening of rare large pit anomalies on polished surfaces using a minimalist sampling scheme

Lawrence Livermore National Laboratory (LLNL) has made significant strides in generating clean energy through its inertial confinement fusion (ICF) experiments. These experiments rely on high-density carbon (HDC) coated shells to encapsulate the fusion fuel. The success of these experiments is heavily dependent on the surface quality of these shells, as even minor imperfections, such as deep pits, can negatively impact fusion yield. Ensuring the required smoothness involves an extensive surface-finishing process that spans approximately 20 stages, making it both time-intensive and resource-demanding. A critical challenge in this process is the need for high-resolution scans to detect rare deep pits, which can be costly and impractical if performed on every shell. This highlights the necessity of developing more efficient scanning methods to optimize time and cost without compromising accuracy. To address these challenges, we introduce a novel approach that employs the multivariate Dvoretzky–Kiefer–Wolfowitz (DKW) inequality to provide a probabilistic upper bound on the error in estimating pit distribution characteristics via a Kernel Density Estimator (KDE). This error bound enables efficient and reliable estimation of pit distribution characteristics at a specified statistical confidence level using a minimal number of surface scans. The integrated DKW-KDE approach was validated through surface-finishing experiments across two batches of HDC-coated shells, demonstrating consistent and robust performance across multiple stages of the surface-finishing experiments. The validation studies suggest that the integrated DKW-KDE approach achieves comparable accuracy in estimating the risk of deleterious large pits with six scans, thus conserving time and resources. Further evaluations show that performance remains consistent across batches and over multiple polishing stages. In conclusion, based on these findings, one can leverage the minimal-scan insights to strategically improve the bottleneck inspection process, thus enhancing the productivity and quality of shell polishing and similar challenging manufacturing processes.

Inertial confinement fusion↗
DOE OSTI2024

A Nonintrusive Optical Approach to Characterize Heliostats in Utility-Scale Power Tower Plants: Camera Position Sensitivity Analysis

Optics plays a major role in the effectiveness of concentrating solar power (CSP) technologies. The nonintrusive optical (NIO) approach developed by the National Renewable Energy Laboratory uses uncrewed aircraft system (UAS)-based imaging to survey heliostats in a commercial-scale power tower CSP plant and characterize their optical errors. The image processing algorithm uses photogrammetry to calculate the camera position for each image frame, and the accuracy of the estimated optical errors is highly sensitive to the calculated camera position accuracy. In this study, we simulate a series of case studies in python to examine the impact of different parameters of the sensitivity of the camera calculation, including the number of facet corners used as control points for the photogrammetric calculation, precision error in the detected pixel locations of the facet corners in the image, and precision error of the canting and mounting positions of the facets of the heliostat. The case studies consider heliostat geometry based on three commercial designs to serve as representative examples of different possible sizes of heliostats that the NIO method could be applied to. The results show that increasing the number of control points can improve accuracy for heliostats with many facets, pixel precision has a significantly larger impact on camera calculation accuracy than facet canting and mounting errors, and the camera distance and focal length must be chosen to ensure adequate pixel accuracy on the heliostat surface depending on the size of heliostat. In conclusion, based on the results, recommendations for the allowable values of each parameter are provided to achieve the required NIO optical error estimation accuracy depending on the size of heliostat.

14 SOLAR ENERGY↗
DOE OSTI2024

A Novel Noise-Aware Classical Optimizer for Variational Quantum Algorithms

A key component of variational quantum algorithms (VQAs) is the choice of classical optimizer employed to update the parameterization of an ansatz. It is well recognized that quantum algorithms will, for the foreseeable future, necessarily be run on noisy devices with limited fidelities. Thus, the evaluation of an objective function (e.g., the guiding function in the quantum approximate optimization algorithm (QAOA) or the expectation of the electronic Hamiltonian in variational quantum eigensolver (VQE)) required by a classical optimizer is subject not only to stochastic error from estimating an expected value but also to error resulting from intermittent hardware noise. Model-based derivative-free optimization methods have emerged as popular choices of a classical optimizer in the noisy VQA setting, based on empirical studies. However, these optimization methods were not explicitly designed with the consideration of noise. In this work we adapt recent developments from the “noise-aware numerical optimization” literature to these commonly used derivative-free model-based methods. We introduce the key defining characteristics of these novel noise-aware derivative-free model-based methods that separate them from standard model-based methods. In conclusion, we study an implementation of such noise-aware derivative-free model-based methods and compare its performance on demonstrative VQA simulations to classical solvers packaged in scikit-quant.

classical optimizers↗
DOE OSTI2025

Dephasing and error dynamics affecting a singlet-triplet qubit during coherent spin shuttling

Quantum information transport over micron to millimeter scale distances is critical for the operation of practical quantum processors based on spin qubits. One method of achieving a long-range interaction is by coherent electron spin shuttling through an array of silicon quantum dots. In order to execute many shuttling operations with high fidelity, it is essential to understand the dynamics of qubit dephasing and relaxation during the shuttling process in order to mitigate them. However, errors arising after many repeated shuttles are not yet well documented. Here, we probe decay dynamics contributing to dephasing and relaxation of a singlet-triplet qubit during coherent spin shuttling over many N repeated shuttle operations, in a small external magnetic field B 0 ≈ 0−10 mT, and in the absence of a micromagnet. We find that losses are dominated by magnetic dephasing, most visible for small N < 10 3 . However, incoherent spin-flip type shuttle errors become evident for large N > 10 3 . Additionally, we estimate shuttle error rates below 10 −4 out to at least N = 10 3 , representing an encouraging figure for future implementations of spin shuttling to entangle distant qubits.

quantum dots↗
DOE OSTI2024

Blue Canyon Dome Emplacement Artifacts

The Blue Canyon Dome (BCD) experiments are a series of intermediate scale subsurface chemical explosions designed to explore the geophysical relationship between scaled yields and formation damage. Phase 2 of this experiment involved explosions designed to evaluate the effect of damage, demonstrating the effect of geologic features and previous damage on the near field seismic signatures, not just in amplitude, but also in spectral content. Based on the emplacement effect on the seismic signature, the moment-based yield estimates showed error between 10s-1000s of percent. This report details an alternative approach to emplacement-based filtering techniques of the seismic signatures to improve yield estimates for explosions in damaged rock. The results will be used to develop an improved inverse yield model that leverages the information in the seismic signature for each source based on emplacement type.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2024

Estimating the Contributions to Human Error Probability from the Convolution of the Distribution of Time Available and Time Required

As part of their duties, Human Reliability Analysis must often evaluate if crews in nuclear power plants (NPPs) can complete tasks associated with a human-failure event within time limits. For example, the time required in NPP scenarios is determined by systematic and structured walkthroughs, feasibility studies, recorded times from training exercises, and interviews with experienced operators and experts. Typically, a point estimate is derived for the estimate (mean, maximum, or 95th percentile of time required). Using point-estimate values can mask the risk associated with variability among crews, plant conditions and set-up, environmental conditions, and other impact factors under which these actions are executed. While point estimates for time required and time available have served the industry well, without considering the uncertainty they could lead to biased understanding about the risk. The Integrated Human Event Analysis System - General Methodology (IDHEAS-G) model (developed by the US Nuclear Regulatory Commission, NRC) for human error probability calculates human error probability by summing two probabilities: insufficient time and cognitive error. As such, the model takes a more holistic approach by considering the full distributions for time required and time available to calculate the human error probability because the time available to complete the task is insufficient. In this study, we expand on the work of the NRC and discuss methods for estimating these time considerations. For example, for the time required, the impact of Performance Influencing Factors (PIFs) on the distribution was divided into impacts that are aleatory in nature, such as crew-to-crew variability, and those that are epistemic (i.e., the PIFs). Starting with the factors that introduce aleatory uncertainty, a first-order distribution was developed from a large set of time required (i.e., NPP task completion times) data for the range of operator actions that occur in the NPP control room under simulated accident conditions. The first-order distribution can then be adjusted to account for epistemic uncertainty using research associated with the impact of applicable PIFs on the time required. We also develop guidance for analysts to address the probability distributions for the time available. The guidance we developed on how to estimate time required and time available distributions is based on the identification of pertinent research and data, data analyses, and expert knowledge elicitation.

human error probability, human performance, time e↗
DOE OSTI2026

ZTF-SEDm Type Ia supernova sample for Twins Embedding spectrophotometric standardization

Aims. This paper has two aims: the first aim is to build a large homogeneous spectrophotometric sample of Type Ia supernovae (SNe Ia) from the second data release of the Zwicky Transient Facility (ZTF DR2). We used the spectrum sample from the low-resolution ( R ∼ 100) SEDmachine (SEDm) Integral Field Spectrograph (IFS) that gathered 3069 spectra. This is one of the largest samples of such collections that can attempt to reproduce the Twins Embedding (TE) spectrophotometric standardization method. This is our second objective. The method was developed based on high-quality spectra from 200 SNe Ia of the Nearby Supernova factory (SNfactory) and led to an exceptionally low value of 0.073 mag for the intrinsic scatter. Methods. As the SEDm is not designed as a spectrophotometric instrument, we first improved the flux-calibration accuracy of the SN Ia spectrum sample using the ZTF photometric data, which were calibrated at the percent level. We corrected the spectra for second-order polynomials, fitted by comparing the synthetic photometry in the ZTF g , r , i filters with the light-curve (LC) data. We then applied the three steps of the TE parameterization to a subset of 783 ZTF SN spectra near maximum light while comparing results from SNfactory and ZTF. We finally analyzed the standardization methods based on the TE parameters. Results. The precision of the phase-correction model, which is the first step of the TE, is estimated at 0.01 mag in g band based on ZTF data. Despite the challenge posed by the spectrum-extraction pipeline associated with the SEDm (flux calibration, leftover host signal, low signal-to-noise ratio, and low resolution), we applied a first standardization in color based on the second step of the TE, called read between the lines (RBTL), to the ZTF sample. We reached a Hubble residual scatter of 0.153 mag, all in normalized median absolute deviation, which is to be compared to the ∼0.11 mag obtained with the SNfactory data. The SALT color and stretch standardization reaches a scatter of 0.164 mag for the same ZTF SN Ia sample, and its host steps are ∼0.1 mag and zero for RBTL. When considering the scatter due to the redshift error and flux calibration error, we estimated a RBTL scatter of ∼0.129 mag for this ZTF sample as an upper limit because we identified an additional contribution from a systematic error in color. We tested the standardization based on the nonlinear TE parameters, and, as expected from the low spectrum quality, it did not improve the overall dispersion. Conclusions. We release 1897 flux calibrated spectra of 1607 SNe Ia with an estimated photometric accuracy of 0.07 mag. We further demonstrate that some amount of spectrophotometric SN Ia standardization can be done with limited-quality spectra. The RBTL standardization is more efficient than that of SALT with one parameter less, and the resulting host steps are consistent with zero. This makes it less prone to astrophysical bias. For future spectroscopic surveys, targeting the extraction pipeline for a thorough flux calibration and good signal-to-noise ratio would enable us to compute the full TE standardization, which would further reduce the scatter in the distance estimate.

Ganot, C↗
DOE OSTI2025

Measuring the thermal conductivity of hydrogels with a bidirectional 3w method

Hydrogels are soft, water-absorbing polymer materials with diverse applications in biomedicine and agriculture. Recently, hydrogels have been proposed to encapsulate water-soluble phase change materials which store energy in their latent heat of solidification. In these applications, the thermal conductivity of these materials affects their performance. Few methods exist for measuring the thermal conductivity of small quantities of hydrogels. Here, we describe an implementation of the bidirectional 3w technique to measure the thermal conductivity of hydrogels with particular attention to their moisture content. Our implementation of the technique can probe sample volumes as little as ~20 mL and yields the thermal conductivity without requiring fitting of additional thermal parameters. We numerically simulate 3w sensor designs with frequency-domain 3-D models to quantify and reduce errors introduced by the choice of substrate and insulation layer thickness. Frequencies in the ~1−20 Hz range yield less error for the materials considered here. We verify our setup with measurements on water and report values for polyacrylamide and poly(2-acrylamido-2-methylpropane sulfonic acid) (PAMPS) hydrogels. Our swollen hydrogels exhibited thermal conductivities nearly equivalent to water, 0.6 W m-1 K-1, and we estimate thermal conductivities of 0.43 and 0.42 W m-1 K-1 for neat polyacrylamide and PAMPS, respectively. Finally, we estimate an error of ±7%, consistent with other 3ω methods, with the largest error coming from the sensor calibration. We find our implementation of the bidirectional 3w method gives reasonable results and can be employed for prototyping soft materials relevant for thermal storage.

3-omega, thermal conductivity, hydrogel, moisture ↗
DOE OSTI2025

An Open-Source Python Package for CFD Solution Verification

Informed decision-making using computational fluid dynamics (CFD) results requires quantifying the errors and uncertainties of a simulation. Verification, validation, and uncertainty quantification (VVUQ) methods were developed to address this need and have matured. However, these VVUQ analyses are often non-trivial and require CFD analysts and practitioners to have specific skill sets. This has led to the uneven adoption of VVUQ analyses, in part, based on the availability of software tools to aid CFD analysts and practitioners. Solution verification, a procedure to evaluate the accuracy of a simulation by estimating potential errors arising from the computational model and computing the uncertainties without comparing to results from a physical system, is one of the lagging VVUQ analyses as the absence of software has forced CFD analysts and practitioners to develop their own codes or piece together incomplete software from across the internet. This work presents an opensource Python package, CFDverify, to lower the barrier of entry and fill in the technological gap in solution verification. CFDverify also provides a streamlined framework to remove some potential errors in post-processing CFD results. The hope is that CFDverify can improve the quality and quantity of CFD solution verification in scientific and research studies and attract interest in developing a communal tool. This paper describes the design, features, and an example use of CFDverify.

Weinmeister, Justin [ORNL] (ORCID:0000000160090237↗
DOE OSTI2024

Taylor-Expansion-Based Robust Power Flow in Unbalanced Distribution Systems: A Hybrid Data-Aided Method

Traditional power flow methods often adopt certain assumptions designed for passive balanced distribution systems, thus lacking practicality for unbalanced operation. moreover, their computation accuracy and efficiency are heavily subject to unknown errors and bad data in measurements or prediction data of distributed energy resources (ders). to address these issues, this paper proposes a hybrid data-aided robust power flow algorithm in unbalanced distribution systems, which combines taylor series expansion knowledge with a data-driven regression technique. the proposed method initiates a linearization power flow model to derive an explicitly analytical solution by modified taylor expansion. to mitigate the approximation loss that surges due to the der integration and bad data, we further develop a data-aided robust support vector regression approach to estimate the errors efficiently. comparative analysis in the 13-bus and 123-bus ieee unbalanced feeders shows that the proposed hybrid algorithm achieves superior computational efficiency, with guaranteed accuracy and robustness against outliers.

data-driven↗
DOE OSTI2025

Atacama Cosmology Telescope: DR6 gravitational lensing and SDSS BOSS cross-correlation measurement and constraints on gravity with the 𝐸 𝐺 statistic

We derive new constraints on the 𝐸 𝐺 statistic as a test of gravity, combining the cosmic microwave background (CMB) lensing map estimated from Data Release 6 (DR6) of the Atacama Cosmology Telescope with Sloan Digital Sky Survey III Baryon Oscillation Spectroscopic Survey (SDSS BOSS) CMASS and LOWZ galaxy data. We develop an analysis pipeline to measure the cross-correlation between CMB lensing maps and galaxy data, following a blinding policy and testing the approach through null and consistency checks. By testing the equivalence of the spatial and temporal gravitational potentials, the 𝐸 𝐺 statistic can distinguish Λ⁢ CDM from alternative models of gravity. We find 𝐸 𝐺 ⁡(𝑧 eff = 0.555) = 0.3⁢1$^{+0.06}_{−0.05}$ for Atacama Cosmology Telescope (ACT) and CMASS data at 68.28% confidence level, and 𝐸 𝐺 ⁡(𝑧 eff = 0.316) = 0.4⁢9$^{+0.14}_{−0.11}$ for the ACT and LOWZ. Systematic errors are estimated to be 3% and 4%, respectively. Including CMB lensing information from Planck PR4 results in 𝐸 𝐺 ⁡(𝑧 eff = 0.555) = 0.3⁢4$^{+0.05}_{−0.05}$ with CMASS and 𝐸 𝐺 ⁡(𝑧 eff = 0.316) = 0.4⁢3$^{+0.11}_{−0.09}$ with LOWZ. These are consistent with predictions for the Λ⁢ CDM model that best fits the Planck CMB anisotropy and SDSS BOSS baryon acoustic oscillations (BAO), where 𝐸$^{GR}_{𝐺⁡}$(𝑧 eff =0.555) =0.401 ± 0.005 for CMB lensing combined with CMASS and 𝐸$^{GR}_{𝐺}$⁡(𝑧 eff = 0.316) = 0.452 ± 0.005 combined with LOWZ. We also find 𝐸 𝐺 to be scale independent, with probability to exceed >5%, as predicted by general relativity. The methods developed in this work are also applicable to improved future analyses with upcoming spectroscopic galaxy samples and CMB lensing measurements.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Assessing the Performance of Exchange‐Correlation Density Functionals in Describing the Iron‐Catalyzed Ammonia Synthesis System

Density functional theory (DFT) has been widely employed for elucidating mechanistic aspects of heterogeneous catalysis. However, the accuracy of DFT calculations relies heavily on selecting an exchange-correlation (XC) functional that correctly describes the electronic structure of materials involved in the reactions. This study assesses the accuracy of several XC density functionals for modeling the iron-catalyzed ammonia synthesis system. In the assessment of functional accuracy, experimental references are compared to DFT-calculated values for the formation energy of gas-phase ammonia and nitrogen, bulk Fe/Fe-nitride (γ′-Fe 4 N) lattice constants and cohesive/formation energies, and nitrogen and ammonia binding energies on Fe(100), Fe(111), Fe(110), and γ′-Fe 4 N(111). It is observed that the experimental value for each of these descriptors is accurately modeled by at least one functional. RPBE alone provides reliable estimates for both the lattice constant and cohesive energy of Fe and γ′-Fe 4 N. Temperature-programmed desorption experiments led to estimates for N and NH 3 adsorption across several Fe-based facets that are best captured by RPBE. These results highlight the importance of choosing an appropriate XC functional that accurately describes experimental systems and offer insights into effectively modeling the interactions between nitrogen and ammonia on Fe-based surfaces.

adsorption↗
DOE OSTI2024

On the Representativity of Electrode Microstructure Parameters and Their Electrochemical Response for Lithium Ion Batteries

Lithium-ion battery electrochemical models require an accurate description of the electrodes microstructures to be predictive that can be achieved through nanoscale imaging. Such observations are however limited by their field of view (FOV), as they provide only a subset of the whole electrode volume that does not necessarily represent the whole electrode microstructure heterogeneity, and therefore can bias the microstructure analysis. A microstructure scale electrochemical model was used to investigate lithium plating onset, material non-uniform utilization, and in-plane heterogeneities for an NMC-graphite full cell. To evaluate the representativeness, and thus relevance, of these model predictions, a coupled representativity analysis has been performed on the microstructure parameters and, in a novel way, on the full cell electrochemical response. Electrode microstructure parameters representativeness has been first quantified using the representative volume element (RVE) methodology. The RVE major flaw is that ultimately it can only conclude if a FOV contains representative subvolumes of the FOV, but not if the FOV itself is representative of the electrode volume. Analysis can conclude negatively ('FOV is not representative'), but not positively ('FOV is representative'). One major contribution of this work was to quantify the convergence of the RVE size with the FOV, to actually investigate the FOV representativeness and thus partly remedy this intrinsic limitation. The analysis determined that performing a standard RVE calculation, without exploring its FOV convergence, is likely to strongly underestimate the actual RVE size. The new RVE methodology has been automated in the NREL open-source Microstructure Analysis Toolbox (MATBOX) and is available to the battery community. Representativeness of microstructure parameters is however only an intermediate step, as the end-results of an electrochemical model are performances predictions. Indeed, what is the practical consequence of a given deviation for a microstructure parameter? The microstructure parameter deviation propagations to the 3D microstructure scale electrochemical response have been then quantified for different charge rates. This defines a threshold for the microstructure parameters FOV for a desired maximum deviation of the electrochemical response. Such deviation propagation analysis is analogous to error propagation analysis and is necessary to determine the relevance of microstructure scale model predictions for macroscale predictions. Electrochemical model shows cell representative section areas are increasing with C-rate, due to higher in-plane heterogeneities, indicating larger FOVs are required specifically for fast charge modeling. Therefore, we introduced the novel concept of electrochemical RVE (eRVE) that is a function of the operating conditions (thus defined as a dynamic RVE), with an increasing dependence with the C-rate. Representativity analysis of the investigated cell determined a FOV of 144.4 x 54.4 m2 is large enough to establish a convergence on the representative section areas for low to intermediate C-rate (=2.5C), but not large enough to conclude for higher rates. This work aims to emphasize the importance of representativity analysis for LIB electrode microstructures, as it is required to estimate the error, and thus the relevance, of microstructure parameters intended to be used in macroscale models. The methodology and results can help researchers to select the relevant imaging and associated FOV required to provide accurate enough microstructure parameters.

ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2024

PINN surrogate of Li-ion battery models for parameter inference, Part II: Regularization and application of the pseudo-2D model

Bayesian parameter inference is useful to improve Li-ion battery diagnostics and can help formulate battery aging models. However, it is computationally intensive and cannot be easily repeated for multiple cycles, multiple operating conditions, or multiple replicate cells. To reduce the computational cost of Bayesian calibration, numerical solvers for physics-based models can be replaced with faster surrogates. A physics-informed neural network (PINN) is developed as a surrogate for the pseudo-2D (P2D) battery model calibration. For the P2D surrogate, additional training regularization was needed as compared to the PINN single-particle model (SPM) developed in Part I. Both the PINN SPM and P2D surrogate models are exercised for parameter inference and compared to data obtained from a direct numerical solution of the governing equations. A parameter inference study highlights the ability to use these PINNs to calibrate scaling parameters for the cathode Li diffusion and the anode exchange current density. By realizing computational speed-ups of ~2250x for the P2D model, as compared to using standard integrating methods, the PINN surrogates enable rapid state-of-health diagnostics. Finally, in the low-data availability scenario, the testing error was estimated to ~2 mV for the SPM surrogate and ~10 mV for the P2D surrogate which could be mitigated with additional data.

25 ENERGY STORAGE↗