- About this Journal ·
- Abstracting and Indexing ·
- Aims and Scope ·
- Annual Issues ·
- Article Processing Charges ·
- Author Guidelines ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Table of Contents
Advances in Mechanical Engineering
Volume 2012 (2012), Article ID 746537, 7 pages
Numerical Analysis of Oscillation Death in Coupled Self-Excited Elastic Beams
1Departamento de Materiales, Universidad Autonoma Metropolitana-Azcapotzalco, Avenida San Pablo 180, Col. Reynosa-Tamaulipas, 02200 México, DF, Mexico
2Centro de Investigacion e Innovacion Tecnologica-IPN, Cerrada Cecati s/n, Col. Santa Catarina, 02250 México, DF, Mexico
Received 3 February 2012; Accepted 30 June 2012
Academic Editor: A. Seshadri Sekhar
Copyright © 2012 Miguel A. Barron et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
The emergence of the oscillation death phenomenon in a ring of four coupled self-excited elastic beams is numerically explored in this work. The beams are mathematically represented through partial differential equations which are solved by means of the finite differences method. A coupling scheme based on shared boundary conditions at the roots of the beams is assumed, and as initial conditions, zero velocity of the first beam and three normal vibration modes of a linear elastic beam are employed. The influence of the self-exciting constant on the ring dynamics is analyzed. It is observed that oscillation death arises as result of the singularity of the coupling matrix.
In the past years the collective behavior of coupled nonlinear oscillators has been widely studied in many disciplines, for example, physics , biology , ecology , chemistry , and mechanics . A wide diversity of nonlinear dynamic phenomena such as locking , partial synchronization , full synchronization , antiphase synchronization , and clustering  have been reported in coupled oscillators. Many coupling schemes have also been tested: local , nearest , global , diffusive , adaptive , delayed , hierarchical , and so on. An interesting behavior of coupled oscillators is amplitude death and oscillation death, which are steady states where the coupled oscillators stop their oscillation in a permanent way and become frozen in time [17–19]. Sometimes this cessation of oscillations in time is named quenching . Amplitude death arises through a Hopf bifurcation mechanism in coupled oscillators with an important parameter mismatch or in identical oscillators with time delays . An already existing unstable steady state with zero amplitude is transformed by the coupling into a stable one allowing its observation; that is, the coupling induces stability at the origin of the phase space. On the other hand, oscillation death occurs through a saddle-node bifurcation mechanism allowing the emergence of new fixed points: a new stable steady state with nonzero amplitude is created by the coupling [19, 21]. Frequently, in the literature amplitude death is confused with oscillation death [22–27]. Even the famous finding of Lord Rayleigh  related to the quenching of two organ pipes standing side by side is indistinctly considered as amplitude death or oscillation death . To date, in spite of the significant conceptual and technical differences between amplitude death and oscillation death, there is not yet a clear distinction between these phenomena. Fortunately, a precise distinction between amplitude death and oscillation death concepts can be found in .
The elimination of oscillations in a population of synthetic genetic clocks with phase-repulsive coupling is studied in . These authors report that as the coupling strength is increased, silencing of oscillations is found to occur via the appearance of an inhomogeneous limit cycle, followed by oscillation death. In  the coupled behavior of two different oscillators, namely, a Rössler oscillator and a linear one, is analyzed. The region of oscillation quenching is analytically obtained using a proper coupling strength and parameter region, and, in general, it is concluded that even though two different oscillators are coupled with each other, they exhibit oscillation quenching, just as coupled oscillators with a large parameter mismatch. The coupled behavior of a set of five oscillators that mimics phenomena in cardiovascular systems, each one with its own characteristic frequency and amplitude, is studied in . The authors derive analytic conditions that allow the prediction of oscillation death through the two aforementioned bifurcation routes and conclude that oscillation death occurs not only by changes in couplings but also by changes in the oscillator frequencies or amplitudes. In  the induction of Hopf bifurcation and oscillation death by time delays in two coupled networks, each one with four nodes, is investigated. Two-way coupling between a single node in one network with one in another is assumed. It is concluded that a system with two networks that involve delayed shortcut connections exhibits oscillation death that is caused by delay coupling. Oscillation death in a pair of coupled oscillators with diffusive and direct couplings is reported in . For larger systems, these authors report the emergence of Turing-like pattern structures and find that death states in coupled oscillators in the absence of any effects of parameter mismatch and delayed time are possible.
By far, the bulk of the published reports on the collective behavior of coupled oscillators is related with oscillators described by ordinary differential equations or maps, for example, [17, 18, 33]. In the last decade the published works on the coupled behavior of spatially extended systems described through partial differential equations are growing; however, the majority of them are focused on synchronization [34–37]. Currently, there are few reports on oscillation death in oscillators described by partial differential equations. In a recent paper by one of the present authors  the synchronization phenomenon in a ring of four self-excited elastic beams governed by partial differential equations with nearest neighbor coupling is analyzed. The self-excited elastic beams are intended to represent turbine blades coupled through the shaft and subject to vortex induced vibration. In the above paper it is shown that if the coupling constant is near to zero, the motions of the beam tips are very complex and exhibit multiple frequencies and amplitudes. As the coupling constant is increased to values just below for which the coupling matrix is singular, the motion becomes more regular, and finally the coupled beams exhibit a single frequency and amplitude. In the present paper the emergence of oscillation death in a ring of four coupled self-excited elastic beams is explored. As in , the self-excited oscillations of the beams are assumed to occur due to vortex-induced vibration. Numerical simulations show that this phenomenon occurs for the case in which the coupling matrix becomes singular. The paper is organized as follows: Section 2 presents the mathematical model and the coupling scheme. In Section 3 the mathematical conditions that allow the emergence of oscillation death are derived and explained. Section 4 describes the numerical method and the parameters employed in the solution of the four coupled partial differential equations. Section 5 presents the numerical results for the coupled behavior of the ring. The influence of the initial conditions on the emergence of oscillation death is analyzed. Besides, some numerical results on the nontrivial fixed point stability are discussed. Finally, the concluding remarks are shown in Section 6.
2. Mathematical Model and Coupling Scheme
The motion of an elastic beam of homogeneous section and properties is governed by the dimensionless equation  where , , and are dimensionless variables which correspond to the distance from the root, time, and the transverse displacement, respectively. The van der Pol oscillator  has been employed in the literature since long time ago to represent the time-varying forces on a cylinder due to vortex shedding [41–43]. Then, to simulate the self-excited oscillations in the beam due to vortex-induced vibration, the self-exciting term of the van der Pol oscillator can be added to (1), and the following equation arises : where is the self-exciting constant. As is justified in , (2) contains two key issues of a fluid-structure problem: the elastic response of the structure and the nonlinear fluid-structure interaction. These issues are responsible for the complex vibration behavior of rotating blades in turbomachinery.
A clamped-free beam is considered, which has one end fixed to the shaft whilst the other end is freely moving. This kind of beam is subject to the following boundary conditions : where and correspond to the root and the free end of the beam, respectively.
It is assumed that there are four coupled beams forming a ring, as is seen in Figure 1. Four beams are chosen for study because this is the smallest number of beams for which there is at least one beam that is not an immediate neighbor. Each beam has its own transverse displacement , . Then, the ring is mathematically represented in this way:
As in , coupling through the roots of the beams is assumed. This corresponds to transmission of vibrations to the nearest neighbors through the shaft. Transverse displacements are zero at the bases of the beams, and the coupling occurs through the slopes of the beams at that point: where at and is a coupling constant whose value determines the strength of the coupling.
From (5) the fourth-order linear system is constructed, where and the coupling matrix is given by  Clearly, is a fourth order circulant matrix whose determinant and singular values decomposition vector are, respectively,
3. Emergence of Oscillation Death
Oscillation death arises through a saddle-node bifurcation mechanism which allows the emergence of new fixed points; that is, new stable steady states with nonzero amplitude are created by the coupling [19, 21]. The trivial solution of the Van der Pol oscillator is an unstable equilibrium point located at the origin. Due to the self-excited nature of the Van der Pol oscillator, nontrivial solutions are time dependent and oscillatory; that is, solutions that start close to the equilibrium point move away from it and are transformed into closed trajectories that become a limit cycle. One would expect a similar behavior for the ring of self-excited beams represented by (4). On the other hand, fixed points require the existence of time-independent solutions of (4).
From (7), becomes singular if . This is verified for whenever . If the coupling matrix is singular, the boundary conditions of (5) employed in the coupling are linearly dependent, and the ring of coupled beams exhibits an infinite number of solutions. This implies that stable nonzero equilibrium points, that is, oscillation death, can emerge in the ring of beams. On the contrary, if the coupling matrix is nonsingular, the solutions of the coupled system represented by (4) are oscillatory in nature, as is shown in . Then, emergence of oscillation death is only possible for . Both the value of the steady-state deflection and its stability depend on the initial condition and on the value of the self-exciting constant .
From (8), for , and therefore . Then, when is singular, this matrix has one linearly dependent row or column which can be bypassed, and the ring dynamics is analyzed in a reduced space of three dimensions.
4. Numerical Solution
Many mathematical tools have been employed in the past to solve vibration problems in beams, for example, separation of variables , spectral method , modal perturbation , and so on. Given that analytical solutions of the four coupled nonlinear partial differential equations here considered are difficult to obtain, in this work a numerical solution is selected. The explicit discretization schemes suggested in  are applied. Convergence tests on time and space were carried out using dimensionless time steps and dimensionless spatial steps . Given that similar results were obtained with the previously mentioned discretization parameters, in the computer runs the following parameters were employed: and .
Recently, blade tip-timing techniques have been employed to measure the transient deflection of blades in order to analyze turbomachinery vibration . A tip-timing system allows the derivation of characteristic vibration parameters, such as blade displacement, velocity, and acceleration from the measured data in order to describe the vibration properties of the blade assembly. Therefore, in this work the tip deflection of every beam, that is, the free end transverse displacement, is monitored and tracked for analysis of the ring dynamics. The transient and steady-state behaviors of the coupled beams are analyzed from the corresponding time series of the tip deflections.
5. Numerical Results
is nonsingular for , and in this situation the collective behavior of the coupled beams evolves from a complex motion with multiple frequencies and amplitudes for around zero to full synchronization as is increased to a value near to but less than 0.5 . In the numerical simulations of the present work, a value of is kept constant in order to explore the ring dynamics when is singular. The influence of the initial conditions and the value of the self-exciting constant on the ring dynamics are studied too.
The four beams are assumed to be identical; that is, they have the same value of the self-exciting constant. The tested initial conditions were as follows: (i) the first three normal vibration modes  of a linear elastic beam ruled by (1) and zero velocity for all beams and (ii) nonzero velocity for the first beam and zero positions and velocities for the rest of the beams. For initializing the ring with modes the procedure was as follows: the first beam is initialized in certain mode, and the rest are initialized in a different one.
5.1. Short-Term Behavior
Figure 2 shows a short-term time series of the first beam tip deflections employing and with mode 1-mode 2 as initial conditions; this means that the first beam is initially in the first normal vibration mode and the rest of the beams are in the second one. In this figure a transient death of oscillations is observed for the first beam; however, as time elapses the first beam and the ring finally become oscillatory.
Figure 3 depicts the short-term dynamic behavior of the coupled beams for using mode 1-mode 3 as initial conditions. In this case there is not transient cessation of oscillations, and the dynamics response of the ring is always oscillatory. Numerical results for several combinations of modes as initial conditions show similar behavior of that shown in Figures 2 and 3; that is, permanent oscillatory response is exhibited by the ring.
The effect of the initial velocity of first beam on the short-term tip deflection is appreciated in Figure 4 using and initial velocity of the first beam equal to . The rest of the beams are initially at rest. As in Figure 2, the ring exhibits a transient death of oscillations but the final dynamic behavior is oscillatory. For and an initial velocity of the first beam equal to 2, Figure 5 shows that oscillation death is permanent. Computer run of Figure 5 was carried out for a dimensionless time long enough to detect an eventual reemergence of oscillations, as is explained below. A closer look to the time series of Figure 5 is shown in Figure 6. Steady state is reached when the tip deflection of the first beam converges to the tip deflection of the rest of the beams; this occurs in around 3 dimensionless time units. When permanent oscillation death emerges, numerical simulations show that the time-independent solutions are identical, that is, = = = = . This means that all beams in the ring converge to the same nonzero equilibrium point when oscillation death is present.
5.2. Long-Term Behavior and Oscillation Death
The length of the time to reach the steady state depends on the value of the beams self-oscillation constants and the chosen initial conditions. The long-term stability of oscillation death depends on some kind of eigenvalues of the coupled ring. These eigenvalues would indicate if the final steady state is a fixed point or a permanent oscillating one. Determining the eigenvalues of the ring lies beyond the scope of this numerical study and should constitute the topic of future work.
Emergence of permanent oscillation death when is singular is numerically detected by studying the long-term dynamic performance of the coupled beams until a steady state is reached. Computer runs were carried out for a dimensionless time of . For this time no significant change in the tip deflection of every beam is appreciated. The results of the long-term computer runs are shown in Figure 7 using several values of the self-exciting constant . From this figure one can observe that the permanent tip deflection of the beams increases as the initial velocity of the first beam is increased; however this effect is flattened as the value of is increased. Besides, numerical simulations show that oscillation death in the ring arises for values of the first beam initial velocity above a critical value, as is depicted in Figure 7. If the initial velocity of the first beam is below this critical value, a permanent oscillatory behavior of the ring is exhibited.
5.3. Stability of the Nontrivial Steady State
In order to study the stability of the oscillation death, a perturbation was introduced in the ring once the nonoscillating state is reached. Figure 8 shows the numerical results for a perturbation introduced at using and an initial velocity of the first beam equal to 2. The steady-state deflection is suddenly changed at from the aforementioned value to a lower value of 1.9711 during one time step. In this situation the ring evolves from an oscillation death state to a permanent oscillatory one. When the above steady state is perturbed from at to an upper value of 2.5343, oscillation death remains; however the new steady state becomes . Figure 9 depicts a sketch of the ball of Lyapunov stability of the ring. For a given fixed point located at , oscillation death is maintained within a ball of radius . If , a permanent oscillatory state is exhibited by the ring. Numerical simulations show that the value of depends both on and on the value of the self-exciting constant .
Short-term and long-term dynamic performance of a ring of four coupled self-excited elastic beams has been studied in this work considering singularity of the coupling matrix. Long-term results show that the coupled beams exhibit oscillation death when the coupling matrix is singular. Oscillation death is not observed when the modes of a linear beam are employed as initial conditions. However, oscillation death arises when a nonzero velocity of the first beam is assumed as initial condition. The amplitude of the permanent tip deflection depends on the values of the initial velocity and the self-exciting constant. Oscillation death emerges for values of the initial velocity above a critical value. Below this critical value, oscillation death disappears, and permanent oscillatory behavior is exhibited by the ring.
- H. Erzgräber and S. Wieczorek, “Locking behavior of three coupled laser oscillators,” Physical Review E, vol. 80, no. 2, Article ID 026212, 6 pages, 2009.
- R. E. Mirollo and S. H. Strogatz, “Synchronization of pulse-coupled biological oscillators,” SIAM Journal on Applied Mathematics, vol. 50, no. 6, pp. 1645–1662, 1990.
- B. Blasius, A. Huppert, and L. Stone, “Complex dynamics and phase synchronization in spatially extended ecological systems,” Nature, vol. 399, no. 6734, pp. 354–359, 1999.
- K. Fujii, D. Hayashi, O. Inomoto, and S. Kai, “Noise-induced entrainment between two coupled chemical oscillators in Belouzov-Zhabotinsky reactions,” Forma, vol. 15, no. 2, pp. 219–225, 2000.
- P. Perlikowski, A. Stefański, and T. Kapitaniak, “1:1 Mode locking and generalized synchronization in mechanical oscillators,” Journal of Sound and Vibration, vol. 318, no. 1-2, pp. 329–340, 2008.
- I. A. Heisler, T. Braun, Y. Zhang, G. Hu, and H. A. Cerdeira, “Experimental investigation of partial synchronization in coupled chaotic oscillators,” Chaos, vol. 13, no. 1, pp. 185–194, 2003.
- R. Carareto, F. M. Orsatti, and J. R. C. Piqueira, “Optimized network structure for full-synchronization,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 6, pp. 2536–2541, 2009.
- M. A. Barrón and M. Sen, “Synchronization of four coupled van der Pol oscillators,” Nonlinear Dynamics, vol. 56, no. 4, pp. 357–367, 2009.
- K. Czolczynski, P. Perlikowski, A. Stefanski, and T. Kapitaniak, “Clustering and synchronization of n Huygens' clocks,” Physica A, vol. 388, no. 24, pp. 5013–5023, 2009.
- H. Hong, H. Park, and M. Y. Choi, “Collective phase synchronization in locally coupled limit-cycle oscillators,” Physical Review E, vol. 70, no. 4, Article ID 045204, 4 pages, 2004.
- H. F. El-Nashar, P. Muruganandam, F. F. Ferreira, and H. A. Cerdeira, “Transition to complete synchronization in phase-coupled oscillators with nearest neighbor coupling,” Chaos, vol. 19, no. 1, Article ID 013103, 5 pages, 2009.
- M. Ataka and T. Ohta, “Anomalous synchronization of integrate-and-fire oscillators with global coupling,” Progress of Theoretical Physics, supplement 161, pp. 156–160, 2006.
- E. H. Park, Z. Feng, and D. M. Durand, “Diffusive coupling and network periodicity: a computational study,” Biophysical Journal, vol. 95, no. 3, pp. 1126–1137, 2008.
- F. Sorrentino, “Adaptive coupling for achieving stable synchronization of chaos,” Physical Review E, vol. 80, no. 5, Article ID 056206, 6 pages, 2009.
- W. He and J. Cao, “Exponential synchronization of hybrid coupled networks with delayed coupling,” IEEE Transactions on Neural Networks, vol. 21, no. 4, pp. 571–583, 2010.
- E. D. Lumer and B. A. Huberman, “Hierarchical dynamics in large assemblies of interacting oscillators,” Physics Letters A, vol. 160, no. 3, pp. 227–232, 1991.
- K. Bar-Eli, “On the stability of coupled chemical oscillators,” Physica D, vol. 14, no. 2, pp. 242–252, 1985.
- D. G. Aronson, G. B. Ermentrout, and N. Kopell, “Amplitude response of coupled oscillators,” Physica D, vol. 41, no. 3, pp. 403–449, 1990.
- W. Zou, X. G. Wang, Q. Zhao, and M. Zhan, “Oscillation death in coupled oscillators,” Frontiers of Physics in China, vol. 4, no. 1, pp. 97–110, 2009.
- Y. P. Singh, “Quenching in a system of coupled van der Pol oscillators,” Journal of Sound and Vibration, vol. 73, no. 1, pp. 73–78, 1980.
- J. J. Suarez-Vargas, J. A. González, A. Stefanovska, and P. V. E. McClintock, “Diverse routes to oscillation death in a coupled-oscillator system,” Europhysics Letters, vol. 85, no. 3, Article ID 38008, 2009.
- L. L. Rubchinsky, M. M. Sushchik, and G. V. Osipov, “Patterns in networks of oscillators formed via synchronization and oscillator death,” Mathematics and Computers in Simulation, vol. 58, no. 4–6, pp. 443–467, 2002.
- A. Pikovsky, M. Rosenblum, and J. Kurths, Synchronization a Universal Concept in Nonlinear Sciences, Cambridge University Press, Cambridge, Mass, USA, 2003.
- Z. Hou and H. Xin, “Oscillator death on small-world networks,” Physical Review E, vol. 68, no. 5, Article ID 055103, 4 pages, 2003.
- L. Basnarkov and V. Urumov, “Critical exponents of the transition from incoherence to partial oscillation death in the Winfree model,” Journal of Statistical Mechanics, vol. 2009, no. 10, Article ID P10014, 2009.
- A. P. Kuznetsov and J.P. Roman, “Synchronization of coupled anisochronous auto-oscillating systems,” Nonlinear Phenomena in Complex Systems, vol. 12, no. 1, pp. 54–60, 2009.
- A. Koseska, E. Volkov, A. Zaikin, and J. Kurths, “Inherent multistability in arrays of autoinducer coupled genetic oscillators,” Physical Review E, vol. 75, no. 3, Article ID 031916, 8 pages, 2007.
- J. W. S. Rayleigh, The Theory of Sound, Dover Publications, New York, NY, USA, 1945.
- D. S. Lee, J. W. Ryu, Y. J. Park, W. H. Kye, M. S. Kurdoglyan, and C. M. Kim, “Stabilization of a chaotic laser and quenching,” Applied Physics Letters, vol. 86, no. 18, Article ID 181104, pp. 1–3, 2005.
- E. Ullner, A. Zaikin, E. I. Volkov, and J. García-Ojalvo, “Multistability and clustering in a population of synthetic genetic oscillators via phase-repulsive cell-to-cell communication,” Physical Review Letters, vol. 99, no. 14, Article ID 148103, 4 pages, 2007.
- J. W. Ryu, D. S. Lee, Y. J. Park, and C. M. Kim, “Oscillation quenching in coupled different oscillators,” Journal of the Korean Physical Society, vol. 55, no. 2, pp. 395–399, 2009.
- C. Y. Cheng, “Induction of Hopf bifurcation and oscillation death by delays in coupled networks,” Physics Letters A, vol. 374, no. 2, pp. 178–185, 2009.
- L. Junge and U. Parlitz, “Control and synchronization of spatially extended systems,” in Proceedings of the International Symposium on Nonlinear Theory and its Applications (NOLTA '98), Crans-Montana, 1998.
- L. Junge and U. Parlitz, “Synchronization and control of coupled Ginzburg-Landau equations using local coupling,” Physical Review E, vol. 61, no. 4, pp. 3736–3742, 2000.
- G. W. Wei, “Synchronization of single-side locally averaged adaptive coupling and its application to shock capturing,” Physical Review Letters, vol. 86, no. 16, pp. 3542–3545, 2001.
- J. Bragard, E. Montbrió, C. Mendoza, S. Boccaletti, and B. Blasius, “Defect-enhanced anomaly in frequency synchronization of asymmetrically coupled spatially extended systems,” Physical Review E, vol. 71, no. 2, Article ID 025201, 4 pages, 2005.
- S. Jie, E. M. Bollt, and T. Nishikawa, “Constructing generalized synchronization manifolds by manifold equation,” SIAM Journal on Applied Dynamical Systems, vol. 8, no. 1, pp. 202–221, 2009.
- M. A. Barron and M. Sen, “Synchronization of coupled self-excited elastic beams,” Journal of Sound and Vibration, vol. 324, no. 1-2, pp. 209–220, 2009.
- K. L. Graff, Wave Motion in Elastic Solids, Ohio State University Press, Belfast, Ireland, 1975.
- B. van der Pol and J. van der Mark, “The heartbeat considered as a relaxation oscillation, and an electrical model of the heart,” Philosophical Magazine, vol. 6, pp. 763–775, 1928.
- R. E. D. Bishop and A. Y. Hassan, “The lift and drag forces on a circular cylinder in a flowing fluid,” Proceedings of the Royal Society Series, vol. A277, pp. 32–50, 1963.
- Y. S. Lee, A. F. Vakakis, L. A. Bergman, and D. M. McFarland, “Suppression of limit cycle oscillations in the van der Pol oscillator by means of passive non-linear energy sinks,” Structural Control and Health Monitoring, vol. 13, no. 1, pp. 41–75, 2006.
- R. T. Hartlen and I. G. Currie, “Lift-oscillator model of vortex induced vibration,” Journal of the Engineering Mechanics, vol. 96, pp. 577–591, 1970.
- M. A. Barron-Meza, “Vibration analysis of a self-excited elastic beam,” Journal of Applied Research and Technology, vol. 8, no. 2, pp. 227–239, 2010.
- M. Lou, Q. Duan, and G. Chen, “Modal perturbation method for the dynamic characteristics of Timoshenko beams,” Shock and Vibration, vol. 12, no. 6, pp. 425–434, 2005.
- W. T. Thomson, Theory of Vibration with Applications, Prentice Hall, Englewood Cliffs, NJ, USA, 1981.
- S. Heath and M. Imregun, “A survey of blade tip-timing measurement techniques for turbomachinery vibration,” Journal of Engineering for Gas Turbines and Power, vol. 120, no. 4, pp. 784–790, 1998.