Next Article in Journal
An Experimental Study of Grouping Mutation Operators for the Unrelated Parallel-Machine Scheduling Problem
Previous Article in Journal
Analysis of Multi-Stacked Dielectric Resonator Antenna with Its Equivalent R-L-C Circuit Modeling for Wireless Communication Systems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Stability Results for a Weakly Dissipative Viscoelastic Equation with Variable-Exponent Nonlinearity: Theory and Numerics

by
Adel M. Al-Mahdi
1,2,*,
Mohammad M. Al-Gharabli
1,2,
Maher Noor
3 and
Johnson D. Audu
1
1
The Preparatory Year Program, King Fahd University of Petroleum and Minerals, Dhahran 31261, Saudi Arabia
2
The Interdisciplinary Research Center in Construction and Building Materials, King Fahd University of Petroleum & Minerals, Dhahran 31261, Saudi Arabia
3
DCC-Math, King Fahd University of Petroleum and Minerals, Dhahran 31261, Saudi Arabia
*
Author to whom correspondence should be addressed.
Math. Comput. Appl. 2023, 28(1), 5; https://doi.org/10.3390/mca28010005
Submission received: 3 December 2022 / Revised: 28 December 2022 / Accepted: 29 December 2022 / Published: 4 January 2023

Abstract

:
In this paper, we study the long-time behavior of a weakly dissipative viscoelastic equation with variable exponent nonlinearity of the form u t t + Δ 2 u 0 t g ( t s ) Δ u ( s ) d s + a | u t | n ( · ) 2 u t Δ u t = 0 , where n ( . ) is a continuous function satisfying some assumptions and g is a general relaxation function such that g ( t ) ξ ( t ) G ( g ( t ) ) , where ξ and G are functions satisfying some specific properties that will be mentioned in the paper. Depending on the nature of the decay rate of g and the variable exponent n ( . ) , we establish explicit and general decay results of the energy functional. We give some numerical illustrations to support our theoretical results. Our results improve some earlier works in the literature.

1. Introduction

In this paper, we consider the following weakly dissipative viscoelastic equation with variable exponent nonlinearity
u t t + Δ 2 u 0 t g ( t s ) Δ u ( s ) d s + a | u t | n ( · ) 2 u t Δ u t = 0 in Ω × ( 0 , + ) ,
subject to the following conditions
u = Δ u = 0 , on Ω , u ( x , 0 ) = u 0 ( x ) , u t ( x , 0 ) = u 1 ( x ) , in Ω ,
where Ω is a bounded domain of R d with a smooth boundary Ω , a is a positive constant, n ( . ) is a continuous function satisfying some assumptions, g is a general relaxation function satisfying some conditions and ( u 0 , u 1 ) are the given initial data.
For the stabilization of weakly dissipative second-order systems, Rivera et al. [1] considered the following abstract integro-differential equation
u t t + A u + β u 0 t g ( t s ) A α u d s = 0 ,
where A is a strictly positive self-adjoint linear operator and established a polynomial decay result with the interpolating cases α ( 0 , 1 ) and g decays exponentially to zero. In [2], Hassan and Messaoudi generalized the result of [1] by considering a very general assumption on the relaxation function g and obtained a new general decay rate. In [3], Anaya and Messaoudi discussed the following problem
u t t + Δ 2 u 0 t g ( t s ) Δ u ( s ) d s + h ( u t ) = 0 ,
with a general weak damping term, and they derived general decay rate estimates under certain restrictions for the relaxation function g, and the function h is a nondecreasing C 1 ( R ) function satisfying h ( 0 ) = 0 , s h ( s ) 0 and for c 1 , c 2 > 0 ,
c 1 | s | | h ( s ) | c 2 | s | , s R .
In recent years, there has been increasing interest in treating equations with variable exponents of nonlinearity. Some models from physical phenomena such as flows of electro-rheological fluids or fluids with temperature-dependent viscosity, filtration processes in porous media, nonlinear viscoelasticity, and image processing give rise to such problems. This great interest is motivated by the applications to the mathematical modeling of non-Newtonian fluids. One of these fluids is the electro-rheological fluids which have the ability to drastically change when applying some external electromagnetic field. The variable exponent of nonlinearity is a given function of density, temperature, saturation, electric field, etc. For more information about the mathematical model of electro-rheological fluids, we refer to [4,5]. However, there are few available works in the literature including nonlinearities of variable-exponent type. For example, Antontsev [6,7] considered the following problem
u t t d i v ( a ( x ) | u | p ( x , t ) 2 u ) α Δ u t = b ( x , t ) | u | σ ( x , t ) 2 u ,
where α is a positive constant and the author established local, global existence and blow-up results under some conditions on the functions a , b , p , σ . In [8], Problem (6) was also considered, and the authors proved various blow-up results for the solutions with positive initial energy. Furthermore, in [9], the authors studied the following equation
u t t d i v ( | u | r ( x ) 2 u ) + a | u t | n ( x ) 1 u t = b | u | p ( x ) 2 u ,
and proved the existence of a unique weak solution. The authors of [9] also established a finite-time blow up of the solutions. In [10], Messaoudi et al. considered Problem ( ) , where b = 0 , and established decay estimates for the solutions under suitable assumptions on the initial data and the variable exponents.
Recently, in [11], Messaoudi investigated the following problem
u t t d i v ( a ( x ) | u | r ( . ) 2 u ) Δ u t + | u t | n ( . ) 2 u t = 0 ,
and established exponential and polynomial decay results under some conditions on the variable exponents n and r. The existence, stability and blow up of solutions of other problems with variable exponents such as Petrovsky, Kirchhoff and other viscoelastic equations can be found in the references [12,13,14,15,16].
In the present work, we are interested in establishing explicit and general decay results for the system (1)–(2) by using the energy method and some properties related to the variable exponents. Then, we give some numerical examples to illustrate our results. Our decay results depend on the nature of the decay rate of the relaxation function g and the variable exponent n ( . ) . The results improve the recent works in [1,2,3], where the authors considered the usual constant exponents.
The remainder of this paper is organized as follows: In Section 2, we outline some preliminaries. In Section 3, we state and prove some essential lemmas which are needed in the proofs of the decay results. In Section 4, we establish the general decay results of our problem (1)–(2) and present some examples. Finally, in Section 5, we provide some numerical examples to illustrate our theoretical results.

2. Preliminaries

In this work, L 2 ( Ω ) denotes the standard Lesbesgue space, and we define following
H ( Ω ) : = { u H 3 ( Ω ) : u = Δ u = 0 on Ω }
is a Sobolev space with the usual scalar products and norms. Throughout this paper, k is a generic positive constant. Details on Lebesgue and Sobolev spaces with variable exponents can be found in [17,18,19]). Here, we provide some basic definitions. Let q : Ω [ 1 , ] be a measurable function, where Ω is a domain of R d . The Lebesgue space with a variable exponent q ( · ) is given by
L q ( · ) ( Ω ) : = u : Ω R ; measurable in Ω : ρ q ( · ) ( μ u ) < , for some μ > 0 ,
where
ρ q ( · ) ( u ) = Ω | u ( x ) | q ( x ) d x .
The space L q ( · ) ( Ω ) is Banach (see [18]) when equipped with the norm
u q ( · ) : = inf μ > 0 : Ω | u ( x ) μ | q ( x ) d x < , u L q ( · ) ( Ω ) .
Furthermore, L q ( · ) ( Ω ) is both reflexive and separable for q 2 ( q 1 , ) , provided q ( · ) is bounded with
q 1 : = essinf x Ω q ( x ) , q 2 : = esssup x Ω q ( x ) .
The variable-exponent Sobolev space defined by
W 1 , q ( · ) ( Ω ) = u L q ( · ) ( Ω ) such that u exists and | u | L q ( · ) ( Ω )
is also Banach when equipped with the norm u W 1 , q ( · ) ( Ω ) = u q ( · ) + u q ( · ) . What is more is that W 1 , q ( · ) ( Ω ) is separable and reflexive for 1 < q 1 q 2 < , provided q ( · ) is bounded.
Lemma 1
([18]). Let Ω be a bounded domain in R d with a smooth boundary Ω . Assume that p , q C ( Ω ¯ ) such that for all x Ω ¯ ,
1 < p 1 p ( x ) p 2 < + , 1 < q 1 q ( x ) q 2 < + , ,
and q ( x ) < p * ( x ) Ω ¯ with p * ( x ) =
d p * ( x ) d p * ( x ) , i f p 2 < d , + , i f p 2 d ,
then W 1 , p ( . ) ( Ω ) is continuously and compactly embedded in L q ( . ) ( Ω ) . Consequently, there exists a constant k e > 0 such that
u q k e u W 1 , p ( . ) , u W 1 , p ( . ) ( Ω ) .
In this work, we assume the following:
( C 1 )
The function g : R + R + is a C 1 non-increasing function satisfying
g ( 0 ) > 0 , 1 ω 0 0 g ( s ) d s = > 0
where
u 2 2 ω 0 Δ u 2 2 , u H ( Ω ) ,
and there exists a C 1 function G : ( 0 , ) ( 0 , ) which is strictly increasing and a strictly convex C 2 function on ( 0 , r ] , r g ( 0 ) , with G ( 0 ) = G ( 0 ) = 0 , such that
g ( t ) ξ ( t ) G ( g ( t ) ) , t 0 ,
where ξ is a positive non-increasing differentiable function.
( C 2 )
n : Ω ¯ [ 1 , ) is a continuous function such that
n 1 : = essinf x Ω n ( x ) , n 2 : = esssup x Ω n ( x ) .
and 1 < n 1 n ( x ) n 2 , where
n 2 < , d = 1 , 2 ; n 2 2 d d 2 , d 3 .
Furthermore, the exponent n ( · ) satisfies the log-Hölder continuity condition; that is
| n ( x ) n ( y ) | A log | x y | , for all x , y Ω , with | x y | < δ ,
for any A , δ > 0 .
The energy functional associated to (1)–(2) is given by
E ( t ) = 1 2 u t 2 2 + Δ u 2 2 0 t g ( s ) d s u 2 2 + ( g u ) ( t ) 1 2 u t 2 2 + 1 ω 0 0 t g ( s ) d s Δ u 2 2 + ( g u ) ( t ) 0 ,
for any t 0 , where for v L l o c 2 ( R + ; L 2 ( Ω ) ) ,
( g v ) ( t ) : = 0 t g ( t s ) v ( t ) v ( s ) 2 2 d s .
Lemma 2.
For any t 0 , the energy functional E ( t ) satisfies
E ( t ) = 1 2 ( g u ) ( t ) 1 2 g ( t ) u 2 2 a Ω | u t | n ( x ) d x u t 2 2 0 .
Proof. 
In view of the boundary conditions (2), we can deduce (14) simply by multiplying (1) by u t and integrate over Ω . □
Remark 1
([20]). Using ( C 1 ) , one can show that for any t [ 0 , t 0 ] and for some κ > 0 ,
g ( t ) ξ ( t ) G ( g ( t ) ) κ ξ ( t ) = κ g ( 0 ) ξ ( t ) g ( 0 ) κ g ( 0 ) ξ ( t ) g ( t )
and hence,
ξ ( t ) g ( t ) g ( 0 ) κ ( t ) , t [ 0 , t 0 ] .
Moreover, if G is a strictly increasing and strictly convex C 2 function on ( 0 , r ] , with G ( 0 ) = G ( 0 ) = 0 , then there is a strictly convex and strictly increasing C 2 function G : [ 0 , + ) [ 0 , + ) which is an extension of G . For instance, we can define G ¯ , for any t > r , by
G ¯ ( t ) : = G ( r ) 2 t 2 + G ( r ) G ( r ) r t + G ( r ) + G ( r ) 2 r 2 G ( r ) r .

3. Technical Lemmas

We establish some important lemmas in this section.
Lemma 3
([21]). Under Assumption ( C 1 ) , we have for any
v L l o c 2 [ 0 , + ) ; L 2 ( 0 , L ) ,
0 L 0 t g ( t s ) ( v ( t ) v ( s ) ) d s 2 d x K ε ( h ε v ) ( t ) , t 0 .
where
K ε : = 0 g 2 ( s ) ε g ( s ) g ( s ) d s and h ε ( t ) : = ε g ( t ) g ( t ) ,
for any 0 < ε < 1 .
Lemma 4.
Under the assumptions ( C 1 ) and ( C 2 ) , and for 0 < δ < 1 , the functional
I ( t ) : = Ω u u t d x
satisfies the estimates:
I ( t ) u t 2 2 + 2 ω 0 u t 2 2 4 Δ u 2 2 + k K ε ( h ε u ) ( t ) + Ω k δ ( x ) | u t | n ( x ) d x , n 1 2 .
and if 1 < n 1 < 2 ,
I ( t ) u t 2 2 + 3 ω 0 u t 2 2 4 Δ u 2 2 + k K ε ( h ε u ) ( t ) + k Ω ( 1 + k δ ( x ) ) | u t | n ( x ) d x + Ω | u t | n ( x ) n 1 1 ,
where h ε is given in Lemma 3 and k δ ( x ) will be given in the proof below.
Proof. 
Differentiating I and using (1)–(2), we obtain
I ( t ) = Ω u t 2 d x Ω u . u t d x Δ u 2 2 Ω u ( t ) 0 t g ( t s ) Δ u ( s ) d s d x a Ω u | u t | n ( x ) 2 u t d x .
Using Young’s inequality, we have
Ω u . u t d x δ 0 ω 0 Δ u 2 2 + 1 4 δ 0 u t 2 2 .
Now, recalling (10), applying Young’s inequality, using Lemma 3 and (10), we obtain
Ω u ( t ) 0 t g ( t s ) Δ u ( s ) d s d x = Ω u ( t ) 0 t g ( t s ) Δ u ( s ) Δ u ( t ) d s d x Ω u ( t ) 0 t g ( t s ) Δ u ( t ) d s d x = Ω u ( t ) . 0 t g ( t s ) u ( s ) u ( t ) d s d x + 0 t g ( s ) d s u 2 2 . 0 t g ( s ) d s u 2 2 + 2 ω 0 Ω | u ( t ) | 2 + ω 0 2 Ω 0 t g ( t s ) u ( s ) u ( t ) d s 2 d x . ω 0 0 t g ( s ) d s Δ u 2 2 + 2 ω 0 u ( t ) 2 2 + k K ε ( h ε u ) ( t ) . ( 1 ) Δ u 2 2 + 2 Δ u ( t ) 2 2 + k K ε ( h ε u ) ( t ) . 1 2 Δ u 2 2 + k K ε ( h ε u ) ( t ) .
Applying Young’s inequality with p ( x ) = n ( x ) n ( x ) 1 and p ( x ) = n ( x ) , we can estimate the last term in (19) as follows
| u t | n ( x ) 2 u t u δ | u | n ( x ) + k δ ( x ) | u t | n ( x ) , x Ω ,
where
k δ ( x ) = δ 1 n ( x ) ( n ( x ) ) n ( x ) ( n ( x ) 1 ) n ( x ) 1 .
Hence,
Ω u | u t | n ( x ) u t d x δ Ω | u | n ( x ) d x + Ω k δ ( x ) | u t | n ( x ) d x .
To establish (17), we set
Ω + = { x Ω : | u ( x , t ) | 1 } and Ω = { x Ω : | u ( x , t ) | < 1 } .
Then, using (10), (13), (14), (22) and Lemma 1, we obtain
Ω | u | n ( x ) d x Ω + | u | n ( x ) d x + Ω | u | n ( x ) d x Ω + | u | n 2 d x + Ω | u | n 1 d x Ω | u | n 2 d x + Ω | u | n 1 d x k e n 1 | | u | | 2 n 1 + k e n 2 | | u | | 2 n 2 k e n 1 ω 0 ρ n 1 | | u | | 2 n 1 + k e n 2 ω 0 n 2 | | u | | 2 n 2 k e n 1 ω 0 n 1 | | u | | 2 n 1 2 + k e n 2 ω 0 n 2 | | u | | 2 n 2 2 | | u | | 2 2 k e n 1 ω 0 n 1 2 E ( 0 ) n 1 2 + k e n 2 ω 0 n 2 2 E ( 0 ) n 2 2 | | u | | 2 2 k 0 | | u | | 2 2 ,
where k 0 = k e n 1 ω 0 n 1 2 E ( 0 ) n 1 2 + k e n 2 ω 0 n 2 2 E ( 0 ) n 2 2 .
By fixing δ 0 = 8 ω 0 and δ = 8 k 0 , k δ ( x ) is still bounded; hence, we obtain the required estimate (17). To prove (18), we re-write the fourth term in (19) as follows
a Ω u | u t | n ( x ) 2 u t d x = a Ω 1 u | u t | n ( x ) 2 u t d x a Ω 2 u | u t | n ( x ) 2 u t d x ,
where
Ω 1 = { x Ω : n ( x ) < 2 } and Ω 2 = { x Ω : n ( x ) 2 } .
We notice that on Ω 1 , we have
2 n ( x ) 2 < n ( x ) , and 2 n ( x ) 2 2 n 1 2 .
Therefore, using Young’s and Poincaré’s inequalities and (25), we obtain
Ω 1 u | u t | n ( x ) 2 u t d x η Ω 1 | u | 2 d x + 1 4 η Ω 1 | u t | 2 n ( x ) 2 d x η k ρ 2 | | u | | 2 2 + k Ω 1 + | u t | 2 n ( x ) 2 d x + Ω 1 | u t | 2 n ( x ) 2 d x η k ρ 2 ω 0 | | Δ u | | 2 2 + k Ω 1 + | u t | n ( x ) d x + Ω 1 | u t | 2 n 1 2 d x η k ρ 2 ω 0 | | Δ u | | 2 2 + k Ω | u t | n ( x ) d x + Ω 1 | u t | 2 d x n 1 1 η k ρ 2 ω 0 | | Δ u | | 2 2 + k Ω | u t | n ( x ) d x + Ω 1 | u t | n ( x ) d x n 1 1 η k ρ 2 ω 0 | | Δ u | | 2 2 + k Ω | u t | n ( x ) d x + Ω | u t | n ( x ) d x n 1 1 ,
where
Ω 1 + = { x Ω 1 : | u t ( x , t ) | 1 } and Ω 1 = { x Ω 1 : | u t ( x , t ) | < 1 } .
Fixing η = 8 c ρ 2 ω 0 , (26) becomes
Ω 1 u | u t | n ( x ) 2 u t d x 8 | | Δ u | | 2 2 + c Ω | u t | n ( x ) d x + Ω | u t | n ( x ) d x n 1 1 .
Similarly, we set
Ω 2 + = { x Ω 2 : | u t ( x , t ) | [ 1 , ) } and Ω 2 = { x Ω 2 : | u t ( x , t ) | [ 0 , 1 ) } .
Therefore,
a Ω 2 u | u t | n ( x ) 2 u t d x δ Ω 2 | u | n ( x ) d x + Ω 2 k δ ( x ) | u t | n ( x ) d x δ Ω 2 + | u | n ( x ) d x + δ Ω 2 | u | n ( x ) d x + Ω k δ ( x ) | u t | n ( x ) d x δ Ω 2 + | u | n 2 d x + δ Ω 2 | u | n 1 d x + Ω k δ ( x ) | u t | n ( x ) d x δ Ω | u | n 2 d x + δ Ω | u | n 1 d x + Ω 2 k δ ( x ) | u t | n ( x ) d x δ k e n 1 | | u | | 2 n 1 + k e n 2 | | u | | 2 n 2 + Ω k δ ( x ) | u t | n ( x ) d x δ k e n 1 ω 0 n 1 | | u | | 2 n 1 2 + k e n 2 ω 0 n 2 | | u | | 2 n 2 2 | | u | | 2 2 + Ω k δ ( x ) | u t | n ( x ) d x δ k e n 1 ω 0 n 1 2 E ( 0 ) n 1 2 + k e n 2 ω 0 n 2 2 E ( 0 ) n 2 2 | | u | | 2 2 + Ω k δ ( x ) | u t | n ( x ) d x δ k 0 | | u | | 2 2 + Ω k δ ( x ) | u t | n ( x ) d x .
Combining (24)–(30) with (19) and fixing δ = 8 k 0 , we obtain (18). □
Lemma 5.
Suppose that the assumptions ( C 1 ) and ( C 2 ) are satisfied; then, the functional
J ( t ) : = 0 t f ( t s ) u ( s ) 2 2 d s
satisfies the following estimate:
J ( t ) 3 ( 1 ) Δ u 2 2 1 2 ( g u ) ( t ) ,
where f ( t ) = t g ( s ) d s .
Proof. 
Applying Young’s inequality, ( C 1 ) and the fact that f ( t ) f ( 0 ) = 1 ω 0 , we obtain, for any t 0 ,
J ( t ) = f ( 0 ) u ( t ) 2 2 0 t g ( t s ) u ( s ) 2 2 d s = f ( t ) u ( t ) 2 2 0 t g ( t s ) u ( s ) u ( t ) 2 2 d s 2 Ω 0 t g ( t s ) u ( t ) . u ( s ) u ( t ) d s d x f ( 0 ) u ( t ) 2 2 g u ( t ) + 2 ω 0 ( 1 ) u ( t ) 2 2 + 1 2 g u ( t ) = 3 ω 0 ( 1 ) u ( t ) 2 1 2 g u ( t ) 3 ( 1 ) Δ u ( t ) 2 2 1 2 g u ( t ) .
Lemma 6.
Assume that n 1 2 . Then, the functional L defined by
L ( t ) : = N E ( t ) + ε 1 I ( t )
satisfies the following equivalence relation
L E ,
and for some positive constants N , ε 1 , the functional satisfies the following estimate
L ( t ) k u t 2 2 4 ( 1 ) Δ u 2 2 + 1 4 ( g u ) ( t ) , t [ 0 , ) .
Proof. 
The proof of (32) is completed in [22]. For the proof of (33), combining (14) and (17), and using the fact that g ( t ) : = ε g ( t ) h ε ( t ) yields
L ( t ) N ε 1 k p 3 ε 1 ω 0 u t 2 2 l 4 ε 1 Δ u 2 2 + N ε 2 ( g u ) ( t ) N 2 k K ε ( h ε Δ u ) ( t ) Ω a N k ε 1 | u t | n ( x ) d x .
First, we select ε 1 = N 2 k ρ + 3 ω 0 ; then, we choose
N > max { ε 1 k p + 3 ε 1 ω 0 , k ε 1 a } .
Finally, we set ε = 1 4 N . Thus, (33) is established. □
Lemma 7.
Assume that n 1 ( 1 , 2 ) ; then, the functional L defined by
L ( t ) : = N E ( t ) + ε 1 I ( t )
satisfies
L E ,
for a suitable choice of the positive constants N , ε 1 , the functional satisfies the following estimate
L ( t ) k u t 2 2 4 ( 1 ) Δ u 2 2 + 1 4 ( g u ) ( t ) + k E ( t ) n 1 1 , t [ 0 , ) .
Proof. 
The proof is similar to the above arguments. □
Lemma 8.
Assume that ( C 1 ) and ( C 2 ) hold; then, we have, for n 1 2 ,
0 E ( s ) d s < .
Proof. 
Using Lemmas 5 and 7 we see that the functional L 1 defined by
L 1 ( t ) : = L ( t ) + J ( t ) ,
satisfies, for any t [ 0 , ) and for some positive constant k 1 ,
L 1 ( t ) k 1 u t 2 2 ( 1 ) Δ u 2 2 1 4 g u ( t ) k E ( t ) .
We then conclude that
0 E ( s ) d s < + .
Lemma 9.
Assume that ( C 1 ) and ( C 2 ) hold; then, for 1 < n 1 < 2 , we obtain
0 E 1 n 1 1 ( s ) d s < + .
Moreover,
t 0 E ( s ) d s k ( t t 0 ) 2 n 1 , t t 0 .
Proof. 
Using Lemmas (5) and (7), we obtain that the functional L 1 defined by
L 1 ( t ) : = L ( t ) + J ( t )
is non-negative and satisfies, for some t [ t 0 , ) and for some positive constants k 0 , k 1 ,
L 1 ( t ) k 1 u t 2 2 ( 1 ) Δ u 2 2 1 4 g u ( t ) + K Ω | u t | n ( x ) d x n 1 1 k 0 E ( t ) + k 1 E ( t ) n 1 1 .
Multiplying (40) by E q ( t ) , q = 2 n 1 n 1 1 > 0 , and using Young’s inequality, we obtain:
E q ( t ) L 1 ( t ) k 0 E q + 1 ( t ) + k 1 E q ( t ) E ( t ) n 1 1 k 0 ( 1 ε ) E q + 1 ( t ) + k ( ε ) E ( t ) .
Taking ε to be small enough and the fact that E is non-increasing, (41) becomes:
E q + 1 ( t ) k L 1 ( t ) ,
where L 1 ( t ) = E q ( t ) L 1 ( t ) + k E ( t ) . From this estimate, we conclude that
0 E q + 1 ( s ) d s < + .
Since q + 1 = 1 n 1 1 , we end up with
0 E 1 n 1 1 ( s ) d s < + .
Furthermore, from (44) and Hölder’s inequality, we can deduce that
t 0 t E ( s ) d s ( t t 0 ) q q + 1 t 0 t E q + 1 ( s ) d s 1 q + 1 k ( t t 0 ) q q + 1 = k ( t t 0 ) 2 n 1 , t t 0 .

4. Decay Results

We begin with the statement of our main theorem.
Theorem 1
(The case: n 1 2 ). Assume that ( C 1 ) and ( C 2 ) hold and the initial data ( u 0 , u 1 ) H ( Ω ) × H 0 1 ( Ω ) ; then, there exist constants λ 1 , λ 2 ( 0 , ) , such that the energy functional E ( t ) associated to Problem (1)–(2) satisfies the estimate
E ( t ) λ 2 G 0 1 λ 1 t 0 t ξ ( s ) d s , t > t 0 ,
where G 0 ( τ ) = τ G ( τ ) .
Proof. 
Since ξ is non-increasing, with the identity (14) and inequality (15),
0 t 0 g ( s ) ( u ( t ) u ( t s ) ) 2 2 d s 1 ξ ( t 0 ) 0 t 0 ξ ( s ) g ( s ) ( u ( t ) u ( t s ) ) 2 2 d s g ( 0 ) a ξ ( t 0 ) 0 t 0 ( s ) ( u ( t ) u ( t s ) ) 2 2 d s k ˜ E ( t ) , t [ t 0 , ) .
Combining (47) and (33), we obtain
L ( t ) n E ( t ) k E ( t ) + k t 0 t g ( s ) ( u ( t ) u ( t s ) ) 2 2 d s , t [ t 0 , ) .
Now, define a functional η as
η ( t ) : = γ t 0 t ( u ( t ) u ( t s ) ) 2 2 d s , t t 0 .
Using the inequality,
u 2 2 ω 0 Δ u 2 2 , u H ( Ω ) ,
we deduce that
E ( t ) 2 Δ u ( t ) 2 2 and E ( t ) 2 ω 0 u ( t ) 2 2 , t 0 .
These estimates and (37) yield
t 0 t ( u ( t ) u ( t s ) ) 2 2 d s 2 t 0 t u ( t ) 2 2 + u ( t s ) 2 2 d s 4 ω 0 t 0 t E ( t ) + E ( t s ) d s 8 ω 0 t 0 E ( s ) d s < + , t t 0 .
So, with γ ( 0 , 1 ) , we obtain
η ( t ) ( 1 , ) , t [ t 0 , ) .
Let θ be another functional defined by
θ ( t ) : = t 0 t g ( s ) ( u ( t ) u ( t s ) ) 2 2 d s .
In view of estimate (14), we can observe that E ( t ) 1 2 ( g u ) ( t ) . Therefore,
θ ( t ) k E ( t ) , t [ t 0 , ) .
Next, the facts that G is strictly convex and G ( 0 ) = 0 give that
G ( s τ ) s G ( τ ) , for s [ 0 , 1 ] and τ ( 0 , r ] .
In view of the assumptions ( C 1 ) , ( C 2 ) , (50) and Jensen’s inequality, we obtain
θ ( t ) = 1 η ( t ) t 0 t η ( t ) ( s ) ( u ( t ) u ( t s ) ) 2 2 d s 1 η ( t ) t 0 t η ( t ) ξ ( s ) G ( g ( s ) ) ( u ( t ) u ( t s ) ) 2 2 d s ξ ( t ) η ( t ) t 0 t G ( η ( t ) g ( s ) ) ( u ( t ) u ( t s ) ) 2 2 d s ξ ( t ) γ G γ t 0 t g ( s ) ( u ( t ) u ( t s ) ) 2 2 d s ξ ( t ) γ G ¯ γ t 0 t g ( s ) ( u ( t ) u ( t s ) ) 2 2 d s , t [ t 0 , ) .
This yields, for any t t 0 ,
t 0 t g ( s ) ( u ( t ) u ( t s ) ) 2 2 d s 1 γ G ¯ 1 γ θ ( t ) ξ ( t ) .
Therefore, (48) becomes
F ( t ) n E ( t ) + k γ G ¯ 1 γ θ ( t ) ξ ( t ) , t [ t 0 , )
where F : = L + k E .
Let r 1 ( 0 , r ) , and define a functional F 1 by
F 1 ( t ) : = G ¯ r 1 E ( t ) E ( 0 ) F ( t ) , t t 0 .
Then, using the facts that E 0 , G > 0 and G > 0 together with estimate (52) lead to
F 1 ( t ) = r 1 E ( t ) E ( 0 ) G ¯ r 1 E ( t ) E ( 0 ) F ( t ) + G ¯ r 1 E ( t ) E ( 0 ) F ( t ) n E ( t ) G ¯ r 1 E ( t ) E ( 0 ) + k γ G ¯ r 1 E ( t ) E ( 0 ) G ¯ 1 γ θ ( t ) ξ ( t ) , t t 0 .
G ¯ * is the convex conjugate of G ¯ (see ([23], pp. 61–64)), that is
G ¯ * ( s ) = s G ¯ 1 ( s ) G ¯ G ¯ 1 ( s )
and satisfies the generalized Young inequality
A B G ¯ * ( A ) + G ¯ ( B ) .
Now, set
A = G ¯ r 1 E ( t ) E ( 0 ) and B = G ¯ 1 γ θ ( t ) ξ ( t ) ,
then, in view of (53) and (55), we arrive at
F 1 ( t ) n E ( t ) G ¯ r 1 E ( t ) E ( 0 ) + k γ G ¯ * G ¯ r 1 E ( t ) E ( 0 ) + k θ ( t ) ξ ( t ) n ( E ( 0 ) k r 1 ) E ( t ) E ( 0 ) G ¯ r 1 E ( t ) E ( 0 ) + k θ ( t ) ξ ( t ) , t t 0 .
Fixing r 1 , we obtain
F 1 ( t ) n 1 E ( t ) E ( 0 ) G ¯ r 1 E ( t ) E ( 0 ) + k θ ( t ) ξ ( t ) , t t 0 ,
where n 1 > 0 . Multiplying both sides of (56) by ξ ( t ) r 1 E ( t ) E ( 0 ) < r and using inequality (51), yields
ξ ( t ) F 1 ( t ) n 1 E ( t ) E ( 0 ) G r 1 E ( t ) E ( 0 ) ξ ( t ) + k θ ( t ) n 1 E ( t ) E ( 0 ) G r 1 E ( t ) E ( 0 ) ξ ( t ) k E ( t ) , t t 0 .
Let F 2 = ξ F 1 + k E ; then, we obtain from the non-increasing property of ξ that
n 1 E ( t ) E ( 0 ) G r 1 E ( t ) E ( 0 ) ξ ( t ) F 2 ( t ) , t t 0 .
The map
t E ( t ) G ε 1 E ( t ) E ( 0 )
is non-increasing, since G > 0 and E non-increasing. As a result, integrating (57) over ( t 0 , t ) yields
n 1 E ( t ) E ( 0 ) G r 1 E ( t ) E ( 0 ) t 0 t ξ ( s ) d s t 0 t E ( s ) E ( 0 ) G r 1 E ( s ) E ( 0 ) ξ ( s ) d s F 2 ( t 0 ) F 2 ( t ) F 2 ( t 0 ) , t t 0 .
Finally, we set G 0 ( τ ) = τ G ( τ ) . Then, we obtain for some positive constants λ 1 and λ 2 the following estimate
E ( t ) λ 2 G 0 1 λ 1 t 0 t ξ ( s ) d s , t > t 0 .
We consider the following examples:
Example 1.
(1) Take g ( t ) = λ e β t , t [ 0 , ] , λ , β ( 0 , ) as constants. The constant λ is carefully chosen to satisfy the assumption ( C 1 ) . Consequently,
g ( t ) = β G ( g ( t ) ) , ξ ( t ) = β and G ( s ) = s .
Hence, in view of Theorem 1, we conclude that for some constant K , t 0 ( 0 , ) ,
E ( t ) K t t 0 , t [ t 0 , ) .
(2) 
Let g ( t ) = λ e ( 1 + t ) γ , for t [ 0 , ) , γ ( 0 , 1 ) and λ selected such that ( C 1 ) is satisfied. Then,
g ( t ) = ξ ( t ) G ( g ( t ) ) , ξ ( t ) = γ ( 1 + t ) γ 1 and G ( s ) = s .
In view of Theorem 1, we deduce that for some constant K , t 0 ( 0 , ) ,
E ( t ) K ( 1 + t ) γ , for t large enough .
(3) 
For γ ( 1 , ) , let
g ( t ) = λ ( 1 + t ) γ , t [ 0 , )
and λ carefully chosen so that ( C 1 ) is valid. Then,
g ( t ) = β G ( g ( t ) ) , ξ ( t ) = β and G ( s ) = s p ,
with p = 1 + γ γ ( 1 , 2 ) , and β is a positive constant. It follows from Theorem 1 that for some constants K , t 0 ( 0 , ) ,
E ( t ) K ( 1 + t ) γ / ( γ + 1 ) , t > t 0 .
Theorem 2
(The case: 1 < n 1 < 2 ). Assume that hypotheses ( C 1 ) and ( C 2 ) hold and the data ( u 0 , u 1 ) H ( Ω ) × H 0 1 ( Ω ) . Then, there exist positive constants λ 1 , λ 2 such that the energy functional associated to Problem (1)–(2) satisfies the estimate
E ( t ) λ 2 ( t t 0 ) 2 n 1 G 0 1 λ 1 ( t t 0 ) 2 n 1 n 1 1 t 0 t ξ ( s ) d s , t > t 0 ,
where G 0 ( τ ) = τ G ( τ ) .
Proof. 
Similar to the proof of Theorem 1, we use (47) and (33) to obtain
L ( t ) n E ( t ) k E ( t ) + k t 0 t g ( s ) ( u ( t ) u ( t s ) ) 2 2 d s + k E ( t ) n 1 1 , t [ t 0 , ) .
We then define another functional η as
η ( t ) : = γ ( t t 0 ) 2 n 1 t 0 t ( u ( t ) u ( t s ) ) 2 2 d s , t t 0 .
Using (39), we conclude that
η ( t ) k γ ,
then choosing 0 < γ < 1 small enough so that
η ( t ) < 1 , t t 0 .
Combining this with the hypotheses ( C 1 ) , ( C 2 ) , Jensen’s inequality and (61), we obtain
θ ( t ) = 1 η ( t ) t 0 t η ( t ) ( s ) ( u ( t ) u ( t s ) ) 2 2 d s 1 η ( t ) t 0 t η ( t ) ξ ( s ) G ( g ( s ) ) ( u ( t ) u ( t s ) ) 2 2 d s ξ ( t ) η ( t ) t 0 t G ( η ( t ) g ( s ) ) ( u ( t ) u ( t s ) ) 2 2 d s ( t t 0 ) 2 n 1 ξ ( t ) γ G γ ( t t 0 ) 2 n 1 t 0 t g ( s ) ( u ( t ) u ( t s ) ) 2 2 d s ( t t 0 ) 2 n 1 ξ ( t ) γ G ¯ γ ( t t 0 ) 2 n 1 t 0 t g ( s ) ( u ( t ) u ( t s ) ) 2 2 d s , t ( t 0 , ) ,
where G ¯ is a C 2 extension of G which is strictly increasing and strictly convex on ( 0 , ) . This yields, for any t t 0 ,
t 0 t g ( s ) ( u ( t ) u ( t s ) ) 2 2 d s ( t t 0 ) 2 n 1 γ G ¯ 1 γ θ ( t ) ( t t 0 ) 2 n 1 ξ ( t )
and (59) becomes
F ( t ) n E ( t ) + ( t t 0 ) 2 n 1 γ G ¯ 1 γ θ ( t ) ( t t 0 ) 2 n 1 ξ ( t ) + E ( t ) n 1 1 , t t 0 ,
where F : = L + k E . Let 0 < r 1 < r ; then, define a functional F 1 by
F 1 ( t ) : = G ¯ r 1 ( t t 0 ) 2 n 1 . E ( t ) E ( 0 ) F ( t ) , t t 0 .
Since E 0 , G > 0 and G > 0 , we obtain
F 1 ( t ) = ( n 2 ) r 1 n ( t t 0 ) n + 2 n E ( t ) E ( 0 ) + r 1 ( t t 0 ) 2 n 1 E ( t ) E ( 0 ) G ¯ r 1 ( t t 0 ) 2 n 1 E ( t ) E ( 0 ) F ( t ) + G ¯ r 1 ( t t 0 ) 2 n 1 E ( t ) E ( 0 ) F ( t ) G ¯ r 1 ( t t 0 ) 2 n 1 . E ( t ) E ( 0 ) F ( t ) , t t 0 .
Estimates (63) and (64) imply that
F 1 ( t ) n E ( t ) G ¯ r 1 ( t t 0 ) 2 n 1 . E ( t ) E ( 0 ) + ( t t 0 ) 2 n 1 γ G ¯ 1 γ θ ( t ) ( t t 0 ) 2 n 1 ξ ( t ) G ¯ r 1 ( t t 0 ) 2 n 1 . E ( t ) E ( 0 ) + c G ¯ r 1 ( t t 0 ) 2 n 1 . E ( t ) E ( 0 ) E ( t ) n 1 1 , t t 0 .
Let G ¯ * be defined as in (54) and satisfy (55). Set
A = G ¯ r 1 ( t t 0 ) 2 n 1 . E ( t ) E ( 0 ) and B = G ¯ 1 γ θ ( t ) ( t t 0 ) 2 n 1 ξ ( t ) ;
then, it follows from a combination of (55) and (64) that
F 1 ( t ) n E ( t ) G ¯ r 1 ( t t 0 ) 2 n 1 . E ( t ) E ( 0 ) + k ( t t 0 ) 2 n 1 γ G ¯ * G ¯ r 1 ( t t 0 ) 2 n 1 . E ( t ) E ( 0 ) + θ ( t ) ξ ( t ) + k G ¯ r 1 ( t t 0 ) 2 n 1 . E ( t ) E ( 0 ) E ( t ) n 1 1 .
By definition of G ¯ * and since G ¯ > 0 , we have:
F 1 ( t ) n E ( t ) G ¯ r 1 ( t t 0 ) 2 n 1 . E ( t ) E ( 0 ) + k r 1 E ( t ) E ( 0 ) G ¯ r 1 ( t t 0 ) 2 n 1 . E ( t ) E ( 0 ) + θ ( t ) ξ ( t ) + k G ¯ r 1 ( t t 0 ) 2 n 1 . E ( t ) E ( 0 ) E ( t ) n 1 1 ,
and it can be written as
F 1 ( t ) n ( E ( 0 ) k r 1 ) E ( t ) E ( 0 ) G ¯ r 1 ( t t 0 ) 2 n 1 . E ( t ) E ( 0 ) + θ ( t ) ξ ( t ) + k G ¯ r 1 ( t t 0 ) 2 n 1 . E ( t ) E ( 0 ) E ( t ) n 1 1 .
After fixing r 1 , we arrive at
F 1 ( t ) n 1 E ( t ) E ( 0 ) G ¯ r 1 ( t t 0 ) 2 n 1 . E ( t ) E ( 0 ) + θ ( t ) ξ ( t ) + k G ¯ r 1 ( t t 0 ) 2 n 1 . E ( t ) E ( 0 ) E ( t ) n 1 1 ,
where n 1 > 0 . Multiplying both sides of (66) by ξ ( t ) E 2 n 1 n 1 1 ( t ) , we reach
ξ ( t ) E 2 n 1 n 1 1 ( t ) F 1 ( t ) n 1 E ( t ) E ( 0 ) ξ ( t ) E 2 n 1 n 1 1 ( t ) G ¯ r 1 ( t t 0 ) 2 n 1 . E ( t ) E ( 0 ) k 1 ( E ( t ) ) E 2 n 1 n 1 1 ( t ) + k ξ ( t ) G ¯ r 1 ( t t 0 ) 2 n 1 . E ( t ) E ( 0 ) E 2 n 1 n 1 1 ( t ) E ( t ) n 1 1 .
Using Young’s inequality with p = 1 n 1 1 and q = 1 2 n 1 , we obtain
ξ ( t ) E 2 n 1 n 1 1 ( t ) F 1 ( t ) ( k ε n 1 E ( 0 ) ) ξ ( t ) E 2 n 1 n 1 1 ( t ) G ¯ r 1 ( t t 0 ) 2 n 1 . E ( t ) E ( 0 ) k 1 ( E ( t ) + E 2 n 1 n 1 1 ( t ) ) k 2 E ( t ) .
Let F 2 = ξ F 1 E 2 n 2 n 2 + k 1 E E 2 n 2 n 2 + k 2 E ( t ) ; then, we obtain, from the non-increasing property of ξ and the fact that E 0 and for ε small enough that
F 2 ( t ) n 2 ξ ( t ) E 2 n 1 n 1 1 ( t ) G ¯ r 1 ( t t 0 ) 2 n 1 . E ( t ) E ( 0 ) ,
for some n 2 > 0 . Then, we have for n 3 = n 2 E ( 0 ) and r 1 ( t t 0 ) 2 n 1 . E ( t ) E ( 0 ) < r for small r 1 ,
n 3 E 2 n 1 n 1 1 ( t ) E ( 0 ) G r 1 ( t t 0 ) 2 n 1 . E ( t ) E ( 0 ) ξ ( t ) F 2 ( t ) .
Integrating (69) over ( t 0 , t ) yields
t 0 t n 3 E 2 n 1 n 1 1 ( s ) E ( 0 ) G r 1 ( s t 0 ) 2 n 1 . E ( s ) E ( 0 ) ξ ( s ) d s t 0 t F 2 ( s ) d s F 2 ( t 0 ) .
It follows from the fact that G > 0 and non-increasing property of E ( t ) that the map
t E 2 n 1 n 1 1 ( s ) E ( 0 ) r 1 ( s t 0 ) 2 n 1 . E ( s ) E ( 0 )
is non-increasing. Consequently, we have
n 3 E 2 n 1 n 1 1 ( t ) E ( 0 ) G r 1 ( t t 0 ) 2 n 1 . E ( t ) E ( 0 ) t 0 t ξ ( s ) d s t 0 t n 3 E 2 n 1 n 1 1 ( s ) E ( 0 ) G r 1 ( s t 0 ) 2 n 1 . E ( s ) E ( 0 ) ξ ( s ) d s t 0 t F 2 ( s ) d s F 2 ( t 0 ) = n 4 .
To finish the proof of Theorem 2, we multiply (71) by 1 ( t t 0 ) 2 n 1 n 1 1 to obtain
n 5 E ( t ) E ( 0 ) ( t t 0 ) 2 n 1 n 1 1 G r 1 ( t t 0 ) 2 n 1 . E ( t ) E ( 0 ) t 0 t ξ ( s ) d s n 4 t t 0 2 n 1 n 1 1 .
Next, we set G 0 ( τ ) = τ 1 n 1 1 G ( τ ) which is strictly increasing; then, we obtain for two positive constants λ 1 and λ 2
E ( t ) λ 2 ( t t 0 ) 2 n 1 G 0 1 λ 1 ( t t 0 ) 2 n 1 n 1 1 t 0 t ξ ( s ) d s , t > t 0 .
Example 2.
(1) Let g ( t ) = λ e β ( 1 + t ) γ , t [ 0 , ) , λ , β ( 0 , ) and γ ( 0 , 1 ) and λ is selected such that ( C 1 ) is satisfied; then, g ( t ) = β G ( g ( t ) ) with ξ ( t ) = γ ( 1 + t ) γ 1 and G ( s ) = s . So, in view of Theorem 2, we conclude that the solution of (1)–(2) satisfies the energy estimate
E ( t ) K ( t t 0 ) n 1 1 , t > t 1 .
(2) 
Suppose
g ( t ) = λ ( 1 + t ) γ , γ ( 1 , ) ,
and λ is chosen so that ( C 1 ) is satisfied. Then, for a positive constant β,
g ( t ) = β G ( g ( t ) ) , ξ ( t ) = β and G ( s ) = s p , p = 1 + γ γ .
Similar to the arguments in Example 1, we deduce that the solution of (1)–(2) satisfies the energy estimate
E ( t ) K ( t t 0 ) μ ,
where μ = ( n 1 1 ) ( n 1 + γ 2 ) n 1 + γ 1 ( 0 , ) , and for sufficiently large t and some constant K , t 0 ( 0 , ) .

5. Numerical Results

We produce some numerical experiments in this section to demonstrate the theoretical findings of Theorems 1 and 2. For this reason, we discretize our Problem (1) in the time–space domain ( 0 , 1 ] × [ 0 , 1 ] using the second-order finite difference method (FDM) in time and fourth-order in space. The time interval ( 0 , T ) is split into N = 10,000 subintervals with a time step Δ t = T N , and the spatial interval ( 0 , 1 ) is divided into 50 subintervals.
The homogeneous Dirichlet boundary condition for Problem (1) is stated, and Δ u = 0 at boundary. Based on the relaxation function G and the initial conditions u ( x , 0 ) = s i n ( π x ) and u t ( x , 0 ) = 0 , we contrast the following numerical two tests.
  • Test 1: We show the exponential decay of the solution of the problem (1) and the energy function given by (13), using g ( t ) = e t and n ( x ) = 2 + 1 x + 1 .
  • Test 2: In the second test, we let g ( t ) = 1 ( t + 1 ) 2 and n ( x ) = 3 + 1 x + 1 .
In Figure 1 and Figure 2, we show the cross-sections of the approximate solution u at x = 0.3 , x = 0.5 , x = 0.6 , and x = 0.7 for Test 1 and Test 2, respectively. In Figure 3 and Figure 4, we graph the corresponding energy functional (13). In addition, we show the decay behavior of the whole wave over the time interval [ 0 , 1 ] in Figure 5 and Figure 6 for Test 1 and Test 2, respectively.

6. Conclusions

In this work, we considered a weakly dissipative viscoelastic equation with variable-exponent nonlinearity. We showed that the decay rate of the energy is weaker than that of the relaxation function. An open question is: can we obtain a similar or even weaker decay rate in the absences of the damping term Δ u t ?

Author Contributions

Conceptualization, M.M.A.-G. and A.M.A.-M.; methodology, M.M.A.-G. and A.M.A.-M.; software, M.N.; validation, M.M.A.-G., A.M.A.-M. and J.D.A.; formal analysis, M.M.A.-G. and A.M.A.-M.; investigation, M.M.A.-G. and A.M.A.-M.; data curation, M.N.; writing—original draft preparation, J.D.A.; writing—review and editing, M.M.A.-G. and A.M.A.-M.; visualization, A.M.A.-M.; supervision, M.M.A.-G.; project administration, A.M.A.-M.; funding acquisition, A.M.A.-M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by KFUPM grant number SB20101.

Acknowledgments

The authors would like to express their profound gratitude to King Fahd University of Petroleum and Minerals (KFUPM) for its continuous support. The authors also thank the referee for his/her very careful reading and valuable comments. This work was funded by KFUPM under Project #SB201012.

Conflicts of Interest

The authors declare that there is no conflict of interest.

References

  1. Rivera, J.E.M.; Naso, M.G.; Vegni, F.M. Asymptotic behavior of the energy for a class of weakly dissipative second-order systems with memory. J. Math. Anal. Appl. 2003, 286, 692–704. [Google Scholar] [CrossRef] [Green Version]
  2. Hassan, J.H.; Messaoudi, S.A. General decay rate for a class of weakly dissipative second-order systems with memory. Math. Methods Appl. Sci. 2019, 42, 2842–2853. [Google Scholar] [CrossRef]
  3. Anaya, K.; Messaoudi, S.A. General decay rate of a weakly dissipative viscoelastic equation with a general damping. Opuscula Mathematica 2020, 40, 647–666. [Google Scholar] [CrossRef]
  4. Acerbi, E.; Mingione, G. Regularity results for stationary electro-rheological fluids. Arch. Ration. Mech. Anal. 2002, 164, 213–259. [Google Scholar] [CrossRef]
  5. Ruzicka, M. Electrorheological Fluids: Modeling and Mathematical Theory; Springer Science & Business Media: New York, NY, USA, 2000. [Google Scholar]
  6. Antontsev, S. Wave equation with p(x,t)-Laplacian and damping term: Existence and blow-up. Differ. Equ. Appl. 2011, 3, 503–525. [Google Scholar] [CrossRef] [Green Version]
  7. Antontsev, S. Wave equation with p(x,t)-Laplacian and damping term: Blow-up of solutions. Comptes Rendus Mécanique 2011, 339, 751–755. [Google Scholar] [CrossRef]
  8. Guo, B.; Gao, W. Blow-up of solutions to quasilinear hyperbolic equations with p(x,t)-Laplacian and positive initial energy. Comptes Rendus Mécanique 2014, 342, 513–519. [Google Scholar] [CrossRef]
  9. Messaoudi, S.A.; Talahmeh, A.A. A blow-up result for a nonlinear wave equation with variable-exponent nonlinearities. Appl. Anal. 2017, 96, 1509–1515. [Google Scholar] [CrossRef]
  10. Messaoudi, S.A.; Al-Smail, J.H.; Talahmeh, A.A. Decay for solutions of a nonlinear damped wave equation with variable-exponent nonlinearities. Comput. Math. Appl. 2018, 76, 1863–1875. [Google Scholar] [CrossRef]
  11. Messaoudi, S.A. On the decay of solutions of a damped quasilinear wave equation with variable-exponent nonlinearities. Math. Methods Appl. Sci. 2020, 43, 5114–5126. [Google Scholar] [CrossRef]
  12. Antontsev, S.; Shmarev, S. Blow-up of solutions to parabolic equations with nonstandard growth conditions. J. Comput. Appl. Math. 2010, 234, 2633–2645. [Google Scholar] [CrossRef] [Green Version]
  13. Antontsev, S.; Ferreira, J.; Piskin, E. Existence and blow up of solutions for a strongly damped Petrovsky equation with variable-exponent nonlinearities. Electron. J. Differ. Equ. 2021, 2021, 1–18. [Google Scholar]
  14. Abita, R. Existence and asymptotic behavior of solutions for degenerate nonlinear Kirchhoff strings with variable-exponent nonlinearities. Acta Math. Vietnam. 2021, 46, 613–643. [Google Scholar] [CrossRef]
  15. Rahmoune, A.; Benabderrahmane, B. On the viscoelastic equation with Balakrishnan–Taylor damping and nonlinear boundary/interior sources with variable-exponent nonlinearities. Stud. Univ. Babes-Bolyai Math. 2020, 65, 599–639. [Google Scholar] [CrossRef]
  16. Al-Gharabli, M.M.; Al-Mahdi, A.M.; Kafini, M. Global existence and new decay results of a viscoelastic wave equation with variable exponent and logarithmic nonlinearities. AIMS Math. 2021, 6, 10105–10129. [Google Scholar] [CrossRef]
  17. Antontsev, S.; Shmarev, S. Evolution PDEs with Nonstandard Growth Conditions; Atlantis Press: Paris, France, 2015. [Google Scholar]
  18. Diening, L.; Harjulehto, P.; Hästö, P.; Ruzicka, M. Lebesgue and Sobolev Spaces with Variable Exponents; Springer: New York, NY, USA, 2011. [Google Scholar]
  19. Radulescu, V.D.; Repovs, D.D. Partial Differential Equations with Variable Exponents: Variational Methods and Qualitative Analysis; CRC Press: Boca Raton, FL, USA, 2015; Volume 9. [Google Scholar]
  20. Mustafa, M.I. Optimal decay rates for the viscoelastic wave equation. Math. Methods Appl. Sci. 2018, 41, 192–204. [Google Scholar] [CrossRef]
  21. Jin, K.-P.; Liang, J.; Xiao, T.-J. Coupled second order evolution equations with fading memory: Optimal energy decay rate. J. Differ. Equ. 2014, 257, 1501–1528. [Google Scholar] [CrossRef]
  22. Messaoudi, S.A. General decay of solutions of a viscoelastic equation. J. Math. Anal. Appl. 2008, 341, 1457–1467. [Google Scholar] [CrossRef]
  23. Arnol’d, V.I. Mathematical Methods of Classical Mechanics; Springer Science & Business Media: New York, NY, USA, 2013; Volume 60. [Google Scholar]
Figure 1. Test 1: The solution u ( t ) of the problem at fixed values of x.
Figure 1. Test 1: The solution u ( t ) of the problem at fixed values of x.
Mca 28 00005 g001
Figure 2. Test 2: The solution u ( t ) of the problem at fixed values of x.
Figure 2. Test 2: The solution u ( t ) of the problem at fixed values of x.
Mca 28 00005 g002
Figure 3. Test 1: The energy decay.
Figure 3. Test 1: The energy decay.
Mca 28 00005 g003
Figure 4. Test 2: The energy decay.
Figure 4. Test 2: The energy decay.
Mca 28 00005 g004
Figure 5. Test 1: The solution function u ( x , t ) .
Figure 5. Test 1: The solution function u ( x , t ) .
Mca 28 00005 g005
Figure 6. Test 2: The solution function u ( x , t ) .
Figure 6. Test 2: The solution function u ( x , t ) .
Mca 28 00005 g006
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-Mahdi, A.M.; Al-Gharabli, M.M.; Noor, M.; Audu, J.D. Stability Results for a Weakly Dissipative Viscoelastic Equation with Variable-Exponent Nonlinearity: Theory and Numerics. Math. Comput. Appl. 2023, 28, 5. https://doi.org/10.3390/mca28010005

AMA Style

Al-Mahdi AM, Al-Gharabli MM, Noor M, Audu JD. Stability Results for a Weakly Dissipative Viscoelastic Equation with Variable-Exponent Nonlinearity: Theory and Numerics. Mathematical and Computational Applications. 2023; 28(1):5. https://doi.org/10.3390/mca28010005

Chicago/Turabian Style

Al-Mahdi, Adel M., Mohammad M. Al-Gharabli, Maher Noor, and Johnson D. Audu. 2023. "Stability Results for a Weakly Dissipative Viscoelastic Equation with Variable-Exponent Nonlinearity: Theory and Numerics" Mathematical and Computational Applications 28, no. 1: 5. https://doi.org/10.3390/mca28010005

APA Style

Al-Mahdi, A. M., Al-Gharabli, M. M., Noor, M., & Audu, J. D. (2023). Stability Results for a Weakly Dissipative Viscoelastic Equation with Variable-Exponent Nonlinearity: Theory and Numerics. Mathematical and Computational Applications, 28(1), 5. https://doi.org/10.3390/mca28010005

Article Metrics

Back to TopTop