Next Article in Journal
An Efficient Subspace Minimization Conjugate Gradient Method for Solving Nonlinear Monotone Equations with Convex Constraints
Previous Article in Journal
Theory of Functions and Applications
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Class of Fractional Viscoelastic Kirchhoff Equations Involving Two Nonlinear Source Terms of Different Signs

College of Mathematics and Computer Science, Northwest Minzu University, Lanzhou 730030, China
*
Author to whom correspondence should be addressed.
Axioms 2024, 13(3), 169; https://doi.org/10.3390/axioms13030169
Submission received: 21 January 2024 / Revised: 25 February 2024 / Accepted: 2 March 2024 / Published: 5 March 2024
(This article belongs to the Section Mathematical Analysis)

Abstract

:
A class of fractional viscoelastic Kirchhoff equations involving two nonlinear source terms of different signs are studied. Under suitable assumptions on the exponents of nonlinear source terms and the memory kernel, the existence of global solutions in an appropriate functional space is established by a combination of the theory of potential wells and the Galerkin approximations. Furthermore, the asymptotic behavior of global solutions is obtained by a combination of the theory of potential wells and the perturbed energy method.

1. Introduction

In this paper, we deal with the initial boundary value problem for a class of fractional viscoelastic Kirchhoff equations involving two nonlinear source terms of different signs.
u t t + ( a + b [ u ] m 2 p 2 ) ( Δ ) m u 0 t g ( t τ ) ( Δ ) m u ( τ ) d τ + u t = | u | q 2 u | u | r 2 u , x Ω , t > 0 ,
u ( x , 0 ) = u 0 ( x ) , u t ( x , 0 ) = u 1 ( x ) , x Ω ,
u ( x , t ) = 0 , x R N \ Ω , t > 0 ,
where [ u ] m is the Gagliardo seminorm defined by
[ u ] m : = R 2 N | u ( x , t ) u ( y , t ) | 2 | x y | N + 2 m d x d y 1 2 ,
[ u ] m 2 p 2 means that 2 p 2 is a power of [ u ] m , and ( Δ ) m is the fractional Laplacian with 0 < m < 1 , which (up to normalization factors) could be given by
( Δ ) m φ ( x ) : = R N 2 φ ( x ) φ ( x + y ) φ ( x y ) | y | N + 2 m d y , x R N
for all φ C 0 ( R N ) . In addition, a > 0 , b 0 , p > 1 , Ω R N ( N 1 ) is a bounded domain with a Lipschitz boundary, and the functions u 0 ( x ) and u 1 ( x ) are specified later.
The problem (1)–(3) can govern the motion of the viscoelastic string with a fractional length (see Ref. [1]). The unknown function u = u ( x , t ) represents the vertical displacement, 0 t g ( t τ ) ( Δ ) m u ( τ ) d τ is the viscoelastic term, u t is a weakly damped term, | u | q 2 u and | u | r 2 u are two nonlinear source terms of different signs, and q, r, and the memory kernel g are introduced later.
In the case where m = 1 , Equation (1) becomes a classical viscoelastic Kirchhoff wave equation. Concerning this class of equations, Torrejón and Yong [2] investigated
u t t h ( u 2 2 ) Δ u + 0 t g ( t τ ) Δ u ( τ ) d τ = f ,
where h is the Kirchhoff function. They obtained existence, uniqueness, and asymptotic behavior of global solutions. Wu and Tsai [3] treated a viscoelastic Kirchhoff wave equation with nonlinear weak damping and source terms:
u t t h ( u 2 2 ) Δ u + 0 t g ( t τ ) Δ u ( τ ) d τ + f 2 ( u t ) = f 1 ( u ) .
They derived existence and blow-up of local solutions, and also established the estimates on the blow-up time. Liu et al. [4] studied a more general viscoelastic Kirchhoff wave equation:
u t t h ( u 2 2 ) Δ u + 0 t g ( t τ ) Δ u ( τ ) d τ + f 2 ( x , t , u , u t ) = f 1 ( x , t , u ) .
They proved the nonexistence of global solutions.
Fractional partial differential equations arise in continuum mechanics [5], quantum mechanics [6,7], population dynamics [5], anomalous diffusion [8,9], fluid mechanics [10,11], and so on, and they have received considerable attention. Ambrosio and Isernia [12] investigated a fractional stationary Kirchhoff equation:
( a + b ( 1 m ) [ u ] m 2 ) ( Δ ) m u = f ( u ) ,
where a and b are two parameters. By the minimax arguments, they established the multiplicity of solutions, provided b is sufficiently small. Nyamoradi and Ambrosio [13] dealt with a fractional stationary Kirchhoff equation involving two nonlinear source terms:
( a + b [ u ] m 2 θ 2 ) ( Δ ) m u = λ f ( x ) | u | q 2 u + | u | 2 m * 2 u ,
where θ and λ are also parameters, and 2 m * represents the critical exponent of the fractional Sobolev space. Under suitable values of the parameters, they derived the existence and nonexistence of multiple solutions. Fiscella and Mishra [14] studied the following fractional stationary Kirchhoff equation involving two nonlinear source terms:
( a + b [ u ] m 2 θ 2 ) ( Δ ) m u μ u | x | 2 m = λ f 1 ( x ) u p + f 2 ( x ) u 2 m * 1 .
They addressed the existence of at least two positive solutions depending on the parameters by exploiting the Nehari manifold. do Ó et al. [15] studied a fractional stationary Kirchhoff equation involving two nonlinear source terms of different signs:
( a + b [ u ] m 2 θ 2 ) ( Δ ) m u = λ f ( x ) | u | q 2 u | u | 2 m * 2 u .
By using a variational approach based on the Nehari manifold, they obtained the existence of two positive solutions for suitable values of the parameters. Zhang et al. [16] investigated a fractional stationary Kirchhoff equation involving the nonlocal integro-differential operator:
h R 2 N | u ( x ) u ( y ) | 2 K ( x y ) d x d y L K u = f ( x , u ) ,
where K : R N \ { 0 } ( 0 , ) is a function satisfying certain assumptions, and L K denotes the nonlocal integral-differential operator. Equation (4) includes the following fractional Kirchhoff equation as a particular case.
h ( [ u ] m 2 ) ( Δ ) m u = f ( x , u ) .
By computing the critical groups at zero and at infinity, they derived the existence of at least one nontrivial solution via Morse theory. Molica Bisci and Vilasi [17] considered the following fractional stationary Kirchhoff equation involving two nonlinear source terms:
h ( [ u ] m 2 ) L K u = λ f 1 ( x , u ) + μ f 2 ( x , u ) .
By exploiting an abstract critical point theorem for smooth functionals, they obtained the existence of at least three solutions for suitable values of the parameters. Additionally, in the autonomous case, they gave a precise estimate for the range of these parameters by using some properties of the fractional calculus on a specific family of test functions. Concerning the fractional evolution Kirchhoff equations, Xiang et al. [18] studied
u t + h ( [ u ] m 2 ) ( Δ ) m u = | u | p 2 u .
They obtained the existence of non-negative local solutions by using the Galerkin approximations and proved the blow-up of non-negative local solutions with suitable initial data by virtue of a differential inequality technique. Lin et al. [19] studied a fractional evolution Kirchhoff equation of the form
u t t + [ u ] m 2 θ 2 ( Δ ) m u = f ( u ) .
They utilized the concavity arguments to obtain the blow-up of solutions. Pan et al. [20] considered the following fractional evolution Kirchhoff equation with nonlinear weak damping:
u t t + [ u ] m 2 θ 2 ( Δ ) m u + | u t | r 2 u t + u = | u | q 2 u .
They obtained the global existence, vacuum isolation, asymptotic behavior, and blow-up of solutions by using the theory of potential wells. Recently, Xiang and Hu [21] conducted an investigation on the following fractional viscoelastic Kirchhoff equation:
u t t + h ( [ u ] m 2 ) ( Δ ) m u 0 t g ( t τ ) ( Δ ) m u ( τ ) d τ + ( Δ ) s u t = λ | u | q 2 u .
They utilized the Galerkin approximations to establish the existence of local and global solutions, and they employed the concavity arguments to derive the blow-up of solutions.
Inspired by the above works, we deal with the problem (1)–(3). There are two features about our results. On the one hand, our results are independent of any parameters. On the other hand, we focus on the effects of the nonlinear source terms of different signs on solutions. From a physical perspective, the two nonlinear source terms of different signs actually represent the two opposing external forces acting on the viscoelastic string with fractional length. In the absence of the external force | u | q 2 u , it is trivial that the problem (1)–(3) has a global solution with the asymptotic behavior by the arguments similar to Ref. [22]. Our results show that even though the external force | u | q 2 u appears, the problem (1)–(3) can still have a global solution with the asymptotic behavior. In addition, our main technical tool is the theory of potential wells that has been widely employed to analyze the qualitative properties of solutions of evolution equations. In this regard, in addition to the work of Ref. [20] mentioned above, we refer to the following work. Xu et al. [23] used the theory of potential wells to investigate the global existence and blow-up of solutions of a sixth-order nonlinear hyperbolic equation. Cavalcanti and Domingos Cavalcanti [24] modified the theory of potential wells to study the global existence and asymptotic behaviour of solutions of a nonlinear evolution equation. Gazzola and Squassina [25] improved the theory of potential wells to obtain the global existence, asymptotic behaviour, and blow-up of solutions of a damped semilinear wave equation. Liu et al. [26] applied the idea of Ref. [25] to deal with a fourth-order damped nonlinear hyperbolic equation. Xu and Su [27], Luo et al. [28], and Liu and Li [29] introduced a family of potential wells to study the equation under consideration, respectively. Liu et al. [30] introduced a family of potential wells that were different from those in Refs. [27,28,29]. In contrast, the definition of the potential well in this paper differs from those in the above studies.
This paper is organized as follows. Section 1 serves as an introduction. Section 2 provides a comprehensive overview of the relevant functional spaces. Additionally, we define a potential well and give its properties. Section 3 is devoted to the proof of the existence of global solutions by using the theory of potential wells and the Galerkin approximations. In Section 4, we prove the asymptotic behavior of global solutions by using the theory of potential wells and the perturbed energy method [31,32,33]. In Section 5, we summarize our main results.

2. Preliminaries

As in Refs. [34,35], we denote by X the linear space of Lebesgue measurable functions from R N to R such that the restriction to Ω of any function u in X belongs to L 2 ( Ω ) and
Q | u ( x ) u ( y ) | 2 | x y | N + 2 m d x d y < ,
where Q : = R 2 N \ ( C Ω × C Ω ) and C Ω : = R N \ Ω . The space X is endowed with
u X : = u L 2 ( Ω ) + Q | u ( x ) u ( y ) | 2 | x y | N + 2 m d x d y 1 2 .
It is straightforward to verify that · X is a norm on X. We introduce the following closed linear subspace of X,
X 0 : = { u X | u = 0 a . e . in C Ω } .
This is a Hilbert space equipped with the inner product
( u , v ) X 0 : = Q ( u ( x ) u ( y ) ) ( v ( x ) v ( y ) ) | x y | N + 2 m d x d y
and the norm
u X 0 : = Q | u ( x ) u ( y ) | 2 | x y | N + 2 m d x d y 1 2
which is equivalent to u X . Moreover, the embedding X 0 L q ( Ω ) is continuous for any 1 q 2 m * and compact for any 1 q < 2 m * , where
2 m * = 2 N N 2 m if 2 m < N , if 2 m N .
In this paper, the exponents q and r satisfy the following assumption:
( A 1 ) 2 < r < q < 2 m * .
In addition, as in Ref. [33], the memory kernel g satisfies
( A 2 ) g C 1 ( R + ) L 1 ( R + ) , g ( t ) 0 , g ( t ) 0 for all t [ 0 , ) , and κ : = a 0 g ( t ) d t > 0 .
For simplicity, we denote
· p : = · L p ( Ω ) , ( u , v ) : = Ω u v d x
and
( g u ) ( t ) : = 0 t g ( t τ ) u ( t ) u ( τ ) X 0 2 d τ .
Moreover, C > 0 denotes a generic constant.
Definition 1.
A function u L ( 0 , T ; X 0 ) with u t L ( 0 , T ; L 2 ( Ω ) ) is called the weak solution of the problem (1)–(3), if u ( 0 ) = u 0 in X 0 , u t ( 0 ) = u 1 in L 2 ( Ω ) , and
( u t ( t ) , w ) + 0 t ( a + b u ( τ ) X 0 2 p 2 ) ( u ( τ ) , w ) X 0 d τ 0 t 0 s g ( s τ ) ( u ( τ ) , w ) X 0 d τ d s + ( u ( t ) , w ) = ( u 1 , w ) + ( u 0 , w ) + 0 t ( | u | q 2 u , w ) d τ 0 t ( | u | r 2 u , w ) d τ
for any w X 0 and t ( 0 , T ] .
The energy function for the problem (1)–(3) is defined as follows:
E ( t ) : = 1 2 u t ( t ) 2 2 + b 2 p u ( t ) X 0 2 p + 1 2 a 0 t g ( τ ) d τ u ( t ) X 0 2 + 1 2 ( g u ) ( t ) 1 q u ( t ) q q + 1 r u ( t ) r r .
Following the idea from Ref. [36], the potential well and its boundary are defined by
W : = u X 0 u X 0 < 2 q ( q 2 ) κ d 1 2
and
𝜕 W : = u X 0 u X 0 = 2 q ( q 2 ) κ d 1 2 ,
where the depth of the potential well is
d : = q 2 2 q κ q q 2 C 1 2 q q 2 ,
and C 1 is the best Sobolev constant for the embedding X 0 L q ( Ω ) , i.e.,
C 1 : = sup u X 0 \ { 0 } u q u X 0 .
Lemma 1.
Let ( A 1 ) and ( A 2 ) be satisfied. Then,
(i) 
if u W and u X 0 0 , then κ u X 0 2 > u q q ;
(ii) 
if u 𝜕 W , then κ u X 0 2 u q q .
Proof. 
(i) By u W and (5), we obtain
u X 0 < 2 q ( q 2 ) κ d 1 2 .
From (7), we obtain
u X 0 < κ 1 q 2 C 1 q q 2 .
Noting that u X 0 0 , we have
κ u X 0 2 > C 1 q u X 0 q .
Therefore,
κ u X 0 2 > u q q .
(ii) By u 𝜕 W and (6), we obtain
u X 0 = 2 q ( q 2 ) κ d 1 2 .
By adopting similar arguments in the proof of (i), it is evident that
κ u X 0 2 u q q .

3. Existence of Global Solutions

Theorem 1.
Let ( A 1 ) and ( A 2 ) be satisfied. If u 0 W , u 1 L 2 ( Ω ) , and E ( 0 ) < d , then the problem (1)–(3) has a global solution u ( t ) W ¯ : = W 𝜕 W for all t ( 0 , ) .
Proof. 
We divide the proof of this theorem into three steps.
Step I. Galerkin approximations. Let { ω j } j = 1 be the orthogonal basis of X 0 and an orthonormal basis of L 2 ( Ω ) given by eigenfunctions of ( Δ ) m with boundary condition (3) (see Ref. [34], Proposition 9, for details). Denote W n = Span { ω 1 , ω 2 , , ω n } , n = 1 , 2 , . We seek the approximate solutions of the problem (1)–(3)
u n ( t ) = j = 1 n ξ j n ( t ) ω j , n = 1 , 2 , ,
which satisfy
( u n t t ( t ) , w ) + ( a + b u n ( t ) X 0 2 p 2 ) ( u n ( t ) , w ) X 0 0 t g ( t τ ) ( u n ( τ ) , w ) X 0 d τ + ( u n t ( t ) , w ) = ( | u n ( t ) | q 2 u n ( t ) , w ) ( | u n ( t ) | r 2 u n ( t ) , w ) , t > 0 ,
u n ( 0 ) = j = 1 n ξ j n ( 0 ) ω j u 0 in X 0 ,
u n t ( 0 ) = j = 1 n ξ j n ( 0 ) ω j u 1 in L 2 ( Ω ) ,
for any w W n . Let ξ n ( t ) = ( ξ 1 n ( t ) , ξ 2 n ( t ) , , ξ n n ( t ) ) T . Then, the vector function ξ n solves
ξ n ( t ) + ξ n ( t ) + L n ( t , ξ n ( t ) ) = 0 , t > 0 ,
ξ n ( 0 ) = ( ( u 0 , ω 1 ) , ( u 0 , ω 2 ) , , ( u 0 , ω n ) ) T ,
ξ n ( 0 ) = ( ( u 1 , ω 1 ) , ( u 1 , ω 2 ) , , ( u 1 , ω n ) ) T ,
where
L n ( t , ξ n ( t ) ) = ( L 1 n ( t , ξ n ( t ) ) , L 2 n ( t , ξ n ( t ) ) , , L n n ( t , ξ n ( t ) ) ) T ,
L i n ( t , ξ n ( t ) ) = a + b j = 1 n ξ j n ( t ) ω j X 0 2 p 2 j = 1 n ξ j n ( t ) ω j , ω i X 0 0 t g ( t τ ) j = 1 n ξ j n ( τ ) ω j , ω i X 0 d τ + j = 1 n ξ j n ( t ) ω j q 2 j = 1 n ξ j n ( t ) ω j , ω i j = 1 n ξ j n ( t ) ω j r 2 j = 1 n ξ j n ( t ) ω j , ω i .
In light of the standard theory of ordinary differential equations, the problem (12)–(14) has a solution ξ n C 2 [ 0 , T n ) , with T n T . Consequently, u n ( t ) defined by (8) satisfies the problem (9)–(11).
Step II. A priori estimates. Taking w = u n t ( t ) in (9), we obtain
d d t 1 2 u n t ( t ) 2 2 + a 2 u n ( t ) X 0 2 + b 2 p u n ( t ) X 0 2 p 0 t g ( t τ ) ( u n ( τ ) , u n t ( t ) ) X 0 d τ + u n t ( t ) 2 2 = 1 q d d t u n ( t ) q q 1 r d d t u n ( t ) r r .
Concerning the fourth term on the left hand side of (15), we have
0 t g ( t τ ) ( u n ( τ ) , u n t ( t ) ) X 0 d τ = 0 t g ( t τ ) ( u n ( τ ) u n ( t ) , u n t ( t ) ) X 0 d τ + 0 t g ( t τ ) ( u n ( t ) , u n t ( t ) ) X 0 d τ = 1 2 0 t g ( t τ ) d d t u n ( τ ) u n ( t ) X 0 2 d τ + 1 2 0 t g ( t τ ) d d t u n ( t ) X 0 2 d τ = 1 2 d d t ( g u n ) ( t ) 0 t g ( τ ) d τ u n ( t ) X 0 2 + 1 2 ( g u n ) ( t ) 1 2 g ( t ) u n ( t ) X 0 2 .
Substituting this equality into (15) and performing integration with respect to t, we obtain
E n ( t ) + 0 t u n τ ( τ ) 2 2 1 2 ( g u n ) ( τ ) + 1 2 g ( τ ) u n ( τ ) X 0 2 d τ = E n ( 0 )
for all t [ 0 , T ] , where
E n ( t ) = 1 2 u n t ( t ) 2 2 + b 2 p u n ( t ) X 0 2 p + 1 2 a 0 t g ( τ ) d τ u n ( t ) X 0 2 + 1 2 ( g u n ) ( t ) 1 q u n ( t ) q q + 1 r u n ( t ) r r .
Based on the observations from (10) and (11), we conclude that E n ( 0 ) < d and u n ( 0 ) W for a sufficiently large n. We claim that
u n ( t ) W
for all t [ 0 , T ] and a sufficiently large n. Suppose that u n ( t ) W for some 0 < t < T . Then, there exists a time 0 < t 0 < T such that u n ( t 0 ) 𝜕 W and u n ( t ) W for all t [ 0 , t 0 ) . Therefore,
u n ( t 0 ) X 0 = 2 q ( q 2 ) κ d 1 2 .
Utilizing (17) in conjunction with part (ii) in Lemma 1, we are able to derive
E n ( t 0 ) 1 2 κ u n ( t 0 ) X 0 2 1 q u n ( t 0 ) q q = q 2 2 q κ u n ( t 0 ) X 0 2 + 1 q κ u n ( t 0 ) X 0 2 u n ( t 0 ) q q q 2 2 q κ u n ( t 0 ) X 0 2 = d .
In view of (16), we obtain a contradiction with E n ( 0 ) < d .
By (17), assertion (18) and (i) in Lemma 1, we can conclude that
E n ( t ) 1 2 u n t ( t ) 2 2 + 1 2 κ u n ( t ) X 0 2 1 q u n ( t ) q q = 1 2 u n t ( t ) 2 2 + q 2 2 q κ u n ( t ) X 0 2 + 1 q κ u n ( t ) X 0 2 u n ( t ) q q 1 2 u n t ( t ) 2 2 + q 2 2 q κ u n ( t ) X 0 2 ,
which, together with (16), gives
1 2 u n t ( t ) 2 2 + q 2 2 q κ u n ( t ) X 0 2 < d
for all t [ 0 , T ] . Thus,
u n t ( t ) 2 2 < 2 d
and
u n ( t ) X 0 2 < 2 q ( q 2 ) κ d .
Further, it can be derived from (21) that
| u n ( t ) | q 2 u n ( t ) q q 1 = u n ( t ) q q 1 C u n ( t ) X 0 q 1 < C 2 q ( q 2 ) κ d q 1 2
and
| u n ( t ) | r 2 u n ( t ) r r 1 = u n ( t ) r r 1 C u n ( t ) X 0 r 1 < C 2 q ( q 2 ) κ d r 1 2
for all t [ 0 , T ] . Estimates (20)–(23) imply that
{ u n } is bounded in L ( 0 , T ; X 0 ) ,
{ u n t } is bounded in L ( 0 , T ; L 2 ( Ω ) ) ,
{ | u n | q 2 u n } is bounded in L ( 0 , T ; L q q 1 ( Ω ) ) ,
{ | u n | r 2 u n } is bounded in L ( 0 , T ; L r r 1 ( Ω ) ) .
Step III. Passage to the limit. By virtue of (24)–(27), there is a subsequence of { u n } (still represented by { u n } ) and a function u such that as n ,
u n u weakly star in L ( 0 , T ; X 0 ) ,
u n t u t weakly star in L ( 0 , T ; L 2 ( Ω ) ) ,
| u n | q 2 u n | u | q 2 u weakly star in L ( 0 , T ; L q q 1 ( Ω ) ) ,
| u n | r 2 u n | u | r 2 u weakly star in L ( 0 , T ; L r r 1 ( Ω ) ) .
Integrating (9) with respect to t, we obtain
( u n t ( t ) , w ) + 0 t ( a + b u n ( τ ) X 0 2 p 2 ) ( u n ( τ ) , w ) X 0 d τ 0 t 0 s g ( s τ ) ( u n ( τ ) , w ) X 0 d τ d s + ( u n ( t ) , w ) = ( u n t ( 0 ) , w ) + ( u n ( 0 ) , w ) + 0 t ( | u n ( τ ) | q 2 u n ( τ ) , w ) d τ 0 t ( | u n ( τ ) | r 2 u n ( τ ) , w ) d τ .
Taking n , we obtain
( u t ( t ) , w ) + 0 t ( a + b u ( τ ) X 0 2 p 2 ) ( u ( τ ) , w ) X 0 d τ 0 t 0 s g ( s τ ) ( u ( τ ) , w ) X 0 d τ d s + ( u ( t ) , w ) = ( u 1 , w ) + ( u 0 , w ) + 0 t ( | u ( τ ) | q 2 u ( τ ) , w ) d τ 0 t ( | u ( τ ) | r 2 u ( τ ) , w ) d τ .
From (10) and (11), we obtain u ( 0 ) = u 0 in X 0 and u t ( 0 ) = u 1 in L 2 ( Ω ) . Hence, u is a global solution of the problem (1)–(3). Furthermore, from (28), we obtain
u ( t ) X 0 lim inf n u n ( t ) X 0 ,
which, together with (21), gives
u ( t ) X 0 2 q ( q 2 ) κ d 1 2 ,
that is, u ( t ) W ¯ for all t ( 0 , ) . □

4. Asymptotic Behavior of Global Solutions

Theorem 2.
In addition to all the assumptions of Theorem 1, assume that there exists a constant ρ > 0 such that g ( t ) ρ g ( t ) for all t [ 0 , ) . Then,
u ( t ) X 0 2 + u t ( t ) 2 2 α e β t , t [ 0 , ) ,
for some constants α , β > 0 .
Proof. 
For the approximate solutions given in the proof of Theorem 1, we construct
L ( t ) = E n ( t ) + ε Ψ ( t ) , t [ 0 , ) ,
where Ψ ( t ) = ( u n ( t ) , u n t ( t ) ) , and ε > 0 is a constant to be determined later.
Step I. We now claim that there exist two constants γ 1 , γ 2 > 0 , depending on ε , such that
γ 1 E n ( t ) L ( t ) γ 2 E n ( t ) , t [ 0 , ) .
In fact, according to Cauchy’s inequality, we obtain
| Ψ ( t ) | 1 2 u n ( t ) 2 2 + 1 2 u n t ( t ) 2 2 ,
and so,
| Ψ ( t ) | C 2 2 2 u n ( t ) X 0 2 + 1 2 u n t ( t ) 2 2 ,
where C 2 is the best Sobolev constant for the embedding X 0 L 2 ( Ω ) . Combining (33) and (19), we obtain | Ψ ( t ) | M E n ( t ) for some constant M > 0 independent of n, which, together with (31), yields that assertion (32) holds.
Step II. We next claim that
L ( t ) ε η E n ( t )
for sufficiently small η and ε . Indeed, since
E n ( t ) = 1 2 ( g u n ) ( t ) 1 2 g ( t ) u n ( t ) X 0 2 u n t ( t ) 2 2
and
Ψ ( t ) = u n t ( t ) 2 2 b u n ( t ) X 0 2 p a u n ( t ) X 0 2 + 0 t g ( t τ ) ( u n ( τ ) , u n ( t ) ) X 0 d τ ( u n ( t ) , u n t ( t ) ) + u n ( t ) q q u n ( t ) r r ,
it follows from (31) that
L ( t ) = 1 2 ( g u n ) ( t ) 1 2 g ( t ) u n ( t ) X 0 2 u n t ( t ) 2 2 + ε u n t ( t ) 2 2 ε b u n ( t ) X 0 2 p ε a u n ( t ) X 0 2 + ε 0 t g ( t τ ) ( u n ( τ ) , u n ( t ) ) X 0 d τ ε ( u n ( t ) , u n t ( t ) ) + ε u n ( t ) q q ε u n ( t ) r r .
For the seventh term on the right hand side of (35), we can derive from Schwarz’s inequality and Cauchy’s inequality with ϵ 1 > 0 that
0 t g ( t τ ) ( u n ( τ ) , u n ( t ) ) X 0 d τ = 0 t g ( t τ ) u n ( t ) X 0 2 d τ + 0 t g ( t τ ) ( u n ( τ ) u n ( t ) , u n ( t ) ) X 0 d τ 0 t g ( τ ) d τ u n ( t ) X 0 2 + ϵ 1 0 t g ( τ ) d τ u n ( t ) X 0 2 + 1 4 ϵ 1 ( g u n ) ( t ) ( a κ ) u n ( t ) X 0 2 + ϵ 1 ( a κ ) u n ( t ) X 0 2 + 1 4 ϵ 1 ( g u n ) ( t ) .
For the eighth term on the right-hand side of (35), we deduce from Cauchy’s inequality with ϵ 2 > 0 that
( u n ( t ) , u n t ( t ) ) ϵ 2 u n ( t ) 2 2 + 1 4 ϵ 2 u n t ( t ) 2 2 ϵ 2 C 2 2 u n ( t ) X 0 2 + 1 4 ϵ 2 u n t ( t ) 2 2 .
Plugging (36) and (37) into (35), we obtain
L ( t ) ε + ε 4 ϵ 2 1 u n t ( t ) 2 2 ε b u n ( t ) X 0 2 p + ε ( ϵ 1 ( a κ ) + ϵ 2 C 2 2 κ ) u n ( t ) X 0 2 + ε 4 ϵ 1 ρ 2 ( g u n ) ( t ) + ε u n ( t ) q q ε u n ( t ) r r ,
and so,
L ( t ) ε η E n ( t ) + ε + ε 4 ϵ 2 + ε η 2 1 u n t ( t ) 2 2 + ε b η 2 p 1 u n ( t ) X 0 2 p + ε ϵ 1 ( a κ ) + ϵ 2 C 2 2 + η a 2 κ u n ( t ) X 0 2 + ε 4 ϵ 1 + ε η 2 ρ 2 ( g u n ) ( t ) + ε u n ( t ) q q ε u n ( t ) r r ε η q u n ( t ) q q + ε η r u n ( t ) r r ,
where η > 0 is a constant to be determined later. We conclude from (16) and (19) that
E n ( 0 ) q 2 2 q κ u n ( t ) X 0 2 .
This leads to
u n ( t ) X 0 2 q ( q 2 ) κ E n ( 0 ) 1 2 .
Hence,
u n ( t ) q q C 1 q u n ( t ) X 0 q 2 u n ( t ) X 0 2 C 1 q 2 q ( q 2 ) κ E n ( 0 ) q 2 2 u n ( t ) X 0 2 .
The substitution of this inequality into (38) yields
L ( t ) ε η E n ( t ) + ε + ε 4 ϵ 2 + ε η 2 1 u n t ( t ) 2 2 + ε b η 2 p 1 u n ( t ) X 0 2 p + ε ϵ 1 ( a κ ) + ϵ 2 C 2 2 + η a 2 + C 1 q 2 q ( q 2 ) κ E n ( 0 ) q 2 2 κ u n ( t ) X 0 2 + ε 4 ϵ 1 + ε η 2 ρ 2 ( g u n ) ( t ) ε u n ( t ) r r + ε η r u n ( t ) r r .
Since
C 1 q 2 q ( q 2 ) κ E n ( 0 ) q 2 2 < C 1 q 2 q ( q 2 ) κ d q 2 2 = κ ,
the values of ϵ 1 , ϵ 2 , and η can be selected to be sufficiently small, such that η < min { 2 p , r } and
ϵ 1 ( a κ ) + ϵ 2 C 2 2 + η a 2 + C 1 q 2 q ( q 2 ) κ E n ( 0 ) q 2 2 κ < 0 .
Consequently, it follows from (39) that
L ( t ) ε η E n ( t ) + ε + ε 4 ϵ 2 + ε η 2 1 u n t ( t ) 2 2 + ε 4 ϵ 1 + ε η 2 ρ 2 ( g u n ) ( t ) .
Thus, for fixed ϵ 1 , ϵ 2 , and η , we can choose
ε < min 1 M , 4 ϵ 2 4 ϵ 2 + 1 + 2 η ϵ 2 , 2 ρ ϵ 1 1 + 2 η ϵ 1
such that assertion (34) holds.
Step III. We prove (30). To accomplish this, by combining (34) with the second inequality stated in assertion (32), we obtain
L ( t ) ε η γ 2 L ( t ) .
Hence,
L ( t ) C e ε η γ 2 t , t [ 0 , ) .
The first inequality in assertion (32) allows us to further deduce that
E n ( t ) C γ 1 e ε η γ 2 t , t [ 0 , ) .
It can be inferred from (28) and (29) that
u ( t ) X 0 2 + u t ( t ) 2 2 lim inf n u n ( t ) X 0 2 + u n t ( t ) 2 2 ,
which, together with (19) and (40), gives
lim inf n u n ( t ) X 0 2 + u n t ( t ) 2 2 lim inf n C E n ( t ) C γ 1 e ε η γ 2 t ,
and so, (30) is derived, where α = C / γ 1 and β = ε η / γ 2 . Theorem 2 is proved. □
Theorem 2 provides a sufficient condition for the asymptotic behavior of solutions of the problem (1)–(3). For example, we take N = 1 , m = 1 / 2 , a = 1 , b = 0 , q = 4 , r = 3 , and g ( t ) = g ( 0 ) e t with g ( 0 ) < 1 . Then, it is easy to verify that ( A 1 ) , ( A 2 ) , and g ( t ) ρ g ( t ) hold, where κ = 1 g ( 0 ) and ρ = 1 . By virtue of (7) and (5), we obtain
d = 1 4 ( 1 g ( 0 ) ) 2 C 1 4
and
W = u X 0 u X 0 < ( 1 g ( 0 ) ) 1 2 C 1 2 .
Therefore, if u 0 and u 1 satisfy
u 0 X 0 < ( 1 g ( 0 ) ) 1 2 C 1 2
and
1 2 u 1 2 2 + 1 2 u 0 X 0 2 1 4 u 0 4 4 + 1 3 u 0 3 3 < 1 4 ( 1 g ( 0 ) ) 2 C 1 4 ,
then, we conclude from Theorem 2 that the norm of global solutions of the problem (1)–(3) in the phase space X 0 × L 2 ( Ω ) decays exponentially to zero when the time tends to infinity.

5. Conclusions

The initial boundary value problem for a class of fractional viscoelastic Kirchhoff equations involving two nonlinear source terms of different signs is studied by this paper. In the case where the two opposing external forces | u | q 2 u and | u | r 2 u appear simultaneously, the existence and asymptotic behavior of global solutions are derived, namely, Theorems 1 and 2. More specifically, Theorem 1 shows that if the initial data u 0 lie in the potential well, and the initial energy is less than the depth of the potential well, then the problem (1)–(3) has a global solution. Theorem 2 shows us that if the memory kernel decays exponentially, then the global solutions of the problem (1)–(3) also do this.
Our main technical tool is the theory of potential wells, which is different from the classical. We describe the potential well by a sphere, whose radius is expressed by d. Although the depth of the potential well d is smaller than the classical, the spatial structure of the potential well is clearer so that it is not necessary to introduce the Nehari functional and the Nehari manifold.
In the future, we will focus on the study of the qualitative properties of solutions of fractional viscoelastic Kirchhoff equations involving more general Kirchhoff functions and nonlinear source terms.

Author Contributions

L.Z. and Y.L. contributed to each part of this study equally. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by the Talent Introduction Research Project of Northwest Minzu University (Grant No. xbmuyjrc2021008), the Innovation Team Project of Northwest Minzu University (Grant No. 1110130131), the First-Rate Discipline of Northwest Minzu University (Grant No. 2019XJYLZY-02), and the Key Laboratory of China’s Ethnic Languages and Information Technology of Ministry of Education at Northwest Minzu University.

Data Availability Statement

Data sharing is not applicable to this article as no datasets were generated or analyzed during the current study.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Fiscella, A.; Valdinoci, E. A critical Kirchhoff type problem involving a nonlocal operator. Nonlinear Anal. 2014, 94, 156–170. [Google Scholar] [CrossRef]
  2. Torrejón, R.; Yong, J.M. On a quasilinear wave equation with memory. Nonlinear Anal. 1991, 16, 61–78. [Google Scholar] [CrossRef]
  3. Wu, S.T.; Tsai, L.Y. Blow-up of positive-initial-energy solutions for an integro-differential equation with nonlinear damping. Taiwan. J. Math. 2010, 14, 2043–2058. [Google Scholar] [CrossRef]
  4. Liu, G.W.; Hou, C.S.; Guo, X.L. Global nonexistence for nonlinear Kirchhoff systems with memory term. Z. Angew. Math. Phys. 2014, 65, 925–940. [Google Scholar] [CrossRef]
  5. Caffarelli, L. Non-local diffusions, drifts and games. In Proceedings of the Nonlinear Partial Differential Equations The Abel Symposium 2010, Oslo, Norway, 28 September–1 October 2010; Springer: Berlin/Heidelberg, Germany, 2012; pp. 37–52. [Google Scholar] [CrossRef]
  6. Applebaum, D. Lévy processes-from probability to finance and quantum groups. Not. Amer. Math. Soc. 2004, 51, 1336–1347. [Google Scholar]
  7. Laskin, N. Fractional quantum mechanics and Lévy path integrals. Phys. Lett. A 2000, 268, 298–305. [Google Scholar] [CrossRef]
  8. Valdinoci, E. From the long jump random walk to the fractional Laplacian. arXiv 2009, arXiv:0901.3261. [Google Scholar]
  9. Vázquez, J.L. Recent progress in the theory of nonlinear diffusion with fractional laplacian operators. Discret. Contin. Dyn. Syst. 2014, 7, 857–885. [Google Scholar] [CrossRef]
  10. Nguyen, H.T.; Tuan, N.A.; Yang, C. Global well-posedness for fractional Sobolev-Galpern type equations. Discret. Contin. Dyn. Syst. 2022, 42, 2637–2665. [Google Scholar] [CrossRef]
  11. Nguyen, A.T.; Tuan, N.H.; Yang, C. On Cauchy problem for fractional parabolic-elliptic Keller-Segel model. Adv. Nonlinear Anal. 2023, 12, 97–116. [Google Scholar] [CrossRef]
  12. Ambrosio, V.; Isernia, T. A multiplicity result for a fractional Kirchhoff equation in RN with a general nonlinearity. Commun. Contemp. Math. 2018, 20, 1750054. [Google Scholar] [CrossRef]
  13. Nyamoradi, N.; Ambrosio, V. Existence and non-existence results for fractional Kirchhoff Laplacian problems. Anal. Math. Phys. 2021, 11, 125. [Google Scholar] [CrossRef]
  14. Fiscella, A.; Mishra, P.K. Fractional Kirchhoff Hardy problems with singular and critical Sobolev nonlinearities. Manuscr. Math. 2022, 168, 257–301. [Google Scholar] [CrossRef]
  15. do Ó, J.M.; He, X.M.; Mishra, P.K. Fractional Kirchhoff problem with critical indefinite nonlinearity. Math. Nachr. 2019, 292, 615–632. [Google Scholar] [CrossRef]
  16. Zhang, B.L.; Rădulescu, V.D.; Wang, L. Existence results for Kirchhoff-type superlinear problems involving the fractional Laplacian. Proc. R. Soc. Edinb. Sect. A 2019, 149, 1061–1081. [Google Scholar] [CrossRef]
  17. Bisci, G.M.; Vilasi, L. On a fractional degenerate Kirchhoff-type problem. Commun. Contemp. Math. 2017, 19, 1550088. [Google Scholar] [CrossRef]
  18. Xiang, M.Q.; Rădulescu, V.D.; Zhang, B.L. Nonlocal Kirchhoff diffusion problems: Local existence and blow-up of solutions. Nonlinearity 2018, 31, 3228–3250. [Google Scholar] [CrossRef]
  19. Lin, Q.; Tian, X.T.; Xu, R.Z.; Zhang, M.N. Blow up and blow up time for degenerate Kirchhoff- type wave problems involving the fractional Laplacian with arbitrary positive initial energy. Discret. Contin. Dyn. Syst. Ser. S 2020, 13, 2095–2107. [Google Scholar] [CrossRef]
  20. Pan, N.; Pucci, P.; Xu, R.Z.; Zhang, B.L. Degenerate Kirchhoff-type wave problems involving the fractional Laplacian with nonlinear damping and source terms. J. Evol. Equ. 2019, 19, 615–643. [Google Scholar] [CrossRef]
  21. Xiang, M.Q.; Hu, D. Existence and blow-up of solutions for fractional wave equations of Kirchhoff type with viscoelasticity. Discret. Contin. Dyn. Syst. Ser. 2021, 14, 4609–4629. [Google Scholar] [CrossRef]
  22. Cavalcanti, M.M.; Cavalcanti, V.N.D.; Ma, T.F. Exponential decay of the viscoelastic Euler–Bernoulli equation with a nonlocal dissipation in general domains. Differ. Integral Equ. 2004, 17, 495–510. [Google Scholar] [CrossRef]
  23. Xu, R.Z.; Zhang, M.Y.; Chen, S.H.; Yang, Y.B.; Shen, J.H. The initial-boundary value problems for a class of six order nonlinear wave equation. Discret. Contin. Dyn. Syst. 2017, 37, 5631–5649. [Google Scholar] [CrossRef]
  24. Cavalcanti, M.M.; Cavalcanti, V.N.D. Existence and asymptotic stability for evolution problems on manifolds with damping and source terms. J. Math. Anal. Appl 2004, 291, 109–127. [Google Scholar] [CrossRef]
  25. Gazzola, F.; Squassina, M. Global solutions and finite time blow up for damped semilinear wave equations. Ann. Inst. Henri Poincare Anal. Non Lineaire 2006, 3, 185–207. [Google Scholar] [CrossRef]
  26. Liu, Y.; Mu, J.; Jiao, Y.J. A class of fourth order damped wave equations with arbitrary positive initial energy. Proc. Edinb. Math. Soc. 2019, 62, 165–178. [Google Scholar] [CrossRef]
  27. Xu, R.Z.; Su, J. Global existence and finite time blow-up for a class of semilinear pseudo-parabolic equations. J. Funct. Anal. 2013, 264, 2732–2763. [Google Scholar] [CrossRef]
  28. Luo, Y.B.; Xu, R.Z.; Yang, C. Global well-posedness for a class of semilinear hyperbolic equations with singular potentials on manifolds with conical singularities. Calc. Var. Partial. Differ. Equ. 2022, 61, 210. [Google Scholar] [CrossRef]
  29. Liu, Y.; Li, W.K. A class of fourth-order nonlinear parabolic equations modeling the epitaxial growth of thin films. Discret. Contin. Dyn. Syst. Ser. S 2021, 14, 4367–4381. [Google Scholar] [CrossRef]
  30. Liu, Y.; Yu, T.; Li, W.K. Global well-posedness, asymptotic behavior and blow-up of solutions for a class of degenerate parabolic equations. Nonlinear Anal. 2020, 196, 111759. [Google Scholar] [CrossRef]
  31. Haraux, A.; Zuazua, E. Decay estimates for some semilinear damped hyperbolic problems. Arch. Ration. Mech. Anal. 1988, 100, 191–206. [Google Scholar] [CrossRef]
  32. Liu, Y. Global attractors for a nonlinear plate equation modeling the oscillations of suspension bridges. Commun. Anal. Mech. 2023, 15, 436–456. [Google Scholar] [CrossRef]
  33. Liu, Y.; Moon, B.; Rădulescu, V.D.; Xu, R.Z.; Yang, C. Qualitative properties of solution to a viscoelastic Kirchhoff-like plate equation. J. Math. Phys. 2023, 64, 051511. [Google Scholar] [CrossRef]
  34. Servadei, R.; Valdinoci, E. Variational methods for non-local operators of elliptic type. Discret. Contin. Dyn. Syst. 2013, 33, 2105–2137. [Google Scholar] [CrossRef]
  35. Servadei, R.; Valdinoci, E. The Brezis-Nirenberg result for the fractional Laplacian. Trans. Amer. Math. Soc. 2015, 367, 67–102. [Google Scholar] [CrossRef]
  36. Liu, Y.; Lv, P.J.; Da, C.J. Blow-up of a nonlocal p-Laplacian evolution equation with critical initial energy. Ann. Polon. Math. 2016, 117, 89–99. [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

Zhang, L.; Liu, Y. A Class of Fractional Viscoelastic Kirchhoff Equations Involving Two Nonlinear Source Terms of Different Signs. Axioms 2024, 13, 169. https://doi.org/10.3390/axioms13030169

AMA Style

Zhang L, Liu Y. A Class of Fractional Viscoelastic Kirchhoff Equations Involving Two Nonlinear Source Terms of Different Signs. Axioms. 2024; 13(3):169. https://doi.org/10.3390/axioms13030169

Chicago/Turabian Style

Zhang, Li, and Yang Liu. 2024. "A Class of Fractional Viscoelastic Kirchhoff Equations Involving Two Nonlinear Source Terms of Different Signs" Axioms 13, no. 3: 169. https://doi.org/10.3390/axioms13030169

APA Style

Zhang, L., & Liu, Y. (2024). A Class of Fractional Viscoelastic Kirchhoff Equations Involving Two Nonlinear Source Terms of Different Signs. Axioms, 13(3), 169. https://doi.org/10.3390/axioms13030169

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