Dynamics of a small neutrally buoyant sphere in a fluid and targeting in Hamiltonian systems

Competing particle attractee in liquid bridges

Assuming the so-called particle accumulation structures (PAS) in liquid bridges as archetypal systems for the investigation of particle self-assembly phenomena in laminar time-periodic flows, an attempt is made here to disentangle the complex hierarchy of relationships existing between the multiplicity of the loci of aggregation (streamtubes which coexist in the physical space as competing attractee) and the particle structures effectively showing up. While the former depends on purely topological (fluid-dynamic) arguments, the influential factors driving the outcomes of the fluid-particle interaction seem to obey a much more complex logic, which makes the arrangement of particles different from realization to realization. Through numerical solution of the governing Eulerian and Lagrangian equations for liquid and mass transport, we show that for a fixed aspect ratio of the liquid bridge, particles can be gradually transferred from one streamtube to another as the Stokes number and/or the Marangoni number are varied. Moreover, ranges exist where these attractors compete resulting in overlapping or intertwined particle structures, some of which, characterized by a strong degree of asymmetry, have never been reported before. This article is part of the theme issue 'New trends in pattern formation and nonlinear dynamics of extended systems'.

Invariant manifold approach for quantifying the dynamics of highly inertial particles in steady and time-periodic incompressible flows

The dynamics of finite-sized particles with large inertia are investigated in steady and time-dependent flows through the numerical solution of the invariance equation, describing the spatiotemporal evolution of the slow/inertial manifold representing the effective particle velocity field. This approach allows for an accurate reconstruction of the effective particle divergence field, controlling clustering/dispersion features of particles with large inertia for which a perturbative approach is either inaccurate or not even convergent. The effect of inertia on heavy and light particles is quantified in terms of the rate of contraction/expansion of volume elements along a particle trajectory and of the maximum Lyapunov exponent for systems exhibiting chaotic orbits, such as the time-periodic sine-flow on the 2D torus and the time-dependent 2D cavity flow.

Particle Coherent Structures in Confined Oscillatory Switching Centrifugation

A small spherical rigid particle in a cylindrical cavity is considered. The harmonic rotation of the cavity wall drives the background flow characterized by the Strouhal number Str, assumed as the first parameter of our investigation. The particle immersed in the flow (assumed Stokesian) has a Stokes number St=1 and a particle-to-fluid density ratio ϱ which is assumed as the second parameter of this study. Building on the theoretical understanding of the recently discovered oscillatory switching centrifugation for inertial particles in unbounded flows, we investigate the effect of a confinement. For the first time we study how the presence of a wall affects the particle trajectory in oscillatory switching centrifugation dynamics. The emergence of two qualitatively different particle attractors is characterized for particles centrifuged towards the cavity wall. The implication of two such classes of attractors is discussed focusing on the application to microfluidics. We propose some strategies for exploiting the confined oscillatory switching centrifugation for selective particle segregation and for the enhancement of particle interaction events.

The fluid mechanics of poohsticks

The year 2019 marked the bicentenary of George Gabriel Stokes, who in 1851 described the drag-Stokes drag-on a body moving immersed in a fluid, and 2020 is the centenary of Christopher Robin Milne, for whom the game of poohsticks was invented; his father A. A. Milne's The House at Pooh Corner, in which it was first described in print, appeared in 1928. So this is an apt moment to review the state of the art of the fluid mechanics of a solid body in a complex fluid flow, and one floating at the interface between two fluids in motion. Poohsticks pertains to the latter category, when the two fluids are water and air. This article is part of the theme issue 'Stokes at 200 (part 2)'.

Quantitative analysis of the gain in probability of escaping for ideal phototactic swimmers due to chaotic dynamics

We study the dynamics of ideal phototactic swimmers in a steady two-dimensional model flow with transport barriers. We consider a distant light source, in which case the self-propulsion velocity of the swimmers is, at any instant, along a predetermined direction. The probability of transport along that direction emerges from the competing effects of the swimmers' self-propulsion and the flow's transport barriers. For swimmers bounded to have the same time average self-propulsion speed, temporal modulation of that speed increases the probability of escaping due to the formation of a stochastic layer which fosters transport. We use separatrix-map techniques to calculate the gain in the probability of escaping.

Turbulence induces clustering and segregation of non-motile, buoyancy-regulating phytoplankton

Turbulence plays a major role in shaping marine community structure as it affects organism dispersal and guides fundamental ecological interactions. Below oceanographic mesoscale dynamics, turbulence also impinges on subtle physical-biological coupling at the single cell level, setting a sea of chemical gradients and determining microbial interactions with profound effects on scales much larger than the organisms themselves. It has been only recently that we have started to disentangle details of this coupling for swimming microorganisms. However, for non-motile species, which comprise some of the most abundant phytoplankton groups on Earth, a similar level of mechanistic understanding is still missing. Here, we explore by means of extensive numerical simulations the interplay between buoyancy regulation in non-motile phytoplankton and cellular responses to turbulent mechanical cues. Using a minimal mechano-response model, we show how such a mechanism would contribute to spatial heterogeneity and affect vertical fluxes and trigger community segregation.

Point-source dispersion of quasi-neutrally-buoyant inertial particles

We analyze the evolution of the distribution, both in the phase space and in the physical space, of inertial particles released by a spatially-localized (punctual) source and advected by an incompressible flow. The difference in mass density between fluid and particles is assumed to be small, and represents the basic parameter for a regular perturbative expansion. By means of analytical techniques such as Hermitianization, we derive a chain of equations of the advection-diffusion-reaction type, easily solvable at least numerically. Our procedure provides results also for finite particle inertia, away from the over-damped limit of quasi-tracer dynamics.

Particle slip velocity influences inertial focusing of particles in curved microchannels

Size based separation and identification of particles in microfluidics through purely hydrodynamic means has gained significant interest due to a number of possible biomedical applications. Curved micro-channels, particularly spiral micro-channels with rectangular cross-section and the dynamics of particles in such channels have been extensively researched to achieve size based separation of particles. In this paper we present evidence that sheds new light on the dynamics of particles in such curved channels; that a unique particle slip velocity is associated with the focusing positions in the cross sections, which leads to a balance of forces. Our experiments therefore imply that the forces acting on the particle lead to convergence to an attractor in both the physical space (the cross section of the channel) and the slip velocity space.

Transport, diffusion, and energy studies in the Arnold-Beltrami-Childress map

We study the transport and diffusion properties of passive inertial particles described by a six-dimensional dissipative bailout embedding map. The base map chosen for the study is the three-dimensional incompressible Arnold-Beltrami-Childress (ABC) map chosen as a representation of volume preserving flows. There are two distinct cases: the two-action and the one-action cases, depending on whether two or one of the parameters (A,B,C) exceed 1. The embedded map dynamics is governed by two parameters (α,γ), which quantify the mass density ratio and dissipation, respectively. There are important differences between the aerosol (α<1) and the bubble (α>1) regimes. We have studied the diffusive behavior of the system and constructed the phase diagram in the parameter space by computing the diffusion exponents η. Three classes have been broadly classified-subdiffusive transport (η<1), normal diffusion (η≈1), and superdiffusion (η>1) with η≈2 referred to as the ballistic regime. Correlating the diffusive phase diagram with the phase diagram for dynamical regimes seen earlier, we find that the hyperchaotic bubble regime is largely correlated with normal and superdiffusive behavior. In contrast, in the aerosol regime, ballistic superdiffusion is seen in regions that largely show periodic dynamical behaviors, whereas subdiffusive behavior is seen in both periodic and chaotic regimes. The probability distributions of the diffusion exponents show power-law scaling for both aerosol and bubbles in the superdiffusive regimes. We further study the Poincáre recurrence times statistics of the system. Here, we find that recurrence time distributions show power law regimes due to the existence of partial barriers to transport in the phase space. Moreover, the plot of average particle kinetic energies versus the mass density ratio for the two-action case exhibits a devil's staircase-like structure for higher dissipation values. We explain these results and discuss their implications for realistic systems.

Backward Finite-Time Lyapunov Exponents in Inertial Flows

Inertial particles are finite-sized objects that are carried by fluid flows and in contrast to massless tracer particles they are subject to inertia effects. In unsteady flows, the dynamics of tracer particles have been extensively studied by the extraction of Lagrangian coherent structures (LCS), such as hyperbolic LCS as ridges of the Finite-Time Lyapunov Exponent (FTLE). The extension of the rich LCS framework to inertial particles is currently a hot topic in the CFD literature and is actively under research. Recently, backward FTLE on tracer particles has been shown to correlate with the preferential particle settling of small inertial particles. For larger particles, inertial trajectories may deviate strongly from (massless) tracer trajectories, and thus for a better agreement, backward FTLE should be computed on inertial trajectories directly. Inertial backward integration, however, has not been possible until the recent introduction of the influence curve concept, which - given an observation and an initial velocity - allows to recover all sources of inertial particles as tangent curves of a derived vector field. In this paper, we show that FTLE on the influence curve vector field is in agreement with preferential particle settling and more importantly it is not only valid for small (near-tracer) particles. We further generalize the influence curve concept to general equations of motion in unsteady spatio-velocity phase spaces, which enables backward integration with more general equations of motion. Applying the influence curve concept to tracer particles in the spatio-velocity domain emits streaklines in massless flows as tangent curves of the influence curve vector field. We demonstrate the correlation between inertial backward FTLE and the preferential particle settling in a number of unsteady vector fields.

Nonperfect mixing affects synchronization on a large number of chemical oscillators immersed in a chemically active time-dependent chaotic flow

The problem of synchronization of finite-size chemical oscillators described by active inertial particles is addressed for situations in which they are immersed in a reacting nonstationary chaotic flow. Active substances in the fluid will be modeled by Lagrangian particles closely following the fluid streamlines. Their interaction with the active inertial particles as well as the properties of the fluid dynamics will result in modifying the synchronization state of the chemical oscillators. This behavior is studied in terms of the exchange rate between the Lagrangian and inertial particles, and the finite-time Lyapunov exponents characterizing the flow. The coherence of the population of oscillators is determined by means of the order parameter introduced by Kuramoto. The different dynamics observed for the inertial particles (chemical oscillators) and Lagrangian particles (describing chemicals in the flow) lead to nonlinear interactions and patches of synchronized regions within the domain.

Lagrangian coherent structures and inertial particle dynamics

In this work we investigate the dynamics of inertial particles using finite-time Lyapunov exponents (FTLE). In particular, we characterize the attractor and repeller structures underlying preferential concentration of inertial particles in terms of FTLE fields of the underlying carrier fluid. Inertial particles that are heavier than the ambient fluid (aerosols) attract onto ridges of the negative-time fluid FTLE. This negative-time FTLE ridge becomes a repeller for particles that are lighter than the carrier fluid (bubbles). We also examine the inertial FTLE (iFTLE) determined by the trajectories of inertial particles evolved using the Maxey-Riley equations with nonzero Stokes number and density ratio. Finally, we explore the low-pass filtering effect of Stokes number. These ideas are demonstrated on two-dimensional numerical simulations of the unsteady double-gyre flow.

Clustering of inertial particles in compressible chaotic flows

Clustering of inertial particles is analyzed in chaotic compressible flows. A simplified dynamical model for the motion of inertial particles in a compressible flow has been derived. Clustering enhancement has been observed for intermediate Stokes times and characterized in terms of the number of particles with negative finite-time Lyapunov exponents and the Lyapunov dimension of the model attractor. Cluster formation has been observed to depend on the nature of the flow; vortical or shear. The motion of heavy and light particles is analyzed in terms of the compressibility and correlation length of the density field.

Microscale patches of nonmotile phytoplankton

Phytoplankton cells have evolved sophisticated strategies for actively responding to environmental signals, most notably to mechanical stresses of hydrodynamic origin. A largely unanswered question, however, is the significance of these cellular responses for the largely heterogeneous spatial distribution of cells found in the oceans. Motivated by the physiological regulation of buoyancy prevalent in nonmotile phytoplankton species, we solve here a minimal model for "active" sinking that incorporates these cellular responses. Within this model, we show how buoyancy regulation leads to intense patchiness for nonmotile species as compared to passive tracers, resulting in important variations in settling speeds and, as a consequence, determining escape rates to the deep ocean.

Dynamics of impurities in a three-dimensional volume-preserving map

We study the dynamics of inertial particles in three-dimensional incompressible maps, as representations of volume-preserving flows. The impurity dynamics has been modeled, in the Lagrangian framework, by a six-dimensional dissipative bailout embedding map. The fluid-parcel dynamics of the base map is embedded in the particle dynamics governed by the map. The base map considered for the present study is the Arnold-Beltrami-Childress (ABC) map. We consider the behavior of the system both in the aerosol regime, where the density of the particle is larger than that of the base flow, as well as the bubble regime, where the particle density is less than that of the base flow. The phase spaces in both the regimes show rich and complex dynamics with three types of dynamical behaviors--chaotic structures, regular orbits, and hyperchaotic regions. In the one-action case, the aerosol regime is found to have periodic attractors for certain values of the dissipation and inertia parameters. For the aerosol regime of the two-action ABC map, an attractor merging and widening crisis is identified using the bifurcation diagram and the spectrum of Lyapunov exponents. After the crisis an attractor with two parts is seen, and trajectories hop between these parts with period 2. The bubble regime of the embedded map shows strong hyperchaotic regions as well as crisis induced intermittency with characteristic times between bursts that scale as a power law behavior as a function of the dissipation parameter. Furthermore, we observe a riddled basin of attraction and unstable dimension variability in the phase space in the bubble regime. The bubble regime in the one-action case shows similar behavior. This study of a simple model of impurity dynamics may shed light upon the transport properties of passive scalars in three-dimensional flows. We also compare our results with those seen earlier in two-dimensional flows.

The physics of broadcast spawning in benthic invertebrates

Most benthic invertebrates broadcast their gametes into the sea, whereupon successful fertilization relies on the complex interaction between the physics of the surrounding fluid flow and the biological properties and behavior of eggs and sperm. We present a holistic overview of the impact of instantaneous flow processes on fertilization across a range of scales. At large scales, transport and stirring by the flow control the distribution of gametes. Although mean dilution of gametes by turbulence is deleterious to fertilization, a variety of instantaneous flow phenomena can aggregate gametes before dilution occurs. We argue that these instantaneous flow processes are key to fertilization efficiency. At small scales, sperm motility and taxis enhance contact rates between sperm and chemoattractant-releasing eggs. We argue that sperm motility is a biological adaptation that replaces molecular diffusion in conventional mixing processes and enables gametes to bridge the gap that remains after aggregation by the flow.

The role of structured stirring and mixing on gamete dispersal and aggregation in broadcast spawning

Broadcast-spawning benthic invertebrates synchronously release sperm and eggs from separate locations into the surrounding flow, whereupon the process depends on structured stirring by the flow field (at large scales), and sperm motility and taxis (at small scales) to bring the gametes together. The details of the relevant physical and biological aspects of the problem that result in successful and efficient fertilization are not well understood. This review paper includes relevant work from both the physical and biological communities to synthesize a more complete understanding of the processes that govern fertilization success; the focus is on the role of structured stirring on the dispersal and aggregation of gametes. The review also includes a summary of current trends and approaches for numerical and experimental simulations of broadcast spawning.

Memory effects are relevant for chaotic advection of inertial particles

A systematic investigation of the effect of the history force on particle advection is carried out in a paradigmatic model flow of chaotic advection, the von Kármán flow. All investigated properties turn out to heavily depend on the presence of memory when compared to previous studies neglecting this force. We find a weaker tendency for accumulation and for caustics formation. The Lyapunov exponent of transients becomes larger, the escape rates are strongly altered. Attractors are found to be suppressed by the history force, and periodic ones have a very slow, t(-1/2)-type convergence towards the asymptotic form.

Reduced transport of swimming particles in chaotic flow due to hydrodynamic trapping

We computationally study the transport of active, self-propelled particles suspended in a two-dimensional chaotic flow. The pointlike, spherical particles have their own intrinsic swimming velocity, which modifies the dynamical system so that the particles can break the transport barriers present in the carrier flow. Surprisingly, we find that swimming does not necessarily lead to enhanced particle transport. Small but finite swimming speed can result in reduced transport, as swimmers get stuck for long times in traps that form near elliptic islands in the background flow. Our results have implications for models of transport and encounter rates for small marine organisms.

Dynamic particle trapping, release, and sorting by microvortices on a substrate

This paper examines particle trapping and release in confined microvortex flows, including those near a solid surface due to variations in the electrokinetic slip velocity and those at a liquid-gas interface due to an external momentum source. We derive a general analytical solution for a two-dimensional microvortex flow within a semicircular cap. We also use a bifurcation theory on the kinetic equation of particles under various velocity and force fields to delineate the conditions for a vortex trap, a point trap, a limit cycle trap, and the selective sorting of the particles into different traps. In the presence of only divergence-free forces on suspended particles, we find that two parameters, such as those related to Stokes drag, gravity, and flow vorticity, are sufficient to classify all the trap topologies for a given slip velocity distribution. We also show that nondivergence-free forces such as nonuniform repulsion or attraction can capture suspended particles in one trap and selectively sort a binary suspension into different traps.

Clustering criterion for inertial particles in two-dimensional time-periodic and three-dimensional steady flows

We derive an analytic condition that predicts the exact location of inertial particle clustering in three-dimensional steady or two-dimensional time-periodic flows. The particles turn out to cluster on attracting inertial Lagrangian coherent structures that are smooth deformations of invariant tori. We illustrate our results on three-dimensional steady flows, including the Hill's spherical vortex and the Arnold-Beltrami-Childress flow, as well as on a two-dimensional time and space periodic flow that models a meandering jet in a channel.

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.

Measuring segregation of inertial particles in turbulence by a full Lagrangian approach

Preferential concentration of inertial particles in turbulence is studied numerically by evaluating the Lagrangian compressibility of the particle velocity field using the "full Lagrangian method." This is compared with the "mesoscopic Eulerian particle velocity field" both in a direct numerical simulation of turbulence and in a synthetic flow field. We demonstrate that the Lagrangian method, in contrast to the Eulerian, accurately predicts the compressibility of the particle velocity field even when the latter is characterized by singularities. In particular we use the method to evaluate the growth rates of spatial moments of the particle number density which reflect the fractal structure of segregation and the occurrence of singularities.

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.

Transport of finite-sized particles in chaotic flow

By extending traditional particle tracking techniques, we study the dynamics of neutrally buoyant finite-sized particles in a spatiotemporally chaotic flow. We simultaneously measure the flow field and the trajectories of millimeter-scale particles so that the two can be directly compared. While the single-point statistics of the particles are indistinguishable from the flow statistics, the particles often move in directions that are systematically different from the underlying flow. These differences are especially evident when Lagrangian statistics are considered.

Dynamics of tidal synchronization and orbit circularization of celestial bodies

We take a dynamical-systems approach to study the qualitative dynamical aspects of the tidal locking of the rotation of secondary celestial bodies with their orbital motion around the primary. We introduce a minimal model including the essential features of gravitationally induced elastic deformation and tidal dissipation that demonstrates the details of the energy transfer between the orbital and rotovibrational degrees of freedom. Despite its simplicity, our model can account for both synchronization into the 1:1 spin-orbit resonance and the circularization of the orbit as the only true asymptotic attractors, together with the existence of relatively long-lived metastable orbits with the secondary in p:q synchronous rotation.

Particle segregation by Stokes number for small neutrally buoyant spheres in a fluid

It is a commonly observed phenomenon that spherical particles with inertia in an incompressible fluid do not behave as ideal tracers. Due to the inertia of the particle, the planar dynamics are described in a four-dimensional phase space and thus can differ considerably from the ideal tracer dynamics. Using finite-time Lyapunov exponents, we compute the sensitivity of the final position of a particle with respect to its initial velocity, relative to the fluid, and thus partition the relative velocity subspace at each point in configuration space. The computations are done at every point in the relative velocity subspace, thus giving a sensitivity field. The Stokes number, being a measure of the independence of the particle from the underlying fluid flow, acts as a parameter in determining the variation in these partitions. We demonstrate how this partition framework can be used to segregate particles by Stokes number in a fluid. The fluid model used for demonstration is a two-dimensional cellular flow.

Quantifying turbulence-induced segregation of inertial particles

Particles with different density from the advecting turbulent fluids cluster due to the different response of light and heavy particles to turbulent fluctuations. This study focuses on the quantitative characterization of the segregation of dilute polydisperse inertial particles evolving in turbulent flow, as obtained from direct numerical simulation of homogeneous isotropic turbulence. We introduce an indicator of segregation amongst particles of different inertia and/or size, from which a length scale r_{seg}, quantifying the segregation degree between two particle types, is deduced.

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.

Clustering, chaos, and crisis in a bailout embedding map

We study the dynamics of inertial particles in two-dimensional incompressible flows. The particle dynamics is modeled by four-dimensional dissipative bailout embedding maps of the base flow which is represented by 2-d area preserving maps. The phase diagram of the embedded map is rich and interesting both in the aerosol regime, where the density of the particle is larger than that of the base flow, as well as the bubble regime, where the particle density is less than that of the base flow. The embedding map shows three types of dynamic behavior, periodic orbits, chaotic structures, and mixed regions. Thus, the embedding map can target periodic orbits as well as chaotic structures in both the aerosol and bubble regimes at certain values of the dissipation parameter. The bifurcation diagram of the 4-d map is useful for the identification of regimes where such structures can be found. An attractor merging and widening crisis is seen for a special region for the aerosols. At the crisis, two period-10 attractors merge and widen simultaneously into a single chaotic attractor. Crisis induced intermittency is seen at some points in the phase diagram. The characteristic times before bursts at the crisis show power-law behavior as functions of the dissipation parameter. Although the bifurcation diagram for the bubbles looks similar to that of aerosols, no such crisis regime is seen for the bubbles. Our results can have implications for the dynamics of impurities in diverse application contexts.

Curvature fields, topology, and the dynamics of spatiotemporal chaos

The curvature field is measured from tracer-particle trajectories in a two-dimensional fluid flow that exhibits spatiotemporal chaos and is used to extract the hyperbolic and elliptic points of the flow. These special points are pinned to the forcing when the driving is weak, but wander over the domain and interact in pairs at stronger driving, changing the local topology of the flow. Their behavior reveals a two-stage transition to spatiotemporal chaos: a gradual loss of spatial and temporal order followed by an abrupt onset of topological changes.

Interaction of solid particles with a tangle of vortex filaments in a viscous fluid

Homogeneous isotropic turbulence consists of coherent filamentary vortex structures superimposed to a more incoherent background. The question which we address is the effect of these structures on the dynamics of small, neutrally buoyant solid particles. Rather than generating the turbulence by direct numerical simulation (DNS) of the Navier-Stokes equations, we use a model of turbulence based entirely on viscous vortex filaments which interact via inertial forces and reconnect with each other. Using this model, we show that solid particles can become trapped around vortex filaments, something difficult to achieve with DNS. Unlike most studies, we have not neglected inviscid inertial effects. By comparing the Stokes, local, and convective components of the particle's acceleration, we also show that the convective part clearly identifies the trapping.

Fractal dimension in dissipative chaotic scattering

The effect of weak dissipation on chaotic scattering is relevant to situations of physical interest. We investigate how the fractal dimension of the set of singularities in a scattering function varies as the system becomes progressively more dissipative. A crossover phenomenon is uncovered where the dimension decreases relatively more rapidly as a dissipation parameter is increased from zero and then exhibits a much slower rate of decrease. We provide a heuristic theory and numerical support from both discrete-time and continuous-time scattering systems to establish the generality of this phenomenon. Our result is expected to be important for physical phenomena such as the advection of inertial particles in open chaotic flows, among others.

Signatures of fractal clustering of aerosols advected under gravity

Aerosols under chaotic advection often approach a strange attractor. They move chaotically on this fractal set but, in the presence of gravity, they have a net vertical motion downwards. In practical situations, observational data may be available only at a given level, for example, at the ground level. We uncover two fractal signatures of chaotic advection of aerosols under the action of gravity. Each one enables the computation of the fractal dimension D(0) of the strange attractor governing the advection dynamics from data obtained solely at a given level. We illustrate our theoretical findings with a numerical experiment and discuss their possible relevance to meteorology.

Embryonic nodal flow and the dynamics of nodal vesicular parcels

We address with fluid-dynamical simulations using direct numerical techniques three important and fundamental questions with respect to fluid flow within the mouse node and left-right development. First, we consider the differences between what is experimentally observed when assessing cilium-induced fluid flow in the mouse node in vitro and what is to be expected in vivo. The distinction is that in vivo, the leftward fluid flow across the mouse node takes place within a closed system and is consequently confined, while this is no longer the case on removing the covering membrane and immersing the embryo in a fluid-filled volume to perform in vitro experiments. Although there is a central leftward flow in both instances, we elucidate some important distinctions about the closed in vivo situation. Second, we model the movement of the newly discovered nodal vesicular parcels (NVPs) across the node and demonstrate that the flow should indeed cause them to accumulate on the left side of the node, as required for symmetry breaking. Third, we discuss the rupture of NVPs. Based on the biophysical properties of these vesicles, we argue that the morphogens they contain are likely not delivered to the surrounding cells by their mechanical rupture either by the cilia or the flow, and rupture must instead be induced by an as yet undiscovered biochemical mechanism.

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.

Newton's descent method for the determination of invariant tori

We formulate a fictitious-time-flow equation which drives an initial guess torus to a torus invariant under a given dynamics, provided such a torus exists. The method is general and applies in principle to continuous time flows and discrete time maps in arbitrary dimension and to both Hamiltonian and dissipative systems.

Basin topology in dissipative chaotic scattering

Chaotic scattering in open Hamiltonian systems under weak dissipation is not only of fundamental interest but also important for problems of current concern such as the advection and transport of inertial particles in fluid flows. Previous work using discrete maps demonstrated that nonhyperbolic chaotic scattering is structurally unstable in the sense that the algebraic decay of scattering particles immediately becomes exponential in the presence of weak dissipation. Here we extend the result to continuous-time Hamiltonian systems by using the Henon-Heiles system as a prototype model. More importantly, we go beyond to investigate the basin structure of scattering dynamics. A surprising finding is that, in the common case where multiple destinations exist for scattering trajectories, Wada basin boundaries are common and they appear to be structurally stable under weak dissipation, even when other characteristics of the nonhyperbolic scattering dynamics are not. We provide numerical evidence and a geometric theory for the structural stability of the complex basin topology.

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.

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).