Oscillating mushrooms: adiabatic theory for a non-ergodic system

Can elliptic islands contribute to sustained energy growth as parameters of a Hamiltonian system slowly vary with time? In this paper we show that a mushroom billiard with a periodically oscillating boundary accelerates the particle inside it exponentially fast. We provide an estimate for the rate of acceleration. Our numerical experiments confirms the theory. We suggest that a similar mechanism applies to general systems with mixed phase space.Comment: final revisio

An analytical study of transport, mixing and chaos in an unsteady vortical flow

We examine the transport properties of a particular two-dimensional, inviscid incompressible flow using dynamical systems techniques. The velocity field is time periodic and consists of the field induced by a vortex pair plus an oscillating strainrate field. In the absence of the strain-rate field the vortex pair moves with a constant velocity and carries with it a constant body of fluid. When the strain-rate field is added the picture changes dramatically; fluid is entrained and detrained from the neighbourhood of the vortices and chaotic particle motion occurs. We investigate the mechanism for this phenomenon and study the transport and mixing of fluid in this flow. Our work consists of both numerical and analytical studies. The analytical studies include the interpretation of the invariant manifolds as the underlying structure which govern the transport. For small values of strain-rate amplitude we use Melnikov's technique to investigate the behaviour of the manifolds as the parameters of the problem change and to prove the existence of a horseshoe map and thus the existence of chaotic particle paths in the flow. Using the Melnikov technique once more we develop an analytical estimate of the flux rate into and out of the vortex neighbourhood. We then develop a technique for determining the residence time distribution for fluid particles near the vortices that is valid for arbitrary strainrate amplitudes. The technique involves an understanding of the geometry of the tangling of the stable and unstable manifolds and results in a dramatic reduction in computational effort required for the determination of the residence time distributions. Additionally, we investigate the total stretch of material elements while they are in the vicinity of the vortex pair, using this quantity as a measure of the effect of the horseshoes on trajectories passing through this region. The numerical work verifies the analytical predictions regarding the structure of the invariant manifolds, the mechanism for entrainment and detrainment and the flux rate

Soft billiards with corners

We develop a framework for dealing with smooth approximations to billiards with corners in the two-dimensional setting. Let a polygonal trajectory in a billiard start and end up at the same billiard's corner point. We prove that smooth Hamiltonian flows which limit to this billiard have a nearby periodic orbit if and only if the polygon angles at the corner are ''acceptable''. The criterion for a corner polygon to be acceptable depends on the smooth potential behavior at the corners, which is expressed in terms of a {scattering function}. We define such an asymptotic scattering function and prove the existence of it, explain how it can be calculated and predict some of its properties. In particular, we show that it is non-monotone for some potentials in some phase space regions. We prove that when the smooth system has a limiting periodic orbit it is hyperbolic provided the scattering function is not extremal there. We then prove that if the scattering function is extremal, the smooth system has elliptic periodic orbits limiting to the corner polygon, and, furthermore, that the return map near these periodic orbits is conjugate to a small perturbation of the Henon map and therefore has elliptic islands. We find from the scaling that the island size is typically algebraic in the smoothing parameter and exponentially small in the number of reflections of the polygon orbit

Stable motions of high energy particles interacting via a repelling potential

The motion of N particles interacting by a smooth repelling potential and confined to a compact d-dimensional region is proved to be, under mild conditions, non-ergodic for all sufficiently large energies. Specifically, choreographic solutions, for which all particles follow approximately the same path close to an elliptic periodic orbit of the single-particle system, are proved to be KAM stable in the high energy limit. Finally, it is proved that the motion of N repelling particles in a rectangular box is non-ergodic at high energies for a generic choice of interacting potential: there exists a KAM-stable periodic motion by which the particles move fast only in one direction, each on its own path, yet in synchrony with all the other parallel moving particles. Thus, we prove that for smooth interaction potentials the Boltzmann ergodic hypothesis fails for a finite number of particles even in the high energy limit at which the smooth system appears to be very close to the Boltzmann hard-sphere gas

Symmetry breaking perturbations and strange attractors

The asymmetrically forced, damped Duffing oscillator is introduced as a prototype model for analyzing the homoclinic tangle of symmetric dissipative systems with \textit{symmetry breaking} disturbances. Even a slight fixed asymmetry in the perturbation may cause a substantial change in the asymptotic behavior of the system, e.g. transitions from two sided to one sided strange attractors as the other parameters are varied. Moreover, slight asymmetries may cause substantial asymmetries in the relative size of the basins of attraction of the unforced nearly symmetric attracting regions. These changes seems to be associated with homoclinic bifurcations. Numerical evidence indicates that \textit{strange attractors} appear near curves corresponding to specific secondary homoclinic bifurcations. These curves are found using analytical perturbational tools

Parabolic resonances and instabilities in near-integrable two degrees of freedom Hamiltonian flows

When an integrable two-degrees-of-freedom Hamiltonian system possessing a circle of parabolic fixed points is perturbed, a parabolic resonance occurs. It is proved that its occurrence is generic for one parameter families (co-dimension one phenomenon) of near-integrable, t.d.o. systems. Numerical experiments indicate that the motion near a parabolic resonance exhibits new type of chaotic behavior which includes instabilities in some directions and long trapping times in others. Moreover, in a degenerate case, near a {\it flat parabolic resonance}, large scale instabilities appear. A model arising from an atmospherical study is shown to exhibit flat parabolic resonance. This supplies a simple mechanism for the transport of particles with {\it small} (i.e. atmospherically relevant) initial velocities from the vicinity of the equator to high latitudes. A modification of the model which allows the development of atmospherical jets unfolds the degeneracy, yet traces of the flat instabilities are clearly observed

Non-ergodicity of the motion in three dimensional steep repelling dispersing potentials

It is demonstrated numerically that smooth three degrees of freedom Hamiltonian systems which are arbitrarily close to three dimensional strictly dispersing billiards (Sinai billiards) have islands of effective stability, and hence are non-ergodic. The mechanism for creating the islands are corners of the billiard domain.Comment: 6 pages, 8 figures, submitted to Chao

Stickiness in Hamiltonian systems: from sharply divided to hierarchical phase space

We investigate the dynamics of chaotic trajectories in simple yet physically important Hamiltonian systems with non-hierarchical borders between regular and chaotic regions with positive measures. We show that the stickiness to the border of the regular regions in systems with such a sharply divided phase space occurs through one-parameter families of marginally unstable periodic orbits and is characterized by an exponent \gamma= 2 for the asymptotic power-law decay of the distribution of recurrence times. Generic perturbations lead to systems with hierarchical phase space, where the stickiness is apparently enhanced due to the presence of infinitely many regular islands and Cantori. In this case, we show that the distribution of recurrence times can be composed of a sum of exponentials or a sum of power-laws, depending on the relative contribution of the primary and secondary structures of the hierarchy. Numerical verification of our main results are provided for area-preserving maps, mushroom billiards, and the newly defined magnetic mushroom billiards.Comment: To appear in Phys. Rev. E. A PDF version with higher resolution figures is available at http://www.pks.mpg.de/~edugal

Universal behaviour of entrainment due to coherent structures in turbulent shear flow

I suggest a solution to a persistent mystery in the physics of turbulent shear flows: cumulus clouds rise to towering heights, practically without entraining the ambient medium, while apparently similar turbulent jets in general lose their identity within a small distance through entrainment and mixing. From dynamical systems computations on a model chaotic vortical flow, I show that entrainment and mixing due to coherent structures depend sensitively on the relative speeds of different portions of the flow. A small change in these speeds, effected for example by heating, drastically alters the sizes of the KAM tori and the chaotic mixing region. The entrainment rate and, hence, the lifetime of a turbulent shear flow, shows a universal, non-monotone dependence on the heating.Comment: Preprint replaced in order to add the following comment: accepted for publication in Phys. Rev. Let