A 2D numerical model is developed for a laminar flow fuel cell considering ion transport and the electric double layer around the electrodes. The Frumkin-Butler-Volmer equation is used for the fuel cell kinetics. The finite volume method is used to form algebraic equations from governing partial differential equations. The numerical solution was obtained using Newton’s method and a block TDMA solver. The model accounts for the coupling of charged ion transport with the electric field and is able to fully resolve the diffuse regions of the electric double layer in both the stream-wise and cross-channel directions. Different operating phenomena, such as laminar flow separation and the development of the depletion boundary layers and electric double layers are obtained. These numerical results demonstrate the model’s ability to capture the complex behavior of a microfluidic fuel cell which has been ignored in previous 1D models.
