AN EFFICIENT EULERIAN-LAGRANGIAN METHOD FOR SOLVING SOLUTE TRANSPORT PROBLEMS IN STEADY AND TRANSIENT FLOW-FIELDS

A computationally efficient, yet relatively simple Eulerian-Lagrangian method is proposed for solving the one-dimensional convection-dispers ion solute transport equation assuming a steady or transient velocity field. The method uses a modified single-step reverse particle trackin g (MSRPT) technique to handle steep concentration fronts. The scheme u tilizes two weighting factors to control the movement of particles dur ing a backward tracking step. One weighting factor greater than unity is used in the upstream region of the convection front, while another weighting factor less than unity is taken in the downstream region. Th e two factors were related empirically to the grid Peclet and Courant numbers. The MSRPT technique is carried out only within the concentrat ion plume at each time step. For transient flow fields, the weighting factors were determined using an automatically adjustable procedure ba sed on mass balance errors. The MSRPT method maintains the advantages of the traditional single-step reverse particle tracking (SRPT) proced ure, i.e., producing efficient and oscillation-free calculations, but circumvents numerical dispersion introduced by SRPT. A large number of tests against analytical solutions for one-dimensional transport in u niform flow fields indicate that the proposed method can handle the en tire range of Peclet numbers from zero to infinity. Numerical tests al so show that the MSRPT method is a relatively accurate, efficient and mass-conservative algorithm for solute transport in transient flow fie lds. The Courant number at present cannot exceed 1. The MSRPT approach was found especially useful for convection-dominated problems; in fac t, an exact numerical solution may be obtained with MSRPT for pure con vection.