1. Introduction
In this work, we first consider the piezoelectric coupled system given by
where $u=u(x,t)$ represents the longitudinal displacement of the center line, $v=v(x, t)$ is the total load of the electric displacement along the transverse direction at each point, ρ denotes the mass density per unit volume, χ is the elastic stiffness, γ is the piezoelectric coefficient, µ is the magnetic permeability, and β is the water resistance coefficient of the beam and the prescribed voltage on electrodes of the beam. It is a well-known fact that a piezoelectric beam is an elastic beam, which has electrodes at its top and bottom surfaces, and in order to prevent fringing effects, it is insulated at the edges, as well as connected to an external electric circuit. The characteristics mentioned before are the simplest structures on which to study the interaction between the electrical and mechanical energy in these systems, see [Reference Kong, Nonato, Liu, Santos and Raposo10] for more information. Experimentally, it is observed that the magnetic effects are minor in the overall dynamics in polarized ceramics, and therefore these effects are ignored in the piezoelectric beam models. The piezoelectric beam models have a unique characteristic of converting mechanical energy into electrical and magnetic energy, and vice versa and, hence, are a competitive choice for many tasks in industry, particularly those involving control of structures. In fact, piezoelectric materials have been employed in civil, industrial, automotive, aeronautic, and space structures. The coupled system (1.1) has been studied from different points of views. Ramos et al. [Reference Ramos, Gonçalves and Correa Neto15], for example, inserted a mechanical dissipative term $\delta u_{t}$ into Equation (1.1)1 for χ > 0 and an appropriate boundary condition to prove that the energy decays exponentially. Later, Ramos et al. [Reference Ramos, Freitas and Almeida14] showed that under different boundary conditions, the system is exponentially stable, regardless of any relationship between system parameters, and also showed that the exponential stability is equivalent to an exact observability at the boundary. The topic has been thoroughly explained in [Reference Morris and Özer13, Reference Yang17, Reference Yang18] and supported by the multiple references therein.
On the other hand, fractional calculus is a mathematical tool which studies integral and differential operators of arbitrary order, not necessarily of integer order. For the past three decades or so, a number of scientists have shown a growing interest on the study of fractional calculus. Several engineering fields, applied sciences, and mathematical physics have benefited from this ascending wave of growing applications in this area. Moreover, space sciences, fluid mechanics, porous media flows, and viscoelastic and biological processes are areas in which fractional differential equations have become a favoured tool to seek new paths.
Fractional derivatives have been used in the model frequency–dependent damping behaviour of many viscoelastic materials, as well as in modelling many chemical processes. A diversity of problems from the scientific applied areas, including analysis of viscoelastic materials, heat conduction in materials with memory, electrodynamics with memory, signal processing, among others, can be modelled with fractional differential calculus. Models involving fractional derivatives represent an important number of natural phenomena in a more realistic way than models involving classical derivatives. For more information, we would like to direct the reader to [Reference Kilbas, Srivastava, Trujillo and van Mill9, Reference Mbodje11, Reference Samko, Kilbas and Marichev16] and the references therein.
The objective of this paper is to study polynomial stabilization in the case of $\mu_{1}(t) \equiv 0$ in Equation (1.1) by means of a boundary viscoelastic damper, the action of which is to cause a feedback force opposite the direction of the fractional derivative of the position of a limit point. By using an energy method, it is shown that the energy of the system decays polynomially. The method is motivated by the fact that the input–output relationship as generalized diffusion equation can be expressed in terms of fractional derivative operator and following the ideas given by [Reference Achori, Eddine and Benaissa1, Reference Benaissa and Abderrahmane3, Reference Mbodje11, Reference Mbodje and Montseny12]. In order to achieve results of polynomial stability, the dissipation given in the domain was eliminated and placed on the border. The imposed dissipation is of a fractional derivative type. Thus, we have focused on the following system:
where $\ 0 \lt \alpha \lt 1,$ $\eta \geq 0 $ and $u=u(x,t),$ $v=v(x,t)$ are real-valued functions.
Lemma 1.1. For every solution (u, v) of the system (1.2), the total energy ${\cal E}_{1}: \mathbb{R}^{+}\rightarrow \mathbb{R}^{+}$ is given at time t by
along with
and the positive definite quadratic form
Proof. Straightforward calculations.
Remark 1.2. We have seen that $\frac{{\rm d}{\cal E}_{1}(t)}\,{{\rm d}t} \neq 0,$ which follows from Equation (1.3), and it is not energy conserving. Moreover, due to the negativity of Equation (1.3), the system (1.2) is a dissipative kind of energy.
The interest in studying the fractional differential equations arose as a separate topic about 40 years ago. It has been of special interest to prove the existence of solutions of Cauchy-type problems involving the Caputo fractional derivatives. Several articles published during the last decade have developed engineering applications in viscoelastic damping and structural mechanics (see [Reference Caputo5, Reference Mbodje11]). A systematic and rigorous study of problems of this kind involving fractional differential equations can be found in recent literature; for a brief review, refer to [Reference Kilbas and Trujillo7, Reference Kilbas and Trujillo8]. Although numerous theoretical applications of fractional calculus operators have been found during an extended period time, many mathematicians and applied researchers have tried to model real processes using fractional differential equations. However, achieving results has been complex due in part to the fact that many of the properties of the ordinary derivative do not translate analogously to the case of the fractional derivative operator $D^{\alpha},$ such as, for example, a clear geometric or physical meaning, product rules, and chain rules.
In this work, we have obtained polynomial stability for a piezoelectric beam with a boundary dissipation of the fractional derivative type. In order to obtain the partial result, we used a similar approach to Mbodje and Montseny [Reference Mbodje and Montseny12]. This relationship help us achieve our main result in the following:
Theorem 1.3. The semigroup $\{{\cal S}_{\cal A}(t)\}_{t\geq 0}$ is polynomially stable and
This paper is organized as follows. In $\S$ 2, we reformulate the model (1.2) into an augmented system, which couples the system (1.1) with a suitable diffusion equation. In $\S$ 3, we establish the well-posedness of the system (2.10). In $\S$ 4, we show the polynomial stability of the corresponding semigroup.
Throughout this paper, C is a generic constant, not necessarily the same at each occasion. It will change from line to line and depends on the indicated quantities.
2. Augmented Model
This section is concerned with the reformulation of the model (1.2) with boundary conditions (1.2) $_{5,\,6,\,7,\,8}$ into an augmented system. To fulfill this goal, we will use a method based on the input–output relationship in diffusion equations which realize the fractional derivative operator. This tool in fractional calculus allows integrals and derivatives of any positive order. This can be considered a branch of mathematical physics which deals with integro-differential equations, where integrals are of a convolution type and exhibit weakly singular kernels of power law type. There are multiple ways in which fractional derivatives and integrals can be defined. These are not all equivalent to each other, but each of them has its own advantages and disadvantages.
Let $0 \lt \alpha \lt 1.$ We define the Caputo fractional integral of order α > 0 as
where Γ is the gamma function, $f \in L^{1}(0,t),$ and t > 0.
The Caputo fractional derivative operator of order α > 0 is defined by
with $f \in W^{1,1}(0,t)$ and t > 0. If the function f(t) represents the strained history within a viscoelastic material whose relaxation function is $[\Gamma(1 - \alpha)t^{\alpha}]^{-1}$, then the material will experience at any time t a total stress given the expression $D^{\alpha}f(t).$ Also, note that D α is a left inverse of $J^{\alpha},$ but in general, it is not a right inverse. More precisely, we have
For the proof of the above mentioned equalities and more properties of fractional calculus, we suggest the reader to refer [Reference Samko, Kilbas and Marichev16].
We consider slightly different versions of (2.1) and (2.2). For this work, we have focused on [Reference Choi and Maccamy6], where Choi and MacCamy establish the following definition of fractional integro-differential operators with exponential weight.
Let $0 \lt \alpha \lt 1, \eta \ge 0.$ The exponential fractional integral of order α is defined by
with $f \in L^{1}(0,t)$ and t > 0. The exponential fractional derivative operator of order α is given by
where $f \in W^{1, 1}(0,t)$ and t > 0. Note that the two operators (2.2) and (2.3) differ just by their kernels. $D^{\alpha,\eta}$ is merely Caputo’s fractional derivative operator with the exception of its exponential factor.
Theorem 2.1. ([Reference Mbodje11])
Let µ be the function
Then the relation between the Input U and the Output O of the following system
is given by $O = I^{1 - \alpha,\eta}U = D^{\alpha,\eta}U,$ where $\mathfrak{C} = \pi^{-1}\,\sin(\alpha \pi)$ and
Proof. Multiplying Equation (2.5) by ${\rm e}^{t\xi^{2}},$ integrating over $[0,t]$, and performing straightforward calculations, we have
Moreover, replacing Equation (2.4) into Equation (2.7), we get
To continue the process of making a clear relation between output and input, we turn Equation (2.8) into Equation (2.6) having Equation (2.9) as result:
With the change of variable $m = \xi^{2}(t - s)$ in the above equality, we are left with
and considering that $ \mathfrak{C} = \frac{1}{\Gamma(\alpha)\Gamma(1 - \alpha)}.$
Lemma 2.2. If $ \lambda \in D = \{\lambda\in \mathbb{C}: {\rm Re}\lambda + \eta \gt0\} \cup \{\lambda\in \mathbb{C}: {\rm Im}\lambda \neq 0\}. $ Then
Proof. See [Reference Mbodje11].
Remark 2.3. In what remains of this paper, we will write $\mathfrak{C} = \pi^{-1}\,\sin(\alpha \pi).$
Now, we reformulate system (1.2) using Theorem 2.1. Therefore, the system is included into the augmented model
Remark 2.4. The relation between the boundary conditions (1.2) and (2.10) are given by
Lemma 2.5. Let $(u, v, \varphi, \phi)$ be a solution of the system (2.10). Then, the energy functional is defined by
where ${\cal E}_{1}(t)$ and ${\cal E}_{2}(t)$ are the energies associated with the waves coupled systems and diffusion equations, respectively, and $q(u,\ v)$ is the positive definite quadratic form given in Equation (1.5), satisfying
The following theorems are fundamental to the proof of our results:
Theorem 2.6. ([Reference Arendt and Batty2])
Let ${\cal A}$ be the generator of a uniformly bounded C0-semigroup $\{{\cal S}(t)\}_{t\geq 0}$ on a Hilbert space $\cal H$. If:
i) ${\cal A}$ does not have eigenvalues on $\mathbb{R}.$
ii) the intersection of the spectrum $\sigma({\cal A})$ with $i\mathbb{R}$ is at most a countable set.
Then the semigroup $\{{\cal S}(t)\}_{t\geq 0}$ is asymptotically stable, that is,
Theorem 2.7. ([Reference Borichev and Tomilov4])
Let ${\cal A}$ be the generator of a uniformly bounded C0-semigroup $\{{\cal S}(t)\}_{t\geq 0}$ on a Hilbert space $\cal H.$ If
for some $\ell$, then there exist c such that
3. Setting of the Semigroup
In this section, we use results from the semigroup theory of linear operators in order to obtain an existence theorem of the system (2.10). Moreover, afterwards we will use the following standard $L^{2}(0,L)$ space, where the scalar product and the norm are denoted by
In a similar way, let $L^{2}(\mathbb{R})$ be the Hilbert space of all measurable square integrable functions on the real line with the inner product
We define $\displaystyle \mathbb{H}(0, L) = \{\Psi\in H^{1}(0, L):\ \Psi(0)=0\}.$ Then
equipped with the inner product given by
where $ U=(u, v, w, z, \varphi, \phi)^{\rm T}\ $ and $\widetilde{U} =(\widetilde{u}, \widetilde{v}, \widetilde{w}, \widetilde{z}, \widetilde{\varphi}, \widetilde{\phi})^{\rm T}.$ The norm is given by
In order to advance, we will transform the initial boundary value problem (2.10) to an abstract problem in the Hilbert space $\mathcal{H}.$ To achieve this, we will introduce the functions $u_{t} = w,$ and $v_{t} = z$ and rewrite the system (2.10) as the following initial value problem:
where $U=(u, v, w, z, \varphi, \phi)^{\rm T}\ $ and $\ U_{0}=(u_{0}, v_{0}, w_{1}, z_{1}, \varphi_{0}, \phi_{0})^{\rm T}, $ and the operator $\,\mathcal{A}:\mathcal{D}(\mathcal{A})\subset \mathcal{H}\rightarrow \mathcal{H}$ is given by
with domain
Proposition. 3.1. The operator ${\cal A}$ is the infinitesimal generator of the contraction semigroup $\{{\cal S}_{{\cal A}}(t)\}_{t\geq 0}.$
Proof. We will show that ${\cal A}$ is a dissipative operator and $0\in \varrho({\cal A})$ the resolvent set of ${\cal A}.$ Then the conclusion will follow Lumer–Phillips’ theorem. If $U = (u, v, w, z, \varphi, \phi)^{\rm T}\in {\cal D}({\cal A}),$ then from Equations (3.3) and (1.4) and using the energy ${\cal E}(t)$ definition, we get
In fact, let $U = (u, v, w, z, \varphi, \phi) \in {\cal D}({\cal A}),$ and using Equation (2.10) $_{9,10}$, we have
Subsequently, if we integrate by parts, we have
Replacing these integrals in $\langle {\cal A}U,\ U\rangle_{\cal H}$, we obtain
Grouping and simplifying the before mentioned, we are left with
Considering the boundary conditions (2.10) $_{9,\,10}$ and replacing them in the above equality, we have as a result
Consequently,
Taking the real part in the equality mentioned above, it follows properly that ${\cal A}$ is a dissipative operator and the statement (3.5) holds true.
Remark 3.2. From Equation (2.10) $_{9,\,10}$, we have
and equivalently $ A X_x = \varphi,$ where $\varphi= ( \varphi_1, \varphi_2)^{\rm T}$ and
If $\ \mathfrak{J} = \det(A) = \chi \beta - \gamma^{2} \beta^{2} \neq 0\ $, then
and X x is given by
We will prove that the operator $\lambda I - {\cal A}$ is surjective for λ > 0. For this, let $F =(f_1, f_{2}, f_{3}, f_{4}, f_{5}, f_{6})^{T}\in {\cal H},$ we seek $ U=(u, v, w, z, \varphi, \phi)^{\rm T}\in {\cal D}({\cal A}) $ such that $\lambda I - {\cal A}U=F,$ that is,
From Equation (3.8) $_{5,\,6}$, we make clear that
Note that, from Equation (3.8) $_{1,\,2}$, we have that $w, z \in \mathbb{H}(0, L)$, and by replacing w and z into Equation (3.8) $_{3,\,4}$, respectively, it yields that
To solve Equation (3.10), we need to find $u, v\in H^{2}(0, L)\cap \mathbb{H}(0,L)$ such that
for all $\zeta, \aleph \in \mathbb{H}(0,L).$ Integrating by parts Equation (3.11)1, we have
Now, using the same idea with Equation (3.11)2, we obtain
Moreover, using Equation (2.10) $_{9,\,10}$, we have
and
From Equation (3.9) $_{1, 2}$, we have, respectively,
and
Since $w(L, t) = \lambda u(L, t) - f_{1}(L) $ and $ z(L, t) = \lambda v(L, t) - f_{2}(L), $ replacing these functions into Equations (3.12) and (3.13), respectively, we obtain
and
The system in Equations (3.14)–(3.15) represents the equivalent to the problem
where the bilinear form continuous and coercive $\mathfrak{a} \left[\mathbb{H}^{1}(0, L)\times \mathbb{H}^{1}(0,L)\right]^{2} \rightarrow\mathbb{R}$ and the continuous linear form ${\cal L}: \mathbb{H}^{1}(0,L)\times \mathbb{H}^{1}(0,L) \rightarrow\mathbb{R}$ are defined by
and
By applying Lax–Milgram theorem, it follows that for all $(\zeta, \chi)\in \mathbb{H}^{1}(0, L)\times \mathbb{H}^{1}(0, L)$, the problem (3.16) admits a unique solution $(u, v)\in \mathbb{H}^{1}(0, L)\times \mathbb{H}^{1}(0, L).$ Using elliptic regularity, stated in Equations (3.14)–(3.15), it follows that $(u, v)\in \mathbb{H}^{1}(0, L)\times \mathbb{H}^{1}(0, L).$ Therefore, the operator $\lambda\,I - {\cal A}$ is surjective for any λ > 0.
Now, applying the Hille–Yosida theorem, we have the following result
Theorem 3.3. (a) If $U_{0}\in {\cal H},$ then the system (3.3) has a unique mild solution
(b) If $U_{0}\in {\cal D}({\cal A}),$ then the system (3.3) has a unique strong solution
4. Asymptotic behavior
In this section, we study the asymptotic behavior for the solution of our problem. First, we prove that the semigroup $\{{\cal S}_{{\cal A}}(t)\}_{t\geq 0}$ is asymptotically stable. To show this, we apply the Theorem 2.6 of Arend–Batty [Reference Arendt and Batty2]. Next, using Theorem 2.7 of Borichev–Tomilov, we show that $\{{\cal S}_{\cal A}(t)\}_{t\geq 0}$ is polynomially stable.
To achieve our first result, we need the following lemmas:
Lemma 4.1. λ = 0 is not an eigenvalue of ${\cal A}.$
Proof. By applying Equation (3.4), we have as a result $U = (u, v, w, z, \varphi, \phi)^{\rm T}\in \ker({\cal A})$ if and only if ${\cal A}U = 0,$ that is,
From Equation (4.1) $_{1, 2, 5, 6}$, we have $w = z = \varphi = \phi = 0.$ Multiplying Equation (4.1) $_{3, 4}$ by u and v, respectively, and integrating over $x\in (0,\,L)$ and using the definition of $\mathbb{H}(0, L)$, it follows that
Using the boundary conditions (2.10) $_{9, 10}$ and that $\varphi = \phi = 0$, we obtain
Adding Equation (4.3)1 and Equation (4.3)2, we obtain
We write this equality as
where $X^{\rm T} A X$ is the positive definite quadratic form given by Equation (1.5).
From Equation (4.4), we have found that $u_{x} = v_{x} = 0,$ which involves $u(x, t) = C,$ where C is a constant, but $u(0, t) = 0,$ then $u \equiv 0.$ Similarly, we obtain $v \equiv 0.$ Then we conclude that $(u, v, w, z, \varphi, \phi)^{\rm T}=0$ and the Lemma follows.
Lemma 4.2. It verifies that
Proof. We will prove the equality by contradiction. In order to do that, we suppose there exists $\lambda\in\mathbb{R},$ λ ≠ 0, and $U\neq 0,$ such that $\left(i \lambda - {\cal A}\right)U = 0.$ Then we have
According to Equation (3.5), we imply $\varphi(\xi,\,t) = 0$ and $\phi(\xi,\,t) = 0,$ and by replacing these in Equation (4.6) $_{5, 6}$, we obtain $ w(L, t) = z(L, t) = 0;$ hence, from Equation (4.6) $_{1, 2}$ it is implied that $u(L, t) = v(L, t) = 0.$ Moreover, by revisiting Equation (2.10) $_{9, 10} $ it follows that $u_{x}(L, t) = v_{x}(L, t) = 0.$
Now, replacing Equation (4.6) $_{1, 2}$ onto (4.6) $_{3, 4}$, we obtain
To clarify the conclusion, we rewrite Equation (4.7) as the initial value problem:
where $A, X$ are given by Equation (3.6) and A −1 by Equation (3.7) for $\mathfrak{J} = \chi \beta - \gamma^{2} \beta^{2} \not = 0,$ and
Using the Picard theorem for ordinary differential equations, the system in Equation (4.8) has a unique solution $(u, v) = 0$; of these and Equation (4.6) $_{1, 2}$, we have w = 0 and z = 0. Therefore, $(u, v, w, z, \varphi, \phi) = 0$, and we conclude that ${\cal A}$ does not have eigenvalues on $\mathbb{R}.$
In the Theorem 2.6, the condition (ii) holds if we show that $\sigma({\cal A})\cap \{i \mathbb{R}\}$ is at most a countable set.
In fact, we will show that the operator $i \lambda I - {\cal A}$ is surjective for λ ≠ 0. For this purpose, and given $F =(f_{1}, f_{2}, f_{3}, f_{4}, f_{5}, f_{6})^{\rm T}\in {\cal H},$ and $U=(u, v, w, z, \varphi, \phi)^{\rm T}\in {\cal D}({\cal A})$ such that $i \lambda I - {\cal A}U=F,$ that is,
with the initial condition
From Equation (4.10) $_{5, 6}$, we have
Note that from Equation (4.10) $_{1, 2}$, it follws that $w, z \in \mathbb{H}(0, L).$ Now, replacing Equation (4.10) $_{1, 2}$ into Equation (4.10) $_{3, 4}$ we get, respectively,
Solving system (4.13) is equivalent to finding $u, v\in H^{2}(0, L)\cap \mathbb{H}(0, L)$ such that
for all $V = (\zeta, \aleph) \in \mathbb{H}(0, L)\times \mathbb{H}(0, L).$ Integrating Equation (4.14)1 by parts
If we replace Equation (2.10)9 in the above equality, we will get
Now, using Equation (4.12)1, it follows that
By replacing Equation (4.10)1 onto Equation (4.15), we have
Similarly to what we find in Equation (4.16), we obtain, starting from Equation (4.14)2 that
The system (4.16)–(4.17) is equivalent to the problem
with
and in addition to
Using that $L^{2}(0, L) \stackrel{c}{\hookrightarrow} H^{-1}(0, L)$ and $\mathbb{H}(0, L)\stackrel{c}{\hookrightarrow} L^{2}(0, L)$, it follows that the operator ${\cal L}_{\lambda}$ is compact from $L^{2}(0, L)$ into $L^{2}(0, L).$ By using Fredholm alternative, to prove the existence of solution, Equation (4.18) is reduce to show that 1 is not a eigenvalue of ${\cal L}_{\lambda}.$ In fact, by contradiction, if 1 is an eigenvalue, then there exists $U\neq 0,$ such that
In particular, for U = V, we obtain:
Thus, by the above equation, the imaginary terms are equal to zero, then we should have
Consequently, in Equation (4.20), we have
Furthermore, from Equation (4.19), integrating by parts and performing straightforward calculations, we deduce that $ u_{x}(L, t) = 0 \qquad \mbox{and}\ v_{x}(L, t) = 0;$ and, from this premise Equation (4.14), we sustain
Now, we rewrite Equation (4.23) as the initial value problem
where A, X, and A −1 are defined in Remark 3.2, and B has been given by Equation (4.9).
Due to the Picard theorem for ordinary differential equation, we arrive at Equation (4.24), which has a unique solution u = 0 and v = 0. Moreover, from Equation (4.6) equality, we obtain w = 0 and z = 0. Therefore $U=(u, v, w, z, \varphi, \phi)=0.$
Lemma 4.3. If $\eta \neq 0,$ we have that $\ 0\in\varrho({\cal A}).$
Proof. We have to prove that given $F=(f_{1}, f_{2}, f_{3}, f_{4}, f_{5}, f_{6})\in {\cal H}$, there exists a unique $U=(u, v, w, z, \varphi, \phi)$ in ${\cal D}({\cal A}),$ such that $-{\cal A}U=F$, namely
with the initial condition
To achieve this, we use Equations (4.16) and (4.17) with $\lambda = 0,$ and it follows that
and
The system in Equations (4.26) and (4.27) is equivalent to the problem
where the bilinear form is continuous and coercive $\mathfrak{a}_{\eta}: \left[\mathbb{H}(0,\,L)\times \mathbb{H}(0, L)\right]^{2}\rightarrow\mathbb{R},$ and the continuous linear ${\cal L}_{\eta}: \left[\mathbb{H}(0, L)\right]^{2}\rightarrow\mathbb{R}$ are defined by
and
Applying the Lax–Milgran theorem, we have that for all $(\zeta, \aleph)\in \mathbb{H}(0, L)\times \mathbb{H}(0, L)$, the problem (4.28) admits a unique solution $(u, v)\in \mathbb{H}(0, L)\times \mathbb{H}(0, L).$ Then by using elliptic regularity, stated in Equations (4.26) and (4.27), it is deduced that $(u, v) \in H^{2}(0, L)\times H^{2}(0, L).$ Therefore, the operator ${\cal A}$ is surjective.
Lemma 4.4. Let ${\cal A}$ be defined by Equation (3.4), then
with the domain
Proof. It is a straightforward verification to show that $\ \langle {\cal A}G, G\rangle = \langle G, {\cal A}^{*}G\rangle.\ $ In fact, let $G = (u, v, w, z, \varphi, \phi)$ and $\widetilde{G} = (\widetilde{u}, \widetilde{v}, \widetilde{w}, \widetilde{z}, \widetilde{\varphi}, \widetilde{\phi}).$ Then
If we integrate by parts the first eight terms, grouping and replacing on before the equality, we obtain
Revisiting the boundary conditions (2.10) $_{9, 10},$ we have
where ${\cal A}^{*}$ is given in Equation (4.29). Lemma 4.4 follows.
Remark 4.5. An operator may be injective and bounded below, but not necesarily invertible, for example, the unilateral shift on $l^{2}(\mathbb{N}).$ This shift operator is isometry and bounded inferiorly, but it is not invertible since it is not surjective. The set of λ for which $\lambda\,I - {\cal A}$ is injective but does not have a dense range is known as the residual spectrum of ${\cal A}$ and is denoted by $\sigma_{r}({\cal A}).$
Theorem 4.6. The set of residual spectrum verifies that $\sigma_{r}({\cal A}) = \emptyset.$
Proof. We consider that the eigenvalues of ${\cal A}$ are symmetric on the real axis; hence, we will show that $\sigma_{r}({\cal A}) = \sigma_{p}({\cal A}^{*}).$
In fact, let $\lambda\in\mathbb{C}$ such that ${\cal A}^{*}U = \lambda\,U$ and $0\neq U=(u,\,v,\,w,\,z,\,\varphi,\,\phi)\in \mathcal{D}({\cal A}^{*}).$ We recall the definition (4.29) and therefore we have
Replacing Equation (4.30) $_{1, 2}$ into (4.30) $_{3, 4}$, we obtain
with the boundary condition
The system (4.31) and (4.32) matches with the eigenvalue problem of ${\cal A}.$ We conclude that ${\cal A}^{*}$ has the same eigenvalues as that of ${\cal A}$ and the theorem follows.
Considering the previous lemmas and utilizing Theorem 2.6, we have proved the following:
Theorem 4.7. The $C_0-$semigroup of contraction ${\cal S}_{{\cal A}}(t)_{t\geq 0}$ is asymptotically stable, which means
Therefore, we are able to prove the second main theorem of this section.
Theorem 4.8. Let η ≠ 0. The semigroup ${\cal S}_{{\cal A}}(t)_{t\geq 0}$ is polynomially stable and
Proof. We will analyze the resolvent equation $(i \lambda I - {\cal A})U = F,$ with $\lambda\in\mathbb{R},$ which in length is
where $F = (f_{1}, f_{2}, f_{3}, f_{4}, f_{5}, f_{6}).$ Revisiting Equation (3.5), we have
that is,
Furthermore, using $\left||a| - |b| \right| \leq \left|a - b\right|$ into Equation (4.34) $_{1,\,2}$ we have, respectively,
then
Thus,
Subsequently, from Equation (2.10) $_{9, 10}$ and by using the Cauchy–Schwartz inequality, we have
as well as
Ultimately, by considering Remark 3.2, we obtain
Moreover, from Equation (4.34) $_{5, 6}$, we can derive, respectively, that
An by multiplying Equation (4.36)1 by $(\xi^{2} + \eta + i \lambda)^{-1} \mu(\xi)$, we have
Then by applying absolute values and integrating them over $\xi\in \mathbb{R}$, we have the following result:
Once again, by using the Cauchy–Schwartz inequality and having straightforward estimates in front, it follows that
Also by applying power square on both sides of the inequality, as well as using $2 a b \leq a^{2} + b^{2}$, we are able to obtain
Hence, by using Equation (4.35) $_{1},$ we get
Applying Lemma 2.2, it follows that
In a similar way, we extrapolate the term (4.36) $_{2},$ resulting in
Consequently, we get
To conclude the verification of the theorem, we need the following lemma:
Lemma 4.9. Let $q\in H^{1}(0, L).$ Therefore, we have
where $\displaystyle R_{1} = 2 \rho\int_{0}^{L}q {\rm Re}(f_{3} \overline{u}_{x})\, {\rm d}x + 2 \rho\int_{0}^{L}q {\rm Re}(\overline{f}_{1x} w)\, {\rm d}x $ and also
along with $\displaystyle R_{2} = 2 \mu\int_{0}^{L}q {\rm Re}(f_{4} \overline{v}_{x})\, {\rm d}x + 2 \mu\int_{0}^{L}q {\rm Re}(\overline{f}_{2x} z)\, {\rm d}x. $
Note that for $R_{i},$ $(i = 1, 2)$, we have
Proof. To demonstrate Lemma 4.9, we will multiply Equation (4.34)3 by $q \overline{u}_{x}$ and integrate over $(0, L)$, and we have that
From Equation (4.34)1, we have $\overline{i \lambda u_{x}} = \overline{w}_{x} + \overline{f}_{1x}.$ Then
Thus,
Taking the real part of the aforementioned equality, it follows that
Using the precedent of ${\rm Re}(\Phi \overline{\Phi}_{x}) = \frac{1}{2} \frac{\rm d}{{\rm d}x}|\Phi|^{2},$ we have
By integrating by parts the first two integrals, we get
It continues as
where $R_{1} = 2 \rho\int_{0}^{L}q {\rm Re}(f_{3} \overline{u}_{x})\, {\rm d}x + 2 \rho\int_{0}^{L}q {\rm Re}(\overline{f}_{1x} w)\, {\rm d}x.$
Therefore, performing similar calculations onto Equation (4.39), we obtain
with $ R_{2} = 2 \mu\int_{0}^{L}q {\rm Re}(f_{4} \overline{v}_{x})\, {\rm d}x + 2 \mu\int_{0}^{L}q {\rm Re}(\overline{f}_{2x} z)\, {\rm d}x. $
The lemma is corroborated.
Now, we continue with the proof of Theorem 4.8. Taking $q(x) = x$ in Lemma 4.9, we get
where
and
with
Moreover, integrating by parts and bounding the following integral,
By adding Equation (4.41) with Equation (4.42) and using Equations (4.43) and (4.38), we have that
where $\chi \geq \gamma \beta$ and $1\geq \gamma.$ Furthermore, from Equation (4.35), we are able to confirm
Then by adding Equations (4.44) and (4.45), it follows that
From Equation (4.34) $_{1, 2}$, we also have
and
Replacing Equations (4.47) and (4.48) onto Equation (4.46), we get
If $|\lambda|\gt1$, we have
Consequently, it follows that
for a positive constant $C.$ The conclusion is obtained by applying Theorem 2.7.
Funding Statement
V.P. and O.V. are partially financed by project Fondecyt 1191137.