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.
A numerical procedure for geochemical compaction in the presence of discontinuous reactions
AGOSTI, ABRAMO;GIOVANARDI, BIANCA;FORMAGGIA, LUCA;SCOTTI, ANNA
2016-01-01
Abstract
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.File | Dimensione | Formato | |
---|---|---|---|
p1.pdf
Accesso riservato
:
Publisher’s version
Dimensione
2.44 MB
Formato
Adobe PDF
|
2.44 MB | Adobe PDF | Visualizza/Apri |
I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.