Calibrating the exchange coefficient in the modified coupled continuum pipe-flow model for flows in karst aquifers

Summary We investigate the validity of the popular coupled-continuum pipe-flow (CCPF) model for flow in a karst aquifer. The (Navier) Stokes–Darcy model is used as the “true model” for calibrating the exchange coefficient in the CCPF model by minimizing the relative differences between results from the two models or at least by having those differences being below a prescribed threshold value. We find that although the CCPF model is never in perfect agreement with the Stokes–Darcy model, there is an almost universal choice for a nearly optimal exchange coefficient such that the relative error is below one percent. Our numerics suggest that the nearly optimal choice of the exchange coefficient should be sufficiently large instead of being a small quantity that is proportional to the hydraulic conductivity, as suggested in existing literatures. We also show that this nearly optimal choice of exchange coefficient is robust under a wide range of model parameters. This result demonstrates that the CCPF model is a valid approximation for flows in karst aquifers as long as we set the fluid exchange coefficient sufficiently large and at least in the simple two-dimensional setting that we consider.

[1]  Eve L. Kuniansky U.S. Geological Survey Karst Interest Group Proceedings, Bowling Green, Kentucky, May 27-29, 2008 , 2008 .

[2]  M. Gunzburger,et al.  Coupled Stokes-Darcy model with Beavers-Joseph interface boundary condition , 2010 .

[3]  P. Saffman On the Boundary Condition at the Surface of a Porous Medium , 1971 .

[4]  B. Katz,et al.  Changes in the isotopic and chemical composition of ground water resulting from a recharge pulse from a sinking stream , 1998 .

[5]  Georg Teutsch,et al.  Hydraulic boundary conditions as a controlling factor in karst genesis: A numerical modeling study on artesian conduit development in gypsum , 2003 .

[6]  D. Joseph,et al.  Boundary conditions at a naturally permeable wall , 1967, Journal of Fluid Mechanics.

[7]  Yanzhao Cao,et al.  Analysis and finite element approximation of a coupled, continuum pipe‐flow/Darcy model for flow in porous media with embedded conduits , 2011 .

[8]  Paul Williams,et al.  Karst Geomorphology and Hydrology , 1989 .

[9]  Sebastian Bauer,et al.  Documentation of a Conduit Flow Process (CFP) for MODFLOW-2005 , 2007 .

[10]  Donald A. Nield,et al.  Onset of convection in a fluid layer overlying a layer of a porous medium , 1977, Journal of Fluid Mechanics.

[11]  B. Ji,et al.  The media and movement of karst water , 1988 .

[12]  Fei Hua,et al.  Modeling, analysis and simulation of the Stokes -Darcy system with Beavers -Joseph interface condition , 2009 .

[13]  Ivan Yotov,et al.  Coupling Fluid Flow with Porous Media Flow , 2002, SIAM J. Numer. Anal..

[14]  Xiaoming Wang,et al.  On the coupled continuum pipe flow model (CCPF) for flows in karst aquifer , 2009 .

[15]  Georg Teutsch,et al.  Simulation of the development of karst aquifers using a coupled continuum pipe flow model , 2003 .

[16]  A. W. Harbaugh MODFLOW-2005 : the U.S. Geological Survey modular ground-water model--the ground-water flow process , 2005 .

[17]  Bill X. Hu,et al.  Laboratory analog and numerical study of groundwater flow and solute transport in a karst aquifer with conduit and matrix domains. , 2009, Journal of contaminant hydrology.

[18]  M. Sauter Quantification and forecasting of regional groundwater flow and transport in a karst aquifer (Gallusquelle, Malm, SW. Germany) , 1992 .

[19]  Sebastian Bauer,et al.  Modelling of karst development considering conduit-matrix exchange flow , 2000 .

[20]  G. I. Barenblatt,et al.  Basic concepts in the theory of seepage of homogeneous liquids in fissured rocks [strata] , 1960 .

[21]  E. Greene,et al.  Quantitative Approaches in Characterizing Karst Aquifers , 2001 .

[22]  G. Teutsch,et al.  A combined continuum and discrete network reactive transport model for the simulation of karst development , 1996 .

[23]  Sebastian Bauer,et al.  Modeling of karst aquifer genesis: Influence of exchange flow , 2003 .

[24]  E. Sudicky,et al.  On the Incorporation of Drains into Three‐Dimensional Variably Saturated Groundwater Flow Models , 1996 .

[25]  E. Miglio,et al.  Mathematical and numerical models for coupling surface and groundwater flows , 2002 .