Next Article in Journal
Flow of a Self-Similar Non-Newtonian Fluid Using Fractal Dimensions
Next Article in Special Issue
Mixed Convection of Fractional Nanofluids Considering Brownian Motion and Thermophoresis
Previous Article in Journal
A Novel Adaptive Fractional Differential Active Contour Image Segmentation Method
Previous Article in Special Issue
Existence of Solutions to a System of Riemann-Liouville Fractional Differential Equations with Coupled Riemann-Stieltjes Integrals Boundary Conditions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Initial Boundary Value Problem for a Fractional Viscoelastic Equation of the Kirchhoff Type

College of Mathematics and Computer Science, Northwest Minzu University, Lanzhou 730030, China
*
Author to whom correspondence should be addressed.
Submission received: 24 August 2022 / Revised: 25 September 2022 / Accepted: 6 October 2022 / Published: 11 October 2022

Abstract

:
In this paper, we study the initial boundary value problem for a fractional viscoelastic equation of the Kirchhoff type. In suitable functional spaces, we define a potential well. In the framework of the potential well theory, we obtain the global existence of solutions by using the Galerkin approximations. Moreover, we derive the asymptotic behavior of solutions by means of the perturbed energy method. Our main results provide sufficient conditions for the qualitative properties of solutions in time.

1. Introduction

In this paper, we study the following initial boundary value problem for a fractional viscoelastic equation of the Kirchhoff type:
u t t + h ( [ u ] m 2 ) ( Δ ) m u 0 t g ( t τ ) ( Δ ) m u ( τ ) d τ + u t = f ( 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 = R 2 N | u ( x , t ) u ( y , t ) | 2 | x y | N + 2 m d x d y 1 2
is the Gagliardo seminorm, ( Δ ) m is the fractional Laplace operator with 0 < m < 1 , and Ω R N ( N 1 ) is a bounded domain with a Lipschitz boundary. The unknown function u = u ( x , t ) is the vertical displacement of the small-amplitude vibrating viscoelastic string with the fractional length at position x and time t, 0 t g ( t τ ) ( Δ ) m u ( τ ) d τ is the viscoelastic term, u t is the weak damping term, the Kirchhoff function h ( s ) = 1 + s p 1 for all s 0 , p > 1 , and the source term f ( u ) = | u | q 2 u . The exponent q and the memory kernel g will be specified later.
For the classical viscoelastic wave equation of the Kirchhoff type, Wu and Tsai [1] studied the following equation:
u t t h ( u 2 2 ) Δ u + 0 t g ( t τ ) Δ u ( τ ) d τ Δ u t = f ( u ) .
They obtained the local existence, global existence, asymptotic behavior, and blow-up of solutions and provided the estimates on the decay rate of the energy function and the blow-up time of the solutions. Moreover, in [2], they considered the following viscoelastic wave equation of the Kirchhoff type with nonlinear weak damping:
u t t h ( u 2 2 ) Δ u + 0 t g ( t τ ) Δ u ( τ ) d τ + f 2 ( u t ) = f 1 ( u ) .
They obtained the local existence and blow-up of solutions and also derived the estimates of the blow-up times of the solutions.
When we examine the deep properties of real-world problems and extend them to other studies, some concepts usually have their own limitations. In this regard, many researchers pointed out the limitations of integer-order calculus while studying the systems related to non-Markovian mechanisms, hereditary properties, and other factors. In this situation, fractional calculus plays an important role, which is a generalization of classical calculus (see [3]). In recent years, fractional partial differential equations have attracted a great deal of attention due to their wide applicability in continuum mechanics, quantum and statistical mechanics, population dynamics, optimal control, game theory, and so on (see, for instance, [3,4,5,6,7,8,9,10,11] and the references therein). Fiscella and Valdinoci [12] proposed a fractional stationary Kirchhoff equation which models the vibration of a string with a fractional length by considering the nonlocal aspect of the tension. Subsequently, many fractional Kirchhoff equations were widely studied. Autuori et al. [13] investigated
h ( [ u ] m 2 ) L K u = λ f ( x , u ) + | u | q * 2 u ,
where L K is a fractional integro-differential operator, λ is a parameter, and q * is the critical exponent of the fractional Sobolev space H m ( R N ) . They proved the existence and asymptotic behavior of nonnegative solutions. Molica Bisci and Vilasi [14] dealt with
h ( [ u ] m 2 ) L K u = λ f ( x , u ) + μ g ( x , u ) ,
and derived the existence of at least three weak solutions for suitable values of the parameters by the variational approach. Moreover, they provided a concrete estimate for the range of these parameters in the autonomous case. Pucci et al. [15] investigated
h ( [ u ] m , p p ) ( Δ ) p m u + V ( x ) | u | p 2 u = λ ω ( x ) | u | q 2 u ν ( x ) | u | r 2 u ,
where ( Δ ) p m is the fractional p-Laplace operator, which may be defined along any φ C 0 ( R N ) as
( Δ ) p m φ ( x ) = 2 lim ε 0 + R N \ B ε ( x ) | φ ( x ) φ ( y ) | p 2 ( φ ( x ) φ ( y ) ) | x y | N + m p d y
for x R N and
[ u ] m , p = R 2 N | u ( x ) u ( y ) | p | x y | N + m p d x d y 1 p .
By using the variational approach and topological degree theory, they proved the multiplicity results depending on the parameter λ and under the suitable general integrability properties of the ratio between some powers of the weights. Moreover, they obtained the existence of infinitely many pairs of entire solutions by genus theory. Wang et al. [16] studied the following fractional Kirchhoff equation involving Choquard nonlinearity and singular nonlinearity:
( a + b ( [ u ] m , p ( θ 1 ) p ) ) ( Δ ) p m u = λ f 1 ( x ) | u | β + R N f 2 ( y ) | u ( y ) | q | x y | μ d y f 2 ( x ) u q 1 ,
where a, b, θ , λ , β , and μ are constants that meet certain conditions. They obtained the existence and multiplicity of nonnegative solutions by using the Nehari manifold approach combined with the Hardy–Littlehood–Sobolev inequality. Recently, Lin et al. [17] considered the fractional evolution Kirchhoff equation of the form
u t t + [ u ] m 2 ( θ 1 ) ( Δ ) m u = f ( u ) ,
and obtained the finite time blow-up of solutions with arbitrary positive initial energy by the concavity arguments.
Continuum mechanics attempts to describe the motions and equilibrium states of deformable bodies. Two types of materials are usually considered in basic texts on continuum mechanics: elastic materials and viscous fluids. At each material point of an elastic material, the stress at the present time depends only on the present value of the strain. On the other hand, for an incompressible viscous fluid, the stress at a given point is a function of the present value of the velocity gradient at that point (plus an undetermined pressure). Viscoelastic materials have properties between those of elastic materials and viscous fluids. Such materials have memory, where the stress depends not only on the present values of the strain or velocity gradient but also on the entire temporal history of motion (see [18]). Therefore, the research on the vibration of the viscoelastic string with a fractional length has important physical significance and scientific value. More recently, Xiang and Hu [19] investigated the following fractional viscoelastic equation of the Kirchhoff type:
u t t + h ( [ u ] m 2 ) ( Δ ) m u 0 t g ( t τ ) ( Δ ) m u ( τ ) d τ + ( Δ ) s u t = λ | u | q 2 u .
They proved the local and global existence of solutions by the Galerkin approximations and obtained the blow-up of solutions by the concavity arguments. However, to the best of our knowledge, much less effort has been devoted to similar studies.
Motivated by the above works, we would like to deal with the problems in Equations (1)–(3). In suitable functional spaces, we aim to study the global existence and asymptotic behavior of solutions in time. First of all, compared with [19], Equation (1) is non-degenerate due to the expression of the Kirchhoff function. Secondly, although we also evaluate the evolutional properties of solutions, we concentrate on the relationship between the initial data and them. In addition, our main method is the potential well theory that is different from classical ones. In the framework of our potential well theory, it is not necessary to introduce the Nehari functional or the Nehari manifold.
This paper is organized as follows. Section 1 is the introduction. In Section 2, we prepare the preliminary knowledge on the functional space. Applying the idea from [20], we define a potential well and provide its properties. Moreover, we display assumptions and notations corresponding to the problems in Equations (1)–(3). In Section 3, we introduce our main method in detail. In Section 4, we prove the global existence of solutions. Section 5 is devoted to the proof of the asymptotic behavior of the solutions by means of the perturbed energy method [21,22]. In Section 6, we summarize our main results.

2. Preliminaries

In this section, we first recall some necessary definitions and properties (see [23,24,25] for further details).
Let X be 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 easy to check that · X is a norm on X. Moreover, 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 ) * : = ( 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 * : = u X 0 = Q | u ( x ) u ( y ) | 2 | x y | N + 2 m d x d y 1 2 .
Here, u * is equivalent to u X .
The embedding X 0 L q ( Ω ) is continuous for any 1 q q * and compact for any 1 q < q * , where
q * = 2 N N 2 m if 2 m < N , if 2 m N .
In this paper, the exponent q of the source term satisfies the following assumption:
( A 1 ) 2 < q < q * .
Moreover, as in [26], 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 κ : = 1 0 g ( t ) d t > 0 .
For the sake of simplicity, we denote
· p : = · L p ( Ω ) , ( u , v ) : = Ω u v d x ,
and
( g u ) ( t ) : = 0 t g ( t τ ) u ( t ) u ( τ ) * 2 d τ .
Definition 1.
A function u L ( 0 , T ; X 0 ) with u t L ( 0 , T ; L 2 ( Ω ) ) is called a weak solution to Equations (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 h ( u ( τ ) * 2 ) ( u ( τ ) , w ) * d τ 0 t 0 s g ( s τ ) ( u ( τ ) , w ) * d τ d s + ( u ( t ) , w ) = ( u 1 , w ) + ( u 0 , w ) + 0 t ( f ( u ( τ ) ) , w ) d τ
for any w X 0 and t ( 0 , T ] .
We define the total energy function associated with the problems in Equations (1)–(3) as follows:
E ( t ) = 1 2 u t ( t ) 2 2 + 1 2 p u ( t ) * 2 p + 1 2 1 0 t g ( τ ) d τ u ( t ) * 2 + 1 2 ( g u ) ( t ) 1 q u ( t ) q q .
The potential well is
W = u X 0 u * < 2 q ( q 2 ) κ d 1 2
and its boundary is
W = u X 0 u * = 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 .
In addition, C 1 is the best Sobolev constant for the embedding X 0 L q ( Ω ) ; in other words, we have
C 1 = sup u X 0 \ { 0 } u q u * .
Lemma 1. 
Let ( A 1 ) and ( A 2 ) be fulfilled. Then, the following are true:
(i) 
If u W and u * 0 , then κ u * 2 > u q q ;
(ii) 
If u W , then κ u * 2 u q q .
Proof. 
(i) By u W and Equation (4), we have
u * < 2 q ( q 2 ) κ d 1 2 .
It follows from Equation (6) that
u * < κ 1 q 2 C 1 q q 2 .
Noting that u * 0 , we have
κ u * 2 > C 1 q u * q .
Hence, we obtain
κ u * 2 > u q q .
(ii) By u W and Equation (5), we find
u * = 2 q ( q 2 ) κ d 1 2 .
By the similar arguments in the proof of (i), it is easy to see that κ u * 2 u q q . □
The main results of this paper are proven in Section 4 and Section 5.

3. Methods

The potential well was first proposed by Sattinger [27] in order to study the global existence of solutions to a nonlinear hyperbolic equation. Subsequently, it was widely employed to analyze the qualitative properties of the solutions to evolution equations (see, for example, [18,28,29,30,31,32,33,34,35,36,37,38,39] and the references therein), and it has now developed into a theoretical system.
In general, by the energy functional J ( u ) and the Nehari functional I ( u ) , the classical potential well can usually be defined by
W = { u | J ( u ) < d , I ( u ) > 0 } { 0 } .
The critical points of J ( u ) are stationary solutions of the problem under consideration. Under appropriate assumptions, J ( u ) satisfies the Palais–Smale condition, and the problem under consideration admits at least a positive stationary solution whose energy d, namely the depth of the potential well, can be defined by
d = inf u N J ( u ) ,
where the Nehari manifold is
N = { u | I ( u ) = 0 } \ { 0 } .
In the present paper, we describe the potential well as a sphere (see Equation (4)) whose radius is expressed by d (see Equation (6)). Thus, the spatial structure of the potential well is clearer, and it is not necessary to introduce I ( u ) and N . As for the original definition and calculation process of d, we refer interested readers to [20].

4. Global Existence of Solutions

Theorem 1.
Let ( A 1 ) and ( A 2 ) be fulfilled. Assume that u 0 W , u 1 L 2 ( Ω ) , and E ( 0 ) < d . Then, Equations (1)–(3) admit a global solution u ( t ) W ¯ : = W W for all t ( 0 , ) .
Proof. 
Let { ω j } j = 1 be an orthogonal basis of X 0 and an orthonormal basis of L 2 ( Ω ) given by the eigenfunctions of ( Δ ) m with the boundary condition in Equation (3) (see [24] (Proposition 9) for details). Denote W n = Span { ω 1 , ω 2 , , ω n } , n = 1 , 2 , . We seek the approximate solutions to Equations (1)–(3), given by
u n ( t ) = j = 1 n ξ j n ( t ) ω j , n = 1 , 2 , ,
which satisfy
( u n t t ( t ) , w ) + h ( u n ( t ) * 2 ) ( u n ( t ) , w ) * 0 t g ( t τ ) ( u n ( τ ) , w ) * d τ + ( u n t ( t ) , w ) = ( f ( 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 ) ) = F n ( ξ n ( t ) ) , 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 ) ) = h j = 1 n ξ j n ( t ) ω j * 2 j = 1 n ξ j n ( t ) ω j , ω i * 0 t g ( t τ ) j = 1 n ξ j n ( τ ) ω j , ω i * d τ ,
F n ( ξ n ( t ) ) = ( F 1 n ( ξ n ( t ) ) , F 2 n ( ξ n ( t ) ) , , F n n ( ξ n ( t ) ) ) T ,
F i n ( ξ n ( t ) ) = f j = 1 n ξ j n ( t ) ω j , ω i .
In terms of standard theory for ODEs, the Cauchy problem in Equations (11)–(13) admits a solution ξ n C 2 [ 0 , T n ) with T n T . In turn, this gives a solution u n ( t ) defined by Equation (7) and satisfying Equations (8)–(10). The following estimates will allow us to extend the local solution to [ 0 , T ] for any T > 0 .
By using w = u n t ( t ) in Equation (8), we obtain
d d t 1 2 u n t ( t ) 2 2 + 1 2 u n t ( t ) * 2 + 1 2 p u n ( t ) * 2 p 0 t g ( t τ ) ( u n ( τ ) , u n t ( t ) ) * d τ + u n t ( t ) 2 2 = 1 q d d t u n ( t ) q q .
Note that
0 t g ( t τ ) ( u n ( τ ) , u n t ( t ) ) * d τ = 0 t g ( t τ ) ( u n ( τ ) u n ( t ) , u n t ( t ) ) * d τ + 0 t g ( t τ ) ( u n ( t ) , u n t ( t ) ) * d τ = 1 2 0 t g ( t τ ) d d t u n ( τ ) u n ( t ) * 2 d τ + 1 2 0 t g ( t τ ) d d t u n ( t ) * 2 d τ = 1 2 d d t ( g u n ) ( t ) 0 t g ( τ ) d τ u n ( t ) * 2 + 1 2 ( g u n ) ( t ) 1 2 g ( t ) u n ( t ) * 2 .
By substituting this equality into Equation (14) and integrating it with respect to t, we deduce that
E n ( t ) + 0 t u n τ ( τ ) 2 2 1 2 ( g u n ) ( τ ) + 1 2 g ( τ ) u n ( τ ) * 2 d τ = E n ( 0 )
for all t [ 0 , T ] , where
E n ( t ) = 1 2 u n t ( t ) 2 2 + 1 2 p u n ( t ) * 2 p + 1 2 1 0 t g ( τ ) d τ u n ( t ) * 2 + 1 2 ( g u n ) ( t ) 1 q u n ( t ) q q .
In light of Equations (9) and (10), we infer that E n ( 0 ) < d and u n ( 0 ) W for a sufficiently large n. We now 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 ) . Hence, we obtain
u n ( t 0 ) * = 2 q ( q 2 ) κ d 1 2 .
Through Equation (16) and (ii) in Lemma 1, we obtain
E n ( t 0 ) 1 2 κ u n ( t 0 ) * 2 1 q u n ( t 0 ) q q = q 2 2 q κ u n ( t 0 ) * 2 + 1 q κ u n ( t 0 ) * 2 u n ( t 0 ) q q q 2 2 q κ u n ( t 0 ) * 2 = d ,
which contradicts E n ( 0 ) < d according to Equation (15).
From Equation (16), the assertion in Equation (17), and (i) in Lemma 1, it follows that
E n ( t ) 1 2 u n t ( t ) 2 2 + 1 2 κ u n ( t ) * 2 1 q u n ( t ) q q = 1 2 u n t ( t ) 2 2 + q 2 2 q κ u n ( t ) * 2 + 1 q κ u n ( t ) * 2 u n ( t ) q q 1 2 u n t ( t ) 2 2 + q 2 2 q κ u n ( t ) * 2 ,
which, together with Equation (15), gives
1 2 u n t ( t ) 2 2 + q 2 2 q κ u n ( t ) * 2 < d
for all t [ 0 , T ] . Thus, for all t [ 0 , T ] , we find
u n t ( t ) 2 2 < 2 d
and
u n ( t ) * 2 < 2 q ( q 2 ) κ d .
Furthermore, we deduce from Equation (19) that
f ( u n ( t ) ) r r = u n ( t ) q q C 1 q u n ( t ) * q < C 1 q 2 q ( q 2 ) κ d q 2
for all t [ 0 , T ] , where r = q q 1 .
The above estimates mean the following:
{ u n } is bounded in L ( 0 , T ; X 0 ) ,
{ u n t } is bounded in L ( 0 , T ; L 2 ( Ω ) ) ,
{ f ( u n ) } is bounded in L ( 0 , T ; L r ( Ω ) ) .
Therefore, there exist u, χ , and a subsequence of { u n } , still denoted by { u n } , such that as n , the following are true:
u n u weakly star in L ( 0 , T ; X 0 ) ,
u n t u t weakly star in L ( 0 , T ; L 2 ( Ω ) ) ,
f ( u n ) χ weakly star in L ( 0 , T ; L r ( Ω ) ) .
Thus, we have the following:
u n u in L 2 ( 0 , T ; L 2 ( Ω ) ) and a . e . in Ω × [ 0 , T ] .
In terms of [32] (Chapter 1, Lemma 1.3), we have χ = f ( u ) .
Integrating Equation (8) with respect to t yields
( u n t ( t ) , w ) + 0 t h ( u n ( τ ) * 2 ) ( u n ( τ ) , w ) * d τ 0 t 0 s g ( s τ ) ( u n ( τ ) , w ) * d τ d s + ( u n ( t ) , w ) = ( u n t ( 0 ) , w ) + ( u n ( 0 ) , w ) + 0 t ( f ( u n ( τ ) ) , w ) d τ .
Using n , we further obtain
( u t ( t ) , w ) + 0 t h ( u n ( τ ) * 2 ) ( u ( τ ) , w ) * d τ 0 t 0 s g ( s τ ) ( u ( τ ) , w ) * d τ d s + ( u ( t ) , w ) = ( u 1 , w ) + ( u 0 , w ) + 0 t ( f ( u ( τ ) ) , w ) d τ .
By virtue of Equations (9) and (10), we have u ( 0 ) = u 0 in X 0 and u t ( 0 ) = u 1 in L 2 ( Ω ) . Therefore, u is a global solution to Equations (1)–(3). In addition, from Equation (20), we have
u ( t ) * lim inf n u n ( t ) * ,
which, together with Equation (19), tells us that
u ( t ) * 2 q ( q 2 ) κ d 1 2 .
In other words, u ( t ) W ¯ for all t ( 0 , ) . □

5. Asymptotic Behavior of the Solutions

Theorem 2.
In addition to all the assumptions of Theorem 1, suppose that there exists a constant ρ > 0 such that g ( t ) ρ g ( t ) for all t [ 0 , ) . Then, we have
u ( t ) * 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.
We now claim that there exist two constants γ i > 0 ( i = 1 , 2 ) , depending on ε , such that
γ 1 E n ( t ) L ( t ) γ 2 E n ( t ) , t [ 0 , ) .
Indeed, by virtue of Cauchy’s inequality, we find
| Ψ ( t ) | 1 2 u n ( t ) 2 2 + 1 2 u n t ( t ) 2 2 ,
and thus
| Ψ ( t ) | C 2 2 2 u n ( t ) * 2 + 1 2 u n t ( t ) 2 2 ,
where C 2 is the best Sobolev constant for the embedding X 0 L 2 ( Ω ) . By combining Equations (24) and (18), we obtain | Ψ ( t ) | C 1 E n ( t ) for some constant C 1 > 0 independent of n which, together with Equation (22), yields that the assertion in Equation (23) holds.
It can be said that
E n ( t ) = 1 2 ( g u n ) ( t ) 1 2 g ( t ) u n ( t ) * 2 u n t ( t ) 2 2 .
Then, a direct calculation gives
L ( t ) = 1 2 ( g u n ) ( t ) 1 2 g ( t ) u n ( t ) * 2 u n t ( t ) 2 2 + ε u n t ( t ) 2 2 ε u n ( t ) * 2 p ε u n ( t ) * 2 + ε 0 t g ( t τ ) ( u n ( τ ) , u n ( t ) ) * d τ ε ( u n ( t ) , u n t ( t ) ) + ε u n ( t ) q q .
For the seventh term on the right side of Equation (25), it follows from Schwarz’s inequality and Cauchy’s inequality with ϵ 1 > 0 that
0 t g ( t τ ) ( u n ( τ ) , u n ( t ) ) * d τ = 0 t g ( t τ ) u n ( t ) * 2 d τ + 0 t g ( t τ ) ( u n ( τ ) u n ( t ) , u n ( t ) ) * d τ 0 t g ( τ ) d τ u n ( t ) * 2 + ϵ 1 0 t g ( τ ) d τ u n ( t ) * 2 + 1 4 ϵ 1 ( g u n ) ( t ) ( 1 κ ) u n ( t ) * 2 + ϵ 1 ( 1 κ ) u n ( t ) * 2 + 1 4 ϵ 1 ( g u n ) ( t ) .
For the eighth term on the right side of Equation (25), it follows 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 ) * 2 + 1 4 ϵ 2 u n t ( t ) 2 2 .
Hence, we have
L ( t ) ε + ε 4 ϵ 2 1 u n t ( t ) 2 2 ε u n ( t ) * 2 p + ε ( ϵ 1 ( 1 κ ) + ϵ 2 C 2 2 κ ) u n ( t ) * 2 + ε 4 ϵ 1 ρ 2 ( g u n ) ( t ) + ε u n ( t ) q q ,
and so
L ( t ) ε η E n ( t ) + ε + ε 4 ϵ 2 + ε η 2 1 u n t ( t ) 2 2 + ε η 2 p 1 u n ( t ) * 2 p + ε ϵ 1 ( 1 κ ) + ϵ 2 C 2 2 + η 2 κ u n ( t ) * 2 + ε 4 ϵ 1 + ε η 2 ρ 2 ( g u n ) ( t ) + ε u n ( t ) q q ε η q u n ( t ) q q ,
where η > 0 is a constant to be determined later. It follows from Equations (15) and (18) that
E n ( 0 ) q 2 2 q κ u n ( t ) * 2 ,
which leads to
u n ( t ) * 2 q ( q 2 ) κ E n ( 0 ) 1 2 .
Hence, we have
u n ( t ) q q C 1 q u n ( t ) * q 2 u n ( t 0 ) * 2 C 1 q 2 q ( q 2 ) κ E n ( 0 ) q 2 2 u n ( t ) * 2 .
By substituting this inequality into Equation (26), we obtain
L ( t ) ε η E n ( t ) + ε + ε 4 ϵ 2 + ε η 2 1 u n t ( t ) 2 2 + ε η 2 p 1 u n ( t ) * 2 p + ε ϵ 1 ( 1 κ ) + ϵ 2 C 2 2 + η 2 + C 1 q 2 q ( q 2 ) κ E n ( 0 ) q 2 2 κ u n ( t ) * 2 + ε 4 ϵ 1 + ε η 2 ρ 2 ( g u n ) ( t ) .
Note that
C 1 q 2 q ( q 2 ) κ E n ( 0 ) q 2 2 < C 1 q 2 q ( q 2 ) κ d q 2 2 = κ .
We choose a sufficiently small ϵ i ( i = 1 , 2 ) and η such that η < 2 p and
ϵ 1 ( 1 κ ) + ϵ 2 C 2 2 + η 2 + C 1 q 2 q ( q 2 ) κ E n ( 0 ) q 2 2 κ 0 .
Thus, for a fixed ϵ i ( i = 1 , 2 ) and η , we can choose
ε < min 1 C 1 , 4 ϵ 2 4 ϵ 2 + 1 + 2 η ϵ 2 , 2 ρ ϵ 1 1 + 2 η ϵ 1
such that L ( t ) ε η E n ( t ) which, together with the second inequality in the assertion in Equation (23), gives L ( t ) ε η γ 2 L ( t ) . Hence, there exists a constant C 2 > 0 independent of n such that
L ( t ) C 2 e ε η γ 2 t , t [ 0 , ) .
We further conclude from the first inequality in the assertion in Equation (23) that
E n ( t ) C 2 γ 1 e ε η γ 2 t , t [ 0 , ) .
From Equations (20) and (21), it follows that
u ( t ) * 2 + u t ( t ) 2 2 lim inf n u n ( t 0 ) * 2 + u n t ( t ) 2 2 ,
which, combined with Equations (18) and (27), gives the conclusion of Theorem 2. □

6. Conclusions

In this paper, we studied the initial boundary value problem for a fractional viscoelastic equation of the Kirchhoff type. In the framework of the potential well theory, we established the global existence theorem, specifically Theorem 1. Under appropriate assumptions of the exponent of the source term and the memory kernel, it has been shown that if the initial data u 0 lies in the potential well, and the initial energy is less than the depth of the potential well, then the initial boundary value problem admits a global solution that lies in the closure of the potential well. Moreover, we have established the asymptotic behavior theorem, specifically Theorem 2. It is established that as the time variable tends toward infinity, the norm of the solutions in the phase space decays exponentially to zero at the same rate as the memory kernel. In light of the applications, once the initial data and the external force are effectively controlled, the vibration of the string with a fractional length and appropriate viscoelasticity will be stable. In this regard, the methods in [40] may be helpful.

Author Contributions

Investigation, Y.L. and L.Z.; Methodology, Y.L. and L.Z.; Project administration, Y.L.; Validation, Y.L.; Writing—original draft, Y.L. and L.Z.; Writing—review and editing, Y.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Fundamental Research Funds for the Central Universities (Grant No. 31920220062), the Science and Technology Plan Project of Gansu Province in China (Grant No. 21JR1RA200), the Talent Introduction Research Project of Northwest Minzu University (Grant No. xbmuyjrc2021008), and the Key Laboratory of China’s Ethnic Languages and Information Technology of the Ministry of Education at Northwest Minzu University.

Data Availability Statement

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

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wu, S.T.; Tsai, L.Y. On global existence and blow-up of solutions for an integro-differential equation with strong damping. Taiwan. J. Math. 2006, 10, 979–1014. [Google Scholar] [CrossRef]
  2. 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]
  3. Yao, S.W.; Ilhan, E.; Veeresha, P.; Baskonus, H.M. A powerful iterative approach for quintic complex Ginzburg-Landau equation within the frame of fractional operator. Fractals 2021, 29, 2140023. [Google Scholar] [CrossRef]
  4. Applebaum, D. Lévy processes-from probability to finance and quantum groups. Not. Am. Math. Soc. 2004, 51, 1336–1347. [Google Scholar]
  5. Baishya, C.; Veeresha, P. Laguerre polynomial-based operational matrix of integration for solving fractional differential equations with non-singular kernel. Proc. R. Soc. Ser. A 2021, 477, 20210438. [Google Scholar] [CrossRef]
  6. Caffarelli, L. Non-local diffusions, drifts and games. In Nonlinear Partial Differential Equations, Abel Symp., 7; Springer: Berlin/Heidelberg, Germany, 2012; pp. 37–52. [Google Scholar]
  7. Dubey, P.R.S.; Goswami, H.M.; Baskonus, T.; Gomati, A. On the Existence and Uniqueness Analysis of Fractional Blood Glucose-Insulin Minimal Model. Int. J. Model. Simul. Sci. Comput. 2022. [Google Scholar] [CrossRef]
  8. Kumar, A.; Prakash, A.; Baskonus, H.M. The epidemic COVID-19 model via Caputo-Fabrizio fractional operator. Waves Random Complex Media 2022. [Google Scholar] [CrossRef]
  9. Laskin, N. Fractional quantum mechanics and Lévy path integrals. Phys. Lett. A 2000, 268, 298–305. [Google Scholar] [CrossRef] [Green Version]
  10. Valdinoci, E. From the long jump random walk to the fractional Laplacian. Bol. Soc. Esp. Mat. Apl. SeMA 2009, 49, 33–44. [Google Scholar]
  11. 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]
  12. Fiscella, A.; Valdinoci, E. A critical Kirchhoff type problem involving a nonlocal operator. Nonlinear Anal. 2014, 94, 156–170. [Google Scholar] [CrossRef]
  13. Autuori, G.; Fiscella, A.; Pucci, P. Stationary Kirchhoff problems involving a fractional elliptic operator and a critical nonlinearity. Nonlinear Anal. 2015, 125, 699–714. [Google Scholar] [CrossRef]
  14. Bisci, G.M.; Vilasi, L. On a fractional degenerate Kirchhoff-type problem. Commun. Contemp. Math. 2017, 19, 1550088. [Google Scholar] [CrossRef]
  15. Pucci, P.; Xiang, M.; Zhang, B. Existence and multiplicity of entire solutions for fractional p-Kirchhoff equations. Adv. Nonlinear Anal. 2016, 5, 27–55. [Google Scholar] [CrossRef]
  16. Wang, F.; Hu, D.; Xiang, M. Combined effects of Choquard and singular nonlinearities in fractional Kirchhoff problems. Adv. Nonlinear Anal. 2021, 10, 636–658. [Google Scholar] [CrossRef]
  17. Lin, Q.; Tian, X.; Xu, R.; Zhang, M. 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] [Green Version]
  18. Xu, R.; Yang, Y.; Liu, Y. Global well-posedness for strongly damped viscoelastic wave equation. Appl. Anal. 2013, 92, 138–157. [Google Scholar] [CrossRef]
  19. Xiang, M.; Hu, D. Existence and blow-up of solutions for fractional wave equations of Kirchhoff type with viscoelasticity. Discret. Contin. Dyn. Syst. Ser. S 2021, 14, 4609–4629. [Google Scholar] [CrossRef]
  20. Liu, Y.; Lv, P.; Da, C. Blow-up of a nonlocal p-Laplacian evolution equation with critical initial energy. Ann. Polon. Math. 2016, 117, 89–99. [Google Scholar] [CrossRef]
  21. Haraux, A.; Zuazua, E. Decay estimates for some semilinear damped hyperbolic problems. Arch. Ration. Mech. Anal. 1988, 100, 191–206. [Google Scholar] [CrossRef]
  22. Zuazua, E. Exponential decay for the semilinear wave equation with locally distributed damping. Comm. Partial. Differ. Equ. 1990, 15, 205–235. [Google Scholar]
  23. Servadei, R.; Valdinoci, E. Mountain Pass solutions for non-local elliptic operators. J. Math. Anal. Appl. 2012, 389, 887–898. [Google Scholar] [CrossRef] [Green Version]
  24. Servadei, R.; Valdinoci, E. Variational methods for non-local operators of elliptic type. Discret. Contin. Dyn. Syst. 2013, 33, 2105–2137. [Google Scholar] [CrossRef]
  25. Servadei, R.; Valdinoci, E. The Brezis-Nirenberg result for the fractional Laplacian. Trans. Amer. Math. Soc. 2015, 367, 67–102. [Google Scholar] [CrossRef]
  26. Liu, Y. Long-time behavior of a class of viscoelastic plate equations. Elec. Res. Arch. 2020, 28, 549–567. [Google Scholar] [CrossRef] [Green Version]
  27. Sattinger, D.H. On global solution of nonlinear hyperbolic equations. Arch. Ration. Mech. Anal. 1968, 30, 148–172. [Google Scholar] [CrossRef]
  28. Cavalcanti, M.M.; Cavalcanti, V.N.D.; Martinez, P. Existence and decay rate estimates for the wave equation with nonlinear boundary damping and source term. J. Differ. Equ. 2004, 203, 119–158. [Google Scholar] [CrossRef] [Green Version]
  29. Esquivel-Avila, J.A. A characterization of global and nonglobal solutions of nonlinear wave and Kirchhoff equations. Nonlinear Anal. 2003, 52, 1111–1127. [Google Scholar] [CrossRef]
  30. Gazzola, F.; Squassina, M. Global solutions and finite time blow up for damped semilinear wave equations. Ann. Inst. H. Poincaré Anal. Non Linéaire 2006, 23, 185–207. [Google Scholar] [CrossRef]
  31. Lian, W.; Xu, R. Global well-posedness of nonlinear wave equation with weak and strong damping terms and logarithmic source term. Adv. Nonlinear Anal. 2020, 9, 613–632. [Google Scholar] [CrossRef]
  32. Lions, J.L. Quelques Méthodes de Résolution des Problèmes aux Limites Non Linéaires; Dunod: Paris, France, 1969. [Google Scholar]
  33. Liu, Y.; Li, W. 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]
  34. Liu, Y.; Mu, J.; Jiao, Y. A class of fourth order damped wave equations with arbitrary positive initial energy. Proc. Edinb. Math. Soc. 2019, 62, 165–178. [Google Scholar] [CrossRef]
  35. Liu, Y.; Yu, T.; Li, W. 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]
  36. Payne, L.E.; Sattinger, D.H. Sadle points and instability of nonlinear hyperbolic equations. Israel J. Math. 1975, 22, 273–303. [Google Scholar] [CrossRef]
  37. Tsutsumi, M. On solutions of semilinear differential equations in a Hilbert space. Math. Japon. 1972, 17, 173–193. [Google Scholar]
  38. Vitillaro, E. Global nonexistence theorems for a class of evolution equations with dissipation. Arch. Ration. Mech. Anal. 1999, 149, 155–182. [Google Scholar] [CrossRef]
  39. Xu, R.; Zhang, M.; Chen, S.; Yang, Y.; Shen, J. 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]
  40. Veeresha, P. Analysis of the spread of infectious diseases with the effects of consciousness programs by media using three fractional operators. In Methods of Mathematical Modelling: Infectious Diseases; Academic Press: Cambridge, MA, USA, 2022; pp. 113–135. [Google Scholar]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Liu, Y.; Zhang, L. Initial Boundary Value Problem for a Fractional Viscoelastic Equation of the Kirchhoff Type. Fractal Fract. 2022, 6, 581. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract6100581

AMA Style

Liu Y, Zhang L. Initial Boundary Value Problem for a Fractional Viscoelastic Equation of the Kirchhoff Type. Fractal and Fractional. 2022; 6(10):581. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract6100581

Chicago/Turabian Style

Liu, Yang, and Li Zhang. 2022. "Initial Boundary Value Problem for a Fractional Viscoelastic Equation of the Kirchhoff Type" Fractal and Fractional 6, no. 10: 581. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract6100581

Article Metrics

Back to TopTop