We consider an Klein–Gordon relativistic equation with a boundary dissipation of fractional derivative type. We study of stability of the system using semigroups theory and classical theorems over asymptotic behavior.
The classical works of Klein (1927) and Gordon (1926) derived a relativistic equation for a charged particle in an electromagnetic field. This equation is of conservative dispersive type and has played an important role in the study of elementary particles. In general, the Klein–Gordon equation is used to describe dispersive wave phenomena. In a series of papers [5,6], Mainardi has considered the Klein–Gordon equation with dissipation. These works are relevant for provide interesting examples of normal dispersion (usually met in the absence of dissipation) and anomalous dispersion (always present on viscoelasticity waves). In this work, we consider the following problem
This equation can be thought as a simple prototype of two regimes of dispersion that depend on the relative weight of their two independent parameters γ and β. For , the equation is the linear version of the Klein–Gordon equation, and when . the equation gives the telegraph equation. Moreover, if the occurrence of normal dispersion is observed and, the case the anomalous dispersion is obtained, in fact, this case turns out to be similar to viscoelastic waves in view of the common anomalous dispersion in their full range of frequencies. If , this equation corresponds to the distortion-less wave propagation, where there is attenuation with no dispersion. For the special case , it reproduces the Maxwell model of viscoelasticity.
In [8], Mbodje and Montseny shown that the fractional derivative forces the system to become dissipative and the solution approaches the equilibrium state. Hence, when dissipation is applied on the boundary, we can consider them as controllers helping to reduce vibrations.
In this article we are interested in studying the stabilization in the case in (1.1). In order to achieve this, we have eliminated the dissipation given in the domain and placed the dissipation on the border. It is importante to note that dissipation imposed is of a fractional derivative type and thus we focused on the following system:
where is a real-value function, denote the order of the exponential fractional derivative, and the weight in the exponential factor of this derivative (see (2.4)).
In recent decades, the equations that involving real order derivatives have assumed an important role in modeling the anomalous dynamics of many processes related to complex systems in the most diverse areas of sciences and engineering. The interest in the specific topic of fractional calculus surged only at the end of the last century, however, the study of fractional differential equations as a separate topic arose some 40 years ago. Questions about the existence of solutions to Cauchy type problems involving the Caputo fractional derivatives and integral. Several articles published during the last decade have developed engineering applications in viscoelastic damping and structural mechanics (see [7] and references therein). A systematic and rigorous study of some problems of this kind involving fractional differential equations can be found in recent literature [4]. Although numerous theoretical applications of fractional calculation operators have been found during their long history, many applied researchers have tried to model processes using this theory but it has not been easy. This may be due, in part, to the fact that many of the useful properties of the entera derivative are not transferred analogously to fractional case, such as, for example, a clear geometric or physical meaning, product rules, chain rules, and so on.
In this work, we have obtained the polynomial stability of a thermoelastic mixture with a boundary dissipation of fractional derivative type. In order to obtain the partial result we used a similar approach that Mbodje and Montseny [8], which is concerned on the fact that the input-output relationship in a certain diffusion equation realizes the fractional derivative operator.
This paper is organized as follows. In Section 2 briefly outlines the notation and also we reformulate the model (1.2) into an augmented system, coupling the system (1.2) with a suitable diffusion equation. In Section 3, we establish the well-posedness of the system (2.9). In Section 4 we show the lack of exponential stability. In Section 5 we obtain the polynomial stability of the corresponding semigroup.
Finally, 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.
Augmented model
This section is concerned with the reformulation of the model (1.2) with boundary conditions (1.2)3, 4 into an augmented system. Method based on the fact the input-output relationship in a certain diffusion equation realizes the fractional derivative operator. We review some basic properties of fractional integrals and derivatives, which we will use later in the analysis of our problem. Fractional calculus, in allowing integrals and derivatives of any positive order, can be considered a branch of mathematical physics which deals with integro-differential equations, where integrals are of convolution type and exhibit weakly singular kernels of power law type. There are multiplies 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 . We define the Caputo fractional integral of order as
where Γ is the well-known gamma function and , and .
Moreover, the Caputo fractional derivative operator of order is defined by
with , and . Note that Caputo definition of fractional derivative does possess a very simple but interesting interpretation: if the function represents the strain history within a viscoelastic material whose relaxation function is then the material will experience at any time t a total stress given the expression . Also, it easy to show that is a left inverse of , but in general it is not a right inverse. More precisely, we have
For the proof of above equalities and more properties of fractional calculus see [10].
In this work, we consider slightly different version of (2.1) and (2.2). In [3], Choi and MacCamy establish the following definition of fractional integro-differential operators with weight exponential. Let , , the exponential fractional integral of order α is defined by
with and . The exponential fractional derivative operator of order α is defined by
with and . Note that .
Very little attention has been given to this type of non-local fractional feedback. In this article, we work with fractional derivatives that involve an integrable kernel. However, in future works, we are interested in considering feedback that involve singular and non-integrable kernels, this leads to substantial mathematical difficulties because the methods developed for terms of convolution with integrable kernel are no longer valid.
Let u be a solution of the system (
1.2
). Then, the energy functionalis given byand satisfies
In fact, multiplying (1.2)1 by , integrating by part in x over [0,1] and performing straightforward calculations we obtain (2.5), deriving with respect to t we have
Replacing (1.2)4 the result follows. □
Now, we given a reformulation of the model (1.2) with boundary conditions (1.2)3, 4 into an augmented system. For this, we introduce the equation
where . Multiplying the above equation by and integrating we get
and
using Fubini theorem and recalling the definition of the Gamma function we get that
For more details of this deduction see [7]. The next technical lemma will be very useful for our purposes.
With the representation (2.8) for the fractional derivative, the system (1.2) can be written as
From (2.8), the boundary conditions in (1.2) and (2.9) note that .
The energy associated with to previous system is given by
and satisfies
Setting of the semigroup
In this section we use results of the semigroup theory of linear operators to obtain an existence theorem of the system (2.9). We will use the following standard space, we are the scalar product and the norm are denoted by
In a similar way, let be the Hilbert space of all measurable square integrable functions on the real line with the inner product
We define . Then
equipped with the inner product given by
where and . The norm is given by
Now, we wish to transform the initial boundary value problem (2.9) to an abstract problem in the Hilbert space . We introduce the functions and rewrite the system (2.9) as the following initial value problem
where and , and the operator is given by
with domain
For the linear operator we define the resolvent set by
and the spectrum of is .
Note that an operator may be injective, even bounded below, but not invertible. The unilateral shift on is such an example. This shift operator is an isometry, therefore bounded below by 1. But it is not invertible as it is not surjective. The set of λ for which is injective but does not have dense range is known as the residual spectrum or compression spectrum of and is denoted by .
The operatoris the infinitesimal generator of a contraction semigroup.
We will show that is a dissipative operator and that the operator is surjective for any . Then our conclusion follows using the well known the Lumer–Phillips theorem. We observe that if , then using (2.11) and (3.3), we get
In fact,
Integrating by parts we have
Taking the real part, it follows that is a dissipative operator and (3.5) follows.
Next, we will prove that the operator is surjective for . Let and such that , that is,
From (3.6)3 we have
and by (3.6)1 we obtain
Replacing (3.8) into (3.6)2 yields
To solve (3.9) is equivalent to finding such that
for all . Then
Integrating by parts we have
Using (2.9)5 it follows that
From (3.7) we have
Then
Replacing (3.8) into (3.11) we obtain
This equality is equivalent at the problem
where the bilinear form continuous and coercive and the continuous linear form are defined by
and
It follows by the Lax–Milgram theorem that for all the problem (3.12) admits a unique solution . Using elliptic regularity, it follows from (3.12) that . therefore, the operator is surjective for any . □
Now, using the Hille–Yosida theorem, we have the following results.
If, then the system (
3.3
) has a unique mild solution
If, then the system (
3.3
) has a unique strong solution
The lack of exponential stability
In this section, we show the lack of exponential stability. For this, the following results are fundamental.
Letbe a-semigroup of contractions on Hilbert spacewith generator. The following statements are equivalent:
The semigroupis exponentially stable:
The setand
Equipped with the foregoing theorem, we going to prove our claim.
The semigroupgenerated by the operatoris not exponentially stable.
Case. We will show that .
In fact, we consider , then denoting by the image of by , we can see that . But for . This way .
We are going to conclude the proof by contradiction. Suppose that , this implies that is invertible and for all . From the equality
we have
for all . If then the left side of the equality belongs to and the right side does not belong to which it is a contradiction. This implies that is not in the resolvent set of .
Case. Using the Theorem 4.2, it suffices to prove that is unbounded for .
In fact, let , , and be such that:
Replacing (4.1)1 into (4.1)2 and performing straightforward calculations we have
Writing for a good choice of β. Let , we have the equality
The solution to the homogeneous problem associated with (4.2) is given by
Using (4.1)4 it follows that
Since the wronskian is , we have that the particular solution of equation (4.2) is
where
Hence
From (4.3) and (4.6) we obtain that the general solution of (4.2) is given by
Moreover,
and
From (4.1)3 we have , hence
Replacing (4.7), (4.9) and (4.10) into (4.1)3 we obtain
where
Let be the sequence defines by
with , . Remembering that and solving ω we have , hence if .
Now, replacing (4.13) into (4.6) and deriving, we can to obtain
Hence
Thus,
On the other hand, for all and , we have
Then, into (4.15)
where
integrating respect to x on and remembering that , , we obtain
and the following conditions holds:
Moreover
which, together with (4.17), implies that
Replacing (4.13) in (4.12) and (4.11) and using (4.16)
which combined with Lemma 4.1 we obtain that as . This fact together with (4.19), (4.18) and (4.14), we reach to
the result follow of Theorem 4.2. □
Asymptotic behavior
In this section, we show the asymptotic behavior of the solution of problem. First, we prove that the semigroup is asymptotically stable, for this we use Theorem 5.1 due to Arend–Batty [1]. Next, we show that is polynomially stable using the Theorem 5.7 of Borichev–Tomilov, see [2].
To achieve our first proof we use the following result.
Letbe the generator of a uniformly bounded C0-semigroupon a Hilbert space. Ifis at most a countable set andthen the semigroupis asymptotically stable, that is,
In order to obtain our result, we will test the following lemmas.
We have.
By contradiction. We suppose that there , and , such that , that is, . Then
From (3.5) we have , and by (5.1)3 we obtain . Hence, from (5.1)1 and (2.9)5 it follows that and . On the other hand, replacing (5.1)1 into (5.1)2 we obtain
We can rewrite (5.2) as the initial value problem
In the Theorem 5.7, the condition (2) holds if we show that any point is at most a countable set.
In fact, we will prove that the operator is surjective for . For this purpose, let , we seek such that , that is,
with the initial condition
Suppose that we have v, from (5.4)1 it follows that . Thus, . On the other hand, replacing (5.4)1 into (5.4)2 we get
this equation is equivalent to find such that
for all . We estimate (5.6),
Clearing ϕ in (5.4)3 and using (2.9)3, we get
From (5.4)1, we have
Rewrite this equality as
with
and .
Using that and it follows that the operator is compact from into . This way, by Fredholm alternative, proving the existence of U solution of (5.8) reduces to show that 1 is not a eigenvalue of . In fact, if 1 is an eigenvalue, then there exists , such that
In particular, for we obtain
and, of this .
Moreover, from (5.9), integrating by parts and performing straightforward calculations we can to deduce that , and from (5.6) we have
We can rewrite the above as the initial value problem
Using the Picard theorem (ordinary differential equations), (5.10) has a unique solution . It follows from (5.1) that . Therefore, . □
Ifthen.
From (3.4) we have that if and only if , that is,
From (5.11) we have . Multiplying (5.11)2 by u, integrating over and using the definition of it follows that
Moreover, using that into (2.9)5 we obtain . This way, (5.12) is given by
Hence and and is injective.
Given , we must show that there exists a unique in , such that , namely,
with the initial condition
From (5.7) it follows that
The equation (5.15) is equivalent to the problem
where the bilinear form continuous and coercive and, the continuous linear form are defined by
and
Applying the Lax–Milgran theorem, we have that for all the problem (5.16) admits an unique solution . Using elliptic regularity, it follows from (5.15) that . Therefore, the operator is surjective. □
Since , the proof will be successful if we can show that . This is because we have considering that the eigenvalues of are symmetric on the real axis. In fact, we will consider the eigenvalue problem for and , that is, from (5.17)
Replacing (5.18)1 into (5.18)2 we obtain
with the following boundary condition
However, the system (5.19)–(5.20) is exactly the eigenvalue problem of . Thus, has the same eigenvalues with . The lemma follows. □
With the previous lemmas and using the Theorem 5.1 we have proved that.
The-semigroup of contractionis asymptotically stable, that is,
To achieve our second proof we use the following result.
Letbe a C0-semigroup of contractions on Hilbert space. Iffor some ℓ, then there exist c such that
Now, we will prove the second main theorem of this section.
Ifthen the semigroupis polynomially stable and
We will study the resolvent equation , . That is,
where . Taking the inner product in with U and using (3.5) we have
that is,
Moreover, from (5.22)1 we have , then
On the other hand, from (2.9)5 and using the Cauchy-Schwartz inequality we have
Hence
From (5.22)3 we obtain
Multiplying this equality by we have
Applying absolute values and integrating over we have
Applying the Cauchy-Schwartz inequality and straightforward estimates it follows that
Applying power squared on both sides of the inequality and using we obtain
Hence, using (5.23)
From Lemma 2.2, it follows that
To conclude the proof of the theorem we need the following lemma:
Let. Then we havewhere
For R, we have
Proof of Lemma Multiplying (5.22)2 by and integrating over we have
From (5.22)1 we have . Then
Thus
Then integrating by parts
Taking the real part it follows that
On the other hand, . Hence into (5.32) we have
Then
It follows that
and, the lemma follows. □
Now, we return to the proof of the Theorem 5.8. Taking in Lemma 5.9 we have
Using the Young inequality and (5.25) we get
Moreover, from (5.23)
Hence, adding (5.35) and (5.36) it follows that
From (5.22)1 for , we have . Then for ,
If we get . It follows that
for a positive constant C. The conclusion follows applying the Theorem 5.7. □
Footnotes
Acknowledgements
The first author was partially supported by CNPq-Brazil: grants 158706/2014-5, 164793/2015-1 and 402689/2012-7. The authors were partially financed by project Fondecyt 1191137.
References
1.
W.Arendt and C.J.K.Batty, Tauberian theorems and stability of one-parameter semigroups, Trans. Amer. Math. Soc.306(2) (1988), 837–852. doi:10.1090/S0002-9947-1988-0933321-3.
2.
A.Borichev and Y.Tomilov, Optimal polynomial decay of functions and operator semigroups, Math. Ann.347(2) (2010), 455–478. doi:10.1007/s00208-009-0439-0.
3.
J.Choi and R.Maccamy, Fractional order Volterra equations with applications to elasticity, J. Math. Anal. Appl.139 (1989), 448–464. doi:10.1016/0022-247X(89)90120-0.
4.
A.Kilbas and J.Trujillo, Differential equation of fractional order: Methods, results and problems, I Appl. Anal.78(2) (2002), 435–493, II Appl. Anal., 81, 1-2, 2001, 153–192. doi:10.1080/00036810108840931.
5.
F.Mainardi, On signal velocity of anomalous dispersive waves, Il Nuovo Cimento. B74 (1983), 52–58. doi:10.1007/BF02721684.
6.
F.Mainardi, Linear dispersive waves with dissipation, in: Wave Phenomena: Modern Theory and Applications, C.Rogers and T.B.Moodie, eds, North-Holland, Amsterdam, 1984, pp. 307–4317. doi:10.1016/S0304-0208(08)71274-X.
7.
B.Mbodje, Wave energy decay under fractional derivative controls, IMA Journal of Mathematical Control and Information.23 (2006), 237–257. doi:10.1093/imamci/dni056.
8.
B.Mbodje and G.Montseny, Boundary fractional derivative control of the wave equation, IEEE Trans. Autom. Control.40 (1995), 368–382.
9.
J.Prüss, On the spectrum of C0-semigroup, Trans. Am. Math. Soc.248 (1984), 847–867. doi:10.2307/1999112.
10.
S.Samko, A.Kilbas and O.Marichev, Integral and Derivatives of Fractional Order, Gordon Breach, New York, 1993.