1
|
A Primer on the Dynamical Systems Approach to Transport in Porous Media. Transp Porous Media 2022. [DOI: 10.1007/s11242-022-01811-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
2
|
Trefry MG, Lester DR, Metcalfe G, Wu J. Lagrangian Complexity Persists with Multimodal Flow Forcing in Compressible Porous Systems. Transp Porous Media 2020. [DOI: 10.1007/s11242-020-01487-w] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
3
|
Cho MS, Zhao Z, Thomson NR, Illman WA. Use of steady-state hydraulic tomography to inform the selection of a chaotic advection system. JOURNAL OF CONTAMINANT HYDROLOGY 2020; 229:103559. [PMID: 31784037 DOI: 10.1016/j.jconhyd.2019.103559] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/07/2019] [Revised: 10/08/2019] [Accepted: 10/09/2019] [Indexed: 06/10/2023]
Abstract
The concept of chaotic advection is a novel approach that has the potential to overcome some of the challenges associated with mixing of reagents that commonly occur when injection based in situ treatment techniques are used. The rotated potential mixing (RPM) flow system is one configuration which has been theorized to achieve chaotic advection in porous media, and enhance reagent mixing by periodically re-oriented dipole pumping at a series of radial wells. Prior to field implementation of chaotic advection, the selection of an RPM flow protocol will likely require a numerical model that can adequately represent groundwater flow within the zone of interest. As expected, the hydraulic conductivity (K) field is the most critical input requirement for the selected groundwater flow model. Hydraulic tomography (HT) is an innovative characterization approach that has shown potential to provide information on a K field. In this investigation, we explored whether the same well system required to invoke chaotic advection can also be applied in a HT analysis, and evaluated the use of the generated K tomogram for the selection of RPM flow parameters that can enhance reagent mixing. A series of dipole pumping tests were conducted within an area of interest as defined by the limits of the circular network of eight injection/extraction wells used to invoke chaotic advection. Hydraulic head data collected from independent dipole pumping tests were used in an inverse model to perform steady-state hydraulic tomography (SSHT) analysis to generate a K tomogram. Both the K tomogram and an effective parameter approach (i.e., a single K value assigned across the entire spatial domain as determined by single well pumping and slug tests) produced estimates of hydraulic head that closely resembled those observed due to the relative homogeneous nature of the aquifer and the small spatial scale of the area of interest. In contrast, particle tracking results showed that incorporating a heterogeneous K field significantly enhanced the spatial distribution of particle trajectories indicative of reagent mixing. These findings support the hypothesis that the same well system used to invoke chaotic advection can be combined with SSHT analysis as a viable site characterization tool for delineating the spatial variability of K. Incorporating this K tomogram in a groundwater flow model with a particle tracking engine can be used as a design tool to aid in the selection of a site-specific RPM flow protocol to achieve enhanced reagent mixing.
Collapse
Affiliation(s)
- Michelle S Cho
- Department of Civil and Environmental Engineering, University of Waterloo, 200 University Avenue West, Waterloo, ON N2L 3G1, Canada; Now at Geosyntec Consultants, Inc., Toronto, ON, Canada.
| | - Zhanfeng Zhao
- Department of Earth and Environmental Sciences, University of Waterloo, 200 University Avenue West, Waterloo, ON N2L 3G1, Canada; Now at The Key Laboratory of Water Cycle and Related Land Surface Processes, Institute of Geographic Science and Natural Resources Research, Chinese Academy of Sciences, Beijing, China
| | - Neil R Thomson
- Department of Civil and Environmental Engineering, University of Waterloo, 200 University Avenue West, Waterloo, ON N2L 3G1, Canada
| | - Walter A Illman
- Department of Earth and Environmental Sciences, University of Waterloo, 200 University Avenue West, Waterloo, ON N2L 3G1, Canada
| |
Collapse
|
4
|
|
5
|
|
6
|
Turuban R, Lester DR, Le Borgne T, Méheust Y. Space-Group Symmetries Generate Chaotic Fluid Advection in Crystalline Granular Media. PHYSICAL REVIEW LETTERS 2018; 120:024501. [PMID: 29376725 DOI: 10.1103/physrevlett.120.024501] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/12/2017] [Indexed: 06/07/2023]
Abstract
The classical connection between symmetry breaking and the onset of chaos in dynamical systems harks back to the seminal theory of Noether [Transp. Theory Statist. Phys. 1, 186 (1918)10.1080/00411457108231446]. We study the Lagrangian kinematics of steady 3D Stokes flow through simple cubic and body-centered cubic (bcc) crystalline lattices of close-packed spheres, and uncover an important exception. While breaking of point-group symmetries is a necessary condition for chaotic mixing in both lattices, a further space-group (glide) symmetry of the bcc lattice generates a transition from globally regular to globally chaotic dynamics. This finding provides new insights into chaotic mixing in porous media and has significant implications for understanding the impact of symmetries upon generic dynamical systems.
Collapse
Affiliation(s)
- R Turuban
- Geosciences Rennes, UMR 6118, Université de Rennes 1, CNRS, 35042 Rennes, France
| | - D R Lester
- School of Engineering, RMIT University, 3000 Melbourne, Australia
| | - T Le Borgne
- Geosciences Rennes, UMR 6118, Université de Rennes 1, CNRS, 35042 Rennes, France
| | - Y Méheust
- Geosciences Rennes, UMR 6118, Université de Rennes 1, CNRS, 35042 Rennes, France
| |
Collapse
|
7
|
Koshel KV, Ryzhov EA, Zhmur VV. Effect of the vertical component of diffusion on passive scalar transport in an isolated vortex model. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2015; 92:053021. [PMID: 26651793 DOI: 10.1103/physreve.92.053021] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/04/2015] [Indexed: 06/05/2023]
Abstract
On the basis of the ellipsoidal vortex model and a Monte-Carlo-type diffusion simulation, we examine the flux and ensuing distribution of passive fluid particles through the boundary of an idealized geophysical vortex. Our focus is on features that the horizontal and vertical diffusion components introduce into the fluid particle transport. We examine the concurrent effect of both components, and we compare it with the only horizontal diffusion impact. We analyze the ellipsoid vortex model in two cases: (i) the steady state when the ellipsoid is motionless, i.e., there is no variation in its axes' lengths, and consequently the exterior fluid is not being stirred; (ii) the perturbed case when the ellipsoid rotates periodically, varying it axes' lengths, which results in the appearance of stirred fluid outside the ellipsoid. Influenced by diffusion, a fluid particle is now permitted to move to another vertical horizon, thus there is an increased possibility that the particle will eventually be located in the exterior stirred region rather than in the ellipsoid vortex with the regular dynamics. This is because the area of the horizontal section of the ellipsoid vortex decreases with depth, but the region of stirred exterior fluid extends significantly deeper. Numerical calculations show that factoring in the vertical component of diffusion significantly affects scalar spreading in the horizontal plane in the perturbed case, while in the steady state the vertical component of diffusion only induces dispersion linear growth according to a Gaussian distribution.
Collapse
Affiliation(s)
- Konstantin V Koshel
- V. I. Il'ichev Pacific Oceanological Institute of RAS, 43, Baltiyskaya Street, Vladivostok 690041, Russia and Far Eastern Federal University, 8, Sukhanova Street, Vladivostok 690950, Russia
| | - Evgeny A Ryzhov
- V. I. Il'ichev Pacific Oceanological Institute of RAS, 43, Baltiyskaya Street, Vladivostok 690041, Russia
| | - Vladimir V Zhmur
- P. P. Shirshov Institute of Oceanology of RAS, 36, Nakhimovski prospect, Moscow 117997, Russia and Moscow Institute of Physics and Technology, 9, Institutskiy Pereulok, Dolgoprudnyi, Moscow region 141700, Russia
| |
Collapse
|
8
|
Schlick CP, Isner AB, Umbanhowar PB, Lueptow RM, Ottino JM. On Mixing and Segregation: From Fluids and Maps to Granular Solids and Advection–Diffusion Systems. Ind Eng Chem Res 2015. [DOI: 10.1021/acs.iecr.5b01268] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Conor P. Schlick
- Department of Engineering Sciences
and Applied Mathematics, ‡Department of Chemical
and Biological Engineering, ¶Department of Mechanical Engineering,
and §The Northwestern Institute
on Complex Systems (NICO), Northwestern University, Evanston, Illinois 60208, United States
| | - Austin B. Isner
- Department of Engineering Sciences
and Applied Mathematics, ‡Department of Chemical
and Biological Engineering, ¶Department of Mechanical Engineering,
and §The Northwestern Institute
on Complex Systems (NICO), Northwestern University, Evanston, Illinois 60208, United States
| | - Paul B. Umbanhowar
- Department of Engineering Sciences
and Applied Mathematics, ‡Department of Chemical
and Biological Engineering, ¶Department of Mechanical Engineering,
and §The Northwestern Institute
on Complex Systems (NICO), Northwestern University, Evanston, Illinois 60208, United States
| | - Richard M. Lueptow
- Department of Engineering Sciences
and Applied Mathematics, ‡Department of Chemical
and Biological Engineering, ¶Department of Mechanical Engineering,
and §The Northwestern Institute
on Complex Systems (NICO), Northwestern University, Evanston, Illinois 60208, United States
| | - Julio M. Ottino
- Department of Engineering Sciences
and Applied Mathematics, ‡Department of Chemical
and Biological Engineering, ¶Department of Mechanical Engineering,
and §The Northwestern Institute
on Complex Systems (NICO), Northwestern University, Evanston, Illinois 60208, United States
| |
Collapse
|
9
|
Lester DR, Metcalfe G, Trefry MG. Is chaotic advection inherent to porous media flow? PHYSICAL REVIEW LETTERS 2013; 111:174101. [PMID: 24206492 DOI: 10.1103/physrevlett.111.174101] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/03/2013] [Indexed: 06/02/2023]
Abstract
We show that chaotic advection is inherent to flow through all types of porous media, from granular and packed media to fractured and open networks. The basic topological complexity inherent to all porous media gives rise to chaotic flow dynamics under steady flow conditions, where fluid deformation local to stagnation points imparts a 3D fluid mechanical analog of the baker's map. The ubiquitous nature of chaotic advection has significant implications for the description of transport, mixing, chemical reaction and biological activity in porous media.
Collapse
Affiliation(s)
- D R Lester
- CSIRO Mathematics, Informatics and Statistics, P.O. Box 56, Highett, Victoria 3190, Australia
| | | | | |
Collapse
|
10
|
Trefry MG, Lester DR, Metcalfe G, Ord A, Regenauer-Lieb K. Toward enhanced subsurface intervention methods using chaotic advection. JOURNAL OF CONTAMINANT HYDROLOGY 2012; 127:15-29. [PMID: 21600670 DOI: 10.1016/j.jconhyd.2011.04.006] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/25/2010] [Revised: 04/14/2011] [Accepted: 04/28/2011] [Indexed: 05/30/2023]
Abstract
Many intervention activities in the terrestrial subsurface involve the need to recover/emplace distributions of scalar quantities (e.g. dissolved phase concentrations or heat) from/in volumes of saturated porous media. These scalars can be targeted by pump-and-treat methods or by amendment technologies. Application examples include in-situ leaching for metals, recovery of dissolved contaminant plumes, or utilizing heat energy in geothermal reservoirs. While conventional pumping methods work reasonably well, costs associated with maintaining pumping schedules are high and improvements in efficiency would be welcome. In this paper we discuss how transient switching of the pressure at different wells can intimately control subsurface flow, generating a range of "programmed" flows with various beneficial characteristics. Some programs produce chaotic flows which accelerate mixing, while others create encapsulating flows which can isolate fluid zones for lengthy periods. In a simplified model of an aquifer subject to balanced pumping, chaotic flow topologies have been predicted theoretically and verified experimentally using Hele-Shaw cells. Here, a survey of the key characteristics of chaotic advection is presented. Mathematical methods are used to show how these characteristics may translate into practical situations involving regional flows and heterogeneity. The results are robust to perturbations, and withstand significant aquifer heterogeneity. It is proposed that chaotic advection may form the basis of new efficient technologies for groundwater interventions.
Collapse
|