Natural Convection in Redox Electrochemistry: Model, Simulations and Experiments

A natural convection model has been implemented on redox species TMPD with supporting electrolyte TBAP in acetonitrile solvent. The cross section geometry used for the two-dimensional simulations is 20 mm (W) x 10 mm (H). The working electrode is a 2 mm wide strip at the center of the cell floor and the counter electrode is the entire cell ceiling. For the three-dimensional simulations, a 3mm diameter working electrode in the semi-infinite configuration was used. A uniform magnetic field was applied in the direction perpendicular the working electrode surface. The first case we have studied has an initial concentration of 10.3 mM TMPD - and no oxidized form of TMPD. The supporting electrolyte concentration was 0.5 M. The cell was operated in the potential step mode and the applied potential was high enough for operation in the diffusion-limited regime. Simulations were run for two values of the magnetic field intensity. The electrical conductivity of the bulk solution was 0.625 S/m. The density was calculated using the formulation in our natural convection model based on electroneutrality, where non-electroactive ions migrate to the diffusion layer to neutralize the charge buildup due to the redox reaction. Cases with gravity not present, and gravity applied in the +x, -x, +y and -y directions in 2D were investigated. For the 3D simulations, gravity was applied in the direction normal to the electrode, toward and away from the electrode surface. To study the effect of pure natural convection, cases were run without applying the magnetic field. For this electrochemical system, the density changes from the initial bulk density (785 kg/m 3 ) is very small with a maximum value of 0.0231 kg/m 3 , which gives a difference of 0.00294% between the highest and the lowest density. Even under this small difference, the evolution of the flow field is quite remarkable. With increasing time from the start, the effect of natural convection on the flow field becomes greater. The maximum velocity increased by an order of magnitude in ~30s from the start. At larger times, the maximum velocity tends toward an asymptotic value along with the electrode current also leveling off and reaching an asymptotic value. The behavior of mixed convection characterized by competing forces of buoyancy and Lorentz force was also investigated, where two values of the magnetic field were chosen. The natural convection-only cases show that the effect is stronger when the established flow is aw