The Nernst-Planck-Poisson Reactive Transport Model for Concrete Carbonation and Chloride Diffusion in Carbonated and Non-carbonated Concrete

Abstract

The intrusion of chlorides and carbon dioxide into a reinforced concrete (RC) structure can initiate corrosion of the reinforcing steel, which, due to its expansive nature, can damage the structure and adversely affects its serviceability and safety. Corrosion will initiate if at the steel surface the concrete free chloride concentration exceeds a defined limit, or its pH falls below a critical level. Hence, determination of the time to reaching these critical limits is key to the assessment of RC structures durability and service life. Due to the ionic nature of the chlorides and the bicarbonate anion (HCO3-) formed by the CO2 in the multi-ionic pore solution, the transport of both species is driven by Fickian diffusion combined with electromigration and ionic activity, which can be mathematically expressed by the Nernst-Planck-Poisson (NPP) equations. For a complete representation of the phenomenon, however, the NPP equations must be supplemented by the relevant chemical equilibrium equations to ensure chemical balance among the various species within the concrete pore solution. The combination of NPP with the chemical equilibrium equations is often termed the NPP reactive transport model. In this study, such a model is developed, coded into the MATLAB platform, validated by available experimental data, and applied to analyze the time-dependent concrete carbonation and the movement of chlorides in carbonated and non-carbonated concrete. The results of these analyses can be used to predict the time to corrosion initiation. The transient one-dimensional governing equations of NPP are numerically solved using the Galerkin’s finite element formulation in space and the backward (implicit) Euler scheme in the time domain. The associated system of chemical equilibrium equations accounts for the key homogeneous and heterogeneous chemical reactions that take place in the concrete during carbonation and chlorides transport. At each stage of the analysis, the effects of these reactions on the changes in the pore solution chemical composition, pH, cement chloride binding capacity, concrete porosity, and the hydrated cement solids volumetric ratio are determined. The study demonstrates that given accurate input data, the presently developed NPP reactive transport model can accurately simulate the complex transport processes of chlorides and CO2 in concrete as a reactive porous medium, and the ensuing physical and chemical changes that occur due to the reaction of these species with the pore solution and the other cement hydration products. This conclusion is supported by the good agreement between results of the current analyses with the corresponding available experimental data from physical tests involving carbonation, and chloride diffusion in non-carbonated and carbonated concrete.