Open Access Article
This Open Access Article is licensed under a
Creative Commons Attribution 3.0 Unported Licence

Correction: Nafion film transport properties in a low-Pt PEM fuel cell: impedance spectroscopy study

Tatyana Reshetenkoa and Andrei Kulikovsky*b
aHawaii Natural Energy Institute, University of Hawaii, Honolulu, Hawaii 96822, USA. E-mail: tatyanar@hawaii.edu
bTheory and Computation of Energy Materials (IEK-13), Forschungszentrum Jülich GmbH, Institute of Energy and Climate Research, D-52425 Jülich, Germany. E-mail: A.Kulikovsky@fz-juelich.de

Received 22nd January 2021 , Accepted 22nd January 2021

First published on 8th February 2021


Abstract

Correction for ‘Nafion film transport properties in a low-Pt PEM fuel cell: impedance spectroscopy study’ by Tatyana Reshetenko et al., RSC Adv., 2019, 9, 38797–38806, DOI: 10.1039/C9RA07794D.


The authors would like to propose a modified version of the model presented in their article, which is consistent with the idea of the 1d + 1d approach reported in ref. 1. This adjustment makes the model more sensitive to the Nafion film parameters.

In this Correction article, equation numbers and notations referred to in the text are those used in the original article. On the right side of eqn (3) and its dimensionless version, eqn (12), the oxygen flux is calculated at the metal surface. A more accurate approximation is to use the flux NN,p at the open pore radius Rp instead of NN,m, i.e.,

 
image file: d1ra90053f-t1.tif(1)
and its dimensionless variant should be substituted into eqn (3) and (12), respectively. It follows that in eqn (23) the perturbed flux image file: d1ra90053f-t2.tif should be replaced by image file: d1ra90053f-t3.tif and eqn (23) takes the form
 
image file: d1ra90053f-t4.tif(2)

The flux image file: d1ra90053f-t5.tif is obtained from solving eqn (24), resulting in

 
image file: d1ra90053f-t6.tif(3)
where
 
image file: d1ra90053f-t7.tif(4)
 
image file: d1ra90053f-t8.tif(5)
 
image file: d1ra90053f-t9.tif(6)
and
 
image file: d1ra90053f-t10.tif(7)
Here, J0, J1 are the Bessel functions of the first kind, and K0, K1 are the modified Bessel functions of the second kind. This modification leads to marginal changes in the calculation results reported in the original article. In the range of current densities 100 to 800 mA cm−2, the Nafion film thickness changes between 8.5 and 12 nm and the film resistivity decreases from 0.6 s cm−1 down to 0.25 s cm−1, in agreement with the results in the original article.

Qualitatively, the pore radius is about an order of magnitude larger than the Nafion film thickness, meaning that the shape of the small film element is close to planar and hence the real part of the fluxes image file: d1ra90053f-t11.tif and image file: d1ra90053f-t12.tif are almost equal. Eqn (23) takes into account the phase shift due to transport of image file: d1ra90053f-t13.tif perturbation through the film. However, replacing NN,m by NN,p in this equation increases the effect of the film on the phase shift. In general, this makes the model more sensitive to the film parameters.

The Royal Society of Chemistry apologises for these errors and any consequent inconvenience to authors and readers.

References

  1. A. A. Kulikovsky, Analytical Modelling of Fuel Cells, Elsevier, Amsterdam, 2nd edn, 2019 Search PubMed .

This journal is © The Royal Society of Chemistry 2021