Skip to Main content Skip to Navigation

Improvement of the numerical capacities of simulation tools for reactive transport modeling in porous media

Abstract : Reactive transport modeling in porous media involves the simulation of several physico‑chemical processes: flow of fluid phases, transport of species, heat transport, chemical reactions between species in the same phase or in different phases. The resolution of the system of equations that describes the problem can be obtained by a fully coupled approach or by a decoupled approach. Decoupled approaches can simplify the system of equations by breaking down the problem into smaller parts that are easier to handle. Each of the smaller parts can be solved with suitable integration techniques. The decoupling techniques might be non‑iterative (operator splitting methods) or iterative (fixed‑point iteration), having each its advantages and disadvantages. Non‑iterative approaches have an error associated with the separation of the coupled effects, and iterative approaches might have problems to converge. In this thesis, we develop an open‑source code written in MATLAB ( in order to model the problematic of concrete atmospheric carbonation for an intermediate‑level long‑lived nuclear waste package in a deep geological repository. The code uses a decoupled approach. Classical operator splitting approaches, such as sequential, alternating or Strang splitting, and less classical splitting approaches, such as additive or symmetrically weighted splitting, have been implemented. Besides, two iterative approaches based on an specific formulation (SIA CC, and SIA TC) have also been implemented. The code has been interfaced in a generic way with different transport solvers (COMSOL, pdepe MATLAB, FVTool, FD scripts) and geochemical solvers (iPhreeqc, PhreeqcRM). In order to validate the implementation of the different approaches, a series of classical benchmarks in the field of reactive transport have been solved successfully and compared with analytical and external numerical solutions. Since the associated error due to the combination of operator splitting and numerical techniques may be complex to assess, we explore the existing mathematical tools used to evaluate it. Finally, we frame the atmospheric carbonation problem and run preliminary simulations, stating the relevant problems and future steps to follow.
Document type :
Complete list of metadatas

Cited literature [218 references]  Display  Hide  Download
Contributor : Abes Star :  Contact
Submitted on : Tuesday, December 12, 2017 - 1:00:46 AM
Last modification on : Thursday, March 5, 2020 - 5:15:46 PM


Version validated by the jury (STAR)


  • HAL Id : tel-01661536, version 1


Daniel Jara Heredia. Improvement of the numerical capacities of simulation tools for reactive transport modeling in porous media. Earth Sciences. Université Rennes 1, 2017. English. ⟨NNT : 2017REN1S036⟩. ⟨tel-01661536⟩



Record views


Files downloads