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

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.


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][2][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][7][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.

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 ,  0 ,  1 are the displacements,  is the damping parameter,  is the nonlinear coefficient,  and  are the stiffness coefficients of the linear oscillators, and  ≪ 1 the mass of the nonlinear attachment.
After applying the linear singular transformation V =  −1/2  0 +  1/2 ,  =  −1/2  0 −  −1/2  the system assumes the form: By just expanding the solution of V and  1 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 ( 1 , 1) and ( 2 , 1), with while the nonhomogeneous part appears in ().
Notice that for  = 0 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, V =  1  1 +  2  2 ,  1 =  1 +  2 , and assume zero initial displacements and nonzero initial velocities ż 1 (0) =  1  10 , ż 2 (0) =  2  20 .We replace the solutions of the nonhomogeneous linear oscillators in the first equation of (2) and after performing a time transformation t →  20 , where  2 =  20 + , 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  20 , 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 t →  10  and continue the analysis in the same way with the resonance being between  10 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 =  0 ( 0 ,  1 , . ..) +  1 ( 0 ,  1 , . ..) +  () ,  () =  ( 0 ,  1 , . ..) =  0 ( 0 ,  1 , . ..) +  1 ( 0 ,  1 , . ..) +  () , where  0 = t is the fast time and  1 =  t is the slow time.By keeping (1) 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  0 , we examine the limit of the dynamics as  0 → ∞ and impose the conditions  N0 / 0 = 0,  0 / 0 = 0.This will provide us with the long-term behavior of the dynamics in the limit of large values of the fast time scale.Therefore t 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  t =   then system (9) becomes   =  (, ,   ) +  () ,   =  (, ,   ) +  () (15) and in the singularity limit we get  =  = 0 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][5][6].

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 ener The energy that is stored in the nonlinear oscillator is given, as a fraction of the total energy, by

Numerical Simulations
The structure of the SIM depends on the parameters , , the initial conditions ẋ 0 (0), ẋ 1 (0) 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 ẋ 0 (0), ẋ 1 (0) 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 t and the figures depict the steady state  2 0 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 ( 2 ) 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.
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.
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.
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.
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).