Operator-splitting errors in coupled reactive transport codes for transient variably saturated flow and contaminant transport in layered soil profiles

Research output: Contribution to journalArticlepeer-review

Authors

Institutes & Expert groups

Documents & links

Abstract

This paper evaluates the accuracy of the operator-splitting approach for typical soil environmental problems involving transient atmospheric boundary conditions (precipitation, evapotranspiration) and layered soil profiles. The recently developed HP1 code was used to solve the coupled transport and chemical equations. For steady-state flow conditions, the accuracy was found to be mainly a function of the adopted spatial discretization and to a lesser extent of the temporal discretization. For transient flow situations, the accuracy depended in a complex manner on grid discretization, time stepping and the main flow conditions. Whereas a finer grid size reduced the numerical errors during steady-state flow or the main infiltration periods, the errors sometimes slightly increased (generally less than 50%) when a finer grid size was used during periods with a high evapotranspiration demand (leading to high pressure head gradients near the soil surface). This indicates that operator-splitting errors are most significant during periods with high evaporative boundary conditions. The operator-splitting errors could be decreased by constraining the time step using the performance index (the product of the grid Peclet and Courant numbers) during infiltration, or the maximum time step during evapotranspiration.

Details

Original languageEnglish
Pages (from-to)197-218
JournalJournal of Contaminant Hydrology
Volume88
Issue number3-4
DOIs
Publication statusPublished - Dec 2006

Keywords

  • Coupled reactive transport, Unsaturated flow modelling, Geochemical reactions, Operator-splitting errors

ID: 227833