A Dual Continuum Model for Groundwater Flow in Karst Aquifers

R. Painter[1], J. Harris[1], T. Byl[1], L. Sharpe[1]
[1]Tennessee State University, Nashville, TN, USA
Published in 2014

In karst terrain, the interaction between the fast moving, highly turbulent conduit flow and the slow moving, laminar flow through porous media is largely unknown. Yet the dual nature of the karst system impacts the residence time of solutes in the groundwater supplies of over 25% of the world’s population. This project attempts to address the need for a better understanding of the dual flow patterns in karst aquifer systems. Because of the highly disparate nature of the two flow regimes, a dual-continuum model was produced to predict storage in porous media as a result of the increased flow due to storm events. COMSOL Multiphysics® software and specifically the Forchheimer Flow Module. We extended the Forchheimer Flow Module by transforming its geometry to reflect an actual karst aquifer conduit and by the coupling the existing model with the physics for contaminant transport. The model was developed as a staged solution in that it solves for the steady state velocity and pressure fields and then uses these fields as inputs for the second stage solution of the advection dispersion equation describing contaminant transport. The steady state velocity field is depicted in Figure 1. Initial testing of the model demonstrates that it is capable of modeling a wide variety of karst tracer response behavior ranging from near Gaussian or symmetric responses to bimodal responses.