Numerical modeling of coupled hydrological phenomena using the Modified Eulerian-Langrangian method

Abstract

The chapter presents the advantages of using the Modified Eulerian-Lagrangian (MEL) method for solving the transport of extensive quantities in a porous medium. In this modified formulation, the material derivatives are written in terms of modified velocities. These are the velocities at which the various phase and component variables propagate in the domain, along their respective characteristic curves. It is shown that these velocities depend on the heterogeneity of various solid matrix and fluid properties. Accordingly, an extension to the Peclet number is presented that also accounts for governing equations that may be advective dominant with no reference to the fluid velocity or even when this velocity is not introduced.A mathematical analysis proves that for coupled partial differential equations (PDEs), unlike the Eulerian implicit finite difference scheme the MEL method unconditionally guarantees the absent of spurious oscillations.The MEL formulation is demonstrated for a coupled set of PDEs concerning the problem of saltwater intrusion, with heat transfer.A numerical example suggests that the MEL scheme produces better resolution compared to the Eulerian and the Eulerian-Lagrangian ones.

Numerical Models of Groundwater Flow and Transport

“…Following Neuman and Sorek (1982) for transport problems (e.g. see also Neuman, 1984;Sorek, 1988a,b), for flow problems (Sorek, 1985;Sorek and Braester, 1988) and a modified Eulerian-Lagrangian method for coupled flow and transport model (Bear et al, 1997;Sorek et al, 2000), a technique consisting of the following two steps is used:…”