An isogeometric variational multiscale method for large-eddy simulation of coupled multi-ion transport in turbulent flow

Electrochemical processes, such as electroplating of large items in galvanic baths, are often coupled to turbulent flow. In this study, we propose an isogeometric residual-based variational multiscale finite element method for multi-ion transport in dilute electrolyte solutions under turbulent flow conditions. In other words, this means that the concepts of isogeometric discretization and variational multiscale methods are successfully combined for developing a method capable of simulating the challenging problem of coupled multi-ion transport in turbulent flow. We present a comprehensive three-dimensional computational method taking into account, among others, coupled convection-diffusion-migration equations subject to an electroneutrality constraint in combination with phenomenological electrode-kinetics modeling. The electrochemical subproblem is one-way coupled to turbulent incompressible flow via convection. Ionic mass transfer in turbulent Taylor-Couette flow is investigated, representing an important model problem for rotating-cylinder-electrode configurations. Multi-ion transport as considered here is an example for mass transport at high Schmidt number (Sc=1389). An isogeometric discretization is especially advantageous for the present problem, since (i) curved boundaries can be represented exactly, and (ii) it has been proven to provide very accurate solutions for flow quantities when being applied in combination with residual-based variational multiscale modeling. We demonstrate that the method is robust and provides results which are in good agreement with direct numerical simulation results as well as empirical mass-transfer correlations reported in literature.

[1]  F. Walsh,et al.  The rotating cylinder electrode: a review of development , 1983 .

[2]  Y. Saad,et al.  GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems , 1986 .

[3]  Thomas J. R. Hughes,et al.  Isogeometric Analysis: Toward Integration of CAD and FEA , 2009 .

[4]  D. Gabe.  The rotating cylinder electrode , 1974 .

[5]  T. Hughes,et al.  Variational multiscale residual-based turbulence modeling for large eddy simulation of incompressible flows , 2007 .

[6]  DNS of passive scalar transport in turbulent channel flow at high Schmidt numbers , 2007 .

[7]  Pavel B. Bochev,et al.  On the Finite Element Solution of the Pure Neumann Problem , 2005, SIAM Rev..

[8]  Frank C. Walsh,et al.  Numerical simulation of the current, potential and concentration distributions along the cathode of a rotating cylinder Hull cell , 2007 .

[9]  J. Talbot,et al.  A Model of Electrocodeposition on a Rotating Cylinder Electrode , 2007 .

[10]  Wolfgang A. Wall,et al.  A computational approach for the simulation of natural convection in electrochemical cells , 2013, J. Comput. Phys..

[11]  J. Deconinck,et al.  Multi-ion transport and reaction simulations in turbulent parallel plate flow , 2004 .

[12]  S. Collis,et al.  Monitoring unresolved scales in multiscale turbulence modeling , 2001 .

[13]  D. Lincot,et al.  Mean and unsteady hydrodynamic and mass transport properties at a rotating cylinder electrode: From laminar to transitional flow regime , 2006 .

[14]  Isabelle Calmet,et al.  Large-eddy simulation of high-Schmidt number mass transfer in a turbulent channel flow , 1997 .

[15]  Victor M. Calo,et al.  The role of continuity in residual-based variational multiscale modeling of turbulence , 2007 .

[16]  T. Hughes,et al.  Large Eddy Simulation and the variational multiscale method , 2000 .

[17]  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..

[18]  Gregory C. Burton,et al.  The nonlinear large-eddy simulation method applied to Sc≈1 and Sc⪢1 passive-scalar mixing , 2008 .

[19]  D. R. Gabe,et al.  The rotating cylinder electrode: its continued development and application , 1998 .

[20]  M. A. Patrick,et al.  Turbulent fluid flow and electrochemical mass transfer in an annular duct with an obstruction , 2009 .

[21]  I. Akkerman,et al.  Large eddy simulation of turbulent Taylor-Couette flow using isogeometric analysis and the residual-based variational multiscale method , 2010, J. Comput. Phys..

[22]  R. Lui,et al.  Global existence of solutions of a strongly coupled quasilinear parabolic system with applications to electrochemistry , 2003 .

[23]  T. Hughes.  Multiscale phenomena: Green's functions, the Dirichlet-to-Neumann formulation, subgrid scale models, bubbles and the origins of stabilized methods , 1995 .

[24]  Gerald Farin,et al.  NURBS: From Projective Geometry to Practical Use , 1999 .

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

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

[27]  Les A. Piegl,et al.  The NURBS Book , 1995, Monographs in Visual Communication.

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

[30]  G. Hulbert,et al.  A generalized-α method for integrating the filtered Navier–Stokes equations with a stabilized finite element method , 2000 .

[31]  Suchuan Dong,et al.  Direct numerical simulation of turbulent Taylor–Couette flow , 2007, Journal of Fluid Mechanics.

[32]  Volker Gravemeier,et al.  The variational multiscale method for laminar and turbulent flow , 2006 .

[33]  M. Eisenberg,et al.  Ionic Mass Transfer and Concentration Polarization at Rotating Electrodes , 1953 .

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

[35]  Wolfgang A. Wall,et al.  Towards a taxonomy for multiscale methods in computational mechanics: building blocks of existing methods , 2007 .

[36]  R. Codina,et al.  Time dependent subscales in the stabilized finite element approximation of incompressible flow problems , 2007 .

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

[38]  Wolfgang A. Wall,et al.  An algebraic variational multiscale-multigrid method for large-eddy simulation: generalized-α time integration, Fourier analysis and application to turbulent flow past a square-section cylinder , 2011 .

[39]  Wolfgang A. Wall,et al.  Variational Multiscale Methods for incompressible flows , 2007, Int. J. Comput. Sci. Math..

[40]  Xi-yun Lu,et al.  Large eddy simulation of turbulent channel flow with mass transfer at high-Schmidt numbers , 2003 .

[41]  Jong-Yeon Hwang,et al.  Direct Numerical Simulation of Turbulent Flow Around a Rotating Circular Cylinder , 2007 .

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

[43]  Wolfgang A. Wall,et al.  Time-dependent subgrid scales in residual-based large eddy simulation of turbulent channel flow , 2010 .

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

[45]  Effects of Schmidt Number on Turbulent Mass Transfer Around a Rotating Circular Cylinder , 2011 .