Thermokinetic Model of Borosilicate Glass Dissolution: Contextual Affinity

Abstract

ABSTRACTShort and long-term geochemical interactions of R7T7 nuclear glass with water at 100°C were simulated with the DISSOL thermokinetic computer code. Both the dissolved glass quantity and the resulting water composition, saturation states and mineral quantities produced were calculated as a function of time. The rate equation used in the simulation was first proposed by Aagaard and Hegelson: v = k+.S.a(H+)-n(l - e-(A/RT)). It simulates a gradually diminishing dissolution rate as the reaction affinity diminishes. The best agreement with 1-year experimental data was obtained with a reaction affinity calculated from silica activity (Grambow's hypothesis) rather than taking into account the activity of all the glass components as proposed by Jantzen and Plodinec. The concept of residual affinity was introduced by Grambow to express the fact that the glass dissolution rate does not cease. We prefer to replace the term “residual affinity” by “contextual affinity”, which expresses the influence on the dissolution rate of three factors: the solution chemistry, the metastability of SiO2(m), and the possible precipitation of certain aluminosilicates such as zeolites.