Chaotic diffusion caused by close encounters with several massive asteroids The (4) Vesta case

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Chaotic diffusion caused by close encounters with several massive asteroids: II. the regions of (10) Hygiea, (2) Pallas, and (31) Euphrosyne

Context. Close encounters with (1) Ceres and (4) Vesta, the two most massive bodies in the main belt, are known to be a mechanism of dynamical mobility able to significantly alter proper elements of minor bodies, and they are the main source of dynamical mobility for medium-sized and large asteroids (D > 20 km, approximately). Recently, it has been shown that drift rates caused by close encounters with massive asteroids may change significantly on timescales of 30 Myr when different models (i.e., different numbers of massive asteroids) are considered. Aims. So far, not much attention has been given to the case of diffusion caused by the other most massive bodies in the main belt: (2) Pallas, (10) Hygiea, and (31) Euphrosyne, the third, fourth, and one of the most massive highly inclined asteroids in the main belt, respectively. Since (2) Pallas is a highly inclined object, relative velocities at encounter with other asteroids tend to be high and changes in proper elements are therefore relatively small. It was thus believed that the scattering effect caused by highly inclined objects in general should be small. Can diffusion by close encounters with these asteroids be a significant mechanism of long-term dynamical mobility? Methods. By performing simulations with symplectic integrators, we studied the problem of scattering caused by close encounters with (2) Pallas, (10) Hygiea, and (31) Euphrosyne when only the massive asteroids (and the eight planets) are considered, and the other massive main belt asteroids and non-gravitational forces are also accounted for. Results. By finding relatively small values of drift rates for (2) Pallas, we confirm that orbital scattering by this highly inclined object is indeed a minor effect. Unexpectedly, however, we obtained values of drift rates for changes in proper semi-major axis a caused by (10) Hygiea and (31) Euphrosyne larger than what was previously found for scattering by (4) Vesta. These high rates may have repercussions on the orbital evolution and age estimate of their respective families. © 2013 ESO.

Chirikov diffusion in the asteroidal three-body resonance (5,-2,-2)

The theory of diffusion in many-dimensional Hamiltonian system is applied to asteroidal dynamics. The general formulation developed by Chirikov is applied to the NesvornA1/2-Morbidelli analytic model of three-body (three-orbit) mean-motion resonances (Jupiter-Saturn-asteroid). In particular, we investigate the diffusion along and across the separatrices of the (5, -2, -2) resonance of the (490) Veritas asteroidal family and their relationship to diffusion in semi-major axis and eccentricity. The estimations of diffusion were obtained using the Melnikov integral, a Hadjidemetriou-type sympletic map and numerical integrations for times up to 10(8) years.

Moonlets wandering on a leash-ring

Since the Voyager flybys, embedded moonlets have been proposed to explain some of the surprising structures observed in Saturn's narrow F ring. Experiments conducted with the Cassini spacecraft support this suggestion. Images of the F ring show bright compact spots, and seven occultations of stars by the F ring, monitored by ultraviolet and infrared experiments, revealed nine events of high optical depth. These results point to a large number of such objects, but it is not clear whether they are solid moonlets or rather loose particle aggregates. Subsequent images suggested an irregular motion of these objects so that a determination of their orbits consistent with the F ring failed. Some of these features seem to cross the whole ring. Here we show that these observations are explained by chaos in the F ring driven mainly by the 'shepherd' moons Prometheus and Pandora. It is characterized by a rather short Lyapunov time of about a few hundred orbital periods. Despite this chaotic diffusion, more than 93 per cent of the F-ring bodies remain confined within the F ring because of the shepherding, but also because of a weak radial mobility contrasted by an effective longitudinal diffusion. This chaotic stirring of all bodies involved prevents the formation of 'propellers' typical of moonlets, but their frequent ring crossings explain the multiple radial 'streaks' seen in the F ring. The related 'thermal' motion causes more frequent collisions between all bodies which steadily replenish F-ring dust and allow for ongoing fragmentation and re-accretion processes (ring recycling).

Controle de dinâmica caótica com toros robustos

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Analytical Studies on Diffusion and lntermittency in Chaotic Maps

The study of simple chaotic maps for non-equilibrium processes in statistical physics has been one of the central themes in the theory of chaotic dynamical systems. Recently, many works have been carried out on deterministic diffusion in spatially extended one-dimensional maps This can be related to real physical systems such as Josephson junctions in the presence of microwave radiation and parametrically driven oscillators. Transport due to chaos is an important problem in Hamiltonian dynamics also. A recent approach is to evaluate the exact diffusion coefficient in terms of the periodic orbits of the system in the form of cycle expansions. But the fact is that the chaotic motion in such spatially extended maps has two complementary aspects- - diffusion and interrnittency. These are related to the time evolution of the probability density function which is approximately Gaussian by central limit theorem. It is noticed that the characteristic function method introduced by Fujisaka and his co-workers is a very powerful tool for analysing both these aspects of chaotic motion. The theory based on characteristic function actually provides a thermodynamic formalism for chaotic systems It can be applied to other types of chaos-induced diffusion also, such as the one arising in statistics of trajectory separation. It was noted that there is a close connection between cycle expansion technique and characteristic function method. It was found that this connection can be exploited to enhance the applicability of the cycle expansion technique. In this way, we found that cycle expansion can be used to analyse the probability density function in chaotic maps. In our research studies we have successfully applied the characteristic function method and cycle expansion technique for analysing some chaotic maps. We introduced in this connection, two classes of chaotic maps with variable shape by generalizing two types of maps well known in literature.

Clustering and diffusion in a symplectic map lattice with non-local coupling

We study a symplectic chain with a non-local form of coupling by means of a standard map lattice where the interaction strength decreases with the lattice distance as a power-law, in Such a way that one can pass continuously from a local (nearest-neighbor) to a global (mean-field) type of coupling. We investigate the formation of map clusters, or spatially coherent structures generated by the system dynamics. Such clusters are found to be related to stickiness of chaotic phase-space trajectories near periodic island remnants, and also to the behavior of the diffusion coefficient. An approximate two-dimensional map is derived to explain some of the features of this connection. (C) 2008 Elsevier Ltd. All rights reserved.

Mean Motion in Stochastic Plasmas With a Space-Dependent Diffusion Coefficient

Radial transport in the tokamap, which has been proposed as a simple model for the motion in a stochastic plasma, is investigated. A theory for previous numerical findings is presented. The new results are stimulated by the fact that the radial diffusion coefficients is space-dependent. The space-dependence of the transport coefficient has several interesting effects which have not been elucidated so far. Among the new findings are the analytical predictions for the scaling of the mean radial displacement with time and the relation between the Fokker-Planck diffusion coefficient and the diffusion coefficient from the mean square displacement. The applicability to other systems is also discussed. (c) 2009 WILEY-VCH GmbH & Co. KGaA, Weinheim

Blocking Radial Diffusion in a Double-Waved Hamiltonian Model

A non-twist Hamiltonian system perturbed by two waves with particular wave numbers can present Robust Tori, barriers created by the vanishing of the perturbing Hamiltonian at some defined positions. When Robust Tori exist, any trajectory in phase space passing close to them is blocked by emergent invariant curves that prevent the chaotic transport. We analyze the breaking up of the RT as well the transport dependence on the wave numbers and on the wave amplitudes. Moreover, we report the chaotic web formation in the phase space and how this pattern influences the transport.

Scaling invariance of the diffusion coefficient in a family of two-dimensional Hamiltonian mappings

We consider a family of two-dimensional nonlinear area-preserving mappings that generalize the Chirikov standard map and model a variety of periodically forced systems. The action variable diffuses in increments whose phase is controlled by a negative power of the action and hence effectively uncorrelated for small actions, leading to a chaotic sea in phase space. For larger values of the action the phase space is mixed and contains a family of elliptic islands centered on periodic orbits and invariant Kolmogorov-Arnold-Moser (KAM) curves. The transport of particles along the phase space is considered by starting an ensemble of particles with a very low action and letting them evolve in the phase until they reach a certain height h. For chaotic orbits below the periodic islands, the survival probability for the particles to reach h is characterized by an exponential function, well modeled by the solution of the diffusion equation. On the other hand, when h reaches the position of periodic islands, the diffusion slows markedly. We show that the diffusion coefficient is scaling invariant with respect to the control parameter of the mapping when h reaches the position of the lowest KAM island. © 2013 American Physical Society.

A family of dissipative two-dimensional mappings: Chaotic, regular and steady state dynamics investigation

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

On the chaotic rotation of a liquid-filled gyrostat via the Melnikov-Holmes-Marsden integral

The non-linear motions of a gyrostat with an axisymmetrical, fluid-filled cavity are investigated. The cavity is considered to be completely filled with an ideal incompressible liquid performing uniform rotational motion. Helmholtz theorem, Euler's angular momentum theorem and Poisson equations are used to develop the disturbed Hamiltonian equations of the motions of the liquid-filled gyrostat subjected to small perturbing moments. The equations are established in terms of a set of canonical variables comprised of Euler angles and the conjugate angular momenta in order to facilitate the application of the Melnikov-Holmes-Marsden (MHM) method to investigate homoclinic/heteroclinic transversal intersections. In such a way, a criterion for the onset of chaotic oscillations is formulated for liquid-filled gyrostats with ellipsoidal and torus-shaped cavities and the results are confirmed via numerical simulations. (c) 2006 Elsevier Ltd. All rights reserved.

A bistable reaction-diffusion system in a stretching flow

We examine the evolution of a bistable reaction in a one-dimensional stretching flow, as a model for chaotic advection. We derive two reduced systems of ordinary differential equations (ODEs) for the dynamics of the governing advection-reaction-diffusion partial differential equations (PDE), for pulse-like and for plateau-like solutions, based on a non-perturbative approach. This reduction allows us to study the dynamics in two cases: first, close to a saddle-node bifurcation at which a pair of nontrivial steady states are born as the dimensionless reaction rate (Damkoehler number) is increased, and, second, for large Damkoehler number, far away from the bifurcation. The main aim is to investigate the initial-value problem and to determine when an initial condition subject to chaotic stirring will decay to zero and when it will give rise to a nonzero final state. Comparisons with full PDE simulations show that the reduced pulse model accurately predicts the threshold amplitude for a pulse initial condition to give rise to a nontrivial final steady state, and that the reduced plateau model gives an accurate picture of the dynamics of the system at large Damkoehler number. Published in Physica D (2006)

Chaotic mixing of a competitive-consecutive reaction

The evolution of a competitive-consecutive chemical reaction is computed numerically in a two-dimensional chaotic fluid flow with initially segregated reactants. Results from numerical simulations are used to evaluate a variety of reduced models commonly adopted to model the full advection-reaction-diffusion problem. Particular emphasis is placed upon fast reactions, where the yield varies most significantly with Peclet number (the ratio of diffusive to advective time scales). When effects of the fluid mechanical mixing are strongest, we find that the yield of the reaction is underestimated by a one-dimensional lamellar model that ignores the effects of fluid mixing, but overestimated by two other lamellar models that include fluid mixing.

Reduced models of chemical reaction in chaotic flows

We describe and evaluate two reduced models for nonlinear chemical reactions in a chaotic laminar flow. Each model involves two separate steps to compute the chemical composition at a given location and time. The “manifold tracking model” first tracks backwards in time a segment of the stable manifold of the requisite point. This then provides a sample of the initial conditions appropriate for the second step, which requires solving one-dimensional problems for the reaction in Lagrangian coordinates. By contrast, the first step of the “branching trajectories model” simulates both the advection and diffusion of fluid particles that terminate at the appropriate point; the chemical reaction equations are then solved along each of the branched trajectories in a second step. Results from each model are compared with full numerical simulations of the reaction processes in a chaotic laminar flow.