The impact of inertial effects on solute dispersion in a channel with periodically varying aperture

Abstract

We investigate solute transport in channels with a periodically varying aperture, when the flow is still laminar but sufficiently fast for inertial effects to be nonnegligible. The flow field is computed for a two-dimensional setup using a finite element analysis, while transport is modeled using a random walk particle tracking method. Recirculation zones are observed when the aspect ratio of the unit cell and the relative aperture fluctuations are sufficiently large; under non-Stokes flow conditions, the flow in non-reversible, which is clearly noticeable by the horizontal asymmetry in the recirculation zones. After characterizing the size and position of the recirculation zones as a function of the geometry and Reynolds number, we investigate the corresponding behavior of the longitudinal effective diffusion coefficient. We characterize its dependence on the molecular diffusion coefficient D m , the Péclet number, the Reynolds number, and the geometry. The proposed relation is a generalization of the well-known Taylor-Aris relationship relating the longitudinal dispersion coefficient to D m and the Péclet number for a channel of constant aperture at sufficiently low Reynolds number. Inertial effects impact the exponent of the Péclet number in this relationship; the exponent is controlled by the relative amplitude of aperture fluctuations. For the range of parameters investigated, the measured dispersion coefficient always exceeds that corresponding to the parallel plate geometry under Stokes conditions; in other words, boundary fluctuations always result in increased dispersion. The transient approach to the asymptotic regime is also studied and characterized quantitatively. We show that the measured characteristic time to attain asymptotic conditions is controlled by two competing effects: (i) the trapping of particles in the near-immobile zone and, (ii) the enhanced mixing in the central zone where most of the flow takes place (mainstream), due to its thinning. C 2012 American Institute of Physics. [http://dx

Pore-Scale Modeling of Mixing-Induced Reaction Transport through a Single Self-Affine Fracture

“…The main flow channel becomes narrow due to the expansion of the eddy domain. Several previous works [28][29][30][31][32] showed similar results for the growth of eddies due to the increasing Reynolds number. At high Reynolds number (Re=170), it can be seen clearly that the eddy-controlled domain occupies a large portion of the fracture.…”

Tail shortening with developing eddies in a rough‐walled rock fracture

“…However, unlike previous numerical studies, our experiments showed that particles that entered eddies were relatively quickly cast back into the main flow channel. That is, there was no well-defined separation stream surface delineating so-called immobile (recirculation) and mobile zones as is predicted by purely two-dimensional simulations [Boutt et al, 2006;Cardenas et al, 2007;Bouquain et al, 2012]. This phenomenon became very clear at Re = 17.13 (Movie S14 in the supporting information).…”

Mass Transfer Between Recirculation and Main Flow Zones: Is Physically Based Parameterization Possible?

“…Anomalous or non-Fickian transport can be due to the following reasons: (1) the transport has not reached the Fickian regime due to limited spatial and temporal scales (Bouquain et al, 2012;Dentz et al, 2011;Wang et al, 2012;Wang & Cardenas, 2014); (2) pronounced heterogeneity such as manifested in preferential flows (Fiori, 2014;Kang et al, 2016;Molinari et al, 2015;Nordqvist et al, 1996); and, of course, (3) diffusion into and out from slow-flowing, recirculation, or stagnation zones including the rock matrix or grains (Bolster et al, 2014;Boutt et al, 2006;Cardenas et al, 2007;Grisak & Pickens, 1981;Zou et al, 2017). Among these causes for anomalous transport in fractures, the understanding of the role of RZs has remained mostly incomplete and sometimes controversial (Lee et al, 2015).…”