An Aggregation-Based Domain Decomposition Preconditioner for Groundwater Flow

We consider theoretical and computational issues associated with an aggregation-based domain decomposition preconditioner applied to a Bi-CGSTAB iterative solver used to solve both Laplace's equation and an important nonlinear model from hydrology used to simulate unsaturated flow, Richards' equation. Theoretical results for Laplace's equation provide estimates of the condition number and the rate of convergence for a two-level Schwarz domain decomposition preconditioner. Computational results for Laplace's equation and Richards' equation show excellent scalability, although no theory is yet available to support the results for the difficult nonlinear problem.

[1]  Henk A. van der Vorst,et al.  Bi-CGSTAB: A Fast and Smoothly Converging Variant of Bi-CG for the Solution of Nonsymmetric Linear Systems , 1992, SIAM J. Sci. Comput..

[2]  L. A. Richards.  Capillary conduction of liquids through porous mediums , 1931 .

[3]  Van Genuchten,et al.  A closed-form equation for predicting the hydraulic conductivity of unsaturated soils , 1980 .

[4]  Christopher E. Kees,et al.  C++ implementations of numerical methods for solving differential-algebraic equations: design and optimization considerations , 1999, TOMS.

[5]  Barry F. Smith,et al.  Domain Decomposition: Parallel Multilevel Methods for Elliptic Partial Differential Equations , 1996 .

[6]  Stacy E. Howington,et al.  Newton-Krylov-Schwarz methods for Richards' equation , 1999 .

[7]  Jinchao Xu,et al.  Some Nonoverlapping Domain Decomposition Methods , 1998, SIAM Rev..

[8]  J. Mandel.  Balancing domain decomposition , 1993 .

[9]  Y. Mualem.  A New Model for Predicting the Hydraulic Conductivity , 1976 .

[10]  Linda R. Petzold,et al.  Using Krylov Methods in the Solution of Large-Scale Differential-Algebraic Systems , 1994, SIAM J. Sci. Comput..

[11]  Linda R. Petzold,et al.  Numerical solution of initial-value problems in differential-algebraic equations , 1996, Classics in applied mathematics.

[12]  Marian Brezina,et al.  Robust Iterative Methods On Unstructured Meshes , 1997 .

[13]  Marian Brezina,et al.  A Black-Box Iterative Solver Based on a Two-Level Schwarz Method , 1999, Computing.

[14]  Alan C. Hindmarsh,et al.  Description and use of LSODE, the Livermore Solver for Ordinary Differential Equations , 1993 .

[15]  Roland W. Freund,et al.  A Transpose-Free Quasi-Minimal Residual Algorithm for Non-Hermitian Linear Systems , 1993, SIAM J. Sci. Comput..

[16]  G. M.,et al.  Partial Differential Equations I , 2023, Applied Mathematical Sciences.

[17]  J. Mandel,et al.  Hybrid Domain Decomposition with Unstructured Subdomains , .

[18]  Jun Zou,et al.  Boundary Treatments for Multilevel Methods on Unstructured Meshes , 1999, SIAM J. Sci. Comput..

[19]  C. Kelley.  Iterative Methods for Linear and Nonlinear Equations , 1987 .

[20]  Jinchao Xu,et al.  Iterative Methods by Space Decomposition and Subspace Correction , 1992, SIAM Rev..

[21]  Stacy E. Howington,et al.  A two-level aggregation-based newton-Krylov-Schwarz method for hydrology , 1999 .

[22]  Stacy E. Howington,et al.  A model to simulate the interaction between groundwater and surface water , 1999 .

[23]  Y. Saad,et al.  GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems , 1986 .