Numerical and analytical studies of the electrical conductivity of a concentrated colloidal suspension

J Phys Chem B. 2006 Mar 30;110(12):6179-89. doi: 10.1021/jp057030e.

Abstract

In the past few years, different models and analytical approximations have been developed facing the problem of the electrical conductivity of a concentrated colloidal suspension, according to the cell-model concept. Most of them make use of the Kuwabara cell model to account for hydrodynamic particle-particle interactions, but they differ in the choice of electrostatic boundary conditions at the outer surface of the cell. Most analytical and numerical studies have been developed using two different sets of boundary conditions of the Neumann or Dirichlet type for the electrical potential, ionic concentrations or electrochemical potentials at that outer surface. In this contribution, we study and compare numerical conductivity predictions with results obtained using different analytical formulas valid for arbitrary zeta potentials and thin double layers for each of the two common sets of boundary conditions referred to above. The conductivity will be analyzed as a function of particle volume fraction, phi, zeta potential, zeta, and electrokinetic radius, kappaa (kappa(-1) is the double layer thickness, and a is the radius of the particle). A comparison with some experimental conductivity results in the literature is also given. We demonstrate in this work that the two analytical conductivity formulas, which are mainly based on Neumann- and Dirichlet-type boundary conditions for the electrochemical potential, predict values of the conductivity very close to their corresponding numerical results for the same boundary conditions, whatever the suspension or solution parameters, under the assumption of thin double layers where these approximations are valid. Furthermore, both analytical conductivity equations fulfill the Maxwell limit for uncharged nonconductive spheres, which coincides with the limit kappaa --> infinity. However, some experimental data will show that the Neumann, either numerical or analytical, approach is unable to make predictions in agreement with experiments, unlike the Dirichlet approach which correctly predicts the experimental conductivity results. In consequence, a deeper study has been performed with numerical and analytical predictions based on Dirichlet-type boundary conditions.