Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “battery modeling and simulation”

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

Control-Oriented Modeling of Lithium-Ion Batteries

Battery management systems (BMSs), which monitor and optimize performance while ensuring safety, require control-oriented models, i.e., models tailored to the design and implementation of estimation and control algorithms. Physics-based electrochemical models describe detailed battery phenomena, but are too computationally intensive for use in estimation and control applications. Single particle models (SPMs), which retain some of the physics of electrochemical models, are often used for control-oriented battery modeling since they are computationally efficient; however, they are only valid over very low frequency ranges and C-rates. Empirical equivalent circuit models (ECMs) are also used for control-oriented battery modeling since they are computationally efficient and can describe battery behavior over wide frequency ranges; however, they provide no physical understanding of the battery and, therefore, have limited applicability. Further, fractional order terms (e.g., Warburg impedances) are often employed, making the models unwieldy for use in the time domain. This work provides a control-oriented battery model that combines the benefits of SPM and ECM models, while overcoming their limitations. The proposed model incorporates some of the battery physics found in electrochemical models, can easily be used in both the time and frequency domains, and describes battery behavior over its entire frequency range. A linearized single particle model, which incorporates key electrochemical parameters, is used for modeling battery physics at very low frequencies. For low frequencies, integer-order linear systems are used to approximate diffusion physics described by Warburg impedances, and high frequency behavior is modeled by the double layer capacitance effect. The proposed battery model is more computationally efficient than full electrochemical models since it does not require the solution of PDEs, is accurate for a wider frequency range than the SPM considered in this paper, and does not suffer from the unwieldiness and limited applicability of empirical ECMs. Finally, the model is validated in the time and frequency domains via a comparison to pseudo-two-dimensional (P2D) model simulations and experimental data.

25 ENERGY STORAGE↗

Thermodynamic Modeling of Complex Solid Solutions in the Lu-H-N System via Graph Neural Network Accelerated Monte Carlo Simulations

Metal hydrides are important across diverse applications, such as hydrogen storage, batteries, gas sensors, nuclear reactions, and high-temperature superconductivity. Previous computational studies of metal hydrides under extreme pressures, e.g., 𝑂⁡(10 2 ) ⁢GPa, usually treat them as stoichiometric compounds without considering interstitial lattice disorder. As pressures become more moderate in the 𝑂⁡(10 0 ) ⁢GPa and below range, hydrogen disorder at interstitial lattice sites becomes prominent, e.g., in the N-doped Lu hydride that was recently claimed superconducting near 1 GPa. Further adding compositional complexity from alloying and/or multielement interstitial occupation makes elucidating pressure- and temperature-dependent observables intractable by first-principles calculations alone. We therefore propose a lattice graph neural-network surrogate modeling approach to predict configuration- and pressure-dependent equation-of-state properties. Their efficiency permits Monte Carlo simulations to calculate Gibbs energies and pressure-dependent phase diagrams, thereby revealing insights into the synthesis conditions required for achieving desired phase equilibria. We demonstrate this concept for the compositionally complex cubic Lu(H,N,Va) 3 system where three constituents (hydrogen, nitrogen and vacancy) have disordered multielement interstitial occupancies and insights into pressure-dependent phase equilibria are critically needed, e.g., N-doping levels can significantly lower dehydrogenation temperatures and provide a new strategy to optimize hydrogen-storage alloys. This work can improve the thermodynamic understanding of the Lu-H-N system and help rational synthesis of N-doped Lu hydrides, but more generally demonstrates an efficient approach to model pressure-dependent thermodynamics of multicomponent solid solutions.

Monte Carlo methods↗

Impact Modeling and Testing of Pouch and Prismatic Cells

Understanding battery response under impact is critical to improve the safety of electrified vehicles. Nevertheless, predicting the impact behaviors of batteries is not straightforward since a battery cell usually contains hundreds of thin layers with dramatically different material properties and multiple physical processes occur simultaneously during cell deformation. Here we utilized both empirical tests and numerical models to capture the failure process of pouch and prismatic cells in various impact scenarios. In each test, a cell was hit once by an indenter dropped from a certain height. During which the cell penetration, loading force, voltage and temperature were monitored to characterize the cell’s response. Meanwhile, numerical models were developed to capture the coupled mechanical, electrical, electrochemical and thermal responses of batteries. In these models, the cell bulk was treated as a homogeneous part to achieve computational efficiency required by large-scale simulations, and it was represented by the geologic cap model that allows both shear and compaction deformation. Simulation results showed agreement with experimental data in essential features of cell behaviors during impact. Details of the test setup, model development and cell failure behaviors are presented in this paper. Additionally, capabilities, limitations and future improvement of the battery safety modeling are discussed.

25 ENERGY STORAGE↗

Dynamics of Intra-Cell Thermal Front Propagation in Lithium-Ion Battery Safety Issues

Thermal runaway (TR), a critical failure mode in lithium-ion batteries (LIBs), poses significant safety risks and hinders wider application of LIBs. TR typically begins at a localized heat source and spreads across the cell. Understanding thermal front propagation (TFP) characteristics, such as front and velocity, is crucial for assessing energy release and temperature distribution for battery hazardous estimation. Recent studies assume that TR within cells propagates at a near-constant velocity, based on the reaction kinetics and thermal properties. Here, in this study, an intra-battery TR model is further proposed and it indicates that TFP velocity stabilizes when the front is distanced from the heat source. Theoretical estimates for propagation velocity and front are developed and validated through numerical simulations and experimental tests from the NREL Battery Failure Databank. The energy release rate during TFP and the impact of preheating based on a point heat source are explored. This work clarifies the long-standing clouds of the thermal font propagation behaviors within the single cell, highlights the power and beauty of mathematics modeling to describe the complicated thermal behaviors, and provides important guidelines for thermal hazardous understanding for next-generation batteries.

25 ENERGY STORAGE↗

LIM1TR: Lithium-Ion Modeling with 1-D Thermal Runaway (V.1.0)

LIM1TR (Lithium-Ion Modeling with 1-D Thermal Runaway) is an open-source code that uses the finite volume method to simulate heat transfer and chemical kinetics on a quasi 1-D domain. The target application of this software is to simulate thermal runaway in systems of lithium-ion batteries. The source code for LIM1TR can be found at https://github.com/ajkur/lim1tr. This user guide details the steps required to create and run simulations with LIM1TR starting with setting up the Python environment, generating an input file, and running a simulation. Additional details are provided on the output of LIM1TR as well as extending the code with custom reaction models. This user guide concludes with simple example analyses of common battery thermal runaway scenarios. The corresponding input files and processing scripts can be found in the “Examples” folder in the on-line repository, with select input files included in the appendix of this document.

25 ENERGY STORAGE↗

Mathematical Model for Li-S Cell with Shuttling-Induced Capacity Loss Approximation

Lithium sulfur (Li-S) batteries have the potential to outperform the current lithium ion batteries and transform the technology of the future. However, dissolution, diffusion, and shuttling of the dissolved polysulfides result in parasitic reactions and substantial capacity loss. To provide a better understanding of the shuttling process, a 1D porous electrode mathematical model has been developed in this paper. An approximation method is used to account for the shuttling-induced capacity loss by adding an extra source/sink term in the material balance equations for the species involved in the parasitic reactions. Shuttling constants used in the source terms can be determined by fitting the model predictions to the experimental measurements. Here, the results showed that by including the approximation method, the model was able to predict the active material loss and the continuous decrease of volume fractions of Li 2 S on the cathode surface. The model sheds light on the capacity loss mechanism occurring inside the cell as a result of the shuttling of polysulfides.

25 ENERGY STORAGE↗

Contributions of Solar + Storage to Future Power Needs

The costs of solar photovoltaics (PV) have been dropping in recent years, leading to increasing installations of solar PV systems and growing interest in how the technology will impact the electric grid. Additionally, as battery costs decline, there is increasing interest in understanding the benefits and limitations of battery storage for the grid as well as potential co-benefits of coupling battery storage and photovoltaics. This paper presents recent methodological developments to more accurately represent the value and limitations of coupled PV and battery systems in capacity expansion models. These are demonstrated using the Resource Planning Model (RPM), which co-optimizes capacity investments, transmission investments, and reduced-order dispatch through 2045. We simulate the evolution of the generation and transmission system in the Western Interconnection both with and without a coupled PV and battery technology option under two core scenarios - a baseline scenario and a high renewable penetration scenario - coupled with sensitivities assuming low and mid PV and battery cost projections. We analyze the degree to which a variety of configurations of coupled PV + Battery systems are able to meet future electricity needs, under which system conditions this coupled technology becomes beneficial to the electric grid and identify the displaced technologies when coupled solar and storage is introduced into the model. We find that capturing the benefits of co-locating these systems becomes increasingly important as the penetration of solar PV rises, and that in particular it is important to accurately capture the ability of these technologies to shift energy, provide firm capacity, and reduce expected curtailment.

14 SOLAR ENERGY↗

MPC4CLR (Model-Predictive-Control-for-Critical-Load-Restoration-in-Power-Distribution-Systems) [SWR-22-24]

Model predictive control (MPC) is a system or process control technique for making decisions under uncertainty via rolling look-ahead optimizations at each control step where only the current step decisions are applied, and the rest are discarded. In this work, we developed an MPC for a critical load restoration (CLR) in power distribution systems to recover system service (electricity delivery) following an extreme event-triggered substation outage. The method considers the problem of controlling distributed energy resources (DERs) of the distribution system with the objective of achieving maximum load pick up while satisfying distribution network flow and voltage constraints. A linearized optimal power flow (OPF) model is employed to represent the physics of the network. The problem formulation is augmented with a ramping (up) reserve product for the DERs to ensure improved and upward monotonic load restoration as time evolves. Simulation analysis and performance tests are performed using a modified IEEE 13-bus test feeder integrated with wind, solar, microturbine, and energy storage battery. The software is developed using various software packages in Julia and Python. The MPC model is implemented using the JuMP optimization language in Julia while the data analytics including renewable generation and load demand forecasts, running the MPC simulation and visualizations is performed in Python.

Eseye, Abinet Tesfaye↗

Modeling Lithium Diffusion in Battery Cathodes Considering Chemo-Mechanically Induced Damage

This talk will present a 3D, continuum-level damage model for simulating Lithium diffusion within generated Li_{x}Ni_{0.5}Mn_{0.3}Co_{0.2} (NMC 532) secondary cathode particles. The primary motivation of the particle-level model is to inform cathode-particle design and determine charging profiles that reduce cathode fracture. The model considers NMC 532 secondary particles containing an agglomeration of anisotropic, randomly oriented grains. The model predicts that secondary-particle fracture is primarily due to non-ideal grain interactions with slight dependence on high-rate charge demands. The model predicts that small secondary-particles with large grains develop significantly less damage than larger secondary particles with small grains. Finally, the model predicts most of the chemo-mechanical damage accumulates in the first high-rate cycles. This chemo-mechanical “damage saturation” effect indicates that initial secondary-particle fracture occurs within the first few cycles, while long-term cathode degradation is not solely chemo-mechanically induced.

47 OTHER INSTRUMENTATION↗

Modeling Methods for Capturing System Interactions of Combined Technologies: A Study of PV + Battery

The costs of solar photovoltaics (PV) have been dropping in recent years, leading to increasing installations of solar PV systems and growing interest in how the technology will impact the electric grid at higher penetrations. Additionally, as battery costs decline it becomes important to understand the benefits and limitations of battery storage for the grid as well as potential benefits of co-locating battery storage and PV, particularly within the realm of future system planning. However, it is non-trivial to capture these potential benefits within a linearized capacity expansion model (CEM). This paper presents methodological developments to more fully represent the value and limitations of coupled PV and battery systems (PV + Battery) in CEMs using the Resource Planning Model (RPM), which co-optimizes capacity investments, transmission investments, and reduced-order dispatch in the Western Interconnection of North America through 2045. We use the model to simulate the evolution of the generation and transmission system under two core scenarios - a baseline scenario and a high renewable penetration scenario - coupled with sensitivities assuming low and midline PV and battery cost projections. When incorporating PV + battery, we find that it is important for CEMs to capture the ability of the coupled technology to provide firm capacity and reduce expected curtailment, compared to independent systems. These interactions can have even more dramatic impacts at higher solar penetrations.

14 SOLAR ENERGY↗

Grid Resiliency with a 100% Renewable Microgrid

San Diego Gas & Electric Company (SDG&E) installed America’s first and largest utility-scale microgrid in Borrego Springs in 2013. The first generation Borrego Springs Microgrid utilized diesel generators to form and stabilize the microgrid island, with support from grid-scale batteries and local solar photovoltaic (PV) generation. In this project, SDG&E in partnership with National Renewable Energy Laboratory (NREL) demonstrated through modeling, simulation and utility field testing that blackstart and islanding of the microgrid can be led with 100% renewable, inverter based resources (IBRs), to help reduce community reliance on conventional generation resources. Through equipment upgrades, grid-forming island leader capability was transitioned to a battery IBR instead of the Borrego Springs Microgrid diesel generators. A new microgrid controller was integrated to the microgrid and programmed to control and manage multiple energy storage systems. Synchrophasor and other power quality data verified autonomous, high-speed response of the IBRs through blackstart, islanding, and load step testing. Results of project field evaluations provide distribution systems operators (DSO) with increased confidence that renewable, IBR can replace traditional generators to blackstart and island microgrids and rapidly establish stable island frequency with rapid changes in peak power demand. Importantly, the project validated the integration feasibility of a distributed energy resource management system (DERMS) controller that manages multiple grid-forming and grid-following IBRs, establishing a standard design interface to reduce the complexity of integrating new DERs in the future and supporting replication by the industry. As a result of learnings in this project, SDG&E has implemented the microgrid controller strategy at multiple other microgrid sites, thereby validating the replicability of the solution. Hardware-in-the-loop (HIL) simulations including power and controller HIL hardware — along with electromagnetic transient (EMT) simulations of Borrego Springs Microgrid —informed adjustments to inverter parameters and were important to characterize the performance of the IBRs in relevant operating conditions before deployment. The EMT and HIL simulations of islanding the entire community are important contributions in providing confidence in IBR performance prior to future islanding of the community in the field. High-fidelity EMT and/or HIL simulation of IBRs can de-risk field operations, and its relevance and importance as a tool is increasing as distribution grids and microgrids become more complex and dynamic with an increasing proportion of renewable generation, distributed energy storage, and two-way power and energy flows.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Advanced Distributed Wind Turbine Controls Series: Part 2-Wind Energy in Isolated Grids – Microgrids, Infrastructure Resilience, and Advanced Controls Launchpad (MIRACL)

In an isolated grid, wind turbines are typically deployed to provide energy to maximize energy production, reduce diesel-fuel consumption, reduce carbon emissions, and reduce costs for energy and fuel transportation. However, in addition to solely providing energy to the power system, wind turbines contain rotating masses and inverter-based controls that can enable various reliability and resilience services through advance controls. As part of the Microgrids, Infrastructure Resilience, and Advanced Controls Launchpad (MIRACL) this paper demonstrates, through desktop simulations, advanced wind turbine controls that can be employed to support higher contributions of wind in isolated grids, and to demonstrate ways that wind can play a role in supporting stability of an isolated grid. This isolated grid used in these desktop simulations is comprised of a wind turbine (600 kW), solar PV (430 kW), battery energy storage system (1 MW/1MWh), a simulated diesel generator (2 MW) and various types of loads (critical, dynamic). We developed a model of the subsystems in MATLAB/Simulink and validated them with available data from real-world components on NREL's Flatirons Campus. These validated models are then configured for various case studies. We compare the output of the desktop simulation with a baseline case with the diesel generator. Active and reactive power control of the wind turbine can help improve frequency and voltage responses in the isolated grid, respectively. By utilizing a small integrated battery energy storage system in the DC-link of the wind turbine, we also demonstrate that wind turbines can help blackstart a critical load comparable to its rated power and support other renewables (e.g. solar PV) come online and pick up an additional load. This report illustrates some of these reliability and resilience services a wind turbine can provide in an isolated grid.

17 WIND ENERGY↗

Scenario-based analysis of electric vehicle adoption in the United States: Technology, infrastructure, and electricity pricing

This work investigates the impact of battery technology advancement, charging infrastructure development, and time-of-use (TOU) electricity pricing on vehicle adoption by 6 powertrain types in the United States through 2050. Using the Market Acceptance of Advanced Automotive Technologies (MA3T) model, we simulate 15 scenarios, examining individual cost factors and their combinations. We assess outcomes through market share, consumer surplus, and energy consumption. Results show that battery cost reductions are the strongest driver of EV adoption, increasing 2050 battery electric vehicle (BEV) share by 27 percentage points over baseline, raising annual consumer surplus by $511 per household, and reducing cumulative energy consumption by 16,610 trillion Btu. These gains are two to five times larger than those from other individual factors. Reducing home charging installation costs produces moderate impact, while TOU pricing alone yields only small gains, raising 2050 BEV market share by 1–2 percentage points. However, when cost factor improvements are combined, their effects are amplified beyond simple additivity. Pairing modest battery cost reductions with charging installation cost reductions and TOU pricing results in the largest 2050 BEV sales combined impact. The analysis demonstrates that moderate progress targeting multiple cost barriers may be more impactful than focusing on any single barrier.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Editors’ Choice—Perspective—Challenges in Moving to Multiscale Battery Models: Where Electrochemistry Meets and Demands More from Math

There has been significant recent interest in studying multiscale characteristics of current and next-generation batteries, including lithium-metal and lithium-sulfur batteries. Advances in computing power make researchers believe that the detailed multiscale models can be efficiently simulated to arrive at the insights for the degradation and performance loss; however, this is not true and special attention needs to be paid to local singularities, boundary layers, moving boundaries, etc. This article presents 2D examples that illustrate the importance of grid convergence studies, provides well-defined detailed models to test the efficiency of numerical schemes, and discusses the associated simulation challenges.

Shah, Krishna (ORCID:0000000278026361)↗

SIRENOpt.jl

SAND2026-22945O SIRENOpt.jl is a Julia software package for prototype hybrid power, storage, and platform dynamics. It integrates solar, wind, wave, hydrokinetic, diesel, generator, converter, battery, hydrogen, desalination, mooring, and floating-platform model interfaces in an automatic-differentiation-friendly simulation framework. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Michelen Strofer, Carlos [Sandia National Lab. (SN↗

Discrete-Element and Material-Point Method (DEM and MPM) Based Solvers for Sustainable Technologies

We present the use of discrete element method (DEM) and material point method (MPM) in three relevant green technology applications that include biomass feedstock handling, lithium-ion battery manufacturing, and high-pressure reverse osmosis. Our open-source DEM and MPM solvers are developed using performance portable grid and particle management library, AMReX, thus enabling superior performance on NVIDIA and AMD GPUs with > 100 million particles. Our DEM solver resolves the motion of individual particles in a granular system and includes a bonded sphere method for modeling non-spherical particles along with Hertzian and liquid bridge-based contact models. We simulate highly variable biomass feedstock flows in large-scale hoppers for biofuel production and electrode calendering in battery manufacturing using DEM. Our simulations predict flow blockage in large scale biomass hoppers and electrode microstructure variations, thus providing valuable information for biofuel and battery manufacturers, respectively. The second half of the talk will be on MPM and its application towards pore resolved simulations of reverse osmosis membranes under compressive loads. We present a validation study of our MPM simulations with membrane microscopy imaging thus providing useful insights on membrane stability under high pressure conditions. We also present a spectral stability analysis of using linear hat, quadratic and cubic spline basis in MPM indicating regions of numerical stability.

BIOMASS FUELS,MATHEMATICS AND COMPUTING↗

Modeling assisted synthesis of Zr-doped Li 3-x In 1-x Zr x Cl 6 with ultrahigh ionic conductivity for lithium-ion batteries

All-solid-state lithium-ion batteries (ASSLBs) are an important milestone for the future of energy storage because of their capability of impressive energy density and outstanding safety. However, oxide and sulfide solid-state electrolytes (SSEs) suffer from either low ionic conductivity or poor chemical stability. In contrast, halide-based SSEs, are promising as candidate materials owing to high conductivity, good stability, and broad cathode compatibility. Though element doping of the SSEs is an effective and common approach to further improve their electrochemical properties, dopant exploration and optimization through solely experimental trials are both costly and time-consuming. For this aspect, computational simulations for dopant element and concentration screening are adopted in this research and zirconium is selected as a suitable dopant for Li 3 InCl 6 . Further, the synthesized Li 2.75 In 0.75 Zr 0.25 Cl 6 exhibited Li ionic conductivity of 5.82 x 10 -3 Scm -1 at room temperature, which is the highest among reported halide SSEs. The ASSLB formed with Li 2 CoO 2 -Li 2.75 In 0.75 Zr 0.25 Cl 6 -Li/In delivers a high initial capacity of 129.3 mAh∙g -1 . Conclusively, this work provides an effective approach which combines computational modeling and experimental verification for the development of halide SSEs with improved stability and conductivity. The successful design approach and compelling results provide further possibilities and capabilities in future SSE research.

25 ENERGY STORAGE↗

Simulating hindered grain boundary diffusion using the smoothed boundary method

Abstract Grain boundaries can greatly affect the transport properties of polycrystalline materials, particularly when the grain size approaches the nanoscale. While grain boundaries often enhance diffusion by providing a fast pathway for chemical transport, some material systems, such as those of solid oxide fuel cells and battery cathode particles, exhibit the opposite behavior, where grain boundaries act to hinder diffusion. To facilitate the study of systems with hindered grain boundary diffusion, we propose a model that utilizes the smoothed boundary method to simulate the dynamic concentration evolution in polycrystalline systems. The model employs domain parameters with diffuse interfaces to describe the grains, thereby enabling solutions with explicit consideration of their complex geometries. The intrinsic error arising from the diffuse interface approach employed in our proposed model is explored by comparing the results against a sharp interface model for a variety of parameter sets. Finally, two case studies are considered to demonstrate potential applications of the model. First, a nanocrystalline yttria-stabilized zirconia solid oxide fuel cell system is investigated, and the effective diffusivities are extracted from the simulation results and are compared to the values obtained through mean-field approximations. Second, the concentration evolution during lithiation of a polycrystalline battery cathode particle is simulated to demonstrate the method’s capability.

Materials Science↗