Numerical Solutions By Method of Lines Approach for Fluid Flow in a Modified Rotating Disk Electrode Apparatus

Abstract

Studying the reaction kinetics of the electrochemical redox reactions on a rotating disk electrode (RDE) is a well-known technique to access kinetic information under various competing levels of kinetic and diffusion dominated conditions. Reactions of flow batteries are typically carried out on graphite felt electrode surfaces at both cell and stack level. But the kinetics of the both positive (VO2 +/VO2+) and negative (V3+/V2+) redox couple reactions are not well understood, as there is no clear agreement of redox kinetics information available in literature1-2for these electrolytes. Porous rotating disk electrode (PRDE) is a commonly used technique to study the kinetics of graphite felt electrodes. A semi-infinite condition is typically used to arrive at well-known Levich profiles for velocity and limiting current analysis for rotating disk electrodes3. However, this type of analysis didn’t account for the data observed for felt electrodes in a RDE set up. To obtain a numerical solution for fluid flow around a rotating disk electrode in a modified RDE setup is proposed in this work as shown in Figure 1. In this setup, the boundary conditions at the stationary walls, rotating walls and axial symmetry are well-defined. This method can be used for obtaining primary and secondary distributions effectively in electrochemical systems efficiently and also for simulating flow problems4. Numerical solutions with well-defined boundary conditions will help in enhancing the reliability of the simulated results. In this 2D fluid flow model, an axially symmetric cylindrical geometry is used and fluid flow was assumed to vary in radial (r) and axial (z) directions, whereas azimuthal (θ) dependence of fluid flow is ignored. Three components of Navier-Stokes equations (r, θ and z) and a continuity equation represent the fluid flow outside the rotating disk electrode. The velocities (ur, u θ and uz ), pressure and radial and axial distances were non-dimensionalized through a method, using fluid properties and angular velocity, as reported in literature5. The 2D model equations were simulated using a new technique in which numerical method is applied in the z-direction and the resulting equations are integrated using efficient Boundary value problem (BVP) solvers in r. Efficiency and accuracy of the proposed method compared with different finite difference and discretization schemes in r will be presented. Acknowledgement: This work was supported by a contract from The American Chemical Society- Petroleum Research Fund (ACS-PRF) Grant PRF# 52588-ND10.