Next Article in Journal
Novel Aspects of Cilia-Driven Flow of Viscoelastic Fluid through a Non-Darcy Medium under the Influence of an Induced Magnetic Field and Heat Transfer
Previous Article in Journal
Nonparametric Threshold Estimation for Drift Function in Jump–Diffusion Model of Interest Rate Using Asymmetric Kernel
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Perturbed Cauchy Viscoelastic Problem in an Exterior Domain

1
Department of Mathematics, College of Science, King Saud University, P.O. Box 2455, Riyadh 11451, Saudi Arabia
2
Department of Mathematics and Computer Science, University of Palermo, Via Archirafi 34, 90123 Palermo, Italy
*
Author to whom correspondence should be addressed.
Mathematics 2023, 11(10), 2283; https://doi.org/10.3390/math11102283
Submission received: 6 April 2023 / Revised: 11 May 2023 / Accepted: 11 May 2023 / Published: 14 May 2023
(This article belongs to the Section Difference and Differential Equations)

Abstract

:
A Cauchy viscoelastic problem perturbed by an inverse-square potential, and posed in an exterior domain of R N , is considered under a Dirichlet boundary condition. Using nonlinear capacity estimates specifically adapted to the non-local nature of the problem, the potential function and the boundary condition, we establish sufficient conditions for the nonexistence of weak solutions.
MSC:
46F25; 46G05; 46N20; 35A01

1. Introduction

Let B 1 be the open unit ball of R N ( N 2 ). In this paper, we study the following perturbed viscoelastic equation:
u t t L λ u + δ Γ ( α ) 0 t ( t s ) α 1 L λ u ( s , x ) d s + u t = | u | p , t > 0 , x R N \ B 1 .
In this problem, the elliptic operator L λ = Δ λ | x | 2 for λ R , is the Laplacian with Hardy perturbation, u = u ( t , x ) , Γ ( · ) is the gamma function, and we assume the following hypotheses on the data
λ N 2 2 2 , δ 0 , α > 0 , p > 1 .
We are interested in the range of values of N, λ (see also the parameter λ N in Equation (8)) and p for which Equation (1) does not admit weak solutions, under the Dirichlet boundary condition
u ( t , x ) = w ( x ) , t > 0 , x B 1 ,
where w L 1 ( B 1 ) , and imposing the initial conditions
u ( 0 , x ) = u 0 ( x ) , u t ( 0 , x ) = u 1 ( x ) , x R N \ B 1 .
We shall discuss separately the cases w = 0 and w 0 . We now briefly recall some known results related to perturbed viscoelastic Equation (1). When δ = λ = 0 , Equation (1) reduces to the semilinear damped wave equation
u t t Δ u + u t = | u | p , t > 0 , x R N \ B 1 .
Ogawa-Takeda [1] investigated Equation (4) under the boundary condition (2) with w = 0 . Hence, for compactly supported initial data, they showed that there is a non-negative global solution whenever 1 < p < 1 + 2 N . The approach used in [1] is based on the Kaplan–Fujita method (see the works of Kaplan [2] and Fujita [3] for more information). Later, Fino-Ibrahim-Wehbe [4] proved that the value p = 1 + 2 N belongs to the blow-up case. In a recent paper, Jleli-Samet [5] considered Equation (4) under the boundary condition (2) in the case when w is a non-negative nontrivial function. Hence, they obtained the following results:
(i)
If N = 2 , then for all p > 1 , Equation (4) under the boundary condition (2) admits no global weak solution;
(ii)
If N 3 , then for all 1 < p < 1 + 2 N 2 , Equation (4) under the boundary condition (2) admits no global weak solution;
(iii)
If N 3 and p > 1 + 2 N 2 , then Equation (4) under the boundary condition (2) admits suitable solutions for some w > 0 .
When dealing with problem (4) posed in the whole space R N , we mention the works of Kirane-Qafsaoui [6] (m-iterated Laplacian equation, m 1 ), Todorova-Yordanov [7] and Zhang [8] (global existence, blow-up and asymptotic behavior of global solutions); see also the references therein.
On the other hand, the issue of nonexistence and blow-up in finite time for viscoelastic wave equations of the form
u t t Δ u + 0 t g ( t s ) Δ u ( s , x ) d s + h ( u t ) = f ( u ) , t > 0 , x Ω
is present in many publications. We mention the works of Haraux-Zuazua [9] (hyperbolic problems), Kafini-Messaoudi [10,11] (nonlinear viscoelastic system and equation, respectively), and Messaoudi [12] (blow-up of solutions with negative initial energy). For instance, in [11], the authors investigated (5) in Ω = R N with h ( u t ) = u t and f ( u ) = | u | p 1 u . Namely, under a certain condition on the kernel function g, it was shown that, if 1 p < N N 2 , N 3 ; or p 1 , N = 1 , 2 , and
0 t g ( s ) d s < 2 p 2 2 p 1 ,
then for any initial data ( u 0 , u 1 ) H 1 ( R N ) × L 2 ( R N ) with compact support, satisfying
E ( 0 ) = 1 2 u 1 2 2 + 1 2 u 0 2 2 1 p + 1 u 0 p + 1 p + 1 0 , R N u 0 ( x ) u 1 ( x ) d x 0 ,
the corresponding solution blows up in finite time. We point out that the approach in [11] is based on the energy method.
Now we recall some references in the literature on evolution equations and inequalities perturbed by the Hardy potential λ | x | 2 . We refer to the works of Abdellaoui-Miri-Peral-Touaoula [13] (p-Laplacian equation), Abdellaoui-Peral-Primo [14,15] (Laplacian equations), Jleli-Samet-Vetro [16] (inhomogeneous wave inequalities) and again the work of Haraux-Zuazua [9]. However, to the best of our knowledge, problems of type (1) have not been previously studied in the literature. The motivation to consider Equation (1) originates from the idea to combine the effects of viscoelastic behavior and singular Hardy potential into a single wave equation. Referring to a physical context, viscoelastic materials (i.e., polymers) exhibit both the behavior of a liquid (viscous case) and of a solid (elastic case). For instance, first a suitable tension produces some elastic deformation, then (time-dependent) viscous stress occurs, hence there are material properties leading to so-called memory effects. Now, the degree of viscoelasticity can be controlled by a parameter varying in an appropriate range (see also Chapter 7 of Mills-Jenkins-Kukureka [17]). From a mathematical perspective, the effects of memory are linked to the kernel function in the integral term of the equation (i.e., the function g in (5)), hence it is interesting to show the behavior of solutions to classes of viscoelastic wave equations, under minimal (or specific) assumptions on g. We mention, for example, the works of Cavalcanti et al. [18] (using the multiplier method together with a lemma about convergent and divergent series for establishing the uniform decay of the energy of the solution) and Wu [19] (using the perturbed energy technique for establishing the uniform decay of the energy of the solution to the system of viscoelastic wave equations). On the other hand, the singular Hardy potential is recognized as a suitable prototype to analyze the critical behavior of different nonlinear problems in physics, hence in dealing with the existence and stability of solutions (for more details and information, we refer to the comprehensive book of Alonso-de Diego [20]). It makes more sense to study how the parametric Hardy potential λ | x | 2 for λ R , affects instantaneous and complete blow-up of solutions to (1) (i.e., nonexistence phenomenon).
In order to define weak solutions to (1) under conditions (2) and (3), we recall below some notions from fractional calculus (see the comprehensive book of Kilbas-Srivastava-Trujillo [21] for more details), hence we fix notation.
Let T > 0 be fixed. Given f L 1 ( [ 0 , T ] ) and β > 0 , the left-sided and right-sided Riemann–Liouville fractional integrals of order β of f, are defined, respectively, by
I 0 β f ( t ) = 1 Γ ( β ) 0 t ( t s ) β 1 f ( s ) d s
and
I T β f ( t ) = 1 Γ ( β ) t T ( s t ) β 1 f ( s ) d s
for almost everywhere t [ 0 , T ] . The following property can be found in Kilbas-Srivastava-Trujillo ([21], Lemma 2.7).
Lemma 1. 
Let β > 0 , m , q 1 , and 1 m + 1 q 1 + β ( m 1 and q 1 if 1 m + 1 q = 1 + β ). If f L m ( [ 0 , T ] ) and g L q ( [ 0 , T ] ) , then
0 T g ( t ) I 0 β f ( t ) d t = 0 T f ( t ) I T β g ( t ) d t .
Let F : [ 0 , T ] × R N \ B 1 R be a given function. The left-sided and right-sided Riemann–Liouville fractional integrals of order β > 0 of F with respect to the time-variable t, are denoted, respectively, by I 0 β F and I T β F , namely we have
I 0 β F ( t , x ) = I 0 β F ( · , x ) ( t )
and
I T β F ( t , x ) = I T β F ( · , x ) ( t ) .
Using the above notations, the nonlocal term in Equation (1) can be written in the form
1 Γ ( α ) 0 t ( t s ) α 1 L λ u ( s , x ) d s = I 0 α L λ u ( t , x ) .
For some contributions related to the applications of fractional derivative operators in diffusion processes, we mention the works of Hurtado-Salvatierra [22], Jleli [23], Nashine et al. [24], Villagran et al. [25], and the references therein. Precisely, [22] deals with non-local diffusion problems driven by the fractional p-Laplacian differential operator in the Heisenberg group. The approach is based on the theory of monotone operators and pullback attractors. In [23], the author investigates the existence of solutions to a Dirichlet problem for the Kohn Laplacian on the Heisenberg group too, using partial ordering methods. In [24], the authors study generalized fractional integral equations, using fixed-point arguments in Banach spaces. In [25], the authors investigate stability for a system of wave equations. They establish well-posedness and polynomial stability using the semigroup theory and certain sharp results. For all T > 0 , we denote
Q T = [ 0 , T ] × R N \ B 1 , Γ T = [ 0 , T ] × B 1 ,
hence, Γ T Q T . Let Φ T be the set of functions φ = φ ( t , x ) satisfying the following properties:
(P 1 )
φ C 2 ( Q T ) , supp x ( φ ) R N \ B 1 , φ 0 ;
(P 2 )
φ ( T , · ) = φ t ( T , · ) = 0 ;
(P 3 )
φ | Γ T = 0 .
Using standard integrations by parts, together with Lemma 1 and (6), we define weak solutions to problem (1)–(3) as follows.
Definition 1. 
Let u 0 , u 1 L loc 1 R N \ B 1 and w L 1 ( B 1 ) . We say that
u L loc p ( [ 0 , ) × R N \ B 1 )
is a weak solution to (1)–(3) if
Q T | u | p φ d x d t + R N \ B 1 u 0 ( x ) φ ( 0 , x ) φ t ( 0 , x ) d x + R N \ B 1 u 1 ( x ) φ ( 0 , x ) d x Γ T φ ν δ ( I T α φ ) ν w ( x ) d S x d t = Q T u φ t t L λ φ + δ L λ ( I T α φ ) φ t d x d t
for every T > 0 and φ Φ T , where ν denotes the outward unit normal vector on B 1 , relative to R N \ B 1 .
For λ N 2 2 2 , we introduce the parameter
λ N = N 2 2 + N 2 2 2 + λ
and the truncation function
H λ ( x ) = | x | λ N 1 | x | 2 N 2 λ N if λ > N 2 2 2 , | x | λ N ln | x | if λ = N 2 2 2 .
Our aim in this work is to establish sufficient conditions for the nonexistence of weak solutions to problem (1)–(3). Therefore, we need to find ways to deal with the nonlocal nature of the problem, the elliptic operator L λ and the boundary condition (2). We come up with an approach based on nonlinear capacity estimates specifically adapted to our needs.
The rest of the paper is organized as follows. In Section 2, we obtain some preliminary estimates. Namely, we first prove an a priori estimate for problem (1)–(3), then we construct a family of test functions belonging to the functional space Φ T , and provide some useful estimates involving such functions. In Section 3, we provide the proofs of Theorems 1 and 2.

2. Preliminaries

In this section, we give the mathematical background necessary to establish our results. Here, the symbols C , C i denote always generic positive constants, which are independent of the scaling parameters T , R and the solution u. Their values could be changed from one line to another. First, we impose the following hypotheses on the data: N 2 , λ N 2 2 2 , α > 0 , δ 0 , p > 1 , u 0 , u 1 L loc 1 R N \ B 1 and w L 1 ( B 1 ) .

2.1. A Priori Estimate

For T > 0 and φ ϕ T , we introduce the integral terms
J 1 ( φ ) = supp ( φ ) φ 1 p 1 | φ t | p p 1 d x d t ,
J 2 ( φ ) = supp ( φ ) φ 1 p 1 | φ t t | p p 1 d x d t ,
J 3 ( φ ) = supp ( φ ) φ 1 p 1 | L λ φ | p p 1 d x d t ,
J 4 ( φ ) = supp ( φ ) φ 1 p 1 | L λ ( I T α φ ) | p p 1 d x d t .
Hence, we establish the following a priori estimate.
Lemma 2. 
If u L loc p ( [ 0 , ) × R N \ B 1 ) is a weak solution to problem (1)–(3), then we have the estimate
R N \ B 1 u 0 ( x ) φ ( 0 , x ) φ t ( 0 , x ) d x + R N \ B 1 u 1 ( x ) φ ( 0 , x ) d x Γ T φ ν δ ( I T α φ ) ν w ( x ) d S x d t C i = 1 4 J i ( φ )
for every T > 0 and φ Φ T , provided that J i ( φ ) < , j = 1 , 2 , 3 , 4 .
Proof. 
Let u L loc p ( [ 0 , ) × R N \ B 1 ) be a weak solution to problem (1)–(3). For T > 0 and φ Φ T satisfying the conditions J i ( φ ) < , j = 1 , 2 , 3 , 4 , using (7) (i.e., Definition 1), we obtain
Q T | u | p φ d x d t + R N \ B 1 u 0 ( x ) φ ( 0 , x ) φ t ( 0 , x ) d x + R N \ B 1 u 1 ( x ) φ ( 0 , x ) d x Γ T φ ν δ ( I T α φ ) ν w ( x ) d S x d t Q T | u | | φ t | d x d t + Q T | u | | φ t t | d x d t + Q T | u | | L λ φ | d x d t + | δ | Q T | u | | L λ ( I T α φ ) | d x d t .
Considering each one of the integrals on the right hand side of (15) separately, by means of Young’s inequality, we first obtain
Q T | u | | φ t | d x d t = Q T | u | φ 1 p φ 1 p | φ t | d x d t 1 4 Q T | u | p φ d x d t + C J 1 ( φ ) ,
and similarly we obtain the other bounds
Q T | u | | φ t t | d x d t 1 4 Q T | u | p φ d x d t + C J 2 ( φ ) ,
Q T | u | | L λ φ | d x d t 1 4 Q T | u | p φ d x d t + C J 3 ( φ ) ,
| δ | Q T | u | | L λ ( I T α φ ) | d x d t 1 4 Q T | u | p φ d x d t + C J 4 ( φ ) .
After combining together inequalities (16)–(19) with the principal inequality (15), we obtain the desired estimate in (14). □

2.2. Construction of a Family of Functions Belonging to Φ T

Let T > 0 . For sufficiently large , we introduce the function
ι T ( t ) = T ( T t ) , 0 t T .
Next, for sufficiently large R, we consider a family of cut-off functions { ξ R } R fulfilling the following properties:
0 ξ R 1 , ξ R C ( R N \ B 1 ) , ξ R = 1 i f 1 | x | R , ξ R = 0 i f | x | 2 R ,
but also
| ξ R | C R 1 , | Δ ξ R | C R 2 .
We now introduce a test function of the form
φ ( t , x ) = ι T ( t ) ψ R ( x ) , ( t , x ) Q T ,
where
ψ R ( x ) = H λ ( x ) ξ R ( x ) , x R N \ B 1
and H λ ( · ) is the truncation in (9).
We immediately find that the test function in (23) possesses the properties (P 1 )–(P 3 ), hence φ Φ T .
Lemma 3. 
For sufficiently large T , R and ℓ, the function φ defined by (23) belongs to Φ T .
Proof. 
Property (P 1 ) follows immediately from (9) and (20)–(24). Moreover, by (9) and (20), one has
ι T ( T ) = ι T ( T ) = 0 , H λ | B 1 = 0 ,
which shows by (23) that properties (P 2 ) and (P 3 ) are satisfied too. □

2.3. Estimates of J i ( φ )

Now we consider again the integral terms (10)–(13) to establish their bounds from above. For sufficiently large T , R and , let φ be the test function defined by (23). The first result holds for J i ( φ ) , with j = 1 , 2 , and is given in the following lemma.
Lemma 4. 
For i = 1 , 2 and λ N given in (8), we have
J i ( φ ) C T 1 i p p 1 R λ N + N ln R .
Proof. 
Starting from the formulas (10) and (11), and involving the cut-off function (21) and the test function (23) (recall (24) too), for i = 1 , 2 , we obtain
J i ( φ ) = 0 T ι T 1 p 1 ( t ) | ι T ( i ) ( t ) | p p 1 d t 1 < | x | < 2 R H λ ( x ) ξ R ( x ) d x ,
where ι T ( i ) = d i ι T d t i . On the other hand, by (20), for all 0 < t < T one has
| ι T ( i ) ( t ) | = C T ( T t ) i ,
which yields
ι T 1 p 1 ( t ) | ι T ( i ) ( t ) | p p 1 = C T ( T t ) i p p 1 .
Integrating over ( 0 , T ) this equation, we obtain
0 T ι T 1 p 1 ( t ) | ι T ( i ) ( t ) | p p 1 d t = C T 1 i p p 1 .
Moreover, we use the truncation (9) together with the appropriate property in (21) (namely 0 ξ R 1 ) to deduce that
H λ ( x ) ξ R ( x ) H λ ( x ) | x | λ N ln | x | , 1 < | x | < 2 R .
Integrating over ( 1 , 2 R ) , we obtain (notice that λ N + N 1 > 0 )
1 < | x | < 2 R H λ ( x ) ξ R ( x ) d x 1 < | x | < 2 R | x | λ N ln | x | ln R r = 1 2 R r λ N + N 1 d r C ln R R λ N + N 1 R = C R λ N + N ln R .
Combining (26)–(28), we conclude the estimate (25). □
The next result follows by elementary calculations, hence we avoid the proof of this lemma.
Lemma 5. 
The function H λ defined by (9) satisfies the following property:
L λ H λ ( x ) = 0 , x R N \ B 1 .
Now, we consider the integral term J 3 ( φ ) and establish the following estimate.
Lemma 6. 
The following estimate holds:
J 3 ( φ ) C T R λ N + N 2 p p 1 ( ln R ) p p 1 .
Proof. 
Starting from the formula (12), and involving the cut-off function (21) and the test function (23), we obtain
J 3 ( φ ) = 0 T ι T ( t ) d t 1 < | x | < 2 R ψ R 1 p 1 | L λ ψ R | p p 1 d x .
On the other hand, by (20), one has
0 T ι T ( t ) d t = T 0 T ( T t ) d t = C T .
Next, by (24), for 1 < | x | < 2 R , we obtain
L λ ψ R ( x ) = L λ ( H λ ξ R ) ( x ) = Δ ( H λ ( x ) ξ R ( x ) ) λ | x | 2 H λ ( x ) ξ R ( x ) = ξ R ( x ) Δ H λ ( x ) + H λ ( x ) Δ ( ξ R ( x ) ) + 2 H λ ( x ) · ( ξ R ( x ) ) λ | x | 2 H λ ( x ) ξ R ( x ) = ξ R ( x ) L λ H λ ( x ) + H λ ( x ) Δ ( ξ R ( x ) ) + 2 H λ ( x ) · ( ξ R ( x ) ) ,
where · denotes the inner product in R N . Then, by Lemma 5, we deduce that
L λ ψ R ( x ) = H λ ( x ) Δ ( ξ R ( x ) ) + 2 H λ ( x ) · ( ξ R ( x ) ) ,
which implies, by (21), that
1 < | x | < 2 R ψ R 1 p 1 | L λ ψ R | p p 1 d x = R < | x | < 2 R ψ R 1 p 1 | L λ ψ R | p p 1 d x .
On the other hand, by using the truncation (9), we obtain
C 1 R λ N H λ ( x ) C 2 R λ N ln R , | H λ ( x ) | C R λ N 1 ln R , R < | x | < 2 R .
Moreover, using the properties (22) together with
Δ ξ R = ( 1 ) ξ R 2 | ξ R | 2 + ξ R 1 Δ ξ R ,
we deduce that (recall 0 ξ R 1 )
| Δ ( ξ R ( x ) ) | C R 2 ξ R 2 ( x ) , | ( ξ R ( x ) ) | C R 1 ξ R 2 ( x ) , R < | x | < 2 R .
Hence, by Cauchy–Schwarz inequality, it follows from (32), (34) and (35) that
| L λ ψ R ( x ) | C R λ N 2 ln R ξ R 2 ( x ) , R < | x | < 2 R .
Furthermore, using (24), (34) and (36), we obtain
ψ R 1 p 1 ( x ) | L λ ψ R ( x ) | p p 1 C R λ N 2 p p 1 ( ln R ) p p 1 ξ R 2 p p 1 ( x ) , R < | x | < 2 R .
Integrating appropriately this inequality and using (33), we obtain
1 < | x | < 2 R ψ R 1 p 1 | L λ ψ R | p p 1 d x C R λ N 2 p p 1 ( ln R ) p p 1 R < | x | < 2 R ξ R 2 p p 1 ( x ) d x C R N + λ N 2 p p 1 ( ln R ) p p 1 .
Finally, (29) follows from (30), (31) and (37). Hence, the estimate of J 3 ( φ ) is reached. □
Now, we consider the integral term J 4 ( φ ) and prove the following result.
Lemma 7. 
The following estimate holds:
J 4 ( φ ) C T 1 + α p p 1 R λ N + N 2 p p 1 ( ln R ) p p 1 .
Proof. 
Starting from the formula (13), and using the test function (23), we obtain
J 4 ( φ ) = 0 T ι T 1 p 1 ( t ) | I T α ι T ( t ) | p p 1 d t 1 < | x | < 2 R ψ R 1 p 1 | L λ ψ R | p p 1 d x .
On the other hand, by (20), for all 0 < t < T , one has
I T α ι T ( t ) = 1 Γ ( α ) t T ( s t ) α 1 ι T ( s ) d s = T Γ ( α ) t T ( s t ) α 1 ( T s ) d s = T Γ ( α ) t T ( T t ) ( T s ) α 1 ( T s ) d s = T ( T t ) α 1 Γ ( α ) t T 1 T s T t α 1 ( T s ) d s .
Then, by the change of variable ϑ = T s T t , we obtain
I T α ι T ( t ) = T ( T t ) σ + Γ ( α ) 0 1 ( 1 ϑ ) α 1 ϑ ( + 1 ) 1 d ϑ = T ( T t ) α + Γ ( α ) B ( α , + 1 ) ,
where B ( · , · ) is the Beta function. Hence, there holds
ι T 1 p 1 ( t ) | I T α ι T ( t ) | p p 1 = C T ( T t ) + α p p 1 .
Integrating this equation over ( 0 , T ) , we obtain
0 T ι T 1 p 1 ( t ) | I T α ι T ( t ) | p p 1 d t = C T 1 + α p p 1 .
Combining (37), (39) and (41), we conclude the estimate (38). □

3. Main Results

Our main results are stated in the following theorems. As already mentioned, we first consider the case w = 0 .
Theorem 1. 
Let N 2 , λ N 2 2 2 , α > 0 , δ 0 , w = 0 and u 0 , u 1 L loc 1 R N \ B 1 . If
( u 0 + u 1 ) H λ L 1 R N \ B 1 , R N \ B 1 ( u 0 ( x ) + u 1 ( x ) ) H λ ( x ) d x > 0 ,
and
1 < p < p * ( α , λ , N ) : = 1 + 2 ( α + 1 ) ( λ N + N ) ,
then problem (1)–(3) admits no weak solution.
Proof. 
We use the contradiction argument. Namely, let us suppose that u L loc p ( [ 0 , ) × R N \ B 1 ) is a weak solution to problem (1)–(3). Then, by Lemma 2 (with w = 0 ) and Lemma 3, for sufficiently large T , R and , we have
R N \ B 1 u 0 ( x ) φ ( 0 , x ) φ t ( 0 , x ) + u 1 ( x ) φ ( 0 , x ) d x C i = 1 4 J i ( φ ) ,
where φ is the function defined by (23). On the other hand, by (23), for all x R N \ B 1 , one has
φ ( 0 , x ) = ι T ( 0 ) ψ R ( x ) = ψ R ( x ) = H λ ( x ) ξ R ( x ) ,
as well as
φ t ( 0 , x ) = ι T ( 0 ) ψ R ( x ) = T 1 H λ ( x ) ξ R ( x ) .
Hence, we deduce that
R N \ B 1 u 0 ( x ) φ ( 0 , x ) φ t ( 0 , x ) + u 1 ( x ) φ ( 0 , x ) d x = R N \ B 1 u 0 ( x ) 1 + T 1 + u 1 ( x ) H λ ( x ) ξ R ( x ) d x .
Using Lemmas 4, 6 and 7, together with (44) and (45), we obtain
R N \ B 1 u 0 ( x ) 1 + T 1 + u 1 ( x ) H λ ( x ) ξ R ( x ) d x C T 1 p p 1 R λ N + N ln R + T 1 + α p p 1 R λ N + N 2 p p 1 ( ln R ) p p 1 .
Taking T = R θ , where
θ = 2 α + 1 ,
one obtains
T 1 p p 1 R λ N + N = T 1 + α p p 1 R λ N + N 2 p p 1 = R κ ,
with
κ = λ N + N 2 ( α + 1 ) ( p 1 ) .
Consequently, the inequality (46) reduces to the following one
R N \ B 1 u 0 ( x ) 1 + R θ + u 1 ( x ) H λ ( x ) ξ R ( x ) d x C R κ ln R + ( ln R ) p p 1 .
The hypotheses (42) and properties (21) of cut-off functions, together with the dominated convergence theorem, lead to
lim R R N \ B 1 u 0 ( x ) 1 + R θ + u 1 ( x ) H λ ( x ) ξ R ( x ) d x = R N \ B 1 ( u 0 ( x ) + u 1 ( x ) ) H λ ( x ) d x > 0 .
Moreover, the hypothesis (43) gives us
κ < 0 .
Next, passing to the limit as R in (48), by the limit in (49) and the sign condition (50), we obtain the following contradiction:
0 < R N \ B 1 ( u 0 ( x ) + u 1 ( x ) ) H λ ( x ) d x 0 .
It follows that problems (1)–(3) admit no weak solutions. This completes the proof of Theorem 1. □
Following similar arguments of proof, in the inhomogeneous case w 0 , we conclude the following theorem.
Theorem 2. 
Let N 2 , λ N 2 2 2 , α > 0 , δ 0 , w L 1 ( B 1 ) and u 0 , u 1 L loc 1 R N \ B 1 , u 0 , u 1 0 . Assuming
δ B 1 w ( x ) d S x < 0 ,
we have the following:
(I) 
If λ N + N 2 = 0 , then for all p > 1 , problems (1)–(3) admit no weak solutions.
(II) 
If λ N + N 2 0 , then for all
1 < p < p * * ( α , λ , N ) : = 1 + 2 ( α + 1 ) ( λ N + N 2 )
problems (1)–(3) admit no weak solutions.
Proof. 
We use the contradiction argument. Let us suppose that u L loc p ( [ 0 , ) × R N \ B 1 ) is a weak solution to the problem (1)–(3). Then, by Lemmas 2 and 3, for sufficiently large T , R and , we have
R N \ B 1 u 0 ( x ) φ ( 0 , x ) φ t ( 0 , x ) d x + R N \ B 1 u 1 ( x ) φ ( 0 , x ) d x Γ T φ ν δ ( I T α φ ) ν w ( x ) d S x d t C i = 1 4 J i ( φ ) ,
where φ is the function defined by (23). On the other hand, since u 0 , u 1 0 , it follows from (45) that
R N \ B 1 u 0 ( x ) φ ( 0 , x ) φ t ( 0 , x ) d x + R N \ B 1 u 1 ( x ) φ ( 0 , x ) d x 0 .
Then, we have
Γ T φ ν δ ( I T α φ ) ν w ( x ) d S x d t C i = 1 4 J i ( φ ) .
Further, by using (21) and (23), we obtain
Γ T φ ν δ ( I T α φ ) ν w ( x ) d S x d t = 0 T I T α ι T ( t ) 1 δ ι T ( t ) d t B 1 H λ ν ( x ) δ w ( x ) d S x .
Involving (20) and (40), for all 0 < t < T , we deduce that
I T α ι T ( t ) 1 δ ι T ( t ) = B ( α , + 1 ) Γ ( α ) T ( T t ) α + T δ ( T t ) .
Integrating this equation over ( 0 , T ) , we obtain
0 T I T α ι T ( t ) 1 δ ι T ( t ) d t = T α + 1 B ( α , + 1 ) Γ ( α ) ( α + + 1 ) 1 δ ( + 1 ) T α .
Since
lim T B ( α , + 1 ) Γ ( α ) ( α + + 1 ) 1 δ ( + 1 ) T α = B ( α , + 1 ) Γ ( α ) ( α + + 1 ) > 0 ,
then, for sufficiently large T, one has
0 T I T α ι T ( t ) 1 δ ι T ( t ) d t C T α + 1 .
On the other hand, with respect to the truncation function (9), for all x B 1 , we obtain
H λ ν ( x ) = 2 N 2 λ N < 0 i f λ > N 2 2 2 , 1 i f λ = N 2 2 2 .
This shows that
B 1 H λ ν ( x ) δ w ( x ) d S x = C δ B 1 w ( x ) d S x .
By Lemmas 4, 6 and 7, together with formulas (53)–(56) and hypothesis (51), we obtain the following inequality:
T α + 1 δ B 1 w ( x ) d S x C T 1 p p 1 R λ N + N ln R + T 1 + α p p 1 R λ N + N 2 p p 1 ( ln R ) p p 1 ,
that is,
δ B 1 w ( x ) d S x C T α p p 1 R λ N + N ln R + T α p 1 R λ N + N 2 p p 1 ( ln R ) p p 1 .
Taking T = R θ , where θ is given by (47), then (57) reduces to
δ B 1 w ( x ) d S x C R μ ln R + ( ln R ) p p 1 ,
where
μ = ( α + 1 ) ( λ N + N 2 ) p ( α + 1 ) ( λ N + N 2 ) + 2 ( α + 1 ) ( p 1 ) .
Observe that in the case (I), that is λ N + N 2 = 0 , one has μ < 0 . Similarly, in the case (II), that is λ N + N 2 0 and p satisfies (52), we have μ < 0 . Hence, passing to the limit as R in (58), we obtain
δ B 1 w ( x ) d S x 0 ,
which contradicts hypothesis (51). We conclude that problem (1)–(3) admits no weak solution. This completes the proof of Theorem 2. □

4. Conclusions

In this paper, we have obtained the nonexistence of weak solutions for the viscoelastic equation (1) in the presence of both homogeneous and inhomogeneous Dirichlet boundary conditions. Then, we have constructed the proofs over corresponding estimates of integral terms in the definition of weak solutions to (1). These estimates help in the analysis of the behavior of solutions to viscoelastic equation (1) in comparison with classical (damped) wave equations of physical interest. For instance, we point out the following two facts:
(i)
Referring to Theorem 1, in the limit case α 0 + we note that
p * ( α , λ = 0 , N ) 1 + 2 N ,
which is the critical exponent for Equation (4) under the boundary condition (2) with w = 0 .
(ii)
Referring to Theorem 2, in the limit case α 0 + we note that
p * * ( α , λ = 0 , N ) 1 + 2 N 2 ( N 3 ) ,
which is the critical exponent for Equation (4) under the boundary condition (2) with w 0 .
Regularization and general decay of energy of solutions for different viscoelastic equations are interesting topics that can be further studied under different data and boundary conditions. We mention the work of Han-Wang [26] (positive decaying kernel function in the memory term) and Thanh Binh et al. [27] (strongly damped wave equation involving statistical discrete data).

Author Contributions

Investigation, B.S. and C.V.; Writing—original draft, B.S. and C.V. All authors have read and agreed to the published version of the manuscript.

Funding

The first author is supported by Researchers Supporting Project number (RSP2023R4), King Saud University, Riyadh, Saudi Arabia. The second author is supported by the research fund of University of Palermo: “FFR 2022 Calogero Vetro”.

Data Availability Statement

Data sharing not applicable to this paper as no data sets were generated or analyzed during the current study.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ogawa, T.; Takeda, H. Non-existence of weak solutions to nonlinear damped wave equations in exterior domains. Nonlinear Anal. 2009, 70, 3696–3701. [Google Scholar] [CrossRef]
  2. Kaplan, S. On the growth of solutions of quasi-linear parabolic equations. Commun. Pure Appl. Math. 1963, 16, 305–330. [Google Scholar] [CrossRef]
  3. Fujita, H. On the blowing up of solutions of the Cauchy problem for utu+u1+α. J. Sci. Univ. Tokyo Sect. I 1996, 13, 109–124. [Google Scholar]
  4. Fino, A.Z.; Ibrahim, H.; Wehbe, A. A blow-up result for a nonlinear damped wave equation in exterior domain: The critical case. Comput. Math. Appl. 2017, 73, 2415–2420. [Google Scholar] [CrossRef]
  5. Jleli, M.; Samet, B. New blow-up results for nonlinear boundary value problems in exterior domains. Nonlinear Anal. 2019, 178, 348–365. [Google Scholar] [CrossRef]
  6. Kirane, M.; Qafsaoui, M. Fujita’s exponent for a semilinear wave equation with linear damping. Adv. Nonlinear Stud. 2002, 2, 41–49. [Google Scholar] [CrossRef]
  7. Todorova, G.; Yordanov, B. Critical exponent for a nonlinear wave equation with damping. J. Differ. Equ. 2001, 174, 464–489. [Google Scholar] [CrossRef]
  8. Zhang, Q.S. A blow-up result for a nonlinear wave equation with damping: The critical case. C. R. Acad. Sci. Paris Ser. I 2001, 333, 109–114. [Google Scholar] [CrossRef]
  9. Haraux, A.; Zuazua, E. Decay estimates for some semilinear damped hyperbolic problems. Arch. Ration. Mech. Anal. 1988, 150, 191–206. [Google Scholar] [CrossRef]
  10. Kafini, M.; Messaoudi, S.A. A blow up result for a viscoelastic system in RN. Electron. J. Differ. Equ. 2006, 2006, 1–7. [Google Scholar]
  11. Kafini, M.; Messaoudi, S.A. A blow up result in a Cauchy viscoelastic problem. Appl. Math. Lett. 2008, 21, 549–553. [Google Scholar] [CrossRef]
  12. Messaoudi, S.A. Blow up and global existence in nonlinear viscoelastic wave equations. Math. Nachrichten 2003, 260, 58–66. [Google Scholar] [CrossRef]
  13. Abdellaoui, B.; Miri, S.E.; Peral, I.; Touaoula, T.M. Some remarks on quasilinear parabolic problems with singular potential and a reaction term. Nonlinear Differ. Equ. Appl. 2014, 21, 453–490. [Google Scholar] [CrossRef]
  14. Abdellaoui, B.; Peral, I.; Primo, A. Influence of the Hardy potential in a semi-linear heat equation. Proc. R. Soc. Edinburgh. Sect. A 2009, 139, 897–926. [Google Scholar] [CrossRef]
  15. Abdellaoui, B.; Peral, I.; Primo, A. Strong regularizing effect of a gradient term in the heat equation with the Hardy potential. J. Funct. Anal. 2010, 258, 1247–1272. [Google Scholar] [CrossRef]
  16. Jleli, M.; Samet, B.; Vetro, C. On the critical behavior for inhomogeneous wave inequalities with Hardy potential in an exterior domain. Adv. Nonlinear Anal. 2021, 10, 1267–1283. [Google Scholar] [CrossRef]
  17. Mills, N.; Jenkins, M.; Kukureka, S. Plastics: Microstructure and Engineering Applications, 4th ed.; Butterworth-Heinemann: Oxford, UK, 2020. [Google Scholar]
  18. Cavalcanti, M.M.; Domingos Cavalcanti, V.N.; Martinez, P. General decay rate estimates for viscoelastic dissipative systems. Nonlinear Anal. 2008, 68, 177–193. [Google Scholar] [CrossRef]
  19. Wu, S.-T. General decay of solutions for a nonlinear system of viscoelastic wave equations with degenerate damping and source terms. J. Math. Anal. Appl. 2013, 406, 34–48. [Google Scholar] [CrossRef]
  20. Alonso, I.P.; de Diego, F.S. Elliptic and Parabolic Equations Involving the Hardy–Leray Potential; Walter de Gruyter GmbH & Co. KG: Berlin, Germany, 2021. [Google Scholar]
  21. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations. In North-Holland Mathematics Studies; Mill, J.V., Ed.; Elsevier: Amsterdam, The Netherlands, 2006; Volume 204. [Google Scholar]
  22. Hurtado, E.J.; Salvatierra, A.P. A stability result of a fractional heat equation and time fractional diffusion equations governed by fractional fluxes in the Heisenberg group. Rend. Circ. Mat. Palermo II. Ser. 2023. [Google Scholar] [CrossRef]
  23. Jleli, M. On the existence of nonnegative radial solutions for Dirichlet exterior problems on the Heisenberg group. Demonstr. Math. 2023, 56, 20220193. [Google Scholar] [CrossRef]
  24. Nashine, H.K.; Ibrahim, R.W.; Agarwal, R.P.; Can, N.H. Existence of local fractional integral equation via a measure of non-compactness with monotone property on Banach spaces. Adv. Differ. Equ. 2020, 2020, 697. [Google Scholar] [CrossRef]
  25. Villagran, O.P.V.; Nonato, C.A.; Raposo, C.A.; Ramos, A.J.A. Stability for a weakly coupled wave equations with a boundary dissipation of fractional derivative type. Rend. Circ. Mat. Palermo II. Ser. 2023, 72, 803–831. [Google Scholar] [CrossRef]
  26. Han, X.; Wang, M. General decay of energy for a viscoelastic equation with nonlinear damping. Math. Methods Appl. Sci. 2009, 32, 346–358. [Google Scholar] [CrossRef]
  27. Thanh Binh, T.; Huu Can, N.; Quoc Nam, D.H.; Ngoc Thach, T. Regularization of a two-dimensional strongly damped wave equation with statistical discrete data. Math. Methods Appl. Sci. 2020, 43, 4317–4335. [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

Samet, B.; Vetro, C. A Perturbed Cauchy Viscoelastic Problem in an Exterior Domain. Mathematics 2023, 11, 2283. https://doi.org/10.3390/math11102283

AMA Style

Samet B, Vetro C. A Perturbed Cauchy Viscoelastic Problem in an Exterior Domain. Mathematics. 2023; 11(10):2283. https://doi.org/10.3390/math11102283

Chicago/Turabian Style

Samet, Bessem, and Calogero Vetro. 2023. "A Perturbed Cauchy Viscoelastic Problem in an Exterior Domain" Mathematics 11, no. 10: 2283. https://doi.org/10.3390/math11102283

APA Style

Samet, B., & Vetro, C. (2023). A Perturbed Cauchy Viscoelastic Problem in an Exterior Domain. Mathematics, 11(10), 2283. https://doi.org/10.3390/math11102283

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop