Testing Accuracy of Finite Element and Random Walk Schemes in Prediction of Pollutant Dispersion in Coastal Waters

Distribution of pollutants in coastal waters is usually represented by depth averaged twodimensional convection-dispersion equation. Under very specific conditions this equation can be solved analytically. Although such a solution is restricted to simplified situations it provides a very useful case for testing the performance of various numerical solution techniques currently available for the simulation of convective-dispersion of pollutants in natural water systems. In this paper the analytical solution of the convective dispersion equation is used as a benchmark against which the accuracy of other techniques are assessed. These assessments are based on quantitative comparisons between the results of the solution of two-dimensional convection-dispersion equation by the deterministic finite element and stochastic random walk methods. Both Eulerian and Lagrangian frameworks are employed to obtain the finite element solution of the convection-dispersion problem. It has been shown that the Lagrange–Galerkin finite element scheme yields the most accurate results for the case under study. However, computational costs of the Lagrange–Galerkin method can be relatively high and under certain conditions it may be reasonable to use a less accurate but cost effective random walk scheme to make water quality management decisions.