Formulation and analysis of a DPG discretization for a simplified electrochemical model

Jaime Mora-Paz

2024-07-10

Abstract:We present a simplified model consisting on two linear elliptic boundary-value problems that represent a single step and single fixed-point iteration in an electrochemical battery model. The main variables are the concentration and the electric potential, whose equation is assigned a Robin BC with a very important physical interpretation. The solvability of both equations is studied in different funcional settings, to finally prove the well-posedness of a broken mixed variational formulation. The latter formulation opens the opportunity of performing discretization and numerical solution via the Discontinuous Petrov-Galerkin (DPG) method, which guarantees discrete stability thanks to optimal test functions. With only the usual assumptions on the data and the discretization, we show that the method herein proposed is convergent. This analytical effort complements other recent works on batttery multiphysics, that have been more modeling and computing oriented, and establishes the DPG approach as a valuable tool to contribute toward an efficient analysis and design of batteries.

Numerical Analysis,Mathematical Physics

What problem does this paper attempt to address?

### Problems the paper attempts to solve

This paper aims to solve the problems of discretization and numerical solution of two linear elliptic boundary - value problems in simplified electrochemical models. Specifically, it focuses on the modeling and analysis of the concentration field \( c \) and the electric potential field \( \phi \) in the microstructure of lithium - ion batteries.

#### 1. **Concentration problem**
- **Description**: The concentration problem involves the volume concentration \( c: \Omega \to \mathbb{R} \) of species in the medium, where \( \Omega \subset \mathbb{R}^d \) (\( d \in \{2,3\} \)). This problem can be expressed as a reaction - diffusion partial differential equation (PDE) with Neumann boundary conditions.
- **Equation**:
\[
\begin{cases}
c-\Delta t D \nabla\cdot\nabla c = c_{\text{prev}} & \text{in } \Omega\\
(-D \nabla c)\cdot n = J & \text{on } \partial\Omega
\end{cases}
\]
- **Physical meaning**: This equation describes the diffusion and reaction processes of species in the medium, as well as the flux on the boundary.

#### 2. **Electric potential problem**
- **Description**: The electric potential problem involves the electric potential field \( \phi: \Omega \to \mathbb{R} \), which takes the form of Poisson's equation or Laplace's equation, with mixed boundary conditions, especially Robin boundary conditions.
- **Equation**:
\[
\begin{cases}
-\kappa \nabla\cdot\nabla \phi=\nabla\cdot S & \text{in } \Omega\\
\phi = 0 & \text{on } \Gamma_D\\
(-\kappa \nabla \phi)\cdot n = I & \text{on } \Gamma_N\\
(-\kappa \nabla \phi)\cdot n-\beta \phi = R & \text{on } \Gamma_R
\end{cases}
\]
- **Physical meaning**: The electric potential equation describes the distribution of the electric field in the medium. In particular, through the Robin boundary condition, it relates the current density and the electric potential.

#### 3. **Main research contents**
- **Well - posedness of variational forms**: The author studied the solvability of these two problems in different function spaces and proved the well - posedness of a broken mixed variational form. This provides a theoretical basis for discretization and numerical solution using the discontinuous Petrov - Galerkin (DPG) method.
- **Application of the DPG method**: The DPG method ensures discrete stability by choosing optimal test functions, making this method efficient and accurate in dealing with complex electrochemical models.

#### 4. **Motivation and background**
- **Electrochemical model**: These problems are proposed based on the electrochemical model of the microstructure of lithium - ion batteries. This model considers the interactions among three materials: the negative electrode, the positive electrode, and the electrolyte, as well as the spatio - temporal distribution and evolution of the concentration field and the electric potential field among them.
- **Butler - Volmer reaction kinetics**: The electrochemical reaction on the electrode - electrolyte interface is the main mechanism of the coupled physical processes, and the Butler - Volmer formula describes the relationship between the current density and the concentration and the electric potential.

### Summary
By introducing the DPG method, this paper solves the problems of discretization and numerical solution of the concentration and electric potential problems in the simplified electrochemical model, providing new tools and methods for understanding and designing efficient battery models.