Using a three-dimensional computational fluid dynamics model and an agglomerate model to investigate the effect of varying agglomerate parameters and output voltages on proton exchange membrane fuel cell performance

Abstract

Modeling of proton exchange membrane (PEM) fuel cells is attracting more attention as fuel cell technology continues to develop. In this study, we considered a hybrid model that combines an agglomerate model based on the agglomeration of catalyst particles and the coverage-dependent kinetic equation of platinum oxide for ORR, and another 3D numerical model of a PEM fuel cell based on computational fluid dynamics (CFD). The obtained results from our developed models were validated with experimental results from literature. In fact, we investigated the effects of changing the agglomerate radius <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mrow><mml:mo>(</mml:mo> <mml:msub><mml:mi>R</mml:mi> <mml:mrow><mml:mi>a</mml:mi> <mml:mi>g</mml:mi> <mml:mi>g</mml:mi></mml:mrow> </mml:msub> <mml:mo>)</mml:mo></mml:mrow> </mml:math> , the ionomer volume fraction within the agglomerate <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow><mml:mrow><mml:mo>(</mml:mo> <mml:msub><mml:mi>Y</mml:mi> <mml:mrow><mml:mi>i</mml:mi> <mml:mo>,</mml:mo> <mml:mi>a</mml:mi> <mml:mi>g</mml:mi> <mml:mi>g</mml:mi></mml:mrow> </mml:msub> <mml:mo>)</mml:mo></mml:mrow> <mml:mtext>,</mml:mtext></mml:mrow> </mml:math> the effective agglomerate surface area <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mrow><mml:mo>(</mml:mo> <mml:msub><mml:mi>A</mml:mi> <mml:mrow><mml:mi>i</mml:mi> <mml:mo>,</mml:mo> <mml:mi>a</mml:mi> <mml:mi>g</mml:mi> <mml:mi>g</mml:mi></mml:mrow> </mml:msub> <mml:mo>)</mml:mo></mml:mrow> </mml:math> , the distribution of the gases and the temperature on the cell performances. The results revealed that the cell performances are strongly influenced by changing <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow><mml:msub><mml:mi>R</mml:mi> <mml:mrow><mml:mi>a</mml:mi> <mml:mi>g</mml:mi> <mml:mi>g</mml:mi></mml:mrow> </mml:msub> </mml:mrow> </mml:math> and <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow><mml:msub><mml:mi>Y</mml:mi> <mml:mrow><mml:mi>i</mml:mi> <mml:mo>,</mml:mo> <mml:mi>a</mml:mi> <mml:mi>g</mml:mi> <mml:mi>g</mml:mi></mml:mrow> </mml:msub> </mml:mrow> </mml:math> for medium and high current densities: The activation loss increases with increasing <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow><mml:msub><mml:mi>R</mml:mi> <mml:mrow><mml:mi>a</mml:mi> <mml:mi>g</mml:mi> <mml:mi>g</mml:mi></mml:mrow> </mml:msub> </mml:mrow> </mml:math> and decreasing <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow><mml:msub><mml:mi>Y</mml:mi> <mml:mrow><mml:mi>i</mml:mi> <mml:mo>,</mml:mo> <mml:mi>a</mml:mi> <mml:mi>g</mml:mi> <mml:mi>g</mml:mi></mml:mrow> </mml:msub> </mml:mrow> </mml:math> . Also, <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow><mml:msub><mml:mi>A</mml:mi> <mml:mrow><mml:mi>i</mml:mi> <mml:mo>,</mml:mo> <mml:mi>a</mml:mi> <mml:mi>g</mml:mi> <mml:mi>g</mml:mi></mml:mrow> </mml:msub> </mml:mrow> </mml:math> increases with decreasing <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow><mml:msub><mml:mi>R</mml:mi> <mml:mrow><mml:mi>a</mml:mi> <mml:mi>g</mml:mi> <mml:mi>g</mml:mi></mml:mrow> </mml:msub> </mml:mrow> </mml:math> and increasing <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow><mml:msub><mml:mi>Y</mml:mi> <mml:mrow><mml:mi>i</mml:mi> <mml:mo>,</mml:mo> <mml:mi>a</mml:mi> <mml:mi>g</mml:mi> <mml:mi>g</mml:mi></mml:mrow> </mml:msub> </mml:mrow> </mml:math> . In addition, the PEM fuel cell's power output is significantly enhanced when <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow><mml:msub><mml:mi>R</mml:mi> <mml:mrow><mml:mi>a</mml:mi> <mml:mi>g</mml:mi> <mml:mi>g</mml:mi></mml:mrow> </mml:msub> </mml:mrow> </mml:math> is decreased and <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow><mml:msub><mml:mi>Y</mml:mi> <mml:mrow><mml:mi>i</mml:mi> <mml:mo>,</mml:mo> <mml:mi>a</mml:mi> <mml:mi>g</mml:mi> <mml:mi>g</mml:mi></mml:mrow> </mml:msub> </mml:mrow> </mml:math> is increased, the optimal power being obtained for values of <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow><mml:msub><mml:mi>R</mml:mi> <mml:mrow><mml:mi>a</mml:mi> <mml:mi>g</mml:mi> <mml:mi>g</mml:mi></mml:mrow> </mml:msub> <mml:mo>=</mml:mo> <mml:mn>100</mml:mn> <mml:mi>n</mml:mi> <mml:mi>m</mml:mi></mml:mrow> </mml:math> and <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow><mml:msub><mml:mi>Y</mml:mi> <mml:mrow><mml:mi>i</mml:mi> <mml:mo>,</mml:mo> <mml:mi>a</mml:mi> <mml:mi>g</mml:mi> <mml:mi>g</mml:mi></mml:mrow> </mml:msub> </mml:mrow> </mml:math> = 0.6. The numerical results also showed that decreasing the output voltage from 0.95V to 0.35V can accelerate the electrochemical reaction process.