A numerical solution of the Nernst-Planck and Poisson equations is presented. The equations are discretized in a finite difference scheme using the method of lies on variable spatial and temporal grids. A Gear’s stiffly stable predictmector integration procedure which automatically adjusts the order of the predictor and corrector equations (and the step size) to ensure the accuracy of the results is incorporated. Some advantages of our approach over more classical ones are discussed. The numerical solution is applied to the study of the equilibrium and nonequilibrium diffuse electrical double layer (EDL) at the metal electrode/electrolyte solution interface. The prescription of this layer is similar to that used in the classical Gouy-Chapman theory. Electrode kinetics are described by the Butler-Volmer equation. Concentration, faradaic and displacement electric current densities, and electric potential profdea as functions of time a- the cell thickness, and partiahly in the EDL regions at the metal electrode/solution interfaces, are obtained. Two physical problems are studied: (i) the formation of the equilibrium EDL, and (ii) the transient response of the system to an electrical perturbation. Thew examplea illustrate the potential applications of the numerical method. Numerical techniques for the solution of transport problems in electrochemistry have contributed significantly to the analysis of many complex processes which are difficult to deal with using more conventional approaches. Though other techniques like the boundary element method are experiencing an increasing popularity,’ fmite differences2J has been one of the most widely used in electrochemistry since the pioneering work by Feldberg!.s The need for numerical solutions appears in a multitude of problems of practical interest, e.g., transport phenomena described by diffusional type or second-order parabolic PDEs (partial differential equations) subjected to temporal boundary conditions,*J diffusion-migration situations involving the coupled, nonlinear Nemst-Planck and Poisson equations,b10 convective diffusion problems in electrochemical cell~,~J etc. Modern numerical solutions to the Nemst-Planck and Poisson equations in electrochemical cells including electrical double layer (EDL) effects seem to be lacking.’ We present, here, a numerical solution of these equations for a system consisting of three ions
Read full abstract