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

被引:55
作者
Healy, RW
Russell, TF
机构
[1] US Geological Survey, MS 413, Box 25046, Lakewood
[2] Department of Mathematics, University of Colorado at Denver, Campus Box 170, Denver, CO 80217-3364
基金
美国国家科学基金会;
关键词
finite-volume method; Eulerian-Lagrangian method; advection-dispersion equation; mass conservation; method of characteristics; forward tracking; backtracking;
D O I
10.1016/S0309-1708(96)00033-4
中图分类号
TV21 [水资源调查与水利规划];
学科分类号
081501 ;
摘要
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. Published by Elsevier Science Ltd.
引用
收藏
页码:11 / 26
页数:16
相关论文
共 24 条
[1]  
[Anonymous], NUMER METHODS PARTIA
[2]   A CHARACTERISTICS-MIXED FINITE-ELEMENT METHOD FOR ADVECTION-DOMINATED TRANSPORT PROBLEMS [J].
ARBOGAST, T ;
WHEELER, MF .
SIAM JOURNAL ON NUMERICAL ANALYSIS, 1995, 32 (02) :404-424
[3]   A finite volume Eulerian-Lagrangian localized adjoint method for solution of the contaminant transport equations in two-dimensional multiphase flow systems [J].
Binning, P ;
Celia, MA .
WATER RESOURCES RESEARCH, 1996, 32 (01) :103-114
[4]  
BINNING PJ, 1994, THESIS PRINCETON U
[5]  
BOULOUTAS ET, 1988, P 7 INT C COMP METH, V2, P15
[6]  
CARRANO CS, 1994, COMPUTATIONAL METHOD, V10, P191
[7]   AN EULERIAN-LAGRANGIAN LOCALIZED ADJOINT METHOD FOR THE ADVECTION-DIFFUSION EQUATION [J].
CELIA, MA ;
RUSSELL, TF ;
HERRERA, I ;
EWING, RE .
ADVANCES IN WATER RESOURCES, 1990, 13 (04) :187-206
[8]  
CELIA MA, 1994, COMPUTATIONAL METHOD, V10, P207
[9]  
CELIA MA, 1990, COMPUTATIONAL METHOD, V8, P383
[10]   A NEW TOTAL VARIATION DIMINISHING SCHEME FOR THE SOLUTION OF ADVECTIVE-DOMINANT SOLUTE TRANSPORT [J].
COX, RA ;
NISHIKAWA, T .
WATER RESOURCES RESEARCH, 1991, 27 (10) :2645-2654