Topological Kinematics of Mixing Jean-Luc Thiffeault Matthew Finn Emmanuelle Gouillart http://www.ma.imperial.ac.uk/ jeanluc Department of Mathematics Imperial College London Topological Kinematics of Mixing p.1/29
Experiment of Boyland et al. [P. L. Boyland, H. Aref, and M. A. Stremler, J. Fluid Mech. 403, 277 (2000)] [Movies by Matt Finn: boyland1 boyland2 ] Topological Kinematics of Mixing p.2/29
s 1 Generators of the n-braid Group s i σ i i 1 i i + 1 i + 2 σ 1 i i 1 i i + 1 i + 2 A generator σ i, i = 1,..., n 1 is the clockwise interchange of the i th and (i + 1)th rod. The generators obey the presentation σ i+1 σ i σ i+1 = σ i σ i+1 σ i σ i σ j = σ j σ i, i j > 1 These generators are used to characterise the motion of the rods. Topological Kinematics of Mixing p.3/29
The Two BAS Stirring Protocols σ 1 σ 2 protocol σ 1 1 σ 2 protocol [P. L. Boyland, H. Aref, and M. A. Stremler, J. Fluid Mech. 403, 277 (2000)] Topological Kinematics of Mixing p.4/29
Rod Trajectories as Braids σ 1 σ 2 protocol σ 1 1 σ 2 protocol [P. L. Boyland, H. Aref, and M. A. Stremler, J. Fluid Mech. 403, 277 (2000)] Topological Kinematics of Mixing p.5/29
Matrix Representation of the σ i We use the convenient Burau representation 1 τ 0 σ i = I i 2 0 τ 0 I n i 2 0 1 1 where τ C. The matrices are (n 1) (n 1). The Burau matrices satisfy the presentation (of course), but for n > 4 they do not provide a faithful representation. This is of no great consequence for our purposes. Topological Kinematics of Mixing p.6/29
Topological Entropy of a Braid Practically speaking, the topological entropy of a braid is a lower bound on the line-stretching exponent of the flow! This is Eminently Reasonable TM : I II σ 2 = II The Burau representation has an awesome property: if Σ is the Burau representation of a braid word, I topological entropy of braid sup spr Σ τ =1 [Kolev] where spr denotes the spectral radius (the magnitude of the largest eigenvalue). Topological Kinematics of Mixing p.7/29
The Difference between BAS s Two Protocols The matrices associated the generators have eigenvalues on the unit-circle. The first (bad) stirring protocol has eigenvalues on the unit circle The second (good) protocol has largest eigenvalue (3 + 5)/2 = 2.6180. So for the second protocol the length of a line joining the rods grows exponentially! That is, material lines have to stretch by at least a factor of 2.6180 each time we execute the protocol σ 1 1 σ 2. This is guaranteed to hold in some neighbourhood of the rods (Thurston Nielsen theorem). Topological Kinematics of Mixing p.8/29
Freely-moving Rods in a Cavity Flow [A. Vikhansky, Physics of Fluids 15, 1830 (2003)] Topological Kinematics of Mixing p.9/29
Particle Orbits are Topological Obstacles Choose any fluid particle orbit (green dot). Material lines must bend around the orbit: it acts just like a rod! The idea: pick any n fluid particles and follow them. How do they braid around each other? Topological Kinematics of Mixing p.10/29
Detecting Braiding Events σ i σ 1 i σ 1 i ements PSfrag replacements σ i In the second case there is no net braid: the two elements cancel each other. Topological Kinematics of Mixing p.11/29
Random Sequence of Braids We end up with a sequence of braids, with matrix representation Σ (N) = σ (N) σ (2) σ (1) where σ (µ) {σ i, σi 1 } and N is the number of braiding events detected after a time t. The largest eigenvalue of Σ (N) is a measure of the complexity of the braiding motion, called the braiding factor. Random matrix theory says that the braiding factor can grow exponentially! We call the rate of exponential growth the braiding Lyapunov exponent or just braiding exponent. Topological Kinematics of Mixing p.12/29
Non-braiding Motion ents First consider the motion of of three points in concentric circles with irrationally-related frequencies. y 1 0.5 0 0.5 1 1 0 1 x Eigenvalue of Σ (N) 20 15 10 5 0 0 100 200 300 400 500 The braiding factor grows linearly, which means that the braiding exponent is zero. Notice that the eigenvalue often returns to unity. t Topological Kinematics of Mixing p.13/29
Blinking-vortex Flow To demonstrate good braiding, we need a chaotic flow on a bounded domain (a spatially-periodic flow won t do). Aref s blinking-vortex flow is ideal. Active Vortex Inactive Vortex First half of period Second half of period The only parameter is the circulation Γ of the vortices. Topological Kinematics of Mixing p.14/29
Blinking Vortex: Non-braiding Motion ents For Γ = 0.5, the blinking vortex has only small chaotic regions. y 1 0.5 0 0.5 1 1 0 1 x Eigenvalue of Σ (N) 80 60 40 20 0 0 50 100 150 200 One of the orbits is chaotic, the other two are closed. t Topological Kinematics of Mixing p.15/29
Blinking Vortex: Braiding Motion ents For Γ = 13, the blinking vortex is globally chaotic. y 1 0.5 0 0.5 1 1 0 1 x Eigenvalue of Σ (N) 10 20 10 10 10 0 0 50 100 150 200 The braiding factor now grows exponentially. In the same time interval as for Γ = 0.5, the final value is now of order 10 20 rather than 80! t Topological Kinematics of Mixing p.16/29
Averaging over many Triplets 400 Log of eigenvalue of Σ (N) replacements 350 300 250 200 150 100 50 slope = 0.187 Γ = 13 0 0 500 1000 1500 2000 Averaged over 100 random triplets. t Topological Kinematics of Mixing p.17/29
Comparison with Lyapunov Exponents Lyapunov exponent g replacements 2.6 2.4 2.2 2 1.8 1.6 1.4 Line stretching exponent Lyapunov exponent 0.15 0.2 0.25 Braiding Lyapunov exponent Γ varies from 8 to 20. Topological Kinematics of Mixing p.18/29
Beyond Three Particles 1 0.8 even odd Braiding exponent eplacements 0.6 0.4 0.2 Γ = 16.5 0 0 5 10 15 20 n Topological Kinematics of Mixing p.19/29
But does it Saturate? Braiding exponent eplacements 1.5 1 0.5 even odd Γ = 16.5 Well, it really should... 0 0 20 40 60 n Topological Kinematics of Mixing p.20/29
One Rod Mixer: The Kenwood Chef Topological Kinematics of Mixing p.21/29
Poincaré Section Topological Kinematics of Mixing p.22/29
Stretching of Lines Topological Kinematics of Mixing p.23/29
Motion of Islands Make a braid from the motion of the rod and the periodic islands. Most (74%) of the topological entropy is accounted for by this braid. Topological Kinematics of Mixing p.24/29
Motion of Islands and Unstable Periodic Orbits Now we also include unstable periods orbits as well as the stable ones (islands). Almost all (99%) of the topological entropy is accounted for by this braid. Topological Kinematics of Mixing p.25/29
Periodic Orbits as Rods [movie: sf_periodic.avi] Topological Kinematics of Mixing p.26/29
Blowup of the Tip 0.52 0.51 0.5 0.49 0.48 0.47 0.06 0.065 0.07 0.075 0.08 0.085 0.09 0.095 0.1 Topological Kinematics of Mixing p.27/29
Curvature of the Tip 12000 10000 8000 slope = 786.6 curvature 6000 4000 2000 0 0 5 10 15 periods [Preliminary result: not sure this is generic.] Topological Kinematics of Mixing p.28/29
Conclusions Topological chaos involves moving obstacles in a 2D flow, which create nontrivial braids. The complexity of a braid can be represented by the largest eigenvalue of a product of matrices the braiding factor. Any collection of n particles can potentially braid. The complexity of the braid is a good measure of chaos. No need for infinitesimal separation of trajectories or derivatives of the velocity field. Many issues to investigate: faithfulness of representation, improved lower bound for topological entropy, reducibility of braids, size of ghost rod for periodic orbits... See J-LT, Phys. Rev. Lett. 94, 084502 (2005). Topological Kinematics of Mixing p.29/29