Reactive transport modelling of iron bentonite interaction in the FEBEX in situ experiment

The evolution of steel/bentonite interfaces in engineered barrier systems of radioactive waste disposal is governed by the corrosion of the steel and the interaction of corrosion products with the bentonite. In the early post-closure phase of the repository, the transition from aerobic to anaerobic corrosion in combination with evolving temperature conditions and chemical gradients due to the hydration of the bentonite take place. In situ tests in underground laboratories provide a high degree of representativeness for the early transient evolution of the interfaces. Complex Fe-bentonite interaction patterns have been observed, which include a concave Fe accumulation front and distinct coloured halos around the corroding steel. A reactive transport model for the FEBEX in situ test is presented, which describes the transition from aerobic to anaerobic corrosion, the transport of O2 in the gas and liquid phase and the chemical evolution of the steel/bentonite interface during the 18 years of the experiment. The model successfully reproduces the major steps of the interface evolution: an aerobic corrosion phase characterized by accumulation of goethite in the corrosion layer, and an anaerobic corrosion phase, where Fe(II) and Fe(II)/Fe(III) corrosion products form, including a gradual re-dissolution of the aerobic corrosion product. In the anaerobic corrosion phase, some of the Fe(II) diffuses into the bentonite, where it partly reacts with remaining O2 to form Fe(III) precipitates or interacts with the clay minerals via surface complexation or cation exchange. Two different approaches for the implementation of a potential electron transfer from sorbed Fe(II) to structural Fe(III) are presented, but the lack of experimental data does not allow evaluation of their representativeness for the FEBEX in situ experiment. The model calculations qualitatively reproduce the concave shape of the Fe accumulation front in the bentonite and a distinct zonation with an iron-oxide precipitation dominated zone adjacent to the steel and a Fe(II) sorption dominated zone further into the bentonite. Sensitivity cases with respect to the hydration of the bentonite and the transport parameters point to the importance of the fast O2 diffusion in the gas phase for the formation of these characteristic interface features. The calculated thickness of the Fe-bentonite interaction zone varies from 1 to 4 cm, which is smaller than observed in some areas of the experiment, where locally this zone extended > 10 cm into the bentonite. It is possible that this difference is due to a more complex flow pattern in the in situ experiment that is not captured by the model or due to an additional decoupled electron transport across oxide rich layers.