Simulation of the Deformation for Cycling Chemo-Mechanically Coupled Battery Active Particles with Mechanical Constraints

R. Schoof,G. F. Castelli,W. Dörfler

DOI: https://doi.org/10.1016/j.camwa.2023.08.027

2023-09-06

Abstract:Next-generation lithium-ion batteries with silicon anodes have positive characteristics due to higher energy densities compared to state-of-the-art graphite anodes. However, the large volume expansion of silicon anodes can cause high mechanical stresses, especially if the battery active particle cannot expand freely. In this article, a thermodynamically consistent continuum model for coupling chemical and mechanical effects of electrode particles is extended by a change in the boundary condition for the displacement via a variational inequality. This switch represents a limited enlargement of the particle swelling or shrinking due to lithium intercalation or deintercalation in the host material, respectively. For inequality constraints as boundary condition a smaller time step size is need as well as a locally finer mesh. The combination of a primal-dual active set algorithm, interpreted as semismooth Newton method, and a spatial and temporal adaptive algorithm allows the efficient numerical investigation based on a finite element method. Using the example of silicon, the chemical and mechanical behavior of one- and two-dimensional representative geometries for a charge-discharge cycle is investigated. Furthermore, the efficiency of the adaptive algorithm is demonstrated. It turns out that the size of the gap has a significant influence on the maximal stress value and the slope of the increase. Especially in two dimension, the obstacle can cause an additional region with a lithium-poor phase.

Numerical Analysis

What problem does this paper attempt to address?

The problem that this paper attempts to solve is: the behavior simulation of silicon anode materials in lithium - ion batteries during the charging and discharging process when the volume expansion and contraction caused by lithium intercalation and de - intercalation are mechanically constrained. Specifically, the article focuses on how to improve the existing thermodynamically consistent continuum models by introducing inequality boundary conditions (such as obstacle problems) to more accurately describe the deformation behavior of electrode particles under the coupling of chemical and mechanical effects.

### Problem Background

1. **High Energy Density Requirement**: Next - generation lithium - ion batteries use silicon anodes, which have a higher energy density compared to traditional graphite anodes.
2. **Volume Expansion Problem**: However, silicon anodes can experience up to 300% volume expansion during the charging and discharging process, which may lead to severe mechanical stress, ultimately causing particle fracture, shortening battery life and accelerating the aging process.
3. **Limitations of Existing Models**: In previous simulation studies, it is usually assumed that electrode particles can expand freely, ignoring the restrictions of the external environment on particle expansion (such as battery casings, current collectors or the contact of other electrode particles).

### Research Objectives

To better understand these problems, this paper proposes an improved model that can:
- Simulate the expansion and contraction behavior of electrode particles in a limited space;
- Take into account the influence of mechanical constraints (such as obstacles) on particle deformation;
- Efficiently solve complex problems with inequality boundary conditions by numerical methods.

### Main Contributions

1. **Introduction of Obstacle Boundary Conditions**: Obstacle boundary conditions are introduced through variational inequalities to simulate the situation where particles are in contact with obstacles during expansion.
2. **Efficient Numerical Algorithm**: Combining the primal - dual active - set algorithm (interpreted as the semi - smooth Newton method) and the space - time adaptive algorithm, the effective application of the finite element method is realized.
3. **Simulation of Multi - Dimensional Geometric Structures**: Through the simulation of one - dimensional and two - dimensional representative geometric structures, the chemical and mechanical behaviors of particles under different conditions are demonstrated.

### Formula Representation

To ensure the correctness and readability of the formulas, the following are some key formulas involved in the paper:

1. **Total Deformation Gradient Tensor**:
\[
F = F_{\text{ch}} F_{\text{el}}
\]
where \(F_{\text{ch}}\) represents the chemical part and \(F_{\text{el}}\) represents the elastic part.

2. **Free - Energy Density**:
\[
\psi(c, \nabla_0 u) = \psi_{\text{ch}}(c) + \psi_{\text{el}}(c, \nabla_0 u)
\]

3. **Chemical Potential**:
\[
\mu = \frac{\partial \psi}{\partial c} = -F U_{\text{OCV}} - \frac{v_{\text{pm}}}{3 \lambda_{\text{ch}}^5} (F^T F) : C[E_{\text{el}}]
\]

4. **Obstacle Contact Conditions**:
\[
\begin{cases}
u \leq g & \text{on } (0, t_{\text{end}}) \times \partial \Omega, \\
\sigma \cdot n \leq 0 & \text{on } (0, t_{\text{end}}) \times \partial \Omega, \\
[\sigma \cdot n] [u - g] = 0 & \text{on } (0, t_{\text{end}}) \times \partial \Omega.
\end{cases}
\]

The introduction of these formulas and models enables researchers to more accurately predict and understand the behavior of silicon anode materials under actual battery working conditions, thus providing theoretical support for the design of more durable and efficient lithium - ion batteries.