Experimental observation of batchelor dispersion of passive tracers

Statistical independence of the initial conditions in chaotic mixing

Experimental evidence of the scalar convergence towards a global strange eigenmode independent of the scalar initial condition in chaotic mixing is provided. This convergence, underpinning the independent nature of chaotic mixing in any passive scalar, is presented by scalar fields with different initial conditions casting statistically similar shapes when advected by periodic unsteady flows. As the scalar patterns converge towards a global strange eigenmode, the scalar filaments, locally aligned with the direction of maximum stretching, as described by the Lagrangian stretching theory, stack together in an inhomogeneous pattern at distances smaller than their asymptotic minimum widths. The scalar variance decay becomes then exponential and independent of the scalar diffusivity or initial condition. In this work, mixing is achieved by advecting the scalar using a set of laminar flows with unsteady periodic topology. These flows, that resemble the tendril-whorl map, are obtained by morphing the forcing geometry in an electromagnetic free surface 2D mixing experiment. This forcing generates a velocity field which periodically switches between two concentric hyperbolic and elliptic stagnation points. In agreement with previous literature, the velocity fields obtained produce a chaotic mixer with two regions: a central mixing and an external extensional area. These two regions are interconnected through two pairs of fluid conduits which transfer clean and dyed fluid from the extensional area towards the mixing region and a homogenized mixture from the mixing area towards the extensional region.

Topologies of velocity-field stagnation points generated by a single pair of magnets in free-surface electromagnetic experiments

The velocity fields generated by a static pair of magnets in free-surface electromagnetically forced flows are analyzed for different magnet attitudes, ionic currents, and brine depths. A wide range of laminar velocity fields is obtained despite the forcing simplicity. The velocity fields are classified according to their temporal mean flow topology, which strongly depends on the forcing geometry but barely on its strength, even through the bifurcation to unsteady regimes. The mean flow topology possesses a major influence on the critical Reynolds number Rec under which the steady velocity fields remain stable. The qualitative comparison of the dependence of Rec on the topology is in agreement with previous works. The unsteady configurations evidence the advection of smaller flow structures by the largest scales, commonly known as "sweeping."

Lamination and mixing in three fundamental flow sequences driven by electromagnetic body forces

This article pursues the idea that the degree of striations, called lamination, could be engineered to complement stretching and to design new sequential mixers. It explores lamination and mixing in three new mixing sequences experimentally driven by electromagnetic body forces. To generate these three mixing sequences, Lorentz body forces are dynamically controlled to vary the flow geometry produced by a pair of local jets. The first two sequences are inspired from the "tendril and whorl" and "blinking vortex" flows. The third novel sequence is called the "cat's eyes flip." These three mixing sequences exponentially stretch and laminate material lines representing the interface between two domains to be mixed. Moreover, the mixing coefficient (defined as 1-σ(2)/σ(0)(2) where σ(2)/σ(0)(2) is the rescaled variance) and its rate grow exponentially before saturation. This saturation of the mixing process is related to the departure of the mixing rate from an exponential growth when the striations' thicknesses reach the diffusive length scale of the measurements or species and dyes. Incidentally, in our experiments, for the same energy or forcing input, the cat's eyes flip sequence has higher lamination, stretching, and mixing rates than the tendril and whorl and the blinking vortex sequences. These features show that bakerlike in situ mixers can be conceived by dynamically controlling a pair of local jets and by integrating lamination during stirring stages with persistent geometries. Combined with novel insights provided by the quantification of the lamination, this paper should offer perspectives for the development of new sequential mixers, possibly on all scales.

Moving walls accelerate mixing

Mixing in viscous fluids is challenging, but chaotic advection in principle allows efficient mixing. In the best possible scenario, the decay rate of the concentration profile of a passive scalar should be exponential in time. In practice, several authors have found that the no-slip boundary condition at the walls of a vessel can slow down mixing considerably, turning an exponential decay into a power law. This slowdown affects the whole mixing region, and not just the vicinity of the wall. The reason is that when the chaotic mixing region extends to the wall, a separatrix connects to it. The approach to the wall along that separatrix is polynomial in time and dominates the long-time decay. However, if the walls are moved or rotated, closed orbits appear, separated from the central mixing region by a hyperbolic fixed point with a homoclinic orbit. The long-time approach to the fixed point is exponential, so an overall exponential decay is recovered, albeit with a thin unmixed region near the wall.

Three-dimensional flow in electromagnetically driven shallow two-layer fluids

Recent experiments on a freely evolving dipolar vortex in a homogeneous shallow fluid layer have clearly shown the existence and evolution of complex three-dimensional (3D) flow structures. The present contribution focuses on the 3D structures of a dipolar vortex evolving in a stable shallow two-layer fluid. Experimentally, Stereoscopic Particle Image Velocimetry is used to measure instantaneously all three components of the velocity field in a horizontal plane and 3D numerical simulations provide the full 3D velocity and vorticity fields over the entire flow domain. Remarkably, the experimental results, supported by the numerical simulations, show to a large extent the same 3D structures and evolution as in the single-layer case. The numerical simulations indicate that the so-called frontal circulation in the two-layer fluid is due to deformations of the internal interface. The 3D flow structures will also affect the distribution of massless passive particles released at the free surface. With numerical studies it is shown that these passive particles tend to accumulate or deplete locally where the horizontal velocity field is not divergence-free. This is in contrast with pure two-dimensional incompressible flows where the divergence of the velocity field is zero by definition.

Mechanism to explore lamination rate

A mechanism amenable to laminate and fold flows is identified and quantified. This laminating mechanism follows from a physical and experimental approach relying on the interlaced structure of velocity and Lagrangian acceleration. The Lagrangian acceleration being the resultant of the forces applied on particle fluids, the component of acceleration perpendicular to the velocity vector allows the quantification of a rate of change of the velocity's direction, i.e., the local angular Lagrangian velocity, theta. The spatial variation in theta is then used to predict and measure the lamination and folding rate. To support and illustrate this approach, three basic experimental flows, driven by electromagnetic forces, are discussed and compared. Folding rate intensities are extracted for different characteristic length scales. Also, good agreement is found between grid deformation and the prediction of lamination rate. This quantification of lamination rate opens new avenues for the design of mixers, in particular at low Reynolds numbers.

Slow decay of concentration variance due to no-slip walls in chaotic mixing

Chaotic mixing in a closed vessel is studied experimentally and numerically in different two-dimensional (2D) flow configurations. For a purely hyperbolic phase space, it is well known that concentration fluctuations converge to an eigenmode of the advection-diffusion operator and decay exponentially with time. We illustrate how the unstable manifold of hyperbolic periodic points dominates the resulting persistent pattern. We show for different physical viscous flows that, in the case of a fully chaotic Poincaré section, parabolic periodic points at the walls lead to slower (algebraic) decay. A persistent pattern, the backbone of which is the unstable manifold of parabolic points, can be observed. However, slow stretching at the wall forbids the rapid propagation of stretched filaments throughout the whole domain, and hence delays the formation of an eigenmode until it is no longer experimentally observable. Inspired by the baker's map, we introduce a 1D model with a parabolic point that gives a good account of the slow decay observed in experiments. We derive a universal decay law for such systems parametrized by the rate at which a particle approaches the no-slip wall.

Walls inhibit chaotic mixing

We report on experiments of chaotic mixing in a closed vessel, in which a highly viscous fluid is stirred by a moving rod. We analyze quantitatively how the concentration field of a low-diffusivity dye relaxes towards homogeneity, and we observe a slow algebraic decay of the inhomogeneity, at odds with the exponential decay predicted by most previous studies. Visual observations reveal the dominant role of the vessel wall, which strongly influences the concentration field in the entire domain and causes the anomalous scaling. A simplified 1D model supports our experimental results. Quantitative analysis of the concentration pattern leads to scalings for the distributions and the variance of the concentration field consistent with experimental and numerical results.

Pair dispersion and doubling time statistics in two-dimensional turbulence

Experimental measurements of pair separation statistics in two-dimensional turbulence are reported for an electromagnetically forced stratified-layer system with simultaneous ranges of direct-enstrophy and inverse-energy transfer separated by a well-defined spatial injection scale. Data for pair separation as a function of time are analyzed to determine the dependence of separation statistics in both regimes. Using doubling-time statistics, we show how the measured scalings of the mean quantities are consistent with exponential behavior in the enstrophy range and power-law behavior in the inverse-energy range. Exponential scaling of the doubling-time probability distribution function agrees well with theoretical predictions. Finite size effects are shown to play an important role in the interpretation of the data.

Chaotic mixing in a steady flow in a microchannel

We report experiments on mixing of a passively advected fluorescent dye in a low Reynolds number flow in a microscopic channel. The channel is a chain of repeating segments with a custom designed profile that generates a steady three-dimensional flow with stretching and folding, and chaotic mixing. A few statistical characteristics of mixing in the flow are studied and are all found to agree with theoretical and experimental results for the flows in the Batchelor regime of mixing that are chaotic in time. The proposed microchannel provides fast and efficient mixing and is simple to fabricate.

Batchelor scaling in fast-flowing soap films

The dynamics of a passive scalar such as a dye in the far dissipative range of fluid turbulence is a central problem in nonlinear physics. An important prediction for this problem was made by Batchelor over 40 years ago and is known as Batchelor's scaling law. We here present strong evidence in favor of this law for the thickness fluctuations in the flow of a soap film past a flat plate. The results also capture the dissipative range of the scalar which turns out to have universal features. The probability density function of the scalar increments and their structure functions come out in nice agreement with theoretical predictions.

Simple ideas on mixing and fragmentation

The mechanisms building the overall concentration distribution in a scalar mixture, and the drops in a spray, are examined successively. In both cases, the distributions belong to a unique family of distributions stable by self-convolution, the signature of the aggregation process from which they originate.

Accurate discretization of advection-diffusion equations

We present an exact mathematical transformation which converts a wide class of advection-diffusion equations into a form allowing simple and direct spatial discretization in all dimensions, and thus the construction of accurate and more efficient numerical algorithms. These discretized forms can also be viewed as master equations which provide an alternative mesoscopic interpretation of advection-diffusion processes in terms of diffusion with spatially varying hopping rates.

Mixing by polymers: experimental test of decay regime of mixing

By using high molecular weight fluorescent passive tracers with different diffusion coefficients and by changing the fluid velocity we study the dependence of a characteristic mixing length on the Peclet number, Pe, which controls the mixing efficiency. The mixing length is found to be related to Pe by a power law, L(mix) proportional, variant Pe0.26+/-0.01, and increases faster than expected for an unbounded chaotic flow. The role of the boundaries in the mixing length abnormal growth is clarified. The experimental findings are in good quantitative agreement with recent theoretical predictions.

Mixing as an aggregation process

Experiments show how a stirred scalar mixture relaxes towards uniformity through an aggregation process. The elementary bricks are stretched sheets whose rates of diffusive smoothing and coalescence build up the overall mixture concentration distribution. The cases studied, in particular, include mixtures in two and three dimensions, with different stirring protocols and Reynolds numbers which all lead to a unique family of concentration distributions stable by self-convolution, the signature of the aggregation mechanism from which they originate.

Boundary effects on chaotic advection-diffusion chemical reactions

A theory of a fast binary chemical reaction, A+B-->C, in a statistically stationary bounded chaotic flow at large Peclet number Pe and large Damköhler number Da is described. The first stage correspondent to formation of the developed lamellar structure in the bulk part of the flow is terminated by an exponential decay, proportional, variant exp((-lambdat) (where lambda is the Lyapunov exponent of the flow), of the chemicals in the bulk. The second and the third stages are due to the chemicals remaining in the boundary region. During the second stage, the amounts of A and B decay proportional, variant 1/sqrt[t], whereas the decay law during the third stage is exponential, proportional, variant exp((-gammat), where gamma approximately lambda/sqrt[Pe].

Decay of scalar turbulence revisited

We demonstrate that at long times the rate of passive scalar decay in a turbulent, or simply chaotic, flow is dominated by regions where mixing is less efficient. We examine two situations. The first is of a spatially homogeneous stationary turbulent flow with both viscous and inertial scales present. It is shown that at large times scalar fluctuations decay algebraically in time at all spatial scales. The second example explains chaotic stationary flow in a disk/pipe. The boundary region of the flow controls the long-time decay, which is algebraic at some transient times, but becomes exponential, with the decay rate dependent on the scalar diffusion coefficient, at longer times.

Decay of passive scalars under the action of single scale smooth velocity fields in bounded two-dimensional domains: from non-self-similar probability distribution functions to self-similar eigenmodes

We examine the decay of passive scalars with small, but nonzero, diffusivity in bounded two-dimensional (2D) domains. The velocity fields responsible for advection are smooth (i.e., they have bounded gradients) and of a single large scale. Moreover, the scale of the velocity field is taken to be similar to the size of the entire domain. The importance of the initial scale of variation of the scalar field with respect to that of the velocity field is strongly emphasized. If these scales are comparable and the velocity field is time periodic, we see the formation of a periodic scalar eigenmode. The eigenmode is numerically realized by means of a deterministic 2D map on a lattice. Analytical justification for the eigenmode is available from theorems in the dynamo literature. Weakening the notion of an eigenmode to mean statistical stationarity, we provide numerical evidence that the eigenmode solution also holds for aperiodic flows (represented by random maps). Turning to the evolution of an initially small scale scalar field, we demonstrate the transition from an evolving (i.e., non-self-similar) probability distribution function (pdf) to a stationary (self-similar) pdf as the scale of variation of the scalar field progresses from being small to being comparable to that of the velocity field (and of the domain). Furthermore, the non-self-similar regime itself consists of two stages. Both stages are examined and the coupling between diffusion and the distribution of the finite time Lyapunov exponents is shown to be responsible for the pdf evolution.

Polymers in 2D turbulence: suppression of large scale fluctuations

Small quantities of a long chain molecule or polymer affect two-dimensional turbulence in unexpected ways. Their presence inhibits the transfers of energy to large scales causing their suppression in the energy density spectrum. This also leads to the change of the spectral properties of a passive scalar which turns out to be highly sensitive to the presence of energy transfers.

Small-scale structure of nonlinearly interacting species advected by chaotic flows

We study the spatial patterns formed by interacting biological populations or reacting chemicals under the influence of chaotic flows. Multiple species and nonlinear interactions are explicitly considered, as well as cases of smooth and nonsmooth forcing sources. The small-scale structure can be obtained in terms of characteristic Lyapunov exponents of the flow and of the chemical dynamics. Different kinds of morphological transitions are identified. Numerical results from a three-component plankton dynamics model support the theory, and they serve also to illustrate the influence of asymmetric couplings. (c) 2002 American Institute of Physics.

Thickness fluctuations in turbulent soap films

Rapidly flowing soap films provide a simple and attractive system to study two-dimensional hydrodynamics and turbulence. By measuring the rapid fluctuations of the thickness of the film in the turbulent regime, we find that the statistics of these fluctuations closely resemble those of a passive scalar field in a turbulent flow. The scalar spectra are well described by Kolmogorov-like scaling while the high-order moments show clear deviations from regular scaling just like dye or temperature fluctuations in 3D turbulent flows.

Scalings of scalar structure functions in a velocity field with coherent vortical structures

In planar turbulence modeled as an isotropic and homogeneous collection of two-dimensional noninteracting compact vortices, the structure functions S(p)(r) of a statistically stationary passive scalar field have the following scaling behavior in the limit where the Péclet number Pe-->infinity: S(p)(r) approximately const+ln(r/L Pe(-1/3)) for L Pe(-1/3)<<r<<L, S(p)(r) approximately (r/L Pe(-1/3))(6(1-D)) for L Pe(-1/2)<<r<<L Pe(-1/3), where L is a large scale and D is the fractal codimension of the spiral scalar structures generated by the vortices (1/2 < or =D < 2/3). Note that L Pe(-1/2) is the scalar Taylor microscale that stems naturally from our analytical treatment of the advection-diffusion equation. The essential ingredients of our theory are the locality of interscale transfer and Lundgren's time average assumption. A phenomenological theory explicitly based only on these two ingredients reproduces our results and a generalization of this phenomenology to spatially smooth chaotic flows yields (k ln k)(-1) generalized power spectra for the advected scalar fields.

Intermittency of a passive tracer in the inverse energy cascade

We report an experimental study of the dispersion of a passive tracer in the two-dimensional inverse energy cascade, which shows that a nonintermittent velocity field can sustain a strongly intermittent concentration field. The experiment suggests the exponents of the intermittent concentration field saturate at large orders towards xi(infinity) approximately 1.2. These observations are in excellent agreement with a recent numerical work [A. Celani, A. Lanotte, A. Mazzino, and M. Vergassola, Phys. Rev. Lett. 84, 2385 (2000)] and theoretical expectations [E. Balkovsky and V. Lebedev, Phys. Rev. E 58, 5776 (1998); V. Yakhot, ibid. 55, 329 (1997)].

Critical geometry of two-dimensional passive scalar turbulence

Passive scalars advected by a magnetically driven two-dimensional turbulent flow are analyzed using methods of statistical topography. The passive tracer concentration is interpreted as the height of a random surface, and the scaling properties of its contour loops are analyzed. Various exponents that describe the loop ensemble are measured and compared to a scaling theory. This leads to a geometrical criterion for the intermittency of scalar fluctuations.