The Existence of Global Attractor for Kirchhoff-Type Strongly Damped Wave Equation with Nonlinear Memory

Abstract

This paper addresses the existence of a global attractor for Kirchhoff-type strongly damped wave equation with nonlinear memory effects. The key innovation of our work lies in reformulating the historical memory term as a convolution integral involving the memory kernel μ and a nonlinear power-law function | u( s ) | β u( s ) . First, we rigorously establish the existence, uniqueness, and regularity of solutions for Equation (1.2) through a systematic application of a priori energy estimates and the Faedo-Galerkin approximation method, while simultaneously demonstrating the presence of a bounded absorbing set. To analyze the asymptotic dynamics, we decompose the solution semigroup S( t ) into two components: S 1 ( t ) , governed by higher-order regularity, and S 2 ( t ) , capturing dissipative effects. The compactness of S 1 ( t ) is established via operator regularity analysis combined with the compact sobolev embedding theorem, while the uniform exponential decay of S 2 ( t ) is proven through refined energy estimation techniques. By synthesizing these results, we conclusively demonstrate the existence of a global attractor for the system under study, thereby extending the theoretical framework for nonlinear wave equations with memory-driven dissipation.

Share and Cite:

Xu, J. and Wang, X. (2025) The Existence of Global Attractor for Kirchhoff-Type Strongly Damped Wave Equation with Nonlinear Memory. Journal of Applied Mathematics and Physics, 13, 1212-1231. doi: 10.4236/jamp.2025.134064.

1. Introduction

The wave equation, as a cornerstone of continuum mechanics, plays a pivotal role in modeling dissipative wave propagation phenomena across diverse physical systems. These include but are not limited to: acoustic attenuation in viscoelastic media (e.g., polymer melts), stress relaxation in aging composite materials, and energy dissipation in nonlinear metamaterials. The mathematical formulation of dissipative wave equations traces back to foundational work by Hale et al. [1], with the general form:

{ u tt +g( u t )Δu+f( u )=h( x ), xΩ,t>0, u( x,0 )= u 0 ( x ), u t ( x,0 )= u 1 ( x ), xΩ, u| Ω =0. (1.1)

where the dissipative term g( u t ) governs energy loss mechanisms. Physically, the strong damping term g( u t )=Δ u t (as studied in [2]-[4]) models highly viscous systems where the damping force depends on the Laplacian of the velocity field—a feature critical for simulating wave attenuation in frequency-dependent media such as biological tissues or rheological fluids.

Recent advances in materials science have highlighted the necessity of incorporating memory effects to describe time-history-dependent responses. For instance, in viscoelastic polymers subjected to cyclic loading, the stress-strain relationship depends nonlinearly on both current deformations and accumulated microstructural damage—a phenomenon that cannot be captured by linear memory models. This has motivated studies on memory-driven wave equations of the form:

u tt Δuσ( u 2 )Δ u t +ρ | u( t ) | γ u( t ) 0 μ ( ts ) | u( s ) | β u( s )ds+f( u )=h( x )

where the nonlinear memory term 0 μ ( ts ) | u( s ) | β u( s )ds introduces a power-law dependence on historical states β>0 , enabling the modeling of fatigue damage accumulation in alloys or creep behavior in geomaterials. In particular, wave equations with linear memory terms have emerged as powerful tools for characterizing the time-dependent evolution of viscoelastic materials, where the current state depends intrinsically on its historical configurationsa framework rigorously explored in [5] [6], see related literature [7]-[9]. Such equations typically incorporate nonlocal operators expressed as convolution integrals between the unknown function u( t ) and a prescribed memory kernel μ . The kernel, often assumed to decay monotonically at infinity, encodes the systems fading memory properties and governs its asymptotic dissipation dynamics. These memory-driven models hold profound implications across multiple disciplines, including solid mechanics, fluid dynamics, and heat transfer theory, where they provide critical insights into energy dissipation mechanisms and long-term stability. Meanwhile, recent advances in dynamical systems theory and global attractor analysis have spurred significant interest in mathematical models incorporating memory effects, as evidenced by foundational studies in [2] [6] [10].

In their seminal 2003 work [11], Yang Meihua and Sun Chunyou established the existence of a global attractor for a strongly damped wave equation with critically growing nonlinearity in the phase space H 0 1 ( Ω )× H 0 1 ( Ω ) . They further demonstrated that this attractor is bounded in the refined space ( H 2 ( Ω ) H 0 1 ( Ω ) )×( H 2 ( Ω ) H 0 1 ( Ω ) ) , exhibiting attraction properties for bounded sets H 0 1 ( Ω )× L 2 ( Ω ) under its natural norm. Subsequent advances include Pata’s 2008 study [4], which derived conditions for the existence of global attractors in strongly damped wave equations with memory, requiring the memory kernel μ( s ) to satisfy the exponential decay criterion μ( s+σ )K e δσ μ( s ) ( K1,δ>0 ). Further contributions emerged in 2012 with Han Yinghao’s analysis of nonlinear memory effects in damped wave equations [12], followed by D’Abbicco’s 2014 investigations into weakly damped and structurally damped wave equations with nonlinear memory in [3] [13], respectively. Parallel developments in Kirchhoff-type damped wave equations have also been extensively explored, as evidenced by works such as [14]-[17]. Despite these advancements, the literature remains sparse regarding systems combining nonlinear memory with Kirchhoff-type damping. Motivated by the methodologies in [3] [12]-[14], this paper addresses this gap by investigating the existence of global attractors for a Kirchhoff-type strongly damped wave equation with nonlinear memory. The governing equation under consideration is formulated as follows:

{ u tt Δuσ( u 2 )Δ u t +ρ | u( t ) | γ u( t ) 0 μ ( ts ) | u( s ) | β u( s )ds+f( u )=h( x ), x Ω,t>0, u( x,t )=0, x Ω,t0, u( x,t )= u 0 ( x ), u t ( x,t )= u 1 ( x ), x Ω,t0 (1.2)

where Ω N be a bounded domain with sufficiently smooth boundary Ω , and γ and β be positive constants, σ( u 2 )Δ u t is a Kirchhoff-type strong damping. The damping coefficient σ( ) C 1 ( + ) monotonically decreases, and for all s + there exists σ( s )>0 . Moreover, h( ) L 2 ( Ω ) .

This paper is organized as follows. Section 2 presents fundamental definitions and key theorems in dynamical systems theory. In Section 3, we rigorously establish the existence and uniqueness of solutions for system (1.2), thereby demonstrating that the system generates a well-defined dynamical system. Finally, Section 4 constructs a compact and connected global attractor for the dynamical system associated with Equation (1.2).

2. Preliminaries

In this Section, we introduce some notations and the useful theorem, which is the key technique to establish the well-posedness and existence of global attractors, which are from [1] [16] [18] [19].

Let L p ( Ω ) , 1p is the usual Lebesgue spaces. The norm and inner-product of L 2 ( Ω ) will be denoted by and ( , ) , respectively. Let C be an arbitrary positive constant, which may be different from line to line and even in the same line. We set H= V 0 = L 2 ( Ω ) , V= V 1 = H 0 1 ( Ω ) . Let V s =D( A s 2 ) , and A=Δ is a positive definite operator on H with discrete spectrum, with

A e i = λ i e i ,0< λ 1 λ 2 , lim k λ k =.

Therefore, { e k } forms a normalized orthogonal basis of space H and simultaneously serves as a normalized orthogonal basis in V . The inner product and norm for space V s are defined as follow

( u,v ) V s =( A s 2 u, A s 2 v ), u V s = A s u .

Using H and V respectively represent the dual spaces of H and V , and in this paper, we use C to represent different constants. λ 1 >0 is the first eigenvalue of the operator Δ .

We assume that the memory kernel μ is subject to the following hypotheses:

( h 1 ) μ C 1 ( + ) L 1 ( + ) , + =( 0, ) ;

( h 2 ) μ( s )0 , μ ( s )0 , s + ;

( h 3 ) 0 μ ( s )ds=1 ;

( h 4 ) m 1 μ( s ) μ ( s ) m 2 μ( s ) , s + , where m 1 , m 2 >0 .

In addition, the memory kernel μ satisfying assumptions ( h 1 ) - ( h 4 ) exists, such as μ( s )=m e ms .

We assume that the nonlinearity fC( ) and f( 0 )=0 , is subject to the following hypotheses: there exists positive constants ν0 , γβ and

ρ> γ+2 β+2 | Ω | 1 γ 1 β such that

( f 1 ) b s 2 A( x )F( s )sf( s )+θ ;

( f 2 ) | f ( s ) | c 1 ( 1+ | s | ν ) , ν0 ,

where 0<b< λ 1 2 , A( ) L 1 ( Ω ) , and A( )θ , F( s )= 0 s f ( y )dy .

Definition 2.1. Let 0 be an open subset of a matric space X , and S( t ) be a semigroup on X . We say that 0 is an absorbing set of S( t ) if the orbit of any bounded set enters into 0 after a certain time (which may depend on the set X ). Namely, for any bounded set , there exists t 1 ( ) , such that S( t ) 0 ,t t 1 ( ).

Definition 2.2. Let X be a metric space. The Kuratowski measure of a noncompact set in X is defined by

K X ( )=inf{ d:hasafinitecoverofopenballsofXofdiameterlessthand }

Theorem 2.1. [18] Let X be a Banach space, and let S( t )( t0 ) be a continuous operator semigroup on X . Suppose there exist an open set UX and a bounded set U that is absorbing in U , and one of the following conditions holds:

i) Uniform Compactness for Large t : For every bounded set , there exists t= t 0 ( ) such that

t t 0 S( t )

is relatively compact in X .

ii) Decomposition into Compact and Decaying Parts: S( t )= S 1 ( t )+ S 2 ( t ) , where

S 1 ( t ) : is uniformly compact for sufficiently large t ,

S 2 ( t ) : XX is a continuous mapping, and for every bounded set X ,

sup φ S 2 ( t )φ X 0.

Then, the ω -limit set A=ω( ) is a compact attractor that attracts all bounded sets in U , is the maximal bounded attractor in U , and is connected if U is both convex and connected.

Theorem 2.2. [18] Let X be a complete metric space, and S( t ) be a continuous semigroup on X . Suppose S( t ) admits a bounded absorbing set UX and satisfies that for every bounded subset X , the Kuratowski measure of noncompactness α( S( t ) )0 as t . Then, the ω -limit set A=ω( U ) is a compact maximal attractor for S( t ) .

3. Well-Posedness of Solutions

We have the existence and uniqueness of solutions which are obtained by the Faedo-Galerkin approximation and a compactness method.

Theorem 3.1. If ( h 1 ) - ( h 4 ) and ( f 1 ) - ( f 2 ) hold, and γ,β,ν satisfy the following conditions

{ 0ν,β,γ<,βγ, whenn=1,2, 0ν,β,γ<2,βγ, whenn=3, ν=β=γ=0, whenn4,

then for any u 0 V, u 1 H and hH , Equation (1.2) has a unique global solution u such that

uC( [ 0, );V ), u t C( [ 0, );H ), u tt C( [ 0, ); V ).

Proof. To prove the existence of solutions, we use the standard Galerkin method. Let V m =span{ e 1 , e 2 , e 3 ,, e m } , be the subspace of H . We consider

u m = i=1 m a im ( t ) e i ,t0, (3.1)

which is the approximate solution of Faedo-Galerkin of order m , that is, we get a system of ODEs in the variables a im ( t ) ( i=1,2,,m ) of the form

( tt u m , e i )( Δ u m , e i )( σ( u m 2 ) t Δ u m , e i )+( ρ | u m ( t ) | γ u m ( t ), e i ) ( 0 μ ( ts ) | u m ( s ) | β u m ( s )ds, e i )+( f( u m ), e i )=( h( x ), e i ), (3.2)

subject to the initial conditions

u m ( 0 )= u m0 = i=1 m a im ( 0 ) e i u 0 , (3.3)

t u m ( 0 )= u m1 = t i=1 n a im ( 0 ) e i u 1 . (3.4)

It is well-known that the above finite-dimensional system of ordinary functional differential equations is well-posed at least locally (see for example [19]). Indeed, for fixed m, the system (3.2) defines a linear system of differential equations on Ω N . Then we can apply differential equations theory for local existence and uniqueness of solutions to the system (3.2). Hence, the system of (3.2) possesses a unique local solution a im ( t ) defined in [ 0,T ] , with 0<T . Also, let us prove that a priori estimate for the Faedo-Galerkin approximate solutions u m ( t ) . Multiplying Equation of (3.2) by t a im and ε a im , respectively, then summing over i and adding the results, we obtain

( tt u m , t u m +ε u m )( Δ u m , t u m +ε u m )( σ( u m 2 ) t Δ u m , t u m +ε u m ) +( ρ | u m ( t ) | γ u m ( t ), t u m +ε u m )( 0 μ ( ts ) | u m ( s ) | β u m ( s )ds, t u m +ε u m ) +( f( u m ), t u m +ε u m )=( h( x ), t u m +ε u m ). (3.5)

Define also

( μ u m )( t )= 0 μ ( ts ) | | u m ( s ) | β 2 u m ( s ) | u m ( s ) | β 2 u m ( t ) | 2 ds; (3.6)

( μ u m )( t )= 0 μ ( ts ) | | u m ( s ) | β 2 u m ( s ) | u m ( s ) | β 2 u m ( t ) | 2 ds; (3.7)

1 2 d dt ( 0 μ ( ts ) u m ( s ) β+2 β+2 ds ) = 1 2 μ( 0 ) u m ( t ) β+2 β+2 + 1 2 0 μ ( ts ) u m ( s ) β+2 β+2 ds. (3.8)

Combining (3.6) - (3.8), we get

( 0 μ ( ts ) | u m ( s ) | β u m ( s )ds, t u m ) = 1 2 ( μ u m )( t ) 1 2 d dt ( μ u m )( t ) 1 2 μ( 0 ) u m ( t ) β+2 β+2 + 1 2 d dt ( 0 μ ( ts ) | u m ( s ) β 2 u m ( t ) | 2 ds ) 1 2 0 μ ( ts ) | u m ( s ) β 2 u m ( t ) | 2 ds. (3.9)

Let v= t u m +ε u m , we substitute (3.8) and (3.9) in (3.5), by direct calculations, one gets

( h,v )= 1 2 d dt [ v 2 +( 1+εσ( u m 2 ) ) u m 2 + 2ρ γ+2 u m ( t ) γ+2 +( μ u m )( t ) 0 μ ( ts ) | u m ( s ) β 2 u m ( t ) | 2 ds+ 0 μ ( ts ) u m ( s ) β+2 β+2 ds +2 Ω ( F( u m )+A( x ) )dx ]ε v 2 + ε 2 ( u m ,v )+ε u m 2 +σ( u m 2 ) t u m 2 1 2 u m 2 ε d dt [ σ( u m 2 ) ] +ερ u m ( t ) γ+2 1 2 ( μ u m )( t ) 1 2 0 μ ( ts ) u m ( s ) β+2 β+2 ds + 1 2 0 μ ( ts ) | u m ( s ) β 2 u m ( t ) | 2 ds+ε Ω f ( u m ) u m dx ε( 0 μ ( ts ) | u m ( s ) | β u m ( s )ds, u m ( t ) ). (3.10)

By Hölder inequality and Young’s inequality, we obtain

( 0 μ ( ts ) | u m ( s ) | β u m ( s )ds, u m ( t ) ) = Ω 0 μ ( ts ) | u m ( s ) | β u m ( s ) u m ( t )dsdx 0 μ ( ts ) u m ( s ) β 2 +1 | u m ( s ) | β 2 u m ( t )ds 1 2 0 μ ( ts ) u m ( s ) β+2 β+2 ds+ 1 2 0 μ ( ts ) | u m ( s ) β 2 u m ( t ) | 2 ds. (3.11)

On the one hand, since σ( ) monotonically decreases and strictly positive in + , there exists a constant σ 0 >0 such that t0 , one gets

σ( u 2 ) σ 0 , (3.12)

and d dt σ( u 2 )<0 . On the other hand, by the following inequality

v 2 = t u m +ε u m 2 2( t u m 2 + ε 2 u m 2 ),

we have

ε v 2 + ε 2 ( u m ,v )+ε u m 2 +σ( u m 2 ) t u m 2 1 2 u m 2 ε d dt [ σ( u m 2 ) ] ε v 2 +ε u m 2 + σ 0 λ 1 ( 1 2 v 2 ε 2 λ 1 u m 2 ) ε 2 λ 1 ( λ 1 v 2 + u m 2 4 λ 1 ) [ σ 0 λ 1 2 ε( 1+ε ) ] v 2 +ε( 1ε( σ 0 + 1 4 λ 1 ) ) u m 2 . (3.13)

Meanwhile, by hypothesis ( f 1 ) , one has

ε Ω f( u m ) u m ε Ω ( F( u m )θ )dx =ε Ω ( F( u m )+A( x ) )dx ε Ω ( A( x )+θ )dx . (3.14)

Then, taking into account hypothesis ( f 1 ) and (3.11) - (3.14) in (3.10), we obtain

( h,v ) 1 2 d dt [ v 2 +( 1+ε σ 0 ) u m 2 + 2ρ γ+2 u m ( t ) γ+2 +( μ u m )( t ) 0 μ ( ts ) | u m ( s ) β 2 u m ( t ) | 2 ds+ 0 μ ( ts ) u m ( s ) β+2 β+2 ds +2 Ω ( F( u m )+A( x ) )dx ]+[ σ 0 λ 1 2 ε( 1+ε ) ] v 2 +ε( 1ε( σ 0 + 1 4 λ 1 ) ) u m 2 +ερ u m ( t ) γ+2 + m 1 2 ( μ u m )( t ) m 2 +ε 2 0 μ ( ts ) | u m ( s ) β 2 u m ( t ) | 2 ds + m 1 ε 2 0 μ ( ts ) u m ( s ) β+2 β+2 ds+ε Ω ( F( u m )+A( x ) )dx ε Ω ( A( x )+θ )dx . (3.15)

Besides, since

( h,v ) h v 1 σ 0 λ 1 h 2 + σ 0 λ 1 4 v 2 ,

from (3.15), we get

1 2 d dt [ v 2 +( 1+ε σ 0 ) u m 2 + 2ρ γ+2 u m ( t ) γ+2 +( μ u m )( t ) 0 μ ( ts ) | u m ( s ) β 2 u m ( t ) | 2 ds+ 0 μ ( ts ) u m ( s ) β+2 β+2 ds +2 Ω F ( u m )+A( x )dx ]+[ σ 0 λ 1 4 ε( 1+ε ) ] v 2 +ε[ 1ε( σ 0 + 1 4 λ 1 ) ] u m 2 +ερ u m ( t ) γ+2 + m 1 2 ( μ u m )( t ) m 2 +ε 2 0 μ ( ts ) | u m ( s ) β 2 u m ( t ) | 2 ds + m 1 ε 2 0 μ ( ts ) u m ( s ) β+2 β+2 ds+ε Ω ( F( u m )+A( x ) )dx ε Ω ( A( x )+θ )dx + 1 σ 0 λ 1 h 2 . (3.16)

Choosing ε small enough such that σ 0 λ 1 4 ε( 1+ε ) , ε( 1ε( σ 0 + 1 4 λ 1 ) ) , m 1 ε 2 , are positive constants, and let

δ=2min{ σ 0 λ 1 4 ε( 1+ε ), ε( 1ε( σ 0 + 1 4 λ 1 ) ) 1+ε σ 0 , ε 2 , ε( γ+2 ) 2 , m 1 2 , m 1 ε 2 , m 2 +ε 2 },

thus, (3.16) becomes

d dt E m ( t )+δ E m ( t )2ε Ω ( A( x )+θ )dx + 2 σ 0 λ 1 h 2 , (3.17)

where

E m ( t )= v 2 +( 1+ε σ 0 ) u m 2 + 2ρ γ+2 u m ( t ) γ+2 +( μ u m )( t ) 0 μ ( ts ) | u m ( s ) β 2 u m ( t ) | 2 ds+ 0 μ ( ts ) u m ( s ) β+2 β+2 ds +2 Ω ( F( u m )+A( x ) )dx .

Using Young’s inequality and hypothesis ( h 3 ) , we get that

0 μ ( ts ) | | u m ( s ) | β 2 u m ( t ) | 2 ds β β+2 0 μ ( ts ) u m ( s ) β+2 β+2 ds+ 2 β+2 0 μ ( ts ) u m ( t ) β+2 β+2 ds β β+2 0 μ ( ts ) u m ( s ) β+2 β+2 ds+ 2 β+2 | Ω | 1 γ 1 β u m ( t ) β+2 . (3.18)

Because of ρ> γ+2 β+2 | Ω | 1 γ 1 β , putting (3.18) in E m ( t ) , and by means of the

hypothesis ( f 1 ) , it immediately gives that

E m ( t ) v 2 +( 1+ε σ 0 ) u m 2 +( μ u m )( t )+2 Ω ( F( u m )+A( x ) )dx + 2 β+2 0 μ ( ts ) u m ( s ) β+2 β+2 ds+( 2ρ γ+2 2 β+2 | Ω | 1 γ 1 β ) u m ( t ) β+2 v 2 +( 1+ε σ 0 ) u m 2 +( μ u m )( t )+2 Ω ( F( u m )+A( x ) )dx v 2 +( 1+ε σ 0 ) u m 2 +( μ u m )( t ) 2b λ 1 u m 2 v 2 +ε σ 0 u m 2 .

Applying Gronwall’s inequality in (3.17), we deduce that

E m ( t ) E m ( 0 ) e δt +( 2 Ω ( A( x )+θ )dx + 2 ε σ 0 λ 1 h 2 )( 1 e δt ). (3.19)

According to the hypothesis ( f 2 ) , when u( 0 ) and u t ( 0 ) are bounded, Ω F ( u( 0 ) )dx is also bounded, it is obvious that E m ( 0 ) is bounded, therefore, we claim that

R= sup ( u 0 , u 1 )B E m ( 0 )<,

if we set μ 0 2 =2 Ω ( A( x )+θ )dx + 2 ε σ 0 λ 1 h 2 , by using (3.19), for all t[ 0, ) , we have that

t u m +ε u m 2 +ε σ 0 u m 2 = v 2 +ε σ 0 u m 2 E m ( t )R e δt + μ 0 2 . (3.20)

Consequently, we can get from the above estimations that

u m isboundeduniformlyin L ( [ 0, ); H 0 1 );

t u m isboundeduniformlyin L ( [ 0, ); L 2 ).

And because of the weak compactness, we can conclude that there exists a subsequence of solutions of the Galerkin approximations, denoted also by { u m } m , such that, for some u ,

u m * uweaklystarin L ( [ 0, ); H 0 1 );

t u m * t uweaklystarin L ( [ 0, ); L 2 ).

Meanwhile, since H 0 1 ↪↪ L 2γ+2 L 2 , by Aubin-Lions theorem, u m u is convergent strongly in L 2 ( [ 0,T ]; L 2 ) . Therefore, u m u is true almost everywhere in Ω×[ 0,T ] .

In the following, we will verify the convergence of the term ρ | u m | γ u m . Since u m u is almost everywhere in Ω×[ 0,T ] , and | x | a x( a>0 ) is a continuous function, we have

ρ | u m | γ u m ρ | u | γ uisconvergentalmosteverywhereinΩ×[ 0,T ].

By the embedding H 0 1 L 6 ( Ω ) L 2γ+2 ( Ω ) and (3.20), we obtain

0 Ω | ρ | u m | γ u m |dxdt ρ 0 T u m 2γ+2 2γ+2 dt C. (3.21)

Using (3.20) and the Lions lemma we also have that

ρ | u m | γ u m ρ | u | γ uin L 2 ( Ω×[ 0,T ] ).

Next we establish the convergence of the nonlinear term f . It follows from u m u is almost everywhere in Ω×[ 0,T ] and from f( u m ) is continuous, one deduces that

f( u m )f( u )isconvergentalmosteverywhereinΩ×[ 0,T ],

we apply the hypothesis ( f 2 ) and the Differential mean value theorem, there exists a constant θ( 0,1 ) , such that

f( u m ) 2 = f( u m )f( 0 ) 2 = f ( θ u m ) u m 2 Ω | C( 1+ | u m | ν )| u m | | 2 dx C Ω ( | u m | 2 + | u m | 2ν+2 )dx C( u m 2 + u m 2ν+2 2ν+2 ).

Due to the embedding H 0 1 L 6 ( Ω ) L 2ν+2 ( Ω ) , one implies that

0 t f( u m ) 2 dt C( u m 2 + u m 2ν+2 2ν+2 )C. (3.22)

Moreover, by the Lions lemma, it holds that

f( u m )f( u )in L 2 ( Ω×[ 0,T ] ).

By the above convergence results, integrating (3.2) over [ 0,T ] and taking the limit m , it is clear that we obtain the solutions of (1.2) exist in L ( 0,T; V ) .

Now, we will prove the continuous dependence of the solutions on the initial data. To do this, we consider u and v are two solutions of (1.2) corresponding to the initial conditions ( u 0 , u 1 ) and ( v 0 , v 1 ) respectively. Set w=uv , therefore, it has w( 0 )= u 0 v 0 , w t ( 0 )= u 1 v 1 . Substituting u and v into (1.2) respectively, then subtracting the two formulas, we readily obtain that

w tt Δw( σ( u 2 )Δ u t σ( v 2 )Δ v t )+( ρ | u( t ) | γ u( t )ρ | v( t ) | γ v( t ) ) ( 0 μ ( ts ) | u( s ) | β u( s )ds 0 μ ( ts ) | v( s ) | β v( s )ds )+( f( u )f( v ) )=0. (3.23)

Multiplying this Equation (3.23) by w t and integrating over Ω , we derive that

1 2 d dt w t 2 + 1 2 d dt w 2 ( σ( u 2 )Δ u t σ( v 2 )Δ v t , w t ) =( ρ | v( t ) | γ v( t )ρ | u( t ) | γ u( t ), w t )+( 0 μ ( ts ) | u( s ) | β u( s )ds 0 μ ( ts ) | v( s ) | β v( s )ds, w t )+( f( v )f( u ), w t ). (3.24)

In the following, we will estimate the third item on the left and the items on the right of (3.24), respectively, where

( σ( u 2 )Δ u t σ( v 2 )Δ v t , w t ) =( σ( u 2 ) u t σ( v 2 ) v t , w t ) σ 0 ( u t v t , w t )= σ 0 w t 2 σ 0 λ 1 . (3.25)

By applying the Mazur inequality, the Young’s inequality, the Hölder inequality, and the embedding H 0 1 L 6 ( Ω ) L γ+2 ( Ω ) , we get that

( ρ | v( t ) | γ v( t )ρ | u( t ) | γ u( t ), w t ) ρ Ω | | v( t ) | γ v( t ) | u( t ) | γ u( t ) || w t |dx ρ( γ+1 ) Ω ( | v( t ) | γ + | u( t ) | γ )| w || w t |dx ρ( γ+1 )( v( t ) γ + u( t ) γ ) w w t Cρ( γ+1 )( v( t ) γ + u( t ) γ ) w w t Cρ( γ+1 )( w 2 + w t 2 ). (3.26)

Similarly, we conclude that

( 0 μ ( ts ) | u( s ) | β u( s )ds 0 μ ( ts ) | v( s ) | β v( s )ds, w t ) ( β+1 ) Ω 0 μ ( ts )| | u( s ) | β + | v( s ) | β | | w | | w t |dsdx C( β+1 ) | u( s ) | β + | v( s ) | β β+2 w w t C( β+1 )( u( s ) β + v( s ) β ) w w t C( β+1 )( w 2 + w t 2 ). (3.27)

From the hypothesis ( f 2 ) and the Differential mean value theorem, combining the Hölder inequality, the Young inequality and the embedding H 0 1 L 6 ( Ω ) L ν+2 ( Ω ) , yield that

( f( u )f( v ), w t ) =( f( θ u+( 1θ )v )( w ), w t ) C Ω ( 1+ | u | ν + | v | ν )| w || w t |dx C( 1+ u ν+2 ν + v ν+2 ν ) w w t C( 1+ u ν+2 ν + v ν+2 ν )( 1 2 w 2 + 1 2 w t 2 ) C( w 2 + w t 2 ). (3.28)

Plugging (2.25) - (2.28) into (2.24), and combining the Poincaré inequality, we show that

C d dt ( w t 2 + w 2 )C[ ρ( γ+1 )+( β+1 )+1 ] w 2 +[ C( ρ( γ+1 )+( β+1 )+1 ) σ 0 λ 1 ] w t 2 , (3.29)

where C min{ 1 2 , λ 1 2 } . Consequently, thanks to the Gronwall lemma, we deduce

w t 2 + w 2 C( w t ( 0 ) 2 + w( 0 ) 2 ),t0. (3.30)

hence, the continuous dependence of the solutions on the initial data is turned out. Meanwhile, the formula (3.30) suggests that the solutions of (1.2) are unique.

4. Existence of Global Attractor

By Theorem 3.1 problem (1.2) generates an evolution semigroup S( t ) in the space V×H by the formula S( t )y=( u( t ); u t ( t ) ) where y=( u 0 ; u 1 )V×H and u( t ) is a weak solution to (1.2).

In this Section, we first show the existence of a bounded absorbing set for S( t ) in V×H , then it follows from Theorem 2.1 that we prove the ω -limit set of 0 is the connected and compact global attractor of S( t ) .

As in the proof of the priori estimates of Theorem 3.1 using the formula (3.20) one can see that there exists a constant t 0 = t 0 ( ) such that for all t t 0

t u( t 0 ) 2 + u( t 0 ) 2 R e δt + μ 0 2 .

Let 0 = t0 S( t ) 0 , where

0 ={ ( u 0 , u 1 )V×H: u 0 2 + u 1 2 2 μ 0 2 },

this implies that 0 is the bounded absorbing set of semigroup S( t ) in V×H . Therefore, the bounded dissipativity of the semigroup S( t ) can be obtained as follows:

Theorem 4.1. Let the forcing term hH . Assume that the nonlinear term f( ) C 1 ( ) satisfies conditions ( f 1 ) - ( f 2 ) , and further assume that conditions ( h 1 ) - ( h 4 ) hold. Let ( u, u t ) denote the solution to equation (1.2). corresponding to the initial data ( u 0 , u 1 ) . Then there exists a constant R 0 >0 such that the closed ball 0 =B( 0, R 0 ) is a bounded absorbing set in V×H for the solution semigroup S( t ) generated by Problem (1.2). Specifically, for every bounded subset V×H , there exists a critical time t 0 >0 ensuring the inclusion S( t ) 0 holds for all t t 0 .

Next, we will decompose the semigroup S( t ) to demonstrate its asymptotic compactness using the Kuratowski measure and α noncompactness measure. Specifically, the semigroup S( t ) can be decomposed into two components S 1 ( t ) and S 2 ( t ) , where:

S 1 ( t ) : Uniformly compact in V×H ,

S 2 ( t ) : Exponential decay.

Consequently, for any bounded set V×H , as t ,

α( S( t ) )α( S 1 ( t ) )+α( S 2 ( t ) )=α( S 2 ( t ) )0.

Decompose u into u ¯ + u ˜ , where u ¯ is the unique solution to the following problem:

{ u ¯ tt Δ u ¯ σ( u 2 )Δ u ¯ t = 0 μ ( ts ) | u( s ) | β u( s )dsρ | u( t ) | γ u( t )f( u )+h( x ), x Ω,t>0, u ¯ ( x,t )=0, x Ω,t0, u ¯ ( x,0 )=0, u ¯ t ( x,0 )=0, x Ω,t0 (4.1)

and, u ˜ is the unique solution to the complementary problem:

{ u ˜ tt Δ u ˜ σ( u 2 )Δ u ˜ t =0, xΩ,t>0, u ˜ ( x,t )=0, xΩ,t0, u ˜ ( x,0 )= u 0 , u ˜ t ( x,0 )= u 1 , xΩ,t0. (4.2)

Lemma 4.1. Define the abstract equation:

{ u tt +Au+σ( u 2 )A u t =h( x ), xΩ,t>0, u( x,t )=0, xΩ,t0, u( x,0 )= u 0 , u t ( x,0 )= u 1 , xΩ,t0. (4.3)

If hH , ( u 0 , u 1 )V×H , u is the solution to the equation (4.3), then ( u, u t ) C b ( + ;V×H ) such that the solution u satisfies the relation

u t 2 + u 2 ( u 1 2 + u 0 2 ) e α 1 t + C α 1 ( 1 e α 1 t ),

where α 1 >0 .

Proof. Taking the inner product in H of both sides of equation (4.3) with u t , we derive

1 2 d dt ( u t 2 + u 2 )+σ( u 2 ) u t 2 =( h, u t ),

since σ( ) C 1 ( + ) is monotonically decreasing and satisfies σ( s )>0 for all s + , then there exists a constant α 1 >0 such that σ( u 2 )> α 1 >0 . By applying the Cauchy-Schwarz inequality and Young’s inequality, we deduce that

( h, u t ) h u t 1 2 α 1 h 2 + α 1 2 u t . (4.4)

Consequently, we derive the following differential inequality:

d dt ( u t 2 + u 2 )+ α 1 ( u t 2 + u 2 )C, (4.5)

by applying Gronwall’s inequality, we obtain the following estimate:

u t 2 + u 2 ( u 1 2 + u 0 2 ) e α 1 t + C α 1 ( 1 e α 1 t ). (4.6)

We now proceed to demonstrate the uniform compactness of the semigroup S 1 ( t ) in the product space V×H .

Lemma 4.2. Under the assumptions of Theorem 3.1 on parameters γ , β , ν , there exists σ >0 such that for any ( u, u t ) C b ( + ,V×H ) , the following regularity holds:

t ( x,t ) L ( + , V σ 1 ). (4.7)

Furthermore, for any bounded set in V×H , there exists a constant C( )>0 such that whenever ( u 0 , u 1 ) , the following uniform bound is satisfied:

sup uV t ( x,t ) L ( + , V σ 1 ) C( ), (4.8)

where

( x,t )= 0 μ ( ts ) | u( s ) | β u( s )dsρ | u( t ) | γ u( t )f( u )+h( x ). (4.9)

Proof. By direct computation, we obtain the following expression for the time derivative of ( x,t ) :

t ( x,t )= 0 μ ( ts ) | u( s ) | β u( s )ds+μ( 0 ) | u( t ) | β u( t ) ρ( γ+1 ) | u( t ) | γ u t ( t ) f ( u ) u t . (4.10)

Under the assumptions of Theorem 3.1, we analyze the problem by considering three cases for the dimension n .

Case 1: For n=1,2 , we choose 1 σ >0 such that the following inequality holds:

nτ 2τ+4 <1 σ < n 2 , (4.11)

where τ=max{ γ,ν } . Consequently, we obtain

β+2γ+2τ+2 2n n2( 1 σ ) ,

which implies the following chain of continuous embeddings:

L 2n n2( 1 σ ) ( Ω ) L τ+2 ( Ω ) L γ+2 ( Ω ) L β+2 ( Ω ). (4.12)

Moreover, by the Sobolev embedding theorem and interpolation theory, we have

V 1 σ H 1 σ ( Ω ) L 2n n2( 1 σ ) ( Ω ). (4.13)

Since all embeddings are continuous, it follows from (4.12) and (4.13) that

V 1 σ L τ+2 ( Ω ) L γ+2 ( Ω ) L β+2 ( Ω ). (4.14)

Clearly, for any 1q< , H 0 1 ( Ω ) L q ( Ω ) . If we set q= nν 1 σ , then for all

u H 0 1 ( Ω ) with u M , the following estimate holds:

u( t ) γn 1 σ C u( t ) , (4.15)

finally, by assumption ( f 2 ) , we conclude that

f ( u ) L q ν C( M ). (4.16)

Case 2: For n=3 , the Sobolev embedding theorem yields

H 0 1 ( Ω ) L 6 ( Ω ); H s ( Ω ) L 6 32s ( Ω );0s< 3 2 , (4.17)

choosing 1 σ = τ 2 , we have

0<2( 1 σ )=τ<2< 6 2 , (4.18)

and

τ+2 2n n2( 1 σ ) . (4.19)

Thus, we recover (4.13) and (4.14). Since

γn 1 σ = 2γn τ <2n, (4.20)

we can deduce (4.15). This completes the proof that (4.13) - (4.14) hold for 1n3 .

First, let w( t ) V 1 σ . Taking the inner product of w( t ) with (4.10) in L 2 ( Ω ) , we obtain

( t ( x,t ),w( t ) )=( 0 μ ( ts ) | u( s ) | β u( s )ds,w( t ) ) +( μ( 0 ) | u( t ) | β u( t ),w( t ) ) ( ρ( γ+1 ) | u( t ) | γ u t ( t ),w( t ) )( f ( u ) u t ,w( t ) ). (4.21)

We now estimate each term on the right-hand side of (4.21). By Hölder inequality and Young’s inequality, we get

( μ( 0 ) | u( t ) | β u( t ),w( t ) ) μ( 0 ) ( Ω ( | u( t ) | β+1 ) β+2 β+1 dx ) β+1 β+2 w( t ) β+2 μ( 0 ) u( t ) β+2 β+1 w( t ) β+2 μ( 0 )( β+1 ) β+2 u( t ) β+2 β+2 w( t ) β+2 + μ( 0 ) β+2 w( t ) β+2 . (4.22)

Similarly, applying Hölder inequality, Young’s inequality, and ( h 4 ) , we derive

( 0 μ ( ts ) | u( s ) | β u( s )ds,w( t ) ) m 1 0 μ ( ts ) u( s ) β+2 β+1 w( t ) β+2 ds m 1 ( β+1 ) β+2 0 μ ( ts ) u( s ) β+2 β+2 w( t ) β+2 ds+ m 1 β+2 w( t ) β+2 . (4.23)

Using Hölder inequality ( 2( 1 σ ) 2n + n 2n + n2( 1 σ ) 2n =1 ) and (4.15), we obtain

( ρ( γ+1 ) | u( t ) | γ u t ( t ),w( t ) ) ρ( γ+1 ) u( t ) γn/ ( 1 σ ) γ u t ( t ) 2 w( t ) 2n/ ( n2( 1 σ ) ) Cρ( γ+1 ) u( t ) γ u t ( t ) w( t ) 2n/ ( n2( 1 σ ) ) . (4.24)

Similarly, using Hölder inequality ( 2( 1 σ ) 2n + n 2n + n2( 1 σ ) 2n =1 ) and (4.16), we derive

( f ( u ) u t ,w( t ) ) f ( u ) n/ ( 1 σ ) u t ( t ) 2 w( t ) 2n/ ( n2( 1 σ ) ) C( M ) u t ( t ) w( t ) 2n/ ( n2( 1 σ ) ) . (4.25)

Substituting (4.22) - (4.25) into (4.21) and combining with (4.12) and (4.13), we obtain

( t ( x,t ),w( t ) ) ( μ( 0 )( β+1 )C β+2 u( t ) β+2 β+2 +Cρ( γ+1 ) u( t ) γ u t ( t ) + μ( 0 )C β+2 + m 1 C β+2 + m 1 ( β+1 )C β+2 0 μ ( ts ) u( s ) β+2 β+2 ds+C( M ) u t ( t ) ) w( t ) V 1 σ C( ) w( t ) V 1 σ , (4.26)

where

C( )= sup ( u 0 , u 1 ),t + { μ( 0 )( β+1 )C β+2 u( t ) β+2 β+2 +Cρ( γ+1 ) u( t ) γ u t ( t ) + μ( 0 )C β+2 + m 1 C β+2 + m 1 ( β+1 )C β+2 0 μ ( ts ) u( s ) β+2 β+2 ds+C( M ) u t ( t ) }.

In summary, combining (4.19) and (4.16), we conclude that C( ) is bounded. This implies that t ( x,t ) V 1 σ belongs to the dual space V σ 1 of V 1 σ , and its norm is bounded in V σ 1 . Consequently, (4.7) and (4.8) hold.

Case 3: For n4 , it follows from (4.9) and (4.10) that

( x,t )=f( u ), t ( x,t )= f ( u ) u t .

Under the assumptions of Theorem 3.1, for all uH , we have f ( u ) L ( Ω ) . Thus, the conclusion holds.

Therefore, for all n1 , Lemma 4.3 holds.

Lemma 4.3. For sufficiently large t , the operator S 1 ( t ) is uniformly compact in V×H .

Proof. Let u ¯ be a solution to equation (4.1). By differentiating equation (4.1) with respect to t , we conclude that w= u ¯ t satisfies the following equation:

{ w tt Δwσ( u 2 )Δ w t +2 σ ( u 2 )( Δu, u t )Δw = 0 μ ( ts ) | u( s ) | β u( s )ds+μ( 0 ) | u( t ) | β u( t ) ρ( γ+1 ) | u( t ) | γ u t ( t ) f ( u ) u t ,xΩ,t>0, w( 0 )=0, w t ( 0 )=hf( u 0 ),xΩ,t0. (4.27)

From Lemma 4.2, we obtain

t ( x,t ) C b ( + ; V σ 1 ),

and thus

w= u ¯ t C b ( + ; V σ ), w t = u ¯ tt C b ( + ; V σ 1 ), w tt = u ¯ ttt C b ( + ; V σ ).

From (4.1), it follows that

0 μ ( ts ) | u( s ) | β u( s )dsρ | u( t ) | γ u( t )f( u )+h( x ) C b ( + , L 2 ( Ω ) ),

and

A u ¯ C b ( + , V σ'1 ).

Thus, we have u ¯ C b ( + , V σ +1 ) . Furthermore, by Lemma 4.1, u ¯ t is bounded, and u ¯ t C( + , V σ ) , which implies Δw=Δ u ¯ t C b ( + ; V σ ) .

Consequently,

( u ¯ , u ¯ t )C( + ; V σ +1 × V σ ),

and

t t 0 S( t ) V σ +1 × V σ .

By the embedding theorem, for any σ 1 > σ 2 , we have V σ 1 V σ 2 . Since V σ +1 × V σ V×H , it follows that t t 0 S( t ) is compact in V×H . Therefore, S 1 ( t ) is uniformly compact in V×H , and the conclusion is proved.

Theorem 4.2. For any bounded set , there exists t 0 >0 such that t t 0 S( t ) is relatively compact in V×H . Therefore, S( t ) is asymptotically compact in V×H .

Proof. We decompose S( t ) into S 1 ( t ) and S 2 ( t ) , i.e., S( t )= S 1 ( t )+ S 2 ( t ) , where S 1 ( t )=( t ) and S 2 ( t )=S( t )( t ) . For any solution u ˜ to equation (4.2) with initial condition u 0 , the following inequality holds:

u ˜ t 2 + u ˜ 2 C( )exp( α 1 t ).

Thus, as t , we have

α( S 2 ( t )( ) )0.

By Lemma 4.3, the operator S 1 ( t ) is uniformly compact in V×H . Therefore, as t ,

α( S( t ) )α( S 1 ( t ) )+α( S 2 ( t ) )=α( S 2 ( t ) )0.

This completes the proof.

Theorem 4.3. If the external force term h L 2 ( Ω ) , the nonlinear term f C 1 ( ) satisfies ( f 1 ) - ( f 2 ) , and ( h 1 ) - ( h 4 ) hold, then the solution semigroup S( t ) associated with Equation (1.2) possesses a compact and connected global attractor in the space V×H .

Proof. Following the same reasoning, Theorem 4.3 can be derived by applying Theorem 4.1 and Theorem 4.2 within the framework of Theorem 2.1.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

References

[1] Hale, J.K. (1971) Introduction. In: Applied Mathematical Sciences, Springer, 1-10.
https://doi.org/10.1007/978-1-4615-9968-5_1
[2] Duong Toan, N. (2022) Time-dependent Global Attractors for Strongly Damped Wave Equations with Time-Dependent Memory Kernels. Dynamical Systems, 37, 466-492.
https://doi.org/10.1080/14689367.2022.2072710
[3] D’Abbicco, M. (2014) The Influence of a Nonlinear Memory on the Damped Wave Equation. Nonlinear Analysis: Theory, Methods & Applications, 95, 130-145.
https://doi.org/10.1016/j.na.2013.09.006
[4] Pata, V., Di Plinio, F. and Zelik, S. (2008) On the Strongly Damped Wave Equation with Memory. Indiana University Mathematics Journal, 57, 757-780.
https://doi.org/10.1512/iumj.2008.57.3266
[5] Bae, J.J., Park, J.Y. and Jeong, J.M. (2003) On Uniform Decay of Solutions for Wave Equation of Kirchhoff Type with Nonlinear Boundary Damping and Memory Source Term. Applied Mathematics and Computation, 138, 463-478.
https://doi.org/10.1016/s0096-3003(02)00163-7
[6] Kang, J. (2016) Asymptotic Behavior of the Thermoelastic Suspension Bridge Equation with Linear Memory. Boundary Value Problems, 2016, Article No. 206.
https://doi.org/10.1186/s13661-016-0707-3
[7] Ma, Q., Wang, J. and Liu, T. (2019) Time-dependent Attractor of Wave Equations with Nonlinear Damping and Linear Memory. Open Mathematics, 17, 89-103.
https://doi.org/10.1515/math-2019-0008
[8] Ma, Q., Wang, J. and Liu, T. (2018) Time-Dependent Asymptotic Behavior of the Solution for Wave Equations with Linear Memory. Computers & Mathematics with Applications, 76, 1372-1387.
https://doi.org/10.1016/j.camwa.2018.06.031
[9] Mei, X. and Zhu, K. (2023) Asymptotic Behavior of Solutions for Hyperbolic Equations with Time-Dependent Memory Kernels. Discrete and Continuous Dynamical Systems, 28, 1855-1885.
https://doi.org/10.3934/dcdsb.2022150
[10] Feng, B., Pelicer, M.L. and Andrade, D. (2016) Long-Time Behavior of a Semi-Linear Wave Equation with Memory. Boundary Value Problems, 2016, Article No. 37.
https://doi.org/10.1186/s13661-016-0551-5
[11] Yang, M. and Sun, C. (2009) Attractors for Strongly Damped Wave Equations. Nonlinear Analysis: Real World Applications, 10, 1097-1100.
https://doi.org/10.1016/j.nonrwa.2007.12.001
[12] Han, Y.H., Zhen, G.Y. and Zheng, G.J. (2012) Global Attractor for Damped Wave Equations with Nonlinear Memory. Journal of Mathematical Research with Applications, 32, 213-222.
[13] D’Abbicco, M. (2014) A Wave Equation with Structural Damping and Nonlinear Memory. Nonlinear Differential Equations and Applications, 21, 751-773.
https://doi.org/10.1007/s00030-014-0265-2
[14] Zhou, C. and Sun, C. (2022) Global Attractor for a Wave Model with Nonlocal Nonlinear Damping. Journal of Mathematical Analysis and Applications, 507, Article 125818.
https://doi.org/10.1016/j.jmaa.2021.125818
[15] Kirchhoff, G. (1883) Vorlesungenüber Mechanik. Teubner.
[16] Chueshov, I. (2012) Long-time Dynamics of Kirchhoff Wave Models with Strong Nonlinear Damping. Journal of Differential Equations, 252, 1229-1262.
https://doi.org/10.1016/j.jde.2011.08.022
[17] Narciso, V. (2020) On a Kirchhoff Wave Model with Nonlocal Nonlinear Damping. Evolution Equations & Control Theory, 9, 487-508.
https://doi.org/10.3934/eect.2020021
[18] Temam, R. (1988) Infinite-Dimensional Dynamical Systems in Mechanics and Physics. Springer-Verlag.
[19] Babin, A.V. and Vishik, M.I. (1992) Attractors of Evolution Equations. North Holland.

Copyright © 2025 by authors and Scientific Research Publishing Inc.

Creative Commons License

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.