TY - GEN

T1 - Numerical simulation of geochemical compaction with discontinuous reactions

AU - Agosti, Abramo

AU - Formaggia, Luca

AU - Giovanardi, Bianca

AU - Scotti, Anna

PY - 2015

Y1 - 2015

N2 - The present work deals with the numerical simulation of porous media subject to the coupled effects of mechanical compaction and reactive flows that can significantly alter the porosity due to dissolution, precipitation or transformation of the solid matrix. These chemical processes can be effectively modelled as ODEs with discontinuous right hand side, where the discontinuity depends on time and on the solution itself. Filippov theory can be applied to prove existence and to determine the solution behaviour at the discontinuities. From the numerical point of view, tailored numerical schemes are needed to guarantee positivity, mass conservation and accuracy. In particular, we rely on an event-driven approach such that, if the trajectory crosses a discontinuity, the transition point is localized exactly and integration is restarted accordingly.

AB - The present work deals with the numerical simulation of porous media subject to the coupled effects of mechanical compaction and reactive flows that can significantly alter the porosity due to dissolution, precipitation or transformation of the solid matrix. These chemical processes can be effectively modelled as ODEs with discontinuous right hand side, where the discontinuity depends on time and on the solution itself. Filippov theory can be applied to prove existence and to determine the solution behaviour at the discontinuities. From the numerical point of view, tailored numerical schemes are needed to guarantee positivity, mass conservation and accuracy. In particular, we rely on an event-driven approach such that, if the trajectory crosses a discontinuity, the transition point is localized exactly and integration is restarted accordingly.

KW - Compaction processes

KW - Discontinuous RHS ODE

KW - Mineral dissolution/precipitation

UR - http://www.scopus.com/inward/record.url?scp=84938704025&partnerID=8YFLogxK

M3 - Conference contribution

AN - SCOPUS:84938704025

T3 - COUPLED PROBLEMS 2015 - Proceedings of the 6th International Conference on Coupled Problems in Science and Engineering

SP - 300

EP - 311

BT - COUPLED PROBLEMS 2015 - Proceedings of the 6th International Conference on Coupled Problems in Science and Engineering

A2 - Onate, Eugenio

A2 - Papadrakakis, Manolis

A2 - Schrefler, Bernhard A.

PB - International Center for Numerical Methods in Engineering

T2 - 6th International Conference on Computational Methods for Coupled Problems in Science and Engineering, COUPLED PROBLEMS 2015

Y2 - 18 May 2015 through 20 May 2015

ER -