Simulation and Visualization of Cell Performance Distribution of SOFCs
Authors: Takino Keisuke , Tachikawa Yuya , Shiratori Yusuke , Taniguchi Shunsuke , Sasaki Kazunari 
Published: February 2020

Introduction

It has been reported for SOFC systems that high electrical efficiency beyond 70% could be achieved theoretically in a recent study[1]. It is important to understand the internal characteristics to design highly efficient SOFC. However, it is difficult to visualize internal reactions experimentally. From these reasons, it is desired to develop a simulation technique for visualizing various phenomena within fuel cell stacks. Here, we focus on numerical analysis on electrode kinetics. Exchange current density is an important phenomenological parameter in considering electrochemical reaction kinetics at SOFC electrodes and in simulating cell performance. Here, the exchange current density of fully pre-reformed methane-fueled SOFC anodes is measured and phenomenological equations of anode exchange current density are derived. Furthermore, we combine numerical analysis with the phenomenological relations to simulate the spatial distribution of anode overvoltages.

Experimental conditions

In this study, electrolyte-supported cells were used to measure the anode exchange current density. The electrolyte material was ScSZ (ScSZ:10mol% Sc2O3–1mol% CeO2–89mol% ZrO2) and the electrode prepared using NiO-ScSZ was applied to the anode. The electrode prepared using LSM（LSM: (La0.8Sr0.2)0.98MnO3）and LSM-ScSZ was applied to the cathode. By changing operation temperature, fuel utilization (Uf), and steam-to-carbon ratio (S/C), anode overvoltage was measured at each current density (0.25, 0.3, 0.35, and 0.4 cm2) for methane-fueled pre-reformed SOFCs by the current interrupt method. The anode exchange current density was then derived from the current density and anode overvoltage values. In simulation, three-dimensional cell model was created using a numerical analysis software. A parallel co-flow model was considered. The inlet gas temperature was 800°C and the anode gas composition was fixed with the S/C of 2.5.

Results and discussion

The equation of anode exchange current density based on hydrogen-fueled model[2] was revised to consider the equilibrium reforming reaction, the water gas shift reaction, and the conservation law of anode gas partial pressure. It has been reported that the electrochemical reaction, where the adsorbed hydrogen atom (Had) in hydrogen-based fuels is rate-limiting[2], an Arrhenius-type equation shown in Eq. 1 has been proposed based on the assumption that the electrochemical reaction of Had is rate- limiting:

i
0,a = a·exp(b/RT)·PH2(A)

0.5·PH2O(A)

0.5 (1)
The partial pressure terms were modified considering the equilibrium reforming reaction and the water gas shift reaction. Pre-exponential factor a and activation energy b were fixed using the experimental results. Fig. 1 shows the experimental results of anode exchange current density. The phenomenological values derived from Eq. 1 are shown in Fig. 1. As the fuel utilization varied, the exchange current density increased and reached the maximum at around 40% of Uf. The experimental results could be well fitted by Eq.1. The SOFC cell performance simulation was performed using the three-dimensional SOFC model. Fig. 2 shows the anode overvoltage distribution along the fuel flow direction from the left side to the right side. The overvoltage decreased from the inlet towards the outlet. This is partly because oxygen partial pressure on the anode increased to the outlet and the electromotive force and the exchange current decreased.

References

[1] Y. Matsuzaki et al., Sci. Rep., 5, 12640 (2015).
[2] T. Hosoi et al., J. Electrochem. Soc., 162 (1), F136-F152 (2015).

Figure 1