Passive scalar evolution in peripheral regions

Skipping the Boundary Layer: High-Speed Droplet-Based Immunoassay Using Rayleigh Acoustic Streaming

Acoustic mixing of droplets is a promising way to implement biosensors that combine high speed and minimal reagent consumption. To date, this type of droplet mixing is driven by a volume force resulting from the absorption of high-frequency acoustic waves in the bulk of the fluid. Here, we show that the speed of these sensors is limited by the slow advection of analyte to the sensor surface due to the formation of a hydrodynamic boundary layer. We eliminate this hydrodynamic boundary layer by using much lower ultrasonic frequencies to excite the droplet, which drives a Rayleigh streaming that behaves essentially like a slip velocity. At equal average flow velocity in the droplet, both experiment and three-dimensional simulations show that this provides a three-fold speedup compared to Eckart streaming. Experimentally, we further shorten a SARS-CoV-2 antibody immunoassay from 20 min to 40 s taking advantage of Rayleigh acoustic streaming.

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.

Active mixing of complex fluids at the microscale

Mixing of complex fluids at low Reynolds number is fundamental for a broad range of applications, including materials assembly, microfluidics, and biomedical devices. Of these materials, yield stress fluids (and gels) pose the most significant challenges, especially when they must be mixed in low volumes over short timescales. New scaling relationships between mixer dimensions and operating conditions are derived and experimentally verified to create a framework for designing active microfluidic mixers that can efficiently homogenize a wide range of complex fluids. Active mixing printheads are then designed and implemented for multimaterial 3D printing of viscoelastic inks with programmable control of local composition.

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.

Transport phenomena in chaotic laminar flows

In many important chemical processes, the laminar flow regime is inescapable and defines the performance of reactors, separators, and analytical instruments. In the emerging field of microchemical process or lab-on-a-chip, this constraint is particularly rigid. Here, we review developments in the use of chaotic laminar flows to improve common transport processes in this regime. We focus on four: mixing, interfacial transfer, axial dispersion, and spatial sampling. Our coverage demonstrates the potential for chaos to improve these processes if implemented appropriately. Throughout, we emphasize the usefulness of familiar theoretical models of transport for processes occurring in chaotic flows. Finally, we point out open challenges and opportunities in the field.

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.

Scalar fluctuations from a point source in a turbulent boundary layer

The downstream development of the concentration probability distribution along the mean-plume centerline of a dispersing plume in the wake of a ground-level continuous point source in a neutrally stratified wall-shear layer is studied. It is shown that the concentration distribution is well described by a family of one-parameter gamma distributions, as first suggested by Villermaux and Duplat [Phys. Rev. Lett. 91, 184501 (2003)] in the context of confined mixing. A prediction of the downstream evolution of the parameter k (which specifies the gamma distribution) is obtained. This prediction includes explicitly the effects of mean shear on the mean-square concentration.

Scalar decay in a three-dimensional chaotic flow

The decay of a passive scalar in a three-dimensional chaotic flow is studied using high-resolution numerical simulations. The (volume-preserving) flow considered is a three-dimensional extension of the randomized alternating sine flow employed extensively in studies of mixing in two dimensions. It is used to show that theoretical predictions for two-dimensional flows with small diffusivity carry over to three dimensions even though the stretching properties differ significantly. The variance decay rate, scalar field structure, and time evolution of statistical moments confirm that there are two distinct regimes of scalar decay: a locally controlled regime, which applies when the domain size is comparable to the characteristic length scale of the velocity field, and a globally controlled regime, which applies when the domain is larger. Asymptotic predictions for the variance decay rate in both regimes show excellent agreement with the numerical results. Consideration of both the forward flow and its time reverse makes it possible to compare the scalar evolution in flows with one or two expanding directions; simulations confirm the theoretical prediction that the decay rate of the scalar is the same in both flows, despite the very different scalar field structures.

Scaling laws of peripheral mixing of passive scalar in a wall-shear layer

The scaling laws governing the concentration moments of a passive scalar released from a ground-level localized source in a neutrally stratified wall-shear layer are investigated using a theoretical framework recently formulated by Lebedev and Turitsyn [Phys. Rev. E 69, 036301 (2004)]. For the current application, this theoretical framework is generalized from the smooth random velocity field applicable in the viscous sublayer to the nonsmooth random velocity field that applies to the bulk of the wall-shear layer. Theoretical relationships for the passive scalar concentration moments are compared to a water-channel simulation of turbulent diffusion from a ground-level source in a wall-shear layer. The diffusion measurements in the wall-shear layer are shown to be consistent with the theoretical description and also imply the robustness of the identified scaling laws for the scalar concentration moments.

Streaming by leaky surface acoustic waves

Acoustic streaming, the generation of mean flow by dissipating acoustic waves, provides a promising method for flow pumping in microfluidic devices. In recent years, several groups have been experimenting with acoustic streaming induced by leaky surface waves: (Rayleigh) surface waves excited in a piezoelectric solid interact with a small volume of fluid where they generate acoustic waves and, as result of the viscous dissipation of these waves, a mean flow. We discuss the computation of the corresponding Lagrangian mean flow, which controls the trajectories of fluid particles and hence the mixing properties of the flows generated by this method. The problem is formulated using the averaged vorticity equation which extracts the dominant balance between wave dissipation and mean-flow dissipation. Particular attention is paid to the thin boundary layer that forms at the solid/liquid interface, where the flow is best computed using matched asymptotics. This leads to an explicit expression for a slip velocity, which includes the effect of the oscillations of the boundary. The Lagrangian mean flow is naturally separated into three contributions: an interior-driven Eulerian mean flow, a boundary-driven Eulerian mean flow and the Stokes drift. A scale analysis indicates that the latter two contributions can be neglected in devices much larger than the acoustic wavelength but need to be taken into account in smaller devices. A simple two-dimensional model of mean flow generation by surface acoustic waves is discussed as an illustration.

Scaling laws of passive tracer dispersion in the turbulent surface layer

Experimental results for passive tracer dispersion in the turbulent surface layer under stable conditions are presented. In this case, the dispersion of tracer particles is determined by the interplay of three mechanisms: relative dispersion (celebrated Richardson's mechanism), shear dispersion (particle separation due to variation of the mean velocity field) and specific surface-layer dispersion (induced by the gradient of the energy dissipation rate in the turbulent surface layer). The latter mechanism results in the rather slow (ballistic) law for the mean squared particle separation. Based on a simplified Langevin equation for particle separation we found that the ballistic regime always dominates at large times. This conclusion is supported by our extensive atmospheric observations. Exit-time statistics are derived from the experimental data set and show a reasonable match with the simple dimensional asymptotes for different mechanisms of tracer dispersion, as well as predictions of the multifractal model and experimental data from other sources.

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.

Near-wall turbulent transport of large-Schmidt-number passive scalars

Turbulent diffusion of a passive scalar with a large-Schmidt number (Sc >> 1) is considered in the viscous sublayer of a turbulent channel flow. Close to the wall, the corresponding eddy diffusivity coefficient is expanded as a power series in terms of the viscous distance to the wall y . The coefficients of the series depend on the Schmidt number and the analysis of recent numerical results allows to conclude that in the close vicinity of the wall (y << Sc(-1/3)) , the y(3) term is the dominant term; whereas, at distances relatively large from the wall (Sc(-1/3)<< y <<1) , the y(4) term becomes dominant. Accordingly, in this region the turbulent Schmidt number is not constant but follows a hyperbolic law in terms of the distance to the wall that matches the values taken in the vicinity of the wall, on the order of Sc(-1/3) , with the values of order unity in the rest of the viscous layer. The implications of this behavior on the surface-transfer coefficient are analyzed.

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.

Role of elastic stress in statistical and scaling properties of elastic turbulence

The role of elastic stress in statistical and scaling properties of elastic turbulence in a polymer solution flow between two disks is discussed. The analogy with a small-scale magnetodynamics and a passive scalar turbulent advection in the Batchelor regime is used to explain the experimentally observed statistical properties, the flow structure, and the scaling of elastic turbulence. The emergence of a new length scale, namely, the boundary layer thickness, is observed and studied.

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.

Diffusion of passive scalar in a finite-scale random flow

We consider a solvable model of the decay of scalar variance in a single-scale random velocity field. We show that if there is a separation between the flow scale k(-1 )(flow ) and the box size k(-1 )(box ) , the decay rate lambda proportional, variant ( k(box) / k(flow) )(2) is determined by the turbulent diffusion of the box-scale mode. Exponential decay at the rate lambda is preceded by a transient powerlike decay (the total scalar variance approximately t(-5/2) if the Corrsin invariant is zero, t(-3/2) otherwise) that lasts a time t approximately 1/lambda . Spectra are sharply peaked at k= k(box) . The box-scale peak acts as a slowly decaying source to a secondary peak at the flow scale. The variance spectrum at scales intermediate between the two peaks ( k(box) <<k<< k(flow) ) is approximately k+a k(2) +... (a>0) . The mixing of the flow-scale modes by the random flow produces, for the case of large Péclet number, a k(-1+delta) spectrum at k>> k(flow) , where delta proportional lambda is a small correction. Our solution thus elucidates the spectral make up of the "strange mode," combining small-scale structure and a decay law set by the largest scales.