Entropy Increasing Numerical Methods for Prediction of Reversible and Irreversible Heating in Supercapacitors

Jie Ding,Xiang Ji,Shenggao Zhou

2024-09-15

Abstract:Accurate characterization of entropy plays a pivotal role in capturing reversible and irreversible heating in supercapacitors during charging/discharging cycles. However, numerical methods that can faithfully capture entropy variation in supercapacitors are still in lack. This work develops first-order and second-order finite-volume schemes for the prediction of non-isothermal electrokinetics in supercapacitors. Semi-implicit discretization that decouples temperature from ionic concentrations and electric potential results in an efficient first-order accurate scheme. Its numerical analysis theoretically establishes the unique solvability of the nonlinear scheme with the existence of positive ionic concentrations and temperature at discrete level. To obtain an entropy-increasing second-order scheme, a modified Crank-Nicolson approach is proposed for discretization of the logarithm of both temperature and ionic concentrations, which is employed to enforce numerical positivity. Moreover, numerical analysis rigorously demonstrates that both first-order and second-order schemes are able to unconditionally preserve ionic mass conservation and original entropy increase for a closed, thermally insulated supercapacitor. Extensive numerical simulations show that the proposed schemes have expected accuracy and robust performance in preserving the desired properties. Temperature oscillation in the charging/discharging processes is successfully predicted, unraveling a quadratic scaling law of temperature rising slope against voltage scanning rate. It is also demonstrated that the variation of ionic entropy contribution, which is the underlying mechanism responsible for reversible heating, is faithfully captured. Our work provides a promising tool in predicting reversible and irreversible heating in supercapacitors.

Numerical Analysis

What problem does this paper attempt to address?

The problem that this paper attempts to solve is how to accurately predict reversible and irreversible heating phenomena during the charging and discharging process of supercapacitors. Specifically, the author points out that the accurate characterization of entropy is crucial for capturing reversible and irreversible heating in supercapacitors. However, at present, there is a lack of numerical methods that can faithfully reflect entropy changes. Therefore, in this paper, first - order and second - order finite - volume schemes for predicting non - isothermal electrokinetics are developed to solve this problem.

### Main contributions
1. **First - order and second - order finite - volume schemes**:
- Developed first - order and second - order finite - volume schemes for predicting non - isothermal electrokinetics in supercapacitors.
- The first - order scheme adopts a semi - implicit discretization method to decouple temperature from ion concentration and electric potential, thereby improving computational efficiency.
- The second - order scheme discretizes the logarithms of temperature and ion concentration by modifying the Crank - Nicolson method to maintain numerical positivity.

2. **Numerical analysis**:
- Theoretically proved the unique solvability of the nonlinear equations of the first - order scheme at the discrete level, and the existence of positive ion concentration and temperature.
- Proved that both the first - order and second - order schemes can unconditionally maintain ion mass conservation and the original entropy increase.

3. **Numerical simulation**:
- Extensive numerical simulations show that the proposed schemes have the expected accuracy and robust performance and can effectively retain the required physical properties.
- Successfully predicted the temperature oscillations during the charging and discharging process and revealed the quadratic scaling relationship between the temperature rise slope and the voltage scanning rate.
- Accurately captured the change in the ion entropy contribution, which is the basic mechanism of reversible heating.

### Key formulas
- **Entropy function**:
\[
S(V, t)=-\int_V\left(\sum_{l = 1}^M c_l\log c_l - C_T\log T\right)dx
\]
where \( c_l \) is the concentration of the \( l \) - th ion, \( T \) is the temperature, and \( C_T \) is a constant related to the ionic heat capacity.

- **Temperature evolution equation**:
\[
C_T k_B\frac{\partial T}{\partial t}=k\Delta T - k_B T\sum_{l = 1}^M c_l\nabla\cdot\mathbf{u}_l+\sum_{l = 1}^M\nu_l c_l|\mathbf{u}_l|^2
\]
where \( k \) is the thermal conductivity, \( \nu_l \) is the dynamic viscosity of the \( l \) - th ion, and \( \mathbf{u}_l \) is the velocity of the \( l \) - th ion.

- **Entropy increase theorem**:
\[
\frac{dS}{dt}=\int_V\left(k\left\|\frac{\nabla T}{T}\right\|^2+\sum_{l = 1}^M\nu_l c_l\frac{|\mathbf{u}_l|^2}{T}\right)dx\geq0
\]

### Conclusion
The method proposed in this paper provides a promising tool for predicting reversible and irreversible heating in supercapacitors. Through numerical simulation, the effectiveness and robustness of the proposed method are verified, and it can accurately capture the temperature oscillations and entropy changes during the charging and discharging process. These results are of great significance for the design of high - performance supercapacitors.