Modélisation du transport réactif en milieux poreux non saturés avec une méthode ELLAM en maillage variable

Modelling contaminant transfer with biological/chemical/radioactive processes needs appropriate numerical methods able to reproduce sharp concentration fronts. In this work, we develop a new Eulerian–Lagrangian Localized Adjoint Method (ELLAM) for solving the reactive transport equation with non-constant coefficients. To avoid interpolation (leading to errors), we use a moving grid to define the solution and test functions. The method is used to simulate first the infiltration of solute into a column of unsaturated porous medium and second the multispecies transport. The developed ELLAM gives accurate results without non-physical oscillations or numerical diffusion, even when using large time steps. To cite this article: A. Younes, C. R. Geoscience 336 (2004).