The nonuniformity of the potential and current distributions at electrodes was recognized long ago by Wagner 1 and others. Studies have focused on the solution of the diffusion-migration equation based on the transport law. One comprehensive description of the equation was given by Newman. 2 Mathematical modeling progress before the early 1980s was summarized in several review papers, e.g. Ref. 3 by Ibl. Newman's initial work 4 and further extensions 5 indicated that the current density at a disk electrode increased with increasing distance from the electrode center, reaching its maximum value at the disk edge, infinite in the case of the primary distribution and finite in the case of the secondary or tertiary distributions. Later work by Shih and Pickering, 6 following the approach by Wagner, 1 employed the Green theorem and a special polarization parameter to compute the secondary distribution of a square electrode. Their threedimensional results were consistent with Newman's and showed further that the highest currents occurred at the corners. Cahan et al. 7 developed an iterative boundary integral equation method to calculate the potential and current distributions. The method was capable of handling problems with complex boundary geometry and nonlinear boundary conditions and hence could treat the secondary distribution problems. Dimpault-Darcy and White, 8 on the other hand, used a commercial finite element method (FEM) software to solve problems with Butler-Volmer-type boundary conditions. Their method may handle the noncontinuity within the system and convection-type conditions at the boundary. The method was further developed and used for the design of an electroplating cell. 9 A review of the distribution literature in the electrodeposition discipline was conducted by Dukovic. 10 In addition to the numerous theoretical/mathematical contributions, quite a few experimental works on the nonuniform potential and current distributions were reported. Papers of Remppel and Exner, 11 Dinan et al., 12 and Ponthiaux et al. 13,14 addressed this issue from different angles. Recently, Kranc and Sagüés 15 published their work of modeling the time-dependent response to external polarization of a corrosion macrocell on steel in concrete.When modeling a corrosion system, the primary distributions are easily obtained by routine methods of numerical analysis. However, the results are usually far from reality due to the electrode polarization that occurs immediately after the circuit is switched on. The polarization changes the initial boundary conditions used for the computation. More accurate secondary distribution results may be obtained after the polarization effect is taken into consideration. A trial and error technique plus iteration can be used for this purpose to get accurate results. One description of the method was given by Ibl. 3 The iteration method used in the past was usually for systems in which the electrodes were only in the active state. For a partially passivated electrode, the computation is more di...