Mixing and transport processes are extremely important topics with regard to the water quality of environmental flows. Examples include transport and dispersion of contaminants, distribution of mine tailings in tailings ponds, growth and decay of micro-organisms in reservoirs or lakes, motions of suspended sediment in rivers, lakes and estuaries, distribution of heavy metal particles near outfalls, and so on. Recent research has shown that mixing in non-turbulent flows can be greatly enhanced by complicated particle behaviour caused by chaotic advection. Chaotic advection, firstly introduced by Hassan Aref in early 1980s, is a state-of-the-art concept derived from nonlinear dynamics and is widely used as an approach to investigate transport and mixing problems in fluid flows.
Consider a passive particle tracer in a two-dimensional flow field
,
which
are the so-called advection equations and provide a Lagrangian description of
the fluid motion.
The above two-dimensional advection equations are integrable if the
flow is steady, while for unsteady flow they may be non-integrable.
Integrable solutions of the advection
equations lead to regular advection, and the non-integrable cases are
characterized by chaotic advection
which
can be described as particle motion sensitive to initial conditions, i.e.,
initially nearby trajectories diverge at an exponential rate. To produce
chaotic advection, the Eulerian velocity field in the advection equations
is
not necessarily very complicated
|
BV flow: Poincaré sections formed by recording particle positions at the end of each period, indicating that particle motions change from regular to chaotic as the governing parameter μ increases its magnitude. |
|
BV flow: advection of a particle line by the flow with μ = 0.01 (characterised by regular particle motions). |
|
BV flow: advection of a particle line by the flow with μ = 0.5. |
|
BV flow: advection of a small particle patch originally located in the middle of the two vortices by the flow with μ = 0.3. |
|
BV flow: advection of a small particle patch originally located in the middle of the two vortices by the flow with μ = 0.5. |
|
BV flow: advection of a small particle patch originally located in the middle of the two vortices by the flow with μ = 2.0. |
Since it was introduced in early 1980s, chaotic advection has been widely used by researchers in diverse areas of fluid mechanics to explain a variety of experimental or real world problems. However the chaotic advection problems considered in most of the published work are restricted to use idealised analytical, numerical solutions or experimental data as input flow fields. Few investigations have considered mixing and transport processes involving chaotic advection in natural shallow water bodies.
We investigate wind-induced passive mixing and transport in shallow water basins. A shallow water basin, like a lake or reservoir, can often be defined as an enclosed system, in which wind shear is normally the major source of momentum while other factors such as inflows and outflows can only locally modify the flow pattern. Kranenburg (Kranenburg 1992, J of Hydraulic Research 30) studied the chaotic advection caused by oscillatory wind-induced circulation in a circular lake with the intention of applying the concept of chaotic advection to environmental hydraulics. Using a simplified analytical model for the wind-induced hydrodynamics, Kranenburg demonstrated that inert particle motions became chaotic when the surface wind stress periodically and abruptly changed direction even though its magnitude was constant.
|
Wind-induced advection in a circular shallow basin: Poincaré sections corresponding to particle tracers released along the x-axis for different μ, where μ represents dimensionless storm duration (two storm durations form one period). |
|
(a) t = 0 T (b) t = 3 T (c) t = 10 T (d) t = 16 T Wind-induced advection in a circular shallow basin: advection of particle line by the flow with μ = 0.28. |
|
(a) t = 2 T (b) t = 4 T (c) t = 7 T (d) t = 16 T Wind-induced advection in a circular shallow basin: advection of particle line by the flow with μ = 0.84. |
Kranenburg’s method is not general since an idealised analytical velocity field is not available for natural shallow water flows, where the domain geometries and flow patterns are usually very complicated. Herein we examine the chaotic mixing of particles in shallow flows using numerical methods that are applicable to natural cases. The wind-induced flow field is numerically predicted using an adaptive quadtree grid based shallow flow solver that can be refined according to the boundary geometry and local flow features. Wind-induced mixing and transport of passive re-suspended material in a natural tailings pond is then considered. And the wind field is assumed to be sinusoidally changing , to provide a more realistic representation for natural cases.
|
(a) (b)
Mine tailings
pond:
|
|
Mine tailings pond: Poincaré sections corresponding to 10 particle tracers for different μ (dimensionless storm duration). |
|
Mine tailings pond: advection of a small particle patch (indicated by the tiny black square at t = 0) by the flow with μ = 1.125. |
In certain cases, the particles themselves impose dynamical effects in addition to being advected by the background fluid flows. Examples include the growth and decay of algae in lakes, and chemical reactions due to mine tailing re-suspension. The dynamical behavour of active particles is significantly different from that of passive particle tracers in chaotic hydrodynamical flows. Herein, the advected active particles are assumed to undergo certain chemical or biological changes but do not affect the background fluid flow, i.e., the advected substances are chemically or biologically active but hydrodynamically passive. This topic has attracted considerable attentions, and the derived theories have been applied to a wide range of chemical and ecological problems.
|
(a)
(b) (c)
Chaotic
advection of reactive particles (A
+ B → 2B)
in
von Kármán vortex street: (a) flow field predicted by the quadtree based shallow
water equation solver, visualised in terms of vorticity; (b) active particle (species
A
|