Codimension-two bifurcations induce hysteresis behavior and multistabilities in delay-coupled Kuramoto oscillators

Hysteresis phenomena and multistability play crucial roles in the dynamics of coupled oscillators, which are now interpreted from the point of view of codimension-two bifurcations. On the Ott–Antonsen’s manifold, two-parameter bifurcation sets of delay-coupled Kuramoto model are derived regarding coupling strength and delay as bifurcation parameters. It is rigorously proved that the system must undergo Bautin bifurcations for some critical values; thus, there always exists saddle-node bifurcation of periodic solutions inducing hysteresis loop. With the aid of center manifold reduction method and the MATLAB package DDE-BIFTOOL, the location of Bautin and double Hopf points and detailed dynamics are theoretically determined. We find that, near these critical points, four coherent states (two of which are stable) and a stable incoherent state may coexist and that the system undergoes Neimark–Sacker bifurcation of periodic solutions. Finally, the clear scenarios about the synchronous transition in delayed Kuramoto model are depicted.     

Constructing Galerkin's approximations of invariant tori using MACSYMA

Invariant tori of solutions for nonlinearly coupled oscillators are generalizations of limit cycles in the phase plane. They are surfaces of aperiodic solutions of the coupled oscillators with the property that once a solution is on the surface it remains on the surface. Invariant tori satisfy a defining system of nonlinear partial differential equations. This case study shows that with the help of a symbolic manipulation package, such as MACSYMA, approximations to the invariant tori can be developed by using Galerkin's variational method. The resulting series must be manipulated efficiently, however, by using the Poisson series representation for multiply periodic functions, which makes maximum use of the list processing techniques of MACSYMA. Three cases are studied for the single van der Pol oscillator with forcing parameter =0.5, 1.0, 1.5, and three cases are studied for a pair of nonlinearly coupled van der Pol oscillators with forcing parameters =0.005, 0.5, 1.0. The approximate tori exhibit good agreement with direct numerical integrations of the systems.Contribution of the National Institute of Standards and Technology, a Federal agency. Not subject to copyright.     

Dynamics of Nonlinear Cyclic Systems with Structural Irregularity

The dynamics of weakly coupled, nonlinear cyclic assemblies are investigated in the presence of weak structural mistuning. The method of multiple scales is used to obtain a set of nonlinear algebraic equations which govern the steady-state, synchronous (modal-like) motions for the structures. Considering a degenerate assembly of uncoupled oscillators, spatially localized modes are computed corresponding to motions during which vibrational energy is spatially confined to one oscillator (strong localization) or a subset of oscillators (weak localization). In the limit of weak substructural coupling, asymptotic solutions are obtained which correspond to (i) spatially extended, (ii) strongly localized, and (iii) weakly localized modes for fully coupled systems. Throughout the analysis, the influence of structural mistunings on the resulting solutions are discussed. Additionally, numerical solutions (including linearized stability characteristics) are obtained for prototypical two- and three-degree-of-freedom (DoF) systems with various structural mistunings. The numerical results provide insight into the strong influence of structural irregularities on the dynamical behavior of nonlinear cyclic systems, and demonstrate that the combined influences of structural mistunings and nonlinearities do not lead to uniform improvement of motion confinement characteristics.     

Knotted periodic orbits in suspensions of Smale's horseshoe: Torus knots and bifurcation sequences

We construct a suspension of Smale's horseshoe diffeomorphism of the two-dimensional disc as a flow in an orientable three manifold. Such a suspension is natural in the sense that it occurs frequently in periodically forced nonlinear oscillators such as the Duffing equation. From this suspension we construct a knot-hòlder or template—a branched two-manifold with a semiflow—in such a way that the periodic orbits are isotopic to those in the full three-dimensional flow. We discuss some of the families of knotted periodic orbits carried by this template. In particular we obtain theorems of existence, uniqueness and non-existence for families of torus knots. We relate these families to resonant Hamiltonian bifurcations which occur as horseshoes are created in a one-parameter family of area preserving maps, and we also relate them to bifurcations of families of one-dimensional quadratic like maps which can be studied by kneading theory. Thus, using knot theory, kneading theory and Hamiltonian bifurcation theory, we are able to connect a countable subsequence of one-dimensional bifurcations with a subsequence of area-preserving bifurcations in a two parameter family of suspensions in which horseshoes are created as the parameters vary. One implication is that infinitely many bifurcation sequences are reversed as one passes from the one dimensional to the area-preserving family: there are no universal routes to chaos!     

Regular oscillations,chaos, and multistability in a system of two coupled van der Pol oscillators: numerical and experimental studies

In this paper, the dynamics of a system of two coupled van der Pol oscillators is investigated. The coupling between the two oscillators consists of adding to each one’s amplitude a perturbation proportional to the other one. The coupling between two laser oscillators and the coupling between two vacuum tube oscillators are examples of physical/experimental systems related to the model considered in this paper. The stability of fixed points and the symmetries of the model equations are discussed. The bifurcations structures of the system are analyzed with particular attention on the effects of frequency detuning between the two oscillators. It is found that the system exhibits a variety of bifurcations including symmetry breaking, period doubling, and crises when monitoring the frequency detuning parameter in tiny steps. The multistability property of the system for special sets of its parameters is also analyzed. An experimental study of the coupled system is carried out in this work. An appropriate electronic simulator is proposed for the investigations of the dynamic behavior of the system. Correspondences are established between the coefficients of the system model and the components of the electronic circuit. A comparison of experimental and numerical results yields a very good agreement.     

Dynamic Buckling of Autonomous Systems Having Potential Energy Universal Unfoldings of Cuspoid Catastrophe

This work deals with dynamic buckling universal solutions of discrete nondissipative systems under step loading of infinite duration. Attention is focused on total potential energy functions associated with universal unfoldings of cuspoid type catastrophes with one active coordinate. The fold, dual cusp and tilted cusp catastrophes under statically applied loading occurring via limit points, asymmetric/symmetric bifurcations and nondegenerate hysteresis points are extended to the case of dynamic loading. Catastrophe manifolds of these types showing imperfection sensitivity under both types of loading are fully assessed. Important findings regarding dynamic buckling of imperfect systems generated by perfect systems associated with imperfect bifurcations are explored. The analysis is supplemented by a numerical application of a system exhibiting imperfect bifurcation when it is perfect as well as a hysteresis point associated with a tilted cusp catastrophe, when it becomes imperfect.     

Exciting dynamical behavior in a network of two coupled rings of Chen oscillators

We study exotic patterns appearing in a network of coupled Chen oscillators. Namely, we consider a network of two rings coupled through a “buffer” cell, with \(\mathbf {Z}_3\times \mathbf {Z}_5\) symmetry group. Numerical simulations of the network reveal steady states, rotating waves in one ring and quasiperiodic behavior in the other, and chaotic states in the two rings, to name a few. The different patterns seem to arise through a sequence of Hopf bifurcations, period-doubling, and halving-period bifurcations. The network architecture seems to explain certain observed features, such as equilibria and the rotating waves, whereas the properties of the chaotic oscillator may explain others, such as the quasiperiodic and chaotic states. We use XPPAUT and MATLAB to compute numerically the relevant states.     

The Smooth and Nonsmooth Travelling Wave Solutions in a Nonlinear Wave Equation

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｅｘｉｓｔｅｎｃｅｏｆｐｅａｋｏｎＴＷＳｏｆａｎｏｎｌｉｎｅａｒｗａｖｅｅｑｕａｔｉｏｎρｔ =ｂｕｘ 12 [(ｕ2 ±ｕ2 ｘ) ρ] ｘ,　ρ =ｕ±ｕｘｘ ( 1 )ｗａｓｃｏｎｓｉｄｅｒｅｄｂｒｅｉｆｌｙｂｙＰ .Ｒｏｓｅｎａｕ (ｓｅｅ [1 ] ) .Ｅｑ.( 1 )ｉｓｆｏｕｎｄｂｙ“ｒｅｓｈｕｆｆｌｉｎｇ”Ｈａｍｉｌｔｏｎｉａｎｏｐｅｒａｔｏｒｏｆｂｉ＿ＨａｍｉｌｔｉｏｎｉａｎｓｔｒｕｃｔｕｒｅｉｎＫｄＶａｎｄｍＫｄＶｅｑｕａｔｉｏｎ (ｓｅｅ [2 ] ) .ＢｅｃａｕｓｅＥｑ.( 1 )ｈａｓｓｔｒｏｎ…     

Several new types of bounded wave solutions for the generalized two-component Camassa–Holm equation

Li and Qiao studied the bifurcations and exact traveling wave solutions for the generalized two-component Camassa–Holm equation $$\begin{aligned} \left\{ \begin{array}{l} m_{t}+\sigma um_{x}-Au_{x}+2m \sigma u_{x}+3(1-\sigma )uu_{x}\\ \quad +\rho \rho _{x}=0, \\ \rho _{t} +(\rho u)_{x}=0, \end{array} \right. \end{aligned}$$ \(m=u-u_{xx}, A>0\) . They showed that there exist solitary wave solutions, cusp wave solutions, and periodic wave solutions for the equation, and their analysis focused on the bifurcations when \(\sigma >0\) . In this paper, we first complement the bifurcations when \(\sigma <0\) by following the same procedure as that of Li, and then show the existence and implicit expressions of several new types of bounded wave solutions, including solitary waves, periodic waves, compacton-like waves, and kink-like waves. In addition, the numerical simulations of the bounded wave solutions are given to show the correctness of our results.     

On the Influence of Viscosity on Riemann Solutions

We show how the existence and uniqueness of Riemann solutions are affected by the precise form of viscosity which is used to select shock waves admitting a viscous profile. We study a complete list of codimension-1 bifurcations that depend on viscosity and distinguish between Lax shock waves with and without a profile. These bifurcations are the saddle–saddle heteroclinic bifurcation, the homoclinic bifurcation, and the nonhyperbolic periodic orbit bifurcation. We prove that these influence the existence and uniqueness of Riemann solutions and affect the number and type of waves comprising a Riemann solution. We present generic situations in which viscous Riemann solutions differ from Lax solutions.     

Bifurcations of Poiseuille flow between parallel plates: Three-dimensional solutions with large spanwise wavelength

The spatial dynamics approach is applied to the analysis of bifurcations of the three-dimensional Poiseuille flow between parallel plates. In contrast to the classical studies, we impose time periodicity as well as spatial periodicity with period 2/ in the streamwise direction. However, we make no assumptions on the behavior in the spanwise direction, except the uniform closeness of the bifurcating solution to the basic flow. In an abstract setting it is shown how the dimension of the critical eigenspace of the spatial dynamics analysis can be uniquely determined from the classical linear stability problem. For the three-dimensional Poiseuille problem we are able to find all relevant coefficients from the analysis of the purely two-dimensional problem. Moreover, we are able to analyze precisely the influence of a spanwise pressure gradient and the associated spanwise mass flux. The study of the reduced problem shows that there are two different kinds of solutions (spirals and ribbons) which are 2p/ periodic in the spanwise direction, as in the Couette-Taylor problem, and both of them bifurcate in the same direction.     

Chaotic thermal convection of couple-stress fluid layer

In this work we study the pattern of bifurcations and intermittent-chaos of non-Newtonian couple-stress shallow fluid layer subject to heating from below. The couple-stress parameter delays onset of convection, synchronizes chaotic behavior, and decreases the heat transfer . Some global aspects of the dynamics such as homoclinic bifurcations and transition to chaos are explored. The effects of particle size on the intermittent-chaos regime at particular normalized Rayleigh number, say \(r=166.1\), are investigated. With the increase in couple-stress parameter, the present Lorenz-like system synchronizes to a steady state via a series of periodic solutions interspersed with intervals of chaotic behaviors.     

Homoclinic Trajectories in Discontinuous Systems

We study bifurcations of bounded solutions from homoclinic orbits for time-perturbed discontinuous systems. Functional analytic method is used. An illustrative example of a periodically perturbed piecewise linear differential equation in is presented.      

The non-uniqueness of the atomistic stress tensor and its relationship to the generalized Beltrami representation

The non-uniqueness of the atomistic stress tensor is a well-known issue when defining continuum fields for atomistic systems. In this paper, we study the non-uniqueness of the atomistic stress tensor stemming from the non-uniqueness of the potential energy representation. In particular, we show using rigidity theory that the distribution associated with the potential part of the atomistic stress tensor can be decomposed into an irrotational part that is independent of the potential energy representation, and a traction-free solenoidal part. Therefore, we have identified for the atomistic stress tensor a discrete analog of the continuum generalized Beltrami representation (a version of the vector Helmholtz decomposition for symmetric tensors). We demonstrate the validity of these analogies using a numerical test. A program for performing the decomposition of the atomistic stress tensor called MDStressLab is available online at http://mdstresslab.org.     

The measurement of the material parameters of viscoelastic fluids using a rotating sphere and a rheogoniometer

Summary In this work, measurement of the flow field around a rotating sphere has been used to obtain the material parameters of a second-order Rivlin-Ericksen fluid. Experiments were carried out with a Laser-Doppler anemometer to obtain the velocity distribution and usingGiesekus' analysis, the material parameters for the second-order fluid were obtained.

---

Zusammenfassung In dieser Untersuchung wird die Ausmessung des Strömungsfeldes um eine rotierende Kugel dazu verwendet, um die Stoffparameter einer Rivlin-Ericksen-Flüssigkeit zweiter Ordnung zu erhalten. Die Experimente zur Bestimmung der Geschwindigkeitsverteilung werden mit einem Laser-Doppler-Anemometer durchgeführt, und zur Auswertung der Parameter der Flüssigkeit zweiter Ordnung wird eine Analyse vonGiesekus benutzt.

---

Notations A ₁,A₂ Rivlin-Ericksen tensor - A ² Parameter used in eq. [12] - a Radius of the sphere - B Parameter used in eq. [12] - I Unit tensor - m ₀(₁ – ₂)/a², parameter used by ref. (8) - N ₁,N₂ First and second normal stress difference - p Isotropic pressure - Radial distance from the centre of the rotating body - S ₁,S₂ Stress tensor - v _r,v,v Velocity components in a spherical coordinate system - ₀,₁,₂ Material parameters used in eq. [2] - Shear rate - _a Apparent voscosity - ₀ Zero-shear viscosity - Angle measured from the axis of rotation - Fluid density - Stream function - Shear stress - Angular velocity With 3 figures     

Some general rheological properties of dilute polymer solutions predicted from intrinsic viscosity measurements

Expressions for the rheological properties of dilute polymer solutions at low and moderate deformation rates are established through the computation and analysis of exact Zimm's eigenvalues. It is shown that they can be expressed in terms of measurable parameters from intrinsic viscosity data. Under moderate deformation rates one needs to introduce a slippage between the solvent and the smoothed polymer to be able to describe shear-thinning behaviour.     

Bifurcation and coalescence of a plethora of homoclinic orbits for a Hamiltonian system

This is a further study of the set of homoclinic solutions (i.e., nonzero solutions asymptotic to 0 as ¦x¦) of the reversible Hamiltonian systemu ^iv +Pu +u–u ²=0. The present contribution is in three parts. First, rigorously for P –2, it is proved that there is a unique (up to translation) homoclinic solution of the above system, that solution is even, and on the zero-energy surface its orbit coincides with the transverse intersection of the global stable and unstable manifolds. WhenP=–2 the origin is a node on its local stable and unstable manifolds, and whenP(–2,2) it is a focus. Therefore we can infer, rigorously, from the discovery by Devaney of a Smale horseshoe in the dynamics on the zero energy set, there are infinitely many distinct infinite families of homoclinic solutions forP(–2, –2+) for some>0. Buffoni has shown globally that there are infinitely many homoclinic solutions for allP(–2,0], based on a different approach due to Champneys and Toland. Second, numerically, the development of the set of symmetric homoclinic solutions is monitored asP increases fromP=–2. It is observed that two branches extend fromP=–2 toP=+2 where their amplitudes are found to converge to 0 asP 2. All other symmetric solution branches are in the form of closed loops with a turning point betweenP=–2 andP=+2. Numerically it is observed that each such turning point is accompanied by, though not coincident with, the bifurcation of a branch of nonsymmetrical homoclinic orbits, which can, in turn, be followed back toP=–2. Finally, heuristic explanations of the numerically observed phenomena are offered in the language of geometric dynamical systems theory. One idea involves a natural ordering of homoclinic orbits on the stable and unstable manifolds, given by the Horseshoe dynamics, and goes some way to accounting for the observed order (in terms ofP-values) of the occurrence of turning points. The near-coincidence of turning and asymmetric bifurcation points is explained in terms of the nontransversality of the intersection of the stable and unstable manifolds in the zero energy set on the one hand, and the nontransversality of the intersection of the same manifolds with the symmetric section in ⁴ on the other. Some conjectures based on present understanding are recorded.     

Calculating the parameters of an electron beam that drifts across a strong homogeneous magnetic field

Electron drift in specified fields has been examined in [1] and, as applied to a magnetron, in [2–4] with the averaging method. In [1,2], a first- and in [3,4] in a second-order approximation of the small parameter ) E/²L was used. Here and below, E and H=(c/) are the field strengths, L is the characteristic dimension of the field heterogeneity, is the charge-mass ratio of an electron (>0), and c is the velocity of light. An attempt to construct similar approximations for a drifting electron beam with allowance for the space-charge field, within the framework of the averaging method, involves considerable mathematical difficulties. This paper describes an attempt to solve the latter problem for a stationary monoenergetic beam that drifts under the influence of a plane electric field with potential (x,y) across a strong homogeneous magnetic field H_z H=const. Solutions are constructed by the method of successive approximations, in powers of the parameter =h/L, where h is the Larmor electron radius for narrow beams with a width on the order of 2h.I thank A. N. Ievlevu for assistance in the computational and graphical work, V. Ya. Kislov for a discussion of the results, and L. A. Vainshtein for suggesting the problem examined in §3 and for critical comments.     

A Degenerate Bifurcation Structure in the Dynamics of Coupled Oscillators with Essential Stiffness Nonlinearities

We study the degenerate bifurcations of the nonlinear normal modes(NNMs) of an unforced system consisting of a linear oscillator weaklycoupled to a nonlinear one that possesses essential stiffnessnonlinearity. By defining the small coupling parameter , we study thedynamics of this system at the limit 0. The degeneracy in the dynamics ismanifested by a 'bifurcation from infinity' where a bifurcation point isgenerated at high energies, as perturbation of a state of infiniteenergy. Another (nondegenerate) bifurcation point is generated close tothe point of exact 1:1 internal resonance between the linear andnonlinear oscillators. The degenerate bifurcation structure can bedirectly attributed to the high degeneracy of the uncoupled system inthe limit 0, whose linearized structure possesses a double zero, and aconjugate pair of purely imaginary eigenvalues. First we construct localanalytical approximations to the NNMs in the neighborhoods of thebifurcation points and at other energy ranges of the system. Then, we`connect' the local approximations by global approximants, and identifyglobal branches of NNMs where unstable and stable mode and inverse modelocalization between the linear and nonlinear oscillators take place fordecreasing energy.