Practical finite analytic methods for simulation of advection-dominated solute transport

Abstract

The practical finite-analytic (PFA) method was applied to the solution of the one-dimensional advection–dispersion equation (ADE) for solute transport in porous media under advection-dominated (high Peclet number, Pe) conditions. Several PFA spatial-temporal computational molecules were developed for Cauchy and pulse loading boundary conditions. The PFA solutions were compared with solutions from the upwind method and quadratic upwind differencing (QUICK) scheme. For all boundary conditions the trapezoidal explicit PFA (EPFA) computational molecule gave the most accurate results at very high Pe number as long as the Courant number (Cr) was close to one. Stability analysis shows that the PFA molecules are always stable for high Pe number.