A numerical procedure for geochemical compaction in the presence of discontinuous reactions

The process by which rocks are formed from the burial of a fresh sediment involves the coupled effects of mechanical compaction and geochemical reactions. Both of them affect the porosity and permeability of the rock and, in particular, geochemical reactions can significantly alter them, since dissolution and precipitation processes may cause a structural transformation of the solid matrix. Often, the differential problems that arise from the modeling of these chemical reactions may present a discontinuous right hand side, where the discontinuity depends on the solution itself. In this work we have developed a numerical model to simulate this complex multi-physics problem by treating the discontinuous right hand side with a specially tailored event-driven numerical scheme. We show the performance of this strategy in terms of positivity and mass conservation, also in comparison with a more traditional approach that relies on a regularization of the discontinuous terms.