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

Fluid-Solid Reaction in Porous Media as a Chaotic Restart Process

Chemical and biological reactions at fluid-solid interfaces are central to a broad range of porous material applications and research. Pore-scale solute transport limitations can reduce reaction rates, with marked consequences for a wide spectrum of natural and engineered processes. Recent advances show that chaotic mixing occurs spontaneously in porous media, but its impact on surface reactions is unknown. We show that pore-scale chaotic mixing significantly increases reaction efficiency compared to nonchaotic flows. We find that reaction rates are well described in terms of diffusive first-passage times of reactants to the solid interface subjected to a stochastic restart process resulting from Lagrangian chaos. Under chaotic mixing, the shear layer at no-slip interfaces sets the restart rate and leads to a characteristic scaling of reaction efficiency with Péclet number, in excellent agreement with numerical simulations. Reaction rates are insensitive to the flow topology as long as flow is chaotic, suggesting the relevance of this process to a broad range of porous materials.

Geometric mixing

Mixing fluids often involves a periodic action, like stirring one's tea. But reciprocating motions in fluids at low Reynolds number, in Stokes flows where inertia is negligible, lead to periodic cycles of mixing and unmixing, because the physics, molecular diffusion excepted, is time reversible. So how can fluid be mixed in such circumstances? The answer involves a geometric phase. Geometric phases are found everywhere in physics as anholonomies, where after a closed circuit in the parameters, some system variables do not return to their original values. We discuss the geometric phase in fluid mixing: geometric mixing. This article is part of the theme issue 'Stokes at 200 (part 2)'.

Efficiency of laminar and turbulent mixing in wall-bounded flows

A turbulent flow mixes in general more rapidly a passive scalar than a laminar flow does. From an energetic point of view, for statistically homogeneous or periodic flows, the laminar regime is more efficient. However, the presence of walls may change this picture. We consider in this investigation mixing in two-dimensional laminar and turbulent wall-bounded flows using direct numerical simulation. We show that for sufficiently large Schmidt number, turbulent flows more efficiently mix a wall-bounded scalar field than a chaotic or laminar flow does. The mixing efficiency is shown to be a function of the Péclet number, and a phenomenological explanation yields a scaling law, consistent with the observations.

Deceleration of one-dimensional mixing by discontinuous mappings

We present a computational study of a simple one-dimensional map with dynamics composed of stretching, permutations of equally sized cells, and diffusion. We observe that the combination of the aforementioned dynamics results in eigenmodes with long-time exponential decay rates. The decay rate of the eigenmodes is shown to be dependent on the choice of permutation and changes nonmonotonically with the diffusion coefficient for many of the permutations. The global mixing rate of the map M in the limit of vanishing diffusivity approximates well the decay rates of the eigenmodes for small diffusivity, however this global mixing rate does not bound the rates for all values of the diffusion coefficient. This counterintuitively predicts a deceleration in the asymptotic mixing rate with an increasing diffusivity rate. The implications of the results on finite time mixing are discussed.

Short-time behavior of advecting-diffusing scalar fields in Stokes flows

This article addresses the short-term decay of advecting-diffusing scalar fields in Stokes flows. The analysis is developed in two main subparts. In the first part, we present an analytic approach for a class of simple flow systems expressed mathematically by the one-dimensional advection-diffusion equation w(y)∂(ξ)φ=ε∂(y)(2)φ+iV(y)φ-ε'φ, where ξ is either time or axial coordinate and iV(y) an imaginary potential. This class of systems encompasses both open- and closed-flow models and corresponds to the dynamics of a single Fourier mode in parallel flows. We derive an analytic expression for the short-time (short-length) decay of φ, and show that this decay is characterized by a universal behavior that depends solely on the singularity of the ratio of the transverse-to-axial velocity components V(eff)(y)=V(y)/w(y), corresponding to the effective potential in the imaginary potential formulation. If V(eff)(y) is smooth, then ||φ||(L(2))(ξ)=exp(-ε'ξ-bξ(3)), where b>0 is a constant. Conversely, if the effective potential is singular, then ||φ||(L(2))(ξ)=1-aξ(ν) with a>0. The exponent ν attains the value 5/3 at the very early stages of the process, while for intermediate stages its value is 3/5. By summing over all of the Fourier modes, a stretched exponential decay is obtained in the presence of nonimpulsive initial conditions, while impulsive conditions give rise to an early-stage power-law behavior. In the second part, we consider generic, chaotic, and nonchaotic autonomous Stokes flows, providing a kinematic interpretation of the results found in the first part. The kinematic approach grounded on the warped-time transformation complements the analytical theory developed in the first part.

Rate of chaotic mixing and boundary behavior

We discuss rigorous results on the rate of mixing for an idealized model of a class of fluid mixing device. These show that the decay of correlations of a scalar field is governed by the presence of boundaries in the domain, and in particular by the behavior of the modeled fluid at such boundaries.

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.

Rotation shields chaotic mixing regions from no-slip walls

We report on the decay of a passive scalar in chaotic mixing protocols where the wall of the vessel is rotated, or a net drift of fluid elements near the wall is induced at each period. As a result the fluid domain is divided into a central isolated chaotic region and a peripheral regular region. Scalar patterns obtained in experiments and simulations converge to a strange eigenmode and follow an exponential decay. This contrasts with previous experiments [Gouillart, Phys. Rev. Lett. 99, 114501 (2007)] with a chaotic region spanning the whole domain, where fixed walls constrained mixing to follow a slower algebraic decay. Using a linear analysis of the flow close to the wall, as well as numerical simulations of Lagrangian trajectories, we study the influence of the rotation velocity of the wall on the size of the chaotic region, the approach to its bounding separatrix, and the decay rate of the scalar.