Finite-element method in modeling geologic transport processes

Deterministic mathematical modeling of complex geologic transport processes may require the use of odd boundary shapes, time dependency, and two or three dimensions. Under these circumstances the governing transport equations must be solved by numerical methods. For a number of transport phenomena a general form of the convective-dispersion equation can be employed. The solution of this equation for complicated problems can be solved readily by the finite-element method. Using quadrilateral isoparametric elements or triangular elements and a computational algorithm based on Galerkin's procedure, solutions to unsteady heat flux from a dike and seawater intrusion in an aquifer have been obtained. These examples illustrate that the finite-element numerical procedure is well suited for solving boundary-value problems resulting from modeling of complex physical phenomena.