The Time-Dependent Attractors for the Wave Equation with Fading Memory and Structural Damping

Abstract

In this article, the asymptotic behavior of the solutions to the wave equation with fading memory and structural damping t 2 uΔu+γ ( Δ ) θ t u 0 μ( s )Δ η t ( s )ds +f( u )=g( x ) is considered. First of all, when the growth exponent of nonlinear terms f satisfies 2p3+2θ , the well-posedness of solutions is obtained by applying Faedo-Galerkin approximation method and energy estimation; secondly, the asymptotic compactness of the solution process is proved via the method of contraction function; finally, the existence of time-dependent global attractor is obtained in the natural energy space H 0 1 ( Ω )× L 2 ( Ω )× L μ 2 ( + ; H 0 1 ( Ω ) ) .

Share and Cite:

Wang, X. and Yan, L. (2026) The Time-Dependent Attractors for the Wave Equation with Fading Memory and Structural Damping. Journal of Applied Mathematics and Physics, 14, 669-693. doi: 10.4236/jamp.2026.142036.

1. Introduction

In this paper, we investigate the long-time behavior of solutions to the wave equation with fading memory and structural damping:

t 2 uΔu+γ ( Δ ) θ t u 0 μ( s )Δ η t ( s )ds +f( u ) =g( x ),( x,t )Ω×[ τ,+ ), (1.1)

u| Ω =0,u( x,τ )= u τ ( x ), t u( x,τ )= t u τ ( x ),xΩ, (1.2)

where θ( 1 2 ,1 ) , γ>0 is the damping coefficient, and Ω 3 is a bounded domain with smooth boundary.

Suppose that the nonlinear function f satisfies the following conditions:

(M1) f C 2 ( ) , f( 0 )=0 . For any s , f satisfies

liminf | s |+ f( s ) s > λ 1 , (1.3)

as well as the growth condition:

| f ( s ) | C 0 ( 1+ | s | p2 ),2p p θ =3+2θ, (1.4)

where the constant C 0 >0 , and λ 1 is the first eigenvalue of the operator Δ under Dirichlet boundary conditions.

Following the ideas in [1] [2], the memory kernel function μ and the forcing term g satisfy the following conditions:

(M2) μ C 1 ( + ) L 1 ( + ) and

0 μ( s )ds = k 0 <. (1.5)

And there exists a positive constant k such that

μ ( s )kμ( s )0,s0 .(1.6)

(M3) g L 2 ( Ω ) .

Remark 1. From (1.3), there exists a constant β 0 satisfying 0< β 0 <1 such that

F( s ),1 ( 1 β 0 ) λ 1 2 s 2 C β 0 ,

f( s ),s ( 1 β 0 ) λ 1 s 2 C β 0 ,s,

holds, where F( s )= 0 s f( r )dr .

In recent years, wave equations with damping and memory terms have extensive applications in the viscoelastic dynamical systems while viscoelastic material serves as the medium for energy transmission, as can be found in [1] [2].

When Equation (1.1) does not contain a memory term (that is, when μ is a Dirac measure at a fixed time or μ takes the zero value), and the damping dissipation exponent satisfies θ=0 or θ=1 , Equation (1.1) reduces to a weakly damped or strongly damped wave equation. For this model, when the nonlinear term satisfies subcritical or critical growth conditions, Pata et al. discussed the well-posedness of solutions and long-time dynamical behaviors in [3] [4]. In addition to strong damping and weak damping, there exists another type of damping called structural damping, which is usually expressed in the form of a fractional power of the operator Δ , namely ( Δ ) θ u t with 0<θ<1 . Its dissipation strength lies between weak damping and strong damping, and it can more essentially reflect the internal friction effect of materials. For wave equations with structural damping, when 0<θ< 1 2 , Savostianov investigated the existence of attractors and exponential attractors for wave equations with structural damping under the subcritical growth of the nonlinear term in [5]. When 1 2 <θ<1 , Wang Xuan et al. studied a Kirchhoff wave equation model with structural damping in [6]:

ϵ( t ) t 2 uM( u 2 )Δu+ ( Δ ) γ t u+f( u )=g( x ),γ( 1 2 ,1 ).

They gained the well-posedness and regularity of solutions to Kirchhoff-type wave equations with time-dependent coefficients and structural damping, and proved the existence of time-dependent global attractors by using the method of contraction functions. In addition, many scholars have carried out extensive research on equations with structural damping, leading to a large number of research results on this model, as shown in references (see [5]-[9] and related literature).

For dissipative evolution equations with fading memory, Dafermos systematically established the theoretical framework of memory kernels for viscoelastic models in [2], established the well-posedness of linear viscoelastic wave equations with exponentially fading memory kernels within this framework. In [10], Lions combined the Faedo-Galerkin method with compactness theorems to solve the existence problem of weak solutions for nonlinear wave equations with fading memory terms. In [11], Chueshov systematically expounded the global attractor theory of wave equations with fading memory terms and established a unified framework for analyzing the long-term dynamical behaviors of such systems. On the basis of this theory, Ma Qiaozhen et al. discussed the asymptotic behaviors of solutions to wave equations with linear memory in time-dependent spaces in [12].

Inspired by the aforementioned research findings, this paper investigates the wave equation with fading memory terms and structural damping. To the best of our knowledge, the long-time dynamical behavior of solutions to wave equations with damping dissipation exponent 1 2 <θ<1 and fading memory terms in time-dependent spaces has not yet been explored. Meanwhile, the structural damping term, nonlinear term, and fading memory term in the equation bring essential difficulties to the derivation of dissipative estimates for solutions, the verification of the existence of bounded absorbing sets, and the proof of asymptotic compactness of the solution process. When the damping dissipation exponent satisfies 1 2 <θ<1 , classical Sobolev embeddings fail to guarantee compactness; in addition, the fading memory term introduces a history-dependent component, making it difficult for traditional energy functionals to simultaneously characterize the instantaneous dissipation of fractional damping and the cumulative effect of memory terms. To address these challenges, we construct a memory-coupled energy functional that can simultaneously describe the effects of structural damping and fading memory. By combining energy estimation techniques, asymptotic regularity estimates, as well as the contraction function tailored to the memory-damping coupling and the relevant theory of time-dependent attractors, we overcome these technical obstacles. Furthermore, we establish the well-posedness of solutions to problem (1.1) - (1.2) and verify their Lipschitz continuity in the space H 0 1 ( Ω )× L 2 ( Ω )× L μ 2 ( + ; H 0 1 ( Ω ) ) . Subsequently, we confirm the asymptotic compactness of the solution process, and finally prove the existence of time-dependent global attractors for problem (1.1) - (1.2) in the space H 0 1 ( Ω )× L 2 ( Ω )× L μ 2 ( + ; H 0 1 ( Ω ) ) .

The content and structure of this paper are arranged as follows: In Section 2, we review the preliminary results; in Section 3, we discuss the well-posedness of weak solutions; in Section 4, we prove the asymptotic compactness of the process by using the method of contraction functions, and then obtain the existence of the time-dependent attractor.

In this paper, the symbol C denotes a positive constant. Each occurrence of C in different formulas represents the corresponding positive constant. We also use C i ,i to denote different positive constants, and C( , ) denotes a constant depending on the parameters in the parentheses.

2. Notations and Preliminary Results

Following the ideas in [2] [13] [14], we introduce the history function η= η t ( x,s ) of u , which satisfies

t η t = s η t + t u (2.1)

with the corresponding initial value conditions

{ u( x,t )=0, xΩ,t>τ, η t ( x,s )=0, ( x,s )Ω× + ,t>τ, u( x,t,τ )= u τ ( x ), xΩ,tτ, t u( x,t,τ )= u τ ( x ), xΩ,tτ, η τ ( x,s )= η τ ( x,s ), ( x,s )Ω× + .

We set A=Δ with domain D( A )= H 0 1 ( Ω ) H 2 ( Ω ) .

Consider the family of Hilbert spaces D( A s 2 ) for s , equipped with their respective inner products and norms:

, D( A s 2 ) = A s 2 , A s 2 , u D( A s 2 ) 2 = A s 2 u 2 ,

where , and denote the inner product and norm in L 2 ( Ω ) , respectively.

For the sake of convenience, we introduce the notation V s =D( A s 2 ) for s , with the inner product and norm expressed as follows:

u,v s = Ω A s 2 u( x ) A s 2 v( x )dx , u s 2 = Ω | A s 2 u( x ) | 2 dx ,u,v V s .

Then, V 0 = L 2 ( Ω ) , D( A s 2 )= V s and D( A s 2 )= V s .

By virtue of the Sobolev embedding theorem, we obtain the compact embedding

V s 1 ↪↪ V s 2 ,for s 1 > s 2 , (2.2)

as well as the continuous embedding

V s L 2n n2s . (2.3)

Therefore, the problem (1.1) - (1.2) can be rewritten in the following form:

t 2 u+Au+γ ( Δ ) θ t u+ 0 μ( s )A η t ( s )ds +f( u ) =g( x ),( x,t )Ω×[ τ,+ ), (2.4)

u( x,t )=0,xΩ,t>τ, (2.5)

η τ ( x,s )=0,( x,s )Ω× + ,t>τ, (2.6)

u( x,t,τ )= u τ ( x ), t u( x,t,τ )= t u τ ( x ), η τ ( x,s )= η τ ( x,s ),xΩ,tτ. (2.7)

Define the family of Hilbert spaces:

t 1+ϑ = V 1+ϑ × V ϑ × L μ 2 ( + ; V 1+ϑ ),

equipped with the corresponding inner product and norm

z( t ) t 1+ϑ 2 = ( u( t ), t u( t ), η t ( s ) ) t 1+ϑ 2 = u( t ) 1+ϑ 2 + t u( t ) ϑ 2 + η t ( s ) μ,1+ϑ 2 , (2.8)

In particular, for ϑ=0 , the family of Hilbert spaces t 1 is defined as follows:

t 1 = V 1 × L 2 ( Ω )× L μ 2 ( + ; V 1 ), (2.9)

with its norm given by

z( t ) t 1 2 = ( u( t ), t u( t ), η t ( s ) ) t 1 2 = u( t ) 1 2 + t u( t ) 2 + η t ( s ) μ,1 2 . (2.10)

Furthermore, for ϑ>0 , we have the compact embedding

t 1+ϑ ↪↪ t 1 .

Next, we will review the following concepts and some abstract results, which will be used to study the long-time dynamical behavior of solutions in time-dependent spaces.

Definition 2.1. ([15]) Let { X t } t be a family of normed spaces. A two-parameter family of operators { U( t,τ ): X τ X t ,tτ,τ } is said to be a process, if for any τ ,

1) U( τ,τ )=Id is the identity operator on X τ ;

2) U( t,s )U( s,τ )=U( t,τ ) , τst .

Let { X t } t be a family of normed spaces. For every t , the R -ball of X t is defined by:

B t ( R )={ z X t | z X t R }.

Definition 2.2. ([15]) A family = { C t } t of bounded sets C t X t is called uniformly bounded, if there exists a constant R>0 such that C t B t ( R ) , for all t .

Definition 2.3. ([15]) A uniformly bounded family B t = { B t ( R 0 ) } t is called a time-dependent absorbing set for the process U( t,τ ) , if for any R>0 , there exist a t 0 = t 0 ( R )t and R 0 >0 such that

τt t 0 U( t,τ ) B τ ( R ) B t ( R 0 ).

A process is said to be dissipative if it possesses a time-dependent absorbing set.

Lemma 2.4. ([6]) Let { x n } be a bounded sequence and also let ψC( ) be a monotonic function. Then

ψ( liminf n x n ) liminf n ψ( x n ).

Lemma 2.5. ([16] [17]) Let X , B and Y be three Banach spaces. For any T>0 , if X ↪↪ B Y , and

W 1 ={ u L p ( [ 0,T ];X )| t u L r ( [ 0,T ];Y ),r>1,1p< },

W 2 ={ u L ( [ 0,T ];X )| t u L r ( [ 0,T ];Y ),r>1 }.

Then the following compact embeddings hold:

W 1 ↪↪ L p ( [ 0,T ];B ), W 2 ↪↪ C( [ 0,T ];B ).

Theorem 2.6. ([15]) If U( t,τ ) is asymptotically compact, that is, the set

K={ K= { K t } t | K t X t iscompactandKisattracting }

is nonempty, then the time-dependent attractor A exists and is unique.

Definition 2.7. ([15] [18] [19]) A time-dependent attractor A= { A t } t is invariant, if for all τt ,

U( t,τ ) A τ = A t .

Theorem 2.8. ([20]) Let U( , ) be a process acting on a family of Banach spaces { X t } t . Then U( , ) has a time-dependent global attractor A * = { A t * } t satisfying

A t * = st τs U( t,τ ) B τ ( R ) ¯

if and only if

1) U( , ) has a time-dependent absorbing family B= { B t ( R 0 ) } t ;

2) U( , ) is asymptotically compact.

Definition 2.9. ([21] [22]) Let { X t } t be a family of Banach spaces and let = { C t } t be a uniformly bounded family of subsets of { X t } t . A function Φ τ t ( , ) defined on X t × X t is called a contraction function on C τ × C τ , if for any fixed t and any sequence { x n } n=1 C τ , there exists a subsequence { x n k } k=1 { x n } n=1 such that

lim k lim l Φ τ t ( x n k , x n l )=0,

where τt .

We denote by ( C t ) the set of all contraction functions on C t × C t .

Theorem 2.10. ([20]) Let U( , ) be a process on { X t } t that possesses a time-dependent absorbing family of sets B t = { B t ( R 1 ) } t . If for any ε>0 , there exists a subsequence T( ε )t , Φ T t ( B T ( R ) ) such that

U( t,T )xU( t,T )y ε+ Φ T t ( x,y ),x,y B T ( R )

for any fixed t , then U( , ) is asymptotically compact.

3. Well-Posedness of Solutions

First, we give the following definition of the solution to the problem (2.4) - (2.7).

Definition 3.1. For any τ , a triple z=( u, t u, η t ) is called a weak solution to the problem (2.4) - (2.7) on the interval [ τ,T ] , if

u L ( [ τ,T ]; V 1 ), t u L ( [ τ,T ]; L 2 ( Ω ) ) L 2 ( [ τ,T ]; V θ ), η t L ( [ τ,t ]; L μ 2 ( + ; V 1 ) )

and it satisfies

t 2 u,ω + Au,ω +γ A θ t u,ω + 0 μ( s )A η t ( s )ds ,ω + f( u ),ω = g( x ),ω

for all τt and all ω V 1 .

Theorem 3.2. Suppose that (M1) - (M3) hold. Then for each T>τ and θ( 1 2 ,1 ) , the problem (2.4) - (2.7) admits a weak solution y=( u, t u, η t )C( [ τ,T ]; t 1 ) with t 2 u L ( [ τ,T ]; V 2θ ) , which satisfies

u( t ) 1 2 + t u( t ) 2 + η t μ,1 2 + t 2 u( t ) 2θ 2 + τ T t u( t ) θ dt C( R, β 0 , g , λ 1 , C 2 ),tτ. (3.1)

Furthermore, the weak solution satisfies the following properties:

1) (Dissipativity) There exists a positive constant R independent of θ( 1 2 ,1 ) such that

( u, t u, η t ) t 1 R 0 ,t t 0 ( R ), (3.2)

where τt t 0 ( R ) and t 0 ( R ) is a moment depending on the positive constant R .

2) (Energy equality) For each τtT , the following energy equality holds

E( u( t ), t u( t ), η t ( s ) )+2γ τ T t u( t ) θ 2 dt + 1 2 τ T 0 μ( s ) d ds η t ( s ) 1 2 dsdt =E( u τ , t u τ , η τ ( s ) ). (3.3)

Here,

E( u( t ), t u( t ), η t ( s ) ) = u( t ) 1 2 + t u( t ) 2 + η t ( s ) μ,1 2 +2 F( u( t ) ),1 2 g,u( t ) . (3.4)

3) (Weak lipschitz stability) The solution y=( u, t u, η t ) is Lipschitz continuous in the space V 1θ × V θ × L μ 2 ( + ; V 1θ ) , that is,

u( t ) 1θ 2 + t u( t ) θ 2 + η t ( s ) μ,1θ 2 m 3 m 2 e C ˜ ( tτ ) ( u τ 1θ 2 + t u τ θ 2 + η τ ( s ) μ,1θ 2 )+ C ˜ 0 ( tτ ) m 2 e C ˜ ( tτ ) . (3.5)

where z ˜ =( u ˜ , t u ˜ , η ˜ t )= z 1 z 2 , and z i =( u i , t u i , η i t )( i=1,2 ) are two weak solutions to the problem (2.4) - (2.7) corresponding to the initial values ( u τ i , t u τ i , η τ i )( i=1,2 ) , respectively. Moreover

C ˜ =C( R, β 0 , g , λ 1 , C 1 , C 2 ,γ,k ),

C ˜ 0 =C( R, β 0 , g , λ 1 , C 2 ,γ ).

Proof.

1) (Existence of Weak Solutions) First, we establish some a priori estimates for the solutions to the problem (2.4) - (2.7). Taking the inner product of Equation (2.4) with t u in L 2 ( Ω ) , we obtain

d dt E( u( t ), t u( t ), η t ( s ) )+2γ t u θ 2 + 0 μ( s ) d ds η t ( s ) 1 2 ds =0,

where

E( u( t ), t u( t ), η t ( s ) ) = u( t ) 1 2 + t u( t ) 2 + η t ( s ) μ,1 2 +2 F( u( t ) ),1 2 g,u( t ) . (3.6)

Integrating the above identity over the interval [ s,t ] , we deduce that (3.3) holds.

Since

0 μ( s ) d ds η t ( s ) 1 2 ds = 0 μ ( s ) η t ( s ) 1 2 ds k 0 μ( s ) η t ( s ) 1 2 ds =k η t ( s ) μ,1 2 ,

we have

E( u( t ), t u( t ), η t ( s ) )+2γ τ T t u( t ) θ 2 dt E( u τ , t u τ , η τ ). (3.7)

By virtue of (1.4) and the compact embedding V 1 ↪↪ L p+1 ( Ω ) , it follows that

2 F( u ),1 2 C 1 ( u 2 + u L p+1 ( Ω ) p+1 ) C 2 ( u 1 2 + u 1 p+1 ). (3.8)

From (M3), we know that

2| g,u | β 0 4 u 1 2 + 4 β 0 λ 1 g 2 . (3.9)

Therefore,

E( u τ , t u τ , η τ ) = u τ 1 2 + t u τ 2 + η τ μ,1 2 +2 F( u τ ),1 2 g, u τ u τ 1 2 + t u τ 2 + η τ μ,1 2 + C 2 ( u τ 1 2 + u τ 1 p+1 )+ β 0 4 u τ 1 2 + 4 β 0 λ 1 g 2 m 0 ( u τ 1 2 + t u τ 2 + η τ μ,1 2 )+ C 2 ( u τ 1 p+1 )+C C( R, β 0 , g , λ 1 , C 2 ),

where m 0 =max{ 1,1+ C 2 + β 0 4 } and C= 4 β 0 λ 1 g 2 .

From Remark 1.1, we have

2 Ω F( u )dx ( 1 β 0 ) u( t ) 1 2 2 C β 0 ( β 0 1 ) u( t ) 1 2 2 C β 0 . (3.10)

By virtue of estimates (3.8) - (3.9), we obtain

m 1 ( u( t ) 1 2 + t u( t ) 2 + η t μ,1 2 )C E( u( t ), t u( t ), η t ( s ) )C( R, β 0 , g , λ 1 , C 2 ). (3.11)

where m 1 =min{ 1, 3 β 0 4 } and C= 4 β 0 λ 1 g 2 +2C( β 0 ) .

According to (3.7) and (3.11), we conclude that

τ T t u( t ) θ 2 dt C( R, β 0 , g , λ 1 , C 2 ). (3.12)

Using the embeddings L 1+ 1 p ( Ω ) ↪↪ V 1 ↪↪ V 2θ and Equation (2.3), we get

t 2 u( t ) 2θ 2 u( t ) 22θ 2 +γ t u( t ) 2 + η t ( s ) μ,22θ 2 + f( u ) 2θ 2 + g 2θ 2 C( R, β 0 , g , λ 1 , C 2 )( u( t ) 22θ 2 +γ t u( t ) 2 + η t ( s ) μ,22θ 2 + f( u ) L 1+ 1 p 2 + g 2 ) C( R, β 0 , g , λ 1 , C 2 )( u( t ) 1 2 + u( t ) 1 2p +γ t u( t ) 2 + η t ( s ) μ,1 2 + g 2 ) C( R, β 0 , g , λ 1 , C 2 ). (3.13)

Combining (3.7), (3.11) - (3.13), we deduce that (3.1) holds.

Let us prove the existence of solutions for the problem (2.4) - (2.7) in the space C( [ τ,T ]; t 1 ) . Suppose that { w j } j=1 is an orthonormal basis of V 1 with A w j = λ j w j for j=1,2, . Let { ς j } j=1 be an orthonormal basis of L μ t 2 ( + ; V 1 ) satisfying A ς j = λ j ς j for j=1,2, . For each n , there exist finite-dimensional subspaces

V n =span{ w 1 ,, w n } V 1 , M n =span{ ς 1 ,, ς n } L μ 2 ( + ; V 1 ).

Define P n : V 1 V n as the orthogonal projection onto V n and Q n : L μ 2 ( + ; V 1 ) M n as the orthogonal projection onto M n .

For each n , let z n ( t )=( u n , t u n , η n t ) be an approximate solution to the problem, where u n = j=1 n T j n ( t ) w j with T j n C 1 ( [ τ,T ] ) , and η n t = j=1 n Λ j n ( t ) ς j with Λ j n C 1 ( [ τ,T ] ) . Then for any φ V n and each t[ τ,T ] , z n ( t )=( u n , t u n , η n t ) satisfies

t 2 u,φ + Au,φ +γ A θ t u,φ + 0 μ( s )A η t ( s )ds ,φ + f( u ),φ = g( x ),φ , (3.14)

together with

η n t ( s )={ u n ( t ) u n ( ts ), 0stτ, η τ n ( st+τ )+ u n ( t ) u τ n , s>tτ.

Multiplying Equation (3.14) by ψ C 0 ( [ τ,T ] ) and integrating over [ τ,T ] , we obtain

τ T ψ ( t 2 u,φ + Au,φ +γ A θ t u,φ + 0 μ( s )A η t ( s )ds ,φ + f( u ),φ g( x ),φ )dt=0. (3.15)

Then we have the following results:

u n isboundedin L ( [ τ,T ]; V 1 ),

t u n isboundedin L ( [ τ,T ]; L 2 ( Ω ) ) L 2 ( [ τ,T ]; V θ ),

η n t isboundedin L ( [ τ,T ]; L μ 2 ( + , V 1 ) ),

t 2 u n isboundedin L ( [ τ,T ]; V 2θ ).

By applying the Galerkin approximation method, there exists z=( u, t u, η t ) L ( [ τ,T ]; t 1 ) such that

u n u weakly * in L ( [ τ,T ]; V 1 ),

t u n t u weakly * in L ( [ τ,T ]; L 2 ( Ω ) ),

t u n t uweaklyin L 2 ( [ τ,T ]; V θ ),

η n t η t weakly * in L ( [ τ,T ]; L μ 2 ( + , V 1 ) ),

t 2 u n t 2 u weakly * in L ( [ τ,T ]; V 2θ ).

Applying Lemma 2.2 and Alaoglu Theorem, for 0<α1 , we deduce that

u n uinC( [ τ,T ]; V 1α ),

t u n t uinC( [ τ,T ]; V α ),

u n uin L 2 ( [ τ,T ]; V 1 )and u n ( x,t )u( x,t ),a.e.( x,t )Ω×[ τ,T ],

t u n t uin L 2 ( [ τ,T ]; L 2 ( Ω ) ),

f( u n )f( u )weaklyin L 1+ 1 p ( [ τ,T ]; L 1+ 1 p ( Ω ) ).

Since for arbitrary φ V n , we have

τ T f( u n )f( u ),φ dt C 1 τ T ( 1+ | u n | p1 + | u | p1 )| u n u || φ |dt C 1 ( τ T ( 1+ | u n | p1 + | u | p1 ) p+1 p1 dt ) p1 p+1 ( Ω | u n u | p+1 dt ) 1 p+1 ( Ω | φ | p dt ) 1 p+1 C 1 ( 1+ u n L p+1 p1 + u L p+1 p1 ) u n u L p+1 φ L p+1 C 2 ( 1+ u n 1 p1 + u 1 p1 ) u n u 1 φ 1 C( R, β 0 , g , λ 1 , C 2 ) u n u L 2 ( [ τ,T ]; V 1 ) 0,

Since f( u n ) and f( u ) are bounded in L 1+ 1 p ( [ τ,T ]; L 1+ 1 p ( Ω ) ) , by the Lebesgue Dominated Convergence Theorem, we obtain

τ T ψ f( u n )f( u ),φ dy 0.

Furthermore, for arbitrary φ V n , we have

τ T A u n Au,φ dt τ T A 1 2 ( u n ( t )u( t ) ) A 1 2 φ dt τ T u n ( t )u( t ) 1 φ 1 dt 0.

In addition, since u n and u are bounded in V 1 , an application of the Lebesgue Dominated Convergence Theorem yields

τ T ψ A u n Au,φ dt 0.

We further set

η ¯ n t = η n t η t , u ¯ n = u n u, η ¯ τ n = η τ n η τ , u ¯ τ n = u τ n u τ .

Consider the mapping p n :[ τ,T ] defined by t p n ( t )= u ¯ n ( t ),φ 1 , where p n 0 in L 2 ( [ τ,T ]; V 1 ) .

From Equation (2.2), we have

η ¯ n t μ,1 2 = 0 tτ μ( s ) u ¯ n ( ts ) u ¯ ( ts ) 1 2 ds + tτ μ( s ) η ¯ τ n ( st+τ )+ u ¯ n ( t ) u ¯ τ n 1 2 ds 0 tτ μ( s ) u ¯ n ( ts ) u ¯ ( ts ) 1 2 ds +2 tτ μ( s ) η ¯ τ n ( st+τ ) 1 2 ds +2 tτ μ( s ) u ¯ n ( t ) u ¯ τ n 1 2 ds ,

where

tτ μ( s ) η ¯ τ n ( st+τ ) 1 2 ds = 0 μ( s+tτ ) η ¯ τ n ( s ) 1 2 ds 0 0 μ( s ) η ¯ τ n ( s ) 1 2 ds 0.

As n , u ¯ n 0 in L 2 ( [ τ,T ]; V 1 ) , hence

u ¯ n ( t ) u ¯ ( ts ) u ¯ n ( t ) + u ¯ n ( ts ) 0.

Combining with (1.5) in (M1), we obtain

0 tτ μ( s ) u ¯ n ( t ) u ¯ ( ts ) 1 2 ds 0.

Similarly,

tτ μ( s ) u ¯ n ( t ) u ¯ τ n 1 2 ds 0.

Therefore,

η ¯ t ( s ) μ,1 2 0,t[ τ,T ].

Since

η n t ( s )={ u n ( t ) u n ( ts ), 0stτ, η τ n ( st+τ )+ u n ( t ) u τ n , s>tτ,

we have

0 μ ( s ) η ¯ n t ( s ),φ 1 ds = 0 tτ μ ( s ) u ¯ n ( t ),φ 1 ds 0 tτ μ ( s ) u ¯ n ( ts ),φ 1 ds + tτ μ ( s ) u ¯ n ( t ),φ 1 ds+ tτ μ ( s ) η ¯ τ n ( st+τ ) u ¯ τ n ,φ 1 ds.

By applying (1.5) in (M1) again, we obtain

0 μ( s ) u ¯ n ( t ),φ 1 ds = 0 μ( s ) p n ( t )ds = k 0 p n ( t )0,

0 tτ μ( s ) u ¯ n ( ts ),φ 1 ds = τ t μ( ts ) p n ( s )ds 0,

tτ μ( s ) η ¯ τ n ( st+τ ) u ¯ τ n ,φ 1 ds φ 1 0 μ( s+tτ )( η ¯ τ n ( s ) 1 + u ¯ τ n 1 )ds φ 1 η ¯ τ n ( s ) L μ 2 ( + ; V 1 ) + φ 1 u ¯ τ n L 2 ( [ τ,T ]; V 1 ) 0,

Therefore, we have

lim n 0 μ( s ) η ¯ n t ( s ),φ 1 ds =0,t[ τ,T ].

By applying the Lebesgue Dominated Convergence Theorem, we deduce that

lim n τ T ψ( t ) 0 μ( s ) η ¯ n t ( s ),φ 1 dsdt =0.

As a result, letting n in Equation (3.15), we conclude that z=( u, t u, η t ) is a weak solution to the problem and satisfies the estimate (3.1).

Next, we shall prove that the solution z=( u, t u, η t ) to the problems (2.4) - (2.7) belongs to C( [ τ,T ]; t 1 ) .

Since ( u, t u )C( [ τ,T ]; V 1α × V α ) and ( u, t u, η t ) L ( τ,T; t 1 ) , it follows that ( u, t u, η t ) C w ( [ τ,T ]; t 1 ) and

( u, t u, η t ) t 1 liminf τt ( u τ , t u τ , η τ ) τ 1 .

For any t[ τ,T ] , by (3.3) we have

lim τt E( u τ , t u τ , η τ )=E( u( t ), t u( t ), η t ( s ) ). (3.16)

From (3.16), we deduce that u( x,τ )u( x,t ) a.e. xΩ as τt . Applying Lemma 2.4, Remark 1.1 and Fatou Lemma, we have

lim τt 2 g, u τ =2 g,u( t ) ,

( u( t ), t u( t ), η t ) t 1 2 liminf τt ( u τ , t u τ , η τ ) τ 1 2 ,

Ω ( 2F( u( t ) )+( 12 C β 0 ) λ 1 | u( t ) | 2 +2 C β 0 )dx liminf τt Ω ( 2F( u τ )+( 12 C β 0 ) λ 1 | u τ | 2 +2 C β 0 )dx liminf τt Ω 2F( u τ )dx +( 12 C β 0 ) λ 1 u 2 +2 C β 0 | Ω |,

That is,

Ω 2F( u( t ) )dx liminf τt Ω 2F( u τ )dx .

From the above estimates and (3.16), we get

liminf τt t u τ 2 + liminf τt u τ 1 2 + liminf τt η τ μ,1 2 + liminf τt 2 F( u τ ),1 lim τt [ t u τ 2 + u τ 2 + η τ μ,1 2 +2 F( u τ ),1 ] = t u( t ) 2 + u( t ) 1 2 + η t μ,1 2 +2 F( u( t ) ),1 liminf τt t u τ 2 + liminf τt u τ 1 2 + liminf τt η τ μ,1 2 + liminf τt 2 F( u τ ),1 ,

Hence

t u( t ) 2 = lim τt t u τ 2 . (3.17)

Similarly, we obtain

u( t ) 1 2 = lim τt u τ 1 2 , (3.18)

η t μ,1 2 = lim τt η τ μ,1 2 . (3.19)

By the uniform continuity of the space t 1 , combining (3.17) - (3.19) and ( u, t u, η t ) C w ( [ τ,T ]; t 1 ) , we conclude that ( u, t u, η t )C( [ τ,T ]; t 1 ) .

2) (Weak lipschitz continuity) Let z i ( t ) ( i=1,2 ) be solutions to the problems (2.4) - (2.7) satisfying z i ( τ ) τ 1 R ( i=1,2 ). Then z ˜ =( u ˜ , t u ˜ , η ˜ t )= z 1 z 2 satisfies

t 2 u ˜ +A u ˜ +γ A θ t u ˜ + 0 μ( s )A η ˜ t ( s )ds + f 1 f 2 =0,( x,t )Ω×[ τ, ), (3.20)

u ˜ ( x,t )=0,xΩ,t>τ, (3.21)

η ˜ τ ( x,s )=0,( x,s )Ω× + ,t>τ, (3.22)

u ˜ ( x,t,τ )= u ˜ τ 1 ( x ) u ˜ τ 2 ( x ), t u ˜ ( x,t,τ )= t u ˜ τ 1 ( x ) t u ˜ τ 2 ( x ),

η ˜ τ ( x,s )= η ˜ τ 1 ( x,s ) η ˜ τ 2 ( x,s ),xΩ,tτ. (3.23)

where f i =f( u i ) for i=1,2 .

In the following estimates, we choose δ to be an arbitrarily small positive constant. Taking the inner product of Equation (3.20) with 2 A θ t u ˜ +2δ A θ u ˜ , we obtain

d dt K( u ˜ , t u ˜ , η ˜ t )+2γ t u ˜ 2 2δ t u ˜ θ 2 +2δ u ˜ 1θ 2 +2 0 μ( s ) d ds η ˜ t 1θ 2 ds +2δ 0 μ( s ) η ˜ t , u ˜ 1θ ds = j=1 2 Π j . (3.24)

where

K( u ˜ , t u ˜ , η ˜ t )=2δ t u ˜ , A θ u ˜ + t u ˜ θ 2 + u ˜ 1θ 2 +δγ u ˜ 2 + η ˜ t μ,1θ 2 ,

Π 1 =2 f( u 1 )f( u 2 ), A θ t u ˜ ,

Π 2 =2δ f( u 1 )f( u 2 ), A θ u ˜ .

Since

2| δ t u ˜ , A θ u ˜ |2 δ 2 λ 1 1 u ˜ 1θ 2 + 1 2 t u ˜ θ 2 ,

there exist constants m 2 , m 3 such that

m 2 ( u ˜ ( t ) 1θ 2 + t u ˜ ( t ) θ 2 + η ˜ t μ,1θ 2 )K( u ˜ , t u ˜ , η ˜ t ) m 3 ( u ˜ ( t ) 1θ 2 + t u ˜ ( t ) θ 2 + η ˜ t μ,1θ 2 ), (3.25)

where m 2 =min{ 1 2 ,12 δ 2 λ 1 1 } and m 3 =max{ 3 2 ,1+2 δ 2 λ 1 1 +δγ λ 1 θ1 } .

By the Interpolation Theorem, we deduce that

| Π 1 |2 Ω | f( u 1 )f( u 2 ) || A θ t u ˜ |dx 2 C 1 Ω ( 1+ | u 1 | p1 + | u 2 | p1 )| u ˜ || A θ t u ˜ |dx 2 C 1 ( Ω ( 1+ | u 1 | p1 + | u 2 | p1 ) p+1 p1 dx ) p1 p+1 ( Ω | u ˜ | p+1 dx ) 1 p+1 ( Ω | A θ t u ˜ | p+1 dx ) 1 p+1 C 1 ( 1+ u 1 L p+1 ( Ω ) p1 + u 2 L p+1 ( Ω ) p1 )( u ˜ L p+1 ( Ω ) 2 + A θ t u ˜ L p+1 ( Ω ) 2 ) C( R, β 0 ,g, λ 1 , C 1 , C 2 )( u ˜ 1α 2 + t u ˜ 12θα 2 ) 2γ( u ˜ 1 2 + t u ˜ 2 )+C( R, β 0 , g , λ 1 , C 1 , C 2 ,γ )( u ˜ 1θ 2 + t u ˜ θ 2 ),

| Π 2 |2δ Ω | f( u 1 )f( u 2 ) || A θ u ˜ |dx 2 C 1 δ Ω ( 1+ | u 1 | p1 + | u 2 | p1 )| u ˜ || A θ u ˜ |dx 2 C 1 δ ( Ω ( 1+ | u 1 | p1 + | u 2 | p1 ) p+1 p1 dx ) p1 p+1 ( Ω | u ˜ | p+1 dx ) 1 p+1 ( Ω | A θ u ˜ | p+1 dx ) 1 p+1 C 1 δ( 1+ u 1 L p+1 ( Ω ) p1 + u 2 L p+1 ( Ω ) p1 )( u ˜ L p +1 2 + A θ u ˜ L p +1 2 ) Cδ( 1+ u 1 L p+1 ( Ω ) p1 + u 2 L p+1 ( Ω ) p1 )( u ˜ 1α 2 + A θ u ˜ 1α2θ 2 ) γ( u ˜ 1 2 + u ˜ 2 )+C( R, β 0 , g , λ 1 , C 1 , C 2 ,δ,γ ) u ˜ 1θ 2 γ( 1+ λ 1 1 ) u ˜ 1 2 +C( R, β 0 , g , λ 1 , C 1 , C 2 ,δ,γ ) u ˜ 1θ 2 ,

where we have used the Sobolev embedding V 1α L p+1 ( Ω ) for 0<α1 .

Substituting the above estimates into Equation (3.24), we get

d dt K( u ˜ , t u ˜ , η ˜ t ) C( R, β 0 , g , λ 1 , C 2 ,γ )+C( R, β 0 , g , λ 1 , C 1 , C 2 ,δ,γ,k )K( u ˜ , t u ˜ , η ˜ t ),

where

C ˜ =C( R, β 0 , g , λ 1 , C 1 , C 2 ,δ,γ,k )

C ˜ 0 =C( R, β 0 , g , λ 1 , C 2 ,γ ).

Furthermore, we can derive Equation (3.5).

3) (Dissipativity) Let K 1 ( u, t u, η t )=E( u, t u, η t )+2δ t u,u .

By virtue of the estimate

2| δ u, t u | 2 δ 2 λ 1 u 1 2 + 1 2 t u 2

and estimate (3.1), there exist constants m 4 , m 5 such that

m 4 ( u, t u, η t ) t 1 2 C 3 K 1 ( u, t u, η t ) m 5 ( u, t u, η t ) t 1 2 +C( R, β 0 , g , λ 1 , C 2 ), (3.26)

where m 4 =min{ 1 2 , 3 β 0 4 2 δ 2 λ 1 } , m 5 =max{ 1 2 , 2 δ 2 λ 1 } , C 3 = 4 β 0 λ 1 g 2 +2C( β 0 ) .

Multiplying Equation (2.4) by 2 t u+2δu and integrating over Ω , we obtain

d dt ( K 1 ( u, t u, η t )+ C 3 )+δ( K 1 ( u, t u, η t )+ C 3 )+2δγ A θ t u,u +2δ f( u ),u +2γ t u θ 2 2 δ 2 u, t u 3δ t u 2 +δ u 1 2 +2k η t ( s ) μ,1 2 +2δ 0 μ( s ) A η t ( s ),u ds =δ η t ( s ) μ,1 2 +2δ F( u ),1 +δ C 3 . (3.27)

From estimates (3.1) and (3.8) - (3.9), it follows that

δ η t ( s ) μ,1 2 +2δ F( u ),1 +δ C 3 C( R, β 0 , g , λ 1 , C 2 , C 3 ,δ ).

It is easy to see that

| 2δγ A θ t u,u |γ t u θ 2 + δ 2 γ λ 1 1θ u 1 2 ,

t u θ 2 λ 1 θ t u 2 ,

2δ f( u ),u 2δ( ( 13 β 0 ) ) u 1 2 2δ C β 0 6δ β 0 u 1 2 2δ C β 0 2δ u 1 2 ,

| 2δ 0 μ( s ) A η t ( s ),u ds | δ 4 u 1 2 +4δ η t ( s ) μ,1 2 .

Substituting the above estimates into Equation (3.27), we get

d dt ( K 1 ( u, t u, η t )+ C 3 )+δ( K 1 ( u, t u, η t )+ C 3 )+ϒ( u, t u, η t ) C( R, β 0 , g , λ 1 , C 2 , C 3 ,δ ), (3.28)

where

ϒ( u, t u, η t )=( 6δ β 0 2 δ 3 λ 1 δ 2 γ λ 1 1θ ) u 1 2 +( γ λ 1 θ 7δ 2 ) t u 2 +( 2k4δ ) η t ( s ) μ,1 2 0.

Thus, the dissipativity of the solutions to the problems (2.4) - (2.7) can be achieved.

Theorem 3.3. Assume that Conditions (M1) - (M3) hold. If z 1 =( u 1 , t u 1 , η 1 t ) and z 2 =( u 2 , t u 2 , η 2 t ) are two solutions to the problems (2.4) - (2.7) corresponding to the initial values z 1 ( τ ) and z 2 ( τ ) respectively, then for any τ<T , we have

z 1 ( t ) z 2 ( t ) t 1 2 C ˜ 2 C ˜ 1 e 2 C ˜ 3 C ˜ 1 C ˜ ( tτ ) z 1 ( τ ) z 2 ( τ ) τ 1 2 ,t[ τ,T ]. (3.29)

Proof. Let z ˜ = z 1 z 2 , then z ˜ =( u ˜ , t u ˜ , η ˜ t ) satisfies

t 2 u ˜ +A u ˜ +γ A θ t u ˜ + 0 μ( s )A η ˜ t ( s )ds + f 1 f 2 =0,( x,t )Ω×[ τ, ), (3.30)

u ˜ ( x,t )=0,xΩ,t>τ, (3.31)

η ˜ τ ( x,s )=0,( x,s )Ω× + ,t>τ, (3.32)

u ˜ ( x,t,τ )= u ˜ τ 1 ( x ) u ˜ τ 2 ( x ), t u ˜ ( x,t,τ )= t u ˜ τ 1 ( x ) t u ˜ τ 2 ( x ),

η ˜ τ ( x,s )= η ˜ τ 1 ( x,s ) η ˜ τ 2 ( x,s ),xΩ,tτ. (3.33)

where f i =f( u i ) for i=1,2 .

Taking the inner product of Equation (3.57) with t u ˜ , we obtain

d dt K 2 ( u ˜ , t u ˜ , η ˜ t )+2γ t u ˜ θ 2 + 0 μ( s ) d ds η ˜ t ( s ) 1 2 ds =2 f( u 1 )f( u 2 ), t u ˜ . (3.34)

where

K 2 ( u ˜ , t u ˜ , η ˜ t )= u ˜ 1 2 + t u ˜ 2 + η ˜ t μ,1 2 .

Integrating Equation (3.34) over the interval [ τ,t ] , we get

K 2 ( u ˜ ( t ), t u ˜ ( t ), η ˜ t ( s ) )+2γ τ t t u ˜ ( r ) θ 2 dr +k τ t η ˜ r ( s ) μ,1 2 dr K 2 ( u ˜ τ , t u ˜ τ , η ˜ τ ( s ) ) τ t 2 f( u 1 )f( u 2 ), t u ˜ dr. (3.35)

For 2p p θ =3+2θ , by integration by parts, we have

2 f( u 1 )f( u 2 ), t u ˜ =2 Ω ( f( u 1 )f( u 2 ) ) t u ˜ dx =2 Ω f ( u 1 +c( u 2 u 1 ) ) u ˜ t u ˜ dx = d dt Ω f ( u 1 +c( u 2 u 1 ) ) u ˜ 2 dx + Ω f t ( u 1 +c( u 2 u 1 ) ) u ˜ 2 dx ,

where c[ 0,1 ] . Hence

2 τ t Ω ( f( u 1 )f( u 2 ) ) t u ˜ dxdr = Ω f ( u 1 +c( u 2 u 1 ) ) u ˜ 2 dx | τ t + τ t Ω f t ( u 1 +c( u 2 u 1 ) ) u ˜ 2 dxdr. (3.36)

By virtue of Hölder inequality, we obtain

Ω f ( u 1 ( t )+c( u 2 ( t ) u 1 ( t ) ) ) u ˜ 2 dx C 1 Ω ( 1+ | u 1 ( t ) | p1 + | u 2 ( t ) | p1 ) | u ˜ ( t ) | 2 dx C 1 ( 1+ u 1 ( t ) L p+1 p1 + u 2 ( t ) L p+1 p1 ) u ˜ ( t ) L p+1 2 C 2 ( 1+ u 1 ( t ) 1 p1 + u 2 ( t ) 1 p1 ) u ˜ ( t ) 1α 2 δ u ˜ 1 2 +C( R, β 0 , g , λ 1 , C 2 ,δ ) u ˜ ( t ) 1θ 2 δ u ˜ 1 2 + m 3 m 2 C( R, β 0 , g , λ 1 , C 2 ,δ ) e C ˜ ( tτ ) ( u ˜ τ 1θ 2 + t u ˜ τ θ 2 + η ˜ τ ( s ) μ,1θ 2 ),

where we have used the Sobolev embedding V 1α L p+1 ( Ω ) for 0<α1 and chosen δ as an arbitrarily small positive constant. Similarly, we have

Ω f ( u τ 1 +c( u τ 2 u τ 1 ) ) u ˜ τ 2 dxC( R, C 2 ) u ˜ τ 1 2 .

Again by Hölder inequality, we deduce that

τ t Ω f t ( u 1 +c( u 2 u 1 ) ) u ˜ 2 dxdr C 0 τ t Ω ( 1+| u 1 | p2 + | u 2 | p2 )( | t u 1 |+| t u 2 | ) | u ˜ | 2 dxdr C 0 τ t ( 1+ u 1 L 6 p2 + u 2 L 6 p2 )( t u 1 L 6 6p + t u 2 L 6 6p ) u ˜ L 6 2 dr C 2 τ t ( 1+ u 1 1 p2 + u 2 1 p2 )( t u 1 L 6 32θ + t u 2 L 6 32θ ) u ˜ 1 2 dr C( R, β 0 , g , λ 1 , C 1 , C 2 ) τ t ( t u 1 θ 2 + t u 2 θ 2 ) u ˜ ( r ) 1 2 dr ,

where p2 6 + 6p 6 + 1 3 =1 and 6 6p 6 32θ for p3+2θ .

Substituting the above estimates into Equation (3.36), we get

K 2 ( u ˜ ( t ), t u ˜ ( t ), η ˜ t ( s ) )+2γ τ t t u ˜ ( r ) θ 2 dr K 2 ( u ˜ τ , t u ˜ τ , η ˜ τ ( s ) )+δ u ˜ 1 2 + m 3 m 2 C( R, β 0 , g , λ 1 , C 2 ,δ ) e C ˜ ( tτ ) ( u ˜ τ 1θ 2 + t u ˜ τ θ 2 + η ˜ τ ( s ) μ,1θ 2 ) +C( R, C 2 ) u ˜ ( τ ) 1 2 +C( R, β 0 , g , λ 1 , C 1 , C 2 ) τ t ( t u 1 θ 2 + t u 2 θ 2 ) u ˜ ( r ) 1 2 dr ,

which implies that

( m 4 δ )( u ˜ ( t ) 1 2 + t u ˜ ( t ) 2 + η ˜ t ( s ) μ,1 2 ) m 3 m 2 λ 1 θ C( R, β 0 , g , λ 1 , C 2 ,δ ) e C ˜ ( tτ ) ( u ˜ τ 1 2 + t u ˜ τ 2 + η ˜ τ ( s ) μ,1 2 ) +C( R, C 2 )( u ˜ τ 1 2 + t u ˜ τ 2 + η ˜ τ ( s ) μ,1 2 )+ m 5 ( u ˜ τ 1 2 + t u ˜ τ 2 + η ˜ τ ( s ) μ,1 2 ) +C( R, β 0 , g , λ 1 , C 1 , C 2 ),

τ t ( t u 1 ( r ) θ 2 + t u 2 ( r ) θ 2 )( u ˜ ( r ) 1 2 + t u ˜ ( r ) 2 + η ˜ r ( s ) μ,1 2 )dr C ˜ 2 ( u ˜ τ 1 2 + t u ˜ τ 2 + η ˜ τ ( s ) μ,1 2 ) + C ˜ 3 τ t ( t u 1 ( r ) θ 2 + t u 2 ( r ) θ 2 )( u ˜ ( r ) 1 2 + t u ˜ ( r ) 2 + η ˜ r ( s ) μ,1 2 )dr , (3.37)

where

C ˜ 1 = m 4 δ,

C ˜ 2 =C( R, C 2 )+ m 3 m 2 λ 1 θ C( R, β 0 , g , λ 1 , C 2 ,δ ) e C ˜ ( tτ ) + m 5 ,

C ˜ 3 =C( R, β 0 , g , λ 1 , C 1 , C 2 ).

By applying Gronwall lemma, we can prove inequality (3.29). Meanwhile, we also obtain the uniqueness of the solutions to the problems (2.4) - (2.7) in the space t 1 .

Based on Theorem 3.2 and Theorem 3.3, we can define the process U( t,τ ) for the problems (2.4) - (2.7) as follows:

z( t )=U( t,τ )z( τ ): τ 1 t 1 ,

and this process is continuous from τ 1 to t 1 .

4. The Existence of Time Dependent Attractor

4.1. Time-Dependent Absorbing Set in t 1

Based on Theorem 3.2, we obtain the following result.

Theorem 4.1. Under the assumptions of Theorem 3.2, if for any initial value z( τ ) B τ ( R ) τ 1 , then there exists R 0 >0 such that the process U( t,τ ) corresponding to the problems (2.4) - (2.7) possesses a time-dependent absorbing set, namely the family of sets B t = { B t ( R 0 ) } t .

4.2. A Priori Estimates

Next, we verify the compactness of the solution process U( t,τ ) to the problems (2.4) - (2.7). To this end, we establish the following a priori estimates.

Let z i ( t )=( u i ( t ), t u i ( t ), η i t ( s ) )( i=1,2 ) be solutions to the problems (2.4) - (2.7) corresponding to the initial values ( u τ i , t u τ i , η τ i ) { B τ ( R ) } τ respectively. The difference between the two solutions z ˜ ( t )= z 1 ( t ) z 2 ( t )=( ω( t ), t ω( t ), ζ t ( s ) ) satisfies the following equations:

t 2 ω+Aω+γ A θ t ω+ 0 μ( s )A ζ t ( s )ds + f 1 f 2 =0,( x,t )Ω×[ τ, ), (4.1)

ω( x,t )=0,xΩ,t>τ, (4.2)

ζ τ ( x,s )=0,( x,s )Ω× + ,t>τ, (4.3)

ω( x,t,τ )= u τ 1 ( x ) u τ 2 ( x ), t ω( x,t,τ )= t u τ 1 ( x ) t u τ 2 ( x ),

ζ τ ( x,s )= η τ 1 ( x,s ) η τ 2 ( x,s ),xΩ,tτ. (4.4)

where f i =f( u i ) for i=1,2 .

Define

H( t )= ω( t ) 1 2 + t ω( t ) 2 + ζ t ( s ) μ,1 2 .

We shall carry out the a priori estimates in the following four steps.

Step 1. Multiply Equation (4.1) by 2 t ω and integrate over [ s,t ]×Ω , we obtain

H( t )H( s )+2γ s t Ω | A θ 2 t ω( r ) | 2 dxdr + s t 0 μ( s ) d ds ζ r ( s ) 1 2 dsdr =2 s t Ω ( f( u 1 )f( u 2 ) ) t ω( r )dxdr , (4.5)

where Tst .

From

0 μ( s ) d ds η t ( s ) 1 2 ds k η t ( s ) μ,1 2 ,

s t Ω | A θ 2 t ω( r ) | 2 dxdr λ 1 θ s t Ω | t ω( r ) | 2 dxdr ,

there exists a constant m 8 such that

H( t )H( s )+ m 6 ( s t t ω( r ) 2 dr + s t ζ r ( s ) μ,1 2 dr ) 2 s t Ω ( f( u 1 )f( u 2 ) ) t ω( r )dxdr , (4.6)

where m 6 =min{ 2γ λ 1 θ ,k } .

Then

s t t ω( r ) 2 dr + s t ζ r ( s ) μ,1 2 dr 1 m 6 H( T ) 2 m 6 s t Ω ( f( u 1 )f( u 2 ) ) t ω( r )dxdr . (4.7)

Step 2. Multiply Equation (4.1) by ω and integrate over [ T,t ]×Ω , we get

Ω t ω( t )ω( t )dx Ω t ω( T )ω( T )dx + γ 2 ω( t ) θ 2 γ 2 ω( T ) θ 2 + T t ω( r ) 1 2 dr + T t 0 μ( s ) A ζ r ,ω( r ) dsdr + T t Ω ( f( u 1 )f( u 2 ) )ω( r )dxdr = T t t ω( r ) 2 dr . (4.8)

By virtue of (4.4) and (4.5), we have

T t H( r )dr = T t ( ω( r ) 1 2 + t ω( r ) 2 + ζ r μ,1 2 )dr 1 m 6 H( T ) 2 m 6 s t Ω ( f( u 1 )f( u 2 ) ) t ω( r )dxdr + Ω t ω( T )ω( T )dx Ω t ω( t )ω( t )dx + γ 2 ω( T ) θ 2 + T t t ω( r ) 2 dr T t ζ r ,ω( r ) μ,1 dr T t Ω ( f( u 1 )f( u 2 ) )ω( r )dxdr .

Step 3. Integrate Equation (4.6) with respect to s over [ T,t ] , we obtain

H( t )( tT ) T t H( s )ds 2 T t s t Ω ( f( u 1 )f( u 2 ) ) t ω( r )dxdrds 1 m 6 H( T )+ γ 2 ω( T ) θ 2 + Ω t ω( T )ω( T )dx + T t t ω( s ) 2 ds Ω t ω( t )ω( t )dx T t ζ s ,ω( s ) μ,1 ds 2 T t s t Ω ( f( u 1 )f( u 2 ) ) t ω( r )dxdrds 2 m 6 s t Ω ( f( u 1 )f( u 2 ) ) t ω( s )dxds T t Ω ( f( u 1 )f( u 2 ) )ω( s )dxds .

Step 4. Denote

C( M )= 1 m 6 H( T )+ γ 2 ω( T ) θ 2 + Ω t ω( T )ω( T )dx , (4.9)

and

Φ T t ( ( u 1 ( T ), t u 1 ( T ), ζ 1 T ( s ) ),( u 2 ( T ), t u 2 ( T ), ζ 2 T ( s ) ) )= Ψ 1 + Ψ 2 , (4.10)

where

Ψ 1 = 1 tT ( T t t ω( s ) 2 ds Ω t ω( t )ω( t )dx T t ζ s ,ω( s ) μ,1 ds ),

Ψ 2 = 1 tT ( 2 m 6 T t Ω ( f( u 1 )f( u 2 ) ) t ω( s )dxds + T t Ω ( f( u 1 )f( u 2 ) )ω( s )dxds +2 T t s t Ω ( f( u 1 )f( u 2 ) ) t ω( r )dxdrds ).

Thus

H( t ) 1 tT C M + Φ T t ( ( u 1 ( T ), t u 1 ( T ), η 1 T ( s ) ),( u 2 ( T ), t u 2 ( T ), η 2 T ( s ) ) ). (4.11)

Next, we shall prove the asymptotic compactness of the solution process to the problems (2.4) - (2.7) by using the method of contraction functions.

4.3. Asymptotic Compactness

Theorem 4.2. If the assumptions (M1) - (M3) hold, for any fixed t , any bounded sequence { τ n } n=1 ( ,t ] with τ n as n , and any sequence { x n } n=1 τ n α , the sequence { U( t, τ n ) x n } n=1 has a convergent subsequence.

Proof. For any ε>0 and fixed t , there exists T<t such that C M tT <ε . Thanks to Theorem 2.10, we also need to show that Φ T t ( B T ( R ) ) , for every fixed t .

Let ( u n , t u n , η n t ) be the solutions to the problems (2.4) - (2.7) corresponding to the initial values ( u n 0 , u n 1 , η n 0 ) B T ( R ) . From Theorem 3.2, we know that u n 1 2 + t u n 2 + η n t μ,1 2 is bounded.

By virtue of Alaoglu theorem, Lemma 2.5 and Theorem 3.2, for any τTt , without loss of generality, we assume that

u n u weakly * in L ( [ T,t ]; V 1 ), (4.12)

t u n t uweaklyin L 2 ( [ T,t ]; L 2 ( Ω ) ), (4.13)

t 2 u n t 2 u weakly * in L ( [ T,t ]; V 2θ ), (4.14)

η n t η t weakly * in L ( [ T,t ]; L μ 2 ( + , V 1 ) ), (4.15)

t u n t uweaklyin L 2 ( [ T,t ]; V θ ), (4.16)

u n uin L p+1 ( [ T,t ]; L p+1 ( Ω ) ), (4.17)

u n uin L 2 ( [ T,t ]; V 1 ), (4.18)

u n ( t )u( t )and u n ( T )u( T )in L p+1 ( Ω ), (4.19)

t u n t uin L 2 ( [ T,t ]; L 2 ( Ω ) ). (4.20)

where we have used the compact Sobolev embedding V 1 ↪↪ L p+1 ( Ω ) for p3+2θ .

From (3.29), we obtain that

{ ( u n ( s ), t u n ( s ), η n s ) }C( [ T,t ]; s 1 )isaCauchysequence, (4.21)

and there exists ( u( s ), t u( s ), η s )C( [ T,t ]; s 1 ) such that

( u n ( s ), t u n ( s ), η n s )convergesto( u( s ), t u( s ), η s )inC( [ T,t ]; s 1 ). (4.22)

Next, we deal with each term in (4.11).

First of all, by (4.15) and (4.20), we have

lim n lim m T t t u n t u m 2 ds =0, (4.23)

lim n lim m Ω ( t u n t u m )( u n u m )dx lim n lim m ( t u n t u m u n u m ) lim n lim m ( ( t u n + t u m ) u n u m )=0, (4.24)

lim n lim m T t η n s η m s , u n u m μ,1 ds = lim n lim m T t 0 μ( l ) η n s ( l ) η m s ( l ), u n u m 1 dlds lim n lim m T t 0 μ( l ) η n s ( l ) η m s ( l ) 1 u n u m 1 dlds lim n lim m T t u n u m 1 0 μ( l ) η n s ( l ) η m s ( l ) 1 dlds lim n lim m ( T t u n u m 1 2 ds ) 1 2 ( T t ( 0 μ( l )dl )( 0 μ( l ) η n s ( l ) η m s ( l ) 1 2 dl )ds ) 1 2 k 0 ( T t u n u m 1 2 ds ) 1 2 ( T t η n s ( l ) η m s ( l ) μ,1 2 ds ) 1 2 =0. (4.25)

Combining (4.23) - (4.25), we get

lim n lim m Ψ 1 =0. (4.26)

Secondly, from (1.4) and (4.18), we obtain

lim n lim m T t Ω ( f( u n )f( u m ) )( u n u m )dxds C lim n lim m T t Ω ( 1+ | u n | p1 + | u m | p1 ) | u n u m | 2 dxds C lim n lim m T t ( 1+ u n L p+1 p1 + u m L p+1 p1 ) u n u m L p+1 2 ds C lim n lim m T t ( 1+ u n 1 p1 + u m 1 p1 ) u n u m 1 2 ds C lim n lim m T t u n u m 1 2 ds =0. (4.27)

It is easy to see that

T t Ω ( f( u n )f( u m ) )( t u n t u m )dxds = T t Ω f( u n ) t u n dxds + T t Ω f( u m ) t u m dxds T t Ω f( u m ) t u n dxds T t Ω f( u n ) t u m dxds = Ω F( u n ( t ) )dx Ω F( u n ( T ) )dx + Ω F( u m ( t ) )dx Ω F( u m ( T ) )dx T t Ω f( u m ) t u n dxds T t Ω f( u n ) t u m dxds . (4.28)

By (1.4) and the compact embedding V 1 ↪↪ L p+1 ( Ω ) , we have

| Ω ( F( u n ( t ) )F( u( t ) ) )dx | Ω | f( u( t )+ϑ( u n ( t )u( t ) ) ) || u n ( t )u( t ) |dx C Ω ( 1+ | u n ( t ) | p + | u( t ) | p )| u n ( t )u( t ) |dx C( 1+ u n ( t ) L p+1 p + u( t ) L p+1 p ) u n ( t )u( t ) L p+1 Cε. (4.29)

Because f( u n ) L 2 ( [ τ,T ]; V γ ) and t u m L 2 ( [ τ,T ]; V γ ) as n , m , we have

lim n lim m T t f( u n ), t u m ds = lim n T t f( u n ), t u ds = T t f( u ), t u ds = Ω F( u( t ) )dx Ω F( u( T ) )dx .

Similarly, we have

lim n lim m T t f( u m ), t u n ds = Ω F( u( t ) )dx Ω F( u( T ) )dx .

Therefore,

lim n lim m T t Ω ( f( u n )f( u m ) )( t u n t u m )dxds =0. (4.30)

For each fixed t , the term | s t Ω ( f( u n )f( u m ) )( t u n t u m )dxdr | is bounded. Then, thanks to Lebesgue Dominated Convergence Theorem, we get

lim n lim m T t s t Ω ( f( u n )f( u m ) )( t u n t u m )dxdrds = T t lim n lim m s t Ω ( f( u n )f( u m ) )( t u n t u m )dxdrds = T t 0ds =0. (4.31)

Hence, from (4.28) - (4.31), we obtain

lim n lim m Ψ 2 =0. (4.32)

In conclusion, we have Φ T t ( ( u 1 ( T ), t u 1 ( T ) ),( u 2 ( T ), t u 2 ( T ) ) )( B T ( R ) ) .

4.4. The Time-Dependent Attractors

Theorem 4.3. If the assumptions of Theorem 4.2 hold, then the dynamical system ( U( t,τ ), t 1 ) corresponding to the problems (2.4) - (2.7) possesses an invariant time-dependent attractor A A = { A t } t .

Proof. According to Theorems 3.2, Theorems 3.3, Theorems 4.1 and Theorems 4.2, we obtain that the conclusion of Theorem 4.3 is valid.

5. Conclusion

This paper focuses on the wave equation with fading memory terms and structural damping. When the nonlinear term satisfies the critical exponential growth condition 2p3+2θ ( 1 2 <θ<1 ), we systematically analyze the dynamical behavior of the equation solutions by employing the Faedo-Galerkin approximation method, energy estimation techniques and the contraction function method. We not only rigorously prove the well-posedness, Lipschitz continuity and asymptotic compactness of the solutions to the equation, but also successfully establish the existence of time-dependent attractor in the natural energy space H 0 1 ( Ω )× L 2 ( Ω )× L μ 2 ( + ; H 0 1 ( Ω ) ) .

Funding

National Natural Science Foundation of China (Grant Nos.12561041; 11761062).

Conflicts of Interest

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

References

[1] Giorgi, C. and Pata, V. (2001) Asymptotic Behavior of a Nonlinear Hyperbolic Heat Equation with Memory. Nonlinear Differential Equations and Applications, 8, 157-171.[CrossRef]
[2] Dafermos, C.M. (1970) Asymptotic Stability in Viscoelasticity. Archive for Rational Mechanics and Analysis, 37, 297-308.[CrossRef]
[3] Pata, V. and Zelik, S. (2006) A Remark on the Damped Wave Equation. Communications on Pure & Applied Analysis, 5, 611-616.[CrossRef]
[4] Belleri, V. and Pata, V. (2001) Attractors for Semilinear Strongly Damped Wave Equations on . Discrete & Continuous Dynamical SystemsA, 7, 719-735.[CrossRef]
[5] Savostianov, A. (2015) Strichartz Estimates and Smooth Attractors for a Sub-Quintic Wave Equation with Fractional Damping in Bounded Domains. Advances in Differential Equations, 20, 495-530.[CrossRef]
[6] Wang, X. and Tian, K.H. (2023) Time-Dependent Global Attractors for Kirchhoff-Type Wave Equations with Structural Damping. Chinese Annals of Mathematics, Series A, 44, 163-198.
[7] Luo, X. and Ma, Q. (2022) The Existence of Time-Dependent Attractor for Wave Equation with Fractional Damping and Lower Regular Forcing Term. Discrete and Continuous Dynamical SystemsB, 27, 4817-4835.[CrossRef]
[8] Li, Y., Yang, Z. and Ding, P. (2020) Regular Solutions and Strong Attractors for the Kirchhoff Wave Model with Structural Nonlinear Damping. Applied Mathematics Letters, 104, Article ID: 106258.[CrossRef]
[9] Carvalho, A.N. and Cholewa, J.W. (2002) Attractors for Strongly Damped Wave Equations with Critical Nonlinearities. Pacific Journal of Mathematics, 207, 287-310.[CrossRef]
[10] Lions, J.L. (1969) Quelques Méthodes de Résolution des Problémes aux Limites Non-Linéaires. Dunod.
[11] Chueshov, I. and Lasiecka, I. (2010) Long-Time Dynamics of Second Order Evolution Equations. Springer.
[12] 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.[CrossRef]
[13] Borini, S. and Pata, V. (1999) Uniform Attractors for a Strongly Damped Wave Equation with Linear Memory. Asymptotic Analysis, 20, 263-277.[CrossRef]
[14] Pata, V. and Zucchi, A. (2001) Attractors for a Damped Hyperbolic Equation with Linear Memory. Advances in Mathematical Sciences and Applications, 11, 505-529.
[15] Conti, M., Pata, V. and Temam, R. (2013) Attractors for Processes on Time-Dependent Spaces. Applications to Wave Equations. Journal of Differential Equations, 255, 1254-1277.[CrossRef]
[16] Chepyzhov, V.V. and Vishik, M.I. (2002) Attractors for Equations of Mathematical Physics. vol 49. American Mathematical Society.
[17] Simon, J. (1986) Compact Sets in the Space . Annali di Matematica Pura ed Applicata, 146, 65-96.[CrossRef]
[18] Conti, M. and Pata, V. (2014) Asymptotic Structure of the Attractor for Processes on Time-Dependent Spaces. Nonlinear Analysis: Real World Applications, 19, 1-10.[CrossRef]
[19] Ding, T. and Liu, Y. (2014) Time-dependent Global Attractor for the Nonclassical Diffusion Equations. Applicable Analysis, 94, 1439-1449.[CrossRef]
[20] Meng, F., Yang, M. and Zhong, C. (2015) Attractors for Wave Equations with Nonlinear Damping on Time-Dependent Space. Discrete and Continuous Dynamical SystemsSeries B, 21, 205-225.[CrossRef]
[21] Chueshov, I. and Lasiecka, I. (2007) Long-Time Dynamics of Semilinear Wave Equation with Interior-Boundary Damping and Sources of Critical Exponents. Contemporary Mathematics, 426, 153-192.
[22] Chueshov, I., Lasiecka, I. and Toundykov, D. (2008) Long-Term Dynamics of Semilinear Wave Equation with Nonlinear Localized Interior Damping and a Source Term of Critical Exponent. Discrete & Continuous Dynamical SystemsA, 20, 459-509.[CrossRef]

Copyright © 2026 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.