  We analyze the simplest problem of electrochemical relaxation in more than
one dimension - the response of an uncharged, ideally polarizable metallic
sphere (or cylinder) in a symmetric, binary electrolyte to a uniform electric
field. In order to go beyond the circuit approximation for thin double layers,
our analysis is based on the Poisson-Nernst-Planck (PNP) equations of dilute
solution theory. Unlike most previous studies, however, we focus on the
nonlinear regime, where the applied voltage across the conductor is larger than
the thermal voltage. In such strong electric fields, the classical model
predicts that the double layer adsorbs enough ions to produce bulk
concentration gradients and surface conduction. Our analysis begins with a
general derivation of surface conservation laws in the thin double-layer limit,
which provide effective boundary conditions on the quasi-neutral bulk. We solve
the resulting nonlinear partial differential equations numerically for strong
fields and also perform a time-dependent asymptotic analysis for weaker fields,
where bulk diffusion and surface conduction arise as first-order corrections.
We also derive various dimensionless parameters comparing surface to bulk
transport processes, which generalize the Bikerman-Dukhin number. Our results
have basic relevance for double-layer charging dynamics and nonlinear
electrokinetics in the ubiquitous PNP approximation.
