A parallel mesh-free contaminant transport model based on the Analytic Element and Streamline Methods

This paper introduces a new method for simulating large-scale subsurface contaminant transport that combines an Analytic Element Method (AEM) groundwater flow solution with a split-operator Streamline Method for modeling reactive transport. The key feature of the method is the manner in which the vertically integrated AEM flow solution is used to construct three-dimensional particle tracks that define the geometry of the Streamline Method. The inherently parallel nature of the algorithm supports the development of reactive transport models for spatial domains much larger than current grid-based methods. The applicability of the new approach is verified for cases with negligible transverse dispersion through comparisons to analytic solutions and existing numerical solutions, and parallel performance is demonstrated through a realistic test problem based on the regional-scale transport of agricultural contaminants from spatially distributed sources.