Skip to main content

Pullback attractors for three-dimensional Navier-Stokes-Voigt equations with delays

Abstract

Our aim in this paper is to study the existence of pullback attractors for the 3D Navier-Stokes-Voigt equations with delays. The forcing term g(t,u(tρ(t))) containing the delay is sub-linear and continuous with respect to u. Since the solution of the model is not unique, which is caused by the continuity assumption, we establish the existence of pullback attractors for our problem by using the theory of multi-valued dynamical system.

1 Introduction

Let Ω R 3 be an open, bounded and connected set. We consider the following problem for three-dimensional Navier-Stokes-Voigt (NSV) equations with delays in continuous and sub-linear operators:

{ u t ν u α 2 u t + ( u ) u + p = g ( t , u ( t ρ ( t ) ) ) , in  ( τ , + ) × Ω , div u = 0 , in  ( τ , + ) × Ω , u ( x , t ) = 0 , on  ( τ , + ) × Ω , u ( τ , x ) = u 0 ( x ) , x Ω , u ( τ + t , x ) = φ ( t , x ) , t ( h , 0 ) , x Ω .
(1.1)

Here u=( u 1 (t,x), u 2 (t,x), u 3 (t,x)) is the velocity vector field, ν is a positive constant, α is a characterizing parameter of the elasticity of the fluid, p is the pressure, g is the external force term which contains memory effects during a fixed interval of time of length h>0, ρ(t) is an adequate given delay function, u 0 is the initial velocity field at the initial time τR, φ is the initial datum on the interval (h,0).

Equation (1.1) with α=0 becomes the classical three-dimensional Navier-Stokes (NS) equation. In the past decades, many authors [16] investigated intensively the classical three-dimensional incompressible NS equation. For the sake of direct numerical simulations for NS equations, the NSV model of viscoelastic incompressible fluid has been proposed as a regularization of NS equations.

Equation (1.1) governs the motion of a Klein-Voigt viscoelastic incompressible fluid. Oskolkov [7] was the first to introduce the system which gives an approximate description of the Kelvin-Voigt fluid (see, e.g., [8, 9]). In 2010, Levant et al. [10] investigated numerically the statistical properties of the Navier-Stokes-Voigt model. Kalantarov and Titi [11] studied a global attractor of a semigroup generated by equation (1.1) for the autonomous case. Recently, Luengo et al. [12] obtained asymptotic compactness by using the energy method, and they further got the existence of pullback attractors for three-dimensional non-autonomous NSV equations.

Let us recall some related results in the literature. Yue and Zhong [13] studied the long time behavior of the three-dimensional NSV model of viscoelastic incompressible fluid for system (1.1) by using a useful decomposition method. The authors in [12, 14] deduced the existence of -pullback attractors for 3D non-autonomous NSV equations using the energy method. As we know from [15], delay terms appear naturally. In recent years, Caraballo and Real [1618] developed a fruitful theory of existence, uniqueness, stability of solutions and global attractors for Navier-Stokes models including some hereditary characteristics such as constant, variable delay, distributed delay, etc. However, our present problem has no uniqueness of solutions. To overcome the difficulty, we may cite the results by Ball [19] and by Marín-Rubio and Real [20]. Gal and Medjo [21] proved the existence of uniform global attractors for a Navier-Stokes-Voigt model with memory. As commented before, in comparison with three-dimensional Navier-Stokes equations, there is no regularizing effect. Our result of this paper is to establish the existence of pullback attractors for three-dimensional NSV equations in H 0 1 when the external forcing term g(t,u)H and the function g(t,u) is continuous with respect to u. Another difficulty is to obtain that the multi-valued processes are asymptotically compact. In 2007, Kapustyan and Valero [22] presented a method suitable for verifying the asymptotic compactness. The authors [20] applied this method to 2D Navier-Stokes equations with delays in continuous and sub-linear operators. We shall apply the energy method to prove that our multi-valued processes are asymptotically compact by making some minor modifications caused by the term α 2 Δ u t in (1.1).

This paper is organized as follows. In Section 2, we recall briefly some results on the abstract theory of pullback attractors. In Section 3, we introduce some abstract spaces necessary for the variational statement of the problem and give the proof of the global existence of solutions. In Section 4, we consider the asymptotic behavior of problem (1.1).

2 Basic theory of pullback attractors

By using the framework of evolution processes, thanks to [20, 23, 24], we now briefly recall some theories of pullback attractors and the related results. On the one hand, we have to overcome some difficulties caused by multi-valued processes. On the other hand, since our model is non-autonomous, we should use the related results for classical multi-valued processes in [23, 24], but which are not completely adapted to our model.

Let (X,d) be a metric space, and let P(X) be the class of nonempty subsets of X. As usual, we denote by dist X ( B 1 , B 2 ) the Hausdorff semi-distance in X between B 1 X and B 2 X, i.e.,

dist X ( B 1 , B 2 )= sup x B 1 inf y B 2 d X (x,y)for  B 1 , B 2 X,

where d X (x,y) denotes the distance between two points x and y in X.

We now formulate an abstract result in order to establish the existence of pullback attractors for the multi-valued dynamical system associated with (1.1).

Definition 2.1 A multi-valued process U is a family of mappings U(t,τ):XP(X) for any pair τt of real numbers such that

U(t,τ)xU(t,r)U(r,τ)x,xX,τrt.

If the above relation is not only an inclusion but also an equality, we say that the multi-valued process is strict. For example, the relation generalized by 3D non-autonomous NSV equation (see, e.g., [12]) is an equality, while the relation generalized by 3D NS equations (see, e.g., [22]) is strict.

Definition 2.2 Suppose that D ˆ 0 = { D 0 ( t ) } t R P(X) is a family of sets. A multi-valued process U is said to be D ˆ 0 -asymptotically compact if for any tR, any sequences { τ n } n = 1 with τ n , x n D 0 ( τ n ), and ξ n U(t, τ n ) x n , the sequence { ξ n } is relatively compact in X.

Lemma 2.1 If a multi-valued process U is D ˆ 0 -asymptotically compact, then the sets Λ( D ˆ 0 ,t) are nonempty compact subsets of X, where

Λ( D ˆ 0 ,t)= s t τ s U ( t , τ ) D 0 ( τ ) ¯ X ,tR.

Furthermore, Λ( D ˆ 0 ,t) attracts in a pullback sense to D ˆ 0 at time t, i.e.,

lim τ dist ( U ( t , τ ) D 0 ( τ ) , Λ ( D ˆ 0 , t ) ) =0.

Indeed, it is the minimal closed set with this property.

Definition 2.3 The family of subsets D ˆ 0 = { D 0 ( t ) } t R is said to be pullback-absorbing with respect to a multi-valued process U if for every tR and all bounded subset B of X, there exists a time τ(t,B)t such that

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

Lemma 2.2 Let the family of sets D ˆ 0 = { D 0 ( t ) } t R be pullback-absorbing for the multi-valued process U, and let U be D ˆ 0 -asymptotically compact. Then, for any bounded sets B of X, it holds that

lim τ dist ( U ( t , τ ) B , Λ ( D ˆ 0 , t ) ) =0.

Definition 2.4 Suppose that U is a multi-valued process. A family A= { A ( t ) } t R P(X) is said to be a pullback attractor for a multi-valued process U if the set A(t) is compact for any tR and attracts at time t to any bounded set BX in a pullback sense, i.e.,

lim τ dist ( U ( t , τ ) B , A ( t ) ) =0.

We can see obviously that a pullback attractor does not need to be unique. However, it can be considered unique in the sense of minimal, that is, the minimal closed family with such a property. In this sense, we obtain the following property and the existence of pullback attractors.

Lemma 2.3 [25]

Assume that U is a multi-valued process, and U is D ˆ 0 -asymptotically compact and a family D ˆ 0 is pullback-absorbing for U. Then, for any tR and any bounded subset B of X, the set

Λ(B,t)= s t τ s U ( t , τ ) B ¯ X ,tR,

is a nonempty compact subset contained in Λ( D ˆ 0 ,t), which attracts to B in a pullback sense. In fact, Λ( D ˆ 0 ,t) defined above is the minimal closed set with this property.

Furthermore, for any bounded set B, the set A(t)= B Λ ( B , t ) ¯ X is a pullback attractor. From Definitions 2.3 and 2.4, it is easy to see that A(t)Λ( D ˆ 0 ,t).

If there exists a time TR such that t T D 0 (t) is bounded, for bounded B, then

A(t)= B Λ ( B , t ) ¯ X =Λ( D ˆ 0 ,t),tT.

As we know, for the single-valued processes, the continuity of processes provides invariance; while in the multi-valued processes, the upper semi-continuity (defined below) provides negatively invariance of the omega limit sets Λ(B,t) and the attractor. The following is the definition of the upper semi-continuity of the multi-valued processes (see in [22]).

Definition 2.5 Let U be a multi-valued process on X. It is said to be upper semi-continuous if for all tτ, the mapping U(t,τ) is upper semi-continuous from X into , that is to say, given a converging sequence x n x, for some sequence { y n } such that y n U(t,τ) x n for all n, there exists a subsequence of { y n } converging in X to an element of U(t,τ)x.

Lemma 2.4 [22]

Assume that a multi-valued process U and a family D ˆ 0 . If, in addition, U(t,τ) is D ˆ 0 -asymptotically compact and upper semi-continuous, then the family { Λ ( D ˆ 0 , t ) } t R is negatively invariant, i.e.,

Λ( D ˆ 0 ,t)U(t,τ)Λ( D ˆ 0 ,τ),tτ,

where B is a bounded set B of X. The family { Λ ( B , t ) } t R is also negatively invariant and the family A(t) defined in Lemma  2.3 is also negatively invariant.

Lemma 2.5 [22]

Given a universe , if a multi-valued process U is -asymptotically compact, then, for any tR and for any D ˆ D, the omega limit set Λ( D ˆ ,t) is a nonempty compact set of X that attracts to D ˆ at time t in a pullback sense. Indeed, it is the minimal closed set with this property. If, in addition, the multi-valued process U is upper semi-continuous, then { Λ ( D ˆ , t ) } t R is negatively invariant.

From the above lemmas, we obtain the following results which are rather similar to Theorem 3 in [20]. We only sketch it here.

Lemma 2.6 Suppose that is a universe and D ˆ 0 is pullback -absorbing for a multi-valued process U, which is also D ˆ 0 -asymptotically compact. Then the results in Lemma  2.5 hold. Furthermore, the family A D = { A D ( t ) } t R , where A D (t)=Λ( D ˆ 0 ,t), and the following results hold:

  1. (1)

    For each tR, the set A D (t) defined above is compact.

  2. (2)

    A D attracts pullback to any D ˆ D.

  3. (3)

    Suppose U is upper semi-continuous, then A D is negatively invariant.

  4. (4)

    A D (t)=Λ( D ˆ 0 ,t)= D ˆ D Λ ( D ˆ , t ) ¯ .

  5. (5)

    Assume D ˆ 0 D, the minimal family of closed sets A D attracts pullback to elements of .

  6. (6)

    Assume D ˆ 0 D, each D 0 (t) is closed and the universe is inclusion-closed, then A D D and it is the only family of which satisfies the above properties (1), (2) and (3).

  7. (7)

    If contains the families of fixed bounded sets, then A= { A ( t ) } t R defined in Lemma  2.3 is the minimal pullback attractor of bounded sets, and A(t) A D ( t ) for each tR. In addition, if there exists some TR such that t T D 0 (t) is bounded, then A(t)= A D (t) for all tT.

3 Introduction to some abstract spaces and the existence of solutions

We first recall some notations about the function spaces which will be used later to discuss the regularity of pullback attracting sets. Let us consider the following abstract space:

V= { u ( C 0 ( Ω ) ) 3 : div u = 0 } .

The symbols H, V denote the closures of in L 2 ( Ω ) 3 , H 0 1 ( Ω ) 3 , respectively. In other words, H= the closure of in ( L 2 ( Ω ) ) 3 with the norm || and the inner product (,), where for u,v ( L 2 ( Ω ) ) 3 ,

(u,v)= i = 1 2 Ω u i (x) v i (x)dx.

V= the closure of in ( H 0 1 ( Ω ) ) 3 with the norm associated to the inner product ((,)), where for u,v ( H 0 1 ( Ω ) ) 3 ,

( ( u , v ) ) = i , j = 1 2 Ω u j x i v j x i dx.

We shall use v V to denote the norm of V . The value of a functional from V on an element from V is denoted by brackets ,. It follows that VHH V , and the injections are dense and compact.

Define Au=Pu for all uD(A), where P is the ortho-projector from ( L 2 ( Ω ) ) 3 onto H. Considering the properties of the operator A, we have A:V V as

Au,v:= ( ( u , v ) ) ,u,vV.

We define

b(u,v,w)= i , j = 1 3 Ω u i v j x i w j dx,

for every function u, v, w:Ω R 3 , and the operator B:V×V V as

B ( u , v ) , w =b(u,v,w),u,v,wV.

Obviously, b(u,v,w) is a continuous trilinear form such that

| b ( u , v , w ) | C 1 uvw,u,v,wV,

which yields

B ( u , v ) V C 1 uv,u,vV.

Moreover, b and B satisfy the following:

b ( u , v , w ) = b ( u , w , v ) , u , v , w V , b ( u , v , v ) = 0 , u , v , w V .

Now, we make some assumptions. The given delay function ρ satisfies ρ C 1 ([0,+);[0,h]), and there is a constant ρ 0 independent of t satisfying

ρ (t) ρ 0 <1t0,
(3.1)

where ρ = d ρ d t .

Moreover, we assume that g:[τ,+)×HH satisfies the following assumptions:

(H1) g(,u):[τ,+)H is measurable for all uV.

(H2) For all tτ, g(t,):HH is continuous.

(H3) There exist two functions γ,β:[τ,+)[0,+). The above γ L p (τ,T) and β L 1 (τ,T) for all T>τ, for 1p+, such that g ( t , u ) V 2 γ(t) u 2 +β(t), tτ, uV.

As to the initial datum, we assume

(H4) u 0 V, and φ L 2 q (h,0;V), where 1 p + 1 q =1.

Next, we shall consider the solution of (1.1).

{ u ( t ) + α 2 Au ( t ) + τ t ( ν Au ( s ) + B ( u ( s ) ) ) d s = u 0 + α 2 Au 0 + τ t g ( s , u ( s ρ ( s ) ) ) d s , t τ , u ( τ + t ) = φ ( t ) , a.e.  t ( h , 0 ) .
(3.2)

Definition 3.1 It is said that u is a weak solution to (1.1) if u belongs to u L 2 q (τh,T;V) L (τ,T;V) for all tτ such that u(τ+t) coincides with φ(t) in (h,0) and satisfies equation (3.2) in V for all tτ.

If u is a solution of (3.2), then it is easy to get u(t)+ α 2 Au(t) L 2 (τ,T; V ), and d d t (u(t)+ α 2 Au(t)) L 1 (τ,T; V ) for all T>τ. From the property of the operator A, we have u(t)+ α 2 Au(t)C(τ,+; V ). On the other hand, reasoning as in [12], we have uC([τ,+);V).

Now, we define a functional γ ˜ (t)=γ ζ 1 (t), where ζ:[τ,+)[ρ(τ),+) is a differentiable and nonnegative strictly increasing function given by ζ(s)=sρ(s). We have

τ t g ( t , u ( t ρ ( t ) ) ) V 2 d t τ T γ ( t ) u ( t ρ ( t ) ) 2 d t + τ T β ( t ) d t 1 1 ρ 0 τ ρ ( τ ) T ρ ( T ) γ ˜ ( t ) u ( t ) 2 d t + τ T β ( t ) d t 1 1 ρ 0 ( ρ ( τ ) 0 γ ˜ ( t + τ ) φ 2 d t + τ T γ ˜ ( t ) u ( t ) 2 d t ) + τ T β ( t ) d t .
(3.3)

From the above analysis, taking into account γ ˜ (t) L P (ρ(τ),T) for all T>τ, we obtain that g(t,u(tρ(t))) L 2 (τ,T; V ). Hence, it is clear that u is a weak solution to (1.1) if uC([τ,+);V), u L 2 (τ,T;V) for all T>τ, and satisfies the energy equality

1 2 d d t ( | u ( t ) | 2 + α 2 u ( t ) 2 ) +ν u ( t ) 2 = g ( t , u ( t ρ ( t ) ) ) , u ( t ) ,a.e. t>τ

in the distributional sense on (τ,+).

Theorem 3.1 Suppose that (H1)-(H4) hold. Then there exists a global solution u to (3.2).

Proof We shall prove the result by the Faedo-Galerkin scheme and compactness method. For convenience and without loss of generality, we set the initial time τ=0. As to different value τ, we only proceed by translation.

Consider the Hilbert basis of H formed by the eigenfunctions { v k } of the above operator A, i.e., A v k = λ k v k . In fact, these elements allow to define the operator P m v= k = 1 m ( v k ,v) v k , which is the orthogonal projection of H and V in V m :=span[ v 1 ,, v m ] with their respective norms.

Denote u m (t)= k = 1 m η m k (t) v k , where η m k (t)=( u m (t), v k ), k=1,2,,m, are unknown real functions satisfying the finite-dimensional problem

{ ( u m , v k ) + α 2 ( Au m ( t ) , v k ) + ν 0 t Au m ( s ) , v k d s + 0 t B ( u m ( s ) , u m ( s ) ) , v k d s = ( u 0 , v k ) + α 2 ( Au 0 , v k ) + 0 t ( g ( s , u m ( s ρ ( s ) ) ) , v k ) d s , t > 0 , 1 k m , u m ( t ) = φ m ( t ) , a.e.  t ( h , 0 ) ,
(3.4)

with φ m (t)= P m φ(t). From [12], we can obtain the local well-posedness of this finite-dimensional delay problem. The following provides estimates which imply that the solutions are well defined in the whole [0,+).

By (3.4), we obtain

d d t ( ( u m , v k ) + α 2 ( Au m ( t ) , v k ) ) + ν Au m ( t ) , v k + B ( u m ( t ) , u m ( t ) ) , v k = ( g ( t , u m ( t ρ ( t ) ) ) , v k ) , t > 0 , 1 k m .
(3.5)

Multiplying (3.5) by ( u m (t), v k ), summing from k=1 to k=m, and using the properties of the operator b, we easily get

d d t ( | u m ( t ) | 2 + α 2 u m ( t ) 2 ) +2ν u m ( t ) 2 =2 g ( t , u m ( t ρ ( t ) ) ) , u m ( t ) .
(3.6)

Observing (H3) and using the Young inequality, we have

2 g ( t , u m ( t ρ ( t ) ) ) , u m ( t ) 2 g ( t , u ( t ρ ( t ) ) ) V u m ( t ) 2 ( γ 1 2 ( t ) u m ( t ρ ( t ) ) + β 1 2 ( t ) ) u m ( t ) = 2 γ 1 2 ( t ) u m ( t ρ ( t ) ) u m ( t ) + 2 β 1 2 ( t ) u m ( t ) ν ( 1 ρ 0 ) u m ( t ρ ( t ) ) 2 + γ ( t ) ν ( 1 ρ 0 ) u m ( t ) 2 + β ( t ) + u m ( t ) 2 = ν ( 1 ρ 0 ) u m ( t ρ ( t ) ) 2 + ( 1 + γ ( t ) ν ( 1 ρ 0 ) ) u m ( t ) 2 + β ( t ) .
(3.7)

Considering the above inequality in (3.6) and observing that

0 t u m ( s ρ ( s ) ) 2 d s 1 1 ρ 0 ρ ( 0 ) t ρ ( t ) u m ( s ) 2 d s 1 1 ρ 0 ( h 0 φ ( s ) 2 d s + 0 t u m ( s ) 2 d s ) ,
(3.8)

we easily get

| u m ( t ) | 2 + α 2 u m ( t ) 2 + 2 ν 0 t u m ( s ) 2 d s | u 0 | 2 + α 2 ( t ) u 0 2 + 0 T β ( s ) d s + ν h 0 φ ( s ) 2 d s + 0 t ( 1 + ν + γ ( t ) ν ( 1 ρ 0 ) ) u m ( s ) 2 d s .
(3.9)

From the above inequality and the Gronwall inequality, one has that { u m } is bounded in L 2 (0,T;V), also in L (0,T;V), for any T>0. Observing (3.5), we know that the sequence { d u m d t } m 1 is bounded in L 2 (0,T; V ) for all T>0. The reason is the same as in [12].

By the compactness of the injection of V into H, using the above estimates and the Ascoli-Arzelà theorem, we deduce that there exist a subsequence { u m } m 1 (we relabel the same) and u W 1 , 2 (0,T;V) L 2 q (h,T;V) for any T>0 with u=φ in (h,0). Recalling (3.9) and the above analysis, we have

u m uweakly-star in  L (0,T;V),
(3.10)

u m u strongly in C(0,T;H) for all T>0, and

u m (t)u(t),a.e. t>0 in V.
(3.11)

By the properties of operator A and (3.10), we deduce that Au m Au weakly in L 2 (0,T; V ). Reasoning as in [26] on page 76, we deduce that B( u m )B(u) weakly in L 2 (0,T; V ) for any T>0.

On the other hand, observing that by (3.11), (H2) and the hypothesis on ρ in (3.1), for any time t, we conclude for any T>0

g ( t , u m ( t ρ ( t ) ) ) g ( t , u ( t ρ ( t ) ) ) in V a.e. 0<t<T
(3.12)

and as { u m } is bounded in the space L (0,T; V ) obtained from (3.9), by (H3), we deduce

g ( t , u m ( t ρ ( t ) ) ) V 2 { β ( t ) + γ ( t ) , t > ρ ( t ) , β ( t ) + γ ( t ) φ ( t ρ ( t ) ) 2 , t < ρ ( t ) ,
(3.13)

where C= sup m 1 u m L ( 0 , T ; V ) . Thus we can easily derive from (3.12) and (3.13)

g ( t , u m ( t ρ ( t ) ) ) g ( t , u ( t ρ ( t ) ) ) in  L 2 ( 0 , T ; V ) .
(3.14)

From the above discussion, passing to the limit, we prove that u is a global solution of (3.2) in the sense of Definition 3.1. □

Remark 3.1 We can obtain the uniqueness if there are additional assumptions on the forcing term g. For example, if we suppose that g satisfies (H1), (H3), for an arbitrary T>τ and C T = C 1 max t [ τ , T ] (u(t),v(t))>0, then for the solutions u, v,

g ( t , u ) g ( t , v ) V (ν+ C T )uv,

then we obtain the uniqueness of solutions.

4 Existence of pullback attractors

In this section, we discuss the existence of pullback attractors for the 3D Navier-Stokes-Voight equations with delays in continuous and sub-linear operators. At first, we propose the assumptions for g given in Section 3:

(H1) For all uV, g(,u):RH is measurable.

(H2) For all tR, g(t,):HH is continuous.

(H3) There are two nonnegative functions γ(t),β(t):R[0,+) with γ(t) L l o c p (R) for some 1p+ and β(t) L l o c 1 (R) such that for any uV,

g ( t , u ) V 2 γ(t) u 2 +β(t)for any tR.

To construct a multi-valued process, we introduce symbols C V =C([h,0];V) and S V 2 q =V× L 2 q (h,0;V), where 1 p + 1 q =1 as two phase spaces. Let D(τ, u 0 ,φ) denote the set of global solutions of (1.1) in [τ,+) and the initial datum ( u 0 ,φ) S V 2 q .

By Theorem 3.1, we know there exists a solution to problem (3.2) although we have no discussion on the uniqueness of solutions to problem (3.2). We may define two strict processes, ( C V ,{U(,)}) as

U(t,τ)φ= { u ( t ) : u D ( τ , φ ( 0 ) , φ ) } for any φ C V ,

and ( S V 2 q ,{U(,)}) as

U(t,τ)( u 0 ,φ)= { u ( t ) : u D ( τ , u 0 , φ ) } for any ( u 0 ,φ) S V 2 q .

Considering the regularity of the problem, the asymptotic behavior of the two processes shall be the same, as we shall see in what follows.

In order to simplify the calculation form, we introduce a function κ σ . For any σ>0, we set

κ σ (t,s)= ( d σ + d α 2 ) (ts) e d h α 2 d ( 1 ρ 0 ) s t γ(r)dr,t,sR.
(4.1)

From (4.1), we can find that

κ σ (t,s)= κ σ (0,t) κ σ (0,s),t,sR,
(4.2)

and for any σ: 0<σ< α 2 dd, then

κ σ (0,r) κ σ (0,t)+ ( d σ + d α 2 ) h,r[th,t],
(4.3)

where d=νmin{ λ 1 , 1 α 2 }.

Lemma 4.1 Suppose that (H1′)-(H3′) hold, for any initial datum ( u 0 ,φ) S V 2 q and any uD(τ, u 0 ,φ), it holds

u ( t ) 2 1 α 2 C τ e κ σ ( t , τ ) + 1 α 2 σ 1 τ t e κ σ ( t , s ) β(s)ds,
(4.4)

where C τ =d h 0 e d r φ ( r ) 2 dr+ | u 0 | 2 + α 2 u 0 2 .

Proof Let u be a solution of (3.2), so uD(τ, u 0 ,φ). Multiplying (1.1) by u(t) and using the energy equality and the Poincaré inequality, we have

d d t ( | u ( t ) | 2 + α 2 u ( t ) 2 ) + d ( | u ( t ) | 2 + α 2 u ( t ) 2 ) 2 g ( t , u ( t ρ ( t ) ) ) V u ( t ) V ,

where d=νmin{ λ 1 , 1 α 2 }. Thus

d d t [ e d t ( | u ( t ) | 2 + α 2 u ( t ) 2 ) ] 2 e d t g ( t , u ( t ρ ( t ) ) ) V u ( t ) V 2 e d t ( γ 1 2 ( t ) u ( t ρ ( t ) ) + β 1 2 ( t ) ) u ( t ) 2 e d t γ 1 2 ( t ) u ( t ) u ( t ρ ( t ) ) + 2 e d t β 1 2 ( t ) u ( t ) C 0 1 e d t u ( t ρ ( t ) ) 2 + ( C 0 γ ( t ) + σ ) e d t u ( t ) 2 + σ 1 e d t β ( t ) ,
(4.5)

where we have denoted

C 0 = e d h d ( 1 ρ 0 ) .

Considering

τ t e d s u ( s ρ ( s ) ) 2 d s e d h 1 ρ 0 ( τ h t e d r u ( r ) d r ) = e d h 1 ρ 0 ( e d τ h 0 e d r φ ( r ) 2 d r + τ t e d r u ( r ) 2 d r ) ,
(4.6)

and integrating (4.5) from τ to t, we deduce

e d t | u ( t ) | 2 + α 2 e d t u ( t ) 2 C 0 1 τ t e d s u ( s ρ ( s ) ) 2 d s + τ t ( C 0 γ ( s ) + σ ) e d s u ( s ) 2 d s + σ 1 τ t β ( s ) d s + e d τ | u 0 | 2 + α 2 e d τ u 0 2 C 0 1 × e d h 1 ρ 0 × e d τ h 0 e d r φ ( r ) 2 d r + C 0 1 × e d h 1 ρ 0 τ t e d r u ( r ) 2 d r + τ t ( C 0 γ ( s ) + σ ) e d s u ( s ) 2 d s + σ 1 τ t e d s β ( s ) d s + e d τ | u 0 | 2 + α 2 e d τ u 0 2 = d e d τ h 0 e d r φ ( r ) 2 d r + τ t e d s ( C 0 γ ( s ) + σ + d ) u ( s ) 2 d s + σ 1 τ t e d r β ( r ) d r + e d τ | u 0 | 2 + α 2 e d τ u 0 2 ,
(4.7)

where we set

C τ =d h 0 e d s φ ( s ) 2 ds+ | u 0 | 2 + α 2 u 0 2 .

Observing the above estimates, we easily deduce

e d t | u ( t ) | 2 + α 2 e d t u ( t ) 2 e d τ C τ + τ t ( C 0 γ ( s ) + σ + d ) e d r u ( r ) 2 d r + σ 1 τ t e d r β ( r ) d r .
(4.8)

Applying the Poincaré inequality and the Gronwall inequality to (4.8), we deduce that (4.4) holds. This finishes the proof of this lemma. □

Next, we shall prove that the processes ( C V ,{U(,)}) and ( S V 2 q ,{U(,)}) defined above are pullback-absorbing. To obtain this, we propose the assumptions (H4′)

lim sup t 1 t 0 t γ(s)ds= γ ¯ [0,+),

and the relation among constants σ>0, d defined above, and α in (1.1) satisfies(H5′)

d σ + d α 2 C 0 α 2 >0,

and β(t) satisfies

(H6)

0 e κ σ ( 0 , r ) β(r)dr<+,

where the function κ σ (t,s) is given by (4.1).

Before proving that the two multi-valued processes possess pullback-absorbing sets, we introduce the definition of the two natural tempered universes which shall play the key role for our main purpose.

Definition 4.1 Suppose that R σ (t) is the collection of the sets of all functions r:R[0,+) such that

lim t e κ σ ( 0 , t ) r 2 (t)=0.

Let D S V 2 q be the class of all families D ˆ ={D(t):tR}P( S V 2 q ) such that D(t) B ¯ S V 2 q (0, r D ˆ (t)) for some r D ˆ R σ . In the same way, let D C V σ denote the class of all families D ˆ ={D(t):tR}P( C V ) satisfying D(t) B ¯ C V (0, r D ˆ (t)) for some r D ˆ R σ .

Let B 0 be any fixed bounded subset of S V 2 q . Observing that D C V σ D S V 2 q , which is inclusion-closed, by (H4) and (H5), we deduce that the family B ˆ ={B(t) B 0 ,tR} is contained in D S V 2 q σ . With regard to D C V , we use the same method and obtain a similar conclusion if B 0 is included in C V .

The following lemma provides that there exist pullback-absorbing sets for the two processes mentioned above.

Lemma 4.2 Suppose that (H1)-(H6) hold and the constants α, d, C 0 , γ ¯ satisfy d σ + d α 2 C 0 α 2 γ ¯ >0 and γ ¯ <1.

  1. (1)

    Then, for any tR and any family B ˆ ={B(t):tR}, there exits τ( B ˆ ,t)t such that any initial datum ( u 0 ,φ) S V 2 q and any uD(τ, u 0 ,φ) for any ττ( B ˆ ,t) satisfy that u(t) R V (t), where the positive continuous function R V () is given by

    R V 2 (t)=1+ 1 α 2 σ 1 τ t e κ σ ( t , s ) β(s)ds,tR.
  2. (2)

    Let D ˆ 0 ={ D 0 (t):tR} be included P( C V ) which is given by

    D 0 (t)= B ¯ C V ( 0 , R ˜ V ( t ) ) and R ˜ V (t)= max t h r t R V (r),tR.

Then the set D ˆ 0 D C V σ and is D S V 2 q σ -pullback absorbing for the process ( S V 2 q ,U). Therefore, D ˆ 0 is D C V σ -pullback absorbing for the process ( C V ,U).

Proof Since the proof is a consequence of the definition of D S V 2 q , we only sketch it here. From Lemma 4.1, (H5) and (H6) , we have

u ( t ) 2 1 α 2 C τ e κ σ ( t , τ ) + 1 α 2 σ 1 τ t e κ σ ( t , s ) β ( s ) d s = 1 α 2 C τ e κ σ ( 0 , t ) e κ ( 0 , τ ) + 1 α 2 σ 1 τ t e κ σ ( t , s ) β ( s ) d s = 1 α 2 C τ e κ σ ( 0 , t ) e ( d σ + d α 2 ) τ C 0 α 2 0 τ γ ( r ) d r + 1 α 2 σ 1 τ t e κ σ ( t , s ) β ( s ) d s C τ α 2 e κ σ ( 0 , t ) e ( d σ + d α 2 C 0 α 2 ) τ + 1 α 2 σ 1 τ t e κ σ ( t , s ) β ( s ) d s 1 + 1 α 2 σ 1 t e κ σ ( t , s ) β ( s ) d s .
(4.9)

We complete the proof of (1). Estimate (2) is a consequence of (1). □

For the two processes ( S V 2 q ,U) and ( C V ,U), they possess pullback-absorbing sets. In order to apply Lemma 2.3 to obtain the existence of pullback attractors, it is necessary to prove that the two multi-valued processes are asymptotically compact. This will be done in the following lemma.

Lemma 4.3 Suppose that the assumptions in Lemma  4.2 hold. The two processes ( C V ,U) and ( S V 2 q ,U) are D ˆ 0 -asymptotically compact.

Proof For any t 0 R , a sequence { τ n }(, t 0 2h] with τ n and a sequence { u n } with u n D( τ n , φ n (0), φ n ) with φ n D 0 ( τ n ), we shall prove that the sequence { u t 0 n } is relatively compact in C V . By the properties concerning operator b mentioned in Section 3, we deduce that

( u n ) V + α 2 ( u n ) ν u n + b ( u n , u n , ) V + g ( t , u n ( t ρ ( t ) ) ) V .

It is easy to get the above estimate which is independent of n. The sequences of { u n } and { ( u n ) } possess their subsequence, relabeled the same in suitable spaces such that there exist u L ( t 0 2h;V) and u L 2 ( t 0 h, t 0 ; V ) satisfying

{ u n u weakly-star in  L ( t 0 2 h , t 0 ; V ) , u n u weakly in  L 2 ( t 0 2 h , t 0 ; V ) , d d t u n d d t u weakly in  L 2 ( t 0 h , t 0 ; V ) , u n ( r ) u ( r ) strongly in  V  a.e.  r ( t 0 2 h , t 0 ) .
(4.10)

According to the assumptions on a function g and analogously as in Theorem 3.1, we deduce that g(r, u n (rρ(r)))g(r,u(rρ(r))) strongly in V a.e. r( t 0 h, t 0 ).

By the Lebesgue theorem and the uniform estimate of u n in L ( t 0 2h, t 0 ;V), we deduce that the function g(r, u n (rρ(r))) converges to the function g(r,u(rρ(r))) strongly. Therefore, for any t[ t 0 h, t 0 ], we have uC([ t 0 h, t 0 ];V) and

u ( t ) + α 2 Au ( t ) + t 0 h t ( ν Au ( r ) + B ( u ( r ) ) ) d r = u ( t 0 h ) + α 2 Au ( t 0 h ) + t 0 h t g ( r , u ( r ρ ( r ) ) ) d r .
(4.11)

The uniform estimate of { ( u n ) } in L 2 ( t 0 h, t 0 ; V ) implies that the sequence { u n } is equicontinuous in V for any t 0 ht t 0 . In addition, the sequence { u n } is bounded, which is independent of n in C([ t 0 h, t 0 ];V). Using the Ascoli-Arzelà theorem, we can obtain

u n ustrongly in C ( [ t 0 h , t 0 ] ; V ) .
(4.12)

From the uniform boundedness of the sequence { u n } in C([ t 0 h, t 0 ];V), for any r[ t 0 h, t 0 ], we can also obtain u n (r)u(r), weakly in V.

By the analogous argument, for any compact sequence { r n }[ t 0 h, t 0 ] and { r n }r[ t 0 h, t 0 ], we obtain that the sequence { u n ( r n )} is convergent to u(r) weakly in V. To achieve our result in Lemma 4.3, we only need to prove

u n ustrongly in C ( [ t 0 h , t 0 ] ; V ) .

The proof is slightly different from Proposition 6 in [20] or in [22]. We only sketch it here. We use a contradiction argument. Suppose that it is not true, then there would exist a value ε, a sequence (relabeled the same) { r n }[ r n h, t 0 ], and r [ t 0 h, t 0 ] with r n r satisfying u n ( r n )u( r )ε for all n1. We shall see u n ( r n )u( r ) in V. In order to achieve the last claim, because the sequence { u n ( r n )} is weakly convergent to u(r) in V, we only need the convergence of the norms above. In other words, u n ( r n )u( r ) as n.

From the weak convergence of u n ( r n ) in V, we get

u ( r ) lim n + inf u n ( r ) .

Therefore we have to check that

lim n + sup u n ( r n ) u ( r ) .
(4.13)

From the energy equality, for any t 0 hrt t 0 , we obtain

1 2 | z ( t ) | 2 + 1 2 z ( t ) 2 + ν r t z ( s ) 2 d s = 1 2 | z ( r ) | 2 + 1 2 z ( r ) 2 + r t ( g ( s , z ( s ρ ( s ) ) ) ) d s ,
(4.14)

where z=u or z= u n . For any t[ t 0 h, t 0 ], define the continuous functions J(t) and J n (t) as

J ( t ) = 1 2 | u ( t ) | 2 + 1 2 u ( t ) 2 t 0 h t ( g ( s , u ( s ρ ( s ) ) ) , u ( s ) ) d s , J n ( t ) = 1 2 | u n ( t ) | 2 + 1 2 u n ( t ) 2 t 0 h t ( g ( s , u n ( s ρ ( s ) ) ) , u n ( s ) ) d s .

By (4.14), it is clear that J and J n are non-increasing functions. By the convergence (4.10), for any t( t 0 h, t 0 ), we have that J n (t)J(t). Using the same analysis method as in [22], we can deduce that for nn( κ ε ), J n ( r n )J( r )ε, which gives (4.13) as desired. □

We can apply the technical method for any family in D S V 2 q σ in Lemma 4.3. Suppose that the assumptions in Lemma 4.3 hold. We can deduce that the processes ( C V ,U) and ( S V 2 q ,U) are D C V σ -asymptotically compact and D S V 2 q σ -asymptotically compact.

Lemma 4.4 Suppose that (H1′)-(H6′) hold. The two processes ( C V ,U) and ( S V 2 q ,U) are semi-continuous and that U(t,τ): S V 2 q P( S V 2 q ) and U(t,τ): C V P( C V ) have compact values in their respective topologies.

Proof In fact, the upper semi-continuity of the process ( S V 2 q ,U) can be obtained by similar arguments to those used for the Galerkin sequence in Theorem 3.1.

As to the process ( C V ,U), applying the same energy-procedure in Lemma 4.3, we shall obtain that in [τ,t] any set of solutions possesses a converging subsequence in this process, whence the assertion in Lemma 4.4 follows. □

According to the results in Section 2, the following two theorems shall be obtained, which are our result in this paper. Observing Lemmas 4.3 and 4.4, and applying Lemma 2.6, we obtain the following theorem.

Theorem 4.1 Suppose that (H1′)-(H6′) hold. For any tR, then there exist global pullback attractors A C V ={ A C V (t)} and A D C V σ ={ A D C V σ (t)} for the process ( C V ,U) in the universe of fixed bounded sets and in D C V σ , respectively. Moreover, they are unique in the sense of Lemma  2.5 and negatively and strictly invariant for U respectively, and the following holds:

A C V (t) A D C V σ (t).

The above theorem proves that there exist pullback attractors in the C V framework, while we shall prove that there exist pullback attractors in the S V 2 q framework in the following theorem.

Theorem 4.2 Suppose that the assumptions in Theorem  4.1 hold. For any tR, there exist global pullback attractors A S V 2 q ={ A S V 2 q (t)} and A D S V 2 q σ ={ A D S V 2 q σ (t)} for the process ( S V 2 q ,U) in the universes of fixed bounded sets and in D S V 2 q σ . They are unique in the sense of Lemma  2.5 and negatively and strictly invariant for U, respectively, and we have A S V 2 q (t) A D S V 2 q σ (t). Moreover, the relationship between the attractors for ( S V 2 q ,U) and for ( C V ,U) is as follows:

A S V 2 q (t)=f ( A C V ( t ) ) and A D S V 2 q σ (t)=f ( A D C V σ ( t ) ) ,
(4.15)

where f: C V S V 2 q is the continuous mapping defined by f(φ)=(φ(0),φ).

Proof The proof is rather similar to that of Theorem 5 in [20]. Since the regularity is different from [20], we only sketch the proof of (4.15) here.

By Theorem 3.1, we can conclude that U(t,τ) maps S V 2 q into bounded sets in C V if tτ+h, and also maps bounded sets from S V 2 q into bounded sets of C V .

Noting that A S V 2 q (t) is the minimal closed set, and using Lemma 2.5 and the above arguments, we deduce that f( A C V (t)) also attracts bounded sets in S V 2 q in a pullback sense. Therefore the inclusion A S V 2 q (t)f( A C V (t)) holds.

As to the opposite inclusion of the first identification in (4.15), for any bounded set B, it follows from the continuous injection f( C V ) S V 2 q and the attractor A C V (t)= B Λ C V ( B , t ) ¯ C V . Thus f( A C V (t))= B f ( Λ C V ( B , t ) ) ¯ C V , whence the opposite inclusion of the first identification in (4.15) holds.

Analogously, it is obvious that the second relation in (4.15) holds. □

References

  1. Caraballo T, Real J, Kloeden P: Unique strong solutions and V-attractors of a three dimensional system of globally modified Navier-Stokes equations. Adv. Nonlinear Stud. 2006, 6: 411-436.

    MathSciNet  MATH  Google Scholar 

  2. Kloeden P, Langa J, Real J: Pullback V -attractors of the 3-dimensional globally modified Navier-Stokes equations. Commun. Pure Appl. Anal. 2007, 6: 937-955.

    Article  MathSciNet  MATH  Google Scholar 

  3. Kloeden P, Valero J: The weak connectedness of the attainability set of weak solutions of the three-dimensional Navier-Stokes equations. Proc. R. Soc. Lond., Ser. A, Math. Phys. Eng. Sci. 2007, 463: 1491-1508. 10.1098/rspa.2007.1831

    Article  MathSciNet  MATH  Google Scholar 

  4. Ma T, Wang S World Scientific Series on Nonlinear Science, Series A 53. Bifurcation Theory and Applications 2005.

    Google Scholar 

  5. Ma T, Wang S AMS Monograph and Mathematical Survey Series 119. Geometric Theory of Incompressible Flows with Applications to Fluid, Dynamics 2005.

    Chapter  Google Scholar 

  6. Ma T, Wang S: Stability and Bifurcation of Nonlinear Evolutions Equations. Science Press, Beijing; 2007.

    Google Scholar 

  7. Oskolkov A: The uniqueness and solvability in the large of boundary value problems for the equations of motion of aqueous solutions of polymers. Zap. Nauč. Semin. POMI 1973, 38: 98-116.

    MathSciNet  Google Scholar 

  8. Ladyzhenskaya O: In memory of A. P. Olskolkov. J. Math. Sci. 2000, 99: 799-801. 10.1007/BF02673588

    Article  MathSciNet  Google Scholar 

  9. Oskolkov A: Theory of nonstationary flows of Kelvin-Voigt fluids. J. Math. Sci. 1985, 28: 751-785. 10.1007/BF02112340

    Article  MATH  Google Scholar 

  10. Levant B, Ramos F, Titi ES: On the statistical properties of the 3D incompressible Navier-Stokes-Voigt model. Commun. Math. Sci. 2010, 8: 277-293. 10.4310/CMS.2010.v8.n1.a14

    Article  MathSciNet  MATH  Google Scholar 

  11. Kalantarov V, Titi E: Global attractors and determining modes for the 3D Navier-Stokes-Voigt equations. Chin. Ann. Math., Ser. B 2009, 30: 697-714. 10.1007/s11401-009-0205-3

    Article  MathSciNet  MATH  Google Scholar 

  12. García-Luengo J, Marín-Rubio P, Real J: Pullback attractors for three-dimensional non-autonomous Navier-Stokes-Voigt equations. Nonlinearity 2012, 25: 905-930. 10.1088/0951-7715/25/4/905

    Article  MathSciNet  MATH  Google Scholar 

  13. Yue G, Zhong C: Attractors for autonomous and nonautonomous 3D Navier-Stokes-Voigt equations. Discrete Contin. Dyn. Syst., Ser. B 2011, B16(3):985-1002.

    Article  MathSciNet  MATH  Google Scholar 

  14. Qin, Y, Yang, X: Pullback attractors for a 3D non-autonomous Navier-Stokes-Voigt equations. Acta Math. Appl. Sin. (to appear)

  15. Manitius AZ: Feedback controllers for a wind tunnel model involving a delay: analytical design and numerical simulation. IEEE Trans. Autom. Control 1984, 29: 1058-1068. 10.1109/TAC.1984.1103436

    Article  MathSciNet  Google Scholar 

  16. Caraballo T, Real J: Navier-Stokes equations with delays. R. Soc. Lond. Proc., Ser. A, Math. Phys. Eng. Sci. 2001, 457: 2441-2453. 10.1098/rspa.2001.0807

    Article  MathSciNet  MATH  Google Scholar 

  17. Caraballo T, Real J: Asymptotic behavior for two-dimensional Navier-Stokes equations with delays. R. Soc. Lond. Proc., Ser. A, Math. Phys. Eng. Sci. 2003, 459: 3181-3194. 10.1098/rspa.2003.1166

    Article  MathSciNet  MATH  Google Scholar 

  18. Caraballo T, Real J: Attractors for 2D Navier-Stokes equations with delays. J. Differ. Equ. 2004, 205: 271-297. 10.1016/j.jde.2004.04.012

    Article  MathSciNet  MATH  Google Scholar 

  19. Ball JM: Continuity properties and global attractors of generalized semiflow and the Navier-Stokes equations. J. Nonlinear Sci. 1997, 7: 475-502. 10.1007/s003329900037

    Article  MathSciNet  MATH  Google Scholar 

  20. Marín-Rubio P, Real J: Pullback attractors for 2D-Navier-Stokes equations with delays in continuous and sub-linear operators. Discrete Contin. Dyn. Syst. 2010, 26: 989-1006.

    Article  MathSciNet  MATH  Google Scholar 

  21. Gal CG, Medjo TT: A Navier-Stokes-Voigt model with memory. Math. Methods Appl. Sci. 2013. 10.1002/mma.2771

    Google Scholar 

  22. Kapustyan AV, Valero J: Weak and strong attractors for the 3D Navier-Stokes system. J. Differ. Equ. 2007, 240: 249-278. 10.1016/j.jde.2007.06.008

    Article  MathSciNet  MATH  Google Scholar 

  23. Caraballo T, Łukaszewicz G, Real J: Pullback attractors for asymptotically compact nonautonomous dynamical systems. Nonlinear Anal. TMA 2006, 64: 484-498. 10.1016/j.na.2005.03.111

    Article  MATH  Google Scholar 

  24. Marín-Rubio P, Real J: On the relation between two different concepts of pullback attractors for non-autonomous dynamical systems. Nonlinear Anal. TMA 2009, 71: 3956-3963. 10.1016/j.na.2009.02.065

    Article  MathSciNet  MATH  Google Scholar 

  25. García-Luengo J, Marín-Rubio P, Real J: Pullback attractors in V for non-autonomous 2D-Navier-Stokes equations and their tempered behaviour. J. Differ. Equ. 2012, 252: 4333-4356. 10.1016/j.jde.2012.01.010

    Article  MathSciNet  MATH  Google Scholar 

  26. Lions JL: Quelques Méthodes de Résolution des Problìms aux Limites Non Liníaires. Dunod, Paris; 1969.

    Google Scholar 

Download references

Acknowledgements

This work was supported in part by the NNSF of China (No. 11271066, No. 11031003), by the grant of Shanghai Education Commission (No. 13ZZ048), by Chinese Universities Scientific Fund (No. CUSF-DH-D-2013068).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Yuming Qin.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

Haiyan carried out the main part of this manuscript. Yuming participated in the discussion and corrected the main theorems. All authors read and approved the final manuscript.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 2.0 International License (https://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Li, H., Qin, Y. Pullback attractors for three-dimensional Navier-Stokes-Voigt equations with delays. Bound Value Probl 2013, 191 (2013). https://doi.org/10.1186/1687-2770-2013-191

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1687-2770-2013-191

Keywords