-
Mass crystals and cage effect in vorticity crystals
Authors:
Jean-Régis Angilella
Abstract:
We study the motion of tiny heavy inertial particles advected by a two dimensional inviscid fluid flow composed of $N$ identical point vortices regularly placed on a ring, and forming a crystal. In the limit of weak particle inertia, we show asymptotically that, in the reference frame of the crystal, inertial particles have $N$ asymptotically stable equilibrium positions located outside the crysta…
▽ More
We study the motion of tiny heavy inertial particles advected by a two dimensional inviscid fluid flow composed of $N$ identical point vortices regularly placed on a ring, and forming a crystal. In the limit of weak particle inertia, we show asymptotically that, in the reference frame of the crystal, inertial particles have $N$ asymptotically stable equilibrium positions located outside the crystal, in agreement with numerical observations by Ravichandran et al. (Sadhana 42, 2017). In addition to these "satellite" attracting points, we observe that for $N \ge 3$ the center of the ring, though degenerate, is a stable equilibrium position for inertial particles. This creates a kind of cage effect, where inclusions slowly drift towards the center under the effect of the surrounding vortices. This cage effect is observed to persist even at larger Stokes numbers, in contrast with the satellite attracting points that vanish when the Stokes number is above some critical value.
△ Less
Submitted 5 July, 2024; v1 submitted 2 July, 2024;
originally announced July 2024.
-
Chaotic advection in a steady three-dimensional MHD flow
Authors:
Julien Fontchastagner,
Jean-François Scheid,
Jean-Régis Angilella,
Jean-Pierre Brancher
Abstract:
We investigate a real 3D stationary flow characterized by chaotic advection generated by a magnetic field created by permanent magnets acting on a weakly conductive fluid subjected to a weak constant current. The model under consideration involves the Stokes equations for viscous incompressible fluid at low Reynolds number in which the density forces correspond to the Lorentz force generated by th…
▽ More
We investigate a real 3D stationary flow characterized by chaotic advection generated by a magnetic field created by permanent magnets acting on a weakly conductive fluid subjected to a weak constant current. The model under consideration involves the Stokes equations for viscous incompressible fluid at low Reynolds number in which the density forces correspond to the Lorentz force generated by the magnetic field of the magnets and the electric current through the fluid. An innovative numerical approach based on a mixed finite element method has been developed and implemented for computing the flow velocity fields with the electromagnetic force. This ensures highly accurate numerical results, allowing a detailed analysis of the chaotic behavior of fluid trajectories through the computations of associated Poincaré sections and Lyapunov exponents. Subsequently, an examination of mixing efficiency is conducted, employing computations of contamination and homogeneity rates, as well as mixing time. The obtained results underscore the relevance of the modeling and computational tools employed, as well as the design of the magnetohydrodynamic device used.
△ Less
Submitted 30 May, 2024;
originally announced May 2024.
-
Separatrix crossing due to multiplicative colored noise
Authors:
Jean-Régis Angilella
Abstract:
The effect of weak multiplicative colored noise on the dynamics of a Hamiltonian system is studied by means of asymptotic methods, in the vicinity of homoclinic or heteroclinic trajectories. A general expression for the probability of noise-induced separatrix crossing is obtained, and is illustrated by means of a two-well Duffing oscillator. It shows how weak noise can significantly affect the dyn…
▽ More
The effect of weak multiplicative colored noise on the dynamics of a Hamiltonian system is studied by means of asymptotic methods, in the vicinity of homoclinic or heteroclinic trajectories. A general expression for the probability of noise-induced separatrix crossing is obtained, and is illustrated by means of a two-well Duffing oscillator. It shows how weak noise can significantly affect the dynamics near separatrices. In addition, the influence of the degree of non-linearity of the noise amplitude is discussed.
△ Less
Submitted 14 August, 2020;
originally announced August 2020.
-
Spontaneous oscillations and negative-conductance transitions in microfluidic networks
Authors:
Daniel J. Case,
Jean-Regis Angilella,
Adilson E. Motter
Abstract:
The tendency for flows in microfluidic systems to behave linearly poses a challenge for designing integrated flow control schemes to carry out complex fluid processing tasks. This hindrance has led to the use of numerous external control devices to manipulate flows, thereby thwarting the potential scalability and portability of lab-on-a-chip technology. Here, we devise a microfluidic network exhib…
▽ More
The tendency for flows in microfluidic systems to behave linearly poses a challenge for designing integrated flow control schemes to carry out complex fluid processing tasks. This hindrance has led to the use of numerous external control devices to manipulate flows, thereby thwarting the potential scalability and portability of lab-on-a-chip technology. Here, we devise a microfluidic network exhibiting nonlinear flow dynamics that enable new mechanisms for on-chip flow control. This network is shown to exhibit oscillatory output patterns, bistable flow states, hysteresis, signal amplification, and negative-conductance transitions, all without reliance on dedicated external control hardware, movable parts, flexible components, or oscillatory inputs. These dynamics arise from nonlinear fluid inertia effects in laminar flows that we amplify and harness through the design of the network geometry. We suggest that these results, which are supported by fluid dynamical simulations and theoretical modeling, have the potential to inspire development of new built-in control capabilities, such as on-chip timing and synchronized flow patterns.
△ Less
Submitted 16 June, 2020;
originally announced June 2020.
-
Braess's paradox and programmable behaviour in microfluidic networks
Authors:
Daniel J. Case,
Yifan Liu,
Istvan Z. Kiss,
Jean-Regis Angilella,
Adilson E. Motter
Abstract:
Microfluidic systems are now being designed with precision to execute increasingly complex tasks. However, their operation often requires numerous external control devices due to the typically linear nature of microscale flows, which has hampered the development of integrated control mechanisms. We address this difficulty by designing microfluidic networks that exhibit a nonlinear relation between…
▽ More
Microfluidic systems are now being designed with precision to execute increasingly complex tasks. However, their operation often requires numerous external control devices due to the typically linear nature of microscale flows, which has hampered the development of integrated control mechanisms. We address this difficulty by designing microfluidic networks that exhibit a nonlinear relation between applied pressure and flow rate, which can be harnessed to switch the direction of internal flows solely by manipulating input and/or output pressures. We show that these networks exhibit an experimentally-supported fluid analog of Braess's paradox, in which closing an intermediate channel results in a higher, rather than lower, total flow rate. The harnessed behavior is scalable and can be used to implement flow routing with multiple switches. These findings have the potential to advance development of built-in control mechanisms in microfluidic networks, thereby facilitating the creation of portable systems that may one day be as controllable as microelectronic circuits.
△ Less
Submitted 27 May, 2020;
originally announced May 2020.
-
Probability of noise-induced separatrix crossing for inertial particles in flows
Authors:
Jean-Régis Angilella
Abstract:
The motion of weakly inertial Brownian particles, transported by steady two-dimensional fluid flows, is investigated by means of asymptotic methods. We focus on the phenomenon of noise-induced separatrix crossing, which can force particles to enter or exit recirculation cells in an unpredictable manner. An analytical expression for the probability of separatrix crossing is obtained. It can be appl…
▽ More
The motion of weakly inertial Brownian particles, transported by steady two-dimensional fluid flows, is investigated by means of asymptotic methods. We focus on the phenomenon of noise-induced separatrix crossing, which can force particles to enter or exit recirculation cells in an unpredictable manner. An analytical expression for the probability of separatrix crossing is obtained. It can be applied to a wide variety of flows, provided some elementary kinematical quantities of the fluid flow are known. It does not require to solve particle trajectories.
△ Less
Submitted 7 May, 2019;
originally announced May 2019.
-
The effect of colored noise on heteroclinic orbits
Authors:
Jean-Régis Angilella
Abstract:
The dynamics of a weakly dissipative Hamiltonian system submitted to stochastic perturbations has been investigated by means of asymptotic methods. The probability of noise-induced separatrix crossing, which drastically changes the fate of the system, is derived analytically in the case where noise is an additive Kubo-Anderson process. This theory shows how the geometry of the separatrix, as well…
▽ More
The dynamics of a weakly dissipative Hamiltonian system submitted to stochastic perturbations has been investigated by means of asymptotic methods. The probability of noise-induced separatrix crossing, which drastically changes the fate of the system, is derived analytically in the case where noise is an additive Kubo-Anderson process. This theory shows how the geometry of the separatrix, as well as the noise intensity and correlation time, affect the statistics of crossing. Results can be applied to a wide variety of systems, and are valid in the limit where the noise correlation time scale is much smaller than the time scale of the undisturbed Hamiltonian dynamics.
△ Less
Submitted 27 March, 2019;
originally announced March 2019.
-
Intrinsic viscosity of a suspension of weakly Brownian ellipsoids in shear
Authors:
G. Almondo,
J. Einarsson,
J. R. Angilella,
B. Mehlig
Abstract:
We analyze the angular dynamics of triaxial ellipsoids in a shear flow subject to weak thermal noise. By numerically integrating an overdamped angular Langevin equation, we find the steady angular probability distribution for a range of triaxial particle shapes. From this distribution we compute the intrinsic viscosity of a dilute suspension of triaxial particles. We determine how the viscosity de…
▽ More
We analyze the angular dynamics of triaxial ellipsoids in a shear flow subject to weak thermal noise. By numerically integrating an overdamped angular Langevin equation, we find the steady angular probability distribution for a range of triaxial particle shapes. From this distribution we compute the intrinsic viscosity of a dilute suspension of triaxial particles. We determine how the viscosity depends on particle shape in the limit of weak thermal noise. While the deterministic angular dynamics depends very sensitively on particle shape, we find that the shape dependence of the intrinsic viscosity is weaker, in general, and that suspensions of rod-like particles are the most sensitive to breaking of axisymmetry. The intrinsic viscosity of a dilute suspension of triaxial particles is smaller than that of a suspension of axisymmetric particles with the same volume, and the same ratio of major to minor axis lengths.
△ Less
Submitted 13 June, 2018; v1 submitted 19 May, 2017;
originally announced May 2017.
-
Levitation of heavy particles against gravity in asymptotically downward flows
Authors:
Jean-Regis Angilella,
Daniel J. Case,
Adilson E. Motter
Abstract:
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 magnitud…
▽ More
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.
△ Less
Submitted 7 March, 2017;
originally announced March 2017.
-
Chaotic sedimentation of particle pairs in a vertical channel at low Reynolds number: multiple states and routes to chaos
Authors:
Romuald Verjus,
Sylvain Guillou,
Alexander Ezersky,
Jean-Régis Angilella
Abstract:
The sedimentation of a pair of rigid circular particles in a two-dimensional vertical channel containing a Newtonian fluid is investigated numerically, for terminal particle Reynolds numbers ranging from 1 to 10, and for a confinement ratio equal to 4. While it is widely admitted that sufficiently inertial pairs should sediment by performing a regular DKT oscillation (Drafting-Kissing-Tumbling), t…
▽ More
The sedimentation of a pair of rigid circular particles in a two-dimensional vertical channel containing a Newtonian fluid is investigated numerically, for terminal particle Reynolds numbers ranging from 1 to 10, and for a confinement ratio equal to 4. While it is widely admitted that sufficiently inertial pairs should sediment by performing a regular DKT oscillation (Drafting-Kissing-Tumbling), the present analysis shows in contrast that a chaotic regime can also exist for such particles, leading to a much slower sedimentation velocity. It consists of a nearly horizontal pair, corresponding to a maximum effective blockage ratio, and performing a quasiperiodic transition to chaos under increasing the particle weight. For less inertial regimes, the classical oblique doublet structure and its complex behavior (multiple stable states and hysteresis, period-doubling cascade and chaotic attractor) are recovered, in agreement with previous work [Aidun & Ding, Physics of Fluids 15(6), 2003]. As a consequence of these various behaviors, the link between the terminal Reynolds number and the non-dimensional driving force is complex: it contains several branches displaying hysteresis as well as various bifurcations. For the range of Reynolds number considered here, a global bifurcation diagram is given.
△ Less
Submitted 26 November, 2016;
originally announced November 2016.
-
The role of inertia for the rotation of a nearly spherical particle in a general linear flow
Authors:
F. Candelier,
J. Einarsson,
F. Lundell,
B. Mehlig,
J. -R. Angilella
Abstract:
We analyse the angular dynamics of a neutrally buoyant nearly spherical particle immersed in a steady general linear flow. The hydrodynamic torque acting on the particle is obtained by means of a reciprocal theorem, regular perturbation theory exploiting the small eccentricity of the nearly spherical particle, and assuming that inertial effects are small, but finite.
We analyse the angular dynamics of a neutrally buoyant nearly spherical particle immersed in a steady general linear flow. The hydrodynamic torque acting on the particle is obtained by means of a reciprocal theorem, regular perturbation theory exploiting the small eccentricity of the nearly spherical particle, and assuming that inertial effects are small, but finite.
△ Less
Submitted 11 May, 2015;
originally announced May 2015.
-
Effect of weak fluid inertia upon Jeffery orbits
Authors:
J. Einarsson,
F. Candelier,
F. Lundell,
J. R. Angilella,
B. Mehlig
Abstract:
We consider the rotation of small neutrally buoyant axisymmetric particles in a viscous steady shear flow. When inertial effects are negligible the problem exhibits infinitely many periodic solutions, the "Jeffery orbits". We compute how inertial effects lift their degeneracy by perturbatively solving the coupled particle-flow equations. We obtain an equation of motion valid at small shear Reynold…
▽ More
We consider the rotation of small neutrally buoyant axisymmetric particles in a viscous steady shear flow. When inertial effects are negligible the problem exhibits infinitely many periodic solutions, the "Jeffery orbits". We compute how inertial effects lift their degeneracy by perturbatively solving the coupled particle-flow equations. We obtain an equation of motion valid at small shear Reynolds numbers, for spheroidal particles with arbitrary aspect ratios. We analyse how the linear stability of the \lq log-rolling\rq{} orbit depends on particle shape and find it to be unstable for prolate spheroids. This resolves a puzzle in the interpretation of direct numerical simulations of the problem. In general both unsteady and non-linear terms in the Navier-Stokes equations are important.
△ Less
Submitted 11 April, 2015;
originally announced April 2015.
-
Rotation of a spheroid in a simple shear at small Reynolds number
Authors:
J. Einarsson,
F. Candelier,
F. Lundell,
J. R. Angilella,
B. Mehlig
Abstract:
We derive an effective equation of motion for the orientational dynamics of a neutrally buoyant spheroid suspended in a simple shear flow, valid for arbitrary particle aspect ratios and to linear order in the shear Reynolds number. We show how inertial effects lift the degeneracy of the Jeffery orbits and determine the stabilities of the log-rolling and tumbling orbits at infinitesimal shear Reyno…
▽ More
We derive an effective equation of motion for the orientational dynamics of a neutrally buoyant spheroid suspended in a simple shear flow, valid for arbitrary particle aspect ratios and to linear order in the shear Reynolds number. We show how inertial effects lift the degeneracy of the Jeffery orbits and determine the stabilities of the log-rolling and tumbling orbits at infinitesimal shear Reynolds numbers. For prolate spheroids we find stable tumbling in the shear plane, log-rolling is unstable. For oblate particles, by contrast, log-rolling is stable and tumbling is unstable provided that the aspect ratio is larger than a critical value. When the aspect ratio is smaller than this value tumbling turns stable, and an unstable limit cycle is born.
△ Less
Submitted 9 April, 2015;
originally announced April 2015.
-
Tumbling of asymmetric microrods in a microchannel flow
Authors:
J. Einarsson,
B. M. Mihiretie,
A. Laas,
S. Ankardal,
J. R. Angilella,
D. Hanstorp,
B. Mehlig
Abstract:
We describe results of measurements of the orientational motion of glass microrods in a microchannel flow, following the orientational motion of particles with different shapes. We determine how the orientational dynamics depends on the shape of the particle and on its initial orientation. We find that the dynamics depends so sensitively on the degree to which particle axisymmetry is broken that i…
▽ More
We describe results of measurements of the orientational motion of glass microrods in a microchannel flow, following the orientational motion of particles with different shapes. We determine how the orientational dynamics depends on the shape of the particle and on its initial orientation. We find that the dynamics depends so sensitively on the degree to which particle axisymmetry is broken that it is difficult to find particles that are sufficiently axisymmetric so that they exhibit periodic tumbling ("Jeffery orbits"). The results of our measurements confirm earlier theoretical analysis predicting sensitive dependence on particle shape and its initial orientation. Our results illustrate the different types of orientational dynamics for asymmetric particles predicted by theory.
△ Less
Submitted 18 December, 2015; v1 submitted 10 March, 2015;
originally announced March 2015.
-
Critical Stokes number for the capture of inertial particles by recirculation cells in 2D quasi-steady flows
Authors:
Romuald Verjus,
Jean-Régis Angilella
Abstract:
Inertial particles are often observed to be trapped, temporarily or permanently, by recirculation cells which are ubiquitous in natural or industrial flows. In the limit of small particle inertia, determining the conditions of trapping is a challenging task, as it requires a large number of numerical simulations or experiments to test various particle sizes or densities. Here, we investigate this…
▽ More
Inertial particles are often observed to be trapped, temporarily or permanently, by recirculation cells which are ubiquitous in natural or industrial flows. In the limit of small particle inertia, determining the conditions of trapping is a challenging task, as it requires a large number of numerical simulations or experiments to test various particle sizes or densities. Here, we investigate this phenomenon analytically and numerically in the case of heavy particles (e.g. aerosols) at low Reynolds number, to derive a trapping criterion that can be used both in analytical and numerical velocity fields. The resulting criterion allows to predict the characteristics of trapped particles as soon as single-phase simulations of the flow are performed. Our analysis is valid for two-dimensional particle-laden flows in the vertical plane, in the limit where the particle inertia, the free-fall terminal velocity, and the flow unsteadiness can be treated as perturbations. The weak unsteadiness of the flow generally induces a chaotic tangle near heteroclinic or homoclinic cycles if any, leading to the apparent diffusion of fluid elements through the boundary of the cell. The critical particle Stokes number Stc below which aerosols also enter and exit the cell in a complex manner has been derived analytically, in terms of the flow characteristics. It involves the non-dimensional curvature-weighted integral of the squared velocity of the steady fluid flow along the dividing streamline of the recirculation cell. When the flow is unsteady and St > Stc, a regular motion takes place due to gravity and centrifugal effects, like in the steady case. Particles driven towards the interior of the cell are trapped permanently. In contrast, when the flow is unsteady and St < Stc, particles wander in a chaotic manner in the vicinity of the border of the cell, and can escape the cell.
△ Less
Submitted 19 May, 2016; v1 submitted 24 February, 2015;
originally announced February 2015.
-
Inertial particle trapping in an open vortical flow
Authors:
Jean-Regis Angilella,
Rafael D. Vilela,
Adilson E. Motter
Abstract:
Recent numerical results on advection dynamics have shown that particles denser than the fluid can remain trapped indefinitely in a bounded region of an open fluid flow. Here, we investigate this counterintuitive phenomenon both numerically and analytically to establish the conditions under which the underlying particle-trapping attractors can form. We focus on a two-dimensional open flow composed…
▽ More
Recent numerical results on advection dynamics have shown that particles denser than the fluid can remain trapped indefinitely in a bounded region of an open fluid flow. Here, we investigate this counterintuitive phenomenon both numerically and analytically to establish the conditions under which the underlying particle-trapping attractors can form. We focus on a two-dimensional open flow composed of a pair of vortices and its specular image, which is a system we represent as a vortex pair plus a wall along the symmetry line. Considering particles that are much denser than the fluid, we show that two point attractors form in the neighborhood of the vortex pair provided that the particle Stokes number is smaller than a critical value of order unity. In the absence of the wall, the boundaries of the basins of the attracting points are smooth. When the wall is present, the point attractors describe counter-rotating ellipses in this frame, with a period equal to half the period of one isolated vortex pair. However, their boundaries are shown to become fractal if the distance to the wall is smaller than a critical distance that scales with the inverse square root of the Stokes number. This transformation is related to the breakdown of a separatrix that gives rise to a heteroclinic tangle close to the vortices, which we describe using a Melnikov function. For an even smaller distance to the wall, we demonstrate that a second separatrix breaks down and a new heteroclinic tangle forms farther away from the vortices. Particles released in the open part of the flow can approach the attractors and be trapped permanently provided that they cross the two separatrices. Furthermore, the trapping of heavy particles from the open flow is shown to be robust to the presence of gravity, viscosity, and noise.
△ Less
Submitted 28 March, 2014;
originally announced March 2014.
-
Orientational dynamics of weakly inertial axisymmetric particles in steady viscous flows
Authors:
J. Einarsson,
J. R. Angilella,
B. Mehlig
Abstract:
The orientational dynamics of weakly inertial axisymmetric particles in a steady flow is investigated. We derive an asymptotic equation of motion for the unit axial vector along the particle symmetry axis, valid for small Stokes number St, and for any axisymmetric particle in any steady linear viscous flow. This reduced dynamics is analysed in two ways, both pertain to the case of a simple shear f…
▽ More
The orientational dynamics of weakly inertial axisymmetric particles in a steady flow is investigated. We derive an asymptotic equation of motion for the unit axial vector along the particle symmetry axis, valid for small Stokes number St, and for any axisymmetric particle in any steady linear viscous flow. This reduced dynamics is analysed in two ways, both pertain to the case of a simple shear flow. In this case inertia induces a coupling between precession and nutation. This coupling affects the dynamics of the particle, breaks the degeneracy of the Jeffery orbits, and creates two limiting periodic orbits. We calculate the leading-order Floquet exponents of the limiting periodic orbits and show analytically that prolate objects tend to a tumbling orbit, while oblate objects tend to a log-rolling orbit, in agreement with previous analytical and numerical results. Second, we analyse the role of the limiting orbits when rotational noise is present. We formulate the Fokker-Planck equation describing the orientational distribution of an axisymmetric particle, valid for small St and general Peclet number Pe. Numerical solutions of the Fokker-Planck equation, obtained by means of expansion in spherical harmonics, show that stationary orientational distributions are close to the inertia-free case when Pe St << 1, whereas they are determined by inertial effects, though small, when Pe >> 1/St >> 1.
△ Less
Submitted 10 July, 2013;
originally announced July 2013.
-
Aperiodic tumbling of microrods advected in a microchannel flow
Authors:
J. Einarsson,
A. Johansson,
S. K. Mahato,
Y. Mishra,
J. R. Angilella,
D. Hanstorp,
B. Mehlig
Abstract:
We report on an experimental investigation of the tumbling of microrods in the shear flow of a microchannel (40 x 2.5 x 0.4 mm). The rods are 20 to 30 microns long and their diameters are of the order of 1 micron. Images of the centre-of-mass motion and the orientational dynamics of the rods are recorded using a microscope equipped with a CCD camera. A motorised microscope stage is used to track i…
▽ More
We report on an experimental investigation of the tumbling of microrods in the shear flow of a microchannel (40 x 2.5 x 0.4 mm). The rods are 20 to 30 microns long and their diameters are of the order of 1 micron. Images of the centre-of-mass motion and the orientational dynamics of the rods are recorded using a microscope equipped with a CCD camera. A motorised microscope stage is used to track individual rods as they move along the channel. Automated image analysis determines the position and orientation of a tracked rods in each video frame. We find different behaviours, depending on the particle shape, its initial position, and orientation. First, we observe periodic as well as aperiodic tumbling. Second, the data show that different tumbling trajectories exhibit different sensitivities to external perturbations. These observations can be explained by slight asymmetries of the rods. Third we observe that after some time, initially periodic trajectories lose their phase. We attribute this to drift of the centre of mass of the rod from one to another stream line of the channel flow.
△ Less
Submitted 30 May, 2013; v1 submitted 24 November, 2012;
originally announced November 2012.
-
Asymptotic properties of wall-induced chaotic mixing in point vortex pairs
Authors:
Jean-Regis Angilella
Abstract:
The purpose of this work is to analyze the flow due to a potential point vortex pair in the vicinity of a symmetry line (or "wall"), in order to understand why the presence of the wall, even far from the vortices, accelerates fluid mixing around the vortex pair. An asymptotic analysis, in the limit of large distances to the wall, allows to approximate the wall effect as a constant translation of t…
▽ More
The purpose of this work is to analyze the flow due to a potential point vortex pair in the vicinity of a symmetry line (or "wall"), in order to understand why the presence of the wall, even far from the vortices, accelerates fluid mixing around the vortex pair. An asymptotic analysis, in the limit of large distances to the wall, allows to approximate the wall effect as a constant translation of the vortex pair parallel to the wall, plus a straining flow which induces a natural blinking vortex mechanism with period half the rotation period. A Melnikov analysis of lagrangian particles, in the frame translating and rotating with the vortices, shows that a homoclinic bifurcation indeed occurs, so that the various separatrices located near the vortex pair (and rotating with it) do not survive when a wall is present. The thickness of the resulting inner stochastic layer is estimated by using the separatrix map method, and is shown to scale like the inverse of the squared distance to the wall. This estimation provides a lower-bound to the numerical thickness measured from either Poincaré sections or simulations of lagrangian particles transported by the exact potential velocity field in the laboratory frame. In addition, it is shown that the outer homoclinic cycle, separating the vortices from the external (open) flow, is also perturbed from inside by the rotation of the vortex pair. As a consequence, a stochastic layer is shown to exist also in the vicinity of this cycle, allowing fluid exchange between the vortices and the outer flow. However, the thickness of this outer stochastic zone is observed to be much smaller than the one of the inner stochastic zone near vortices, as soon as the distance to the wall is large enough.
△ Less
Submitted 18 October, 2011; v1 submitted 12 September, 2011;
originally announced September 2011.
-
Note on dust trapping in inviscid vortex pairs with unequal strengths
Authors:
Tatiana Nizkaya,
Jean-Regis Angilella,
Michel Bues
Abstract:
We investigate theoretically the motion of tiny heavy passive particles transported in a plane inviscid flow consisting of two point vortices, in order to understand particle dispersion and trapping during vortex interaction. In spite of their large density, particles are not necessarily centrifugated away from vortices. It is observed that they can have various equilibrium positions in the refer…
▽ More
We investigate theoretically the motion of tiny heavy passive particles transported in a plane inviscid flow consisting of two point vortices, in order to understand particle dispersion and trapping during vortex interaction. In spite of their large density, particles are not necessarily centrifugated away from vortices. It is observed that they can have various equilibrium positions in the reference frame rotating with the vortices, provided the particle response time and the vortex strength ratio lie in appropriate ranges. A stability analysis reveals that some of these points can be asymptotically stable, and can therefore trap particles released in their basin of attraction. A complete trapping diagram is derived, showing that any vortex pair can potentially become a dust trap, provided the vortex strength ratio is different from 0 (single vortex) and -1 (translating symmetrical vortices). Trapping exists for both co-rotating or contra-rotating vortex pairs. In the latter case, particle trapping on a limit cycle is also observed, and confirmed by using Sapsis and Haller's method [Chaos, 20, 017515, 2010] generalized to non-inertial reference frames.
△ Less
Submitted 7 April, 2010; v1 submitted 24 March, 2010;
originally announced March 2010.
-
Asymptotic analysis of chaotic particle sedimentation and trapping in the vicinity of a vertical upward streamline
Authors:
J. R. Angilella
Abstract:
The sedimentation of a heavy Stokes particle in a laminar plane or axisymmetric flow is investigated by means of asymptotic methods. We focus on the occurrence of Stommel's retention zones, and on the splitting of their separatrices. The goal of this paper is to analyze under which conditions these retention zones can form, and under which conditions they can break and induce chaotic particle set…
▽ More
The sedimentation of a heavy Stokes particle in a laminar plane or axisymmetric flow is investigated by means of asymptotic methods. We focus on the occurrence of Stommel's retention zones, and on the splitting of their separatrices. The goal of this paper is to analyze under which conditions these retention zones can form, and under which conditions they can break and induce chaotic particle settling. The terminal velocity of the particle in still fluid is of the order of the typical velocity of the flow, and the particle response time is much smaller than the typical flow time-scale. It is observed that if the flow is steady and has an upward streamline where the vertical velocity has a strict local maximum, then inertialess particle trajectories can take locally the form of elliptic Stommel cells, provided the particle terminal velocity is close enough to the local peak flow velocity. These structures only depend on the local flow topology and do not require the flow to have closed streamlines or stagnation points. If, in addition, the flow is submitted to a weak time-periodic perturbation, classical expansions enable one to write the particle dynamics as a hamiltonian system with one degree of freedom, plus a perturbation containing both the dissipative terms of the particle motion equation and the flow unsteadiness. Melnikov's method therefore provides accurate criteria to predict the splitting of the separatrices of the elliptic cell mentioned above, leading to chaotic particle trapping and chaotic settling. The effect of particle inertia and flow unsteadiness on the occurrence of simple zeros in Melnikov's function is discussed. Particle motion in a plane cellular flow and in a vertical pipe is then investigated to illustrate these results.
△ Less
Submitted 22 March, 2010;
originally announced March 2010.
-
Chaotic particle sedimentation in a rotating flow with time-periodic strength
Authors:
J. R. Angilella
Abstract:
Particle sedimentation in the vicinity of a fixed horizontal vortex with time-dependent intensity can be chaotic, provided gravity is sufficient to displace the particle cloud while the vortex is off or weak. This "stretch, sediment and fold" mechanism is close to the so-called blinking vortex effect, which is responsible for chaotic transport of perfect tracers, except that in the present case…
▽ More
Particle sedimentation in the vicinity of a fixed horizontal vortex with time-dependent intensity can be chaotic, provided gravity is sufficient to displace the particle cloud while the vortex is off or weak. This "stretch, sediment and fold" mechanism is close to the so-called blinking vortex effect, which is responsible for chaotic transport of perfect tracers, except that in the present case the vortex motion is replaced by gravitational settling. In the present work this phenomenon is analyzed for heavy Stokes particles moving under the sole effect of gravity and of a linear drag. The vortex is taken to be a fixed isolated point vortex the intensity of which varies under the effect of either boundary conditions or volume force. When the unsteadiness of the vortex is weak and the free-fall velocity is of the order of the fluid velocity, and the particle response time is small, the particle motion equation can be written asymptotically as a perturbed hamiltonian system the phase portrait of which displays a homoclinic trajectory. A homoclinic bifurcation is therefore likely to occur, and the contribution of particle inertia to the occurrence of this bifurcation is analyzed asymptotically by using Melnikov's method.
△ Less
Submitted 21 March, 2010;
originally announced March 2010.
-
Dust trapping in inviscid vortex pairs
Authors:
Jean-Regis Angilella
Abstract:
The motion of tiny heavy particles transported in a co-rotating vortex pair, with or without particle inertia and sedimentation, is investigated. The dynamics of non-inertial sedimenting particles is shown to be chaotic, under the combined effect of gravity and of the circular displacement of the vortices. This phenomenon is very sensitive to particle inertia, if any. By using nearly hamiltonian…
▽ More
The motion of tiny heavy particles transported in a co-rotating vortex pair, with or without particle inertia and sedimentation, is investigated. The dynamics of non-inertial sedimenting particles is shown to be chaotic, under the combined effect of gravity and of the circular displacement of the vortices. This phenomenon is very sensitive to particle inertia, if any. By using nearly hamiltonian dynamical system theory for the particle motion equation written in the rotating reference frame, one can show that small inertia terms of the particle motion equation strongly modify the Melnikov function of the homoclinic trajectories and heteroclinic cycles of the unperturbed system, as soon as the particle response time is of the order of the settling time (Froude number of order unity). The critical Froude number above which chaotic motion vanishes and a regular centrifugation takes place is obtained from this Melnikov analysis and compared to numerical simulations. Particles with a finite inertia, and in the absence of gravity, are not necessarily centrifugated away from the vortex system. Indeed, these particles can have various equilibrium positions in the rotating reference frame, like the Lagrange points of celestial mechanics, according to whether their Stokes number is smaller or larger than some critical value. An analytical stability analysis reveals that two of these points are stable attracting points, so that permanent trapping can occur for inertial particles injected in an isolated co-rotating vortex pair. Particle trapping is observed to persist when viscosity, and therefore vortex coalescence, is taken into account. Numerical experiments at large but finite Reynolds number show that particles can indeed be trapped temporarily during vortex roll-up, and are eventually centrifugated away once vortex coalescence occurs.
△ Less
Submitted 19 March, 2010;
originally announced March 2010.