Numerical method with convection boundary condition

R.J. Kee and A.A. McKillop. A Numerical Method for Predicting Natural Convection in Horizontal Cylinders with Asymmetric Boundary Conditions. Comput. Fluids, 5 1-14,1977. 

The modeling of the current distribution in a general-geometry cell nearly always requires a numerical solution. The following discussion focuses on the thin boundary layer approximation, with the overpotential components lumped within a thin boundary layer which may be of a varying thickness. The Laplace equation for the potential with nonlinear boundary conditions must be solved. Similar considerations typically apply to the more comprehensive solution of the Nernst-Planck equation (10) however, the need to account for the convective fluid flow in the latter case makes the application of the boundary methods more complex. We focus our brief discussion on the most common methods the finite-difference method, the finite-element method, and the boundary-element method, schematically depicted in Fig. 4. Since the finite-difference method is the simplest to implement and the best known technique, it is discussed in somewhat more detail. 

---