Periodic Orbits Research Articles

B–plane assisted trajectory design with analytical mapping of swing-by parameters

The Artemis program has gathered further attention due to the transfer trajectories to the Moon and beyond. Trajectories such as weak stability boundary and ballistic lunar transfer are known as low-energy transfer toward the vicinity of the Moon. Among these categories, using the Moon swing-by contributes to the reduction of fuel consumption. In this study, we investigate trajectory design with the Moon swing-by as a constraint, and propose a methodology that projects swing-by parameters onto a conventional B–plane to facilitate the efficient analysis and design of post swing-by trajectories. This proposed approach not only enables a comprehensive understanding of all candidate solutions at each swing-by, but also provides flexibility in response to varied swing-by conditions. As illustrative examples, this paper demonstrates transfer problems to periodic orbits in the Earth–Moon system.

Modeling the Nonmonotonic Immune Response in a Tumor–Immune System Interaction

Tumor–immune system interactions are very complicated, being highly nonlinear and not well understood. A large number of tumors can potentially weaken the immune system through various mechanisms such as secreting cytokines that suppress the immune response. In this paper, we propose a tumor–immune system interaction model with a nonmonotonic immune response function and adoptive cellular immunotherapy (ACI). The model has a tumor-free equilibrium and at most three tumor-presence equilibria (low, moderate and high ones). The stability of all equilibria is studied by analyzing their characteristic equations. The consideration of nonmonotonic immune response results in a series of bifurcations such as the saddle-node bifurcation, transcritical bifurcation, Hopf bifurcation and Bogdanov–Takens bifurcation. In addition, numerical simulation results show the coexistence of periodic orbits and homoclinic orbits. Interestingly, along with various bifurcations, we also found two bistable scenarios: the coexistence of a stable tumor-free as well as a high-tumor-presence equilibrium and the coexistence of a stable-low as well as a high-tumor-presence equilibrium, which can show symmetric and antisymmetric properties in a range of model parameters and initial cell concentrations. The new findings indicate that under ACI, patients can possibly reach either a stable tumor-free state or a low-tumor-presence state in the presence of nonmonotonic immune response once the immune system is activated.

Recurrent flow patterns as a basis for two-dimensional turbulence: Predicting statistics from structures

A dynamical systems approach to turbulence envisions the flow as a trajectory through a high-dimensional state space [Hopf, Commun. Appl. Maths 1, 303 (1948)]. The chaotic dynamics are shaped by the unstable simple invariant solutions populating the inertial manifold. The hope has been to turn this picture into a predictive framework where the statistics of the flow follow from a weighted sum of the statistics of each simple invariant solution. Two outstanding obstacles have prevented this goal from being achieved: 1) paucity of known solutions and 2) the lack of a rational theory for predicting the required weights. Here, we describe a method to substantially solve these problems, and thereby provide compelling evidence that the probability density functions (PDFs) of a fully developed turbulent flow can be reconstructed with a set of unstable periodic orbits. Our method for finding solutions uses automatic differentiation, with high-quality guesses constructed by minimizing a trajectory-dependent loss function. We use this approach to find hundreds of solutions in turbulent, two-dimensional Kolmogorov flow. Robust statistical predictions are then computed by learning weights after converting a turbulent trajectory into a Markov chain for which the states are individual solutions, and the nearest solution to a given snapshot is determined using a deep convolutional autoencoder. In this study, the PDFs of a spatiotemporally chaotic system have been successfully reproduced with a set of simple invariant states, and we provide a fascinating connection between self-sustaining dynamical processes and the more well-known statistical properties of turbulence.

Periodic orbits for 2n-dimensional control piecewise smooth dynamical systems

Periodic orbits for 2n-dimensional control piecewise smooth dynamical systems

Dynamics of a Cournot game with bounded rational firms and various scale effects.

In this paper, we focus on a dynamic Cournot game in the market with a nonlinear (isoelastic) demand function. In our model, there are N competing firms featured by nonlinear cost functions, which enhances our model's resemblance to real-world scenarios. Firstly, we focus on the homogeneous case where firms' marginal costs change at equal rates. We establish analytical expressions of the market supply at equilibrium and perform comparative static analysis. In addition, we investigate the local stability under different economies of scale and show that there could be multiple stable equilibria if firms face economies of scale. The heterogeneous case where firms' marginal costs change at distinct rates is much more complex, thus we investigate the duopoly game with only two firms involved. Methods of symbolic computations such as triangular decomposition and partial cylindrical algebraic decomposition are employed in the analytical investigations of the equilibrium, which is nearly impossible by using the pencil-and-paper approach since the closed-form equilibrium is quite complicated. According to the computational results, we derive that two stable positive equilibria may coexist if both firms face economies of scale. Additionally, we conduct preliminary numerical simulations and find two different types of complex dynamics of the model considered in this paper: complex trajectories such as periodic and chaotic orbits may appear; the topological structure of the basins of attraction may be complex.

Data-driven model reduction for pipes conveying fluid via spectral submanifolds

Pipes conveying fluid have non-linearizable dynamics when the flow velocity is large enough. Developing analytical models that accurately capture these highly nonlinear and complex dynamics is challenging. Efficient analysis and predictions of these nonlinear dynamics are essential in the design and control of pipe systems in engineering applications. To address these challenges, here we construct low-dimensional reduced-order models in a data-driven setting. In particular, we perform model reduction on spectral submanifolds (SSMs) via SSMLearn. This open-source package automates the identification of low-dimensional SSMs and their associated reduced-order models (ROMs). Our reduction via SSMs treats pipes with various boundary conditions and flow velocities in a unified way. It achieves a significant reduction because the ROMs are two-dimensional, which enables efficient and even analytical predictions on the nonlinear vibration of the pipes, paving the way for real-time prediction of the non-linearizable dynamics. In addition, our SSM-based reduction shows remarkable extrapolation capability because it is built on invariant manifolds and normal form theory. We demonstrate that the ROMs trained on unforced responses can be directly extrapolated to make reliable predictions on forced vibrations under various excitation frequencies and amplitudes, including periodic and quasi-periodic orbits and their bifurcations.

Distribution in homology classes and discrete fractal dimension

In this note, we examine the proportion of periodic orbits of Anosov flows that lie in an infinite zero density subset of the first homology group. We show that on a logarithmic scale we get convergence to a discrete fractal dimension.

Precessing and periodic timelike orbits and their potential applications in Einsteinian cubic gravity

Einsteinian cubic gravity (ECG) is the most general theory up to cubic order in curvature, which has the same graviton spectrum as the Einstein theory. In this paper, we investigate the geodesic motions of timelike particles around the four dimensional asymptotically flat black holes in ECG, and discuss their potential applications when connecting them with recent observational results. We first explore the effects of the cubic couplings on the marginally bound orbits (MBO), innermost stable circular orbits (ISCO) and on the periodic orbits around the Einsteinian cubic black hole. We find that comparing to Schwarzschild black hole in general relativity, the cubic coupling enhances the energy as well as the angular momentum for all the bound orbits of the particles. Then, we derive the relativistic periastron precessions of the particles and give a preliminary bound on the cubic coupling employing the observational result of the S2 star’s pericenter precession in SgrA*. Finally, after calculating the periodic orbits’ configurations, we preliminarily evaluate the gravitational waveform radiated from several periodic orbits in one complete period of a test object which orbits a supermassive Einsteinian cubic black hole. Our studies could be helpful for us to better understand the gravitational structure of the theory with high curvatures.

Studying the Upper Bounds of the Numbers of Zeros of Abelian Integrals by the Law of Polynomials

For the quadratic reversible systems of genus one, all of their periodic orbits are higher-order algebraic curves. When they are perturbed by polynomials of degree [Formula: see text], the numbers of zeros of their Abelian integrals will change and we study the upper bounds of these numbers by using the methods of Riccati equation and Picard–Fuchs equation. We consider both the highest and lowest degrees of polynomials, and more importantly, we consider the law of polynomials and the range of values for their variables. Consequently, some laws of the polynomials are discovered and many upper bounds are obtained, and these upper bounds are sharper than the results obtained by other techniques.

On non-contractible periodic orbits and bounded deviations

We present a dichotomy for surface homeomorphisms in the isotopy class of the identity. We show that, in the absence of a degenerate fixed point set, either there exists a uniform bound on the diameter of orbits of non-wandering points for the lifted dynamics in the universal covering space, or the map has non-contractible periodic orbits. We then use this new tool to characterize the dynamics of area preserving homeomorphisms of the torus without non-contractible periodic orbits, showing that if the fixed point set is non-degenerate, then either the lifted dynamics is uniformly bounded, or it has a single strong irrational dynamical direction.

Periodic orbits in a galactic potential

We make a numerical study of the periodic orbits of period-1 and their bifurcations in a galactic type potential V=12[α(x2+y2)+x2y2]\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$V=\\frac{1}{2}[\\alpha (x^2+y^2)+x^2y^2]$$\\end{document}, which tends to the Yang–Mills potential V=12x2y2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$V=\\frac{1}{2}x^2y^2$$\\end{document}, when α\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\alpha $$\\end{document} tends to zero. We consider their stability diagrams, the corresponding Poincaré surfaces of section and the forms of the orbits.

Mathematical Analysis for the Influence of Seasonality on Chikungunya Virus Dynamics

In this article, we discuss a mathematical system modelling Chikungunya virus dynamics in a seasonal environment with general incidence rates. We establish the existence, uniqueness, positivity and boundedness of a periodic orbit. We show that the global dynamics is determined using the basic reproduction number denoted by R0 and calculated using the spectral radius of a linear integral operator. We show the global stability of the disease free periodic solution if R0<1 and we show also the persistence of the disease if R0>1 where the trajectories converge to a periodic orbit. Finally, we display some numerical examples confirming the theoretical findings.

Exploration and Maintenance of Homeomorphic Orbit Revs in the Elliptic Restricted Three-Body Problem

A novel station-keeping strategy leveraging periodic revolutions of homeomorphic orbits in the Elliptic Restricted Three-Body Problem within the pulsating frame is presented. A systemic approach founded on arc-length continuation is presented for the discovery, computation, and classification of periodic revolutions that morph from their traditional circular restricted three-body counterparts to build an a priori dataset. The dataset is comprehensive in covering all possible geometric architectures of the restricted problem. Shape similarity is quantified using Hausdorff distance and works as a filter for the station-keeping algorithm in relation to appropriate target conditions. Finally, an efficient scheme to quantify impulsive orbit maintenance maneuvers that minimize the total fuel cost is presented. The proposed approach is salient in its generic applicability across any elliptic three-body system and any periodic orbit family. Finally, average annual station-keeping costs using the described methodology are quantified for selected “orbits of interest” in the cis-lunar and the Sun–Earth systems. The robustness and efficacy of the approach instill confidence in its applicability for realistic mission design scenarios.

Cupolets: History, Theory, and Applications

In chaos control, one usually seeks to stabilize the unstable periodic orbits (UPOs) that densely inhabit the attractors of many chaotic dynamical systems. These orbits collectively play a significant role in determining the dynamics and properties of chaotic systems and are said to form the skeleton of the associated attractors. While UPOs are insightful tools for analysis, they are naturally unstable and, as such, are difficult to find and computationally expensive to stabilize. An alternative to using UPOs is to approximate them using cupolets. Cupolets, a name derived from chaotic, unstable, periodic, orbit-lets, are a relatively new class of waveforms that represent highly accurate approximations to the UPOs of chaotic systems, but which are generated via a particular control scheme that applies tiny perturbations along Poincaré sections. Originally discovered in an application of secure chaotic communications, cupolets have since gone on to play pivotal roles in a number of theoretical and practical applications. These developments include using cupolets as wavelets for image compression, targeting in dynamical systems, a chaotic analog to quantum entanglement, an abstract reducibility classification, a basis for audio and video compression, and, most recently, their detection in a chaotic neuron model. This review will detail the historical development of cupolets, how they are generated, and their successful integration into theoretical and computational science and will also identify some unanswered questions and future directions for this work.

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.

Critical intermittency in rational maps

Intermittent dynamics is characterized by long periods of different types of dynamical characteristics, for instance almost periodic dynamics alternated by chaotic dynamics. Critical intermittency is intermittent dynamics that can occur in iterated function systems, and involves a superattracting periodic orbit. This paper will provide and study examples of iterated function systems by two rational maps on the Riemann sphere that give rise to critical intermittency. The main ingredient for this is a superattracting fixed point for one map that is mapped onto a common repelling fixed point by the other map. We include a study of topological properties such as topological transitivity.

On a theorem by Schlenk

In this paper we prove a generalisation of Schlenk’s theorem about the existence of contractible periodic Reeb orbits on stable, displaceable hypersurfaces in symplectically aspherical, geometrically bounded, symplectic manifolds, to a forcing result for contractible twisted periodic Reeb orbits. We make use of holomorphic curve techniques for a suitable generalisation of the Rabinowitz action functional in the stable case in order to prove the forcing result. As in Schlenk’s theorem, we derive a lower bound for the displacement energy of the displaceable hypersurface in terms of the action value of such periodic orbits. The main application is a forcing result for noncontractible periodic Reeb orbits on quotients of certain symmetric star-shaped hypersurfaces. In this case, the lower bound for the displacement energy is explicitly given by the difference of the two periods. This theorem can be applied to many physical systems including the Hénon–Heiles Hamiltonian and Stark–Zeeman systems. Further applications include a new proof of the well-known fact that the displacement energy is a relative symplectic capacity on R2n\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${\\mathbb {R}}^{2n}$$\\end{document} and that the Hofer metric is indeed a metric.

Dynamics of bouncing convex body

We consider a dynamical system that describes a convex planer body with a smooth boundary, which falls under the influence of gravity and is elastically reflected upon hitting a wall. Under suitable restriction on the system’s parameters, the system is equivalent to a billiard problem within a cylinder. We prove that the dynamics of the system is intimately linked to a variational problem, corresponding to a twist map. Moreover, we utilise the established variational method to investigate the stability of trivial periodic orbits and near-integrable dynamics when this planar body closely approximates a disc with its centre of mass at the geometric centre. Numerical results illustrates our mathematical findings and highlight the system’s intricate dynamics.

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.

Comparison of Homoclinic Bifurcations Between Grid-Following and Grid-Forming Converters

A second-order model is developed to describe the essential behavior of the grid-following converter (GFLC) and the grid-forming converter (GFMC) for studying the synchronization characteristics of grid-connected converter systems. In this paper, a general set of criteria in terms of the manifolds and saddle quantity is derived for studying the homoclinic bifurcation behavior of GFLCs and GFMCs. It is shown that three unstable periodic orbits and two homoclinic bifurcations may exist for the GFLC, while only one stable periodic orbit and one homoclinic bifurcation may exist for the GFMC. It is also found that after the onset of the first homoclinic bifurcation, the stable equilibrium point (SEP) is surrounded by the stable manifolds of the unstable equilibrium point (UEP), which forms the SEP's basin of attraction, for both GFLCs and GFMCs. In this case, the converters may lose their synchronization if the trajectory crosses over the UEP. Therefore, system parameters should be designed to avoid the onset of the aforementioned bifurcation so that the SEP's basin of attraction can be significantly enlarged to contain the UEP. As a result, the GFLC or GFMC can eventually resynchronize with the grid even when its trajectory passes the UEP. Finally, experimental results are provided to verify these theoretical findings.