ABSTRACT
The automated solution of contaminant transport equations inside an aquifer depends on physical and chemical parameters and, most importantly, on the groundwater velocity field. Such a field can be obtained from the numerical solution of a mathematical model describing the hydraulic potential of the saturated flux region. In the pre-processing solution step, geoprocessing techniques were applied to produce computational meshes respecting the geospatial constraints of the initial and boundary conditions to simulate the behavior of the water table around the contaminant. In this work, after the hydraulic potential transference - in terms of velocities field for the transport equation -, a Python numerical code was developed to predict the contaminant migration within an aquifer, describing the amount of pollutant through concentration isolines. Another important contribution of the geoprocessing techniques was observed in the standardization and specification of the georeferenced information from several sources, which provided the geospatial database for the hydrogeologic description of the studied area.
Keywords:
geocomputation; water table; aquifers; dolfin code; ADR model
