Mathematical modelling of Solid Oxide Fuel Cells revisited -- a modified formulation of the problem

Michal Wrobel,Grzegorz Brus

DOI: https://doi.org/10.48550/arXiv.1901.05509

2021-04-20

Abstract:In the paper a mathematical model of the PEN structure (positive electrode - electrolyte - negative electrode) of the Solid Oxide Fuel Cell (SOFC) is analyzed. It is proved that classical formulation of the problem leads inevitably to locally unphysical effects related to negative values of the activation overpotential. Moreover, the active layers' thicknesses are shown to be components of solution and cannot be predefined in an arbitrary way. A modified mathematical formulation of the problem is proposed which includes this novelty alongside an amended definition of the activation overpotential. A dedicated computational scheme is developed for the cathode sub-problem. The accuracy of computations is investigated by means of a newly introduced analytical benchmark example. The numerical results obtained for LSCF cathode are used to discuss certain aspects of the modified formulation and the active layer thickness. The new modelling approach is validated by comparison with experimental data.

Numerical Analysis

What problem does this paper attempt to address?

The problem that this paper attempts to solve is: **The problem of local non - physical effects in the classical solid oxide fuel cell (SOFC) mathematical model, especially the problems related to negative activation overpotential**. Specifically, the author points out that the classical model will produce unreasonable negative activation overpotential and negative charge - transfer current near the boundary of the electrode active layer, which leads to physically infeasible results. In addition, the thickness of the active layer is usually preset in the classical model rather than being part of the solution, which is also inconsistent with the actual situation.

To solve these problems, the author proposes an improved mathematical model. The main improvements include:
1. **Revising the definition of activation overpotential** to avoid negative values.
2. **Taking the thickness of the active layer as part of the solution** instead of a preset parameter, thus more accurately reflecting the actual physical phenomena.

Through these improvements, the author hopes to eliminate local non - physical effects and improve the accuracy and reliability of the model. The improved model has also been verified by numerical calculations and experimental data.

### Summary of Key Formulas

1. **Charge - transfer current equation**:
\[
\frac{d}{dy}\left(\sigma_{el} \frac{d\varphi_{el}}{dy}\right) =
\begin{cases}
0 & \text{for } y \in (-h_3, -h_a) \cup (-h_1, h_1) \cup (h_b, h_2), \\
i_{ct} & \text{for } y \in (-h_a, -h_1), \\
-i_{ct} & \text{for } y \in (h_1, h_b)
\end{cases}
\]
\[
\frac{d}{dy}\left(\sigma_{ion} \frac{d\varphi_{ion}}{dy}\right) =
\begin{cases}
0 & \text{for } y \in (-h_3, -h_a) \cup (-h_1, h_1) \cup (h_b, h_2), \\
-i_{ct} & \text{for } y \in (-h_a, -h_1), \\
i_{ct} & \text{for } y \in (h_1, h_b)
\end{cases}
\]

2. **Specific expression of charge - transfer current**:
\[
i_{ct}^{ano} = i_{tpb,0}^{ano} l_{tpb}^{ano} \left[\exp\left(\frac{2F\eta_{act}^{ano}}{RT}\right) - \exp\left(-\frac{F\eta_{act}^{ano}}{RT}\right)\right]
\]
\[
i_{ct}^{cat} = i_{dpb,0}^{cat} A_{dph}^{cat} \left[\exp\left(\frac{1.2F\eta_{act}^{cat}}{RT}\right) - \exp\left(-\frac{F\eta_{act}^{cat}}{RT}\right)\right]
\]

3. **Definition of activation overpotential**:
\[
\eta_{act}^{ano} = \varphi_{el} - \varphi_{ion} - \eta_{conc}^{ano}
\]
\[
\eta_{act}^{cat} = \varphi_{ion} - \varphi_{el} - \eta_{conc}^{cat}
\]

4. **Definition of concentration overpotential**:
\[