A 3D finite element approach for the coupled numerical simulation of electrochemical systems and fluid flow

A comprehensive finite element method for three‐dimensional simulations of stationary and transient electrochemical systems including all multi‐ion transport mechanisms (convection, diffusion and migration) is presented. In addition, non‐linear phenomenological electrode kinetics boundary conditions are accounted for. The governing equations form a set of coupled non‐linear partial differential equations subject to an algebraic constraint due to the electroneutrality condition. The advantage of a convective formulation of the ion‐transport equations with respect to a natural application of homogeneous flux boundary conditions is emphasized. For one of the numerical examples, an analytical solution for the coupled problem is provided, and it is demonstrated that the proposed computational approach is robust and provides accurate results. Copyright © 2011 John Wiley & Sons, Ltd.