# 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

Authors: Jamal-Odysseas Maaita; Efthymia Meletlidou
Journal: Journal of Nonlinear Dynamics
 (2014)

---

## 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.


---

## Body

## 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 Section2 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(1)ϵy¨+ϵλ(y˙-x˙0)+C(y-x0)3=0,x¨0+d(x0-x1)=ϵλ(y˙-x˙0)+C(y-x0)3,x¨1+ax1+d(x1-x0)=0,
where y, x0, x1 are the displacements, λ is the damping parameter, C is the nonlinear coefficient, a and d are the stiffness coefficients of the linear oscillators, and ϵ≪1 the mass of the nonlinear attachment.After applying the linear singular transformationv=ϵ-1/2x0+ϵ1/2y,w=ϵ-1/2x0-ϵ-1/2y the system assumes the form:
(2)w¨+λ(1+ϵ)w˙+C(1+ϵ)w3=-dv+ϵw1+ϵ+dx1,v¨+d1+ϵv-dx1=-ϵdw1+ϵ,x¨1+(a+d)x1-d1+ϵv=ϵdw1+ϵ.By just expanding the solution ofv and x1 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
(3)ω12=(-4ad(1+ϵ)+((1+ϵ)a+(2+ϵ)d)2(1+ϵ)a+(2+ϵ)dhhh--4ad(1+ϵ)+((1+ϵ)a+(2+ϵ)d)2)×(2(1+ϵ))-1,ω22=(-4ad(1+ϵ)+((1+ϵ)a+(2+ϵ)d)2(1+ϵ)a+(2+ϵ)dhhh+-4ad(1+ϵ)+((1+ϵ)a+(2+ϵ)d)2)×(2(1+ϵ))-1,
and corresponding eigenvectors (K1,1) and (K2,1), with
(4)K1=(-4ad(1+ϵ)+((1+ϵ)a+(2+ϵ)d)2(1+ϵ)a+ϵdhhh--4ad(1+ϵ)+((1+ϵ)a+(2+ϵ)d)2)×(2d)-1,K2=(-4ad(1+ϵ)+((1+ϵ)a+(2+ϵ)d)2(1+ϵ)a+ϵdhhh+-4ad(1+ϵ)+((1+ϵ)a+(2+ϵ)d)2)×(2d)-1,
while the nonhomogeneous part appears in O(ϵ).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=K1z1+K2z2, x1=z1+z2, and assume zero initial displacements and nonzero initial velocities z˙1(0)=ω1z10, z˙2(0)=ω2z20. We replace the solutions of the nonhomogeneous linear oscillators in the first equation of (2) and after performing a time transformation t^→ω20t, where ω2=ω20+ϵB-, we derive the following reduced system:
(5)w′′+λ^w′+C^w3=A^sin(ω1ω20t^)+B^sin(1+ϵB-ω20t^)+O(ϵ),
where
(6)C^=Cω202,λ^=λω20,A^=Aω202,B^=Bω202,A=(d-dK1)z10,B=(d-dK2)z20.Considering the reduced system (5) we apply the complexification-averaging technique [3, 10] by introducing the complex coefficient Ψ=w′+jw. By making the important additional assumption that the nonlinear oscillator is in (1 : 1) resonance with the unperturbed eigenfrequency ω20, that is, Ψ=ϕejt, we derive the averaged complex dynamical system
(7)ϕ′+(λ^2+j2)ϕ-3C^j8|ϕ|2ϕ+J2+jA~2+B^j2ej(ϵB-/ω20)t^+O(ϵ)=0,
where
(8)J=A^ω20ω1π(ω12-ω202)(cos⁡(ω1ω202π)-1),A~=A^ω202π(ω12-ω202)sin(ω1ω202π).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 consideringt^→ω10t 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ϕ=N(t)ejη(t). After separating the real and imaginary parts, we obtain the following two equations for N and η, representing the slow flow of the system:
(9)N′=-λ^N2-J2cos⁡(η)-A~2sin(η)+B^2sin(ϵB-ω20t^-η)+O(ϵ),Nη′=-N2+3C^N38+J2sin(η)-A~2cos⁡(η)-B^2cos⁡(ϵB-ω20t^-η)+O(ϵ).We use multiple scales analysis(10)N(t)=N(t0,t1,…)=N0(t0,t1,…)+ϵN1(t0,t1,…)+O(ϵ),η(t)=η(t0,t1,…)=η0(t0,t1,…)+ϵη1(t0,t1,…)+O(ϵ),
where t0=t^ is the fast time and t1=ϵt^ is the slow time. By keeping O(1) terms in (9) for N and η we derive the equations
(11)∂N0∂t0+λ^N02+J2cos⁡(η0)+A~2sin(η0)-B^2sin(ϵB-t^ω20-η0)=0,N0∂η0∂t0+N02-3C^N038-J2sin(η0)+A~2cos⁡(η0)+B^2cos⁡(ϵB-t^ω20-η0)=0.
To study the steady state dynamics of the above system, in terms of the fast time scale t0, we examine the limit of the dynamics as t0→∞ and impose the conditions ∂N^0/∂t0=0, ∂n0/∂t0=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
(12)λ^N^02=-J2cos⁡(η^0)-A~2sin(η^0)+B^2sin(ϵB-t^ω20-η^0),N^02-3C^N^038=J2sin(η^0)-A~2cos⁡(η^0)-B^2cos⁡(ϵB-t^ω20-η^0).
Manipulating expressions (12) we derive the steady state phase
(13)cos⁡η^0=((A~+B^cos⁡(ϵB-ω20t^))(3CN^034-N^0)hhihh+λ^N^0(B^sin(ϵB-ω20t^)-J))×(J2+A~2+2A~B^cos⁡(ϵB-ω20t^)hhihhh+B^2-2JB^sin(ϵB-ω20t^))-1,
with the steady state amplitude given by
(14)N^06-83C^N^04+16(λ^2+1)9C^2N^02=169C^2(A~2+2A^B^cos⁡(ϵB-ω20t^)+B^2hhhihhhh+J2-2JB^sin(ϵB-ω20t^)).
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^=t′ then system (9) becomes
(15)ϵN′=f(N,n,t′)+O(ϵ),ϵNn′=g(N,n,t′)+O(ϵ)
and in the singularity limit we get f=g=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–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
(16)ener(t)=12(ϵy˙2+x˙02+x˙12)+C4(y-x0)4+d2(x0-x1)2+a2x12.The energy that is stored in the nonlinear oscillator is given, as a fraction of the total energy, by(17)nlnE(t)=(ϵ2y˙2+C4(y-x0)4)×(d2(x0-x1)2+a2x1212(ϵy˙2+x˙02+x˙12)+C4(y-x0)4hhhhh+d2(x0-x1)2+a2x12)-1·100%.


## 4. Numerical Simulations

The structure of the SIM depends on the parametersa, d, the initial conditions x˙0(0),x˙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 parametersa, d there are different structures of the SIM that are related to the different initial conditions x˙0(0),x˙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 N02 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 (N2) 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 Figures1 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.One stable branch of the SIM:a=0.1, d=1.0, λ^=0.8, x0=0.09, and x1=0.05.
(a)
Gray line: SIM (N02), red line: slow flow (N2)
(b)
Instantaneous energy(%) stored in the nonlinear oscillator
(c)
Total energyOne stable branch of the SIM:a=0.2, d=0.2, λ^=0.5, x0=0.01, and x1=0.01.
(a)
Gray line: SIM (N02), red line: slow flow (N2)
(b)
Instantaneous energy(%) stored in the nonlinear oscillator
(c)
Total energyAn important remark here is that, for the cases where the slow flow oscillates rapidly before it follows the SIM (Figures3 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.One stable branch of the SIM:a=0.01, d=1.0, λ^=0.1, x0=0.9, and x1=0.5.
(a)
Gray line: SIM (N02), red line: slow flow (N2)
(b)
Instantaneous energy(%) stored in the nonlinear oscillator
(c)
Total energyOne stable branch of the SIM:a=0.2, d=0.2, λ^=0.05, x0=0.01, and x1=0.01.
(a)
Gray line: SIM (N02), red line: slow flow (N2)
(b)
Instantaneous energy(%) stored in the nonlinear oscillator
(c)
Total energyThe same behavior, as in the case where the SIM has only one stable branch, is seen when the SIM has always three branches (Figure5). 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.Always three branches of the SIM:a=0.01, d=1.5, λ^=0.1, x0=0.09, and x1=0.05.
(a)
Gray line: SIM (N02), red line: slow flow (N2)
(b)
Instantaneous energy(%) stored in the nonlinear oscillator
(c)
Total energyThe next case in our study is where the SIM bifurcates. As we observe in Figures6, 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.Relaxation oscillations of the slow flow:a=4, d=4, λ^=0.15, x0=0.9, and x1=0.5.
(a)
Gray line: SIM (N02), red line: slow flow (N2)
(b)
Instantaneous energy(%) stored in the nonlinear oscillator
(c)
Total energyThe SIM bifurcates:a=6, d=6, λ^=0.5, x0=0.9, and x1=0.04.
(a)
Gray line: SIM (N02), red line: slow flow (N2)
(b)
Instantaneous energy(%) stored in the nonlinear oscillator
(c)
Total energyThe SIM bifurcates:a=6, d=6, λ^=0.15, x0=0.9, and x1=0.5.
(a)
Gray line: SIM (N02), red line: slow flow (N2)
(b)
Instantaneous energy(%) stored in the nonlinear oscillator
(c)
Total energyWe 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 (Figure9). 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).Oscillations of the initial system when the SIM has no bifurcations. Black line:y(t), thick gray line: x0(t), and dashed gray line: x1(t).
(a)
a
=
0.01, d=1, λ^=0.1, x0=0.9, and x1=0.5
(b)
a
=
0.2, d=0.2, λ^=0.05, x0=0.01, and x1=0.01Oscillations of the initial system when the SIM bifurcates. Black line:y(t), thick gray line: x0(t), and dashed gray line: x1(t).
(a)
a
=
4, d=4, λ^=0.15, x0=0.9, and x1=0.5
(b)
a
=
6, d=6, λ^=0.15, x0=0.9, and x1=0.5


## 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 λ^<1/3. 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.



---

*Source: 208171-2014-09-01.xml*