Solution of the advection-dispersion equation in two dimensions by a finite-volume Eulerian-Lagrangian localized adjoint method

We extend the finite-volume Eulerian-Lagrangian localized adjoint method (FVELLAM) for solution of the advection-dispersion equation to two dimensions. The method can conserve mass globally and is not limited by restrictions on the size of the grid Peclet or Courant number. Therefore, it is well suited for solution of advection-dominated ground-water solute transport problems. In test problem comparisons with standard finite differences, FVELLAM is able to attain accurate solutions on much coarser space and time grids. On fine grids, the accuracy of the two methods is comparable. A critical aspect of FVELLAM (and all other ELLAMs) is evaluation of the mass storage integral from the preceding time level. In FVELLAM this may be accomplished with either a forward or backtracking approach. The forward tracking approach conserves mass globally and is the preferred approach. The backtracking approach is less computationally intensive, but not globally mass conservative. Boundary terms are systematically represented as integrals in space and time which are evaluated by a common integration scheme in conjunction with forward tracking through time. Unlike the one-dimensional case, local mass conservation cannot be guaranteed, so slight oscillations in concentration can develop, particularly in the vicinity of inflow or outflow boundaries.

[1]  Richard E. Ewing,et al.  Eulerian‐Lagrangian localized adjoint method: The theoretical framework , 1993 .

[2]  Eliezer J. Wexler,et al.  Analytical solutions for one-, two-, and three-dimensional solute transport in ground-water systems with uniform flow , 1989 .

[3]  T. F. Russell,et al.  An Eulerian-Lagrangian localized adjoint method for the advection-diffusion equation , 1990 .

[4]  Michael A. Celia,et al.  Eulerian-Lagrangian Localized Adjoint Methods for Contaminant Transport Simulations , 1994 .

[5]  D. W. Pollock.  Semianalytical Computation of Path Lines for Finite‐Difference Models , 1988 .

[6]  Thomas F. Russell.  Computational methods in water resources IX , 1992 .

[7]  John L. Wilson,et al.  Efficient and accurate front tracking for two‐dimensional groundwater flow models , 1991 .

[8]  T. F. Russell,et al.  Eulerian-Lagrangian localized adjoint methods for convection-diffusion equations and their convergence analysis , 1994 .

[9]  W. G. Gray,et al.  Computational methods in water resources X , 1994 .

[10]  T. F. Russell,et al.  NUMERICAL METHODS FOR CONVECTION-DOMINATED DIFFUSION PROBLEMS BASED ON COMBINING THE METHOD OF CHARACTERISTICS WITH FINITE ELEMENT OR FINITE DIFFERENCE PROCEDURES* , 1982 .

[11]  M. Celia,et al.  An Analysis of Some Classes of Petrov-Galerkin and Optimal Test Function Methods , 1988 .

[12]  Tracy Nishikawa,et al.  A New Total Variation Diminishing Scheme for the Solution of Advective‐Dominant Solute Transport , 1991 .

[13]  Endre Süli,et al.  Stability of the Lagrange-Galerkin method with non-exact integration , 1988 .

[14]  Ning Lu,et al.  A SEMIANALYTICAL METHOD OF PATH LINE COMPUTATION FOR TRANSIENT FINITE-DIFFERENCE GROUNDWATER FLOW MODELS , 1994 .

[15]  M. Wheeler,et al.  A characteristics-mixed finite element method for advection-dominated transport problems , 1995 .

[16]  M. A. Celia,et al.  A Eulerian-Lagrangian localized adjoint method for reactive transport in groundwater. , 1990 .

[17]  Philip John Binning,et al.  A Finite Volume Eulerian‐Lagrangian Localized Adjoint Method for Solution of the Contaminant Transport Equations in Two‐Dimensional Multiphase flow Systems , 1996 .

[18]  David W. Pollock,et al.  Documentation of computer programs to compute and display pathlines using results from the U.S. Geological Survey modular three-dimensional finite-difference ground-water flow model , 1989 .

[19]  Richard W. Healy,et al.  A finite‐volume Eulerian‐Lagrangian Localized Adjoint Method for solution of the advection‐dispersion equation , 1993 .