Abstract
Concrete is a porous material in which several chemical reactions may develop. One way to handle such a complex set of reactions is to use a modified form of Petersen matrix notation which is described in the paper. This formulation of Petersen matrix is able to handle the case where some of the reactions involved in the set of reactions have no analytical expression, but are described by chemical equilibrium. In order to handle this new situation, the reaction rate is tuned during the simulations in such a way that the chemical equilibrium remains verified. This tuning is performed thanks to a numerical proportional integral (PI) controller technique. This extension of the Petersen matrix is applied to the case of concrete carbonation modelling. The PI coefficients are given for each equilibrium controlled kinetics. The simulation results obtained thanks to a finite-difference solver are compared to experimental results taken from the literature.