Discretization methods with analytical characteristic methods for advection-diffusion-reaction equations and 2d applications

Geiser, Jürgen

ESAIM: Mathematical Modelling and Numerical Analysis , Tome 43 (2009) no. 6, pp. 1157-1183.

Résumé

Our studies are motivated by a desire to model long-time simulations of possible scenarios for a waste disposal. Numerical methods are developed for solving the arising systems of convection-diffusion-dispersion-reaction equations, and the received results of several discretization methods are presented. We concentrate on linear reaction systems, which can be solved analytically. In the numerical methods, we use large time-steps to achieve long simulation times of about $10000$ years. We propose higher-order discretization methods, which allow us to use large time-steps without losing accuracy. By decoupling of a multi-physical and multi-dimensional equation, simpler physical and one-dimensional equations are obtained and can be discretized with higher-order methods. The results can then be coupled with an operator-splitting method. We discuss benchmark problems given in the literature and real-life applications. We simulate a radioactive waste disposals with underlying flowing groundwater. The transport and reaction simulations for the decay chains are presented in 2d realistic domains, and we discuss the received results. Finally, we present our conclusions and ideas for further works.

Détail

DOI : 10.1051/m2an/2009033

Classification : 35K15, 35K57, 47F05, 65M60, 65N30
Mots-clés : advection-diffusion-reaction equation, embedded analytical solutions, operator-splitting methods, characteristic methods, finite-volume methods, multi-physics, simulation of radioactive waste disposals