Comparison of numerical methods for simulating strongly nonlinear and heterogeneous reactive transport problems-the MoMaS benchmark case

Carrayrou J, Hoffmann J, Knabner P, Kräutle S, de Dieuleveult C, Erhel J, Van der Lee J, Lagneau V, Mayer KU, MacQuarrie KTB (2010)

Publication Status: Published

Publication Type: Journal article, Original article

Publication year: 2010

Journal

Computational Geosciences Springer Verlag (Germany)

Publisher: Springer Verlag (Germany)

Book Volume: 14

Pages Range: 483-502

Journal Issue: 3

URI: https://www1.am.uni-erlangen.de/research/publications/Jahr_2010/2010_CarrayrouHoffmannKnKraeutleDieuleveultErhel_ComparisonOfNumericalMethodsForSimulStrongNonlinearAndHetere

DOI: 10.1007/s10596-010-9178-2

Abstract

Although multicomponent reactive transport modeling is gaining wider application in various geoscience fields, it continues to present significant mathematical and computational challenges. There is a need to solve and compare the solutions to complex benchmark problems, using a variety of codes, because such intercomparisons can reveal promising numerical solution approaches and increase confidence in the application of reactive transport codes. In this contribution, the results and performance of five current reactive transport codes are compared for the 1D and 2D subproblems of the so-called easy test case of the MoMaS benchmark (Carrayrou et al., Comput Geosci, 2009, this issue). This benchmark presents a simple fictitious reactive transport problem that highlights the main numerical difficulties encountered in real reactive transport problems. As a group, the codes include iterative and noniterative operator splitting and global implicit solution approaches. The 1D easy advective and 1D easy diffusive scenarios were solved using all codes, and, in general, there was a good agreement, with solution discrepancies limited to regions with rapid concentration changes. Computational demands were typically consistent with what was expected for the various solution approaches. The differences between solutions given by the three codes solving the 2D problem are more important. The very high computing effort required by the 2D problem illustrates the importance of parallel computations. The most important outcome of the benchmark exercise is that all codes are able to generate comparable results for problems of significant complexity and computational difficulty. © Springer Science+Business Media B.V. 2010.

Authors with CRIS profile

Peter Knabner Lehrstuhl für Angewandte Mathematik (Modellierung und Numerik) Serge Kräutle Lehrstuhl für Angewandte Mathematik (Modellierung und Numerik)

Involved external institutions

Mines ParisTech / École Nationale Supérieure des Mines de Paris (ENSMP) FR France (FR) University of New Brunswick (UNB) CA Canada (CA) INRIA Rennes - Bretagne Atlantique FR France (FR) Université de Strasbourg (UDS) FR France (FR) University of British Columbia CA Canada (CA)

How to cite

APA:

Carrayrou, J., Hoffmann, J., Knabner, P., Kräutle, S., de Dieuleveult, C., Erhel, J.,... MacQuarrie, K.T.B. (2010). Comparison of numerical methods for simulating strongly nonlinear and heterogeneous reactive transport problems-the MoMaS benchmark case. Computational Geosciences, 14(3), 483-502. https://doi.org/10.1007/s10596-010-9178-2

MLA:

Carrayrou, Jérôme, et al. "Comparison of numerical methods for simulating strongly nonlinear and heterogeneous reactive transport problems-the MoMaS benchmark case." Computational Geosciences 14.3 (2010): 483-502.

BibTeX: Download