Multiphase flow and reactive transport benchmark for radioactive waste disposal

Compacted bentonite is part of the multi-barrier system of radioactive waste repositories. The assessment of the long-term performance of the barrier requires using reactive transport models. Here we present a multiphase flow and reactive transport benchmark for radioactive waste disposal. The numerical model deals with a 1D column of unsaturated bentonite through which water, dry air and CO2(g)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\hbox {CO}_{2{(g)}}}$$\end{document} may flow and with the following reactions; aqueous complexation, calcite and gypsum dissolution/precipitation, cation exchange and gas dissolution. INVERSE-FADES-CORE V2, DuMuX\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {DuMu}<^>X$$\end{document}, TOUGHREACT and iCP were benchmarked with 6 test cases of increasing complexity, starting with conservative tracer transport under variably unsaturated conditions and ending with water flow, gas diffusion, minerals and cation exchange. The solutions of all codes exhibit similar trends. Small discrepancies are found in conservative tracer transport due to differences in hydrodynamic dispersion. Computed CO2(g)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\hbox {CO}_{2{(g)}}}$$\end{document} pressures agree when a sufficiently refined grid is used. Small discrepancies in CO2(g)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\hbox {CO}_{2{(g)}}}$$\end{document} and pH are found near the no-flow boundary at early times which vanish later. Discrepancies are due differences in the formulations used for gas flow at nearly water-saturated conditions. Computed CO2(g)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\hbox {CO}_{2{(g)}}}$$\end{document} pressures show a fluctuation between 10-4\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$10<^>{-4}$$\end{document} and 10-3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$10<^>{-3}$$\end{document} years which slows down the in-diffusion of CO2(g)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\hbox {CO}_{2{(g)}}}$$\end{document}. This fluctuation is associated with chemical reactions involving CO2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\hbox {CO}_{2}}$$\end{document}. There are discrepancies in solute concentrations due to differences in the Debye-H & uuml;ckel (DH) formulation. They are overcome when all codes use the same DH formulation. The results of this benchmark will contribute to increase the confidence on multiphase reactive transport models for radioactive waste disposal.