Dynamics of a model of two delay-coupled relaxation oscillators

This paper investigates the dynamics of a new model of two coupled relaxation oscillators. The model replaces the usual DDE (differential-delay equation) formulation with a discrete-time approach with jumps. Existence, bifurcation and stability of in-phase periodic motions is studied. Simple periodic motions, which involve exactly two jumps per period, are found to have large plateaus in parameter space. These plateaus are separated by regions of complicated dynamics, reminiscent of the Devil’s Staircase. Stability of motions in the in-phase manifold are contrasted with stability of motions in the full phase space.     

Winding Numbers and Average Frequencies in Phase Oscillator Networks

We study networks of coupled phase oscillators and show that network architecture can force relations between average frequencies of the oscillators. The main tool of our analysis is the coupled cell theory developed by Stewart, Golubitsky, Pivato, and Torok, which provides precise relations between network architecture and the corresponding class of ODEs in R^M and gives conditions for the flow-invariance of certain polydiagonal subspaces for all coupled systems with a given network architecture. The theory generalizes the notion of fixed-point subspaces for subgroups of network symmetries and directly extends to networks of coupled phase oscillators. For systems of coupled phase oscillators (but not generally for ODEs in R^M, where M ≥ 2), invariant polydiagonal subsets of codimension one arise naturally and strongly restrict the network dynamics. We say that two oscillators i and j coevolve if the polydiagonal θ_i = θ_j is flow-invariant, and show that the average frequencies of these oscillators must be equal. Given a network architecture, it is shown that coupled cell theory provides a direct way of testing how coevolving oscillators form collections with closely related dynamics. We give a generalization of these results to synchronous clusters of phase oscillators using quotient networks, and discuss implications for networks of spiking cells and those connected through buffers that implement coupling dynamics.     

The dynamics ofn weakly coupled identical oscillators

Summary We present a framework for analysing arbitrary networks of identical dissipative oscillators assuming weak coupling. Using the symmetry of the network, we find dynamically invariant regions in the phase space existing purely by virtue of their spatio-temporal symmetry (the temporal symmetry corresponds to phase shifts). We focus on arrays which are symmetric under all permutations of the oscillators (this arises with global coupling) and also on rings of oscillators with both directed and bidirectional coupling. For these examples, we classify all spatio-temporal symmetries, including limit cycle solutions such as in-phase oscillation and those involving phase shifts. We also show the existence of “submaximal” limit cycle solutions under generic conditions. The canonical invariant region of the phase space is defined and used to investigate the dynamics. We discuss how the limit cycles lose and gain stability, and how symmetry can give rise to structurally stable heteroclinic cycles, a phenomenon not generically found in systems without symmetry. We also investigate how certain types of coupling (including linear coupling between oscillators with symmetric waveforms) can give rise to degenerate behaviour, where the oscillators decouple into smaller groups.     

Hereditary effects of exponentially damped oscillators with past histories

This paper presents hereditary effects of exponentially damped oscillators with past histories. Unlike the classical viscously damped oscillators, the nonviscously damped ones involve damping forces which depend on time-histories of vibrating motions via convolution integrals. As a result, equations of motion of such systems are a set of coupled second-order Volterra integro-differential equations. In this work, initial value problems for the integro-differential equations are revisited. The initial conditions should contain time-histories of vibrating motions. Then, initialization response of exponentially damped oscillators is obtained. It is used to characterize the hereditary effects on the dynamic response. At last, stability of initialization response is proved from the theoretical viewpoint and verified by numerical simulations. This reveals that the hereditary effects gradually recede with increasing of time.     

Cluster synchronization,switching and spatiotemporal coding in a phase oscillator network

A network of five globally-coupled identical phase oscillators is considered. Cluster states consisting of two synchronized pairs of oscillators and one singleton are investigated. Forcing the system with non-uniform constant inputs results in regular switches between cluster states. The resultant cyclic sequences of switches (spatiotemporal codes) are studied for different initial conditions and input configurations. Implications on information coding in neural systems are briefly discussed. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Integer lattices of the action variables for the generalized Lagrange case

The lattices generated by integer-valued isolines of action variables are constructed for some integrable Hamiltonian systems with two degrees of freedom (the generalized Lagrange case). The monodromy matrices are calculated for critical points of this system.     

Switching dynamics of multiple linear oscillators

In this paper, the switching dynamics of linear oscillators with arbitrary discontinuous forcing are investigated through the concept of switching systems, and such switching systems consist of countable prescribed linear oscillators with different external excitations. The traditional treatments are to smoothen the discontinuity at switching points of two subsystems in a switching system, which can provide an approximate solution only. Therefore, an alternative method is presented to obtain an exact solution of the resultant switching linear system. Under periodic piecewise forcing and random forcing, the corresponding exact solutions and stochastic responses of switching linear systems are developed. For any periodic forcing, the periodic responses and stability of the resultant system composed of multiple linear oscillators in different time intervals are presented. In addition, the resultant switching system consisting of two oscillators are discussed, and the corresponding stability analysis is carried out.     

Geometrical properties of coupled oscillators at synchronization

We study the synchronization of N nearest neighbors coupled oscillators in a ring. We derive an analytic form for the phase difference among neighboring oscillators which shows the dependency on the periodic boundary conditions. At synchronization, we find two distinct quantities which characterize four of the oscillators, two pairs of nearest neighbors, which are at the border of the clusters before total synchronization occurs. These oscillators are responsible for the saddle node bifurcation, of which only two of them have a phase-lock of phase difference equals ± π/2. Using these properties we build a technique based on geometric properties and numerical observations to arrive to an exact analytic expression for the coupling strength at full synchronization and determine the two oscillators that have a phase-lock condition of ± π/2.     

Numerical Analysis of the Synchronization of Coupled Chemical Oscillators

A reaction–diffusion model describing a system of coupled oscillators is constructed and investigated. The oscillators in this study are chemical oscillators that represent an oscillatory heterogeneous catalytic reaction in a granular catalyst layer. The oscillators are arranged serially in the reagent stream and are coupled through the gaseous phase. The dynamic behavior of the system is investigated as a function of the main external parameter — the partial pressure of one of the reagents in the gaseous phase. Existence regions of regular and chaotic oscillations are identified. Synchronization conditions are established for the oscillations in such a chain of coupled chemical oscillators.     

Isotropy of Angular Frequencies and Weak Chimeras with Broken Symmetry

The notion of a weak chimeras provides a tractable definition for chimera states in networks of finitely many phase oscillators. Here, we generalize the definition of a weak chimera to a more general class of equivariant dynamical systems by characterizing solutions in terms of the isotropy of their angular frequency vector—for coupled phase oscillators the angular frequency vector is given by the average of the vector field along a trajectory. Symmetries of solutions automatically imply angular frequency synchronization. We show that the presence of such symmetries is not necessary by giving a result for the existence of weak chimeras without instantaneous or setwise symmetries for coupled phase oscillators. Moreover, we construct a coupling function that gives rise to chaotic weak chimeras without symmetry in weakly coupled populations of phase oscillators with generalized coupling.     

Dynamics of a ring of three coupled relaxation oscillators

The dynamics of a ring of three identical relaxation oscillators is shown to exhibit a variety of periodic motions, including clockwise and counter-clockwise wave-like modes, and a synchronous mode in which all three oscillators are in phase. The model involves individual oscillators which exhibit sudden jumps, modeling the relaxation oscillations of van der Pol oscillators. Methods include (i) numerical integration, (ii) a semi-analytical method involving solving transcendental equations numerically, and (iii) perturbation methods. A variety of bifurcations of the periodic motions are identified. This work is motivated by application to the design of a decision-making machine which can sort initial conditions according to their steady state.     

Synchronization of a chaotic electronic circuit system with cubic term via adaptive feedback control

This paper presents an adaptive feedback control scheme for the synchronization of the chaotic system consisting of Van der Pol oscillators coupled to linear oscillators with cubic term when the parameters of the master system are unknown and different with the those of the slave system. Based on the Lyapunov stability theory, an adaptive control law is derived to make the states of two slightly mismatched chaotic systems asymptotically synchronized. This method is efficient and easy to implement. Numerical simulations results confirming the analytical predictions are shown and pspice simulations are also performed to confirm the efficiency of the proposed control scheme.     

The Spectrum of the Partially Locked State for the Kuramoto Model

We solve a long-standing stability problem for the Kuramoto model of coupled oscillators. This system has attracted mathematical attention, in part because of its applications in fields ranging from neuroscience to condensed-matter physics, and also because it provides a beautiful connection between nonlinear dynamics and statistical mechanics. The model consists of a large population of phase oscillators with all-to-all sinusoidal coupling. The oscillators' intrinsic frequencies are randomly distributed across the population according to a prescribed probability density, here taken to be unimodal and symmetric about its mean. As the coupling between the oscillators is increased, the system spontaneously synchronizes: The oscillators near the center of the frequency distribution lock their phases together and run at the same frequency, while those in the tails remain unlocked and drift at different frequencies. Although this "partially locked" state has been observed in simulations for decades, its stability has never been analyzed mathematically. Part of the difficulty is in formulating a reasonable infinite-N limit of the model. Here we describe such a continuum limit, and prove that the corresponding partially locked state is, in fact, neutrally stable, contrary to what one might have expected. The possible implications of this result are discussed. An erratum to this article is available at .     

Microscopical investigation of wave propagation phenomena in residual saturated porous media

Wave propagation is used in many fields for measurement and characterization. Corresponding multiphase models usually use a continuous approach. Nevertheless, systems like wetted rocks may be saturated residually in certain situations. In such cases, one fluid is distributed as clusters, each different in size and shape. One single, continuous phase cannot account for a variety of fluid clusters, either disconnected from each other or connected only about thin liquid films. Therefore, we present a model that considers a heterogeneous distribution of disconnected fluid clusters in the form of harmonic oscillators. These oscillators are described and distinguished by their mass, damping and eigenfrequency. Hence, the model allows to characterize different clusters and includes an additional damping mechanism due to oscillations of the fluid clusters. (© 2012 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Nonstationary energy localization vs conventional stationary localization in weakly coupled nonlinear oscillators

In this paper we study the effect of nonstationary energy localization in a nonlinear conservative resonant system of two weakly coupled oscillators. This effect is alternative to the well-known stationary energy localization associated with the existence of localized normal modes and resulting from a local topological transformation of the phase portraits of the system. In this work we show that nonstationary energy localization results from a global transformation of the phase portrait. A key to solving the problem is the introduction of the concept of limiting phase trajectories (LPTs) corresponding to maximum possible energy exchange between the oscillators. We present two scenarios of nonstationary energy localization under the condition of 1:1 resonance. It is demonstrated that the conditions of nonstationary localization determine the conditions of efficient targeted energy transfer in a generating dynamical system. A possible extension to multi-particle systems is briefly discussed.     

Synchronization of stochastic coupled systems via feedback control based on discrete-time state observations

In this paper, we formulate and investigate the synchronization of stochastic coupled systems via feedback control based on discrete-time state observations (SCSFD). The discrete-time state feedback control is used in the drift parts of response system. Combining Lyapunov method with graph theory, the upper bound of duration between two consecutive state observations is provided. And a global Lyapunov function of SCSFD is presented, which derives some sufficient criteria to guarantee the synchronization of drive–response systems in the sense of mean-square asymptotical synchronization. In addition, the theoretical results are applied to stochastic coupled oscillators and second-order Kuramoto oscillators. Finally, two numerical examples are given to verify the effectiveness of the theoretical results.     

Homoclinic Tangles, Phase Locking, and Chaos in a Two-Frequency Perturbation of Duffing's Equation

Summary. We study a two-frequency perturbation of Duffing's equation. When the perturbation is small, this system has a normally hyperbolic invariant torus which may be subjected to phase locking. Applying a version of Melnikov's method for multifrequency systems, we detect the occurrence of transverse intersection between the stable and unstable manifolds of the invariant torus. We show that if the invariant torus is not subjected to phase locking, then such a transverse intersection yields chaotic dynamics. When the invariant torus is subjected to phase locking, the situation is different. In this case, there exist two periodic orbits which are created in a saddle-node bifurcation. Using another version of Melnikov's method for slowly varying oscillators, we also give conditions under which the stable and unstable manifolds of the periodic orbits intersect transversely and hence chaotic dynamics may occur. Our results reveal that when the invariant torus is subjected to phase locking, chaotic dynamics resulting from transverse intersection between its stable and unstable manifolds may be interrupted. Received November 18, 1993; final revision received September 9, 1997; accepted October 27,1997     

Analysis of emitter-coupled multivibrators by singularly perturbed systems

Emitter-coupled multivibrators play a decisive role in electrical engineering, especially for phase locked loops which are key-building blocks of analogue RF front-ends. Since multivibrators correspond to relaxation oscillators, in the following the modelling and analysis by the theory of singularly perturbed systems is presented. Models for fast and slow phenomena are derived, and the fast transients of emitter-coupled multivibrators are analysed for the first time. The results of our analysis lead to significant advantages for the design of electrical multivibrators.     

Polynomial functions on upper triangular matrix algebras

There are two kinds of polynomial functions on matrix algebras over commutative rings: those induced by polynomials with coefficients in the algebra itself and those induced by polynomials with scalar coefficients. In the case of algebras of upper triangular matrices over a commutative ring, we characterize the former in terms of the latter (which are easier to handle because of substitution homomorphism). We conclude that the set of integer-valued polynomials with matrix coefficients on an algebra of upper triangular matrices is a ring, and that the set of null-polynomials with matrix coefficients on an algebra of upper triangular matrices is an ideal.     

Passages through resonance in weakly nonlinear systems

A method of determining asymptotic expansions for weakly couplednonlinearly perturbed systems of harmonic oscillators with slowlyvarying frequencies is presented. In an example with two oscillators,each one experiences a separate resonance passage that producesa first-order amplitude change. Simultaneously, second-orderadjustments occur to both oscillators. The determination isachieved by carrying the calculations to third order.