Mathematical model of electromigration allowing the deviation from electroneutrality

The structure of the double layer on the boundary between solid and liquid phases is described by various models, of which the Stern-Gouy-Chapman model is still commonly accepted. Generally, the solid phase is charged, which also causes the distribution of the electric charge in the adjacent diffuse layer in the liquid phase. We propose a new mathematical model of electromigration considering the high deviation from electroneutrality in the diffuse layer of the double layer when the liquid phase is composed of solution of weak multivalent electrolytes of any valence and of any complexity. The mathematical model joins together the Poisson equation, the continuity equation for electric charge, the mass continuity equations, and the modified G-function. The model is able to calculate the volume charge density, electric potential, and concentration profiles of all ionic forms of all electrolytes in the diffuse part of the double layer, which consequently enables to calculate conductivity, pH, and deviation from electroneutrality. The model can easily be implemented into the numerical simulation software such as Comsol. Its outcome is demonstrated by the numerical simulation of the double layer composed of a charged silica surface and an adjacent liquid solution composed of weak multivalent electrolytes. The validity of the model is not limited only to the diffuse part of the double layer but is valid for electromigration of electrolytes in general.