BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Jean-Luc Thiffeault

Jean-Luc Thiffeault appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
physics.flu-dynnlin.CDcond-mat.softmath.DSmath.GTBiological Physicsmath.APphysics.geo-phchao-dyncond-mat.stat-mechmath-phmath.HOmath.MPmath.PRphysics.pop-ph

Catalog footprint

What is connected

35works
15topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

35 published item(s)

preprint2022arXiv

Anisotropic active Brownian particle with a fluctuating propulsion force

The active Brownian particle (ABP) model describes a swimmer, synthetic or living, whose direction of swimming is a Brownian motion. The swimming is due to a propulsion force, and the fluctuations are typically thermal in origin. We present a 2D model where the fluctuations arise from nonthermal noise in a propelling force acting at a single point, such as that due to a flagellum. We take the overdamped limit and find several modifications to the traditional ABP model. Since the fluctuating force causes a fluctuating torque, the diffusion tensor describing the process has a coupling between translational and rotational degrees of freedom. An anisotropic particle also exhibits a mass-dependent noise-induced drift, which does not disappear in the overdamped limit. We show that these effects have measurable consequences for the long-time diffusivity of active particles, in particular adding a contribution that is independent of where the force acts.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Open-flow mixing and transfer operators

We study finite-time mixing in time-periodic open flow systems. We describe the transport of densities in terms of a transfer operator, which is represented by the transition matrix of a finite-state Markov chain. The transport processes in the open system are organized by the chaotic saddle and its stable and unstable manifolds. We extract these structures directly from leading eigenvectors of the transition matrix. We use different measures to quantify the degree of mixing and show that they give consistent results in parameter studies of two model systems.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Reaching for the surface: Spheroidal microswimmers in surface gravity waves

Microswimmers (planktonic microorganisms or artificial active particles) immersed in a fluid interact with the ambient flow, altering their trajectories. In surface gravity waves, a common goal for microswimmers is vertical migration (e.g., to reach the free surface or to dive to deeper depths). By modelling microswimmers as spheroidal bodies with an intrinsic swimming velocity that supplements advection and reorientation by the flow, we investigate how shape and swimming affect vertical transport of microswimmers in waves. We find that it is possible for microswimmers to be initially swimming downwards, but to recover and head back to the surface, and vice versa. This is because the coupling between swimming and flow-induced reorientations introduces a shape dependency in the vertical transport. From a wave-averaged analysis of microswimmer trajectories, we show that each trajectory is bounded by critical planes in the position-orientation phase space that depend only on the shape. We also give explicit solutions to these trajectories and determine the fraction of microswimmers that begin within the water column and eventually reach the surface. For microswimmers that are initially randomly oriented, the fraction that reach the surface increases monotonically as the starting depth decreases, as expected, but also varies with shape and swimming speed. In the limit of small swimming speed, the fraction of highly prolate microswimmers reaching the surface is 0.5, suggesting that these swimmers would be able to choose direction of vertical transport with small changes in swimming behaviour.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The mathematics of burger flipping

What is the most effective way to grill food? Timing is everything, since only one surface is exposed to heat at a given time. Should we flip only once, or many times? We present a simple model of cooking by flipping, and some interesting observations emerge. The rate of cooking depends on the spectrum of a linear operator, and on the fixed point of a map. If the system has symmetric thermal properties, the rate of cooking becomes independent of the sequence of flips, as long as the last point to be cooked is the midpoint. After numerical optimization, the flipping intervals become roughly equal in duration as their number is increased, though the final interval is significantly longer. We find that the optimal improvement in cooking time, given an arbitrary number of flips, is about 29% over a single flip. This toy problem has some characteristics reminiscent of turbulent thermal convection, such as a uniform average interior temperature with boundary layers.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Phase separation in the advective Cahn-Hilliard equation

The Cahn--Hilliard equation is a classic model of phase separation in binary mixtures that exhibits spontaneous coarsening of the phases. We study the Cahn--Hilliard equation with an imposed advection term in order to model the stirring and eventual mixing of the phases. The main result is that if the imposed advection is sufficiently mixing then no phase separation occurs, and the solution instead converges exponentially to a homogeneous mixed state. The mixing effectiveness of the imposed drift is quantified in terms of the dissipation time of the associated advection-hyperdiffusion equation, and we produce examples of velocity fields with a small dissipation time. We also study the relationship between this quantity and the dissipation time of the standard advection-diffusion equation.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Using braids to quantify interface growth and coherence in a rotor-oscillator flow

The growth rate of material interfaces is an important proxy for mixing and reaction rates in fluid dynamics, and can also be used to identify regions of coherence. Estimating such growth rates can be difficult, since they depend on detailed properties of the velocity field, such as its derivatives, that are hard to measure directly. When an experiment gives only sparse trajectory data, it is natural to encode planar trajectories as mathematical braids, which are topological objects that contain information on the mixing characteristics of the flow, in particular through their action on topological loops. We test such braid methods on an experimental system, the rotor-oscillator flow, which is well-described by a theoretical model. We conduct a series of laboratory experiments to collect particle tracking and particle image velocimetry data, and use the particle tracks to identify regions of coherence within the flow that match the results obtained from the model velocity field. We then use the data to estimate growth rates of material interface, using both the braid approach and numerical simulations. The interface growth rates follow similar qualitative trends in both the experiment and model, but have significant quantitative differences, suggesting that the two are not as similar as first seems. Our results shows that there are challenges in using the braid approach to analyze data, in particular the need for long trajectories, but that these are not insurmountable.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Winding of a Brownian particle around a point vortex

We derive the asymptotic winding law of a Brownian particle in the plane subjected to a tangential drift due to a point vortex. For winding around a point, the normalized winding angle converges to an inverse Gamma distribution. For winding around a disk, the angle converges to a distribution given by an elliptic theta function. For winding in an annulus, the winding angle is asymptotically Gaussian with a linear drift term. We validate our results with numerical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2017arXiv

Lyapunov exponents for the random product of two shears

We give lower and upper bounds on both the Lyapunov exponent and generalised Lyapunov exponents for the random product of positive and negative shear matrices. These types of random products arise in applications such as fluid stirring devices. The bounds, obtained by considering invariant cones in tangent space, give excellent accuracy compared to standard and general bounds, and are increasingly accurate with increasing shear. Bounds on generalised exponents are useful for testing numerical methods, since these exponents are difficult to compute in practice.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Chaotic Geodesics

When a shallow layer of inviscid fluid flows over a substrate, the fluid particle trajectories are, to leading order in the layer thickness, geodesics on the two-dimensional curved space of the substrate. Since the two-dimensional geodesic equation is a two degree-of-freedom autonomous Hamiltonian system, it can exhibit chaos, depending on the shape of the substrate. We find chaotic behaviour for a range of substrates.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Finite-Time Braiding Exponents

Topological entropy is a common measure of the rate of mixing in a flow. It can be computed by partition methods, or by estimating the growth rate of material lines or other material elements. This requires detailed knowledge of the velocity field, which is not always available, such as when we only know a few particle trajectories (ocean float data, for example). We propose an alternative approximation to topological entropy,applicable to two-dimensional flows, which uses only a finite number of trajectories as input data. To represent these sparse data sets, we use braids, algebraic objects that record how strands, i.e., trajectories, exchange positions with respect to a projection axis. Material curves advected by the flow are represented as simplified loop coordinates. The exponential rate at which a braid deforms loops over a finite time interval as the strands exchange places is the Finite-Time Braiding Exponent (FTBE) and serves as a proxy for topological entropy of the two-dimensional flow. We demonstrate that FTBEs are robust with respect to the value of numerical time step,details of braid representation, and choice of initial conditions inside the mixing region. We also explore how closely the FTBEs approximate topological entropy depending on the number and length of trajectories used.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Scalar Decay in Chaotic Mixing

I review the local theory of mixing, which focuses on infinitesimal blobs of scalar being advected and stretched by a random velocity field. An advantage of this theory is that it provides elegant analytical results. A disadvantage is that it is highly idealised. Nevertheless, it provides insight into the mechanism of chaotic mixing and the effect of random fluctuations on the rate of decay of the concentration field of a passive scalar.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Estimating topological entropy from the motion of stirring rods

Stirring a two-dimensional viscous fluid with rods is often an effective way to mix. The topological features of periodic rod motions give a lower bound on the topological entropy of the induced flow map, since material lines must `catch' on the rods. But how good is this lower bound? We present examples from numerical simulations and speculate on what affects the 'gap' between the lower bound and the measured topological entropy. The key is the sign of the rod motion's action on first homology of the orientation double cover of the punctured disk.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Topological entropy and secondary folding

A convenient measure of a map or flow's chaotic action is the topological entropy. In many cases, the entropy has a homological origin: it is forced by the topology of the space. For example, in simple toral maps, the topological entropy is exactly equal to the growth induced by the map on the fundamental group of the torus. However, in many situations the numerically-computed topological entropy is greater than the bound implied by this action. We associate this gap between the bound and the true entropy with 'secondary folding': material lines undergo folding which is not homologically forced. We examine this phenomenon both for physical rod-stirring devices and toral linked twist maps, and show rigorously that for the latter secondary folds occur.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Trajectory entanglement in dense granular materials

The particle-scale dynamics of granular materials have commonly been characterized by the self-diffusion coefficient $D$. However, this measure discards the collective and topological information known to be an important characteristic of particle trajectories in dense systems. Direct measurement of the entanglement of particle space-time trajectories can be obtained via the topological braid entropy $\Sbraid$, which has previously been used to quantify mixing efficiency in fluid systems. Here, we investigate the utility of $\Sbraid$ in characterizing the dynamics of a dense, driven granular material at packing densities near the static jamming point $ϕ_J$. From particle trajectories measured within a two-dimensional granular material, we typically observe that $\Sbraid$ is well-defined and extensive. However, for systems where $ϕ\gtrsim 0.79$, we find that $\Sbraid$ (like $D$) is not well-defined, signifying that these systems are not ergodic on the experimental timescale. Both $\Sbraid$ and $D$ decrease with either increasing packing density or confining pressure, independent of the applied boundary condition. The related braiding factor provides a means to identify multi-particle phenomena such as collective rearrangements. We discuss possible uses for this measure in characterizing granular systems.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A low-Reynolds-number treadmilling swimmer near a semi-infinite wall

We investigate the behavior of a treadmilling microswimmer in a two-dimensional unbounded domain with a semi-infinite no-slip wall. The wall can also be regarded as a probe or pipette inserted into the flow. We solve the governing evolution equations in an analytical form and numerically calculate trajectories of the swimmer for several different initial positions and orientations. We then compute the probability that the treadmilling organism can escape the vicinity of the wall. We find that many trajectories in a 'wedge' around the wall are likely to escape. This suggests that inserting a probe or pipette in a suspension of organism may push away treadmilling swimmers.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Detecting coherent structures using braids

The detection of coherent structures is an important problem in fluid dynamics, particularly in geophysical applications. For instance, knowledge of how regions of fluid are isolated from each other allows prediction of the ultimate fate of oil spills. Existing methods detect Lagrangian coherent structures, which are barriers to transport, by examining the stretching field as given by finite-time Lyapunov exponents. These methods are very effective when the velocity field is well-determined, but in many applications only a small number of flow trajectories are known, for example when dealing with oceanic float data. We introduce a topological method for detecting invariant regions based on a small set of trajectories. In the method we regard the two-dimensional trajectory data as a braid in three dimensions, with time being the third coordinate. Invariant regions then correspond to trajectories that travel together and do not entangle other trajectories. We detect these regions by examining the growth of hypothetical loops surrounding sets of trajectories, and searching for loops that show negligible growth.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Stirring up trouble: Multi-scale mixing measures for steady scalar sources

The mixing efficiency of a flow advecting a passive scalar sustained by steady sources and sinks is naturally defined in terms of the suppression of bulk scalar variance in the presence of stirring, relative to the variance in the absence of stirring. These variances can be weighted at various spatial scales, leading to a family of multi-scale mixing measures and efficiencies. We derive a priori estimates on these efficiencies from the advection--diffusion partial differential equation, focusing on a broad class of statistically homogeneous and isotropic incompressible flows. The analysis produces bounds on the mixing efficiencies in terms of the Peclet number, a measure the strength of the stirring relative to molecular diffusion. We show by example that the estimates are sharp for particular source, sink and flow combinations. In general the high-Peclet number behavior of the bounds (scaling exponents as well as prefactors) depends on the structure and smoothness properties of, and length scales in, the scalar source and sink distribution. The fundamental model of the stirring of a monochromatic source/sink combination by the random sine flow is investigated in detail via direct numerical simulation and analysis. The large-scale mixing efficiency follows the upper bound scaling (within a logarithm) at high Peclet number but the intermediate and small-scale efficiencies are qualitatively less than optimal. The Peclet number scaling exponents of the efficiencies observed in the simulations are deduced theoretically from the asymptotic solution of an internal layer problem arising in a quasi-static model.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Using multiscale norms to quantify mixing and transport

Mixing is relevant to many areas of science and engineering, including the pharmaceutical and food industries, oceanography, atmospheric sciences, and civil engineering. In all these situations one goal is to quantify and often then to improve the degree of homogenisation of a substance being stirred, referred to as a passive scalar or tracer. A classical measure of mixing is the variance of the concentration of the scalar, which can be related to the $L^2$ norm of the concentration field. Recently other norms have been used to quantify mixing, in particular the mix-norm as well as negative Sobolev norms. These norms have the advantage that unlike variance they decay even in the absence of diffusion, and their decay corresponds to the flow being mixing in the sense of ergodic theory. General Sobolev norms weigh scalar gradients differently, and are known as multiscale norms for mixing. We review the applications of such norms to mixing and transport, and show how they can be used to optimise the stirring and mixing of a decaying passive scalar. We then review recent work on the less-studied case of a continuously-replenished scalar field --- the source-sink problem. In that case the flows that optimally reduce the norms are associated with transport rather than mixing: they push sources onto sinks, and vice versa.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Hydraulic jumps on an incline

When a fluid jet strikes an inclined solid surface at normal incidence, gravity creates a flow pattern with a thick outer rim resembling a parabola and reminiscent of a hydraulic jump. There appears to be little theory or experiments describing simple aspects of this phenomenon, such as the maximum rise height of the fluid above the impact point, and its dependence on jet velocity and inclination angle. We address this with experiments, and present a simple theory based on horizontal hydraulic jumps which accounts for the rise height and its scaling, though without describing the shape of the parabolic envelope.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Measures of mixing quality in open flows with chaotic advection

We address the evaluation of mixing efficiency in experiments of chaotic mixing inside an open-flow channel. Since the open flow continuously brings new fluid into the limited mixing region, it is difficult to define relevant mixing indices, as fluid particles experience typically very different stretching and mixing histories. The repeated stretching and folding of a spot of dye leads to a persistent pattern. We propose that the normalized standard deviation of this characteristic pattern is a good measure of the mixing quality of the flow. We discuss the link between this measure and mixing of continuously-injected dye, and investigate it using an idealized map.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

On the minimum dilatation of pseudo-Anosov homeomorphisms on surfaces of small genus

We find the minimum dilatation of pseudo-Anosov homeomorphisms that stabilize an orientable foliation on surfaces of genus three, four, or five, and provide a lower bound for genus six to eight. Our technique also simplifies Cho and Ham's proof of the least dilatation of pseudo-Anosov homeomorphisms on a genus two surface. For genus g=2 to 5, the mimimum dilatation is the smallest Salem number for polynomials of degree 2g.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Optimal stirring strategies for passive scalar mixing

We address the challenge of optimal incompressible stirring to mix an initially inhomogeneous distribution of passive tracers. As a quantitative measure of mixing we adopt the $H^{-1}$ norm of the scalar fluctuation field, equivalent to the (square-root of the) variance of a low-pass filtered image of the tracer concentration field. First we establish that this is a useful gauge even in the absence of molecular diffusion: its vanishing as $t --> \infty$ is evidence of the stirring flow's mixing properties in the sense of ergodic theory. Then we derive absolute limits on the total amount of mixing, as a function of time, on a periodic spatial domain with a prescribed instantaneous stirring energy or stirring power budget. We subsequently determine the flow field that instantaneously maximizes the decay of this mixing measure---when such a flow exists. When no such `steepest descent' flow exists (a possible but non-generic situation) we determine the flow that maximizes the growth rate of the $H^{-1}$ norm's decay rate. This local-in-time optimal stirring strategy is implemented numerically on a benchmark problem and compared to an optimal control approach using a restricted set of flows. Some significant challenges for analysis are outlined.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Stirring by squirmers

We analyse a simple 'Stokesian squirmer' model for the enhanced mixing due to swimming micro-organisms. The model is based on a calculation of Thiffeault & Childress [Physics Letters A, 374, 3487 (2010), arXiv:0911.5511], where fluid particle displacements due to inviscid swimmers are added to produce an effective diffusivity. Here we show that for the viscous case the swimmers cannot be assumed to swim an infinite distance, even though their total mass displacement is finite. Instead, the largest contributions to particle displacement, and hence to mixing, arise from random changes of direction of swimming and are dominated by the far-field stresslet term in our simple model. We validate the results by numerical simulation. We also calculate nonzero Reynolds number corrections to the effective diffusivity. Finally, we show that displacements due to randomly-swimming squirmers exhibit PDFs with exponential tails and a short-time superdiffusive regime, as found previously by several authors. In our case, the exponential tails are due to 'sticking' near the stagnation points on the squirmer's surface.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Stirring by swimming bodies

We consider the stirring of an inviscid fluid caused by the locomotion of bodies through it. The swimmers are approximated by non-interacting cylinders or spheres moving steadily along straight lines. We find the displacement of fluid particles caused by the nearby passage of a swimmer as a function of an impact parameter. We use this to compute the effective diffusion coefficient from the random walk of a fluid particle under the influence of a distribution of swimming bodies. We compare with the results of simulations. For typical sizes, densities and swimming velocities of schools of krill, the effective diffusivity in this model is five times the thermal diffusivity. However, we estimate that viscosity increases this value by two orders of magnitude.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Mixing efficiency of open flows

Mixing in open incompressible flows is studied in a model problem with inhomogeneous passive scalar injection on an inlet boundary. As a measure of the efficiency of stirring, the bulk scalar concentration variance is bounded and the bound is shown to be sharp at low Peclet number. Although no specific flow saturating the bound at high Peclet number is produced here, the estimate is conjectured to be approached for flows possessing sufficiently sustained chaotic regions.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Topological optimisation of rod-stirring devices

There are many industrial situations where rods are used to stir a fluid, or where rods repeatedly stretch a material such as bread dough or taffy. The goal in these applications is to stretch either material lines (in a fluid) or the material itself (for dough or taffy) as rapidly as possible. The growth rate of material lines is conveniently given by the topological entropy of the rod motion. We discuss the problem of optimising such rod devices from a topological viewpoint. We express rod motions in terms of generators of the braid group, and assign a cost based on the minimum number of generators needed to write the braid. We show that for one cost function -- the topological entropy per generator -- the optimal growth rate is the logarithm of the golden ratio. For a more realistic cost function,involving the topological entropy per operation where rods are allowed to move together, the optimal growth rate is the logarithm of the silver ratio, $1+\sqrt{2}$. We show how to construct devices that realise this optimal growth, which we call silver mixers.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Braids of entangled particle trajectories

In many applications, the two-dimensional trajectories of fluid particles are available, but little is known about the underlying flow. Oceanic floats are a clear example. To extract quantitative information from such data, one can measure single-particle dispersion coefficients, but this only uses one trajectory at a time, so much of the information on relative motion is lost. In some circumstances the trajectories happen to remain close long enough to measure finite-time Lyapunov exponents, but this is rare. We propose to use tools from braid theory and the topology of surface mappings to approximate the topological entropy of the underlying flow. The procedure uses all the trajectory data and is inherently global. The topological entropy is a measure of the entanglement of the trajectories, and converges to zero if they are not entangled in a complex manner (for instance, if the trajectories are all in a large vortex). We illustrate the techniques on some simple dynamical systems and on float data from the Labrador sea.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Nonlinear dynamics of phase separation in thin films

We present a long-wavelength approximation to the Navier-Stokes Cahn-Hilliard equations to describe phase separation in thin films. The equations we derive underscore the coupled behaviour of free-surface variations and phase separation. We introduce a repulsive substrate-film interaction potential and analyse the resulting fourth-order equations by constructing a Lyapunov functional, which, combined with the regularizing repulsive potential, gives rise to a positive lower bound for the free-surface height. The value of this lower bound depends on the parameters of the problem, a result which we compare with numerical simulations. While the theoretical lower bound is an obstacle to the rupture of a film that initially is everywhere of finite height, it is not sufficiently sharp to represent accurately the parametric dependence of the observed dips or `valleys' in free-surface height. We observe these valleys across zones where the concentration of the binary mixture changes sharply, indicating the formation of bubbles. Finally, we carry out numerical simulations without the repulsive interaction, and find that the film ruptures in finite time, while the gradient of the Cahn--Hilliard concentration develops a singularity.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Open-flow mixing: Experimental evidence for strange eigenmodes

We investigate experimentally the mixing dynamics in a channel flow with a finite stirring region undergoing chaotic advection. We study the homogenization of dye in two variants of an eggbeater stirring protocol that differ in the extent of their mixing region. In the first case, the mixing region is separated from the side walls of the channel, while in the second it extends to the walls. For the first case, we observe the onset of a permanent concentration pattern that repeats over time with decaying intensity. A quantitative analysis of the concentration field of dye confirms the convergence to a self-similar pattern, akin to the strange eigenmodes previously observed in closed flows. We model this phenomenon using an idealized map, where an analysis of the mixing dynamics explains the convergence to an eigenmode. In contrast, for the second case the presence of no-slip walls and separation points on the frontier of the mixing region leads to non-self-similar mixing dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

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 et al., 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.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

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 2-D 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 Poincare 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 1-D 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.

0 citations0 reviewsNo code linkedOpen access
preprint1998arXiv

Invariants and Labels in Lie-Poisson Systems

Reduction is a process that uses symmetry to lower the order of a Hamiltonian system. The new variables in the reduced picture are often not canonical: there are no clear variables representing positions and momenta, and the Poisson bracket obtained is not of the canonical type. Specifically, we give two examples that give rise to brackets of the noncanonical Lie-Poisson form: the rigid body and the two-dimensional ideal fluid. From these simple cases, we then use the semidirect product extension of algebras to describe more complex physical systems. The Casimir invariants in these systems are examined, and some are shown to be linked to the recovery of information about the configuration of the system. We discuss a case in which the extension is not a semidirect product, namely compressible reduced MHD, and find for this case that the Casimir invariants lend partial information about the configuration of the system.

0 citations0 reviewsNo code linkedOpen access