Simulation numérique d'un modèle de transport réactif

Lien vers la page Web de l'orateur Résumé : Le transport de contaminants dans le sous-sol risque de polluer durablement les nappes phréatiques. Pour modéliser et simuler les processus physico-chimiques, il faut coupler une équation de transport à des équations non linéaires de chimie. Une méthode dite globale repose sur une discrétisation spatio-temporelle du modèle couplé. Le schéma implicite en temps utilise un pas de temps adaptatif et une mise à jour adaptative du Jacobien dans les itérations de Newton. Les contraintes de positivité des concentrations peuvent être satisfaites avec des variables logarithmiques, au prix parfois d'un mauvais conditionnement du Jacobien. Dans le cas d'une réaction de précipitation-dissolution, l'apparition et la disparition d'une phase minérale sont modélisées par un problème de complémentarité non linéaire. Des résultats numériques avec le benchmark du groupe Momas illustrent les performances de la méthode.