Finite-size effects on active chaotic advection

Fractal structures in the chaotic advection of passive scalars in leaky planar hydrodynamical flows

The advection of passive scalars in time-independent two-dimensional incompressible fluid flows is an integrable Hamiltonian system. It becomes non-integrable if the corresponding stream function depends explicitly on time, allowing the possibility of chaotic advection of particles. We consider for a specific model (double gyre flow), a given number of exits through which advected particles can leak, without disturbing the flow itself. We investigate fractal escape basins in this problem and characterize fractality by computing the uncertainty exponent and basin entropy. Furthermore, we observe the presence of basin boundaries with points exhibiting the Wada property, i.e., boundary points that separate three or more escape basins.

Levitation of heavy particles against gravity in asymptotically downward flows

In the fluid transport of particles, it is generally expected that heavy particles carried by a laminar fluid flow moving downward will also move downward. We establish a theory to show, however, that particles can be dynamically levitated and lifted by interacting vortices in such flows, thereby moving against gravity and the asymptotic direction of the flow, even when they are orders of magnitude denser than the fluid. The particle levitation is rigorously demonstrated for potential flows and supported by simulations for viscous flows. We suggest that this counterintuitive effect has potential implications for the air-transport of water droplets and the lifting of sediments in water.

Inhomogeneous distribution of water droplets in cloud turbulence

We consider sedimentation of small particles in the turbulent flow where fluid accelerations are much smaller than acceleration of gravity g. The particles are dragged by the flow by linear friction force. We demonstrate that the pair-correlation function of particles' concentration diverges with decreasing separation as a power law with negative exponent. This manifests fractal distribution of particles in space. We find that the exponent is proportional to ratio of integral of energy spectrum of turbulence times the wave number over g. The proportionality coefficient is a universal number independent of particle size. We derive the spectrum of Lyapunov exponents that describes the evolution of small patches of particles. It is demonstrated that particles separate dominantly in the horizontal plane. This provides a theory for the recently observed vertical columns formed by the particles. We confirm the predictions by direct numerical simulations of Navier-Stokes turbulence. The predictions include conditions that hold for water droplets in warm clouds thus providing a tool for the prediction of rain formation.

Influence of the history force on inertial particle advection: gravitational effects and horizontal diffusion

We analyze the effect of the Basset history force on the sedimentation or rising of inertial particles in a two-dimensional convection flow. When memory effects are neglected, the system exhibits rich dynamics, including periodic, quasiperiodic, and chaotic attractors. Here we show that when the full advection dynamics is considered, including the history force, both the nature and the number of attractors change, and a fractalization of their basins of attraction appears. In particular, we show that the history force significantly weakens the horizontal diffusion and changes the speed of sedimentation or rising. The influence of the history force is dependent on the size of the advected particles, being stronger for larger particles.

Dynamics of light particles in oscillating cellular flows

The dynamics of light particles in chaotic oscillating cellular flows is investigated both analytically and numerically by means of Monte Carlo simulations. At level of linear analysis (in the oscillation amplitude) we determined how the known fixed points relative to the stationary cellular flow deform into closed stable trajectories. Once the latter have been analytically determined, we numerically show that they possess the dynamical role of attracting all asymptotic trajectories for a wide range of parameters values. The robustness of the attracting trajectories is tested by adding a white-noise contribution to the particle equation of motion. As a result, attracting trajectories persist up to a critical Péclet number above which an average rising velocity sets in. Possible implications of our results on the long-standing problem related to the explanation of the observed oceanic plankton patchiness will be also discussed.

Coagulation and fragmentation dynamics of inertial particles

Inertial particles suspended in many natural and industrial flows undergo coagulation upon collisions and fragmentation if their size becomes too large or if they experience large shear. Here we study this coagulation-fragmentation process in time-periodic incompressible flows. We find that this process approaches an asymptotic dynamical steady state where the average number of particles of each size is roughly constant. We compare the steady-state size distributions corresponding to two fragmentation mechanisms and for different flows and find that the steady state is mostly independent of the coagulation process. While collision rates determine the transient behavior, fragmentation determines the steady state. For example, for fragmentation due to shear, flows that have very different local particle concentrations can result in similar particle size distributions if the temporal or spatial variation in shear forces is similar.

Transport and diffusion in the embedding map

We study the transport properties of passive inertial particles in two-dimensional (2D) incompressible flows. Here, the particle dynamics is represented by the four-dimensional dissipative embedding map of the 2D area-preserving standard map which models the incompressible flow. The system is a model for impurity dynamics in a fluid and is characterized by two parameters, the inertia parameter alpha and the dissipation parameter gamma . The aerosol regime, where the particles are denser than the fluid, and the bubble regime, where they are less dense than the fluid, correspond to the parameter regimes alpha>1 and alpha<1 , respectively. Earlier studies of this system show a rich phase diagram with dynamical regimes corresponding to periodic orbits, chaotic structures, and mixed regimes. We obtain the statistical characterizers of transport for this system in these dynamical regimes. These are the recurrence time statistics, the diffusion exponent, and the distribution of jump lengths. The recurrence time distribution shows a power-law tail in the dynamical regimes, where there is preferential concentration of particles in sticky regions of the phase space, and an exponential decay in mixing regimes. The diffusion exponent shows behavior of three types-normal, subdiffusive, and superdiffusive, depending on the parameter regimes. Phase diagrams of the system are constructed to differentiate different types of diffusion behavior, as well as the behavior of the absolute drift. We correlate the dynamical regimes seen for the system at different parameter values with the transport properties observed at these regimes and in the behavior of the transients. This system also shows the existence of a crisis and unstable dimension variability at certain parameter values. The signature of the unstable dimension variability is seen in the statistical characterizers of transport. We discuss the implications of our results for realistic systems.

Aggregation and fragmentation dynamics of inertial particles in chaotic flows

Inertial particles advected in chaotic flows often accumulate in strange attractors. While moving in these fractal sets they usually approach each other and collide. Here we consider inertial particles aggregating upon collision. The new particles formed in this process are larger and follow the equation of motion with a new parameter. These particles can in turn fragment when they reach a certain size or shear forces become sufficiently large. The resulting system consists of a large set of coexisting dynamical systems with a varying number of particles. We find that the combination of aggregation and fragmentation leads to an asymptotic steady state. The asymptotic particle size distribution depends on the mechanism of fragmentation. The size distributions resulting from this model are consistent with those found in raindrop statistics and in stirring tank experiments.

Moving finite-size particles in a flow: a physical example of pitchfork bifurcations of tori

The motion of small, spherical particles of finite size in fluid flows at low Reynolds numbers is described by the strongly nonlinear Maxey-Riley equations. Due to the Stokes drag, the particle motion is dissipative, giving rise to the possibility of attractors in phase space. We investigate the case of an infinite cellular flow field with time-periodic forcing. The dynamics of this system are studied in a part of the parameter space. We focus particularly on the size of the particles, whose variations are most important in active physical processes, for example, for aggregation and fragmentation of particles. Depending on their size the particles will settle on different attractors in phase space in the long-term limit, corresponding to periodic, quasiperiodic, or chaotic motion. One of the invariant sets that can be observed in a large part of this parameter region is a quasiperiodic motion in the form of a torus. We identify some of the bifurcations that these tori undergo, as particle size and mass ratio relative to the fluid are varied. In this way we provide a physical example for sub- and supercritical pitchfork bifurcations of tori.

Coexistence of inertial competitors in chaotic flows

We investigate the dynamics of inertial particles immersed in open chaotic flows. We consider the generic problem of competition between different species, e.g., phytoplankton populations in oceans. The strong influence from inertial effects is shown to result in the persistence of different species even in cases when the passively advected species cannot coexist. Multispecies coexistence in the ocean can be explained by the fact that the unstable manifold is different for each advected competitor of different size.

Finite-size effects on open chaotic advection

We study the effects of finite-sizeness on small, neutrally buoyant, spherical particles advected by open chaotic flows. We show that, when observed in the configuration or physical 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-size 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.

Stability of attractors formed by inertial particles in open chaotic flows

Particles having finite mass and size advected in open chaotic flows can form attractors behind structures. Depending on the system parameters, the attractors can be chaotic or nonchaotic. But, how robust are these attractors? In particular, will small, random perturbations destroy the attractors? Here, we address this question by utilizing a prototype flow system: a cylinder in a two-dimensional incompressible flow, behind which the von Kármán vortex street forms. We find that attractors formed by inertial particles behind the cylinder are fragile in that they can be destroyed by small, additive noise. However, the resulting chaotic transient can be superpersistent in the sense that its lifetime obeys an exponential-like scaling law with the noise amplitude, where the exponent in the exponential dependence can be large for small noise. This happens regardless of the nature of the original attractor, chaotic or nonchaotic. We present numerical evidence and a theory to explain this phenomenon. Our finding makes direct experimental observation of superpersistent chaotic transients feasible and it also has implications for problems of current concern such as the transport and trapping of chemically or biologically active particles in large-scale flows.

Continuum description of finite-size particles advected by external flows: the effect of collisions

The equation of the density field of an assembly of macroscopic particles advected by a hydro-dynamic flow is derived from the microscopic description of the system. This equation allows one to recognize the role and the relative importance of the different microscopic processes implicit in the model: the driving of the external flow, the inertia of the particles, and the collisions among them. The validity of the density description is confirmed by comparisons of numerical studies of the continuum equation with direct simulation Monte Carlo simulations of hard disks advected by a chaotic flow. We show that the collisions have two competing roles: a dispersinglike effect and a clustering effect (even for elastic collisions). An unexpected feature is also observed in the system: the presence of collisions can reverse the effect of inertia, so that grains with lower inertia are more clusterized.

Universality in active chaos

Many examples of chemical and biological processes take place in large-scale environmental flows. Such flows generate filamental patterns which are often fractal due to the presence of chaos in the underlying advection dynamics. In such processes, hydrodynamical stirring strongly couples into the reactivity of the advected species and might thus make the traditional treatment of the problem through partial differential equations difficult. Here we present a simple approach for the activity in inhomogeneously stirred flows. We show that the fractal patterns serving as skeletons and catalysts lead to a rate equation with a universal form that is independent of the flow, of the particle properties, and of the details of the active process. One aspect of the universality of our approach is that it also applies to reactions among particles of finite size (so-called inertial particles).

Reactive dynamics of inertial particles in nonhyperbolic chaotic flows

Anomalous kinetics of infective (e.g., autocatalytic) reactions in open, nonhyperbolic chaotic flows are important for many applications in biological, chemical, and environmental sciences. We present a scaling theory for the singular enhancement of the production caused by the universal, underlying fractal patterns. The key dynamical invariant quantities are the effective fractal dimension and effective escape rate, which are primarily determined by the hyperbolic components of the underlying dynamical invariant sets. The theory is general as it includes all previously studied hyperbolic reactive dynamics as a special case. We introduce a class of dissipative embedding maps for numerical verification.

Sand stirred by chaotic advection

We study the spatial structure of a granular material, N particles subject to inelastic mutual collisions, when it is stirred by a bidimensional smooth chaotic flow. A simple dynamical model is introduced where four different time scales are explicitly considered: (i) the Stokes time, accounting for the inertia of the particles, (ii) the mean collision time among the grains, (iii) the typical time scale of the flow, and (iv) the inverse of the Lyapunov exponent of the chaotic flow, which gives a typical time for the separation of two initially close parcels of fluid. Depending on the relative values of these different times, a complex scenario appears for the long-time steady spatial distribution of particles, where clusters of particles may or may not appear.

Advection of finite-size particles in open flows

It is known that small, spherical particles with inertia do not follow the local velocity field of the flow. Here we investigate the motion of such particles and particle ensembles immersed in open, unsteady flows which, in the case of ideal pointlike tracers, generate chaotic Lagrangian trajectories. Due to the extra force terms in the equations of motion (such as Stokes drag, added mass) the inertial tracer trajectories become described by a high-dimensional (2d+1, with d being the flow's dimension) chaotic dynamics, which can drastically differ from the (d+1)-dimensional ideal tracer dynamics. As a consequence, we find parameter regimes (in terms of density and size), where long-term tracer trapping can occur for the inertial particle, even for flows in which no ideal, pointlike passive tracers can be trapped. These studies are performed in a model of a two-dimensional channel flow past a cylindrical obstacle. Since the Lagrangian tracer dynamics is sensitive to the particle density and size parameters, a simple geometric setup in such flows could be used as a (low-density) particle mixture segregator.

Selective sensitivity of open chaotic flows on inertial tracer advection: catching particles with a stick

We investigate the effects of finite size and inertia of a small spherical particle immersed in an open unsteady flow which, for ideal tracers, generates transiently chaotic trajectories. The inertia effects may strongly modify the chaotic motion to the point that attractors may appear in the configuration space. These studies are performed in a model of the two-dimensional flow past a cylindrical obstacle. The relevance to modeling efforts of biological pathogen transport in large-scale flows is discussed. Since the tracer dynamics is sensitive to the particle inertia and size, simple geometric setups in such flows could be used as a particle mixture segregator separating and trapping particles.

Spatial structure of passive particles with inertia transported by a chaotic flow

We study the spatial patterns formed by inertial particles suspended on the surface of a smooth chaotic flow. In addition to the well-known phenomenon of clustering, we show that, in the presence of diffusion and when a steady space-dependent source of particles is considered, the density of particles may show smooth or fractal features in the low density areas. The conditions needed for the appearance of these structures and their characterization with the first order structure function are also calculated.