Formulation généralisée du transport réactif pour les modèles de réseaux de pores saturés en eau

A generalized solution for reactive transport in saturated porous networks

• Auteur / Autrice : Toko Kamtchueng
• Direction : Jean-Louis Rouet, Mohamed Azaroual
• Type : Thèse de doctorat
• Discipline(s) : Physique
• Date : Soutenance le 07/12/2016
• Etablissement(s) : Orléans
• Ecole(s) doctorale(s) : École doctorale Énergie, Matériaux, Sciences de la Terre et de l'Univers (Centre-Val de Loire ; 2012-....)
• Partenaire(s) de recherche : Laboratoire : Institut des sciences de la terre d'Orléans (2010-....)
• Jury : Président / Présidente : Lionel Mercury
• Jury : Examinateurs / Examinatrices : Jean-Louis Rouet, Mohamed Azaroual, Lionel Mercury, Benoît Nœtinger, Marc Prat, Claire Chainais-Hillairet
• Jury : Rapporteurs / Rapporteuses : Benoît Nœtinger, Marc Prat
• DOI : 10.70675/d2fc67c0zf17fz4905z9196zd7c2ef5cc227

Mots clés

Mots clés contrôlés

Écoulement (hydrologie)

Porosité

Polluants

Matériaux poreux

Ressources en eau

Mots clés libres

Transport réactif

Milieu poreux saturé

Modèle de réseau de pores

Schéma numérique

Reactive transport

Saturated porous media

Pore Network Model

Numerical scheme

Résumé

La protection et la remédiation des ressources en eau sont un enjeu sociétal majeur, ainsi il est nécessaire de comprendre l’évolution de solutés, tels les polluants, au sein de la zone saturée et non saturée. Dans ce but, de nombreux travaux ont été consacrés à la modélisation du transport réactif en milieux poreux. Son déroulement à l’échelle de Darcy dépend des hétérogénéités microscopiques du milieu. Les modèles de réseau de pores qui simplifient la géométrie en un ensemble de pores reliés par des liens de sections constantes, permettent de se placer à une échelle mésoscopique, faisant le lien entre l’échelle porale et l’échelle de Darcy. Sur de telles formes géométriques, l’écoulement admet un traitement analytique. En ce qui concerne le transport réactif des solutés, nous proposons une solution analytique dans les liens qui permet de calculer le débit de masse entre pores. Le modèle de transport se formule alors comme un système d’équations de Volterra de secondes espèces dont les noyaux de convolution sont des séries d’exponentielles décroissantes (hormis le premier terme qui est constant). Leurs temps de relaxation sont pilotés essentiellement par le temps de dispersion td. Dans la limite où td tend vers 0 à Péclet constant, les termes transitoires des noyaux se réduisent à un Dirac, débouchant sur un premier modèle simplifié à réponse instantanée c'est-à-dire un modèle de transport quasi-statique. Dans le cas où les volumes des pores sont suffisamment grands, les noyaux se réduisent à leur premier terme. Ces formulations du transport généralisent celles de la littérature. En particulier pour des Péclet petit ou grand on retrouve respectivement les modèles usuels en régime dispersif et convectif. Numériquement, la décroissance exponentielle des noyaux permet d’optimiser le calcul des convolutions avec une précision arbitrairement fixée, réduisant drastiquement le temps de résolution.

Protection and remediation of ground water resources are a major societal challenge. It implies to understand the evolution of solutes as pollutant in the saturated and non-saturated zones. For that purpose numerous studies have been conducted for modeling the reactive transport in a porous media. At Darcy scale, the behavior of solutes depends on microscopic heterogeneity for the media. The Pore Network Models (PNM) simplifies drastically its geometry and considers pores linked by straight throats the section of which is constant. They give a description which is in between the macroscopic and the pore descriptions. With such geometry it is possible to use a Poiseuille flow modeling the flux. With respect to the reactiontransport equation, we seek the analytical solution of the CDE in throats, which in turn allows computing the mass flux in pores. The transport solution consists of a Volterra equation system. Its convolution kernels result in a summation of time function which is decreasing exponentially with time (except the first term which still constant). The time constant is driven by the diffusion time td. As td goes to zero, keeping the Peclet number fixed, each term of the summation reduces to a Dirac. The response of the system is then instantaneous. When the volume of the pore is large enough it is possible to neglect all the term of the kernel except the constant one. In the limit where the Peclet number goes to zero, usual models are recovered. Numerically, the exponential time decreasing of the kernel allow to optimize their computational time up to an arbitrary fixed precision.