Efficient Parallel Implementation of Hybrid Optimization Approaches for Solving Groundwater Inverse Problems

Inverse problems that are constrained by large-scale partial differential equation (PDE) systems demand very large computational resources. Solutions to these problems generally require the solution of a large number of complex PDE systems. Three-dimensional groundwater inverse problems fall under this category. In this paper, we describe the implementation of a parallel simulation-optimization framework for solving PDE-based inverse problems and demonstrate it for the solution of groundwater contaminant source release history reconstruction problems that are of practical importance. The optimization component employs several optimization algorithms, including genetic algorithms (GAs) and several local search (LS) approaches that can be used in a hybrid mode. This hybrid GA-LS optimizer is used to drive a parallel finite-element (FEM) groundwater forward transport simulator. Parallelism is exploited within the transport simulator (fine grained parallelism) as well as the optimizer (coarse grained parallelism) through the exclusive use of the Message Passing Interface (MPI) communication library. Algorithmic and parallel performance results are presented for an IBM SP3 supercomputer. Simulation and performance results presented in this paper illustrate that an effective combination of efficient optimization algorithms and parallel computing can enable solution to three-dimensional groundwater inverse problems of a size and complexity not attempted before.

[1]  Anthony Skjellum,et al.  Using MPI - portable parallel programming with the message-parsing interface , 1994 .

[2]  Mohamed Sayeed,et al.  Hybrid optimization approaches for solving groundwater inverse problems in a parallel computing environment , 2003 .

[3]  Tom V. Mathew,et al.  Transit route network design using parallel genetic algorithm , 2004 .

[4]  Lawrence. Davis,et al.  Handbook Of Genetic Algorithms , 1990 .

[5]  L. L. Rogers,et al.  Optimization of groundwater remediation using artificial neural networks with parallel solute transport modeling , 1994 .

[6]  Mustafa M. Aral,et al.  Identification of Contaminant Source Location and Release History in Aquifers , 2001 .

[7]  Mohamed Sayeed.  A Parallel Optimization Framework for Inverse Problems , 2004 .

[8]  Ganesh Thiagarajan,et al.  Parallelization strategies for element-by-element preconditioned conjugate gradient solver using high-performance Fortran for unstructured finite-element applications on linux clusters , 2002 .

[9]  L. L. Hamm,et al.  Data fusion modeling for groundwater systems , 2000 .

[10]  G. Pinder,et al.  Computational Methods in Subsurface Flow , 1983 .

[11]  N. Sun.  Inverse problems in groundwater modeling , 1994 .

[12]  Z. J. Kabala,et al.  Recovering the release history of a groundwater contaminant using a non-linear least-squares method. , 2000 .

[13]  E. Poeter,et al.  Inverse Models: A Necessary Next Step in Ground‐Water Modeling , 1997 .

[14]  B. H. V. Topping,et al.  Transient Dynamic Nonlinear Analysis Using MIMD Computer Architectures , 2000 .

[15]  Mustafa M. Aral,et al.  Genetic Algorithms in Search of Groundwater Pollution Sources , 1996 .

[16]  Bithin Datta,et al.  Optimal Monitoring Network and Ground-Water–Pollution Source Identification , 1997 .

[17]  Bithin Datta,et al.  Identification of Pollution Sources in Transient Groundwater Systems , 2000 .

[18]  K. R. Bailey,et al.  Estimation of groundwater flow parameters using least squares , 1997 .

[19]  William H. Press,et al.  Numerical recipes in Fortran 90: the art of parallel scientific computing, 2nd Edition , 1996, Fortran numerical recipes.

[20]  Yong Liu,et al.  Full multiscale approach for optimal control of in situ bioremediation , 2004 .

[21]  A. Bagtzoglou,et al.  State of the Art Report on Mathematical Methods for Groundwater Pollution Source Identification , 2001 .

[22]  M. C. Balla,et al.  Distributed Genetic Algorithm Model on Network of Personal Computers , 2000 .

[23]  Bithin Datta,et al.  Optimal Identification of Ground-Water Pollution Sources and Parameter Estimation , 2001 .

[24]  Jack Dongarra,et al.  PVM: Parallel virtual machine: a users' guide and tutorial for networked parallel computing , 1995 .

[25]  G. Mahinthakumar,et al.  A Hybrid Mpi-Openmp Implementation of an Implicit Finite-Element Code on Parallel Architectures , 2002, Int. J. High Perform. Comput. Appl..

[26]  G. Mahinthakumar,et al.  Hybrid Genetic Algorithm—Local Search Methods for Solving Groundwater Source Identification Inverse Problems , 2005 .

[27]  G. Mahinthakumar,et al.  Implementation and Performance Analysis of a Parallel Multicomponent Groundwater Transport Code , 1999, PPSC.

[28]  David E. Goldberg,et al.  Genetic Algorithms in Search Optimization and Machine Learning , 1988 .

[29]  William H. Press,et al.  Numerical Recipes: FORTRAN , 1988 .

[30]  G. Mahinthakumar,et al.  Subsurface Biological Activity Zone Detection Using Genetic Search Algorithms , 1999 .

[31]  G. Mahinthakumar,et al.  Comparison of some parallel Krylov solvers for large scale groundwater contaminant transport simulations , 1997 .

[32]  T. Edgar,et al.  The finite layer method for groundwater flow models , 1992 .

[33]  J. Bear.  Dynamics of Fluids in Porous Media , 1975 .

[34]  Robert D. Morrison,et al.  Application of Forensic Techniques for Age Dating and Source Identification in Environmental Litigation , 2000 .

[35]  Singiresu S. Rao.  Engineering Optimization : Theory and Practice , 2010 .