Research Article | Open Access

# The Effect of Slow Invariant Manifold and Slow Flow Dynamics on the Energy Transfer and Dissipation of a Singular Damped System with an Essential Nonlinear Attachment

**Academic Editor:**Mitsuhiro Ohta

#### Abstract

We study the effect of slow flow dynamics and slow invariant manifolds on the energy transfer and dissipation of a dissipative system of two linear oscillators coupled with an essential nonlinear oscillator with a mass much smaller than the masses of the linear oscillators. We calculate the slow flow of the system, the slow invariant manifold, the total energy of the system, and the energy that is stored in the nonlinear oscillator for different sets of the parameters and show that the bifurcations of the SIM and the dynamics of the slow flow play an important role in the energy transfer from the linear to the nonlinear oscillator and the rate of dissipation of the total energy of the initial system.

#### 1. Introduction

Mechanical structures where the nonlinear attachments have small masses in comparison to the structures to which they are attached and the systems are non conservative act as nonlinear energy sinks (NESs) and absorb, through irreversible transient transfer, energy from the linear parts due to resonance captures in vicinities of resonance manifolds of the underlying conservative systems for certain ranges of parameters and initial conditions [1–3].

Systems, where the masses of the nonlinear attachments are small in comparison to the linear oscillator, are singular and their dynamics are governed by different time scales. Such systems can be treated with the use of methods such as singularity analysis or multiple scale analysis [4, 5]. As it has been shown in previous works the slow flow of the system and the slow invariant manifolds (SIMs) obtained in the singular limit play a very important role in the dynamics under consideration [3, 6–8].

Furthermore, many works [1, 3, 9] have shown that energy transfer from a linear oscillator to its coupled nonlinear oscillator as well as the rate of dissipation of the energy is directly connected to the bifurcations of the SIM and the dynamics of the slow flow of the system.

In this work we study a dissipative system of two linear oscillators coupled with an essential nonlinear oscillator. In previous works [6, 7] we studied the SIM of the system and made a classification of its structure. Specifically, the structure of the SIM may be classified in three cases. First case when the SIM has always one stable branch. Second case when the SIM has always three branches, two of them stable and the one that lies between them unstable. In the above cases the SIM has no bifurcations. The third case is when the SIM Bifurcates. We also showed that the slow flow of the system has rich dynamics. It may oscillate, following the stable branch of the SIM, make relaxation oscillations, or have chaotic orbits.

In this work we study the effect of the slow invariant manifold and the dynamics of the slow flow of the system on the energy transfer from the linear to the nonlinear oscillator and the dissipation of the total energy of the system. From our numerical results we confirm that the SIM, its bifurcations, and the dynamics of the slow flow play an essential role in the energy transfer and dissipation of the system.

The paper is organized as follows. In Section 2 we derive the slow flow and the SIM of the system. In Section 3 we calculate the energy of the system and the instantaneous energy that is stored in the nonlinear oscillator. In Section 4 we present numerical simulations for each case of the SIM and the behavior of the system’s energy transfer and dissipation. Finally we conclude the paper in Section 5.

#### 2. The Slow Flow and the SIM of the System

The initial system considered in this work is composed of two coupled linear oscillators and a nonlinear oscillator interacting through an essential nonlinearity with one of the linear oscillators and a small mass in comparison to the masses of the linear oscillators. The equations of motion are given by where , , are the displacements, is the damping parameter, is the nonlinear coefficient, and are the stiffness coefficients of the linear oscillators, and the mass of the nonlinear attachment.

After applying the linear singular transformation the system assumes the form:

By just expanding the solution of and with respect to and substituting them in the second and third equations of (2) we see that at zeroth order we have only the homogeneous part of the linear equations with eigenfrequencies and corresponding eigenvectors and , with while the nonhomogeneous part appears in .

Notice that for and therefore for small the eigenfrequencies are distinct and can not have a unity ratio. Their ordering is not important.

We introduce modal coordinates for the linear part of (2) through the coordinate transformations, , , and assume zero initial displacements and nonzero initial velocities , . We replace the solutions of the nonhomogeneous linear oscillators in the first equation of (2) and after performing a time transformation , where , we derive the following reduced system: where

Considering the reduced system (5) we apply the complexification-averaging technique [3, 10] by introducing the complex coefficient . By making the important additional assumption that the nonlinear oscillator is in (1 : 1) resonance with the unperturbed eigenfrequency , that is, , we derive the averaged complex dynamical system where

We mention at this point that system (7) is approximate since it takes into account only a single “fast” frequency. Clearly, since the system under consideration is strongly nonlinear, this is only an approximation since higher harmonics will exist in the dynamics. However, we conjecture that there exist regimes where the harmonic components with normalized frequency unity dominate.

We could have chosen another time transformation by considering and continue the analysis in the same way with the resonance being between and the nonlinear oscillator.

Since we are interested in the amplitude of the oscillations of the nonlinear attachment, we use the polar representation of the complex number . After separating the real and imaginary parts, we obtain the following two equations for and , representing the slow flow of the system:

We use multiple scales analysis where is the fast time and is the slow time. By keeping terms in (9) for and we derive the equations To study the steady state dynamics of the above system, in terms of the fast time scale , we examine the limit of the dynamics as and impose the conditions , . This will provide us with the long-term behavior of the dynamics in the limit of large values of the fast time scale. Therefore is not a constant variable and its variation provide us the variation of the SIM. Then, from (11) we find Manipulating expressions (12) we derive the steady state phase with the steady state amplitude given by Equations (13) and (14) represent the SIM of the dynamics of (5) [1, 3].

This is equivalent to singularity analysis of (9). Notice that if we consider then system (9) becomes and in the singularity limit we get which is exactly the SIM of the system, that is, (14).

There are many theorems for this singular manifold which allow us to determine the long-term behavior of (9) [4–6].

#### 3. The Energy of the System

In order to study the effect of the SIM and the slow flow on the energy transfer and dissipation of our initial system we will study the rate of the total energy dissipation of the system. We use as a tool the graph of the logarithm of the total energy versus time for the different cases of the SIM.

Furthermore, in order to study the possible energy transfer from the linear to the nonlinear oscillators we calculate the instantaneous energy that is stored in the nonlinear oscillator.

The total energy of the system is calculated by the Hamiltonian of the initial system (1) without the dissipative terms and it is given by

The energy that is stored in the nonlinear oscillator is given, as a fraction of the total energy, by

#### 4. Numerical Simulations

The structure of the SIM depends on the parameters , , the initial conditions of the initial system (1), and the damping parameter . The SIM may have one branch that is always stable, three branches, two of them stable and the one that lies between them unstable, or it may bifurcate.

In what follows we present numerical simulations for different sets of parameters and initial conditions.

For every set of the parameters , there are different structures of the SIM that are related to the different initial conditions and the structure of the SIM affects the dynamics of the initial system.

Since we are interested only in the amplitude of the oscillation, the SIM is computed by finding the real roots of (14) for every time and the figures depict the steady state with respect to the fast time. The slow flow is computed by solving system (9) numerically, with the use of Runge-Kutta method.

In order to compute the total energy of the system and the amount of energy that is stored in the nonlinear oscillator we solve system (1) numerically, with the use of Runge-Kutta method, calculate the total energy from (16) and the amount of the energy that is stored in the nonlinear oscillator from (17).

Our results are depicted in three types of graphs. The first one gives the SIM and the amplitude of the slow flow () of the system versus time, the second gives the amount of energy that is stored in the nonlinear oscillator versus time, and the third gives the logarithm of the total energy of the system versus time in order to follow the rate of dissipation of the total energy.

First, we study the cases where the SIM has no bifurcations and the slow flow has regular orbits, that is, when the SIM has always one branch or when the SIM has always three branches.

In Figures 1 and 2 we present some examples where the SIM has only one stable branch. In these cases the slow flow follows the SIM. The total energy of the system dissipates smoothly and, as we observe in the diagrams of the instantaneous energy that is stored in the nonlinear oscillator, there is no energy transfer from the linear to the nonlinear oscillator.

**(a) Gray line: SIM (), red line: slow flow ()**

**(b) Instantaneous energy stored in the nonlinear oscillator**

**(c) Total energy**

**(a) Gray line: SIM (), red line: slow flow ()**

**(b) Instantaneous energy stored in the nonlinear oscillator**

**(c) Total energy**

An important remark here is that, for the cases where the slow flow oscillates rapidly before it follows the SIM (Figures 3 and 4), we detect energy transfer to the nonlinear oscillator of the system in the time interval before the slow flow follows the SIM. Furthermore, in these cases, the total energy of the system dissipates faster for the time interval when the slow flow oscillates rapidly around the SIM than the period of time when the slow flow follows the SIM.

**(a) Gray line: SIM (), red line: slow flow ()**

**(b) Instantaneous energy stored in the nonlinear oscillator**

**(c) Total energy**

**(a) Gray line: SIM (), red line: slow flow ()**

**(b) Instantaneous energy stored in the nonlinear oscillator**

**(c) Total energy**

The same behavior, as in the case where the SIM has only one stable branch, is seen when the SIM has always three branches (Figure 5). The amount of energy that is transfered from the linear to the nonlinear oscillator depends on the initial conditions, that is, the initial energy that is given to the system.

**(a) Gray line: SIM (), red line: slow flow ()**

**(b) Instantaneous energy stored in the nonlinear oscillator**

**(c) Total energy**

The next case in our study is where the SIM bifurcates. As we observe in Figures 6, 7, and 8, the SIM has three real roots and then suddenly the two of them disappear (a bifurcation occurs), one real root remains for a period of time, and then two other real roots suddenly reappear (another bifurcation occurs). In this case, energy is transferred to the nonlinear oscillator, irrespectively whether the slow flow performs relaxation oscillations or follows other cases of bifurcations. The rate of the total energy dissipation depends on whether energy is transferred to the nonlinear oscillator or not. Specifically, the system disipates energy faster in the time intervals when there is energy transfer to the nonlinear oscillator.

**(a) Gray line: SIM (), red line: slow flow ()**

**(b) Instantaneous energy stored in the nonlinear oscillator**

**(c) Total energy**

**(a) Gray line: SIM (), red line: slow flow ()**

**(b) Instantaneous energy stored in the nonlinear oscillator**

**(c) Total energy**

**(a) Gray line: SIM (), red line: slow flow ()**

**(b) Instantaneous energy stored in the nonlinear oscillator**

**(c) Total energy**

We must indicate here that the behavior of the initial system, in the case where the slow flow oscillates rapidly until it follows the SIM (with no bifurcations), differs from its behavior when there are bifurcations.

In the first case energy transfer from the linear to the nonlinear oscillator occurs in the beginning of the oscillations, until the time when the system oscillates in a regular way (Figure 9). In the second case energy transfer from the linear to the nonlinear oscillator occurs at a time when the SIM bifurcates. After the energy transfer, the system dissipates its energy and the oscillations fade out (Figure 10).

**(a) , , , , and**

**(b) , , , , and**

**(a) , , , , and**

**(b) , , , , and**

#### 5. Conclusions

From the above study we conclude that the bifurcations of the SIM and the dynamics of the slow flow play an important role in the energy transfer from the linear to the nonlinear oscillator and the dissipation of the total energy of the initial system.

When the SIM has no bifurcations, there is no energy transfer from the linear to the nonlinear oscillator and the energy dissipates smoothly. When the SIM has bifurcations, then energy transfer occurs. Furthermore, when there is energy transfer to the nonlinear oscillator, the rate of the dissipation of the total energy of the system becomes larger.

When the slow flow oscillates rapidly around the SIM, the energy is transfered to the nonlinear oscillator. The amount of energy that transfers is related to the initial energy given to the system. The difference between this case and the previous one is that when the SIM has bifurcations the system dissipates its energy and the oscillations fade out.

The damping parameter determines whether the SIM bifurcates; that is, in order for the SIM to bifurcate, the damping parameter must satisfy the relation . Therefore, from the above analysis, we conclude that the damping parameter determines the ability of the system to transfer energy from the linear to the nonlinear oscillator and plays a role in the rate of the total energy dissipation of the system.

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### References

- O. V. Gendelman, Y. Starosvetsky, and M. Feldman, “Attractors of harmonically forced linear oscillator with attached nonlinear energy sink I: description of response regimes,”
*Nonlinear Dynamics*, vol. 51, no. 1-2, pp. 31–46, 2008. View at: Publisher Site | Google Scholar - O. V. Gendelman, A. F. Vakakis, L. A. Bergman, and D. M. McFarland, “Asymptotic analysis of passive nonlinear suppression of aeroelastic instabilities of a rigid wing in subsonic flow,”
*SIAM Journal on Applied Mathematics*, vol. 70, no. 5, pp. 1655–1677, 2010. View at: Publisher Site | Google Scholar | MathSciNet - A. F. Vakakis, O. V. Gendelman, L. A. Bergman, D. M. McFarland, G. Kerschen, and Y. S. Lee,
*Nonlinear Target Energy Transfer in Mechanical and Structural Systems*, Springer, 2008. - N. Fenichel, “Geometric singular perturbation theory for ordinary differential equations,”
*Journal of Differential Equations*, vol. 31, no. 1, pp. 53–98, 1979. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - F. Verhulst, “Singular perturbation methods for slow-fast dynamics,”
*Nonlinear Dynamics*, vol. 50, no. 4, pp. 747–753, 2007. View at: Publisher Site | Google Scholar | MathSciNet - J. O. Maaita, E. Meletlidou, A. F. Vakakis, and V. Rothos, “The effect of Slow Flow Dynamics on the Oscillations of a singular damped system with an essentially nonlinear attachment,”
*Journal of Applied Nonlinear Dynamics*, vol. 2, no. 4, pp. 315–328, 2013. View at: Publisher Site | Google Scholar - J. O. Maaita, E. Meletlidou, A. F. Vakakis, and V. Rothos, “The dynamics of the slow flow of a singular damped nonlinear system and its Parametric Study,”
*Journal of Applied Nonlinear Dynamics*, vol. 3, no. 1, pp. 37–49, 2014. View at: Publisher Site | Google Scholar - A. F. Vakakis, “Relaxation oscillations, subharmonic orbits and chaos in the dynamics of a linear lattice with a local essentially nonlinear attachment,”
*Nonlinear Dynamics*, vol. 61, no. 3, pp. 443–463, 2010. View at: Publisher Site | Google Scholar | MathSciNet - O. V. Gendelman, “Bifurcations of nonlinear normal modes of linear oscillator with strongly nonlinear damped attachment,”
*Nonlinear Dynamics*, vol. 37, no. 2, pp. 115–128, 2004. View at: Publisher Site | Google Scholar | MathSciNet - L. I. Manevitch, “Complex representation of dynamics of coupled oscillators,” in
*Mathematical Models of Nonlinear Excitations, Transfer Dynamics and Control in Condensed Systems*, pp. 269–300, Kluwer Academic, New York, NY, USA, 1999. View at: Google Scholar

#### Copyright

Copyright © 2014 Jamal-Odysseas Maaita and Efthymia Meletlidou. 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.