No Arabic abstract
We continue our study of chaotic mixing and transport of passive particles in a simple model of a meandering jet flow [Prants, et al, Chaos {bf 16}, 033117 (2006)]. In the present paper we study and explain phenomenologically a connection between dynamical, topological, and statistical properties of chaotic mixing and transport in the model flow in terms of dynamical traps, singular zones in the phase space where particles may spend arbitrary long but finite time [Zaslavsky, Phys. D {bf 168--169}, 292 (2002)]. The transport of passive particles is described in terms of lengths and durations of zonal flights which are events between two successive changes of sign of zonal velocity. Some peculiarities of the respective probability density functions for short flights are proven to be caused by the so-called rotational-islands traps connected with the boundaries of resonant islands (including those of the vortex cores) filled with the particles moving in the same frame. Whereas, the statistics of long flights can be explained by the influence of the so-called ballistic-islands traps filled with the particles moving from a frame to frame.
Cross-jet transport of passive scalars in a kinematic model of the meandering laminar two-dimensional incompressible flow which is known to produce chaotic mixing is studied. We develop a method for detecting barriers to cross-jet transport in the phase space which is a physical space for our model. Using tools from theory of nontwist maps, we construct a central invariant curve and compute its characteristics that may serve good indicators of the existence of a central transport barrier, its strength, and topology. Computing fractal dimension, length, and winding number of that curve in the parameter space, we study in detail change of its geometry and its destruction that are caused by local bifurcations and a global bifurcation known as reconnection of separatrices of resonances. Scenarios of reconnection are different for odd and even resonances. The central invariant curves with rational and irrational (noble) values of winding numbers are arranged into hierarchical series which are described in terms of continued fractions. Destruction of central transport barrier is illustrated for two ways in the parameter space: when moving along resonant bifurcation curves with rational values of the winding number and along curves with noble (irrational) values.
This work reviews the present position of and surveys future perspectives in the physics of chaotic advection: the field that emerged three decades ago at the intersection of fluid mechanics and nonlinear dynamics, which encompasses a range of applications with length scales ranging from micrometers to hundreds of kilometers, including systems as diverse as mixing and thermal processing of viscous fluids, microfluidics, biological flows, and oceanographic and atmospheric flows.
Dynamical and statistical properties of tracer advection are studied in a family of flows produced by three point-vortices of different signs. A collapse of all three vortices to a single point is then possible. Tracer dynamics is analyzed by numerical construction of Poincar{e} sections, and is found to be strongly chaotic: advection pattern in the region around the center of vorticity is dominated by a well developed stochastic sea, which grows as the vortex system approaches the collapse; at the same time, the islands of regular motion around vortices, known as vortex cores, shrink. An estimation of the cores radii from the minimum distance of vortex approach to each other is obtained. Tracer transport was found to be anomalous: for all of the three numerically investigated cases, the variance of the tracer distribution grows faster than a linear function of time, corresponding to a super-diffusive regime. The transport exponent varies with time decades, implying the presence of multi-fractal transport features. Yet, its value is never too far from 3/2, indicating some kind of universality. Statistics of Poincar{e} recurrences is non-Poissonian: distributions have long power-law tails. The anomalous properties of tracer statistics are the result of the complex structure of the advection phase space, in particular, of strong stickiness on the boundaries between the regions of chaotic and regular motion. The role of the different phase space structures involved in this phenomenon is analyzed. Based on this analysis, a kinetic description is constructed, which takes into account different time and space scalings by using a fractional equation.
We study the effects of finite-sizeness on small, neutrally buoyant, spherical particles advected by open chaotic flows. We show that, when projected onto configuration space, the advected finite-size particles disperse about the unstable manifold of the chaotic saddle that governs the passive advection. Using a discrete-time system for the dynamics, we obtain an expression predicting the dispersion of the finite-size particles in terms of their Stokes parameter at the onset of the finite-sizeness induced dispersion. We test our theory in a system derived from a flow and find remarkable agreement between our expression and the numerically measured dispersion.
We discuss a modification to Random Matrix Theory eigenstate statistics, that systematically takes into account the non-universal short-time behavior of chaotic systems. The method avoids diagonalization of the Hamiltonian, instead requiring only a knowledge of short-time dynamics for a chaotic system or ensemble of similar systems. Standard Random Matrix Theory and semiclassical predictions are recovered in the limits of zero Ehrenfest time and infinite Heisenberg time, respectively. As examples, we discuss wave function autocorrelations and cross-correlations, and show that significant improvement in accuracy is obtained for simple chaotic systems where comparison can be made with brute-force diagonalization. The accuracy of the method persists even when the short-time dynamics of the system or ensemble is known only in a classical approximation. Further improvement in the rate of convergence is obtained when the method is combined with the correlation function bootstrapping approach introduced previously.