An assembly model for simulation of large-scale ground water flow and transport.

When managing large-scale ground water contamination problems, it is often necessary to model flow and transport using finely discretized domains--for instance (1) to simulate flow and transport near a contamination source area or in the area where a remediation technology is being implemented; (2) to account for small-scale heterogeneities; (3) to represent ground water-surface water interactions; or (4) some combination of these scenarios. A model with a large domain and fine-grid resolution will need extensive computing resources. In this work, a domain decomposition-based assembly model implemented in a parallel computing environment is developed, which will allow efficient simulation of large-scale ground water flow and transport problems using domain-wide grid refinement. The method employs common ground water flow (MODFLOW) and transport (RT3D) simulators, enabling the solution of almost all commonly encountered ground water flow and transport problems. The basic approach partitions a large model domain into any number of subdomains. Parallel processors are used to solve the model equations within each subdomain. Schwarz iteration is applied to match the flow solution at the subdomain boundaries. For the transport model, an extended numerical array is implemented to permit the exchange of dispersive and advective flux information across subdomain boundaries. The model is verified using a conventional single-domain model. Model simulations demonstrate that the proposed model operated in a parallel computing environment can result in considerable savings in computer run times (between 50% and 80%) compared with conventional modeling approaches and may be used to simulate grid discretizations that were formerly intractable.

[1]  Third-order multi-dimensional Euler/Navier-Stokes solver , 1988 .

[2]  W. Kessels,et al.  A RAM-economizing domain decomposition technique for regional high-resolution groundwater simulation , 2003 .

[3]  Paul P. Wang,et al.  MT3DMS: A Modular Three-Dimensional Multispecies Transport Model for Simulation of Advection, Dispersion, and Chemical Reactions of Contaminants in Groundwater Systems; Documentation and User's Guide , 1999 .

[4]  Andrea Toselli,et al.  Recent developments in domain decomposition methods , 2002 .

[5]  Rolf Hempel,et al.  The MPI Message Passing Interface Standard , 1994 .

[6]  Message P Forum,et al.  MPI: A Message-Passing Interface Standard , 1994 .

[7]  Yu-Shu Wu,et al.  An efficient parallel-computing method for modeling nonisothermal multiphase flow and multicomponent transport in porous and fractured media , 2002 .

[8]  Arlen W. Harbaugh,et al.  A modular three-dimensional finite-difference ground-water flow model , 1984 .

[9]  Stanley A. Leake,et al.  A new ghost-node method for linking different models and initial investigations of heterogeneity and nonmatching grids , 2007 .

[10]  B. P. Leonard,et al.  Cost-effective accurate coarse-grid method for highly convective multidimensional unsteady flows , 1991 .

[11]  Gary B. Lamont,et al.  Parallel real-valued genetic algorithms for bioremediation optimization of TCE-contaminated groundwater , 1999, Proceedings of the 1999 Congress on Evolutionary Computation-CEC99 (Cat. No. 99TH8406).

[12]  T. P. Clement,et al.  Modeling Multispecies Reactive Transport in Ground Water , 1998 .

[13]  Barbara S. Minsker,et al.  Multiscale island injection genetic algorithms for groundwater remediation , 2007 .

[14]  S. Wondzell,et al.  Geomorphic controls on hyporheic exchange flow in mountain streams , 2003 .

[15]  Erik Elmroth,et al.  High Performance Computations for Large Scale Simulations of Subsurface Multiphase Fluid and Heat Flow , 2004, The Journal of Supercomputing.

[16]  Vedat Batu,et al.  Applied Flow and Solute Transport Modeling in Aquifers: Fundamental Principles and Analytical and Numerical Methods , 2005 .

[17]  Mary C. Hill,et al.  Three-dimensional local grid refinement for block-centered finite-difference groundwater models using iteratively coupled shared nodes: A new method of interpolation and analysis of errors , 2004 .

[18]  Arlen W. Harbaugh,et al.  MODFLOW-2000, The U.S. Geological Survey Modular Ground-Water Model - User Guide to Modularization Concepts and the Ground-Water Flow Process , 2000 .

[19]  Brian Sanderson,et al.  Order and resolution for computational ocean dynamics , 1998 .

[20]  H. H. Rachford,et al.  The Numerical Solution of Parabolic and Elliptic Differential Equations , 1955 .

[21]  Roger S. Crouch,et al.  Advantages of using adaptive remeshing and parallel processing for modelling biodegradation in groundwater , 2005 .

[22]  C. T. Kelley,et al.  An Aggregation-Based Domain Decomposition Preconditioner for Groundwater Flow , 2001, SIAM J. Sci. Comput..

[23]  Timothy J. Williams,et al.  A Numerical Simulation of Groundwater Flow and Contaminant Transport on the CRAY T3D and C90 Supercomputers , 1999, Int. J. High Perform. Comput. Appl..

[24]  Onyx Wing-Hong Wai,et al.  An efficient parallel model for coastal transport process simulation , 2000 .

[25]  Steffen Mehl,et al.  Evaluation of a local grid refinement method for steady-state block-centered finite-difference groundwater models , 2002 .

[26]  D. Dudley Williams,et al.  Factors controlling riffle‐scale hyporheic exchange flows and their seasonal changes in a gaining stream: A three‐dimensional groundwater flow model , 2003 .

[27]  Erik Elmroth,et al.  A Parallel Implementation of the TOUGH2 Software Package for Large Scale Multiphase Fluid and Heat Flow Simulations , 1999, SC.