Modeling of cementitious materials exposed to isothermal calcium leaching, considering process kinetics and advective water flow. Part 2: Numerical solution

The second part of the paper presents numerical solutions of the mathematical model of hydro-chemo-mechanical behavior of cementitious materials exposed to contact with deionized water of part 1. The model defines kinetics of the calcium leaching process instead of a direct application of a curve describing equilibrium between solid calcium in the material skeleton and the calcium dissolved in the pore solution. It further takes into account the advective flux of calcium ions. Both aspects are new as compared to previous models. The weak form of the governing equations of the model is derived first using the Galerkin method. Then, the equations are discretized in space with finite elements and in time domain with finite differences, and finally the procedures used for numerical solution of their discretized form are presented. Three numerical examples are solved to test the numerical solution procedure proposed and demonstrate its robustness for solution of 1D and 2D problems concerning fast and slow leaching of cement-based materials. The effect of various factors on the results concerning chemical degradation of structures made of cementitious materials is analyzed as well.