Nonlinear Sciences

In this paper we perform a careful analysis of the forced PP04 model for climate change, in particular the behaviour of the ice-ages. This system models the transition from a glacial to an inter-glacial state through a sudden release of oceanic Carbon Dioxide into the atmosphere. This process can be cast in terms of a Filippov dynamical system, with a discontinuous change in its dynamics related to the Carbon Dioxide release. By using techniques from the theory of non-smooth dynamical systems, we give an analysis of this model in the cases of both no insolation forcing and also periodic insolation forcing. This reveals a rich, and novel, dynamical structure to the solutions of the PP04 model. In particular we see synchronised periodic solutions with subtle regions of existence which depend on the amplitude and frequency of the forcing. The orbits can be created/destroyed in both smooth and discontinuity induced bifurcations. We study both the orbits and the transitions between them and make comparisons with actual climate dynamics.

The magnitudes of the terms in periodic orbit semiclassical trace formulas are determined by the orbits' stability exponents. In this paper, we demonstrate a simple asymptotic relationship between those stability exponents and the phase-space positions of particular homoclinic points.

We consider the system of a two ball automatic dynamic balancer attached to a rotating disc with nonconstant angular velocity. We directly compare the scenario of constant angular velocity with that when the acceleration of the rotor is taken into consideration. In doing so we show that there are cases where one must take the acceleration phase into consideration to obtain an accurate picture of the dynamics. Similarly we identify cases where the acceleration phase of the disc may be ignored. Finally, we briefly consider nonmonotonic variations of the angular velocity, with a view of maximising the basin of attraction of the desired solution, corresponding to damped vibrations.

Covariant Lyapunov vectors or CLVs span the expanding and contracting directions of perturbations along trajectories in a chaotic dynamical system. Due to efficient algorithms to compute them that only utilize trajectory information, they have been widely applied across scientific disciplines, principally for sensitivity analysis and predictions under uncertainty. In this paper, we develop a numerical method to compute the directional derivatives of CLVs along their own directions. Similar to the computation of CLVs, the present method for their derivatives is iterative and analogously uses the second-order derivative of the chaotic map along trajectories, in addition to the Jacobian. We validate the new method on a super-contracting Smale-Williams Solenoid attractor. We also demonstrate the algorithm on several other examples including smoothly perturbed Arnold Cat maps, and the Lorenz attractor, obtaining visualizations of the curvature of each attractor. Furthermore, we reveal a fundamental connection of the CLV self-derivatives with a statistical linear response formula.

The classical three-body problem arose in an attempt to understand the effect of the Sun on the Moon's Keplerian orbit around the Earth. It has attracted the attention of some of the best physicists and mathematicians and led to the discovery of chaos. We survey the three-body problem in its historical context and use it to introduce several ideas and techniques that have been developed to understand classical mechanical systems.

We use the one-dimensional Burgers equation to illustrate the effect of replacing the standard Laplacian dissipation term by a more general function of the Laplacian -- of which hyperviscosity is the best known example -- in equations of hydrodynamics. We analyze the asymptotic structure of solutions in the Fourier space at very high wave-numbers by introducing an approach applicable to a wide class of hydrodynamical equations whose solutions are calculated in the limit of vanishing Reynolds numbers from algebraic recursion relations involving iterated integrations. We give a detailed analysis of their analytic structure for two different types of dissipation: a hyperviscous and an exponentially growing dissipation term. Our results, obtained in the limit of vanishing Reynolds numbers, are validated by high-precision numerical simulations at non-zero Reynolds numbers. We then study the bottleneck problem, an intermediate asymptotics phenomenon, which in the case of the Burgers equation arises when ones uses dissipation terms (such as hyperviscosity) growing faster at high wave-numbers than the standard Laplacian dissipation term. A linearized solution of the well-known boundary layer limit of the Burgers equation involving two numerically determined parameters gives a good description of the bottleneck region.

This work is a generic advance in the study of delocalized (ergodic) to localized (non-ergodic) wave propagation phenomena in the presence of disorder. There is an urgent need to better understand the physics of extreme value process in the context of contemporary climate change. For earth system climate analysis General Circulation Model simulation sizes are rather small, 10 to 50 ensemble members due to computational burden while large ensembles are intrinsic to the study of Anderson localization. We merge universal transport approaches of Random Matrix Theory (RMT), described by the characteristic polynomial of random matrices, with the geometrical universal extremal types max stable limit law. A generic ensemble based random Hamiltonian approach allows a physical proof of state transition properties for extreme value processes. In this work Anderson localization is examined for the extreme tails of the related probability densities. We show that the Generalized Extreme Value (GEV) shape parameter ξ is a diagnostic tool that accurately distinguishes localized from delocalized systems and this property should hold for all wave based transport phenomena.

We show that Hamiltonian nonlinear dispersive wave systems with cubic nonlinearity and random initial data develop, during their evolution, anomalous correlators. These are responsible for the appearance of "ghost" excitations, i.e. those characterized by negative frequencies, in addition to the positive ones predicted by the linear dispersion relation. We use generalization of the Wick's decomposition and the wave turbulence theory to explain theoretically the existence of anomalous correlators. We test our theory on the celebrated β -Fermi-Pasta-Ulam-Tsingou chain and show that numerically measured values of the anomalous correlators agree, in the weakly nonlinear regime, with our analytical predictions. We also predict that similar phenomena will occur in other nonlinear systems dominated by nonlinear interactions, including surface gravity waves. Our results pave the road to study phase correlations in the Fourier space for weakly nonlinear dispersive wave systems.

We propose a new simple three-dimensional continuous autonomous model with two nonlinear terms and observe the dynamical behavior with respect to system parameters. This system changes the stability of fixed point via Hopf bifurcation and then undergoes a cascade of a period-doubling route to chaos. We analytically derive the first Lyapunov coefficient to investigate the nature of Hopf bifurcation and also investigate well-separated regions for different kinds of attractors in two-dimensional parameter space. Next, we introduce a time-scale ratio parameter and calculate the slow manifold using geometric singular perturbation theory. Finally, the chaotic state is annihilated by decreasing the value of the time-scale ratio parameter.

The study of deterministic chaos continues to be one of the important problems in the field of nonlinear dynamics. Interest in the study of chaos exists both in low-dimensional dynamical systems and in large ensembles of coupled oscillators. In this paper, we study the emergence of spatio-temporal chaos in chains of locally coupled identical pendulums with constant torque. The study of the scenarios of the emergence (disappearance) and properties of chaos is done as a result of changes in: (i) the individual properties of elements due to the influence of dissipation in this problem, and (ii) the properties of the entire ensemble under consideration, determined by the number of interacting elements and the strength of the connection between them. It is shown that an increase of dissipation in an ensemble with a fixed coupling force and elements number can lead to the appearance of chaos as a result of a cascade of period doubling bifurcations of periodic rotational motions or as a result of invariant tori destruction bifurcation. Chaos and hyperchaos can occur in an ensemble by adding or excluding one or more elements. Moreover, chaos arises hard, since in this case the control parameter is discrete. The influence of the coupling strength on the occurrence of chaos is specific. The appearance of chaos occurs with small and intermediate coupling and is caused by the overlap of the various out-of-phase rotational modes regions existence. The boundaries of these areas are determined analytically and confirmed in a numerical experiment. Chaotic regimes in the chain do not exist if the coupling strength is strong enough.