Nonlinear Flow Process: A New Package to Compute Nonlinear Flow in MODFLOW

A new MODFLOW package (Nonlinear Flow Process; NLFP) simulating nonlinear flow following the Forchheimer equation was developed and implemented in MODLFOW-2005. The method is based on an iterative modification of the conductance calculated and used by MODFLOW to obtain an effective Forchheimer conductance. The package is compatible with the different layer types, boundary conditions, and solvers as well as the wetting capability of MODFLOW. The correct implementation is demonstrated using four different benchmark scenarios for which analytical solutions are available. A scenario considering transient flow in a more realistic setting and a larger model domain with a higher number of cells demonstrates that NLFP performs well under more complex conditions, although it converges moderately slower than the standard MODFLOW depending on the nonlinearity of flow. Thus, this new tool opens a field of opportunities to groundwater flow simulation with MODFLOW, especially for core sample simulation or vuggy karstified aquifers as well as for nonlinear flow in the vicinity of pumping wells.

[1]  Vassilios A. Tsihrintzis,et al.  Approximate analytical solutions of the Forchheimer equation , 2005 .

[2]  J. Bear.  Hydraulics of Groundwater , 1979 .

[3]  Michael E. Barrett,et al.  Can we simulate regional groundwater flow in a karst system using equivalent porous media models? Case study, Barton Springs Edwards aquifer, USA , 2003 .

[4]  S. Birk,et al.  Understanding changes in the hydrological behaviour within a karst aquifer (Lurbach system, Austria) , 2013, Carbonates and Evaporites.

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

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

[7]  V. Tsihrintzis,et al.  Determination of Forchheimer equation coefficients a and b , 2007 .

[8]  Steffen Birk,et al.  Modifications to the Conduit Flow Process Mode 2 for MODFLOW‐2005 , 2012, Ground water.

[9]  R. Grigg,et al.  A Criterion for Non-Darcy Flow in Porous Media , 2006 .

[10]  T. Geyer,et al.  The significance of turbulent flow representation in single‐continuum models , 2011 .

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

[12]  Haibo Huang,et al.  Evaluation of permeability and non‐Darcy flow in vuggy macroporous limestone aquifer samples with lattice Boltzmann methods , 2013 .

[13]  W Barclay Shoemaker,et al.  Permeameter Data Verify New Turbulence Process for MODFLOW , 2008, Ground water.

[14]  R. Price,et al.  Nonlinear Flow in Karst Formations , 2009, Ground water.

[15]  Guanhua Huang,et al.  Non-Darcian flow to a well in a leaky aquifer using the Forchheimer equation , 2011 .