LBGK method coupled to time splitting technique for solving reaction-diffusion processes in complex systems.

A new approach to numerically solve a reaction-diffusion system is given, specifically developed for complex systems including many reacting/diffusing species with broad ranges of rate constants and diffusion coefficients, as well as complicated geometry of reacting interfaces. The approach combines a Lattice Boltzmann (LB) method with a splitting time technique. In the present work, the proposed approach is tested by focusing on the typical reaction process between a metal ion M and a ligand L, to form a complex ML with M being consumed at an electrode. The aim of the paper is to systematically study the convergence conditions of the associated numerical scheme. We find that the combination of LB with the time splitting method allows us to solve the problem for any value of association and dissociation rate constant of the reaction process. Also, the method can be extended to a mixture of ligands. We stress two main points: (1) the LB approach is particularly convenient for the flux computation of M and (2) the splitting time procedure is very well suited for reaction processes involving association-dissociation rate constants varying on many orders of magnitude.

[1]  A. Fisher,et al.  The boundary element method: applications to steady-state voltammetric simulations within domains extending to infinity , 2001 .

[2]  Shiyi Chen,et al.  Lattice Boltzmann computations for reaction‐diffusion equations , 1993 .

[3]  H. P. van Leeuwen,et al.  Voltammetry of metal complex systems with different diffusion coefficients of the species involved. I. Analytical approaches to the limiting current. , 1987 .

[4]  L. Luo,et al.  Theory of the lattice Boltzmann method: From the Boltzmann equation to the lattice Boltzmann equation , 1997 .

[5]  J. Shadid,et al.  Studies on the accuracy of time-integration methods for the radiation-diffusion equations , 2004 .

[6]  Qing Li,et al.  LBGK Simulations of Turing Patterns in CIMA Model , 2001, J. Sci. Comput..

[7]  David J. Gavaghan,et al.  Two-dimensional implementation of the finite element method with singularity correction for diffusion limited current at an unshielded disc electrode , 1995 .

[8]  Josep Galceran,et al.  Evaluation of the Koutecký-Koryta approximation for voltammetric currents generated by metal complex systems with various labilities , 2002 .

[9]  M. Esteban,et al.  Voltammetry of labile metal-macromolecular systems for any ligand-to-metal ratio, including adsorption phenomena. The role of the stability constant , 1994 .

[10]  Chopard,et al.  Reaction-diffusion cellular automata model for the formation of Leisegang patterns. , 1994, Physical review letters.

[11]  P. Bhatnagar,et al.  A Model for Collision Processes in Gases. I. Small Amplitude Processes in Charged and Neutral One-Component Systems , 1954 .

[12]  Josep Galceran,et al.  Lability and mobility effects on mixtures of ligands under steady-state conditions , 2003 .