Arxiv:Nlin/0311056V1 [Nlin.CD] 26 Nov 2003 Oi Rnnyeblc Nhproi Hoi Scattering, Chaotic Hyperbolic in Hyper Either Nonhyperbolic
Total Page:16
File Type:pdf, Size:1020Kb
[Phys. Rev. E 68, 056307 (2003)] Reactive dynamics of inertial particles in nonhyperbolic chaotic flows Adilson E. Motter,1, ∗ Ying-Cheng Lai,2 and Celso Grebogi3 1Max Planck Institute for the Physics of Complex Systems, N¨othnitzer Strasse 38, 01187 Dresden, Germany 2Departments of Mathematics, Electrical Engineering, and Physics, Arizona State University, Tempe, Arizona 85287, USA 3Instituto de F´ısica, Universidade de S˜ao Paulo, Caixa Postal 66318, 05315-970 S˜ao Paulo, Brazil (Dated: October 31, 2018) Anomalous kinetics of infective (e.g., autocatalytic) reactions in open, nonhyperbolic chaotic flows are im- portant for many applications in biological, chemical, and environmental sciences. We present a scaling theory for the singular enhancement of the production caused by the universal, underlying fractal patterns. The key dynamical invariant quantities are the effective fractal dimension and effective escape rate, which are primarily determined by the hyperbolic components of the underlying dynamical invariant sets. The theory is general as it includes all previously studied hyperbolic reactive dynamics as a special case. We introduce a class of dissipative embedding maps for numerical verification. PACS numbers: 47.70.Fw, 47.53.+n, 47.52.+j, 05.45.-a Many chemical and biological processes in fluids are char- the periodic orbits are unstable and there are no Kolmogorov- acterized by a filamental distribution of active particles along Arnold-Moser (KAM) tori in the phase space, while the non- fractal invariant sets of the advection chaotic dynamics. These hyperbolic counterpart is frequently characterized as having fractal structures act as dynamical catalysts for the reaction, both chaotic and marginally stable periodic orbits. Fundamen- which is relevant for a variety of environmental processes in tal assumptions in such works are: (1) the active particles are open flows, such as ozone depletion in the atmosphere [1] and massless point-like tracers and (2) the advection dynamics of population dynamics of plankton in the oceans [2]. these particles is hyperbolic [5]. However, in realistic situa- The study of active processes in open chaotic flows has tions, the Lagrangian dynamics is typically nonhyperbolic and attracted a great deal of interest from the dynamical system the active particles have finite size and inertia. Indeed, fully community [3, 4]. Most of the studies have been performed hyperbolic systems are quite rare and represent very idealized for time-dependent two-dimensional incompressible flows, in situations as the advection dynamicsof tracers in fluids is usu- the limit of weak diffusion. The flow is nonturbulent, but ally constrained to have a nonhyperbolic character because of the particle dynamics is considered to be chaotic (Lagrangian no-slip boundary conditions at the surface of obstacles. Ob- chaos) and the active particles to interact with one another stacles are at the same time the origin of Lagrangianchaos and without modifying the flow. The advection dynamics of such the origin of nonhyperbolicity. Even away from obstacles and particles can be cast in the context of chaotic scattering, where boundaries, chaotic motions of tracers typically coexist with incoming tracers spend some time in a mixing (scattering) re- regular motions. In addition, the individual active particles are gion before being scattered along the unstable manifold of the often too large to be regarded as noninertial, as is the case for chaotic saddle. As a result, the products of the reaction con- many species of zooplankton in the sea. Therefore, nonhyper- centrate along a fattened-up copy of the unstable manifold, bolic and inertial effects are prevalent in nature and expected giving rise to the observed fractal patterns. to play an important role in most environmental processes. A question of physical importance is then: what happens to the Although filamental patterns have been observed in nature, reactive dynamics when assumptions (1) and (2) are dropped? a clear relation between the observed value of the fractal di- In this article, we present a scaling theory for the reactive arXiv:nlin/0311056v1 [nlin.CD] 26 Nov 2003 mension and the underlying advection dynamics has been dynamics of inertial particles in nonhyperbolic chaotic flows. lacking. For example, in the “flow past a cylinder” system The key concepts in our framework are the effective fractal previously considered [3], the dimension of the unstable man- dimension and effective escape rate, which are respectively ifold is known to be 2 but the relevant dimension governing defined as infective and collisional reactions is about 1.6. This lack of relation, while not reducing the importance of the previous d ln N(ε) D (ε) = − , (1) phenomenological characterization of filamental distributions eff d ln ε of active particles, has led to some skepticism about the merit d ln R(n) κ (ε) = − , (2) of the dynamical system approach to the problem. In general, eff dn the advection dynamics can be characterized as either hyper- bolic or nonhyperbolic. In hyperbolic chaotic scattering, all where N(ε) is the number of ε-squares needed to cover the relevant fractal set, and R(n) is the fraction of particles that takes more than n = n(ε) steps to escape from the mixing region (see below). As a representative application of these ∗Electronic address: [email protected] concepts, we show, for autocatalytic reactions of the form A+ 2 B → 2B, that the area covered by B-particles in the steady ics is robust, nonhyperbolic chaotic scattering typically un- state obeys the following scaling law: dergoes a metamorphosis in the presence of arbitrarily small amount of dissipation. For nonhyperbolic scattering in open − 2 Deff Hamiltonian systems, particles can spend a long time in the A ∼ σ B − , (3) eµeff /(2 Deff ) − 1 neighborhood of KAM tori, resulting in an algebraic decay for the survival probability of particles in the scattering re- where µeff ≡ (κeff +κ ˜)τ, τ is the time interval between gion. As a consequence, the fractal dimension of the invariant successive reactions (time lag), σ is the reaction range, and manifolds is the phase-space dimension [11]. This should be κ˜ is the contraction rate due to dissipation. For nonhyper- contrasted with the hyperbolic case, whose decay is exponen- bolic flows in two dimensions, Deff < 2 and κeff > 0 are tial and fractal dimension is typically smaller. The dissipation, non-trivial functions of the scale ε, which can often be re- however, may convert marginally stable periodic orbits of the garded as constants over a wide interval, even though D = KAM islands into attractors. The survival probability then limε→0 Deff (ε)=2 and κ = limε→0 κeff (ε)=0. We becomes exponential, the dimension of the chaotic saddle be- find, surprisingly, that Deff and κeff are significantly differ- comes fractional, and the overall dynamics of the scattering ent from D and κ, respectively, not only for noninertial but process becomes hyperbolic. also for inertial particles, even though the advection dynamics To understand the meanings of Deff and κeff for fractal of the latter is hyperbolic, meaning that scars of the nonhy- sets arising in the transition from Hamiltonian nonhyperbolic perbolic conservative dynamics are observable in the hyper- to weakly dissipative chaotic scattering, we consider a Cantor bolic dynamics of slightly dissipative systems (κ˜ ≪ 1). The set, which is constructed in the interval [0, 1] according to the previous relations for noninertial particles in hyperbolic fluids rule that in the nth time step, a fraction ∆n = γ/(β + n)+ δ appear as a particular case of our results. is removed from the middle of each one of the N =2n−1 re- The nature of the chaotic scattering arising in the context of maining subintervals, where β, γ, and δ are constants. The particle advection in incompressible fluids may change fun- conservative case corresponds to δ = 0, which is charac- damentally as the mass and size of the particles are increased terized by an algebraic decay with n of the total length re- from zero. Physically, this happens because of the detachment maining, given by R(n) ∼ n−γ for n ≫ β, and by a unity of the particle motion from the local fluid motion. For spher- fractal dimension for the invariant set, D = 1. The re- ≡ ical particles of finite size, the particle velocity v dx/dt moved fraction ∆n decreases at each time step and, as a re- is typically different from the (time dependent) fluid velocity sult, a systematic change of scales is induced, resulting in a u = u(x,t) and, in first order, is governed by the equation [6] non-self-similar invariant set that becomes denser as we go to smaller scales. The relevant consequence is that the box- counting dimension converges slowly to 1, leading to a scale- dv du −1 − α = −a (v − u). (4) dependent effective fractal dimension, Deff ≈ 1 − γ/ ln ε dt dt for small ε. Similarly, the effective escape rate behaves as ≈ ≈ −1 The parameters are α = 3ρ /(ρ +2ρ ) and a = 2 α/St, κeff γ/n γ ln 2/ ln ε , where n = n(ε) is defined as f f p 3 the number of iterations needed to make the length of each where ρf and ρp are the densities of fluid and particle, re- spectively, and St is the Stokes number, which goes to zero remaining subinterval smaller than ε. in the limit of point-like particles [7]. For neutrally buoy- The limiting dynamics changes drastically and acquires ant particles, the mass ratio parameter is α = 1, while for properties of hyperbolic dynamics when a small amount of aerosols and bubbles, we have α < 1 and α > 1, respec- dissipation is allowed, which is modeled by 0 < δ ≪ γ/β. tively. The inertia parameter a determines the contraction rate In particular, the total length of the remaining intervals de- or dissipation in the phase space (x, v), which for incom- cays exponentially, R(n) ∼ (1 − δ)n for n ≫ γ/δ − β, pressible flows can be shown to be −2a.