Journal of Applied Mathematics
Volume 2011 (2011), Article ID 498098, 41 pages
Research Article

Numerical Simulation of Pollutant Transport in Fractured Vuggy Porous Karstic Aquifers

1Department of Mathematics, College of Science, Guizhou University, Huaxi, Guiyang 550025, China
2Computational Transport Phenomena Laboratory (CPTL), Division of Physical Sciences and Engineering (PSE), King Abdullah University of Science and Technology (KAUST), Thuwal 23955-6900, Saudi Arabia

Received 15 December 2010; Accepted 8 January 2011

Academic Editor: Zhangxin Chen

Copyright © 2011 Xiaolin Fan et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


This paper begins with presenting a mathematical model for contaminant transport in the fractured vuggy porous media of a species of contaminant (PCP). Two phases are numerically simulated for a process of contaminant and clean water infiltrated in the fractured vuggy porous media by coupling mixed finite element (MFE) method and finite volume method (FVM), both of which are locally conservative, to approximate the model. A hybrid mixed finite element (HMFE) method is applied to approximate the velocity field for the model. The convection and diffusion terms are approached by FVM and the standard MFE, respectively. The pressure distribution and temporary evolution of the concentration profiles are obtained for two phases. The average effluent concentration on the outflow boundary is obtained at different time and shows some different features from the matrix porous media. The temporal multiscale phenomena of the effluent concentration on the outlet are observed. The results show how the different distribution of the vugs and the fractures impacts on the contaminant transport and the effluent concentration on the outlet. This paper sheds light on certain features of karstic groundwater are obtained.