Eulerian‐Lagrangian Solution of the Convection‐Dispersion Equation in Natural Coordinates

The vast majority of numerical investigations of transport phenomena use an Eulerian formulation for the convenience that the computational grids are fixed in space. An Eulerian-Lagrangian method (ELM) of solution for the convection-dispersion equation is discussed and analyzed. The ELM uses the Lagrangian concept in an Eulerian computational grid system. The values of the dependent variable off the grid are calculated by interpolation. When a linear interpolation is used, the method is a slight improvement over the upwind difference method. At this level of approximation both the ELM and the upwind difference method suffer from large numerical dispersion. However, if second-order Lagrangian polynomials are used in the interpolation, the ELM is proven to be free of artificial numerical dispersion for the convection-dispersion equation. The concept of the ELM is extended for treatment of anisotropic dispersion in natural coordinates. In this approach the anisotropic properties of dispersion can be conveniently related to the properties of the flow field. Several numerical examples are given to further substantiate the results of the present analysis.

[1]  Jeffrey W. Gartner,et al.  Harmonic analysis of tides and tidal currents in South San Francisco Bay, California , 1985 .

[2]  S. P. Neuman,et al.  A Eulerian-Lagrangian numerical scheme for the dispersion-convection equation using conjugate space-time grids , 1981 .

[3]  J. Bredehoeft,et al.  Computer model of two-dimensional solute transport and dispersion in ground water , 1978 .

[4]  R. T. Cheng,et al.  On Lagrangian residual currents with applications in south San Francisco Bay, California , 1982 .

[5]  R. T. Cheng.  EULER-LAGRANGIAN COMPUTATIONS IN ESTUARINE HYDRODYNAMICS. , 1983 .

[6]  O. C. Zienkiewicz,et al.  An ‘upwind’ finite element scheme for two‐dimensional convective transport equation , 1977 .

[7]  T. Mueller,et al.  Numerical solutions of laminar separated flows , 1970 .

[8]  J. P. Huffenus,et al.  The Lagrangian approach of advective term treatment and its application to the solution of Navier—Stokes equations , 1981 .

[9]  Patrick J. Roache,et al.  On artificial viscosity , 1972 .

[10]  K. O'Neill.  Highly efficient, oscillation free solution of the transport equation over long times and large spaces , 1981 .

[11]  F. Holly,et al.  Accurate Calculation of Transport in Two Dimensions , 1977 .

[12]  S. P. Neuman,et al.  Eulerian-Lagrangian Methods for Advection-Dispersion , 1982 .

[13]  T. F. Russell,et al.  NUMERICAL METHODS FOR CONVECTION-DOMINATED DIFFUSION PROBLEMS BASED ON COMBINING THE METHOD OF CHARACTERISTICS WITH FINITE ELEMENT OR FINITE DIFFERENCE PROCEDURES* , 1982 .