MODELLING OF ELUTION TESTS WITH THE ICP INTERFACE BETWEEN COMSOL AND PHREEQC

RÉSUMÉ: L’interface iCP permet de coupler Phreeqc, un code de modélisation géochimique, et COMSOL, un solveur d’éléments finis, pour la modélisation du transport réactif dans les sols. Avec l’interface iCP, les équations décrivant le transport des espèces chimiques dans l’eau interstitielle sont tout d’abord résolues avec COMSOL. Les analyses chimiques pour chaque cellule sont ensuite transmises vers Phreeqc pour des calculs d’équilibre thermodynamique et de cinétique. Dans cet article, l’interface iCP a été testée en modélisant une série d’essais en laboratoire réalisés avec un sol naturel (sable calcaire). Ces essais incluent un essai en bécher durant lequel une solution de nitrates de plomb, cuivre, cadmium et zinc a été mise en contact avec le sol, un essai de traçage au bromure et un essai d’élution avec les mêmes nitrates de métaux qui ont été utilisés pour l’essai en bécher. Pour l’essai en bécher et l’essai d’élution, Phreeqc a été utilisé pour modéliser la dissolution de la calcite et la précipitation des carbonates de métaux. COMSOL a été utilisé pour modéliser le transport des espèces dissoutes pour l’essai de traceur et l’essai en bécher. Cet article est centré sur un couplage entre des phénomènes hydrodynamiques et géochimiques, mais l’interface iCP pourrait aussi être utilisée pour des modèles couplés impliquant des équations mécaniques (ex. cimentation) et thermiques (ex. chaleur produite par des réactions chimiques).

ABSTRACT: The iCP interface allows the coupling of Phreeqc, a geochemical modelling code, and COMSOL, a finite element solver, for the modelling of reactive transport in soils. With the iCP interface, the equations describing the transport of chemical species in pore water are first solved in COMSOL. The updated chemical analyses for each cell are then sent to Phreeqc for thermodynamic equilibrium and kinetics calculations. In this paper, the iCP interface was tested by modelling a series of laboratory tests conducted with a natural soil (calcareous sand). These tests included a batch test where a solution of lead, copper, cadmium and zinc nitrates was put in contact with the soil, a conservative tracer test with bromide and an elution test with the same metal nitrates as for the batch test. For the batch and elution tests, Phreeqc was used to model the dissolution of calcite in the soil and the precipitation of metal carbonates. COMSOL was used to model solute transport for the tracer and elution tests. This paper is centered on a coupling between hydrodynamic and geochemical phenomena, but the iCP interface could also be applied to multiphysics models involving mechanical (e.g., cementation) and thermal (e.g., heat produced by chemical reactions) equations.