Header menu link for other important links
X
Overlapping control volume method for solute transport
, Atul Kumar Verma, V. Eswaran
Published in ASCE, Reston, VA, United States
2000
Volume: 5
   
Issue: 3
Pages: 308 - 316
Abstract
In this paper, an overlapping control volume method is presented for the numerical solution of transient 2D solute transport problems in ground water. The method is applicable for nonorthogonal grids and uses an isoparametric formulation for computing the dispersion and for second-order upwinding. Time integration is performed using an implicit approach. Three test cases are considered for comparing the numerical with analytical solutions. The scheme is second order in space and, when combined with Crank-Nicholson, is also second order in time. The results using Crank-Nicholson and full implicit (first-order) time integration methods are compared for problems with a variety of boundary conditions. For diffusion-dominated flows (PΔ ≤ 2) Crank-Nicholson works well, but for convection-dominated flows it produces spurious oscillations due to numerical dispersion errors. These oscillations are controlled by a flux limiter but only for Courant numbers below unity. It is shown that for high Courant and Peclet numbers a slight weighting of the time stepping toward fully implicit is effective against spurious oscillations and offers an optimum compromise between numerical dissipation and dispersion errors for a wide range of Courant and Peclet numbers. The scheme is shown to work on mildly nonorthogonal grids.In this paper, an overlapping control volume method is presented for the numerical solution of transient 2D solute transport problems in ground water. The method is applicable for nonorthogonal grids and uses an isoparametric formulation for computing the dispersion and for second-order upwinding. Time integration is performed using an implicit approach. Three test cases are considered for comparing the numerical with analytical solutions. The scheme is second order in space and, when combined with Crank-Nicholson, is also second order in time. The results using Crank-Nicholson and full implicit (first-order) time integration methods are compared for problems with a variety of boundary conditions. For diffusion-dominated flows (PΔ ≤ 2) Crank-Nicholson works well, but for convection-dominated flows it produces spurious oscillations due to numerical dispersion errors. These oscillations are controlled by a flux limiter but only for Courant numbers below unity. It is shown that for high Courant and Peclet numbers a slight weighting of the time stepping toward fully implicit is effective against spurious oscillations and offers an optimum compromise between numerical dissipation and dispersion errors for a wide range of Courant and Peclet numbers. The scheme is shown to work on mildly nonorthogonal grids.
About the journal
JournalJournal of Hydrologic Engineering
PublisherASCE, Reston, VA, United States
ISSN10840699
Open AccessYes
Concepts (12)
  •  related image
    Algorithms
  •  related image
    Boundary conditions
  •  related image
    Flow of water
  •  related image
    Integration
  •  related image
    Numerical methods
  •  related image
    Two dimensional
  •  related image
    COURANT NUMBERS
  •  related image
    OVERLAPPING CONTROL VOLUME METHOD
  •  related image
    PECLET NUMBERS
  •  related image
    Solute transport
  •  related image
    Groundwater
  •  related image
    Numerical method