Numerical errors of explicit finite difference approximation for two-dimensional solute transport equation with linear sorption

The numerical errors associated with explicit upstream finite difference solutions of two-dimensional advection-dispersion equation with linear sorption are formulated from a Taylor analysis. The error expressions are based on a general form of the corresponding difference equation. The numerical truncation errors are defined using Peclet and Courant numbers in the X and Y direction, a sink/source dimensionless number and new Peclet and Courant numbers in the XY plane. The effects of these truncation errors on the explicit solution of a two-dimensional advection-dispersion equation with a first-order reaction or degradation are demonstrated by comparison with an analytical solution in uniform flow field. The results show that these errors are not negligible and correcting the finite difference scheme for them results in a more accurate solution. (c) 2004 Elsevier Ltd. All rights reserved.