Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userDetecting chaos, determining the dimensions of tori and predicting slow diffusion in Fermi--Pasta--Ulam lattices by the Generalized Alignment Index method
114
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
The recently introduced GALI method is used for rapidly detecting chaos, determining the dimensionality of regular motion and predicting slow diffusion in multi--dimensional Hamiltonian systems. We propose an efficient computation of the GALI$_k$ indices, which represent volume elements of $k$ randomly chosen deviation vectors from a given orbit, based on the Singular Value Decomposition (SVD) algorithm. We obtain theoretically and verify numerically asymptotic estimates of GALIs long--time behavior in the case of regular orbits lying on low--dimensional tori. The GALI$_k$ indices are applied to rapidly detect chaotic oscillations, identify low--dimensional tori of Fermi--Pasta--Ulam (FPU) lattices at low energies and predict weak diffusion away from quasiperiodic motion, long before it is actually observed in the oscillations.
rate research
Read More
In systems of N coupled anharmonic oscillators, exact resonant interactions play an important role in the energy exchange between normal modes. In the weakly nonlinear regime, those interactions may facilitate energy equipartition in Fourier space. We consider analytically resonant wave-wave interactions for the celebrated Fermi-Pasta-Ulam-Tsingou (FPUT) system. Using a number-theoretical approach based on cyclotomic polynomials, we show that the problem of finding exact resonances for a system of N particles is equivalent to a Diophantine equation whose solutions depend sensitively on the set of divisors of N. We provide an algorithm to construct all possible resonances, based on two methods: pairing-off and cyclotomic, which we introduce to build up explicit solutions to the 4-, 5- and 6-wave resonant conditions. Our results shed some light in the understanding of the long-standing FPUT paradox, regarding the sensitivity of the resonant manifolds with respect to the number of particles N and the corresponding time scale of the interactions leading to thermalisation. In this light we demonstrate that 6-wave resonances always exist for any N, while 5-wave resonances exist if N is divisible by 3 and N > 6. It is known (for finite N) that 4-wave resonances do not mix energy across the spectrum, so we investigate whether 5-wave resonances can produce energy mixing across a significant region of the Fourier spectrum by analysing the interconnected network of Fourier modes that can interact nonlinearly via resonances. The answer depends on the set of odd divisors of N that are not divisible by 3: the size of this set determines the number of dynamically independent components, corresponding to independent constants of motion (energies). We show that 6-wave resonances connect all these independent components, providing in principle a restoring mechanism for full-scale thermalisation.
We study the original $alpha$-Fermi-Pasta-Ulam (FPU) system with $N=16,32$ and $64$ masses connected by a nonlinear quadratic spring. Our approach is based on resonant wave-wave interaction theory, i.e. we assume that, in the weakly nonlinear regime (the one in which Fermi was originally interested), the large time dynamics is ruled by exact resonances. After a detailed analysis of the $alpha$-FPU equation of motion, we find that the first non trivial resonances correspond to six-wave interactions. Those are precisely the interactions responsible for the thermalization of the energy in the spectrum. We predict that for small amplitude random waves the time scale of such interactions is extremely large and it is of the order of $1/epsilon^8$, where $epsilon$ is the small parameter in the system. The wave-wave interaction theory is not based on any threshold: equipartition is predicted for arbitrary small nonlinearity. Our results are supported by extensive numerical simulations. A key role in our finding is played by the {it Umklapp} (flip over) resonant interactions, typical of discrete systems. The thermodynamic limit is also briefly discussed.
The dispersive interacting waves in Fermi-Pasta-Ulam (FPU) chains of particles in textit{thermal equilibrium} are studied from both statistical and wave resonance perspectives. It is shown that, even in a strongly nonlinear regime, the chain in thermal equilibrium can be effectively described by a system of weakly interacting textit{renormalized} nonlinear waves that possess (i) the Rayleigh-Jeans distribution and (ii) zero correlations between waves, just as noninteracting free waves would. This renormalization is achieved through a set of canonical transformations. The renormalized linear dispersion of these renormalized waves is obtained and shown to be in excellent agreement with numerical experiments. Moreover, a dynamical interpretation of the renormalization of the dispersion relation is provided via a self-consistency, mean-field argument. It turns out that this renormalization arises mainly from the trivial resonant wave interactions, i.e., interactions with no momentum exchange. Furthermore, using a multiple time-scale, statistical averaging method, we show that the interactions of near-resonant waves give rise to the broadening of the resonance peaks in the frequency spectrum of renormalized modes. The theoretical prediction for the resonance width for the thermalized $beta$-FPU chain is found to be in very good agreement with its numerically measured value.
We consider a $pi$-mode solution of the Fermi-Pasta-Ulam $beta$ system. By perturbing it, we study the system as a function of the energy density from a regime where the solution is stable to a regime, where is unstable, first weakly and then strongly chaotic. We introduce, as indicator of stochasticity, the ratio $rho$ (when is defined) between the second and the first moment of a given probability distribution. We will show numerically that the transition between weak and strong chaos can be interpreted as the symmetry breaking of a set of suitable dynamical variables. Moreover, we show that in the region of weak chaos there is numerical evidence that the thermostatistic is governed by the Tsallis distribution.
We use the Smaller Alignment Index (SALI) to distinguish rapidly and with certainty between ordered and chaotic motion in Hamiltonian flows. This distinction is based on the different behavior of the SALI for the two cases: the index fluctuates around non--zero values for ordered orbits, while it tends rapidly to zero for chaotic orbits. We present a detailed study of SALIs behavior for chaotic orbits and show that in this case the SALI exponentially converges to zero, following a time rate depending on the difference of the two largest Lyapunov exponents $sigma_1$, $sigma_2$ i.e. $SALI propto e^{-(sigma_1-sigma_2)t}$. Exploiting the advantages of the SALI method, we demonstrate how one can rapidly identify even tiny regions of order or chaos in the phase space of Hamiltonian systems of 2 and 3 degrees of freedom.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
