Effect of dynamical traps on chaotic transport in a meandering jet flow

We continue our study of chaotic mixing and transport of passive particles in a simple model of a meandering jet flow [Prants et al., Chaos 16, 033117 (2006)]. In the present paper we study and phenomenologically explain a connection between dynamical, topological, and statistical properties of chaotic mixing and transport in the model flow in terms of dynamical traps, singular zones in the phase space where particles may spend an arbitrarily long but finite time [Zaslavsky, Phys. D 168-169, 292 (2002)]. The transport of passive particles is described in terms of lengths and durations of zonal flights which are events between two successive changes of sign of zonal velocity. Some peculiarities of the respective probability density functions for short flights are proven to be caused by the so-called rotational-island traps connected with the boundaries of resonant islands (including the vortex cores) filled with the particles moving in the same frame and the saddle traps connected with periodic saddle trajectories. Whereas, the statistics of long flights can be explained by the influence of the so-called ballistic-islands traps filled with the particles moving from a frame to frame.     

Directed transport of two interacting particles in a washboard potential

We study the conservative and deterministic dynamics of two nonlinearly interacting particles evolving in a one-dimensional spatially periodic washboard potential. A weak tilt of the washboard potential is applied biasing one direction for particle transport. However, the tilt vanishes asymptotically in the direction of bias. Moreover, the total energy content is not enough for both particles to be able to escape simultaneously from an initial potential well; to achieve transport the coupled particles need to interact cooperatively. For low coupling strength the two particles remain trapped inside the starting potential well permanently. For increased coupling strength there exists a regime in which one of the particles transfers the majority of its energy to the other one, as a consequence of which the latter escapes from the potential well and the bond between them breaks. Finally, for suitably large couplings, coordinated energy exchange between the particles allows them to achieve escapes — one particle followed by the other — from consecutive potential wells resulting in directed collective motion. The key mechanism of transport rectification is based on the asymptotically vanishing tilt causing a symmetry breaking of the non-chaotic fraction of the dynamics in the mixed phase space. That is, after a chaotic transient, only at one of the boundaries of the chaotic layer do resonance islands appear. The settling of trajectories in the ballistic channels associated with transporting islands provides long-range directed transport dynamics of the escaping dimer.     

Flights in a pseudo-chaotic system

We consider the problem of transport in a one-parameter family of piecewise rotations of the torus, for rotation number approaching 1∕4. This is a zero-entropy system which in this limit exhibits a divided phase space, with island chains immersed in a "pseudo-chaotic" region. We identify a novel mechanism for long-range transport, namely the adiabatic destruction of accelerator-mode islands. This process originates from the approximate translational invariance of the phase space and leads to long flights of linear motion, for a significant measure of initial conditions. We show that the asymptotic probability distribution of the flight lengths is determined by the geometric properties of a partition of the accelerator-mode island associated with the flight. We establish the existence of flights travelling distances of order O(1) in phase space. We provide evidence for the existence of a scattering process that connects flights travelling in opposite directions.     

Chaotic scattering,transport, and fractals in a simple hydrodynamic flow

Advection of passive tracers in an unsteady hydrodynamic flow consisting of a background stream and a vortex is analyzed as an example of chaotic particle scattering and transport. A numerical analysis reveals a nonattracting chaotic invariant set Λ that determines the scattering and trapping of particles from the incoming flow. The set has a hyperbolic component consisting of unstable periodic and aperiodic orbits and a nonhyperbolic component represented by marginally unstable orbits in the particle-trapping regions in the neighborhoods of the boundaries of outer invariant tori. The geometry and topology of chaotic scattering are examined. It is shown that both the trapping time for particles in the mixing region and the number of times their trajectories wind around the vortex have hierarchical fractal structure as functions of the initial particle coordinates. The hierarchy is found to have certain properties due to an infinite number of intersections of the stable manifold in Λ with a material line consisting of particles from the incoming flow. Scattering functions are singular on a Cantor set of initial conditions, and this property must manifest itself by strong fluctuations of quantities measured in experiments.     

Chaotic transmission of waves and "cooling" of signals

Ray dynamics in waveguide media exhibits chaotic motion. For a finite length of propagation, the large distance asymptotics is not uniform and represents a complicated combination of bunches of rays with different intermediate asymptotics. The origin of the phenomena that we call "chaotic transmission," lies in the nonuniformity of the phase space with sticky domains near the boundary of islands. We demonstrate different fractal properties of ray propagation using underwater acoustics as an example. The phenomenon of the kind of Levy flights can occur and it can be used as a mechanism of cooling of signals when the width of spatial spectra dispersion is significantly reduced. (c) 1997 American Institute of Physics.     

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.     

湍流扩散火焰局部熄火和再燃现象的PDF模拟

对一个值班湍流CH4／O2／N2射流扩散火焰(Sandia Flame D)进行了数值模拟研究．所采用的数学物理模型包括双尺度的k—ε湍流模型，标量联合的概率密度函数(PDF)输运方程方法，甲烷氧化的ARM简化化学反应机理(包含16种组分，12步总包反应)和欧几里德最小生成树(EMST)小尺度混合模型．将计算结果和实验数据进行了比较，不仅对于平均量，对于标量的散点分布和条件概率密度分布也是如此．计算结果表明文中采用的模型不仅能够预测宏观的火焰结构，而且预测了湍流燃烧中复杂的局部熄火和再燃过程．     

Hierarchical structures in the phase space and fractional kinetics: I. Classical systems

Hamiltonian chaotic dynamics is not ergodic due to the infinite number of islands imbedded in the stochastic sea. Stickiness of the islands' boundaries makes the wandering process very erratic with multifractal space-time structure. This complication of the chaotic process can be described on the basis of fractional kinetics. Anomalous properties of the chaotic transport become more transparent when there exists a set of islands with a hierarchical structure. Different consequences of the described phenomenon are discussed: a distribution of Poincare recurrences, characteristic exponents of transport, nonuniversality of transport, log periodicity, and chaos erasing. (c) 2000 American Institute of Physics.     

Probing the statistics of transport in the Hénon Map

The phase space of an area-preserving map typically contains infinitely many elliptic islands embedded in a chaotic sea. Orbits near the boundary of a chaotic region have been observed to stick for long times, strongly influencing their transport properties. The boundary is composed of invariant “boundary circles.” We briefly report recent results of the distribution of rotation numbers of boundary circles for the Hénon quadratic map and show that the probability of occurrence of small integer entries of their continued fraction expansions is larger than would be expected for a number chosen at random. However, large integer entries occur with probabilities distributed proportionally to the random case. The probability distributions of ratios of fluxes through island chains is reported as well. These island chains are neighbours in the sense of the Meiss-Ott Markov-tree model. Two distinct universality families are found. The distributions of the ratio between the flux and orbital period are also presented. All of these results have implications for models of transport in mixed phase space.     

Advection of passive particles in the Kolmogorov flow

A statistical analysis of the advection of passive particles in a flow governed by driven two-dimensional Navier-Stokes equations (Kolmogorov flow) is presented. Different regimes are studied, all corresponding to a chaotic behavior of the flow. The diffusion is found to be strongly asymmetric with a very weak transport perpendicular to the forcing direction. The trajectories of the particles are characterized by the presence of traps and flights. The trapping time distributions show algebraic decrease, and strong anomalous diffusion is observed in transient phases. Different regimes lead to different types of diffusion, i.e., no universal behavior of diffusion is observed, and both time and space properties are needed to define anomalous transport. (c) 2001 American Institute of Physics.     

Transient chaotic mixing during a baroclinic life cycle

We discuss how atmospheric eddies affect transport and mixing of tracers at midlatitudes. To this purpose, we study baroclinic life cycles in a simple dynamical model of the atmosphere. We consider the trapping properties of the developing eddies and the characteristics of meridional transport, and we identify regions of increased mixing. Although the flow is in principle three-dimensional, we illustrate how some of the concepts developed in the study of two-dimensional chaotic advection provide useful information on tracer dynamics in more complicated flows. (c) 2000 American Institute of Physics.     

A DNS evaluation of mixing and evaporation models for TPDF modelling of nonpremixed spray flames

In the present work, nonpremixed temporally evolving planar spray jet flames are simulated using both direct numerical simulation (DNS) and the composition transported probability density function (TPDF) method. The objective is to assess the performance of various mixing and evaporation source term distribution models which are required to close the PDF transport equation in spray flames. Quantities which would normally be provided to the TPDF solver by spray models and turbulence models are provided from the DNS: the mean flow velocity, turbulent diffusivity, mixing frequency, and cell-mean evaporation source term. Two cases with different Damköhler numbers (Da) are considered. The low Da case (Da-) features extinction followed by reignition while extinction in the high Da case (Da+) is insignificant. The TPDF modelling considers two mixing models: interaction by exchange with the mean (IEM) and Euclidean minimum spanning trees (EMST). Three models for distribution of the evaporation source terms are considered: EQUAL which distributes them in proportion to notional particles’ mass weight, NEW which creates new particles of pure fuel, and SAT which distributes the sources preferentially to notional particles close to saturation. It is found that the IEM model overpredicts the extinction when used with any evaporation model for both Da- and Da+ cases. The EMST model captures well the trend for extinction and reignition for the Da- case when it is coupled with the EQUAL evaporation model, but it overpredicts the extinction when coupled with the NEW or SAT evaporation model. For the Da+ case, all evaporation models reasonably capture the flame dynamics when coupled with EMST. The flame temperature in the mixture fraction space was examined to further assess the model performance. In general the EMST model results in narrow PDFs with little conditional fluctuation, while the IEM model produces bimodal PDFs with burning and partial extinction branches.     

Ratchet effect and the transporting islands in the chaotic sea

We study directed transport in a classical deterministic dissipative system. We consider the generic case of mixed phase space and show that large ratchet currents can be generated thanks to the presence, in the Hamiltonian limit, of transporting stability islands embedded in the chaotic sea. Because of the simultaneous presence of chaos and dissipation the stationary value of the current is independent of initial conditions, except for initial states with very small measure.     

Transport and aggregation of self-propelled particles in fluid flows

An analysis of the dynamics of prolate swimming particles in laminar flow is presented. It is shown that the particles concentrate around flow regions with chaotic trajectories. When the swimming velocity is larger than a threshold, dependent on the aspect ratio of the particles, all particles escape from regular elliptic regions. For thin rodlike particles the threshold velocity vanishes; thus, the arbitrarily small swimming velocity destroys all transport boundaries. We derive an expression for the minimum swimming velocity required for escape based on a circularly symmetric flow approximation of the regular elliptic regions.     

Effect of particles on the flow structure and dispersion of solid impurities in a two-phase axisymmetric jet

The Euler approach is used for studying the structure of a flow and the propagation of a disperse impurity in a submerged two-phase jet for small values of the mass concentration of particles (M _L1 = 0 to 0.5) upon a variation of the size and material of particles in a wide range. The effect of particles on the propagation of a two-phase jet, gas turbulence, and solid phase dispersion is analyzed. The addition of particles decreases the jet opening angle, increases the jet range, suppresses turbulence, and deteriorates turbulent mixing with the surrounding submerged space. It is shown that at the first stage, particle accumulation effects (pinching) in the axial region of the jet appear upon an increase in the particle size and the density of the particle material. Then, upon an increase in the inertia of particles, pinching changes to intense scattering of the disperse phase in the initial cross sections of the jet. The results are compared with the results of measurements for mono- and polydisperse two-phase jet flows.     

两相壁面射流PDF方程有限分析方法及数值模拟

提出求解位置-速度相空间中高维两相流PDF(probability density function)方程的有限分析方法,将位置-速度相空间颗粒PDF方程约化到速度空间,并解析求解,颗粒的位置PDF用轨道方法求解.对壁面射流两相流动进行数值模拟,并与颗粒雷诺应力轨道方法进行比较计算,结果优于颗粒雷诺应力轨道方法.     

Density and velocity statistics in variable density turbulent mixing

We experimentally study variable–density mixing of miscible gases in an open-circuit wind tunnel using simultaneous particle image velocimetry and planar laser-induced fluorescence. Experiments of a high Atwood number (0.6) and low Atwood number (0.1) are performed to compare non-Boussinesq cases with the Boussinesq limit. The higher density gas is injected into the wind tunnel co-flow using a round jet configuration, and near-field and far-field measurements are performed to examine mixing in both momentum and buoyancy-dominated regimes. The effects of buoyancy are measurable and important in both large-scale mixing features and in turbulence quantities. The low Atwood number PDFs (probability density functions) show fast and uniform mixing. The high Atwood number PDFs of density have skewness towards the larger densities, indicating less mixing of the heavy fluid due to its inertia. The skewness in the density gradient PDFs at high Atwood number displays strong density local variations that can enhance mixing at molecular scales. Turbulent kinetic energy decreases with streamwise distance from the jet for low Atwood number but increases for high Atwood number due to larger buoyancy and density-driven shear. Over 3000 experimental realisations are used to calculate statistical characteristics of the mixing, including valuable and rarely given data such as Favre-averaged turbulent quantities: mass flux velocity, Reynolds stress, turbulent kinetic energy, and density-specific volume correlation. Buoyancy effects are observed in these quantities and the trends are compared qualitatively with direct numerical simulations.     

Resonance phenomena and long-term chaotic advection in volume-preserving systems

Creating chaotic advection is the most efficient strategy to achieve mixing on microscale or in very viscous fluids. In this paper, we present a quantitative theory of the long-time resonant mixing in 3D near-integrable flows. We use the flow between two coaxial elliptic counter-rotating cylinders as a demonstrative model, where multiple scatterings on resonance result in mixing by causing the jumps of adiabatic invariants. We improve the existing estimates of the width of the mixing domain. We show that the resulting mixing both on short and long time scales can be described in terms of a single diffusion-type equation with a diffusion coefficient depending on the averaged effect of multiple passages through resonances. We discuss the exact location of the boundaries of the chaotic domain and show how it affects the properties of mixing.