Next Article in Journal
Investigation and Application of Fractals in Civil Engineering Materials
Next Article in Special Issue
Determination of a Nonlinear Coefficient in a Time-Fractional Diffusion Equation
Previous Article in Journal
Combined Liouville–Caputo Fractional Differential Equation
Previous Article in Special Issue
Compact Difference Schemes with Temporal Uniform/Non-Uniform Meshes for Time-Fractional Black–Scholes Equation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Stabilization of a Non-linear Fractional Problem by a Non-Integer Frictional Damping and a Viscoelastic Term

by
Banan Al-Homidan
and
Nasser-eddine Tatar
*
Interdisciplinary Research Center for Intelligent Manufacturing & Robotics, Department of Mathematics, King Fahd University of Petroleum and Minerals, Dhahran 31261, Saudi Arabia
*
Author to whom correspondence should be addressed.
Fractal Fract. 2023, 7(5), 367; https://doi.org/10.3390/fractalfract7050367
Submission received: 2 January 2023 / Revised: 13 April 2023 / Accepted: 18 April 2023 / Published: 28 April 2023

Abstract

:
In this paper, we consider a fractional equation of order between one and two which may be looked at as an interpolation between the heat and wave equations. The problem is non-linear as it involves a power-type non-linearity. We investigate the possibilities of stabilizing the system by a lower-order fractional term and/or a memory term involving the Laplacian. We prove a global Mittag–Leffler stability result in case a fractional frictional damping is active and a local Mittag–Leffler stability result when the material is viscoelastic in case of small relaxation functions. Unlike the integer-order problems, additional serious difficulties arise in the present case. These difficulties are highlighted clearly in the introduction. They are mainly due to the memory dependence of the fractional derivatives which is the cause of the invalidity of the product rule in particular. We utilize several properties in fractional calculus. Moreover, we introduce new Lyapunov-type functionals in the context of the multiplier technique.

1. Introduction

For many phenomena in engineering and science, fractional calculus offers an appropriate platform for their studies [1,2,3,4,5,6,7]. In fact, for anomalous processes and disordered systems, where the mean square displacement MSD (measure of the deflection of the actual position of a particle with regards to a given reference position) grows in a non-linear fashion, one is forced to consider fractional derivatives. Complex diffusion phenomena are therefore better described by non-integer-order fractional derivatives.
The problem to be discussed here is the following non-linear fractional initial boundary value problem of order between 1 and 2
C D 2 α v ( t , x ) + a C D α v ( t , x ) + v γ v = Δ v ( t , x ) b 0 t k ( t s ) Δ v ( s , x ) d s , t > 0 , x Ω , 1 / 2 < α < 1 , v ( t , x ) = 0 , t > 0 , x Ω , v ( 0 , x ) = v 0 ( x ) , v ( 0 , x ) = v 1 ( x ) , x Ω
where C D α is the Caputo fractional derivative of order α defined below, k is a non-negative function, v 0 ( x ) and v 1 ( x ) are the initial data of the state, a , b 0 will be equal to 1 or 0 depending on whether they are effective or not in the equation, γ > 0 is to be determined later and Ω is a bounded domain in R n with smooth boundary Ω .
It is well known that the relationship between the heat flux vector and the temperature gradient is given by
q = k T
where k is the thermal conductivity. Then, having in mind the conservation of energy
γ T t = div q
where γ is the heat capacity, the basic heat equation
T t = k γ Δ T
is obtained. Taking into account the history dependence between the gradient and the flux
q = k Γ ( μ ) t 0 t ( t s ) μ 1 T ( s ) d s = k R L D μ T , 0 < μ < 1 ,
we are lead to the fractional heat conduction equation
C D α T = δ Δ T
for 0 < α < 1 and δ > 0 (see [5]).
A sub-diffusive problem has been derived in [8] from the Rayleigh–Stokes problem
w t = ( 1 + γ R L D β ) Δ w .
The authors proved that the solutions of Stoke’s first problem follow as a special case.
In case of a material with memory, the relationship between the heat flux vector and the temperature gradient will involve a term of the form
0 t g ( t s ) T ( s ) d s .
The fractional equation
C D 2 α v ( t , x ) = Δ v ( t , x )
is a fractional version interpolating the heat and wave equations. It describes anomalous processes and corresponds to the (fast) super-diffusion (particles move faster than in the normal case) (see [9,10]) when the mean square displacement grows non-linearly in time ( < v ( t ) > ~ t 2 α , 1 < 2 α < 2 ). The importance of super-diffusion is explained and highlighted in many works, see [3,6,7,11] to cite but a few.
A linear growth of the MSD ( < v ( t ) > ~ t ) corresponds to the ordinary case (normal diffusion), when ( < v ( t ) > ~ t α , 0 < α < 1 ), we are in presence of sub-diffusion (particles move slower than in the ordinary case), whilst α = 2 corresponds to the ballistic diffusion. Apart from the normal diffusion and the ballistic diffusion, the other cases cannot be described adequately by the standard heat or wave equations. Therefore, accordingly, the equation
C D 2 α v ( t , x ) = Δ v ( t , x ) b 0 t k ( t s ) Δ v ( s , x ) d s , 1 / 2 < α < 1 ,
interpolates the heat conduction equation [12]
u ( t ) = Δ u ( t ) 0 t k ( t s ) Δ u ( s ) d s
and the viscoelastic equation [13]
w ( t ) = Δ w ( t ) 0 t k ( t s ) Δ w ( s ) d s .
The stability of these problems has been investigated by many researchers: as parabolic problems with memory [14,15] and as viscoelastic problems [16,17,18,19,20,21]. For the viscoelastic problem, it has been proved, roughly, that the memory term produces a weak damping which is able by itself to drive the system to rest. It is the kernel which determines the rate of stability. Roughly, the stability rate is similar to the decay rate of the kernel. Various kinds of stability rates (exponential, polynomial and arbitrary) may be found in the literature [12,13,14,15,16,17,18,19,20,21,22]. It is clear now that these models are not only interesting from the mathematical point of view but also from the point of view of applications. They may describe the vibrations of viscoelastic structures in anomalous media.
Moreover, having in mind the telegraph equation [1,2,22]
2 v t 2 + c v t = Δ v ,
Orsingher and Beghin [4] showed that the law of iterated Brownian motion and the telegraph process with Brownian time are expressed by the fractional telegraph equation
C D 2 α v ( t , x ) + a C D α v ( t , x ) = Δ v ( t , x ) , 0 < α < 1 .
Therefore, our present model (1) is a kind of non-linear fractional telegraph problem with memory.
The telegraph equation in (6) is dissipative. This is due to the presence of a first-order derivative which is often called ’frictional’ damping. In fact, it is easy to see that the energy of the system tends to zero in an exponential manner as time tends to infinity. Hence, it is interesting to know whether the lower-order fractional derivative is a dissipative term for the fractional problem in (1). In the affirmative, what would be the rate of convergence?
While we can find a good number of results on the well-posedness of similar problems to ours for orders between 0 and 1 and for problems without the Laplacian memory term, there are very few papers dealing with the exact form in the problem in (1). We identified three papers: one by El-Sayed and Herzallah [23], one by Ponce [24] (for an infinite history) and one by Agarwal et al. [25], dealing with abstract operators. We mention here the result in [25].
Let ( X , . ) be a Banach space and M , N ( t ) t 0 be closed linear operators defined on domains D ( M ) and D ( N ( t ) ) D ( M ) dense in X. We denote by . 1 the graph norm in D ( M ) , R ( ν , M ) : = ( ν I M ) 1 and ρ ( M ) the resolvent of M . The problem
C D β v ( t , x ) = M v ( t , x ) 0 t N ( t s ) v ( s , x ) d s , 1 < β < 2
with initial data ( v ( 0 ) = v 0 , v ( 0 ) = 0 ) in D ( M ) (see [24,25]), admits a classical solution v C ( ( 0 , ) ; D ( M ) ) C 1 ( ( 0 , ) ; X ) such that t 1 β Γ ( 2 β ) ( v v 0 ) C 2 ( ( 0 , ) ; X ) under the following assumptions
(A1) For some 0 < θ 0 π / 2 and every θ < θ 0 , there is a constant C > 0 such that
Σ 0 , β ω : = { ν C : ν 0 , arg ( ν ) < β ω } ρ ( M )
and R ( ν , M ) < C / ν , ν Σ 0 , β ω with ω = θ + π / 2 .
(A2) For each v D ( M ) , N ( . ) v is strongly measurable on ( 0 , ) . There exists a locally integrable function f ( t ) with Laplace transform f ^ ( ν ) , Re ( ν ) > 0 and N ( t ) v f ( t ) v 1 , t > 0 , v D ( M ) . In addition, N ^ : Σ 0 , π / 2 L ( D ( M ) ; X ) (the space of bounded linear operators from D ( M ) into X) has an analytic extension N ˜ to Σ 0 , ω verifying N ˜ ( ν ) v N ˜ ( ν ) v 1 , v D ( M ) and N ˜ ( ν ) = O 1 ν as ν .
(A3) There exists a subspace S dense in ( D ( M ) , . 1 ) and C ˜ > 0 such that M ( S ) D ( M ) , N ^ ( ν ) ( S ) D ( M ) , M N ^ ( ν ) v C ˜ v , v S and ν Σ 0 , ω .
In the present situation our operator is the well-known Laplacian which fortunately satisfies the above conditions in an appropriate space. Therefore, we can profit from this result to ensure the well-posedness of our problem in (1).
In our case, because of the non-linearity present in the equation, we only obtain a local existence result (see [26]). Assuming that the initial data are in D ( A ) (A is the operator in our problem in (1)), there exists a unique local solution in the space
C * : = v C [ 0 , T ) , D ( A ) C 1 [ 0 , T ) , H 0 1 ( Ω ) : t 1 2 α ( w w 0 ) C 2 [ 0 , T ) , L 2 ( Ω ) ,
for T > 0 , where
D ( A ) : = H 0 1 ( Ω ) H 2 ( Ω ) .
For simplicity and convenience, we shall rather consider the problem
C D α C D α v ( t , x ) + a C D α v ( t , x ) + v γ v = Δ v ( t , x ) b 0 t k ( t s ) Δ v ( s , x ) d s , t > 0 , x Ω , 1 / 2 < α < 1 , v ( t , x ) = 0 , t > 0 , x Ω , v ( 0 , x ) = v 0 ( x ) , x Ω .
Because the fractional composition C D α + β v ( t , x ) = C D α C D β v ( t , x ) is not correct in general, the two problems are not exactly the same. The case α = 1 / 2 corresponds to the parabolic case [12,15]. In the case where α = 1 , we recover the viscoelastic problem studied extensively by many authors [7,14,18,19,20,21,22].
We shall shed some light on the effects of the non-linear term, the fractional lower-order term and the memory term on the stability of the system. The lower-order term and the memory term will be considered separately. It is well known, in the integer-order case, that having both at the same time enhances the dissipation. Thus, the challenging question is to consider only one at a time. We have purposefully written the lower-order term as a half-order fractional derivative of the leading fractional derivative as it corresponds to the one in (7). It would be interesting, at least from a mathematical point of view, to treat other orders as well. We will show below that it is a dissipative term. More precisely, we will prove that this lower-order term is by itself (without the memory term) sufficient to guarantee the stability of the system in the Mittag–Leffler manner (solutions of the system will decay as the Mittag–Leffler functions under certain norms). However, without the lower-order fractional term, the Mittag–Leffler stability is only shown for ’small’ kernels. We are unable to decide on the dissipativity of the system in the early stages. This is justified by the fact that the memory term alone is not capable of producing a considerable amount of damping. The product rule (Proposition 5 below) gives rise to a complicated term. It is defined by three singular integrals involving functions of undefined signs. Moreover, the stability is only of local character due to the difficulty to control the non-linear terms. In fact, this local character is not suitable to fraction, it is also the case in integer-order problems. The objective of this paper is to suggest a way to remedy to these difficulties. Unfortunately, we could not compare our results with existing ones as, to the best of our knowledge, this is the first paper treating this as a non-linear issue.
In case both the lower-order fractional term and memory term are present, the system remains Mittag–Leffler stable for reasonable classes of kernels. In fact, the memory term in this situation does not contribute significantly. On the contrary, it only adds some complexity to the estimations. To keep the size of the paper reasonable, we refrain from treating this case. Instead a simple combination of the first two cases is pursued (see Remark 1 below).
The next section contains some preliminary definitions and results. It is followed by Section 3 on the existence and uniqueness issue. In Section 4 we treat the case where only the lower-order term is present. In Section 5 we determine an appropriate function to work with. Section 6 contains the Mittag–Leffler stability in the presence of the memory term. We end the paper by Section 7.

2. Preliminaries

In this section, we give some definitions, formulas and inequalities involving fractional derivatives [27,28]. We also provide the fractional product rule.
Definition 1.
The Riemann–Liouville fractional integral of order α > 0
I α φ ( t ) = 1 Γ ( α ) 0 t ( t s ) α 1 φ ( s ) d s , α > 0
for any measurable function φ provided that the right-hand side exists. Here, Γ ( α ) is the usual Gamma function.
Definition 2
([28]). The fractional derivative of order α in the Caputo is defined by
C D α φ ( t ) = 1 Γ ( 1 α ) 0 t ( t s ) α φ ( s ) d s , 0 < α < 1 ,
C D α φ ( t ) = 1 Γ ( 2 α ) 0 t ( t s ) 1 α φ ( s ) d s , 1 < α < 2 ,
whereas the fractional derivative of order α in the Riemann–Liouville is defined by
R L D α φ ( t ) = 1 Γ ( 1 α ) d d t 0 t ( t s ) α φ ( s ) d s , 0 < α < 1
provided that the integrals exist. Here, the primes denote time differentiation.
The relationship between both derivatives is
R L D α φ ( t ) = φ ( 0 ) t α Γ ( 1 α ) + C D α φ ( t ) , 0 < α < 1 , t > 0 .
The one-parametric and two-parametric Mittag-Leffler functions are defined by [25]
E α ( z ) : = n = 0 z n Γ ( α n + 1 ) , Re ( α ) > 0 ,
and
E α , β ( z ) : = n = 0 z n Γ ( α n + β ) , Re ( α ) > 0 , Re ( β ) > 0 ,
respectively. Notice that E α , 1 ( z ) E α ( z ) .
The power γ is assumed to satisfy
γ > 0 , n = 1 , 2 , 0 < γ < 2 n 2 , n > 2 .
These conditions are imposed so as to ensure an embedding later in our proofs.
Proposition 1
([29]). If φ ( t ) is a differentiable function satisfying
C D α φ ( t ) γ φ ( t ) , 0 < α < 1
for γ > 0 , then φ ( t ) φ ( 0 ) E α ( γ t α ) , t 0 . In case the derivative is of Riemann–Liouville-type, then the decay is of the form t α 1 E α , α ( γ t α ) .
Proposition 2
([27]). For α , β > 0 , we have
λ t α E α , α + β ( λ t α ) = E α , β ( λ t α ) 1 Γ ( β ) .
Proposition 3
([27], p. 61). For μ , α , β > 0 , it holds that
1 Γ ( μ ) 0 t ( t s ) μ 1 E α , β ( λ s α ) s β 1 d s = t μ + β 1 E α , μ + β ( λ t α ) , t > 0 .
Proposition 4
([28], p. 99). If I 1 α φ ( t ) C 1 ( [ 0 , ) ) , 0 < α < 1 and ψ ( t ) is a continuous function, then
R L D α 0 t φ ( t s ) ψ ( s ) d s = 0 t ψ ( s ) R L D α φ ( t s ) d s + ψ ( t ) lim t 0 + I 1 α φ ( t ) , t > 0 .
A fractional version of the well-known product rule in the integer-order case is
Proposition 5
([30]). Let φ ( t ) and ψ ( t ) be absolutely continuous functions. Then, for t 0 and 0 < α < 1 , we have
φ ( t ) C D α ψ ( t ) + ψ ( t ) C D α φ ( t ) = C D α ( φ ψ ) ( t ) + α Γ ( 1 α ) 0 t d ξ ( t ξ ) 1 α 0 ξ φ ( η ) d η ( t η ) α 0 ξ ψ ( s ) d s ( t s ) α .
From now on, we shall drop the letter “C” referring to “Caputo"” in C D α , for notation convenience.
Proposition 6.
Let v ( t ) be a real-valued differentiable function, then for 0 < α < 1
D α v ρ + 2 ρ + 2 v ρ v D α v .
Proof. 
This is an immediate consequence of a combination of a result in [31]
D α v ρ + 2 ρ + 2 v ρ + 1 D α v
and one in [32]
D α v s g n v D α v .

3. Existence and Uniqueness

Clearly, we can rewrite the above equation in (8) as a system of two equations of order α as follows
C D α v ( t , x ) = w ( t , x ) , C D α w ( t , x ) = Δ v ( t , x ) a v ( t , x ) b 0 t k ( t s ) Δ v ( s , x ) d s v γ v
or in matrix form
C D α v ( t , x ) w ( t , x ) = 0 I Δ a I v ( t , x ) w ( t , x ) 0 t 0 0 b k ( t s ) Δ 0 v ( s , x ) w ( s , x ) d s + 0 v γ ( t , x ) v ( t , x ) .
Notice that, if v ( t ) is continuous near zero then w ( 0 ) = C D α v ( 0 ) = 0 .
It is proven in [13] that the abstract problem
C D γ z ( t , x ) = P z ( t , x ) + 0 t Q ( t s ) z ( s , x ) d s + f ( t , z ( t ) ) , 0 < γ < 1 z ( 0 , x ) = z 0 ( x ) X
admits a unique solution under the above assumptions (A1)–(A3) and f ( t , z ( t ) ) is locally Lipschitz with respect to its second argument. For the fractional case, instead of semigroups, we rather use the notion of resolvents defined next (see also [33,34]).
Definition 3.
The family of bounded linear operators  R γ ( t ) t 0  determines a  γ -resolvent for this problem if
(a) The mapping  R γ ( t ) : [ 0 , ) L ( X ) : = L ( X ; X )  is strongly continuous and  R γ ( 0 ) = I .
(b)  z D ( P ) , R γ ( . ) z C [ 0 , ) ; D ( P ) C γ ( 0 , ) ; X  ( C γ ( 0 , ) ; X  designates the space of continuous functions z such that  C D γ z  exists and is continuous), and
C D γ R γ ( t ) z = P R γ ( t ) z + 0 t Q ( t s ) R γ ( s ) z d s = R γ ( t ) P u + 0 t R γ ( t s ) Q ( s ) z d s , t 0 .
The following family
R γ ( t ) : = 1 2 π i Γ μ γ 1 e μ t μ γ I P Q ^ ( μ ) 1 d μ , t 0
is shown to be a γ -resolvent for a suitable path Γ : = { t e i θ : t r } { r e i ζ : θ ζ θ } { t e i θ : t r } oriented counterclockwise where π / 2 < θ < ϕ and r is determined in the proof. An associated family to R γ ( t ) t 0 is defined by
S γ ( t ) : = t 1 γ 2 π i Γ e μ t μ γ I P Q ^ ( μ ) 1 d μ , t 0 .
Theorem 1
([35]). For z 0 D ( P ) , the function
z ( t ) : = R γ ( t ) z 0 + 0 t ( t s ) γ 1 S γ ( t s ) f ( s , z ( s ) ) d s
belongs to C [ 0 , T ) ; D ( P ) C γ ( 0 , T ) ; X and is a solution to the problem above for T > 0 .
As in our case the Laplacian is a sectorial operator (with open sector spectrum and a resolvent which is uniformly bounded from above outside any large sector) and taking into account the assumptions below on the kernel, the conditions (A1)–(A3) are satisfied. Consequently, there exists a unique local solution, to our problem, in the space
C * : = w C [ 0 , T ) , D ( A ) C 1 [ 0 , T ) , H 0 1 ( Ω ) : t α w t C 1 [ 0 , T ) , L 2 ( Ω ) ,
for some T > 0 , where
D ( A ) : = H 0 1 ( Ω ) H 2 ( Ω ) .

4. The Case a 0 and b = 0

Here, we present and prove our stability result. Without loss of generality, we assume a = 1 and consider the functional
E 0 ( t ) : = 1 2 D α v 2 + v 2 + 1 γ + 2 v ( t ) γ + 2 γ + 2 , t 0 .
From Proposition 5, it appears that
D α E 0 ( t ) = D α v 2 + 1 γ + 2 D α v ( t ) γ + 2 γ + 2 Ω v γ v D α v d x α 2 Γ ( 1 α ) Ω 0 t d ξ ( t ξ ) 1 α 0 ξ D α v ( η ) d η ( t η ) α 2 + 0 ξ v ( η ) d η ( t η ) α 2 d x , t 0 .
We shall assume that our initial data satisfy ( v 0 , v 1 ) H 0 1 ( Ω ) 2 and start the process of looking for a suitable functional to work with. In the next lemma, we introduce the first functional.
Lemma 1.
For the functional
Ψ 0 ( t ) : = Ω v D α v d x , t 0 ,
it holds that
D α Ψ 0 ( t ) 1 + C p 2 D α v 2 1 2 v 2 v ( t ) γ + 2 γ + 2 + α 2 Γ ( 1 α ) Ω 0 t d ξ ( t ξ ) 1 α C p 0 ξ v ( s ) d s ( t s ) α 2 + 0 ξ D α v ( s ) d s ( t s ) α 2 d x , t 0
where C p is the Poincaré constant.
Proof. 
As D α v is absolutely continuous, we may apply Proposition 5 to find, for t 0
D α Ψ 0 ( t ) = D α v 2 + Ω v D α D α v d x α Γ ( 1 α ) Ω 0 t d ξ ( t ξ ) 1 α 0 ξ v ( η ) d η ( t η ) α 0 ξ D α v ( s ) d s ( t s ) α d x = D α v 2 v 2 Ω v D α v d x v ( t ) γ + 2 γ + 2 α Γ ( 1 α ) Ω 0 t d ξ ( t ξ ) 1 α 0 ξ v ( η ) d η ( t η ) α 0 ξ D α v ( s ) d s ( t s ) α d x .
This expression is evaluated as follows
D α Ψ 0 ( t ) 1 + C p 2 D α v 2 1 2 v 2 v ( t ) γ + 2 γ + 2 + α 2 Γ ( 1 α ) Ω 0 t d ξ ( t ξ ) 1 α C p 0 ξ v ( s ) d s ( t s ) α 2 + 0 ξ D α v ( s ) d s ( t s ) α 2 d x , t 0 .
The proof is complete. □
Theorem 2.
The problem in (8) with a = 1 and b = 0 , is Mittag–Leffler stable. This means that there exist positive constants B 0 (depending on E 0 ( 0 ) ) and ν 0 such that
E 0 ( t ) B 0 E α ( ν 0 t α ) , t 0 .
Proof. 
Let
L 0 ( t ) : = N E 0 ( t ) + M Ψ 0 ( t ) , t 0
for some positive constants M and N, to be determined. Then, in view of Proposition 6, we may write
D α L 0 ( t ) N D α v 2 + M 1 + C p 2 D α v 2 M 2 v 2 M v ( t ) γ + 2 γ + 2 α N 2 Γ ( 1 α ) Ω 0 t d ξ ( t ξ ) 1 α 0 ξ D α v ( η ) d η ( t η ) α 2 + 0 ξ v ( η ) d η ( t η ) α 2 d x + α M 2 Γ ( 1 α ) 0 t d ξ ( t ξ ) 1 α C p 0 ξ v ( s ) d s ( t s ) α 2 + 0 ξ D α v ( s ) d s ( t s ) α 2 d x , t 0
or
D α L 0 ( t ) M 1 + C p 2 N D α v 2 M 2 v 2 M v ( t ) γ + 2 γ + 2 + α 2 Γ ( 1 α ) Ω 0 t d ξ ( t ξ ) 1 α × C p M N 0 ξ v ( s ) d s ( t s ) α 2 + M N 0 ξ D α v ( η ) d η ( t η ) α 2 d x
for t 0 . To obtain C 1 E 0 ( t ) for some positive constant C 1 in the right hand side of (10), we need
M 1 + C p 2 < N , C p M N , M N .
This requires
N > max C p , 1 + C p 2 M .
Consequently
D α L 0 ( t ) C 1 E 0 ( t ) , t 0
for some positive constant C 1 and as L 0 ( t ) and E 0 ( t ) are equivalent, we obtain
D α L 0 ( t ) C 2 L 0 ( t ) , t 0
for some positive constant C 2 . The conclusion follows from Proposition 1 and the equivalence of L 0 ( t ) and E 0 ( t ) again. □

5. A Useful Functional

In this section, we determine an appropriate functional before we tackle the case of the memory term. We shall use the notation
( k v ) ( t ) : = Ω 0 t k ( t s ) v ( t ) v ( s ) 2 d s d x , t 0 .
Proposition 7.
If I 1 α k ( t ) C 1 ( [ 0 , ) ) , 0 < α < 1 , then the following identity holds
Ω D α v . 0 t k ( t s ) v ( s ) d s d x = 1 2 ( R L D α k v ) + 0 t k ( s ) d s D α v 2 α 2 Γ ( 1 α ) 0 t d ξ ( t ξ ) 1 α 0 ξ k ( η ) d η ( t η ) α 0 ξ v ( s ) 2 d s ( t s ) α 1 2 D α ( k v ) ( t ) + α Γ ( 1 α ) × Ω 0 t d ξ ( t ξ ) 1 α 0 ξ v ( η ) d η ( t η ) α . 0 ξ ( t s ) α 0 s k ( s τ ) v ( τ ) d τ d s d x , t 0 .
Proof. 
By virtue of Proposition 5 and the identity
Ω 0 t k ( t s ) v ( t ) v ( s ) 2 d s d x = v 2 0 t k ( t s ) d s + 0 t k ( t s ) v ( s ) 2 d s 2 Ω v . 0 t k ( t s ) v ( s ) d s d x , t 0
we see that
D α ( k v ) ( t ) = D α 0 t k ( t s ) d s v 2 + 0 t k ( t s ) d s D α v 2 α Γ ( 1 α ) 0 t d ξ ( t ξ ) 1 α 0 ξ k ( η ) d η ( t η ) α 0 ξ v ( s ) 2 d s ( t s ) α + D α 0 t k ( t s ) v ( s ) 2 d s 2 Ω D α v . 0 t k ( t s ) v ( s ) d s d x 2 Ω v . D α 0 t k ( t s ) v ( s ) d s d x + 2 α Γ ( 1 α ) × Ω 0 t d ξ ( t ξ ) 1 α 0 ξ v ( η ) d η ( t η ) α . 0 ξ ( t s ) α 0 s k ( s τ ) v ( τ ) d τ d s d x , t 0 .
Notice that
( R L D α k v ) ( t ) = v 2 0 t R L D α k ( t s ) d s + 0 t R L D α k ( t s ) v ( s ) 2 d s 2 Ω v . 0 t R L D α k ( t s ) v ( s ) d s d x , t 0
and by Proposition 4
D α 0 t k ( t s ) d s = 0 t R L D α k ( t s ) d s + lim t 0 + I 1 α k ( t ) , t 0 .
Again, by Proposition 4 and the summability of k, we deduce that
D α 0 t k ( t s ) v ( s ) 2 d s = R L D α 0 t k ( t s ) v ( s ) 2 d s t α Γ ( 1 α ) 0 t k ( t s ) v ( s ) 2 d s t = 0 = 0 t R L D α k ( t s ) v ( s ) 2 d s + v ( t ) 2 lim t 0 + I 1 α k ( t ) , t > 0
and
Ω v . D α 0 t k ( t s ) v ( s ) d s d x = Ω v . 0 t R L D α k ( t s ) v ( s ) d s d x + v ( t ) 2 lim t 0 + I 1 α k ( t ) , t > 0 .
Gathering the previous relations (13)–(15) in (11), we find
D α ( k v ) ( t ) = v 2 0 t R L D α k ( t s ) d s + 0 t k ( t s ) d s D α v 2 + 0 t R L D α k ( t s ) v ( s ) 2 d s α Γ ( 1 α ) 0 t d ξ ( t ξ ) 1 α 0 ξ k ( η ) d η ( t η ) α 0 ξ v ( s ) 2 d s ( t s ) α 2 Ω D α v . 0 t k ( t s ) v ( s ) d s d x 2 Ω v . 0 t R L D α k ( t s ) v ( s ) d s d x + 2 α Γ ( 1 α ) × Ω 0 t d ξ ( t ξ ) 1 α 0 ξ v ( η ) d η ( t η ) α . 0 ξ ( t s ) α 0 s k ( s τ ) v ( τ ) d τ d s d x , t 0 .
Therefore, in view of the relation (12) and (16), we obtain for t 0
D α ( k v ) ( t ) = ( R L D α k v ) + 0 t k ( s ) d s D α v 2 α Γ ( 1 α ) 0 t d ξ ( t ξ ) 1 α 0 ξ k ( η ) d η ( t η ) α 0 ξ v ( s ) 2 d s ( t s ) α 2 Ω D α v . 0 t k ( t s ) v ( s ) d s d x + 2 α Γ ( 1 α ) Ω 0 t d ξ ( t ξ ) 1 α 0 ξ v ( η ) d η ( t η ) α . 0 ξ ( t s ) α 0 s k ( s τ ) v ( τ ) d τ d s d x .
This finishes the proof. □
A multiplication of the equation in (8) (with a = 0 and b = 1 ) by D α v gives
Ω D α v D α D α v d x + Ω v γ v D α v d x = Ω D α v . v d x + Ω D α v . 0 t k ( t s ) v ( s , x ) d s d x , t 0
and by the above Proposition 6, we see that
Ω D α v D α D α v d x + Ω v γ v D α v d x + Ω D α v . v d x 1 2 D α v 2 0 t k ( s ) d s + 1 2 D α ( k v ) ( t ) = 1 2 ( R L D α k v ) α 2 Γ ( 1 α ) 0 t d ξ ( t ξ ) 1 α 0 ξ k ( η ) d η ( t η ) α 0 ξ v ( s ) 2 d s ( t s ) α + α Γ ( 1 α ) Ω 0 t d ξ ( t ξ ) 1 α × 0 ξ v ( η ) d η ( t η ) α . 0 ξ ( t s ) α 0 s k ( s τ ) v ( τ ) d τ d s d x , t 0 .
This identity is the main reason behind choosing the functional
E 1 ( t ) : = 1 2 D α v 2 + 1 0 t k ( s ) d s v 2 + k v ( t ) + 1 γ + 2 v ( t ) γ + 2 γ + 2 .

6. The Case a = 0 and b 0

Here, in the absence of the fractional damping, there remains only the viscoelastic term. Without loss of generality, we pick b = 1 . Our assumptions on the kernel k are:
(K) k L l o c 1 ( R + ) is a non-negative continuous function satisfying k ^ ( λ ) = O λ 1 and I 1 α k ( t ) C 1 ( [ 0 , ) ) ,
R L D α k ( t ) C k k ( t ) , t > 0
and
k ¯ : = 0 k ( s ) d s < 1
for some positive constant C k .
To deal with the non-linearity we need to assume that
(Γ) 0 < γ < n n 2 1 = 2 n 2 if n 3 and γ > 0 if n = 1 , 2 .
This assumption will allow us to use the embedding H 0 1 ( Ω ) L q ( Ω ) when 2 q 2 n / ( n 2 ) if n 3 and q 2 if n = 1 , 2 .
It is easy to see from Proposition 1 that k ( t ) K t α 1 E α , α ( C k t α ) , t 0 for some K > 0 and from Propositions 2 and 3 that such kernels are summable and
0 t k ( s ) d s K 0 t s α 1 E α , α ( C k s α ) d s = K t α E α , α + 1 ( C k t α ) K / C k , t > 0 .
In fact, this family of functions t α 1 E α , α ( C k t α ) satisfies (K). Moreover, our proofs below are valid for initial data in H 0 1 ( Ω ) 2 .
Lemma 2.
The functional  E 1 ( t )  satisfies along solutions of (8), for  t 0
D α E 1 ( t ) 1 2 ( R L D α k v ) 1 2 I 1 α k ( t ) v 2 α 2 Γ ( 1 α ) Ω 0 t d ξ ( t ξ ) 1 α 0 ξ D α v ( η ) d η ( t η ) α 2 d x α 4 Γ ( 1 α ) Ω 0 t d ξ ( t ξ ) 1 α 0 ξ v ( η ) d η ( t η ) α 2 d x + α Γ ( 1 α ) Ω 0 t d ξ ( t ξ ) 1 α 0 ξ ( t s ) α 0 s k ( s τ ) v ( τ ) d τ d s 2 d x .
Proof. 
Clearly, a direct fractional differentiation of E 1 ( t ) , using Proposition 5, yields for t > 0
D α E 1 ( t ) = 1 2 D α D α v 2 + 1 0 t k ( s ) d s v 2 + k v ( t ) + 1 γ + 2 D α v ( t ) γ + 2 γ + 2 = Ω D α v D α D α v d x α 2 Γ ( 1 α ) Ω 0 t d ξ ( t ξ ) 1 α 0 ξ D α v ( η ) d η ( t η ) α 2 d x + 1 0 t k ( s ) d s × Ω v . D α v d x α 2 Γ ( 1 α ) Ω 0 t d ξ ( t ξ ) 1 α 0 ξ v ( η ) d η ( t η ) α 2 d x 1 2 D α 0 t k ( s ) d s v 2 + α 2 Γ ( 1 α ) 0 t d ξ ( t ξ ) 1 α 0 ξ v ( s ) 2 d s ( t s ) α 0 ξ k ( η ) d η ( t η ) α + 1 2 D α k v ( t ) + 1 γ + 2 D α v ( t ) γ + 2 γ + 2 .
Along solutions of (8) (see (17)) and in view of Proposition 4, we find
D α E 1 ( t ) = 1 2 ( R L D α k v ) 1 2 I 1 α k ( t ) v 2 + 1 γ + 2 D α v ( t ) γ + 2 γ + 2 Ω v γ v D α v d x α 2 Γ ( 1 α ) Ω 0 t d ξ ( t ξ ) 1 α 0 ξ D α v ( η ) d η ( t η ) α 2 d x α 2 Γ ( 1 α ) Ω 0 t d ξ ( t ξ ) 1 α 0 ξ v ( η ) d η ( t η ) α 2 d x + α Γ ( 1 α ) Ω 0 t d ξ ( t ξ ) 1 α × 0 ξ v ( η ) d η ( t η ) α . 0 ξ ( t s ) α 0 s k ( s τ ) v ( τ ) d τ d s d x , t 0 .
We conclude by using Proposition 6 and the evaluation
Ω 0 t d ξ ( t ξ ) 1 α 0 ξ v ( η ) d η ( t η ) α . 0 ξ ( t s ) α 0 s k ( s τ ) v ( τ ) d τ d s d x 1 4 Ω 0 t d ξ ( t ξ ) 1 α 0 ξ v ( η ) d η ( t η ) α 2 d x + Ω 0 t d ξ ( t ξ ) 1 α 0 ξ ( t s ) α 0 s k ( s τ ) v ( τ ) d τ d s 2 d x .
This finishes the proof. □
The last term in the estimation of D α E 1 ( t ) needs to be controlled somehow. This term does not appear in the integer-order case (second-order) and, whatever evaluation we derive, its weight is going to be considerable. We suggest a new functional in the next lemma.
Lemma 3.
The fractional derivative of the functional
Ψ 1 ( t ) : = Ω 0 t k ( t s ) v ( s ) d s 2 d x , t 0
satisfies
D α Ψ 1 ( t ) 2 λ 1 k ¯ Ω 0 t k ( t s ) v ( s ) 2 d s d x + K λ 1 ( R L D α k v ) ( t ) + K E α ( C k t α ) 2 λ 1 v 2 α Γ ( 1 α ) Ω 0 t d ξ ( t ξ ) 1 α 0 ξ ( t η ) α 0 η k ( η τ ) v ( τ ) d τ d η 2 d x ,
for  λ 1 > 0 , t 0 .
Proof. 
With the help of Proposition 5, we can write
D α Ψ 1 ( t ) = 2 Ω 0 t k ( t s ) v ( s ) d s . D α 0 t k ( t s ) v ( s ) d s d x α Γ ( 1 α ) Ω 0 t d ξ ( t ξ ) 1 α 0 ξ ( t η ) α 0 η k ( η τ ) v ( τ ) d τ d η 2 d x , t 0
and the relation (9) gives
D α Ψ 1 ( t ) = 2 Ω 0 t k ( t s ) v ( s ) d s × . R L D α 0 t k ( t s ) v ( s ) d s d x t α Γ ( 1 α ) 0 t k ( t s ) v ( s ) d s t = 0 d x α Γ ( 1 α ) Ω 0 t d ξ ( t ξ ) 1 α 0 ξ ( t η ) α 0 η k ( η τ ) v ( τ ) d τ d η 2 d x , t 0 .
Next, the summability of k and Proposition 4 imply
D α Ψ 1 ( t ) = 2 Ω 0 t k ( t s ) v ( s ) d s . 0 t R L D α k ( t s ) v ( s ) v ( t ) d s + v ( t ) 0 t R L D α k ( s ) d s + v ( t ) lim t 0 + I 1 α k ( t ) d x α Γ ( 1 α ) Ω 0 t d ξ ( t ξ ) 1 α 0 ξ ( t η ) α 0 η k ( η τ ) v ( τ ) d τ d η 2 d x = 2 Ω 0 t k ( t s ) v ( s ) d s . 0 t R L D α k ( t s ) v ( s ) v ( t ) d s d x + 2 I 1 α k ( t ) Ω 0 t k ( t s ) v ( s ) d s . v ( t ) d x α Γ ( 1 α ) Ω 0 t d ξ ( t ξ ) 1 α 0 ξ ( t η ) α 0 η k ( η τ ) v ( τ ) d τ d η 2 d x , t 0
or, for λ 1 > 0
D α Ψ 1 ( t ) 2 λ 1 k ¯ Ω 0 t k ( t s ) v ( s ) 2 d s d x + 1 λ 1 0 t R L D α k ( s ) d s ( R L D α k v ) ( t ) + I 1 α k ( t ) 2 λ 1 v 2 α Γ ( 1 α ) Ω 0 t d ξ ( t ξ ) 1 α 0 ξ ( t η ) α 0 η k ( η τ ) v ( τ ) d τ d η 2 d x , t 0 .
Observing that (see Proposition 3)
I 1 α k ( t ) K Γ ( 1 α ) 0 t ( t s ) α 1 E α , α ( C k ( t s ) α ) s α d s = K E α ( C k t α )
and
0 t R L D α k ( s ) d s = 0 t D I 1 α k ( s ) d s = I 1 α k ( s ) 0 t = I 1 α k ( 0 ) I 1 α k ( t ) K
it appears that
D α Ψ 1 ( t ) 2 λ 1 k ¯ Ω 0 t k ( t s ) v ( s ) 2 d s d x + K λ 1 ( R L D α k v ) ( t ) + K E α ( C k t α ) 2 λ 1 v 2 α Γ ( 1 α ) Ω 0 t d ξ ( t ξ ) 1 α 0 ξ ( t η ) α 0 η k ( η τ ) v ( τ ) d τ d η 2 d x , t 0
for λ 1 > 0 .
The second functional we introduce is in the next lemma.
Lemma 4.
For the functional
Ψ 2 ( t ) : = Ω v D α v d x , t 0 ,
it holds that
D α Ψ 2 ( t ) D α v 2 1 λ 2 v 2 v ( t ) γ + 2 γ + 2 + k ¯ 4 λ 2 Ω 0 t k ( t s ) v ( s ) 2 d s d x + α λ 3 C p Γ ( 1 α ) Ω 0 t d ξ ( t ξ ) 1 α 0 ξ v ( s ) d s ( t s ) α 2 d x + α 4 λ 3 Γ ( 1 α ) Ω 0 t d ξ ( t ξ ) 1 α 0 ξ D α v ( s ) d s ( t s ) α 2 d x , t 0
for λ 2 , λ 3 > 0 , where C p is the Poincaré constant.
Proof. 
By Proposition 5 we have
D α Ψ 2 ( t ) = D α v 2 + Ω v D α D α v d x α Γ ( 1 α ) Ω 0 t d ξ ( t ξ ) 1 α 0 ξ v ( η ) d η ( t η ) α 0 ξ D α v ( s ) d s ( t s ) α d x = D α v 2 v 2 + Ω v . 0 t k ( t s ) v ( s ) d s d x v ( t ) γ + 2 γ + 2 α Γ ( 1 α ) Ω 0 t d ξ ( t ξ ) 1 α 0 ξ v ( η ) d η ( t η ) α 0 ξ D α v ( s ) d s ( t s ) α d x , t 0 .
Therefore, for λ 2 , λ 3 > 0 and t 0
D α Ψ 2 ( t ) D α v 2 1 λ 2 v 2 + k ¯ 4 λ 2 Ω 0 t k ( t s ) v ( s ) 2 d s d x v ( t ) γ + 2 γ + 2 + α λ 3 C p Γ ( 1 α ) Ω 0 t d ξ ( t ξ ) 1 α 0 ξ v ( s ) d s ( t s ) α 2 d x + α 4 λ 3 Γ ( 1 α ) Ω 0 t d ξ ( t ξ ) 1 α 0 ξ D α v ( s ) d s ( t s ) α 2 d x .
The proof is complete. □
Lemma 5.
The functional
Ψ 3 ( t ) : = 0 t k ( t s ) v ( s ) 2 d s , t 0
satisfies
D α Ψ 3 ( t ) C k Ψ 3 ( t ) + lim t 0 + I 1 α k ( t ) v 2 , t 0 .
Proof. 
This follows immediately from Proposition 4 and (K). □
As we have removed the lower-order fractional term, which represented the damping and was responsible for the appearance of the nice term D α v 2 , we need to find a way to come up with this missing and essential term. The next suggested functional will take care of this matter.
Lemma 6.
If the exponent γ satisfies (Γ), then for the functional
Ψ 4 ( t ) : = Ω D α v 0 t k ( t s ) v ( t ) v ( s ) d s d x , t 0 ,
we have, along solutions of problem (8)
D α Ψ 4 ( t ) = λ 5 v 2 + k ¯ 1 4 λ 5 + C p 4 λ 6 + 1 ( k v ) ( t ) + λ 4 k 0 2 D α v 2 + λ 6 C 1 E 1 γ + 1 ( t ) K C p 2 k 0 ( R L D α k v ) ( t ) + C 4 E α ( C k t α ) λ 4 Ω 0 t d ξ ( t ξ ) 1 α 0 ξ v ( s ) d s ( t s ) α 2 d x + α λ 7 Γ ( 1 α ) Ω 0 t d ξ ( t ξ ) 1 α 0 ξ D α v ( η ) d η ( t η ) α 2 d x + α 4 λ 7 Γ ( 1 α ) Ω 0 t d ξ ( t ξ ) 1 α 0 ξ ( t s ) α 0 s k ( s τ ) v ( s ) v ( τ ) d τ d s 2 d x ,
for λ i > 0 , i = 4 , 5 , 6 , 7 , t t 0 > 0 where k 0 = 0 t 0 k ( s ) d s for some t 0 > 0 .
Proof. 
By Proposition 5, we see that
D α Ψ 4 ( t ) = Ω D α D α v 0 t k ( t s ) v ( t ) v ( s ) d s d x Ω D α v D α 0 t k ( t s ) v ( t ) v ( s ) d s d x + α Γ ( 1 α ) Ω 0 t d ξ ( t ξ ) 1 α 0 ξ D α v ( η ) d η ( t η ) α × 0 ξ ( t s ) α 0 s k ( s τ ) v ( s ) v ( τ ) d τ d s d x , t 0
and along solutions of (8), we find for t 0
D α Ψ 4 ( t ) = Ω Δ v 0 t k ( t s ) Δ v ( s ) d s v γ v 0 t k ( t s ) v ( t ) v ( s ) d s d x Ω D α v D α v ( t ) 0 t k ( t s ) d s D α 0 t k ( t s ) v ( s ) d s d x + α Γ ( 1 α ) × Ω 0 t d ξ ( t ξ ) 1 α 0 ξ D α v ( η ) d η ( t η ) α 0 ξ ( t s ) α 0 s k ( s τ ) v ( s ) v ( τ ) d τ d s d x .
Next, we apply Proposition 4
D α Ψ 4 ( t ) = Ω 1 0 t k ( s ) d s v + 0 t k ( t s ) v ( t ) v ( s ) d s . 0 t k ( t s ) v ( t ) v ( s ) d s d x Ω D α v D α v ( t ) 0 t k ( t s ) d s 0 t R L D α k ( t s ) v ( s ) d s v ( t ) I 1 α k ( 0 ) d x + Ω v γ v 0 t k ( t s ) v ( t ) v ( s ) d s d x + α Γ ( 1 α ) Ω 0 t d ξ ( t ξ ) 1 α × 0 ξ D α v ( η ) d η ( t η ) α 0 ξ ( t s ) α 0 s k ( s τ ) v ( s ) v ( τ ) d τ d s d x
or, for t 0
D α Ψ 4 ( t ) = 1 0 t k ( s ) d s Ω v . 0 t k ( t s ) v ( t ) v ( s ) d s d x + Ω 0 t k ( t s ) v ( t ) v ( s ) d s 2 d x + Ω v γ v 0 t k ( t s ) v ( t ) v ( s ) d s d x Ω D α v D α v ( t ) 0 t k ( t s ) d s 0 t R L D α k ( t s ) v ( s ) d s d x + I 1 α k ( 0 ) Ω v D α v d x + α Γ ( 1 α ) Ω 0 t d ξ ( t ξ ) 1 α 0 ξ D α v ( η ) d η ( t η ) α 0 ξ ( t s ) α 0 s k ( s τ ) v ( s ) v ( τ ) d τ d s d x .
Another application of Proposition 5 gives
D α v ( t ) 0 t k ( t s ) d s = 0 t k ( s ) d s D α v ( t ) + v ( t ) 0 t R L D α k ( t s ) d s + I 1 α k ( 0 ) α Γ ( 1 α ) 0 t d ξ ( t ξ ) 1 α 0 ξ k ( η ) d η ( t η ) α 0 ξ v ( s ) d s ( t s ) α , t 0 .
Therefore, for λ 4 > 0 , t 0
Ω D α v D α v ( t ) 0 t k ( s ) d s d x = Ω v D α v 0 t R L D α k ( t s ) d s + I 1 α k ( 0 ) d x 0 t k ( s ) d s D α v 2 + α Γ ( 1 α ) Ω D α v 0 t d ξ ( t ξ ) 1 α 0 ξ k ( η ) d η ( t η ) α 0 ξ v ( s ) d s ( t s ) α d x Ω v D α v 0 t R L D α k ( t s ) d s d x I 1 α k ( 0 ) Ω v D α v d x 0 t k ( s ) d s D α v 2 + λ 4 D α v 2 + α 2 4 λ 4 Γ 2 ( 1 α ) Ω 0 t d ξ ( t ξ ) 1 α 0 ξ k ( η ) d η ( t η ) α 0 ξ v ( s ) d s ( t s ) α 2 d x .
We can apply the Cauchy-Schwarz inequality to the last term in (19) to get
Ω 0 t d ξ ( t ξ ) 1 α 0 ξ k ( η ) d η ( t η ) α 0 ξ v ( s ) d s ( t s ) α 2 d x Ω 0 t d ξ ( t ξ ) 1 α 0 ξ k ( η ) d η ( t η ) α 2 0 t d ξ ( t ξ ) 1 α 0 ξ v ( s ) d s ( t s ) α 2 d x 0 t d ξ ( t ξ ) 1 α 0 ξ k ( η ) d η ( t η ) α 2 Ω 0 t d ξ ( t ξ ) 1 α 0 ξ v ( s ) d s ( t s ) α 2 d x , t 0 .
In view of assumption (K), Proposition 3 and the fact that E α ( C k t α ) C 2 t α (for some positive constant C 2 ) away from zero
0 t d ξ ( t ξ ) 1 α 0 ξ k ( η ) d η ( t η ) α 2 K Γ ( 1 α ) E α ( C k t α ) 2 0 t d ξ ( t ξ ) 1 α C 3 E α ( C k t α ) ,
for t t 0 > 0 where C 3 : = K 2 Γ 2 ( 1 α ) C 2 . Therefore
Ω D α v D α v ( t ) 0 t k ( s ) d s d x 0 t k ( s ) d s D α v 2 Ω v D α v 0 t R L D α k ( t s ) d s d x + λ 4 D α v 2 I 1 α k ( 0 ) Ω v D α v d x + C 4 E α ( C k t α ) λ 4 Ω 0 t d ξ ( t ξ ) 1 α 0 ξ v ( s ) d s ( t s ) α 2 d x ,
for λ 4 > 0 , t t 0 > 0 where C 4 = α 2 C 3 / 4 Γ 2 ( 1 α ) . Moreover, it is clear that
Ω v . 0 t k ( t s ) v ( t ) v ( s ) d s d x λ 5 v 2 + k ¯ 4 λ 5 ( k v ) ( t ) , λ 5 > 0 , t 0
and the embedding mentioned above ( H 0 1 ( Ω ) L q ( Ω ) when 2 q 2 n / ( n 2 ) if n 3 and q 2 if n = 1 , 2 ) gives
Ω v γ v 0 t k ( t s ) v ( t ) v ( s ) d s d x λ 6 Ω v 2 ( γ + 1 ) d x + k ¯ C p 4 λ 6 ( k v ) ( t ) λ 6 C 1 E 1 γ + 1 ( t ) + k ¯ C p 4 λ 6 ( k v ) ( t ) , λ 6 > 0
where C 1 is the embedding constant. Observe also that
0 0 t R L D α k ( s ) d s = 0 t D I 1 α k ( s ) d s = I 1 α k ( s ) 0 t = I 1 α k ( 0 ) I 1 α k ( t )
implies
Ω D α v 0 t R L D α k ( t s ) v ( s ) d s d x Ω v D α v 0 t R L D α k ( t s ) d s d x = Ω D α v 0 t R L D α k ( t s ) v ( t ) v ( s ) d s d x k 0 2 D α v 2 + C p 2 k 0 0 t R L D α k ( s ) d s ( R L D α k v ) ( t ) k 0 2 D α v 2 C p 2 k 0 I 1 α k ( 0 ) ( R L D α k v ) ( t ) k 0 2 D α v 2 K C p 2 k 0 ( R L D α k v ) ( t ) , t 0 .
Finally, for λ 7 > 0
Ω 0 t d ξ ( t ξ ) 1 α 0 ξ D α v ( η ) d η ( t η ) α 0 ξ ( t s ) α 0 s k ( s τ ) v ( s ) v ( τ ) d τ d s d x λ 7 Ω 0 t d ξ ( t ξ ) 1 α 0 ξ D α v ( η ) d η ( t η ) α 2 d x + 1 4 λ 7 Ω 0 t d ξ ( t ξ ) 1 α 0 ξ ( t s ) α 0 s k ( s τ ) v ( s ) v ( τ ) d τ d s 2 d x .
Taking into account the relations (20)–(24) in (18), we obtain
D α Ψ 4 ( t ) λ 5 v 2 + k ¯ 1 4 λ 5 + C p 4 λ 6 + 1 ( k v ) ( t ) + λ 4 k 0 2 D α v 2 K C p 2 k 0 ( R L D α k v ) ( t ) + λ 6 C 1 E 1 γ + 1 ( t ) + C 4 E α ( C k t α ) λ 4 Ω 0 t d ξ ( t ξ ) 1 α 0 ξ v ( s ) d s ( t s ) α 2 d x + α λ 7 Γ ( 1 α ) Ω 0 t d ξ ( t ξ ) 1 α × 0 ξ D α v ( η ) d η ( t η ) α 2 d x + α 4 λ 7 Γ ( 1 α ) Ω 0 t d ξ ( t ξ ) 1 α × 0 ξ ( t s ) α 0 s k ( s τ ) v ( s ) v ( τ ) d τ d s 2 d x , t t 0 > 0 .
The proof is complete. □
One more term needs to be controlled. To this end we introduce the functional
Ψ 5 ( t ) : = Ω 0 t k ( t s ) v ( t ) v ( s ) d s 2 d x , t 0 .
Lemma 7.
The above functional  Ψ 5 ( t )  fulfills
D α Ψ 5 ( t ) λ 8 D α v 2 + C p λ 9 k ¯ + k ¯ 3 λ 8 + λ 10 α k ¯ Γ ( 1 α ) k v ( t ) + C p K λ 9 R L D α k v ( t ) + C 5 E α ( C k t α ) λ 10 Ω 0 t d ξ ( t ξ ) 1 α 0 ξ v ( s ) d s ( t s ) α 2 d x α Γ ( 1 α ) Ω 0 t d ξ ( t ξ ) 1 α 0 ξ ( t s ) α 0 s k ( s τ ) v ( s ) v ( τ ) d τ d s 2 d x , t 0 .
Proof. 
Clearly,
D α Ψ 5 ( t ) = 2 Ω 0 t k ( t s ) v ( t ) v ( s ) d s D α 0 t k ( t s ) v ( t ) v ( s ) d s d x α Γ ( 1 α ) Ω 0 t d ξ ( t ξ ) 1 α 0 ξ ( t s ) α 0 s k ( s τ ) v ( s ) v ( τ ) d τ d s 2 d x
and
Ω D α 0 t k ( t s ) v ( t ) v ( s ) d s d x = Ω D α v ( t ) 0 t k ( t s ) d s d x Ω D α 0 t k ( t s ) v ( s ) d s d x = 0 t k ( s ) d s D α v ( t ) + v ( t ) 0 t R L D α k ( t s ) d s + I 1 α k ( 0 ) α Γ ( 1 α ) 0 t d ξ ( t ξ ) 1 α 0 ξ k ( η ) d η ( t η ) α 0 ξ v ( s ) d s ( t s ) α Ω D α 0 t k ( t s ) v ( s ) d s d x .
Therefore,
D α Ψ 5 ( t ) = 2 0 t k ( s ) d s Ω D α v ( t ) 0 t k ( t s ) v ( t ) v ( s ) d s d x + 2 Ω 0 t k ( t s ) v ( t ) v ( s ) d s v ( t ) 0 t R L D α k ( t s ) d s + I 1 α k ( 0 ) d x 2 α Γ ( 1 α ) Ω 0 t k ( t s ) v ( t ) v ( s ) d s 0 t d ξ ( t ξ ) 1 α 0 ξ k ( η ) d η ( t η ) α 0 ξ v ( s ) d s ( t s ) α d x 2 Ω 0 t k ( t s ) v ( t ) v ( s ) d s D α 0 t k ( t s ) v ( s ) d s d x α Γ ( 1 α ) Ω 0 t d ξ ( t ξ ) 1 α 0 ξ ( t s ) α 0 s k ( s τ ) v ( s ) v ( τ ) d τ d s 2 d x , t 0
or, using (14), we find
D α Ψ 5 ( t ) = 2 0 t k ( s ) d s Ω D α v ( t ) 0 t k ( t s ) v ( t ) v ( s ) d s d x + 2 Ω 0 t k ( t s ) v ( t ) v ( s ) d s v ( t ) 0 t R L D α k ( t s ) d s + I 1 α k ( 0 ) d x 2 α Γ ( 1 α ) Ω 0 t k ( t s ) v ( t ) v ( s ) d s 0 t d ξ ( t ξ ) 1 α 0 ξ k ( η ) d η ( t η ) α 0 ξ v ( s ) d s ( t s ) α d x 2 Ω 0 t k ( t s ) v ( t ) v ( s ) d s R L D α 0 t k ( t s ) v ( s ) d s + v ( t ) I 1 α k ( 0 ) d x α Γ ( 1 α ) Ω 0 t d ξ ( t ξ ) 1 α 0 ξ ( t s ) α 0 t k ( t s ) v ( t ) v ( s ) d s d s 2 d x , t 0
that is
D α Ψ 5 ( t ) λ 8 D α v 2 + k ¯ 3 C p λ 8 k v ( t ) + 2 Ω 0 t k ( t s ) v ( t ) v ( s ) d s 0 t R L D α k ( t s ) v ( t ) v ( s ) d s d x 2 α Γ ( 1 α ) Ω 0 t k ( t s ) v ( t ) v ( s ) d s 0 t d ξ ( t ξ ) 1 α 0 ξ k ( η ) d η ( t η ) α 0 ξ v ( s ) d s ( t s ) α d x α Γ ( 1 α ) Ω 0 t d ξ ( t ξ ) 1 α 0 ξ ( t s ) α 0 s k ( s τ ) v ( s ) v ( τ ) d τ d s 2 d x , t 0 .
Furthermore, we may evaluate
2 Ω 0 t k ( t s ) v ( t ) v ( s ) d s 0 t R L D α k ( t s ) v ( t ) v ( s ) d s d x λ 9 k ¯ Ω 0 t k ( t s ) v ( t ) v ( s ) 2 d s d x + 1 λ 9 0 t R L D α k ( s ) d s Ω 0 t R L D α k ( t s ) v ( t ) v ( s ) 2 d s d x λ 9 k ¯ C p k v ( t ) + C p λ 9 0 t R L D α k ( s ) d s R L D α k v ( t ) , λ 9 > 0
and as
0 0 t R L D α k ( s ) d s = 0 t D I 1 α k ( s ) d s = I 1 α k ( 0 ) I 1 α k ( t ) K
we deduce that
2 Ω 0 t k ( t s ) v ( t ) v ( s ) d s 0 t R L D α k ( t s ) v ( t ) v ( s ) d s d x λ 9 k ¯ C p k v ( t ) + C p K λ 9 R L D α k v ( t ) .
Therefore,
D α Ψ 5 ( t ) λ 8 D α v 2 + k ¯ 3 C p λ 8 k v ( t ) + λ 9 k ¯ C p k v ( t ) + C p K λ 9 R L D α k v ( t ) 2 α Γ ( 1 α ) Ω 0 t k ( t s ) v ( t ) v ( s ) d s 0 t d ξ ( t ξ ) 1 α 0 ξ k ( η ) d η ( t η ) α 0 ξ v ( s ) d s ( t s ) α d x α Γ ( 1 α ) Ω 0 t d ξ ( t ξ ) 1 α 0 ξ ( t s ) α 0 s k ( s τ ) v ( s ) v ( τ ) d τ d s 2 d x
for λ 8 , λ 9 > 0 , t 0 . Further, using a similar argument to the one used to derive (20), for C 5 = α C 3 Γ ( 1 α ) and λ 8 , λ 9 , λ 10 > 0 , we obtain
D α Ψ 5 ( t ) λ 8 D α v 2 + C p λ 9 k ¯ + k ¯ 3 λ 8 + λ 10 α k ¯ Γ ( 1 α ) k v ( t ) + C p K λ 9 R L D α k v ( t ) + C 5 E α ( C k t α ) λ 10 Ω 0 t d ξ ( t ξ ) 1 α 0 ξ v ( s ) d s ( t s ) α 2 d x α Γ ( 1 α ) Ω 0 t d ξ ( t ξ ) 1 α 0 ξ ( t s ) α 0 s k ( s τ ) v ( s ) v ( τ ) d τ d s 2 d x , t 0 .
This ends the proof. □
Let
L 1 ( t ) : = N E 1 ( t ) + i = 1 5 N i Ψ i ( t ) , t > 0
for some positive constants N , N i , i = 1 , , 5 to be determined later.
Theorem 3.
Assume that k satisfies (K) and γ satisfies (Γ). Then, the solution to the problem in (8) (with a = 0 ) is locally Mittag–Leffler stable, that is there exist positive constants B 1 and ν 1 such that
E 1 ( t ) B 1 E α ( ν 1 t α ) , t 0
for small K (or, alternatively, a large C k ) .
Proof. 
The evaluations obtained in the previous lemmas, with N 1 = N , imply
D α L 1 ( t ) N K E α ( C k t α ) 2 λ 1 + N 3 K + N 4 λ 5 N 2 1 λ 2 v 2 N 2 v ( t ) γ + 2 γ + 2 + N 2 + λ 8 N 5 + N 4 λ 4 k 0 2 D α v 2 + 2 N λ 1 k ¯ + k ¯ N 2 4 λ 2 N 3 C k Ψ 3 ( t ) + λ 6 N 4 C 1 E 1 γ + 1 ( t ) + N 2 K N 4 C p 2 k 0 K N λ 1 C p N 5 K λ 9 ( R L D α k v ) ( t ) + N 4 k ¯ 1 4 λ 5 + C p 4 λ 6 + 1 + k ¯ N 5 C p λ 9 + k ¯ 2 λ 8 + λ 10 α Γ ( 1 α ) k v ( t ) + α 4 Γ ( 1 α ) N 2 λ 3 + 4 λ 7 N 4 2 N Ω 0 t d ξ ( t ξ ) 1 α 0 ξ D α v ( η ) d η ( t η ) α 2 d x + N 4 C 4 λ 4 + N 5 C 5 λ 10 C p E α ( C k t α ) + α 4 λ 3 C p N 2 N 4 Γ ( 1 α ) Ω 0 t d ξ ( t ξ ) 1 α 0 ξ v ( η ) d η ( t η ) α 2 d x + α 4 Γ ( 1 α ) N 4 λ 7 4 N 5 Ω 0 t d ξ ( t ξ ) 1 α 0 ξ ( t s ) α 0 s k ( s τ ) v ( s ) v ( τ ) d τ d s 2 d x .
The strategy consists of choosing the parameters so as to give rise to a term of the form C E 1 ( t ) , C > 0 , in addition to λ 6 N 4 C 1 E 1 γ + 1 ( t ) , on the right-hand side of (25). Observe that, as E α ( C k t α ) decreases towards zero, choosing t 0 large enough allows us to ignore this term. Incidentally, we may also ignore λ 4 and λ 10 . Let λ 1 = 1 , λ 2 = 1 / 2 and recall that k ¯ K / C k . Therefore, for small values of K (or large values of C k ) and consequently small λ 5 , λ 6 , λ 8 , and λ 9 , it suffices that
N 2 λ 3 + 4 λ 7 N 4 2 N , 4 λ 3 C p N 2 N , N 2 < N 4 k 0 2 , N 4 λ 7 4 N 5 .
To this end, we may pick λ 3 = λ 7 = 1 / 4 , N 5 = N 4 , N 2 = k 0 N 4 4 and N N 4 max C p k 0 4 , k 0 + 1 . Finally, we go back and select the remaining parameters to ensure the fractional inequality
D α L 1 ( t ) C 6 E 1 ( t ) + M 4 C 1 E 1 γ + 1 ( t ) , t > t 0
for some C 6 > 0 and by the equivalence of L 1 ( t ) and E 1 ( t )
D α L 1 ( t ) C 7 L 1 ( t ) + C 8 L 1 γ + 1 ( t ) , C 7 , C 8 , t > t 0 .
Let t ¯ , ω > 0 be real numbers (we can assume that t 0 < t ¯ ) such that L 1 ( t 0 ) < ω and C 8 L 1 γ ( t ) ω < C 7 , 0 t t ¯ so that
D α L 1 ( t ) C 7 ω L 1 ( t ) , t 0 t t ¯ .
Consequently,
L 1 ( t ) L 1 ( t 0 ) E α ( C 7 ω t t 0 α ) , t 0 t t ¯ .
Assume that there exists t 0 t 1 t ¯ such that L 1 ( t 1 ) = ω , then as
E α ( C 7 ω t t 0 α ) 1 ,
we have
ν = L 1 ( t 1 ) L 1 ( t 0 ) E α ( χ ω t 1 t 0 α ) < ω .
This is a contradiction. Starting from t ¯ , we can continue the process forever; thus, the proof is complete. □
Remark 1.
The case a 0 and b 0 , results from the above arguments. The functionals in the first case are good enough to obtain an explicit upper bound for the kernels ensuring Mittag–Leffler stability.

7. Conclusions

A half-order fractional derivative of the leading term is found to be capable of driving the system to rest. The obtained and expected rate is of Mittag–Leffler-type which is characteristic of the fractional derivatives. The stability is global in time. Therefore, this term is a damping term just as its integer-order counterpart. In the case of a memory term (viscoelastic), a serious problem arises from the evaluation of the expression
Ω 0 t d ξ ( t ξ ) 1 α 0 ξ ( t η ) α 0 s k ( s τ ) v ( τ ) d τ d η 2 d x
appearing in the energy derivative. To cope with this situation, we introduced a new functional Ψ 1 ( t ) . Regarding non- linearity, we assumed some addition assumptions on the exponent of non-linearity so as to profit from an embedding result. Unfortunately, the obtained fractional inequality differs from the one corresponding to the lower-order fractional term. To extract information about the concerned functional we had to consider small kernels and use an argument to deal with the arising non-linearity. The price to pay is a stability of the local character only.

Author Contributions

Conceptualization, N.-e.T.; Methodology, B.A.-H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by King Fahd University of Petroleum and Minerals through project number SB201014.

Acknowledgments

The authors are very grateful for the financial support and facilities provided by King Fahd University of Petroleum and Minerals through project number SB201014. The authors are also indebted to the anonymous referees whose comments and suggestions helped considerably improve the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Atanackovic, T.M.; Pilipovic, S.; Zorica, D. A diffusion wave equation with two fractional derivatives of different order. J. Phys. A Math. Theor. 2007, 40, 5319–5333. [Google Scholar] [CrossRef]
  2. Chen, J.; Liu, F.; Anh, V. Analytical solution for the time-fractional telegraph equation by the method of separating variables. J. Math. Anal. Appl. 2008, 338, 1364–1377. [Google Scholar] [CrossRef]
  3. Mainardi, F. Fractional Calculus and Waves in Linear Viscoelasticity: An Introduction to Mathematical Models; Imperial College Press: London, UK, 2010. [Google Scholar]
  4. Orsingher, E.; Beghin, L. Time-fractional telegraph equation and telegraph processes with Brownian time. Probab. Theory Relat. Fields 2004, 128, 141–160. [Google Scholar] [CrossRef]
  5. Povstenko, Y. Fractional Thermoelasticity; Springer: New York, NY, USA, 2015. [Google Scholar]
  6. Sandev, T.; Tomovski, Z. The general time fractional wave equation for a vibrating string. J. Phys. A Math. Theor. 2010, 43, 055204. [Google Scholar] [CrossRef]
  7. Schneider, W.R.; Wess, W. Fractional diffusion and wave equations. J. Math. Phys. 2006, 213, 205–213. [Google Scholar] [CrossRef]
  8. Shen, F.; Tan, W.; Zhao, Y.; Masuoka, T. The Raleigh-Stokes problem for a heated generalized second grade fluid with fractional derivative model. Nonlinear Anal. Real World Appl. 2006, 7, 1072–1080. [Google Scholar] [CrossRef]
  9. Barthelemy, P.; Bertolotti, J.; Wiersma, D.S. A Lévy flight for light. Nature 2008, 453, 495–498. [Google Scholar] [CrossRef]
  10. Goychuk, I.; Kharchenko, V.; Metzler, R. Molecular motors pulling cargos in the viscoelastic cytosol: How power strokes beat subdiffusion. Phys. Chem. Chem. Phys. 2014, 16, 16524–16535. [Google Scholar] [CrossRef]
  11. Tarasov, V.E. Fractional Dynamics: Applications of Fractional Calculus to Dynamics of Particles, Fields and Media; Springer: Berlin/Heidelberg, Germany, 2010. [Google Scholar]
  12. Lunardi, A. On the linear heat equation with fading memory. SIAM J. Math. Anal. 1990, 21, 1213–1224. [Google Scholar] [CrossRef]
  13. Fabrizio, M.; Morro, A. Mathematical Problems in Linear Viscoelasticity; Studies in Applied and Numerical Mathematics; Society for Industrial and Applied Mathematics: Philadelphia, PA, USA, 1992. [Google Scholar]
  14. Giorgi, C.; Pata, V. Asymptotic behavior of a nonlinear hyperbolic heat equation with memory. NoDEA Nonlinear Differ. Equ. Appl. 2001, 8, 157–171. [Google Scholar] [CrossRef]
  15. Munteanu, I. Stabilization of semilinear heat equations, with fading memory, by boundary feedbacks. J. Diff. Equ. 2015, 259, 454–472. [Google Scholar] [CrossRef]
  16. Cavalcanti, M.M.; Domingos Cavalcanti, V.N.; Prates Filho, J.S.; Soriano, J.A. Existence and uniform decay rates for viscoelastic problems with nonlocal boundary damping. Diff. Integral Equ. 2001, 14, 85–116. [Google Scholar]
  17. Pata, V. Exponential stability in linear viscoelasticity. Quart. Appl. Math. 2006, LXIV, 499–513. [Google Scholar] [CrossRef]
  18. Tatar, N.-e. Polynomial stability without polynomial decay of the relaxation function. Math. Meth. Appl. Sci. 2008, 31, 1874–1886. [Google Scholar] [CrossRef]
  19. Tatar, N.-e. On a large class of kernels yielding exponential stability in viscoelasticity. Appl. Math. Comp. 2009, 215, 2298–2306. [Google Scholar] [CrossRef]
  20. Tatar, N.-e. Arbitrary decays in viscoelasticity. J. Math. Physics 2011, 52, 013502. [Google Scholar] [CrossRef]
  21. Tatar, N.-e. A new class of kernels leading to an arbitrary decay in viscoelasticity. Mediterr. J. Math. 2013, 10, 213–226. [Google Scholar] [CrossRef]
  22. Tatar, N.-e. Nonexistence results for a fractional problem arising in thermal diffusion in fractal media. Chaos Solit. Fractals 2008, 36, 1205–1214. [Google Scholar] [CrossRef]
  23. El-Sayed, A.M.A.; Herzallah, M.A.E. Continuation and maximal regularity of an arbitrary (fractional) order evolutionary integro-Differential equation. Appl. Anal. 2005, 84, 1–14. [Google Scholar] [CrossRef]
  24. Ponce, R. Bounded mild solutions to fractional integro-differential equations in Banach spaces. Semigroup Forum 2013, 87, 377–392. [Google Scholar] [CrossRef]
  25. Agarwal, R.P.; Dos Santos, J.P.C.; Cuevas, C. Analytic resolvent operator and existence results for fractional integrodifferential equations. J. Abst. Diff. Eqs. Appl. 2012, 2, 26–47. [Google Scholar]
  26. Andrade, B.D.; Dos Santos, J.P.C. Existence of solutions for a fractional neutral integro-differential equation with unbounded delay. Electron. J. Diff. Eqs. 2012, 90, 1–13. [Google Scholar]
  27. Gorenflo, R.; Kilbas, A.A.; Mainardi, F.; Rogosin, S.Y. Mittag–Leffler Functions, Related Topics and Applications; Springer Monographs in Mathematics; Springer: Berlin/Heidelberg, Germeny, 2014. [Google Scholar] [CrossRef]
  28. Podlubny, I. Fractional Differential Equations; Academic Press: San Diego, CA, USA, 1999. [Google Scholar]
  29. Yan, L.; Chen, Y.Q.; Podlubny, I. Stability of fractional-order nonlinear dynamic systems: Lyapunov direct method and generalized Mittag Leffler stability. Comput. Math. Appl. 2010, 59, 1810–1821. [Google Scholar]
  30. Alikhanov, A.A. Boundary value problems for the diffusion equation of the variable order in differential and difference settings. Appl. Math. Comput. 2012, 219, 3938–3946. [Google Scholar] [CrossRef]
  31. Dai, H.; Chen, W. New power law inequalities for fractional derivative and stability analysis of fractional order systems. Nonlinear Dyn. 2017, 87, 1531–1542. [Google Scholar] [CrossRef]
  32. Zhang, S.; Yu, Y.; Wang, H. Mittag–Leffler stability of fractional-order Hopfield neural networks. Nonlinear Anal. Hybrid Syst. 2015, 16, 104–121. [Google Scholar] [CrossRef]
  33. Kukushkin, M.V. Abstract fractional calculus for m-accretive operators. Int. J. Appl. Math. 2001, 34, 1–41. [Google Scholar] [CrossRef]
  34. Shkalikov, A.A. Perturbations of self-adjoint and normal operators with discrete spectrum. Russ. Math. Surveys 2016, 71, 907–964. [Google Scholar] [CrossRef]
  35. Dos Santos, J.P.C. Fractional resolvent operator with α∈(0,1) and applications. Frac. Diff. Calc. 2019, 9, 187–208. [Google Scholar] [CrossRef]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Al-Homidan, B.; Tatar, N.-e. Stabilization of a Non-linear Fractional Problem by a Non-Integer Frictional Damping and a Viscoelastic Term. Fractal Fract. 2023, 7, 367. https://doi.org/10.3390/fractalfract7050367

AMA Style

Al-Homidan B, Tatar N-e. Stabilization of a Non-linear Fractional Problem by a Non-Integer Frictional Damping and a Viscoelastic Term. Fractal and Fractional. 2023; 7(5):367. https://doi.org/10.3390/fractalfract7050367

Chicago/Turabian Style

Al-Homidan, Banan, and Nasser-eddine Tatar. 2023. "Stabilization of a Non-linear Fractional Problem by a Non-Integer Frictional Damping and a Viscoelastic Term" Fractal and Fractional 7, no. 5: 367. https://doi.org/10.3390/fractalfract7050367

APA Style

Al-Homidan, B., & Tatar, N. -e. (2023). Stabilization of a Non-linear Fractional Problem by a Non-Integer Frictional Damping and a Viscoelastic Term. Fractal and Fractional, 7(5), 367. https://doi.org/10.3390/fractalfract7050367

Article Metrics

Back to TopTop