A computational approach for the simulation of natural convection in electrochemical cells

A novel computational approach for the numerical simulation of electrochemical systems influenced by natural convection phenomena is presented. A stabilized finite element framework for multi-ion transport mechanisms including convection, diffusion and migration coupled to an incompressible flow solver is developed. The role of a galvanostatic Butler-Volmer condition including the interaction of ionic concentration at the surface of the electrode and the surface overpotential is emphasized, to obtain a non-uniform surface overpotential distribution. Additionally, a three-dimensional rotationally-symmetric boundary condition is used for modeling rotating cylinder electrodes. The computational framework is tested for various numerical examples exhibiting two- and three-dimensional electrochemical cell configurations including dilute CuSO4 electrolyte solutions with and without excess of supporting H2SO4 electrolyte.

[1]  Y. Kondo,et al.  Concentration Profile of Cu2+ Ion near a Plane Vertical Cathode in Electrolytes Containing CuSO4 and an Excess of H2SO4 as a Supporting Electrolyte , 1983 .

[2]  Wolfgang A. Wall,et al.  Residual‐based variational multiscale methods for laminar, transitional and turbulent variable‐density flow at low Mach number , 2011 .

[3]  A. D. Davydov,et al.  Simulation of ion transfer under conditions of natural convection by the finite difference method , 2003 .

[4]  D. Bograchev,et al.  Determination of Mass Coefficients of Ions in a Quantitative Analysis of the Effect of Natural Convection on Electrochemical Processes , 2005 .

[5]  C. Wagner.  The Role of Natural Convection in Electrolytic Processes , 1949 .

[6]  D. Lincot,et al.  Numerical Simulation of Hydrodynamic and Mass-Transport Properties at a Laminar Rotating Cylinder Electrode , 2006 .

[7]  T. M. Brown,et al.  By Electrochemical methods , 2007 .

[8]  A. Berkovich,et al.  A COMPUTER PROOF OF A POLYNOMIAL IDENTITY IMPLYING A PARTITION THEOREM OF G (cid:127) OLLNITZ , 2022 .

[9]  Wolfgang A. Wall,et al.  An algebraic variational multiscale-multigrid method for large-eddy simulation of turbulent variable-density flow at low Mach number , 2010, J. Comput. Phys..

[10]  S. Kida,et al.  Numerical Calculation of Transient Current Density Distribution along Vertical Plane Electrode in CuSO4 – H2SO4 Electrolyte Solution , 2010 .

[11]  Kenneth E. Jansen,et al.  A stabilized finite element method for the incompressible Navier–Stokes equations using a hierarchical basis , 2001 .

[12]  Yousef Saad,et al.  Iterative methods for sparse linear systems , 2003 .

[13]  S. Kida,et al.  Numerical Simulation of Ionic Mass-Transfer Rates with Natural Convection in CuSO4-H2SO4 Solution , 2009 .

[14]  John Newman,et al.  Electrochemical Systems, 3rd Edition , 2004 .

[15]  M. Moats,et al.  The effect of copper, acid, and temperature on the diffusion coefficient of cupric ions in simulated electrorefining electrolytes , 2000 .

[16]  D. Lincot,et al.  Prediction of the electrodeposition process behavior with the gravity or acceleration value at continuous and discrete scale , 2007 .

[17]  Tayfan E. Tezduyar,et al.  Stabilized Finite Element Formulations for Incompressible Flow Computations , 1991 .

[18]  W. Wall,et al.  A stabilized finite element method for the numerical simulation of multi-ion transport in electrochemical systems , 2012 .

[19]  S. Kida,et al.  Numerical Simulation of Ionic Mass-Transfer Rates with Natural Convection in CuSO4 – H2SO4 Solution I. Numerical Study on the Developments of Secondary Flow and Electrolyte Stratification Phenomena , 2009 .

[20]  Paul T. Lin,et al.  Performance of a parallel algebraic multilevel preconditioner for stabilized finite element semiconductor device modeling , 2009, J. Comput. Phys..

[21]  M. Chung.  A numerical method for analysis of tertiary current distribution in unsteady natural convection multi-ion electrodeposition , 2000 .

[22]  Tom De Mulder,et al.  The role of bulk viscosity in stabilized finite element formulations for incompressible flow: A review , 1998 .

[23]  Thomas J. R. Hughes,et al.  Finite element modeling of blood flow in arteries , 1998 .

[24]  T. Hughes,et al.  Streamline upwind/Petrov-Galerkin formulations for convection dominated flows with particular emphasis on the incompressible Navier-Stokes equations , 1990 .

[25]  J. Selman,et al.  Free‐Convection Mass Transfer with a Supporting Electrolyte , 1971 .

[26]  Y. Kondo,et al.  Studies on the velocity profile in natural convection during copper deposition at vertical cathodes , 1976 .

[27]  T. N. Stevenson,et al.  Fluid Mechanics , 2021, Nature.

[28]  Ray Bert.  Book Review: \IIntegrated Geospatial Technologies: A Guide to GPS, GIS, and Data Logging\N by Jeff THurston, Thomas K. Poiker, and J. Patrick Moore. Hoboken, New Jersey: John Wiley & Sons, Inc., 2003 , 2004 .

[29]  S. Kida,et al.  Transient mass transfer rate of Cu2+ ion caused by copper electrodeposition with alternating electrolytic current , 2010 .

[30]  C. Wallgren,et al.  Mass transport in a weakly stratified electrochemical cell , 1996 .

[32]  Thomas J. R. Hughes,et al.  Multiscale and Stabilized Methods , 2007 .

[33]  K. Eckert,et al.  The concentration field during transient natural convection between vertical electrodes in a small-aspect-ratio cell , 2008 .