Asymptotic analysis on charging dynamics for stack-electrode model of supercapacitors

Lijie Ji,Zhenli Xu,Shenggao Zhou

DOI: https://doi.org/10.1098/rspa.2023.0044

2023-04-21

Abstract:Supercapacitors are promising electrochemical energy storage devices due to their prominent performance in rapid charging/discharging rates, long cycle life, stability, etc. Experimental measurement and theoretical prediction on charging timescale for supercapacitors often have large difference. This work develops a matched asymptotic expansion method to derive the charging dynamics of supercapacitors with porous electrodes, in which the supercapacitors are described by the stack-electrode model. Coupling leading-order solutions between every two stacks by continuity of ionic concentration and fluxes leads to an ODE system, which is a generalized equivalent circuit model for zeta potentials, with the potential-dependent nonlinear capacitance and resistance determined by physical parameters of electrolytes, e.g., specific counterion valences for asymmetric electrolytes. Linearized stability analysis on the ODE system after projection is developed to theoretically characterize the charging timescale. The derived asymptotic solutions are numerically verified. Further numerical investigations on the biexponential charging timescales demonstrate that the proposed generalized equivalent circuit model, as well as companion linearized stability analysis, can faithfully capture the charging dynamics of symmetric/asymmetric electrolytes in supercapacitors with porous electrodes.

Numerical Analysis,Mathematical Physics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is the significant difference between the theoretical prediction and experimental measurement of the charging time scale of supercapacitors. Specifically, the author focuses on the dynamic behavior of supercapacitors with porous electrode structures during the charging process. By introducing a matched asymptotic expansion method, the author aims to derive a mathematical model to describe the charging dynamics of such complex systems.

### Main problems and solutions

1. **Problem background**:
- Supercapacitors have become promising electrochemical energy - storage devices due to their advantages such as fast charge - discharge, long cycle life, and stability.
- However, there is a large difference between the theoretically predicted charging time scale and the actual experimental measurement, especially when considering the porous electrode structure.

2. **Research objectives**:
- Develop a matched asymptotic expansion method to analyze the charging dynamics of supercapacitors with porous electrodes.
- Construct a generalized equivalent circuit model by coupling the continuity conditions of ion concentration and flux between every two stacks.
- Use linearized stability analysis to theoretically characterize the charging time scale.

3. **Specific methods**:
- **Model description**: Use the stack - electrode model to describe the supercapacitor and use the Poisson - Nernst - Planck (PNP) equations to describe ion transport and potential distribution.
- **Asymptotic expansion**: Based on the small parameter \(\epsilon\) (the ratio of the Debye length to the electrode spacing), perform a matched asymptotic analysis to obtain the leading - order solutions of ion concentration, potential distribution, and total diffusive charge.
- **Equivalent circuit model**: Couple the leading - order solutions through the continuity conditions to form a system of ordinary differential equations (ODEs) for the zeta potential. This system can be regarded as a generalized equivalent circuit model, where the nonlinear capacitance and resistance are determined by the physical parameters of the electrolyte.
- **Stability analysis**: Perform a linearized stability analysis on the projected ODE system and solve the eigenvalue problem to determine the charging time scale.

4. **Verification and application**:
- Numerically verify the accuracy of the asymptotic solutions.
- Further numerically study the double - exponential charging time scale and prove that the proposed generalized equivalent circuit model can accurately capture the charging dynamics of symmetric and asymmetric electrolyte systems.

### Formula summary

- **PNP equations**:
\[
\begin{cases}
\frac{\partial c_i}{\partial t}+\nabla\cdot J_i = 0, & i=\pm\\
-\epsilon_0\nabla\cdot(\epsilon\nabla\Phi)=\sum_i z_i e c_i
\end{cases}
\]
where \(J_i=-D_i\left(\nabla c_i+\beta c_i\nabla(z_i e\Phi)\right)\), \(\beta=\frac{1}{k_B T}\).

- **Leading - order solutions**:
\[
\begin{cases}
c_{\pm,0}(x,t)=|z_\mp|\\
\Phi_0(x,t)=j_0(t)x + A(t)
\end{cases}
\]

- **ODE system of the generalized equivalent circuit model**:
\[
C(\zeta)\dot{\zeta}=\alpha(-T\zeta + Y)
\]
where \(\zeta = [\zeta_{n,L},\cdots,\zeta_{1,L},\zeta_{1',R},\cdots\)