Abstract
In this paper a mathematical description and numerical implementation for ion transport in concrete due to current passage is developed, in which the heterogeneous equilibrium between Ca^2+^, OH^β^ and the solid Ca(OH)~2~ is incorporated. The description is based on the NernstβPlanck equation for ion transport, and reaction terms for the dissolution/precipitation of Ca(OH)~2~. This description was implemented in the finite element package Comsol Multiphysics. In this way Ca(OH)~2~ depletion in a zone at a CP anode adjacent to a bulk of concrete with Ca(OH)~2~ could be modelled in one calculation. Drawback of this model is that the kinetic parameters in the reaction terms are not known, and must be chosen high to ensure the dissolution of Ca(OH)~2~ to be in equilibrium. This proved numerically challenging and sometimes caused long calculation times.
The growth rate of the zone without solid depends on the current density applied, concrete cover, the pore liquid composition and the diffusion constants of Ca^2+^ and OH^β^. This rate must be evaluated numerically. This qualitative model of anode acidification shows no participation of Na^+^; therefore transport properties of this ion do not affect the acidification rate of concrete. The same would hold for any other ion included in the model, which is not involved in electrochemical or chemical reactions.