Smoothed Particle Hydrodynamics for anisotropic dispersion in heterogeneous porous media

Abstract

Listen

Translate article

In the context of modeling solute transport through heterogeneous porous media, particle methods possess inherent advantages with respect to mesh-based (Eulerian) methods. In Smoothed Particle Hydrodynamics (SPH), particles represent fluid volumes exchanging concentrations with their neighbors to emulate hydrodynamic dispersion, and advection is simulated by the particles’ displacement. This crucially prevents problems that are otherwise typically associated with Eulerian advection schemes (especially at high grid-Péclet numbers), such as numerical diffusion. Despite the advantages of SPH, modeling dispersion with anisotropic coefficients remains a challenge for the approach, with studies reporting unphysical negative concentrations in conservative problems. This has likely hindered its practical use because dispersion is intrinsically anisotropic in porous media. This article provides a review and numerical evaluation of SPH for simulating dispersion, focusing on three formulations compatible with anisotropic dispersion coefficients. The analysis includes a scheme for which negative concentrations have been formerly reported, plus two more recently developed methods which are applied here for the first time to the problem of anisotropic dispersion in heterogeneous porous media. The SPH schemes are tested under different degrees of dispersion anisotropy for both homogeneous and heterogeneous velocity fields. The results indicate that the newer SPH schemes can produce accurate results without negative concentrations while considering anisotropic dispersion, providing a valid alternative to simulate solute transport through heterogeneous domains.