Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “electrochemical modeling”

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

Modeling the Nucleation and Growth of Lead Sulfate Particles on Lead Electrodes

Lead-acid batteries (LABs) play a pivotal role in the energy storage sector with applications spanning from starting-lighting-ignition batteries to grid energy storage. Passivation of lead negative electrodes by PbSO 4 particles is a fundamental mechanism limiting the performance of LAB. In this regard, an electrochemical model is developed that simulates the nucleation and growth (N&G) dynamics of PbSO 4 particles on a flat lead electrode, responsible for its passivation. The model considers the electrochemical reactions between lead electrode and sulfuric acid, N&G of PbSO 4 particles, passivation of the lead surface, and the ternary transport of PbSO 4 (aq), bisulfate, and protons in H 2 SO 4 electrolyte. The model is validated with a dataset of cyclic voltammetry (CV) responses collected at several scan rates and H 2 SO 4 concentrations. The model shows remarkable qualitative and quantitative agreement with the experimental data including CV peak features, discharge capacity, and particle size. The model was employed to explore key N&G quantities, such as supersaturation, nucleation rate, particle count, growth rate, particle size, and surface coverage, and to examine how their interactions influence electrode utilization. Parametric studies were also conducted to evaluate how scan rates and acid concentrations influence the previously mentioned N&G quantities and, subsequently, the utilization of the electrode.

25 ENERGY STORAGE↗
DOE OSTI2023

Simulation of Impedance Changes with Aging in Lithium Titanate-based Cells Using Physics-Based Dimensionless Modeling

Quantifying aging effects in lithium-ion cells with chemistries that have a flat open circuit potential is challenging. We implement a physics-based electrochemical model to track changes in the electrochemical impedance response of lithium titanate-based cells. Frequency domain equations of a pseudo two-dimensional model are made dimensionless, and the corresponding non-dimensional parameters are estimated using a Levenberg-Marquardt routine. The model weighs the relative contributions of changes in diffusion, ionic conduction within the electrolyte phase against solid phase electronic conduction towards cell aging. Solid-phase diffusion, charge transfer resistance and double layer capacitance at the solid-liquid interface are accounted for in the particle impedance. The estimation routine tracks dimensionless parameters using accelerated cycling data from full cells over 1000 cycles. The model can be deployed within a short time for state estimation using physics-based models without requiring prior knowledge of the battery chemistry, format, or capacity.

25 ENERGY STORAGE↗
DOE OSTI2024

FY24: LEMMs: Long-term, Electrochemical Materials degradation Models

FY24 progress for LEMMs: Long-term, Electrochemical Materials degradation Models LDRD project are presented. Significant progress toward creating validated corrosion models was made in this FY paving success in future FY’s for success with battery model development and validation. Specifically, models for various forms of corrosion were probed for sensitivities showing a significant influence of reactive transport parameters for hydroxide species on the run time of the models. Further, cryo-genic time-of-flight secondary ion mass spectroscopy was utilized to map ions in a frozen corrosion droplet showing the precipitates/precipitate species near corroding locations. Additionally, a realistic corrosion droplet was created that accounted for evaporation and condensation coupled with corrosion. This pushes the boundary of corrosion modeling and will be validated in future FY’s.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Modeling Losses in a Three Electrode System Towards Fast Charge Control

In this work we demonstrate the applicability of a versatile separator-reference mounted electrode in a three-electrode setup to accurately capture the primary current pathway in the battery cell during operation, calibrate a porous electrode model to the individual anode and cathode potential signals, and assign the various resistances to electrochemical phenomena in the battery cell. This calibrated electrochemical model is validated using continuous rate discharges associated with highway driving scenarios in an electric vehicle, and in turn utilized to predict the local anode potential and proximity to lithium plating onset. Finally, we demonstrate the strategy associated with utilization of the model to estimate constant anode potential charging during fast charge scenarios at various rates and starting conditions as a future look to fast charge calibration development and controls.

Garrick, Taylor R. (ORCID:0000000322518129)↗
DOE OSTI2021

Modeling Activities Related to Waste Form Degradation: Progress Report

This report represents the milestone deliverable M4SF-21SN010309021 “Modeling Activities Related to Waste Form Degradation: Progress Report” that describes the progress of R&D activities of ongoing modeling investigations specifically on nuclear waste glass degradation, Density Functional Theory (DFT) studies on clarkeite structure and stability, and electrochemical modeling of spent nuclear fuel (SNF). These activities are part of the newly-created Waste form Testing, Modeling, and Performance work package at Sandia National Laboratories (SNL). This work package is part of the “Inventory and Waste Form Characteristics and Performance” control account that includes various experimental and modeling activities on nuclear waste degradation conducted at Oak Ridge National Laboratory (ORNL), SNL, Argonne National Laboratory (ANL), and Pacific Northwest National Laboratory (PNNL).

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Modeling Activities Related to Waste Form Degradation: Progress Report

This report represents the milestone deliverable M4SF-22SN010309092 “Modeling Activities Related to Waste Form Degradation: Progress Report” that describes the progress of R&D activities of ongoing modeling investigations specifically on nuclear waste glass degradation, Density Functional Theory (DFT) studies on clarkeite structure and stability, and electrochemical model development of spent nuclear fuel (SNF). These activities are part of the Waste Form Testing, Modeling, and Performance work package at Sandia National Laboratories (SNL). This work package is part of the “Inventory and Waste Form Characteristics and Performance” control account that includes various experimental and modeling activities on nuclear waste degradation conducted at Oak Ridge National Laboratory (ORNL), SNL, Argonne National Laboratory (ANL), and Pacific Northwest National Laboratory (PNNL).

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

Continuum Model Development for Electrochemical Systems (Statement of Work: Ilenia Battiato Subcontract)

The subcontractor shall provide the services of qualified multiscale modelers to perform tasks that contribute to reaching the objectives of the LDRD project entitled Automated and Accelerated Continuum Model Development for Electrochemical Systems (Tracking number 24-ERD-051). These tasks relate to the development, deployment, and validation of multiscale models relevant to electrochemical systems.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A Diffuse-Interface Model for Predicting the Evolution of Metallic Negative Electrodes and Interfacial Voids in Solid-State Batteries with Homogeneous and Polycrystalline Solid Electrolyte Separators

Here, this paper presents a novel diffuse-interface electrochemical model that simultaneously simulates the evolution of the metallic negative electrode and interfacial voids during the stripping and plating processes in solid-state batteries. The utility and validity of this model are demonstrated for the first time on a cell with a sodium (Na) negative electrode and a Na-β″-alumina ceramic solid electrolyte (SE) separator. Three examples are simulated. First, stripping and plating with a perfect electrode/electrolyte interface; second, stripping and plating with a single interfacial void at the electrode/electrolyte interface; third, stripping with multiple interfacial voids. Both homogeneous and polycrystalline SEs with low and high-conductivity grain boundaries (GBs) are considered for all three examples. Heterogeneous GB conductivity marginally impacts the thickness of the Na electrode in cases with a perfect electrode/electrolyte interface. Moreover, it results in local changes to void growth due to the interactions between the void edge and the GBs. The void growth rate is a linear function of the flux of Na atoms at the void edge, which in turn depends on the applied current density. We also show that the void coalescence rate increases with applied current density and can be marginally influenced by GB conductivity.

batteries↗
DOE OSTI2021

Modeling the electrical double layer at solid-state electrochemical interfaces

Models of the electrical double layer (EDL) at electrode/liquid-electrolyte interfaces no longer hold for all-solid-state electrochemistry. Here we show a more general model for the EDL at a solid-state electrochemical interface based on the Poisson–Fermi–Dirac equation. By combining this model with density functional theory predictions, the interconnected electronic and ionic degrees of freedom in all-solid-state batteries, including the electronic band bending and defect concentration variation in the space-charge layer, are captured self-consistently. Along with a general mathematical solution, the EDL structure is presented in various materials that are thermodynamically stable in contact with a lithium metal anode: the solid electrolyte Li 7 La 3 Zr 2 O 12 (LLZO) and the solid interlayer materials LiF, Li 2 O and Li 2 CO 3 . Here, the model further allows design of the optimum interlayer thicknesses to minimize the electrostatic barrier for lithium ion transport at relevant solid-state battery interfaces.

25 ENERGY STORAGE↗
DOE OSTI2025

Microstructure Scale Lithium-Ion Battery Modeling: Part II. On In-Plane Heterogeneities and the Mechanisms that Regulate Them

Li-ion batteries performance and degradation are typically modeled at the macroscopic scale, that is neglecting in-plane heterogeneities that can arise from non-uniform electrode microstructures. Herein, a microstructure scale electrochemical model is used to quantify the impact of microstructure heterogeneity on cell performance during fast charging. The model predicts the electrolyte and solid concentration in-plane standard deviation can reach, respectively, ≈200 mol·m −3 and 6–7 kmol·m −3 locally. Further, the intercalation current density in-plane relative standard deviation can reach extremely high values, around 100% in the cathode and well above 100% in the anode graphite. These denote highly non-uniform lithiation rates and material utilization within each slice of the microstructure along the cell thickness. Non-uniform curvatures, at the particle scale (surface roughness) and between particles (size distribution), were found to initiate these in-plane heterogeneities, while an OCP-induced mechanism subsequently regulates them. The present model provides new insights into small length scale heterogeneity impact on battery performance not available with standard macro-scale/P2D modeling.

25 ENERGY STORAGE↗
DOE OSTI2024

Performance Prediction of High‐Entropy Perovskites La 0.8 Sr 0.2 Mn x Co y Fe z O 3 with Automated High‐Throughput Characterization of Combinatorial Libraries and Machine Learning

Perovskite oxides form a large family of materials with applications across various fields, owing to their structural and chemical flexibility. Efficient exploration of this extensive compositional space is now achievable through automated high-throughput experimentation combined with machine learning. In this study, we investigate the composition–structure–performance relationships of high-entropy La 0.8 Sr 0.2 Mn x Co y Fe z O 3±𝞭 perovskite oxides (0 < x, y, z <1; x+y+z≈1) for application as oxygen electrodes in Solid Oxide Cells. Following the deposition of a continuous compositional map using thin-film combinatorial pulsed laser deposition, compositional, structural, and performance properties are characterized using six different techniques with mapping capabilities. Random forests effectively model electrochemical performance, consistently identifying Fe-rich oxides as optimal compounds with the lowest area-specific resistance values for oxygen electrodes at 700 °C. Additionally, the models identify a statistical correlation between oxygen sublattice distortion—derived from spectral analysis of Raman-active modes—and enhanced performance.

high entropy oxides↗
DOE OSTI2021

A Segregated Approach for Modeling the Electrochemistry in the 3-D Microstructure of Li-Ion Batteries and Its Acceleration Using Block Preconditioners

Abstract Battery performance is strongly correlated with electrode microstructure. Electrode materials for lithium-ion batteries have complex microstructure geometries that require millions of degrees of freedom to solve the electrochemical system at the microstructure scale. A fast-iterative solver with an appropriate preconditioner is then required to simulate large representative volume in a reasonable time. In this work, a finite element electrochemical model is developed to resolve the concentration and potential within the electrode active materials and the electrolyte domains at the microstructure scale, with an emphasis on numerical stability and scaling performances. The block Gauss-Seidel (BGS) numerical method is implemented because the system of equations within the electrodes is coupled only through the nonlinear Butler–Volmer equation, which governs the electrochemical reaction at the interface between the domains. The best solution strategy found in this work consists of splitting the system into two blocks—one for the concentration and one for the potential field—and then performing block generalized minimal residual preconditioned with algebraic multigrid, using the FEniCS and the Portable, Extensible Toolkit for Scientific Computation libraries. Significant improvements in terms of time to solution (six times faster) and memory usage (halving) are achieved compared with the MUltifrontal Massively Parallel sparse direct Solver. Additionally, BGS experiences decent strong parallel scaling within the electrode domains. Last, the system of equations is modified to specifically address numerical instability induced by electrolyte depletion, which is particularly valuable for simulating fast-charge scenarios relevant for automotive application.

25 ENERGY STORAGE↗
DOE OSTI2020

Electrochemical Reduction Selectivity of Crotonaldehyde on Copper

Lignocellulosic waste is a potential feedstock for the generation of fuels and commodity chemicals, but existing conversion methods are too cost-intensive to be viable long-term solutions. Electrochemical reductions are promising for decentralized biomass valorization due to their modular scaling and capacity to run intermittently and without high temperatures or pressures. Using crotonaldehyde as a multi-functional model compound for the many partially unsaturated oxygenates found in processed biomass, we here demonstrate the production of butanal, butanol, butene, and butane (variously useful as commodity chemicals and major components of liquified petroleum gas) under ambient conditions by reductive bulk electrolysis with a copper mesh working electrode. We identify an optimum potential for reduced organic production under the reaction conditions and compare product distributions from reductions of intermediate species to further propose branching reaction pathways. Though butanal is typically the most abundant product from crotonaldehyde reduction, most of the butene and butane appear to result from a pathway involving initial reduction of the aldehyde group. We discuss evidence that selectivity is driven by interplay between crotonaldehyde reduction, local pH shifts due to the hydrogen evolution reaction, and changes in site reactivity and availability due to electrode fouling. This demonstration of model electrochemical biomass valorization also serves to inform further exploration into reduction of multi-functional molecules and electrochemical biomass processing in general.

biomass↗
DOE OSTI2026

From Cell to System: Accelerated hpc Simulations of BESS Aging under Frequency Regulation and Arbitrage use cases

Lithium-ion battery energy storage systems (BESS) packs have emerged as a leading solution for grid-scale energy storage, enhancing resiliency and balancing load fluctuations. Yet, experimental characterization of large-format LIB packs-particularly to assess performance and degradation over hundreds of cycles - demands substantial hardware investment and multi-year testing campaigns. In this work, we couple a hierarchical, physics-based modeling framework agnostic to electrode chemistries with high-performance computing to accelerate systems level evaluation by upto two orders of magnitude. Building on the open-source liionpack platform, we implement cell, module, and pack-scale electrochemical models enriched with mechanistic aging mechanisms and deploy them on an HPC cluster to simulate 150−200kWh systems over 500 - 1,000 cycles with in days. We subject these virtual B ESS to both constant-current cycling and realistic grid service profiles spanning frequency regulation, ramp-rate support, and energy arbitrage-and quantify the resulting degradation patterns. Our results reveal that localized cell aging can induce substantial nonuniformity at module and pack levels, with service-specific cycling protocols driving distinct aging modes. This rapid, multiscale modeling approach provides a powerful design-space exploration tool for optimizing electrical architecture, control strategies, and operational schedules to prolong pack lifetime and lower total cost of ownership.

Ayalasomayajula, Surya [ORNL] (ORCID:0009000860788↗
DOE OSTI2025

Microstructure Scale Lithium-Ion Battery Modeling: Part III. When and Where Lithium Plating Occurs and its Correlation with the Electrode Microstructure

Li-ion battery performance and degradation are closely related to the cell’s underlying electrode microstructure. Electrode microstructures are typically characterized with volume-averaged properties that neglect the impact of local heterogeneities. However, local heterogeneities create hot spots that can trigger degradation onset. Herein, a microstructure scale electrochemical model is used to investigate the impact of microstructure heterogeneity on lithium plating. The model predicts lithium plating is not uniform, even when considering a relatively small portion of the electrode (a cross-sectional area of 154×144 µm 2 ), preferring to plate on larger particles as compared to smaller particles. While local heterogeneities control where plating occurs, the model predicts that volume-averaged properties control when plating occurs. Additionally, the model predicts that the active material specific surface area has a linear relationship with the plating onset. However, the linear relationship between increased active material surface area and delayed plating response appears to be sensitive to the microstructure feature used to increase the active interface area. Here, a comparative case-study is explored where the specific surface area is increased by either reducing the active material particle diameter, adding open-porosity cracks, or increasing the active material surface roughness. The model predicts that increasing the specific surface area by reducing the active material particle diameter is the most effective strategy for delaying lithium plating. At 6C, reducing particle size is shown to be 3 and 20 times more effective than, respectively, adding open-porosity cracks and increasing surface roughness. A dual-layer electrode architecture combining gradations both for average properties and uniformities is eventually proposed to improve homogeneous material utilization and reduce degradation at high charge rates.

25 ENERGY STORAGE↗
DOE OSTI2025

Multiscale electro-chemo-mechanical model of high-capacity amorphous silicon anode material in sodium-ion batteries

The growing interest in sodium-ion batteries (SIBs) is fueled by their abundant resources and environmentally friendly nature, with amorphous silicon (a-Si) emerging as a promising anode material for enhancing capacity. However, the key challenge lies in sustaining reversible capacity during cycling. Here, in this work, we developed a multiscale electrochemical model incorporating an a-Si anode to elucidate the performance parameters of SIBs. Additionally, we integrated an electro-chemo-mechanical model at the particle level to capture stress generation, an essential factor in the degradation of high-capacity electrodes. Unlike existing models, our approach accounts for large-deformation chemo-mechanics at the particle scale and includes simulations under varying charge rates to explore multiscale behavior. The results reveal that coupled sodiation significantly prolongs complete cycling times and influences discharge dynamics, indicating that neglecting this coupling leads to an underestimation of actual capacity. Furthermore, we observed pronounced polarization effects at higher charge rates, resulting in heterogeneous stress distributions across the anode. With the identification of critical failure parameters for both active particles and binder materials, offering novel insights for mitigating degradation in high-capacity electrode systems.

Amorphous silicon anode↗
DOE OSTI2025

Static Subspace Approximation for Random Phase Approximation Correlation Energies: Applications to Materials for Catalysis and Electrochemistry

Modeling complex materials using high-fidelity, ab initio methods at low cost is a fundamental goal for quantum chemical software packages. The GW approximation and random phase approximation (RPA) provide a unified description of both electronic structure and total energies using the same physics in a many-body perturbative approach that can be more accurate than generalized-gradient density functional theory (DFT) methods. However, GW/RPA implementations have historically been limited to either specific materials classes or application toward small chemical systems. Here, the static subspace approximation allows for reduced cost full-frequency GW/RPA calculations and has previously been benchmarked thoroughly for GW calculations. Here, we describe our approach to including partial occupations of electronic orbitals in full-frequency GW and RPA calculations for the study of electrocatalysts. We benchmarked RPA total energy calculations using the subspace approximation across a diverse test suite of materials for a variety of computational parameters. The benchmarking quantifies the impact of different extrapolation procedures for representing the static polarizability at infinite screened cutoff, and shows that using screened cutoffs above 20-25 Ryd result in diminishing accuracy returns for predicting RPA total energies. Additionally, for moderately sized electrocatalytic models, 2-3 times fewer computational resources are used to compute RPA total energies by representing the static polarizability with 20-30% of the static subspace basis, with an error of approximately 0.01 eV or better in RPA adsorption energy calculations. Finally, we show that for these electrochemical models RPA can shift DFT adsorption energy shifts by up to 0.5 eV and that GW can frequently shift DFT eigenvalues of surface and adsorbate states by approximately 0.5-1 eV.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1981

Mechanism and models for zinc metal morphology in alkaline media

Based on experimental observations, a mechanism is presented to explain existence of the different morphologies of electrodeposited zinc in alkaline solution. The high current density dendrites appear to be due to more rapid growth on the nonbasal crystallographic planes than on the basal plane. The low current density moss apparently results from dissolution from the nonbasal planes at low cathodic voltages. Electrochemical models were sought which would produce such a phenomenon. The fundamental plating mechanism alone accounts only for different rates on different planes, not for zinc dissolution from a plane in the cathodic region. Fourteen models were explored; two models were in accord with the proposed mechanism. One involves rapid disproportionation of the zinc +1 species on the nonbasal planes. The other involves a redox reaction (corrosion) between the zinc-zincate and hydrogen-water systems.

May, C. E.↗