EN FR
EN FR


Section: New Software and Platforms

GRT3D

Global Reactive Transport in 3D

Keywords: Geochemistry - Dispersion - Scientific calculation - Simulation - Advection

Scientific Description: Participants : Édouard Canot, Jocelyne Erhel [correspondant] .

Version: version 2.0, April 2014

APP: registered

Programming language: C

Abstract: Reactive transport modeling has become an essential tool for understanding complex environmental problems. It is an important issue for MoMaS and C2S@EXA partners (see sections 8.2.5 , 8.2.3 ), in particular Andra. We have developed a method coupling transport and chemistry, based on a method of lines such that spatial discretization leads to a semi-discrete system of algebraic differential equations (DAE system). The main advantage is to use a complex DAE solver, which controls simultaneously the timestep and the convergence of Newton algorithm. The approach SIA uses a fixed-point method to solve the nonlinear system at each timestep, whereas the approach SNIA uses an explicit scheme.

The software suite GRT3D has four executable modules:

SIA1D: Sequential Iterative Approach for 1D domains,

GDAE1D: Global DAE approach for 1D domains,

SNIA3D: Sequential Non Iterative Approach for 1D, 2D or 3D domains.

GDAE3D: Global DAE approach for 1D, 2D or 3D domains. This module has three variants: the original one with logarithms, an optimized one still with logarithms, an optimized one which does not use logarithms.

Current work: extension of the chemistry module and parallelization.

Functional Description: Reactive transport modeling has become an essential tool for understanding complex environmental problems. It is an important issue for MoMaS and C2S@EXA partners, in particular Andra. We have developed a method coupling transport and chemistry, based on a method of lines such that spatial discretization leads to a semi-discrete system of algebraic differential equations (DAE system). The main advantage is to use a complex DAE solver, which controls simultaneously the timestep and the convergence of Newton algorithm. The approach SIA uses a fixed-point method to solve the nonlinear system at each timestep, whereas the approach SNIA uses an explicit scheme.

The software suite GRT3D has four executable modules:

SIA1D: Sequential Iterative Approach for 1D domains,

GDAE1D: Global DAE approach for 1D domains,

SNIA3D: Sequential Non Iterative Approach for 1D, 2D or 3D domains.

GDAE3D: Global DAE approach for 1D, 2D or 3D domains. This module has three variants: the original one with logarithms, an optimized one still with logarithms, an optimized one which does not use logarithms.

  • Participants: Caroline De Dieuleveult, Édouard Canot, Jocelyne Erhel, Nadir Soualem and Souhila Sabit

  • Partner: ANDRA

  • Contact: Jocelyne Erhel