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

Microstructure-Based Modeling of Inner Oxygen Pressure in Solid Oxide Electrolysis Cells

One major degradation mechanism during long-term operation of solid oxide electrolysis cells (SOECs) is delamination of oxygen electrodes (OEs). The driving force for the electrode delamination could be the generated high inner oxygen pressure near the electrode-electrolyte interface during operation. However, the effects of transport properties and electrode thickness on the inner oxygen partial pressure are not well understood. Here, a microstructure-based electrochemical model which includes the conduction of electrons and oxygen ions coupled with Butler-Volmer-type chemical reactions at triple-phase-boundaries (TPBs), is employed to investigate the oxygen pressure in lanthanum strontium manganate (LSM)-based SOECs. The model is applied to both two-dimensional (2D) prototype microstructures and three-dimensional (3D) realistic microstructures, and the oxygen pressure is analyzed as a function of transport properties and electrode thickness under both potentiostatic and galvanostatic operations. The simulation results suggest strategies to suppress electrode delamination. The simulation results are compared to an analytical solution, and the discrepancies are attributed to the Butler-Volmer-type kinetics included in the microstructure-based model.

Xue, Fei↗
DOE OSTI2025

Microstructure-based modeling of inner oxygen pressure in solid oxide electrolysis cells: Analysis of electrode delamination and mitigation

One major degradation mechanism during long-term operation of solid oxide electrolysis cells (SOECs) is delamination of oxygen electrodes (OEs). The driving force for the electrode delamination could be the generated high inner oxygen pressure near the electrode-electrolyte interface during operation. However, the effects of transport properties and electrode thickness on the inner oxygen partial pressure are not well understood. Here a microstructure-based electrochemical model, which includes the conduction of electrons and oxygen ions coupled with Butler-Volmer-type chemical reactions at triple-phase-boundaries (TPBs), is employed to investigate the oxygen pressure in lanthanum strontium manganate (LSM)-based SOECs. The model is applied to both two-dimensional (2D) prototype microstructures and three-dimensional (3D) realistic microstructures, and the oxygen pressure is analyzed as a function of transport properties and electrode thickness under both potentiostatic and galvanostatic operations. The simulation results suggest strategies to suppress electrode delamination. The simulation results are compared to an analytical solution, and the discrepancies are attributed to the Butler-Volmer-type kinetics included in the microstructure-based model.

25 ENERGY STORAGE↗
DOE OSTI2020

Enabling fast charging of lithium-ion batteries through secondary-/dual- pore network: Part II - numerical model

To increase the market share of electric vehicles, it is desirable to reduce the battery charge times, which are significantly limited by poor electrolyte transport. A high rate charging is achievable by using expensive and low energy density cells with thin electrodes. For higher energy density cells, new electrolytes with improved conductivity and diffusivity and/or electrodes with advanced architecture are required to boost the electrolyte transport, leading to a more uniform utilization of active materials. In our previous work, an analytical model was developed to investigate the effect of secondary pore network (SPN) on electrolyte transport and the configuration of SPN was optimized by enforcing equal characteristic diffusion times in through-plane and in-plane directions. Here, to evaluate the effect of SPN on the fast-charging capability of lithium-ion batteries, a 2D physics-based electrochemical model is developed with SPN in either one or both electrodes. Additionally, the effect of SPN on cell energy density and lithium plating is investigated for cells with different loadings and electrode porosities. Combining SPN with elevated charging temperatures, the model predicts that the volumetric discharge energy density of a 3 mA h/cm2 cell can reach 270 Wh/L after a 6C constant-current charging.

25 ENERGY STORAGE↗
DOE OSTI2022

Probing the Role of Multi-scale Heterogeneity in Graphite Electrodes for Extreme Fast Charging

Electrode-scale heterogeneity can combine with complex electrochemical interactions to impede lithium-ion battery performance, particularly during fast charging. This research investigates the influence of electrode heterogeneity at different scales on the lithium-ion battery electrochemical performance under operational extremes. We employ image-based mesoscale simulation in conjunction with a three-dimensional electrochemical model to predict performance variability in 14 graphite electrode X-ray computed tomography data sets. Our analysis reveals that the tortuous anisotropy stemming from the variable particle morphology has a dominating influence on the overall cell performance. Cells with platelet morphology achieve lower capacity, higher heat generation rates, and severe plating under extreme fast charge conditions. On the contrary, the heterogeneity due to the active material clustering alone has minimal impact. Our work suggests that manufacturing electrodes with more homogeneous and isotropic particle morphology will improve electrochemical performance and improve safety, enabling electromobility.

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

Sodium Ion Expansion Power Block for Distributed CSP

The Sodium Ion Expansion Power Block for Distributed CSP was a three-plus-one-year effort under the Concentrating Solar Power: Advanced Projects Offering Low LCOE Opportunities (CSP: APOLLO) funding program within the U.S. Department of Energy Solar Energy Technologies Office. The primary objective of this project is to develop a dual-stage modular sodium thermal electrochemical converter (Na-TEC) heat engine power block, which can be potentially integrated with either a small-scale dish solar or large-scale heliostats and parabolic trough CSP. Na-TEC is a heat engine that generates electricity through the isothermal expansion of sodium ions. The Na-TEC is a closed system that can theoretically achieve conversion efficiencies above 45% when operating between thermal reservoirs at 1150 K and 550 K. However, thermal designs have confined previous single-stage devices to thermal efficiencies below 20%. To mitigate some of these limitations, we consider dividing the isothermal expansion into two stages; one at the evaporator temperature (1150 K) and another at an intermediate temperature (650 K –1050 K). This dual-stage Na-TEC takes advantage of regeneration and reheating, and could be amenable to better thermal management. In light of this, we first designed and developed a thermo-electrochemical model, and thermodynamically demonstrated how the dual-stage device can improve the efficiency by up to 8% points over the best performing single-stage device. We also established an application regime map for the single- and dual-stage Na-TEC in terms of the power density and the total thermal parasitic loss. Moreover, a thermal design of an axisymmetric dual-stage Na-TEC is developed to guide the scale-up and fabrication of sub-components of prototype module. A reduced-order finite-element model is used in conjunction with a Na-TEC thermodynamic model that was developed to determine the total parasitic heat loss of this dual-stage design. A number of simplifications are applied in the reduced-order model to decrease the computational time while maintaining acceptable accuracy. According to this analysis, a maximum efficiency of 29% and a maximum power output of 125 W can be achieved. Ultimately, we were able to demonstrate thermal efficiency improvements of the Na-TEC heat engine from 19% up to 40.3%, in a dual-stage (non-optimized) prototype module that we designed, fabricated, and tested with high temperature stage at 923 K. Furthermore, a cost-performance analysis for this improved dual-stage design was carried out for distributed-CSP systems. A high-level techno-economic analysis (TEA) explores four scenarios where a Na-TEC is used as the heat engine for a distributed-CSP system. Overnight capital cost and levelized cost of electricity (LCOE) are estimated for a system lifetime of 30 years, revealing that overnight capital costs in a range from $3.57 to $17.71 per We are feasible, which equate to LCOEs from 6.9 to 17.2 cents/kWh e -1 . This analysis makes a significant contribution by concurrently quantifying the efficiency and unit costs for a range of multistage configurations, and demonstrating that a Na-TEC may be a promising alternative to Stirling engines for distributed-CSP systems at residential scale of 1–5 kW e .

14 SOLAR ENERGY↗
DOE OSTI2020

Fast Proton Insertion in Layered H 2 W 2 O 7 via Selective Etching of an Aurivillius Phase

H 2 W 2 O 7 , a metastable material synthesized via selective etching of the Aurivillius-related Bi 2 W 2 O 9 , is shown as an electrode for high power proton-based energy storage. Comprehensive structural characterization is performed to obtain a high-fidelity crystal structure of H 2 W 2 O 7 using an iterative approach that combines X-ray diffraction, neutron pair distribution function, scanning transmission electron microscopy, Raman spectroscopy, and density functional theory modeling. Electrochemical characterization shows a capacity retention of ≈80% at 1000 mV s –1 (1.5-s charge/discharge time) as compared to 1 mV s –1 (≈16-min charge/discharge time) with cyclability for over 100 000 cycles. Energetics from density functional theory calculations indicate that proton storage occurs at the terminal oxygen sites within the hydrated interlayer. Last, optical micrographs collected during in situ Raman spectroscopy show reversible, multicolor electrochromism, with color changes from pale yellow to blue, purple, and last, orange as a function of proton content. These results highlight the use of selective etching of layered perovskites for the synthesis of metastable transition metal oxide materials and the use of H 2 W 2 O 7 as an anode material for proton-based energy storage or electrochromic applications.

25 ENERGY STORAGE↗
DOE OSTI2023

Composite Lithium Metal Structure to Mitigate Pulverization and Enable Long‐Life Batteries

In lithium metal batteries, non‐uniform stripping of lithium results in pit formation, which promotes subsequent non‐uniform, dendritic deposition. This viscous cycle leads to pulverization of lithium which promotes cell shorting or capacity degradation, symptoms further exaggerated by high electrode areal loading and lean electrolytes. Here, to address this challenge, a composite lithium metal anode is engineered that contains uniformly distributed, nanometer‐sized carbon particles. This composite lithium is shown to strip more uniformly since the growth of non‐uniform pits is intercepted by the carbon particles. This mechanism is corroborated by a continuum electrochemical model. Subsequent lithium deposition on carbon particles is also found to be more uniform than on the surface with irregular pits. Notably, the pulverization rate of composite lithium is 26 times slower than that of commercial lithium. Moreover, in a Li‐S battery with sulfurized polyacrylonitrile cathode, the use of the composite anode extends the cycle life by three times when the areal capacity is 8 mAh cm −2 . The approach of using an engineered lithium composite structure to address challenges during both stripping and plating can inform future designs of lithium metal anodes for high areal capacity operations.

high areal capacity↗
DOE OSTI2020

Effect of Anode Porosity and Temperature on the Performance and Lithium Plating During Fast-Charging of Lithium-Ion Cells

Twenty-four single-layer approximate to 32 mAh pouch cells are tested to determine the effect of electrode porosity on lithium plating. Twelve cells contain a graphite electrode that is 26% porous, and 47% for the other twelve. The cells are cycled using a 6-C charge and a C/2 discharge protocol at temperatures in the range of 20-50 degrees C. A macro-homogeneous electrochemical model and microstructure analysis tool set are used to help interpret experimental observations for the effect of anode porosity and ambient temperature on fast-charging performance. Additionally, comparison between the two also highlights gaps in current theoretical understanding that need to be addressed. In post-test examination, lithium plating is seen in all cells, regardless of porosity. Elevated temperature is shown to reduce the amount of lithium plating and improve initial fast-charge capacity, but also changes the rate of other, less well-understood degradation mechanisms. Apparent kinetic rate laws, At + Bt 1/2 , where A and B are constants, can be fit to most of the capacity loss and resistance increase data. The relative magnitudes of A and B change with temperature and porosity. The capacity loss data at 50 degrees C from the high-porosity cells are fit by a logistics rate law.

25 ENERGY STORAGE↗
DOE OSTI2025

Optimal sizing of battery energy storage systems for peak shaving and demand response using a degradation-aware Bayesian Optimization-Mixed-Integer Linear Programming framework

The increasing integration of renewable energy and rising electricity demand highlight the importance of battery energy storage systems for peak shaving and demand response. Unlike prior approaches that overlook operational impacts on degradation, this study proposes a Bayesian Optimization–Mixed Integer Linear Programming framework for optimal battery energy storage system sizing. In this framework, Mixed Integer Linear Programming determines short-term scheduling while a calibrated electrochemical model iteratively evaluates degradation. The central hypothesis is that the framework can efficiently identify optimal sizes that yield realistic and economically robust outcomes. The method is tested across three scenarios: peak shaving, peak shaving with energy-reduction demand response, and peak shaving with power-reduction demand response. Results show that the framework converge to the optimum within 20 iterations out of 150 possible sizes. Under baseline conditions, the framework consistently selects the smallest feasible system, minimizing unnecessary degradation costs from oversized storage. Sensitivity analyses reveal that larger systems are favored as demand rates or incentives increase. Comparisons of demand response programs indicate that power-reduction demand response offers greater economic benefits than energy-reduction demand response, although demand savings from peak shaving remain the dominant contributor to overall performance. This study demonstrates that the proposed framework balances computational tractability with degradation fidelity, identifies critical economic thresholds for investment, and offers a practical, flexible tool to guide industrial stakeholders in cost-effective battery energy storage system deployment.

Batteries↗