Numerical methods for simultaneous solution of geofiltration and ge- omigration problems in multilayer systems in contaminant transport investigation have been developed. A multilayer system consists of several aquifers separated by weakly permeable layers. Mathematical models are constructed under the following assump- tion: longitudinal flows predominate in the aquifer and crossflows — in the separation layers. The mass transfer is described by the convection and diffusion equations. The implicit finite volume difference schemes were used. The computational al- gorithm is based on the concept of splitting the schemes on physical processes of two- dimensional filtration and contaminant transport along aquifers and vertical transport through the separation levels. An iterative method of the block Gauss — Seidel type for the realization of implicit finite difference schemes for geofiltration and geomigration problems in the presence of wells and contaminant sources has been constructed. The method is based on the perturbation transfer from cells with sources (wells and sources of contaminant). The iterative process for update of unknowns has been constructed taking into account the position of a cell relative to the calculated cells. The results of the numerical experiment on modeling groundwater flows redistri- bution and contaminant migration in system of two aquifers with operating wells are presented.
© Institute of Mathematics, NAS of Belarus
This article is distributed under the terms of the Creative Commons Attribution Non-Commercial License, which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.