Finite-difference electrochemical kinetic simulations using the Rosenbrock time integration scheme

Abstract Temporal integration based on multistep backward differentiation formulae may not be the best choice for finite-difference electrochemical kinetic simulations. A one-step, highly accurate, unconditionally stable and non-oscillatory implicit time scheme that also provides temporal error estimates, would be more desirable. The Lawson–Morris–Gourlay extrapolation (LMGE) method, suggested by Strutwolf and Schoeller, partially satisfies these requirements, but it is computationally rather expensive, especially for nonlinear governing equations. The Rosenbrock method, thus far unused in electrochemical simulations, is shown here to be competitive. By combining a three-stage ROWDA3 Rosenbrock method for differential–algebraic equations (DAEs), with a finite-difference spatial discretization of the governing equations, an effective simulation algorithm can be obtained. This is demonstrated for three example kinetic models represented by a single reaction–diffusion equation in one space dimension: potential step experiment, linear sweep voltammetry for a first- and second-order irreversible heterogeneous charge transfer, and linear sweep voltammetry for an enzymatic catalytic reaction mechanism. Numerical properties of the algorithm are compared in calculations with those of the finite-difference methods based on first- and second-order backward differentiation formulae, the Crank–Nicolson (CN) method, and the second-order LMGE method. Fixed, uniform grids are used. Under conditions of a large temporal-to-spatial grid step ratio, important for electrochemical simulations, the third-order accurate ROWDA3 integration proves more efficient than the second-order extrapolation method. In one case it is also more efficient than the method based on the second-order (two-step) backward differentiation formula.

[1]  Richard G. Compton,et al.  Modelling electrode transients: the strongly implicit procedure , 1996 .

[2]  Pentti Laasonen,et al.  Über eine Methode zur Lösung der Wärmeleitungs-gleichung , 1949 .

[3]  B. A. Schmitt,et al.  Polynomial preconditioning in Krylov-ROW-methods , 1998 .

[4]  M. Roche,et al.  Rosenbrock methods for Differential Algebraic Equations , 1987 .

[5]  R. D. Richtmyer,et al.  Difference methods for initial-value problems , 1959 .

[6]  Stephen W. Feldberg,et al.  Examination of the behavior of the fully implicit finite-difference algorithm with the Richtmyer modification: behavior with an exponentially expanding time grid , 1995 .

[7]  D. Britz.  Electrochemical digital simulation by runge-kutta integration , 1988 .

[8]  Richard G. Compton,et al.  A General Method for Electrochemical Simulations. 1. Formulation of the Strategy for Two-Dimensional Simulations , 1997 .

[9]  L. Bieniasz.  Use of dynamically adaptive grid techniques for the solution of electrochemical kinetic equations: Part 4: The adaptive moving-grid solution of one-dimensional fast homogeneous reaction-diffusion problems with extremely thin reaction zones away from the electrodes , 1994 .

[10]  Linda R. Petzold,et al.  Numerical solution of initial-value problems in differential-algebraic equations , 1996, Classics in applied mathematics.

[11]  Christopher M.A. Brett,et al.  Electrochemistry: Principles, Methods, and Applications , 1993 .

[12]  G. Fairweather,et al.  A Survey of Higher-order Methods for the Numerical Integration of Semidiscrete Parabolic Problems , 1984 .

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

[14]  Huang Weiguang,et al.  Program Polar 1: Simulator of Voltammograms , 1995, Comput. Chem..

[15]  Dieter Britz.  Time Shift Artifacts and Start-up Protocols with the BDF Method in Electrochemical Digital Simulation , 1998, Comput. Chem..

[16]  Lesław K. Bieniasz.  Use of dynamically adaptive grid techniques for the solution of electrochemical kinetic equations: Part 3. An adaptive moving grid—adaptive time step strategy for problems with discontinuous boundary conditions at the electrodes , 1994 .

[17]  Helmut Podhaisky,et al.  Numerical experiments with Krylov integrators , 1998 .

[18]  P. Eiseman,et al.  Adaptive grid generation , 1987 .

[19]  Richard G. Compton,et al.  Diffusional mass transport to microband electrodes of practical geometries: A simulation study using the strongly implicit procedure , 1995 .

[20]  Stephen W. Feldberg,et al.  The Richtmyer modification of the fully implicit finite difference algorithm for simulations of electrochemical problems , 1994 .

[21]  Ole Østerby,et al.  Numerical Stability of Finite Difference Algorithms for Electrochemical Kinetic Simulations. Matrix Stability Analysis of the Classic Explicit, Fully Implicit and Crank-Nicolson Methods, Extended to the 3- and 4-point Gradient Approximation at the Electrodes , 1995, Comput. Chem..

[22]  Dieter Britz,et al.  Stability of the Backward Differentiation Formula (FIRM) Applied to Electrochemical Digital Simulation , 1997, Comput. Chem..

[23]  P. Delahay,et al.  Theory of Irreversible Waves in Oscillographic Polarography , 1953 .

[24]  Lesław K. Bieniasz,et al.  Use of dynamically adaptive grid techniques for the solution of electrochemical kinetic equations: Part 1. Introductory exploration of the finite-difference adaptive moving grid solution of the one-dimensional fast homogeneous reaction-diffusion problem with a reaction layer , 1993 .

[25]  J. Hale.  Transients in convective systems: I. Theory of galvanostatic, and galvanostatic with current reversal transients, at a rotating disc electrode , 1963 .

[26]  J. Heinze,et al.  Implicit calculation of boundary values in digital simulation applied to several types of electrochemical experiment , 1988 .

[27]  Leon Lapidus,et al.  Solution of Stiff Differential Equations and the Use of Imbedding Techniques , 1978 .

[28]  Alexander Ostermann,et al.  Rosenbrock methods for partial differential equations and fractional orders of convergence , 1993 .

[29]  Bernd Speiser.  Easiest — ein Programmsystem zur Simulation von und Parameterbestimmung aus Elektroanalytischen Experimenten , 1989 .

[30]  Roland W. Lewis,et al.  A comparison of time marching schemes for the transient heat conduction equation , 1975 .

[31]  J. Lang.  Two-dimensional fully adaptive solutions of reaction-diffusion equations , 1995 .

[32]  Joe F. Thompson.  A survey of dynamically-adaptive grids in the numerical solution of partial differential equations , 1984 .

[33]  A. R. Gourlay,et al.  The Extrapolation of First Order Methods for Parabolic Partial Differential Equations, II , 1978 .

[34]  J. Heinze,et al.  The Crank-Nicolson technique—an efficient algorithm for the simulation of electrode processes at macro- and microelectrodes , 1993 .

[35]  Tien D. Bui,et al.  The solution of a combustion problem with Rosenbrock methods , 1986, TOMS.

[36]  P. Schmidt,et al.  Global spline collocation in the simulation of electrochemical diffusion equations , 1980 .

[37]  B. A. Schmitt,et al.  ROWMAP—a ROW-code with Krylov techniques for large stiff ODEs , 1997 .

[38]  David K. Gosser,et al.  Cyclic Voltammetry: Simulation and Analysis of Reaction Mechanisms , 1993 .

[39]  Leslaw K. Bieniasz,et al.  ELSIM - A User-friendly PC Program for Electrochemical Kinetic Simulations. Version 1.0 - Solution of Integral Equations for Linear Scan and Cyclic Voltammetry , 1992, Comput. Chem..

[40]  J. Crank,et al.  A practical method for numerical evaluation of solutions of partial differential equations of the heat-conduction type , 1947 .

[41]  Manfred Rudolph,et al.  A fast implicit finite difference algorithm for the digital simulation of electrochemical processes , 1991 .

[42]  D. Britz.  Investigation of the relative merits of some n-point current approximations in digital simulation , 1987 .

[43]  N. N. Kalitkin,et al.  L q -decreasing monotonic schemes with complex coefficients and applications to complicated PDE systems , 1994 .

[44]  N. Winograd.  An implicit finite difference method simulation of spectro-electrochemical working curves , 1973 .

[45]  R. Compton,et al.  A general computational approach to linear sweep voltammetry at channel electrodes , 1992 .

[46]  Lesław K. Bieniasz,et al.  Use of dynamically adaptive grid techniques for the solution of electrochemical kinetic equations: Part 2. An improved finite-difference adaptive moving grid technique for one-dimensional fast homogeneous reaction-diffusion problems with reaction layers at the electrodes , 1994 .

[47]  Leslaw K. Bieniasz,et al.  ELSIM - a Problem-solving Environment for Electrochemical Kinetic Simulations. Version 3.0-solution of Governing Equations Associated with Interfacial Species, Independent of Spatial Coordinates Or in One-dimensional Space Geometry , 1997, Comput. Chem..

[48]  Ole Østerby,et al.  Some numerical investigations of the stability of electrochemical digital simulation, particularly as affected by first-order homogeneous reactions , 1994 .

[49]  Jürgen Heinze,et al.  Digital simulation of cyclic voltammetric curves by the implicit Crank-Nicolson technique , 1984 .

[50]  Z. Galus.  Fundamentals of electrochemical analysis , 1976 .

[51]  L. Bieniasz.  A method-oriented approach to the formulation of algorithms for electrochemical kinetic simulations. Part 2. Extension to kinetic problems characterized by the simultaneous presence of bulk and interfacial species , 1996 .

[52]  Ole Østerby,et al.  Numerical Stability of Finite Difference Algorithms for Electrochemical Kinetic Simulations: Matrix Stability Analysis of the Classic Explicit, Fully Implicit and Crank-Nicolson Methods and Typical Problems Involving Mixed Boundary Conditions , 1995, Comput. Chem..

[53]  L. Bieniasz.  Automatic derivation of the governing equations that describe a transient electrochemical experiment, given a reaction mechanism of arbitrary complexity. Part 1. Problem parameters and initial conditions , 1996 .

[54]  George W. Batten,et al.  Second-Order Correct Boundary Conditions for the Numerical Solution of the Mixed Boundary Problem for Parabolic Equations , 1963 .

[55]  Dieter Britz,et al.  Brute force digital simulation , 1996 .

[56]  H. H. Rosenbrock,et al.  Some general implicit processes for the numerical solution of differential equations , 1963, Comput. J..

[57]  川口 光年,et al.  R.D.Richtmyer: Difference Methods for Initial-Value Problems. Interscience Pub. Inc. New York, 1957, xii+238頁, 15×23cm, \2,600. , 1958 .

[58]  Bernd Speiser.  EASIEST - a program system for lectroanalytical imulation and parameter imation - I. Simulation of cyclic voltammetric and chronoamperometric experiments , 1990, Comput. Chem..

[59]  Leslaw K. Bieniasz,et al.  A Reaction Compiler for Electrochemical Kinetics , 1996, Comput. Chem..

[60]  Adrian Sandu,et al.  Benchmarking stiff ode solvers for atmospheric chemistry problems II: Rosenbrock solvers , 1997 .

[61]  Jens Lang,et al.  Adaptive FEM for reaction-diffusion equations , 1998 .

[62]  Manfred Rudolph,et al.  Digital simulations with the fast implicit finite difference (FIFD) algorithm: Part II. An improved treatment of electrochemical mechanisms with second-order reactions , 1992 .

[63]  Carlos M. Villa,et al.  Simulation of Complex Electrochemical Reaction Systems , 1995 .

[64]  Klaus Johannsen,et al.  Matrix Stability of the Backward Differentiation Formula in Electrochemical Digital Simulation , 1999, Comput. Chem..

[65]  G. Smith,et al.  Numerical Solution of Partial Differential Equations: Finite Difference Methods , 1978 .

[66]  W. Heineman,et al.  Laboratory techniques in electroanalytical chemistry , 1984 .

[67]  Irma Lavagnini,et al.  Digital simulation of cyclic voltammetric experiments for non-first-order electrode reactions , 1988 .

[68]  Richard G. Compton,et al.  A General Method for Electrochemical Simulations. 2. Application to the Simulation of Steady-State Currents at Microdisk Electrodes: Homogeneous and Heterogeneous Kinetics , 1997 .

[69]  Richard G. Compton,et al.  A comparison of finite difference algorithms for the simulation of microband electrode problems with and without convective flow , 1996 .

[70]  D. Britz.  Electrochemical digital simulation: incorporation of the Crank—Nicolson scheme and n-point boundary expression into the Rudolph algorithm , 1993 .

[71]  A. Lasia.  Applications of implicit finite difference technique to cyclic voltammetry , 1985 .

[72]  Ole Østerby,et al.  The Effect of the Discretization of the Mixed Boundary Conditions on the Numerical Stability of the Crank-Nicolson Algorithm of Electrochemical Kinetic Simulations , 1997, Comput. Chem..

[73]  Andrzej Lasia,et al.  General Model of Electrochemical Hydrogen Absorption into Metals , 1995 .

[74]  Richard G. Compton,et al.  Theory of ECE processes at hemicylinder and band microelectrodes , 1995 .

[75]  Leslaw K. Bieniasz,et al.  ELSIM - A PC Program for Electrochemical Kinetic Simulations. Version 2.0 - solution of the Sets of Kinetic Partial Differential Equations in One-dimensional Geometry, Using Finite Difference, Orthogonal Collocation Methods , 1993, Comput. Chem..

[76]  F. Anson.  Innovations in the Study of Adsorbed Reactants by Chronocoulometry. , 1966 .

[77]  L. Bieniasz.  The von Neumann stability of finite-difference algorithms for the electrochemical kinetic simulation of diffusion coupled with homogeneous reactions , 1993 .

[78]  Allen J. Bard,et al.  Electrochemical Methods: Fundamentals and Applications , 1980 .

[79]  C. Goldstein,et al.  A stability criterion for accurate simulation of electrochemical diffusion-kinetic phenomena at the rotating disk electrode and implications for simulation of diffusion-migration and other problems , 1996 .

[80]  J. Lang.  High-resolution self-adaptive computations on chemical reaction-diffusion problems with internal boundaries , 1996 .

[81]  Brian E. Conway,et al.  Modern Aspects of Electrochemistry , 1974 .

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

[83]  J. Heinze,et al.  Electrochemistry at Ultramicroelectrodes—Simulation of Heterogeneous and Homogeneous Kinetics by an Improved ADI—Technique , 1986 .