A fourth‐order accurate, Numerov‐type, three‐point finite‐difference discretization of electrochemical reaction‐diffusion equations on nonuniform (exponentially expanding) spatial grids in one‐dimensional space geometry

The validity for finite‐difference electrochemical kinetic simulations, of the extension of the Numerov discretization designed by Chawla and Katti [J Comput Appl Math 1980, 6, 189–196] for the solution of two‐point boundary value problems in ordinary differential equations, is examined. The discretization is adapted to systems of time‐dependent reaction‐diffusion partial differential equations in one‐dimensional space geometry, on nonuniform space grids resulting from coordinate transformations. The equations must not involve first spatial derivatives of the unknowns. Relevant discrete formulae are outlined and tested in calculations on two example kinetic models. The models describe potential step chronoamperometry under limiting current conditions, for the catalytic EC, and Reinert‐Berg CE reaction mechanisms. Exponentially expanding space grid is used. The discretization considered proves the most accurate and efficient, out of all the three‐point finite‐difference discretizations on such grids, that have been used thus far in electrochemical kinetics. Therefore, it can be recommended as a method of choice. © 2004 Wiley Periodicals, Inc. J Comput Chem 25: 1515–1521, 2004

[1]  S. Miller.  Polarographic Currents from a Combination of Diffusion and Reaction , 1952 .

[2]  D. Britz,et al.  Higher-order spatial discretisations in digital simulations. Algorithm for any multi-point first- or second derivative on an arbitrarily spaced grid , 2003 .

[3]  Leslaw K. Bieniasz,et al.  Use of the Numerov Method to Improve the Accuracy of the Spatial Discretisation in Finite-difference Electrochemical Kinetic Simulations , 2002, Comput. Chem..

[4]  Manfred Rudolph.  Digital simulations on unequally spaced grids.: Part 1. Critical remarks on using the point method by discretisation on a transformed grid , 2002 .

[5]  Dieter Britz,et al.  Digital Simulation in Electrochemistry , 1981 .

[6]  B. Fornberg.  Generation of finite difference formulas on arbitrarily spaced grids , 1988 .

[7]  Linda R. Petzold,et al.  Algorithms and software for ordinary differential equations and differential-algebraic equations, part II: higher-order methods and software packages , 1995 .

[8]  E. Hairer,et al.  Solving Ordinary Differential Equations II: Stiff and Differential-Algebraic Problems , 2010 .

[9]  M. D. Rooij,et al.  Electrochemical Methods: Fundamentals and Applications , 2003 .

[10]  L. Collatz.  The numerical treatment of differential equations , 1961 .

[11]  W. Cody,et al.  Rational Chebyshev approximations for the error function , 1969 .

[12]  Ole Østerby,et al.  High-order Spatial Discretisations in Electrochemical Digital Simulation. Part 3. Combination with the Explicit Runge-Kutta Algorithm , 2002, Comput. Chem..

[13]  G. L. Stiehl,et al.  Theory of Catalytic Polarographic Currents , 1952 .

[14]  Finite difference methods for a class of two-point boundary value problems with mixed boundary conditions , 1980 .

[15]  R. Bellman,et al.  Quasilinearization and nonlinear boundary-value problems , 1966 .

[16]  M. A. Akanbi,et al.  Numerical solution of initial value problems in differential - algebraic equations , 2005 .

[17]  Dieter Britz,et al.  High-order Spatial Discretisations in Electrochemical Digital Simulation. 1. Combination with the BDF Algorithm , 2000, Comput. Chem..

[18]  Leslaw K. Bieniasz,et al.  High order accurate, one-sided finite-difference approximations to concentration gradients at the boundaries, for the simulation of electrochemical reaction-diffusion problems in one-dimensional space geometry , 2003, Comput. Biol. Chem..

[19]  Stephen W. Feldberg,et al.  Optimization of explicit finite-difference simulation of electrochemical phenomena utilizing an exponentially expanded space grid , 1981 .

[20]  Leslaw K. Bieniasz,et al.  Improving the accuracy of the spatial discretization in finite‐difference electrochemical kinetic simulations, by means of the extended Numerov method , 2004, J. Comput. Chem..

[21]  B. V. Noumerov.  A Method of Extrapolation of Perturbations , 1924 .

[22]  L. Bieniasz.  Comments on the paper by M. Rudolph, entitled “Digital simulations on unequally spaced grids. Part 1. Critical remarks on using the point method by discretisation on a transformed grid” , 2003 .

[23]  D. Britz,et al.  Higher-order spatial discretisations in electrochemical digital simulations. Part 4. Discretisation on an arbitrarily spaced grid , 2003, Comput. Biol. Chem..

[24]  L. Bieniasz.  Automatic derivation of the governing equations that describe a transient electrochemical experiment, given a reaction mechanism of arbitrary complexity. Part 2. Governing equations in one-dimensional geometry , 1996 .

[25]  Jörg Strutwolf,et al.  Digital simulation of potential step experiments using the extrapolation method , 1997 .