Premium
An alternating direction galerkin technique for simulation of contaminant transport in complex groundwater systems
Author(s) -
Daus A. D.,
Frind E. O.
Publication year - 1985
Publication title -
water resources research
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.863
H-Index - 217
eISSN - 1944-7973
pISSN - 0043-1397
DOI - 10.1029/wr021i005p00653
Subject(s) - galerkin method , computation , tridiagonal matrix , finite element method , grid , representation (politics) , hydraulic conductivity , mathematics , mathematical analysis , mathematical optimization , computer science , geometry , algorithm , geology , engineering , soil science , physics , structural engineering , eigenvalues and eigenvectors , quantum mechanics , politics , political science , law , soil water
A general contaminant transport simulation technique based on a Galerkin finite element representation in space and an alternating direction time‐stepping scheme is developed. The formulation, which is in terms of natural coordinates following the principal directions of hydraulic conductivity of the medium, decouples the spatial components of the equations, leading to tridiagonal matrices. The numerical accuracy can be easily controlled in each principal direction through the well‐established Peclet and Courant numbers. Results are found to be comparable to those from the classic finite element scheme, but the cost of computation is generally much less because the computational effort varies only linearly with the grid size. The technique is thus well suited for transport simulations with complex layered anisotropic groundwater systems requiring large and detailed grids.