A Study of Homoclinic/Heteroclinic Bifurcations, Homoclinic-Doubling Bifurcations and Cascades of Dynamical Systems Using a Perturbation-Incremental Method


Student thesis: Doctoral Thesis

View graph of relations


Related Research Unit(s)


Awarding Institution
Award date25 Aug 2017


In this thesis, we will study the application of the perturbation-incremental (PI) method in nonlinear dynamical systems. The thesis consists of three parts.

The first part is about an analytical approximation of heteroclinic connections in the 1:3 resonance problem by the method of nonlinear time transformation. The analysis was carried out considering the slow flow of a self-excited nonlinear Mathieu oscillator corresponding to the normal form near the 1:3 strong resonance. Two types of heteroclinic bifurcation of dimension one occur near the 1:3 resonance: the triangle and clover heteroclinic connections. Using the Hamiltonian system of slow flow of the oscillator, the unperturbed zero-order approximation of the heteroclinic connections is established. Conditions of persistence of connections in the perturbed first-order approximation of the heteroclinic connections provide close analytical approximations of the triangle and clover heteroclinic bifurcation curves, simultaneously. The analytical predictions are compared to the results obtained by numerical simulations for validation.

In the second part, homoclinic-doubling bifurcation and homoclinic-doubling cascade are studied. The normal form near the Takens-Bogdanov bifurcation from a convection problem of a square container is used as the illustrative example. Usually, a geometric approach is employed to investigate the existence of homoclinic-doubling cascades. However, this approach is not able to find analytical expression of homoclinic orbits and construct the bifurcation curves arisen from homoclinic doubling. Here, by using the nonlinear time transformation method, the algebraic approach is able to predict homoclinic-doubling bifurcation and homoclinic-doubling cascade. It has the advantage that the usual geometric consideration such as orbit flip, inclination flip, resonant homoclinic orbit and twisting is avoided. The theoretical prediction is crucial in finding the appropriate parametric values for homoclinic-doubling bifurcation analysis. In the numerical simulation, the interesting scenarios of splitting and merging of several homoclinic orbits are observed. The theoretical results agree very well with numerical simulation.

In the third part, the bifurcations of homoclinic orbits emanating from the merging of a pair of heteroclinic/homoclinic orbits are studied by the Perturbation-Incremental(PI) method. Consider strongly nonlinear self-excited oscillators
x''+g(x)=ε(μ+λ x+δ+c20x2+c11xx'+c02x'2)x',
where g(x) is an odd nonlinear function of its argument, ε is perturbation parameter, μ, λ, δ, c20, c11 and c02 are system parameters. The pair of asymmetric heteroclinic/homoclinic orbits are investigated using the nonlinear time transformation in the perturbation step when λ and δ are non-zero. We predict analytically the bifurcation point at which two branches of the pair of asymmetric orbits coexists. At this point, two new homoclinic orbits appear due to the merging of two heteroclinic/homoclinic orbits. The bifurcation curves of the new-born homoclinic orbits emanating from the bifurcation point are found by using the harmonic balance method in the incremental step. The results agree well with the numerical simulation performed by MatCont.