DDNS Discrete Dynamics in Nature and Society 1607-887X 1026-0226 Hindawi Publishing Corporation 10.1155/2015/721673 721673 Research Article Approximate Kelvin-Voigt Fluid Driven by an External Force Depending on Velocity with Distributed Delay http://orcid.org/0000-0002-1123-1563 Guo Yantao 1,2 Cheng Shuilin 3 http://orcid.org/0000-0001-7445-2817 Tang Yanbin 1 Gori Luca 1 School of Mathematics and Statistics Huazhong University of Science and Technology Wuhan, Hubei 430074 China hust.edu.cn 2 School of Mathematics and Statistics Xuchang University Xuchang, Henan 461000 China xcu.edu.cn 3 School of Statistics and Mathematics Zhongnan University of Economics and Law Wuhan, Hubei 430073 China znufe.edu.cn 2015 11102015 2015 23 03 2015 02 05 2015 11102015 2015 Copyright © 2015 Yantao Guo et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

We consider the approximate 3D Kelvin-Voigt fluid driven by an external force depending on velocity with distributed delay. We investigate the long time behavior of solutions to Navier-Stokes-Voigt equation with a distributed delay external force depending on the velocity of fluid on a bounded domain. By a prior estimate and a contractive function, we give a sufficient condition for the existence of pullback attractor of NSV equation.

1. Introduction

In this paper, we consider 3D Navier-Stokes-Voigt (NSV) equation with a distributed delay external force depending on the velocity of the fluid:(1)ut-α2Δut-νΔu+u·u+p=-h0Gs,ut+sds,divut,x=0,t,xτ,+×Ω,ut,x=0,t,xτ,+×Ω,ut,x=ϕt-τ,x,tτ-h,τ,xΩ,where u=(u1,u2,u3) is the velocity field of the fluid, p is the pressure, ν>0 is the kinematic viscosity, α>0 is the length scale parameter of the elasticity of the fluid, the external force G and initial velocity field ϕ are defined in the interval of time [-h,0], where h is a fixed positive number and Ω is a bounded smooth domain of R3.

The NSV equation was introduced by Oskolkov  to give an approximate description of the Kelvin-Voigt fluid and was proposed as a regularization of 3D Navier-Stokes equation for the purpose of direct numerical simulations in . Since the term -α2Δut changes the parabolic character of the equation, the NSV equation being well posed in 3D, many authors have studied the long time dynamics of this model. Kalantarov and Titi  investigated the existence of the global attractor, the estimation of the upper bounds for the number of determining modes, and the dimension of global attractor of the semigroup generated by the equations. By a useful decomposition method, Yue and Zhong  proved the asymptotic regularity of solution of NSV equation and obtained the existence of the uniform attractor; they also described the structure of the uniform attractor and its regularity. García-Luengo et al.  investigated the existence and relationship between minimal pullback attractor for the universe of fixed bounded sets and universe given by a tempered condition.

Partial differential equations with delays arise from various fields, like physics, control theory, and so on (see, e.g., ); the unknown functions depend on not only present stage but also some past stage. The existence and stability of solution and global attractor for Navier-Stokes equation with discrete delay were established in . The existence of pullback attractors in CH01 and CH01H2 was proved for the processes associated with nonclassical diffusion equations with variable bounded delay in [14, 15]. Delay effect has been considered on an unbounded domain in . The existence of pullback attractor for a Navier-Stokes equation with infinite discrete delay effect was studied in .

The aim of this paper is to investigate the NSV equation with a distributed delay, instead of the discussions with finite delays in the references. Our purpose is twofold. We first show the existence and uniqueness of solution to NSV equation (1) with a distributed delay; then we prove the existence of pullback attractor for the process generated by the NSV equation (1).

This paper is organized as follows. In Section 2, we give some preliminary results and prove existence of solution to NSV equation with a distributed delay. In Section 3, we derive the existence of pullback attractor by prior estimates and contractive functions.

2. Existence of Solutions

In order to prove the existence of solutions to problem (1), we define the function spaces (2)V=uC0Ω3,divu=0.H is the closure of V in (L2(Ω))3 with the inner product (·,·) and associate norm |·|, V is the closure of V in (H01(Ω))3 with scalar product ((·,·)) and associate norm ·, where (3)u,v=i,j=13Ωujxivjxidx,u,vH01Ω3;it follows that VHHV, where the injections are dense and compact. We will use · for the norm in V and ·,· for the duality pairing between V and V.

Define the linear continuous operator A:VV as(4)Au,v=u,v,u,vV.We denote D(A)={uV,AuH}; one has that D(A)=(H2(Ω))3V and Au=-PΔu, for all uD(A) is the Stokes operator, where P is the orthoprojector from (L2(Ω))3 onto H; also denote CH=C0([-h,0];H) and CV=C0([-h,0];V).

Define the trilinear form b on V×V×V by (5)bu,v,w=i,j=13Ωuivjxiwjdx,u,v,wVand the operator B:V×VV as (6)Bu,v,w=bu,v,w,u,v,wV,and denote B(u)=B(u,u).

The trilinear form b satisfies that(7)bu,v,w=-bu,w,v,bu,v,v=0,u,v,wV.We also recall that there exists a constant C depending only on Ω such that(8)bu,v,wCu1/2Au1/2vw,uDA,  vV,  wH,bu,v,wCuvw1/2w1/2,uV,  vV,  wV,bu,v,wCu1/2u1/2vw,uV,  vV,  wV.

For the term containing the time delay, G:R×HH satisfies that

(H1)G(·,u):RH is a measurable function,

(H2)G(t,0)=0 for all tR,

(H3) there exists a positive constant L, such that R>0; if |u|<R and |v|<R, then(9)Gt,u-Gt,v2Lu-v2.

Remark 1.

Hypotheses H2-(H3) imply that |G(t,u)|2L|u|2, so we have |G(t,u)|2L(τ,T) for |u|<R.

Problem (1) can be rewritten as(10)tu+α2Au+νAu+Bu,u=-h0Gs,ut+sds,ut,x=ϕt-τ,x,tτ-h,τ,xΩ;then we get the existence of solution to problem (10).

Theorem 2.

Let ϕCV, let G:R×HH satisfy the hypotheses (H1)(H3), and let τR. Then, T>τ; there exists a unique weak solution to (10) such that(11)ut,xCτ-h,T;V,ut,xtL2τ,T;V.Moreover, if ϕCD(A), then problem (10) admits a strong solution.

Proof.

Consider the Galerkin approximations for problem (10):(12)ddtum+α2Aum+νAum+Bum,um=-h0Gs,umt+sds,where um=j=1mujmej, Aum=j=1mλjujmej, and ej and λj are the corresponding orthonormal eigenfunctions and eigenvalues of operator A, respectively; then, (13)um2=j=1mλjujm2,Aum2=j=1mλj2ujm2,um2=j=1mujm2.

We now derive a prior estimate for the Galerkin approximate solution. Multiplying (12) by ujm, summing from j=1 to m and using the fact (14)Bum,um,um=bum,um,um=0,we obtain that, for a.e. t>τ,(15)ddtum2+α2um2+2νum2=2-h0Gs,umt+sds,um-h0Gs,umt+sds2+um2.Integrating (15) from τ to t, we obtain that (16)umt2+α2umt2+2ντtums2dsumτ2+α2umτ2+τt-h0Gs,umr+sds2dr+τtumr2dr.Remark 1 implies that(17)τt-h0Gs,umr+sds2dr-h0τtGs,umr+s2drds-h0τtLumr+s2drdsL-h0τ+st+sumr2drdsL-h0τ-hτumr2dr+τtumr2drdsLh2ϕCH2+Lhτtumr2dr.Then, t(τ,T) and(18)umt2+α2umt2+2ντtums2dsLh+2τtums2ds+Lh2ϕCH2+ϕCH2+α2ϕCV2.So, we have (19)umt2Cτtums2ds+CϕCH2+ϕCH2+α2ϕCV2.The Gronwall inequality implies that(20)umt2C.Putting (20) into the right-hand side of (18), we have (21)α2umt2+2ντtums2dsC,tτ,T.This implies that(22)um is bounded in  Lτ,T;VL2τ,T;V.

Now, multiplying (12) by tum and integrating over Ω, we have (23)tum2+α2tum2+ν2ddtum2bum,um,tum+-h0Gs,umt+sds,tumcum3/2um1/2tum+12-h0Gs,umt+sds2+12tum2;since(24)um3/2um1/2tumcum3um+α22tum2;then(25)tum2+α2tum2+νddtum2cum3um+-h0Gs,umt+sds2,integrating the above inequality from τ to t, by (17), (20), and (22) we have(26)τttum2+α2tum2ds+νumt2νumτ2+cτtum3umds+τt-h0Gs,umr+sds2drνϕτ2+cτtum3umds+cϕCH2+cτtumr2drνϕCV2+cτtum3umds+cϕCH2+cτtumr2dr.Since {um} is bounded in L(τ,T;V)L2(τ,T;V), we obtain that(27)tum is bounded in  L2τ,T;V.

By the Faedo-Galerkin scheme, for example, see [14, 18], according to the estimates (22) and (27), we can get existence of the weak solution; here we omit the details.

We next consider the uniqueness of solution. Let u,v be two solutions to problem (10) corresponding the initial data ϕ and ψ, respectively.

Denote w=u-v; then, we have(28)tw+α2Aw+νAw+Bu,u-Bv,v=-h0Gs,ut+sds--h0Gs,vt+sds.Multiplying (28) by w and integrating over Ω, we obtain(29)12ddtw2+α2w2+νw2+Bu,u,w-Bv,v,w=-h0Gs,ut+sds--h0Gs,vt+sds,w.Notice that(30)Bu,u,w-Bv,v,w=bu,u,u-v-bv,u,u-v+bv,u,u-v-bv,v,u-v=bu-v,u,u-v-bv,u-v,u-v=bw,u,w.

Substituting (30) into (29) and integrating from τ to t, we get(31)wt2+α2wt2+2ντtws2ds-wτ2-α2wτ2τtbw,u,wds+τtws2ds+τt-h0Gs,ur+sds--h0Gs,vr+sds2dr.(H3) implies that(32)τt-h0Gs,ur+sds--h0Gs,vr+sds2dr-h0τtGs,ur+s-Gs,vr+s2drds-h0τtGs,ur+s-Gs,vr+s2drdsL-h0τtur+s-vr+s2drdsL-h0τ+st+sur-vr2drdsL-h0τ-hτur-vr2dr+τtur-vr2drdsLh2ϕ-ψCH2+Lhτtws2ds.As the property of operator b and Poincaré, we have(33)τtbw,u,wdsCτtw3/2w1/2udsCτtw2uds.Substituting (32) and (33) into (31), we get (34)wt2+α2wt2+2ντtws2dsτtw2uds+Lh2ϕ-ψCH2+Lh+1τtws2ds+wτ2+α2wτ2τtw2uds+Lh2ϕ-ψCH2+Lh+1τtws2ds+ϕ-ψCH2+α2ϕ-ψCV2Cτtws2ds+Cϕ-ψCV2.The last inequality comes from Poincaré inequality and the boundedness of u. Therefore, the Gronwall inequality implies the uniqueness of the solution. The proof is complete.

3. Existence of Pullback Attractor

In this section, we will prove the existence of pullback attractor to problem (10). First we give existence of pullback absorbing set for the process {U(t,τ)} generated by the global solution to problem (10).

Lemma 3.

Assume (H1)(H3) hold and Lν2λ12ɛ/4eɛh; then, the process {U(t,τ)} is pullback dissipative, where 0<ɛ<min{ν/α2,νλ1/4}.

Proof.

Multiplying (10) by u and integrating over Ω, we obtain (35)tu2+α2u2+2νu2ηu2+1η-h0Gs,ut+sds2,where η is a constant determined later.

By Poincaré inequality, we have (36)tu2+α2u2+νu2+νλ1-ηu21η-h0Gs,ut+sds2.Since(37)e-ɛtteɛtu2+α2u2=ɛu2+α2u2+tu2+α2u2,then(38)e-ɛtteɛtu2+α2u2ɛ+η-νλ1u2+ɛα2-νu2+1η-h0Gs,ut+sds2.Integrating (38) from τ to t, we get(39)eɛtut2+α2ut2-eɛτuτ2+α2uτ2τteɛsɛ+η-νλ1u2+ɛα2-νu2ds+1ητteɛr-h0Gs,ur+sds2dr.Assumptions (H1)(H3) imply that(40)τteɛr-h0Gs,ur+sds2drτteɛr-h0Gs,ur+s2dsdrL-h0τteɛrur+s2drdsL-h0τ+st+seɛr-sur2drdsL-h0τ-hteɛr-sur2drdsL-h0e-ɛsτ-hteɛrur2drdsL-h0e-ɛsτ-hτeɛrur2dr+τteɛrur2drdsLeɛhɛ1ɛϕCH2eɛτ+τteɛrur2drLeɛτ+hɛ2ϕCH2+Leɛhɛ·τteɛrur2dr.Substituting (40) into (39), we have (41)eɛtut2+α2ut2-eɛτuτ2+α2uτ2τteɛsɛ+η-νλ1+Leɛh2ηɛu2+ɛα2-νu2ds+Leɛτ+h2ηɛ2ϕCH2.Let η=(1/2)νλ1, choosing 0<ɛ<min{ν/α2,νλ1/4}. Lν2λ12ɛ/4eɛh implies that (42)maxɛ+η-νλ1+η-1Leɛh2ɛ,ɛα2-ν<0;then,(43)eɛtut2+α2ut2Leɛτ+hɛ2νλ1ϕCH2+eɛτuτ2+α2uτ2,which implies (44)ut2+α2ut2Leɛτ+h-tɛ2νλ1ϕCH2+eɛτ-tϕτ2+α2ϕτ2.Now, if we take tτ+h, then, for θ[-h,0], we have (45)ut+θ2+α2ut+θ2Leɛτ+h-t-θɛ2νλ1ϕCH2+eɛτ-t-θϕCH2+α2ϕCV2.

We denote by R the set of all functions r:(-,+)(0,+) such that (46)limτ-eατr2τ=0.Then, the closed ball in CV defined by(47)B=φCV:φCV21is pullback absorbing set for {U(t,τ)}. The proof is complete.

We next prove the asymptotic compactness of solution to problem (10) by contractive functions; see [19, 20].

Let X be a Banach space and let B be a bounded subset of X. We call a function Φ(·,·) which, defined on X×X, is a contractive function on B×B if for any sequence {xn}n=1B there is a subsequence {xnk}k=1{xn}n=1 such that (48)limllimkΦxnk,xnl=0.Denote all such contractive functions on B×B by C(B).

Theorem 4 (see [<xref ref-type="bibr" rid="B9">19</xref>]).

Let {S(t)}tτ be a semigroup on a Banach space (X,·) and have a bounded absorbing set B, ϵ>0; there exist T=T(B,ϵ) and ΦC(B), such that (49)STx-STyϵ+ΦTx,y,x,yB,where ΦT depends on T. Then, {S(t)}tτ is asymptotically compact in X.

Lemma 5.

Assume that (H1)(H3) hold; the process {U(t,τ)}tτ generated by the global solution to problem (10) is asymptotically compact.

Proof.

Let ui(t) be the solution to problem (10) with initial data ϕiB (i=1,2), respectively. Denote v(t)=u1(t)-u2(t); then, v(t) satisfies the equivalent abstract equation(50)tv+α2Av+νAv+Bu1,u1-Bu2,u2=-h0Gs,u1t+s-Gs,u2t+sdswith the initial condition v(t)=ϕ1(t-τ)-ϕ2(t-τ), t[τ-h,τ].

Set an energy function(51)Evt=12Ωv2dx+α22Ωv2dx.Multiplying (50) by v and integrating over [s,T]×Ω with T>t+τ, sτ, we have(52)EvT-Evs+νsTΩv2dxdr+sTΩBu1-Bu2vrdxdr=sTΩ-h0Gσ,u1r+σ-Gσ,u2r+σvrdσdxdr;then we have(53)νsTΩv2dxdrEvs-sTΩBu1-Bu2vrdxdr+sTΩ-h0Gσ,u1r+σ-Gσ,u2r+σvrdσdxdr.Using Poincaré inequality and (51) and (53), we have(54)τTEvsds=12τTΩv2dxdr+α22τTΩv2dxdrCτTΩv2dxdrCEvτ-CτTΩBu1-Bu2vrdxdr+CτTΩ-h0Gσ,u1r+σ-Gσ,u2r+σvrdσdxdr.Integrating (52) from τ to T with respect to s, we obtain(55)TEvT+ντTsTΩv2dxdrds=τTEvsds-τTsTΩBu1-Bu2vrdxdrds+τTsTΩ-h0Gσ,u1r+σ-Gσ,u2r+σvrdσdxdrds.Substituting (54) into (55), we get(56)TEvT+ντTsTΩv2dxdrdsCEvτ-CτTΩBu1-Bu2vrdxdr+CτTΩ-h0Gσ,u1r+σ-Gσ,u2r+σvrdσdxdr-τTsTΩBu1-Bu2vrdxdrds+τTsTΩ-h0Gσ,u1r+σ-Gσ,u2r+σvrdσdxdrds.Set(57)C0=CEvτ=C2Ωϕ10-ϕ202dx+Cα22Ωϕ10-ϕ202dx,(58)Φu1,u2=-CτTΩBu1-Bu2vrdxdr+CτTΩ-h0Gσ,u1r+σ-Gσ,u2r+σvrdσdxdr-τTsTΩBu1-Bu2vrdxdrds+τTsTΩ-h0Gσ,u1r+σ-Gσ,u2r+σvrdσdxdrds;then, by (56), we have(59)EvTC0T+1TΦu1,u2.One has ϵ>0; we choose T>0 large enough such that C0/T<ϵ. By Theorem 4, it suffices to prove that Φ(u1,u2) is a contractive function, since B is a bounded positive invariant set.

If umu (m), we have the limits(60)limnlimmτTΩBun-Bumun-umdxdr=limnlimmτTΩun-um·un-um·un-umun-umdxdr=0,limnlimmτTsTΩBun-Bumun-umdxdrds=limnlimmτTsTΩun-um·un-um·un-umun-umdxdrds=0.By (H1)(H3), we have(61)limnlimmτTΩ-h0Gσ,unr+σ-Gσ,umr+σvrdσdxdrlimnlimmτTΩ-h0Lunr+σ-umr+σunr-umrdσdxdr=0,limnlimmτTsTΩ-h0Gσ,unr+σ-Gσ,umr+σvrdσdxdrdslimnlimmτTsTΩ-h0Lun-umr+σun-umrdσdxdrds=0.Combining (60)-(61) with (58), we have that Φ(u1,u2) is a contractive function.

The proof is complete.

Theorem 6.

Assume that (H1)(H3) hold and Lν2λ12ɛ/4eɛh; then, the process {U(t,τ)} generated by the global solution for problem (10) has a pullback attractor.

Proof.

Lemma 3 implies {U(t,τ)} has a pullback absorbing set B and Lemma 5 implies {U(t,τ)} is asymptotically compact; we obtain the conclusion immediately.

According to Theorem 6, 0<ɛ<min{ν/α2,νλ1/4} and Lν2λ12ɛ/4eɛh. Under the assumptions (H1)(H3), the NSV equation (1) with a distributed delay has a pullback attractor.

Conflict of Interests

The authors declare that they have no competing interests.

Authors’ Contribution

Yantao Guo carried out the long time behavior of solutions. Shuilin Cheng carried out the pullback attractor. Yanbin Tang carried out the distributed delay. All authors read and approved the final paper.

Acknowledgment

This work was supported by the National Natural Science Foundation of China (Grant nos. 11471129 and 11272277).

Oskolkov A. P. The uniqueness and solvability in the large of boundary value problems for the equations of motion of aqueous solutions of polymers Zapiski Nauchnykh Seminarov Leningrad Otdel Mathematics Institute Stekov (LOMI) 1973 38 98 136 Cao Y. Lunasin E. M. Titi E. S. Global well-posedness of the three-dimensional viscous and inviscid simplified Bardina turbulence models Communications in Mathematical Sciences 2006 4 4 823 848 10.4310/CMS.2006.v4.n4.a8 MR2264822 ZBL1127.35034 Kalantarov V. K. Titi E. S. Global attractors and determining modes for the 3D Navier-Stokes-Voight equations Chinese Annals of Mathematics, Series B 2009 30 6 697 714 10.1007/s11401-009-0205-3 MR2570790 2-s2.0-71249118247 Yue G.-C. Zhong C.-K. Attractors for autonomous and nonautonomous 3D Navier-Stokes-Voight equations Discrete and Continuous Dynamical Systems Series B 2011 16 3 985 1002 10.3934/dcdsb.2011.16.985 MR2806333 2-s2.0-80053178469 García-Luengo J. Marín-Rubio P. Real J. Pullback attractors for three-dimensional non-autonomous Navier-Stokes-Voigt equations Nonlinearity 2012 25 4 905 930 10.1088/0951-7715/25/4/905 MR2891552 2-s2.0-84859080811 Hale J. Theory of Functional Differential Equations 1977 New York, NY, USA Springer MR0508721 Tang Y. Wang M. A remark on exponential stability of time-delayed Burgers equation Discrete and Continuous Dynamical Systems. Series B 2009 12 1 219 225 10.3934/dcdsb.2009.12.219 MR2505671 2-s2.0-69249188563 Tang Y. Zhou L. Stability switch and Hopf bifurcation for a diffusive prey-predator system with delay Journal of Mathematical Analysis and Applications 2007 334 2 1290 1307 10.1016/j.jmaa.2007.01.041 MR2338664 2-s2.0-34250634518 Wu J. H. Theory and Applications of Partial Differential Equations 1996 New York, NY, USA Springer Zhou L. Tang Y. B. Hussein S. Stability and Hopf bifurcation for a delay competition diffusion system Chaos, Solitons & Fractals 2002 14 8 1201 1225 10.1016/s0960-0779(02)00068-1 MR1924523 2-s2.0-0036833143 Caraballo T. Real J. Navier-Stokes equations with delays The Royal Society of London. Proceedings. Series A. Mathematical, Physical and Engineering Sciences 2001 457 2014 2441 2453 10.1098/rspa.2001.0807 MR1862662 2-s2.0-34548675437 Caraballo T. Real J. Asymptotic behaviour of two-dimensional Navier-Stokes equations with delays The Royal Society of London. Proceedings. Series A. Mathematical, Physical and Engineering Sciences 2003 459 2040 3181 3194 10.1098/rspa.2003.1166 MR2027360 2-s2.0-1542512445 Caraballo T. Real J. Attractors for 2D-Navier-Stokes models with delays Journal of Differential Equations 2004 205 2 271 297 10.1016/j.jde.2004.04.012 MR2091818 ZBL1068.35088 2-s2.0-4344688565 Hu Z. Wang Y. Pullback attractors for a nonautonomous nonclassical diffusion equation with variable delay Journal of Mathematical Physics 2012 53 7 17 072702 10.1063/1.4736847 MR2985228 2-s2.0-84864771865 Li H. Y. Qin Y. M. Pullback attractors for three-dimensional navier-stokes-voigt equations with delays Boundary Value Problems 2013 2013, article 191 10.1186/1687-2770-2013-191 2-s2.0-84887914763 Marín-Rubio P. Real J. Attractors for 2D-Navier-Stokes equations with delays on some unbounded domains Nonlinear Analysis. Theory, Methods & Applications 2007 67 10 2784 2799 10.1016/j.na.2006.09.035 MR2347436 2-s2.0-34547138177 Marín-Rubio P. Real J. Valero J. Pullback attractors for a two-dimensional Navier-Stokes model in an infinite delay case Nonlinear Analysis: Theory, Methods & Applications 2011 74 5 2012 2030 10.1016/j.na.2010.11.008 MR2764398 2-s2.0-78651378960 Marín-Rubio P. Real J. Pullback attractors for 2D-Navier-Stokes equations with delays in continuous and sub-linear operators Discrete and Continuous Dynamical Systems. Series A 2010 26 3 989 1006 10.3934/dcds.2010.26.989 MR2600726 2-s2.0-77950282483 Khanmamedov A. K. Global attractors for wave equations with nonlinear interior damping and critical exponents Journal of Differential Equations 2006 230 2 702 719 10.1016/j.jde.2006.06.001 MR2269940 ZBL1110.35016 2-s2.0-33748602633 Yang L. Zhong C.-K. Global attractor for plate equation with nonlinear damping Nonlinear Analysis: Theory, Methods & Applications 2008 69 11 3802 3810 10.1016/j.na.2007.10.016 MR2463335 2-s2.0-53949105051