Chaotic Motion Research Articles

Applying Zero-Crossing Method for Frequency Map Analysis of Dynamical Systems

Abstract We present a simple method for estimating fundamental frequencies from orbit’s time series of triaxial dynamical systems based on the zero-crossing technique. The simplicity and high performance of the method make it a viable alternative to more advanced, but rather complex and computationally expensive techniques based on Fourier spectral analysis and iterative torus-mapping with action-angle variables. More importantly, the method provides a simple and quick way to identify irregular orbits and quantify their level of chaoticity, allowing to detect even a subtle onset of chaos. To validate the proposed algorithm, we performed comparison tests with synthetic quasi-periodic signals and with other methods using the triaxial Dehnen potential. We found noticeable differences in frequency maps constructed with our method and methods based on Fourier spectral analysis that we attributed to the way we assign static properties to chaotic orbits. At the same time, our frequency maps turn out to be in very good agreement with the maps obtained using the action-based technique.

Analysis of quasi-periodic and chaotic motion of a dielectric elastomer shell under alternating voltage

Analysis of quasi-periodic and chaotic motion of a dielectric elastomer shell under alternating voltage

Vibration energy harvesting system with cyclically time-varying potential barrier

Nonlinear kinetic energy harvesters are becoming more and more popular as well as advanced and efficient. This paper presents the study of the dynamics of such a system in a wide range of excitation parameters, assuming at the same time the possibility of a cyclical and smooth change of the potential function. We have designed a system that allows to obtain a wide spectrum of potential characteristics, from a single well to a three-well system, and we have analyzed its effectiveness. Next, we checked the influence of parameters characterizing the change of potential using bifurcation diagrams and their comparison with the effective voltage values. We also analyzed the behavior of the system in chaotic and periodic motion zones and presented selected sections of Poincare and Fourier amplitude-frequency spectra of chaotic solutions. The last element of the analysis was the impact of cyclic potential change on coexisting solutions. We have shown that the best effectiveness is achieved when the frequency of the external load is equal to the resonant frequency of the flexible cantilever beam and the change of potential is limited to extreme positions.

Tunable subdiffusion in the Caputo fractional standard map

The Caputo fractional standard map (C-fSM) is a two-dimensional nonlinear map with memory given in action–angle variables (I,θ). It is parameterized by K and α∈(1,2] which control the strength of nonlinearity and the fractional order of the Caputo derivative, respectively. In this work we perform a scaling study of the average squared action I2 along strongly chaotic orbits, i.e. when K≫1. We numerically prove that I2∝nμ with 0≤μ(α)≤1, for large enough discrete times n. That is, we demonstrate that the C-fSM displays subdiffusion for 1<α<2. Specifically, we show that diffusion is suppressed for α→1 since μ(1)=0, while standard diffusion is recovered for α=2 where μ(2)=1. We describe our numerical results with a phenomenological analytical estimation.

Energy flow investigations of Rayleigh-Plesset equation for cavitation simulations

Nonlinear Rayleigh-Plesset Equation (RPE) for cavitation simulations is investigated using energy-flow theory. Nondimensional RPE affected by Reynolds-number, surface-tension, bubble-pressure, and liquid static-dynamic pressure ratio is derived to examine its equilibrium points with stability, possible periodical/chaotic motions by the energy-flow criteria. An example is numerically analyzed to illustrate the developed method. Five cases of the example reveal: 1) It is a damped system, where initial disturbances are gradually reduced with its phase point tends to its stable equilibrium-point; 2) Reynolds-number affects the damping of system, large one corresponds small damping; 3) Bubble-pressure, surface-tension and liquid static-dynamic pressure ratio affect the position of equilibrium point; 4) Periodical orbit appears only in forced vibrations, in which free vibration is reduced with time forward and the system finally shows a stable periodical oscillation; 5) Energy flow criteria for chaotic motions is not reached, and there are no chaotic motions for the cases of example. Numerical simulations confirm the developed energy-flow means with available computer code is effective to investigate generalized RPEs to reveal their inherent characteristics affected by required parameters in engineering cavitation analysis and designs, such as considering mass transports across the boundary of bubble by evaporation or condensation of liquids.

Nonlinear behaviors and chaos control of a fluid-conveying piezoelectric pipe in supercritical regime

The study investigated the steady state response of a piezoelectric pipe conveying fluid under external and internal resonances in the supercritical regime of fluid velocity. Internal resonance conditions between the first and second frequency modes were studied at a specific flow velocity in the supercritical regime. The Galerkin method and the multiple scale method were used to extract the vibration amplitude versus the excitation frequency and amplitude. The steady state solutions lost stability through saddle-node and Hopf bifurcations, leading to periodic, double periodic, multi-periodic, and chaotic behaviors in the forced and frequency response curves. Time response, FFT, phase portrait, and Lyapunov exponents were presented to predict the unstable regions of the piezo-pipe system between the Hopf points. The Wolf algorithm was utilized to evaluate the Lyapunov exponents. Additionally, a fuzzy terminal sliding mode controller was designed in the chaotic region of the piezoelectric pipe conveying fluid, and its parameters were optimized using the genetic algorithm. The controller effectively stabilized the chaotic motion of the system, demonstrating better performance than the sliding mode control in numerical results.

Dynamical analysis of a stay cable with a nonlinear energy sink and time-delayed feedback control

Stay cables exist commonly in engineering and can exhibit sustainable and large vibrations under external excitations. The unexpected oscillations always lead to the fatigue and damage of structures, resulting in the economical loss. This paper studies the vibration control of a stay cable system with a nonlinear vibration absorber and time delay control. The governing equations of motion of the stay cable system under harmonic excitations are given. The stability and different bifurcations of the trivial equilibrium point of the system are analyzed. The bifurcation diagrams with respect to the control parameters are illustrated. The frequency responses of the forced system are obtained by using the harmonic balance method. Numerical simulations are performed to validate the obtained results and various dynamical phenomena are shown, such as the multi-period oscillations, quasi-periodic responses and chaotic motions. The effects of the time-delayed feedback control on the nonlinear dynamics and vibration reduction of the system are extracted. It is shown that the proper values of the time delay and feedback gain can improve the vibration reduction performance of the system, such as the suppression of the resonance peaks and elimination of the complex responses.

Mixed (regular and chaotic) dynamics under the channeling of the high energy positrons in [100] direction of the silicon crystal and Podolskiy-Narimanov distribution

The character of motion (regular or chaotic) on the quantum level manifests itself in the statistics of the set of the system's energy levels. The completely regular case generates the sequence of the levels with exponential (Poisson) level spacing distribution; the completely chaotic one — with Wigner distribution. The most interesting case is the co-existence of the regular and chaotic motion domains in the phase space of the system under consideration. The assumption of independent generation of two level sequencies by one chaotic domain and by all regular domains leads to Berry-Robnik distribution. However, the presence of the chaotic motion domain in the phase space affects the levels produced by the regular domains via the so-called chaos-assisted tunneling (CAT) that leads to Podolskiy-Narimanov distribution of the level spacings. This distribution needs the mean amplitude of the tunnel transition and the relative contribution of the regular domains to the mean level density as the parameters. Using their estimations for the transverse motion of the high energy positrons channeling near [100] direction in the silicon crystal we found that Podolskiy-Narimanov distribution demonstrates the best agreement with the level spacing distribution.

Orbital analysis of stars in the nuclear stellar disc of the Milky Way

Context. While orbital analysis studies were so far mainly focused on the Galactic halo, it is possible now to do these studies in the heavily obscured region close to the Galactic Centre. Aims. We aim to do a detailed orbital analysis of stars located in the nuclear stellar disc (NSD) of the Milky Way allowing us to trace the dynamical history of this structure. Methods. We integrated orbits of the observed stars in a non-axisymmetric potential. We used a Fourier transform to estimate the orbital frequencies. We compared two orbital classifications, one made by eye and the other with an algorithm, in order to identify the main orbital families. We also compared the Lyapunov and the frequency drift techniques to estimate the chaoticity of the orbits. Results. We identified several orbital families as chaotic, z-tube, x-tube, banana, fish, saucer, pretzel, 5:4, and 5:6 orbits. As expected for stars located in a NSD, the large majority of orbits are identified as z-tubes (or as a sub-family of z-tubes). Since the latter are parented by x2 orbits, this result supports the contribution of the bar (in which x2 orbits are dominant in the inner region) in the formation of the NSD. Moreover, most of the chaotic orbits are found to be contaminants from the bar or bulge which would confirm the predicted contamination from the most recent NSD models. Conclusions. Based on a detailed orbital analysis, we were able to classify orbits into various families, most of which are parented by x2-type orbits, which are dominant in the inner part of the bar.

The Evolution of the Phase Space Structure Along Pitchfork and Period-Doubling Bifurcations in a 3D-Galactic Bar Potential

We investigate how the phase space structure of a Three-Dimensional (3D) autonomous Hamiltonian system evolves across a series of successive Two-Dimensional (2D) and 3D pitchfork and period-doubling bifurcations, as the transition of the parent families of Periodic Orbits (POs) from stability to simple instability leads to the creation of new stable POs. Our research illustrates the consecutive alterations in the phase space structure near POs as the stability of the main family of POs changes. This process gives rise to new families of POs within the system, either maintaining the same or exhibiting higher multiplicity compared to their parent families. Tracking such a phase space transformation is challenging in a 3D system. By utilizing the color and rotation technique to visualize the Four-Dimensional (4D) Poincaré surfaces of section of the system, i.e. projecting them onto a 3D subspace and employing color to represent the fourth dimension, we can identify distinct structural patterns. Perturbations of parent and bifurcating stable POs result in the creation of tori characterized by a smooth color variation on their surface. Furthermore, perturbations of simple unstable parent POs beyond the bifurcation point, which lead to the birth of new stable families of POs, result in the formation of figure-8 structures of smooth color variations. These figure-8 formations surround well-shaped tori around the bifurcated stable POs, losing their well-defined forms for energies further away from the bifurcation point. We also observe that even slight perturbations of highly unstable POs create a cloud of mixed color points, which rapidly move away from the location of the PO. Our study introduces, for the first time, a systematic visualization of 4D surfaces of section within the vicinity of higher multiplicity POs. It elucidates how, in these cases, the coexistence of regular and chaotic orbits contributes to shaping the phase space landscape.

Nonlinear combination resonance analysis of parametric-forced excitation for an axially moving piezoelectric rectangular thin plate in thermal- electromechanical coupling field

In this paper, the combination resonance of parametric-forced excitation characteristics for an axially moving rectangular piezoelectric plate under a thermal-electromechanical field is studied. Based on Kirchhoff-Love plate theory and Von Karman theory, the transverse vibration governing equations are derived from Hamilton's principle. The equations are discretized to ordinary differential equations by the Galerkin method. Then, the multiple scales method is applied to solve the system combination resonance equation, two different resonance states and corresponding amplitude-frequency response equations are obtained by eliminating the secular term, respectively. Additionally, the stability of the steady-state responses is analyzed by Lyapunov stability. Based on the numerical analysis, the influence of axial velocity, external voltage, central temperature difference, structural damping, and other parameters on nonlinear combination resonance response are investigated. The effect of parameter variation on period-doubling bifurcation and the chaotic motion of the system are also discussed by the system global bifurcation diagram.

Nonlinear dynamics analysis of electromechanical planetary gear systems with temperature effects considered

PurposeThe electromechanical planetary transmission system has the advantages of high transmission power and fast running speed, which is one of the important development directions in the future. However, during the operation of the electromechanical planetary transmission system, friction and other factors will lead to an increase in gear temperature and thermal deformation, which will affect the transmission performance of the system, and it is of great significance to study the influence of the temperature effect on the nonlinear dynamics of the electromechanical planetary system.Design/methodology/approachThe effects of temperature change, motor speed, time-varying meshing stiffness, meshing damping ratio and error amplitude on the nonlinear dynamic characteristics of electromechanical planetary systems are studied by using bifurcation diagrams, time-domain diagrams, phase diagrams, Poincaré cross-sectional diagrams, spectra, etc.FindingsThe results show that when the temperature rise is less than 70 °C, the system will exhibit chaotic motion. When the motor speed is greater than 900r/min, the system enters a chaotic state. The changes in time-varying meshing stiffness, meshing damping ratio, and error amplitude will also make the system exhibit abundant bifurcation characteristics.Originality/valueBased on the principle of thermal deformation, taking into account the temperature effect and nonlinear parameters, including time-varying meshing stiffness and tooth side clearance as well as comprehensive errors, a dynamic model of the electromechanical planetary gear system was established.

Voting model prediction of nonlinear behavior for double-circumferential-slot air bearing system

Double-circumferential-slot air bearing (DCSAB) systems provide multidirectional supporting forces and have high stiffness, increasing the stability of instruments at high rotational speeds. However, DCSAB systems may exhibit chaotic motion because of a nonlinear pressure distribution within the gas film, supplied gas imbalances, or an inappropriate design. This study investigated the occurrence of nonperiodic motion in a DCSAB system by analyzing the dynamic response of systems with different rotor masses and bearing numbers. The dynamic trajectory, spectral response, bifurcation, Poincaré map, and maximum Lyapunov exponent were analyzed to identify chaotic behavior. Behavior was found to be highly sensitive to rotor mass and bearing number; the system exhibits chaotic behavior when the rotor mass has values in three intervals within 0.1–6.0 kg given a fixed bearing number of Λ = 3.8. To reduce the computational cost of predicting chaotic behavior, the maximum Lyapunov exponent was predicted using various machine learning models; a voting model combining random forest with XGBoost has the highest performance. The results can be used as a guideline for designing of DCASB systems for use in industrial applications.

Electromagnetic field and chaotic charged-particle motion around hairy black holes in Horndeski gravity

The Wald vector potential is an exact solution of the source-less Maxwell equations regarding an electromagnetic field of a vacuum uncharged black hole like the Kerr background black hole in an asymptotically uniform magnetic field. However, it is not if the black hole is a nonvacuum solution in a theory of modified gravity with extra fields or a charged Kerr–Newman spacetime. To satisfy the source-less Maxwell equations in this case, the Wald vector potential must be modified and generalized appropriately. Following this idea, we derive an expression for the vector potential of an electromagnetic field surrounding a hairy black hole in the Horndeski modified gravity theory. Explicit symplectic integrators with excellent long-term behaviour are used to simulate the motion of charged particles around the hairy black hole immersed in the external magnetic field. The recurrence plot method based on the recurrence quantification analysis uses diagonal structures parallel to the main diagonal to show regular dynamics, but adopts no diagonal structures to indicate chaotic dynamics. The method is efficient to detect chaos from order in the curved spacetime, as the Poincaré map and the fast Lyapunov indicator are.

Experimental Research on an Afterburner System Fueled with Hydrogen–Methane Mixtures

A new afterburner installation is proposed, fueled with pure hydrogen (100%H2) or hydrogen–methane mixtures (60% H2 + 40% CH4, 80% H2 + 20% CH4) for use in cogeneration applications. Two prototypes (P1 and P2) with the same expansion angle (45 degrees) were developed and tested. P1 was manufactured by the classic method and P2 by additive manufacturing. Both prototypes were manufactured from Inconel 625. During the tests, analysis of flue gas (CO2, CO, and NO concentration), PIV measurements, and noise measurements were conducted. The flue gas analysis emphasizes that the behavior of the two tested prototypes was very similar. For all three fuels used, the CO2 concentration levels were slightly lower in the case of the additive-manufactured prototype P2. The CO concentration levels were significantly higher in the case of the additive-manufactured prototype P2 when 60% H2/40% CH4 and 80% H2/20% CH4 mixtures were used as fuel. When pure H2 was used as fuel, the measured data suggest that no additional CO was produced during the combustion process, and the level of CO was similar to that from the Garrett micro gas turbine in all five measuring points. The NO emissions gradually decreased as the percentage of H2 in the fuel mixture increased. The NO concentration was significantly lower in the case of the additive-manufactured prototype (P2) in comparison with the classic manufactured prototype (P1). Examining the data obtained from the PIV measurements of the flow within the mixing region shows that the highest axial velocity component value on the centerline was measured for the P1 prototype. The acoustic measurements showed that a higher H2 concentration led to a reduction in noise of approximately 1.5 dB for both afterburner prototypes. The outcomes reveal that the examined V-gutter flame holder prototype flow was smooth, without any perpendicular oscillations, without chaotic motions or turbulent oscillations to the flow direction, across all tested conditions, keeping constant thermal power.

Appearance of transient chaos and energy variability in Ueda oscillator driven by an amplitude‐modulated force

Transient chaos is a characteristic behavior observed in nonlinear systems, where trajectories within a specific region of phase space exhibit chaotic dynamics for a finite duration before transitioning to an external attractor. In many nonlinear systems, transient and chaotic behaviors are closely associated with variations in the system's energy. In this paper, we introduce the concept of an energy variable and establish its connection with the Melnikov integral. We explore the influence of energy variation in an amplitude‐modulated (AM) force‐driven Ueda oscillator through numerical simulations. Our investigation reveals the emergence of transient chaos, chaotic dynamics, and regular behaviors, underscoring the significant role played by the energy variable, denoted as , in the system. We employ various analytical tools, including phase diagrams, Poincaré maps, time series plots, and bifurcation diagrams, to characterize and visualize transient chaos, regular, and chaotic motions within the system.

Dynamical Structures under Nonrestricted Hierarchical Planetary Systems with Different Mass Ratios

Secular dynamics have been extensively studied in both the inner and outer restricted hierarchical three-body systems. In the inner restricted problem, the quadrupole-order resonance (i.e., the well-known Kozai resonance) causes large coupled oscillations of eccentricity and inclination when the maximum inclination is higher than 39.2°, and the octupole-order resonance leads to the behavior of orbital flips. In the outer restricted problem, the behavior of orbital flips is due to the quadrupole-order resonance. Secular dynamics under the inner and outer restricted systems are distinctly different. The mass ratio of inner and outer bodies could change the ratio of circular orbital angular momenta β, which significantly influences dynamical structures of the system. But this influence is still unclear. In this paper, we focus on nonrestricted hierarchical planetary systems where β > 1 and investigate the secular dynamics by changing mass ratios. Dynamical structures are systematically explored from four aspects: periodic orbits, secular resonances, orbital flips, and chaos detection. We find that (a) it tends to lead to more bifurcations in the host family of prograde periodic orbits associated with Kozai resonance with smaller β; (b) with the decrease of β, fewer orbits inside the octupole-order resonance can realize flip; (c) for given initial conditions, the forbidden region appears in the retrograde region and becomes larger as β decreases, meaning that the mutual inclination cannot reach a very high value if β is small; and (d) chaotic orbits are distributed in the low-eccentricity, high-inclination region when β > 1.

Chaotic dynamics of a three-species food chain model

An analysis of a discrete-time food chain model is presented in this paper to investigate the effects of strong prey pressure on the chain. A bifurcation and stability theory is used to analyze equilibrium formation in system. It is shown by numerical simulations that chaos occurs as a result of the bifurcation of invariant curves following Neimark-Sacker (NS) bifurcation rules. The OGY method on stable periodic orbits of period 1 further controls chaos. Stability is achieved through iterations for chaotic motion controlled by different regulator poles. Numerical simulations are performed to examine the effects of the control method along with the results of our theoretical analysis.

Dynamische Systeme

This workshop continues a series of workshops whose current format originated in 1981 under then-organizers Moser and Zehnder, and whose latest iteration took place in July 2023. The general goal of this series of workshops is to discuss the latest developments in the field of dynamical systems, broadly construed, and its connections with neighboring areas of mathematics such as differential geometry, partial differential equations, and more recently contact and symplectic geometry. We continued this tradition, bringing in new participants working in areas of dynamical systems and its connections with other areas of mathematics that are currently highly active and/or showing great promise for future development. Key focus areas for the 2023 workshop include spectral rigidity for planar domains, chaotic and oscillatory motions in celestial mechanics, conformal symplectic dynamics, and relations between dynamics.he workshop by the grant DMS-2230648, “US Junior Oberwolfach Fellows”.

Resonance response and chaotic analysis for an irrational pendulum system

We propose a novel rotating pendulum system with irrational nonlinearity, bi-stable phenomenon and quasi-zero stiffness characteristic, which bears a transition from standard double well dynamics to discontinuous double well dynamics by adjusting a geometric parameter. Nonlinear dynamical behaviors of the proposed irrational pendulum system subjected to viscus damping and periodic excitation are studied. To analytically demonstrate the primary resonance response for this perturbed irrational pendulum system, we introduce an approximate irrational system which is significant parallels with the original pendulum system from the perspective of qualitative analysis and quantitative calculation. Averaging method is applied to investigate the dynamic response of the perturbed irrational pendulum system with quasi-zero stiffness at origin. The effect of the internal and external parameters on the response curves are discussed. Numerical fitting technique and semi-analytical Melnikov method are used to detect the chaotic boundaries of the perturbed irrational pendulum system with two types of homoclinic orbits. Different types of periodic solutions and chaos are clarified in the perturbed irrational pendulum system by using numerical simulations. It is found that a class of chaos with oscillatory and rotational motions bears significant similarities to paroxysmal chaos due to sudden rotational motion and two chaotic motions formed by period-doubling bifurcation merge into one chaos.