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

Résumé élisation du transport réactif avec des processus biologiques, chimiques ou radioactifs nécessite des méthodes numériques capables de reproduire des fronts de concentration très abrupts. Dans ce travail, nous présentons une nouvelle méthode ELLAM (Eulerian–Lagrangian Localized Adjoint Method) pour résoudre lʹéquation de transport réactif avec des coefficients non constants. Afin dʹéviter les interpolations (source dʹerreurs), la méthode utilise des maillages variables pour définir la solution ainsi que les fonctions test. Cette méthode est utilisée pour simuler la propagation dʹun front raide en milieu poreux non saturé ainsi que le transport multi-espèces. Les résultats ne présentent, ni oscillations, ni diffusion numérique, et restent précis, même en utilisant de grands pas de temps. Pour citer cet article : A. Younes, C. R. Geoscience 336 (2004). ing 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).