Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userNatural Groundwater Systems Can Display Chaotic Mixing at the Darcy Scale
151
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Although steady, isotropic Darcy flows are inherently laminar and non-mixing, it is well understood that transient forcing via engineered pumping schemes can induce rapid, chaotic mixing in groundwater. In this study we explore the propensity for such mixing to arise in natural groundwater systems subject to cyclical forcings, e.g. tidal or seasonal influences. Using a conventional linear groundwater flow model subject to tidal forcing, we show that under certain conditions these flows generate Lagrangian transport and mixing phenomena (chaotic advection) near the tidal boundary. We show that aquifer heterogeneity, storativity, and forcing magnitude cause reversals in flow direction over the forcing cycle which, in turn, generate coherent Lagrangian structures and chaos. These features significantly augment fluid mixing and transport, leading to anomalous residence time distributions, flow segregation, and the potential for profoundly altered reaction kinetics. We define the dimensionless parameter groups which govern this phenomenon and explore these groups in connection with a set of well-characterised coastal systems. The potential for Lagrangian chaos to be present near discharge boundaries must be recognized and assessed in field studies.
rate research
Read More
In a recent paper (Trefry et al., 2019) we showed that the interplay of aquifer heterogeneity and poroelasticity can produce complex transport in tidally forced aquifers, with significant implications for solute transport, mixing and reaction. However, what was unknown was how broadly these transport dynamics can arise in natural groundwater systems, and how these dynamics depend upon the aquifer properties, tidal and regional flow characteristics. In this study we answer these questions through parametric studies of these governing properties. We uncover the mechanisms that govern complex transport dynamics and the bifurcations between transport structures with changes in the governing parameters, and we determine the propensity for complex dynamics to occur in natural aquifer systems. These results clearly demonstrate that complex transport structures and dynamics may arise in natural tidally forced aquifers around the world, producing solute transport and mixing behaviour that is very different to that of the conventional Darcy flow picture. Key Points: * Transient Darcy flows can generate complex transport dynamics in heterogeneous compressible aquifers. * This complex transport can trap dispersing solutes for many years. * Global tidal maps indicate widespread potential for complex transport dynamics in coastal zones.
Transport and mixing of scalar quantities in fluid flows is ubiquitous in industry and Nature. Turbulent flows promote efficient transport and mixing by their inherent randomness. Laminar flows lack such a natural mixing mechanism and efficient transport is far more challenging. However, laminar flow is essential to many problems and insight into its transport characteristics of great importance. Laminar transport, arguably, is best described by the Lagrangian fluid motion (`advection) and the geometry, topology and coherence of fluid trajectories. Efficient laminar transport being equivalent to `chaotic advection is a key finding of this approach. The Lagrangian framework enables systematic analysis and design of laminar flows. However, the gap between scientific insights into Lagrangian transport and technological applications is formidable primarily for two reasons. First, many studies concern two-dimensional (2D) flows yet the real world is three dimensional (3D). Second, Lagrangian transport is typically investigated for idealised flows yet practical relevance requires studies on realistic 3D flows. The present review aims to stimulate further development and utilisation of know-how on 3D Lagrangian transport and its dissemination to practice. To this end 3D practical flows are categorised into canonical problems. First, to expose the diversity of Lagrangian transport and create awareness of its broad relevance. Second, to enable knowledge transfer both within and between scientific disciplines. Third, to reconcile practical flows with fundamentals on Lagrangian transport and chaotic advection. This may be a first incentive to structurally integrate the `Lagrangian mindset into the analysis and design of 3D practical flows.
We study numerically joint mixing of salt and colloids by a chaotic velocity field $mathbf{V}$, and how salt inhomogeneities accelerate or delay colloid mixing by inducing a velocity drift $mathbf{V}_{rm dp}$ between colloids and fluid particles as proposed in recent experiments cite{Deseigne2013}. We demonstrate that because the drift velocity is no longer divergence free, small variations to the total velocity field drastically affect the evolution of colloid variance $sigma^2=langle C^2 rangle - langle C rangle^2$. A consequence is that mixing strongly depends on the mutual coherence between colloid and salt concentration fields, the short time evolution of scalar variance being governed by a new variance production term $P=- langle C^2 abla cdot mathbf{V}_{rm dp} rangle/2$ when scalar gradients are not developed yet so that dissipation is weak. Depending on initial conditions, mixing is then delayed or enhanced, and it is possible to find examples for which the two regimes (fast mixing followed by slow mixing) are observed consecutively when the variance source term reverses its sign. This is indeed the case for localized patches modeled as gaussian concentration profiles.
Adjoint-based sensitivity analysis methods are powerful tools for engineers who use flow simulations for design. However, the conventional adjoint method breaks down for scale-resolving simulations like large-eddy simulation (LES) or direct numerical simulation (DNS), which exhibit the chaotic dynamics inherent in turbulent flows. Sensitivity analysis based on least-squares shadowing (LSS) avoids the issues encountered by conventional methods, but has a high computational cost. The following report outlines a new, more computationally efficient formulation of LSS, non-intrusive LSS, and estimates its cost for several canonical flows using Lyapunov analysis.
A new type of instability - electrokinetic instability - and an unusual transition to chaotic motion near a charge-selective surface was studied by numerical integration of the Nernst-Planck-Poisson-Stokes system and a weakly nonlinear analysis near the threshold of instability. Two kinds of initial conditions were considered: (a) white noise initial conditions to mimic room disturbances and subsequent natural evolution of the solution; (b) an artificial monochromatic ion distribution with a fixed wave number to simulate regular wave patterns. The results were studied from the viewpoint of hydrodynamic stability and bifurcation theory. The threshold of electroconvective movement was found by the linear spectral stability theory, the results of which were confirmed by numerical simulation of the entire system. The following regimes, which replace each other as the potential drop between the selective surfaces increases, were obtained: one-dimensional steady solution; two-dimensional steady electroconvective vortices (stationary point in a proper phase space); unsteady vortices aperiodically changing their parameters (homoclinic contour); periodic motion (limit cycle); and chaotic motion. The transition to chaotic motion did not include Hopf bifurcation. Numerical resolution of the thin concentration polarization layer showed spike-like charge profiles along the surface, which could be, depending on the regime, either steady or aperiodically coalescent. The numerical investigation confirmed the experimentally observed absence of regular (near-sinusoidal) oscillations for the overlimiting regimes. There is a qualitative agreement of the experimental and the theoretical values of the threshold of instability, the dominant size of the observed coherent structures, and the experimental and theoretical volt-current characteristics.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
