Submit manuscript...
Open Access Journal of
eISSN: 2641-9335

Mathematical and Theoretical Physics

Review Article Volume 1 Issue 5

Decay of solutions for 2D navier-stokes equations posed on rectangles and on a half-strip

Padilha MV, Larkin NA

Departamento de Matematica, Universidade Estadual de Maring

Correspondence: Larkin NA, Departamento de Matemática,Universidade Estadual de Maringá, Av. Colombo 5790, Agência UEM, 87020-900, Maringá, PR, Brazil

Received: June 26, 2018 | Published: September 27, 2018

Citation: Larkin NA, Padilha MV. Decay of solutions for 2D navier-stokes equations posed on rectangles and on a half-strip. Open Acc J Math Theor Phy. 2018;1(5):203-208. DOI: 10.15406/oajmtp.2018.01.00035

Download PDF

Abstract

Initial-boundary value problems for 2D Navier-Stokes equations posed on rectangles and on a half-strip were considered. The existence and uniqueness of regular global solutions on rectangles and their exponential decay as well as exponential decay of generalized solutions on a half-strip have been established.

Keywords: navier-stokes equations, lipschitz and smooth domains, decay in bounded and unbounded domains

Introduction

The main goal of this work is establishing of sharp estimates for the exponential decay rates of solutions to initial-boundary value problems for the 2D Navier-Stokes equations:

ut+(u)u=νΔup,inΩ×(0,t),ut+(u)u=νΔup,inΩ×(0,t),  (1)

u=0inΩ,u|Ω=0,u=0inΩ,u|Ω=0,  (2)

u(x,y,0)=u0(x,y),u(x,y,0)=u0(x,y),  (3)
where ΩΩ is either a bounded rectangle or a half-strip in 2R2  with the homogeneous Dirichlet condition on the boundary of Ω.Ω.

The question of decay of the energy for generalized solutions had been stated by J.Leray1 and attracts till now attention of many pure and applied mathematicians2‒11. In all of these papers, the decay rate of u(t)L2(Ω)u(t)L2(Ω)  was controlled by the first eigenvalue of the operator A=PΔ,A=PΔ, where PP is the projection operator on solenoidal subspace of L2(Ω).L2(Ω).  Associated with stability questions, problems on dimensions of attractors and nonlinear spectral manifolds also have been studied.2,6,7

It is well-known that solutions of the 2D Navier-Stokes equations posed on smooth bounded domains with the Dirichlet boundary conditions are globally regular.9,11‒14 On the other hand, the question of regularity is not obvious in the case of bounded Lipschitz domains and unbounded Lipschitz and smooth domains. It has been proved that for Lipschitz domains, bounded and unbounded, there exists a unique global generalized solution.9,11,14

u,utL(0,;L2(Ω))L2(0,;H1(Ω)),u,utL(0,;L2(Ω))L2(0,;H1(Ω),)

but it was not clear whether

uL(0,;H2(Ω))uL(0,;H2(Ω))

at least for bounded Lipschitz domains.

In this work, we have established this fact for rectangles making use of ideas.15 The following inequality holds for rectangles

u2H2(Ω)(t)+ut2L2(Ω)(t)Cu02H2(Ω)exp(ν(π2L2+π2B2)t)u2H2(Ω)(t)+ut2L2(Ω)(t)Cu02H2(Ω)exp(ν(π2L2+π2B2)t)

and
u2H1(Ω)(t)+ut2L2(Ω)(t)Cu02H2(Ω)exp(νπ2B2t)u2H1(Ω)(t)+ut2L2(Ω)(t)Cu02H2(Ω)exp(νπ2B2t)

for a half-strip.

Our paper has the following structure: Chapter I is Introduction. Chapter 2 contains notations and auxiliary facts. In Chapter 3, existence and uniqueness of global generalized solutions on either bounded or unbounded Lipschitz domains have been established. In Chapter 4, regularity and decay of solutions on rectangles and on a half-strip have been studied.

Notations and auxiliary facts

Let ΩΩ be a domain in 2.R2.  Define as in:11

Dx=x,Dy=y,Dj=DxjxDyjy=|jx+jy|xjxyjy.Dx=x,Dy=y,Dj=DxjxDyjy=|jx+jy|xjxyjy.

We denote for scalar functions f(x,y,t)f(x,y,t)  by Lp(Ω),1<p<+  the Banach space with the norm

fpLp(Ω)=Ω|f|pdxdy,p(1,+),fL(D)=esssupD|f(x,y)|.

For p=2,L2(Ω  is a Hilbert space with the scalar product

(u,v)=Ωu(x,y)v(x,y)dxdyandthenormu2=Ω|u(x,y)|2dxdy.

The Sobolev space Wm,p(Ω)  is a Banach space with the norm
uWk,p(Ω)=0|α|kDαuLp(Ω).
When p=2,Wm,2(Ω)=Hm(Ω)  is a Hilbert space with the following scalar product and the norm:

((u,v))Hm(Ω)=|j|m(Dju,Djv),u2Hm(Ω)=|j|mDju2.

Let D(Ω) or D(ˉΩ) be the space of C  functions with compact support in Ω  or ˉΩ . The closure of C  functions in Wm,p(Ω)  is denoted by Wm,p0(Ω)  and (Hm0(Ω)whenp=2).

Define the auxiliary spaces which are projections for the solenoidal vector functions,

V={uD(Ω),u=0},V=theclosureofVinH10(Ω),

H=theclosureofVinL2(Ω),

The space H is eqquiped with the natural L2 inner product. The space V will be equipped with the scalar produt
((u,v))=(Dxu,Dxv)+(Dyu,Dyv)

when Ω is bounded. If Ω is unbounded, we define the inner product as the sum of the inner products as following:

                                                      

[[u,v]]=(u,v)+((u,v)).

We use the usual notations of Sobolev spaces Wk,p, Lp and Hk for vector functions and the following notations for the norms:
i) For vector functions u(x,y,t)=(u1(x,y,t),u2(x,y,t)),

upLp(Ω)=Ω(|u1|p+|u2|p)dxdy,

uWk,p(Ω)=0|α|kDαu1Lp(Ω)+Dαu2Lp(Ω),p(1,+).

The closures of V in L2(Ω) and in H10(Ω) are the basic spaces in our study. We denote them by H and V respectively.. Obviously V is a subspace of H10(Ω).

Define the operator

(u)u=(u1u1x+u1u2x+u2u1y+u2u2y).

Lemma 4.1 (The Steklov Inequality)16 Let vH10(0,L).  Then

π2L2v2(t)vx2(t).  (4)

Proof. Let v(t)H10(0,π) , then by the Fourier series,

π0v2t(t)dtπ0v2(t)dt.

Inequality (4) follows by a simple scaling.

Lemma 4.2 (Differential form of the Grownwall Inequality) Let I=[t0,t1]. Suppose that functions a,b:I  are integrable and a function a(t) may be of any sign. Let u:I  be a differentiable function satisfying

u(t)a(t)u(t)+b(t),fortIandu(t0)=u0,  (5)

then
u(t)u0ett0a(t)dt+tt0est0a(r)drb(s)ds.  (6)

Proof. Multiply (5) by the integrating factor est0a(r)dr and integrate from t0 to t.

The next Lemmas will be used in estimates:

Lemma 4.3 (See: 11,14 ) Let vH10(Ω) , then

vL4(Ω)21/4v1/2L2(Ω)v1/2L2(Ω).  (7)

Lemma 4.4 (See: 14 ) Let b(u,v,w)=((u)v,w),  then

|b(u,v,w)|21/2u1/2u1/2VvVw1/2w1/2V

u,v,wH10(Ω) . If uL2(0,;V)L(0,;H) , then we can define the operator Bu such that Bu belongs to L2(0,;V)  and
(Bu,v)=b(u,u,v),

BuL2(0,;V)21/2uL(0,;H)uL2(0,;V).

Existence theorems

Let Ω  be a bounded Lipschitz domain. Given u0H, consider the following problem:

(utνΔu+p+(u)u=0inΩ×(0,t),u=0inΩ×(0,t),u=0onΩ×(0,t),t>0,u(x,y,0)=u0(x,y),inΩ  (8)

equivalent to the variational problem given by,11

(u+Au+Bu=0in(0,t),t>0u(0)=u0,  (9)

where AuV such that (Au,v)=ν((u,v)) for all vV and BuV such that

(Bu,v)=b(u,u,v).  (10)

Theorem 5.1 Given u0H2(Ω)V, there exists a unique generalized solution u to (8) such that for all ΦV,Φ|Ω=0  it satisfies the following integral identity:

Ω{utΦ+ν(uxΦx+uyΦy)u(u)Φ}dxdy=0,  (11)

where

uL(0,;V),utL(0,;H)L2(0;;V).

Proof. The estimates that follow may be established on Gallerkin approximations.14,9 We estimate:

Estimate I - uL(0,;H)L2(0,;V) .

Multiply (9) by u to obtain

(ut,u)(t)+(Au,u)(t)=0.  (12)

It follows from here that

ddtu2(t)+2νu2(t)V=0.  (13)

Integrating (13) over (0,t), , we get

u2(t)+2νt0u2V(s)ds=u02,t>0.  (14)
Hence uL(0,;H)L2(0,;V) .

Estimate II - utL(0,;H)L2(0,;V) .

Derivating (9) and multiplying by ut, we get

ddtut2(t)+2νut2V(t)+2b(ut,u,ut)(t)=0.  (15)

By Lemma 4.4,

2|b(ut,u,ut)(t)|23/2ut(t)utVuV

νut2V+23νu2V(t)ut2(t)

and (15) becomes

ddtut2(t)+νut2V(t)ϕ(t)ut2(t),  (16)

where ϕ(t)=23νu2V(t). Making use of Lemma 4.2, we obtain

ut2(t)ut2(0)et0ϕ(s)ds.  (17)

To prove that ut(0)  is in H, multiply equation (9) by ut(t)  to get

ut2(t)+ν((u,ut))(t)+b(u,u,ut)(t)=0.

In particular, for t=0  we have

ut(0)2=ν(Δu0,ut(0))b(u0,u0,ut(0)),  (18)

where ut(0)=limt0ut(t), [19]. From this

ut(0)νΔu0+Bu0.  (19)

By the Hölder inequality,

|b(u,u,v)|uL4(Ω)(2uL4(Ω)+u2xL4(Ω)+u1yL4(Ω))v
CuVuH2(Ω)v,uH2(Ω),vL2(Ω).  (20)

Hence

|Bu0|Cu0Vu0H2(Ω)Cu02H2(Ω)  (21)

and by (18), ut(0)H.  This and (16) imply that

utL(0,;H)L2(0,;V).

Returning to (12), we calculate

νu2V(t)=(u,ut)(t)u(t)ut(t),  (22)

hence uL(0,;V).  This and (17) prove validity of (11) and consequently the existence part of Theorem 3.1. Uniqueness of the generalized solution, u, utL(0,;H)L2(0,;V)  has been established.9,14

Remark 5.1 Estimates u,utL(0,;HL2(0,;V)  were established first for Lipschitz domains9,14 and were valid also for unbounded domains with a natural condition lim|x|+u(x,y,t)=0.  We repeat them because we will need these estimates while establishing decay of solutions in bounded and unbounded Lipschitz domains.

Regularity and decay on rectangles and on the half-strip

Consider the Poisson problem in a bounded rectangle Ω

(Δu=f(x,y),(x,y)Ω,u|Ω=0,  (23)

Remark 6.1 It has been proved10 that for

Ωπ={x=(x1,...,xn),0<xi<π;i=1,...,n}

the following inequality holds

uW2,p(Ωπ)C(Ω)fLp(Ωπ).

It is easy to generalize this result for any rectangle in 2.  

Theorem 6.1 The problem (23) posed in rectangle Ω={(x,y)2,0<x<L;0<y<B}, where fLp(Ω),   1<p2, has a solution uW2,p(Ω).  Moreover,

uW2,p(Ω)cΩfLp(Ω).  (24)

Returning to the original problem for the Navier-Stokes equations,

(utνΔu+p+(u)u=0inΩ×(0,t)u=0inΩ×(0,t),u=0inΩ×(0,t),t>0,u(x,y,0)=u0(x,y),inΩ,  (25)

where u is a vector function from 2 into 2 and p is a real function from 2 into , and making use of Galerkin approximations, we establish the following result.

Theorem 6.2 Given u0H2(Ω)V,  the problem (25) has a unique solution (u,p)  such that

uL(0,;VH2(Ω)),utL(0,;H)L2(0,;V),
pL(0,;H).       (26)

Moreover,

ut(t)+u(t)H2(Ω)+p(t)Ce12χt,  (27)
where χ=ν(π2L2+π2B2)  and C  depends on u0H2(Ω).

Proof. Decay of L2 Norm

By definition,

u2V(t)=ux2(t)+uy2(t).

Since u|Ω=0, making use of Lemma 4.1, we get

ux2(t)π2L2u2(t),uy2(t)π2B2u2(t).
This implies u2V(t)(π2L2+π2B2)u2(t).  (28)

Returning to (12), we obtain

ddtu2(t)+2ν(π2L2+π2B2)u2(t)0.  (29)

Define χ=ν(π2L2+π2B2). Then (29) implies
u2(t)u02e2χt.  (30)

Decay of H1  Norm
Rewrite (15) in the form

ddtut2(t)+νut2V(t)ϕ(t)ut2(t)0,  (31)

where ϕ=23νu2V(t). Acting similarly to the proof of (29), we obtain

ut2V(t)(π2L2+π2B2)ut2(t).  (32)

Hence (31) reduces to the form

ddtut2(t)+(χϕ(t))ut2(t)0.  (33)

By Lemma 4.2,

ut2(t)ut(0)2et0ϕ(s)dseχt.  (34)

Since uL2(0,,V), then by (14),

t0ϕ(s)ds2νu02,t>0,

and it follows from (13) that

νu2V(t)(u,ut)(t)u(t)ut(t)

u0u0e2νu0eχte12χt.  (35)
Therefore u2V(t)1νu0u0e2νu0e32χt  (36)
and u2H10(Ω)(t)(1νu0u0e2νu0+u02)e32χt.  (37)

Decay of H2 -Norm
In order to estimate uH2(Ω)(t), we will use Theorem 6.1. First write (8) as

Δu=f=1ν(ut+p(u)u).

We estimate

|b(u,u,v)|(t)=|((u),v)(t)|c2u(t)L4(Ω)u(t)H10(Ω)v(t)L4(Ω)
Cu2(t)H10(Ω)v(t)L4(Ω)  (38)

and by (30),

(u)u(t)L4/3(Ω)Ce32χt.

Returning to (9), we obtain

AuL4/3(Ω)(t)BuL4/3(Ω)(t)+utL2(Ω)(t).  (39)

It follows by (38) and (34) that Au(t)L4/3Ce12χt.  By Theorem of de Rham,17 one can check that there exists p  such that11

p=ut+Au+Bu     (40)

and

pL4/3(Ω)(t)utL2(Ω)(t)+AuL4/3(Ω)(t)
+BuL4/3(Ω)(t)Ce12χt.  (41)

Since fL4/3(Ω), due to Theorem 6.1,

u(t)W2,43(Ω)utL2(Ω)(t)+pL4/3(Ω)(t)
+(u)uL4/3(Ω)(t)  (42)

and by (42), we get u(t)W2,4/3(Ω)Ce12χt.  By the Sobolev theorems,

uL(Ω)(t)CuW2,4/3(Ω)(t)Ce12χt.  (43)
This implies

Bu(t)Cu(t)L(Ω)u(t)H10(Ω)L2(Ω).

To prove that the norms utL2, pL2(Ω) and (u)uL2(Ω)  have exponential decay, we use the equality (10)

(u)u(t)=Bu(t),

where BuL2(Ω) such that

Bu,v=Ω(u1u1xv1+u1u1yv2+u2u2xv1+u2u2yv2)dΩ

for every vL2(Ω).  We calculate

|b(u,u,v)|(t)CuL(Ω)(t)uH10(Ω)(t)v(t).  (44)

Since the right-hand side of (44) has exponential decay for every vL2(Ω), it follows

(u)u(t)Ce54χt.  (45)

Returning to (9), we obtain the decay rate for the operator Au

Au(t)Bu(t)+ut(t).

It follows from (34) and (45) that Au(t)Ce12χt.  By (40),

pL2(t)ut(t)+Au(t)+Bu(t)Ce12χt.  (46)

Since now fL2(Ω),  substituting (34), (45), (46) into (24) and making use of Theorem 6,1, we prove
uH2(Ω)(t)Ce12χt.  It means that a unique generalized soliution is regular.

The proof of Theorem 6.2 is complete.

Existence and decay on the half-strip

Theorem 7.1 Consider the half-strip Ω={(x,y)2;0<x,0<y<B}.  Given u0H2(Ω)V,  the following problem:

(utνΔu+p+(u)u=0inΩ×(0,t),u=0inΩ×(0,t),u=0onΩ×(0,t),t>0,limx|u(x,y,t)|=0,t>0,u(x,y,0)=u0(x),inΩ  (47)

has a unique solution (u,p) such that

uL(0,;H10(Ω)),utL(0,;L2(Ω)),
pL(0,;L2(Ω)).  (48)

Moreover,

ut(t)+uH10(Ω)(t)+pL4/3(Ω)(t)C2e12θt,  (49)

where θ=νπ2B2 and C2 depends on ν, u0H2(Ω).

Proof. Obviously, the variational formulation of (47) is also (9). Repeating the proof of Theorem 5.1 (see Remark 3.1), we can proof the existence and uniqueness of the generalized solution18 to problem (47). Note that (14) holds for the problem (47). Using the Steklov inequality with respect to variable y, we obtain

uy2π2B2u2,

hence, similarly to (13),

ddtu2(t)+2νπ2B2u2(t)+2νux2(t)0.  (50)

By Lemma 5, u2(t)u02e2νπ2B2t.  (51)

Since (31) holds for the problem (47), making use of Lemma 4.4, we estimate

 

ddtut2(t)+2νut2(t)V2ut(t)u(t)Vut(t)V  (52)

which we rewrite as

ddtut2(t)+νut2V(t)2νu2V(t)ut2(t)0.  (53)

By Lemma 4.1,

uty2(t)π2B2ut2(t)

and (53) becomes

ddtut2(t)+[νπ2B22νu2V(t)]ut2(t)0.  (54)

By Lemma 5, (54) provides

ut2(t)ut2(0)e23νt0u2V(s)dseνπ2B2t,

hence

ut2(t)ut2(0)e2νu0eνπ2B2t.  (55)

Returning to (35), we estimate

u2V(t)1νut(t)u(t)1νu0ut(0)e2νu0eν3π22B2t.  (56)

Decay for Pressure
In order to obtain decay for pL4/3(Ω)(t),  we start with

(u)uL4/3(Ω)(t)=BuL4(Ω)(t),

where L4(Ω) is the dual of the space L4(Ω).  Since

Au=utBu,

repeating calculations of (38) and making use of (34), we get

Au(t)L4/3(Ω)c1e12θt.  Observing that (40) holds for the problem (47), we obtain

pL4/3(Ω)(t)utL2(Ω)(t)+AuL4/3(Ω)(t)
+BuL4/3(Ω)(t)c2e12θt.  (57)

Jointly (55), (56) and (57) prove (48), (49).

Conclusion

In our work, we tried to respond some questions posed by J. Leray,1 namely, regularity of global solutions of the Navier-Stokes equations and their decay. Therefore, our results can be divided in two parts: the first one concerns decay of global regular solutions of the 2D Navier-Stokes equations posed on rectangles.19 It is known that there exist global regular solutions for the 2D Navier-Stokes equations posed on smooth bounded domains,4,10,11,14 but regularity in nonsmooth (Lipschitz) domains, such as rectangles, is not obvious. For bounded rectangles, we have established the existence of an unique global regular solution which decays exponentially as  We demonstrated that the decay rate is different for different norms, see (26), (30), (36), whereis defined by the geometrical characteristics of a domain

The second part of our work concerns decay of solutions for the 2D Navier-Stokes equations posed on a half-strip. In existing publications,3–11 the decay rate of  is controlled by the first eigenvalue of the operatorwhereis the projection operator on solenoidal subspace of  It is clear that this approach does not work in unbounded domains On the other hand, our approach based on the Steklov inequality with respect toallowed us to estimate the decay rate of a generalized solution for the 2D Navier-Stokes equations posed on a half-strip.

We must emphasize that this estimate is the first in the history which gives an explicit value of the decay rate for unbounded domains. Results established in our work can be used in constructing of numerical schemes for solving initial-boundary value problems for the Navier-Stokes equations appearing in Mechanics of viscous liquid. From the physical point of view, decay estimates show that the decay rate of perturbations of solutions caused by the initial data is bigger for bigger values of viscosityand smaller values of the width and length of the rectangles and the width of a half-strip.

Acknowledgements

None.

Conflict of interest

The authors declare that there are no conflict of interest regarding the publication of this paper.

References

  1. Leray Essai J. Sur le mouvement d'un uide visqueux emplissant l'espace. Acta Math. 1934;63:193–248.
  2.  Brown RM, Perry P, Shen Z. On the dimension of the attractor for the non-homogeneous Navier-Stokes equations in non-smooth domains. Indiana Univ Math Journal. 2000;40(1):81–112.
  3. Edmunds DE. Asymptotic behavior of solutions of the Navier-Stokes equations. Arch Rational Mech Anal. 1966;22(1):15–21.
  4. Foias C, Prodi G. Sur le comportement global des solutions non-stationairos des equations de Navier-Stokes en dimension 2. Rendiconti del Seminario Matemtico della Universit di Padova. 1967;39:1–34.
  5. Foias C, Saut JC. Asynptotic behavior as t! 1of solutions of Navier-stokes equations in nonlinear partial dierential equations and their applications. College of France Seminar. 1983;4:74–86.
  6. Foias C, Temam R. Some analytic and geometric properties of the solutions of the evolution of Navier-Stokes equations. J Math Pures Appl. 1979;58:339–368
  7. Foias C, Saut JC. Asymptotic Behavior, as t! 1 of Solutions of Navier- Stokes Equations and Nonlinear Spectral Manifolds Indiana University Mathematics Journal. 1984;33.10.1512/iumj.1984.33.33025.
  8. Guillop C. Comportement l'in ni des solutions des equations de Navier-Stokes et proprit des ensembles fonctionnels invariants (ou attracteurs). Annales de l'institut Fourier. 1982;32(3):1–37.
  9. Ladyzhenskaya O. Solution in the large of the nonstationary boundary value problem for the navier-stokes system with two space variables. Comm Pur Applied Math. 1959;12(3):427–433.
  10. Prodi G, Qualche risultato riguardo alle equazioni di Navier-Stokes nel caso bidimensional. Rend Sem Mat Univ Padova. 1965;30:1–15.
  11. Temam R. Navier-Stokes Equations. Theory and Numerical Analysis. Noth-Holland, Amsterdam; 1979.
  12. Cattabriga L. Su un problema al contorno relativo al sistema di equazioni di Stokes. Rendiconti del Seminario Matematico della Universit di Padova. 1961;31:308–340.
  13. Constantin P, Foias C. Navier-Stokes Equations. Chicago Lectures in Mathematics; 1988.
  14. Ladyzhenskaya OA. The Mathematical Theory of Viscous Incompressible Flow. 2nd ed. English translation, New York: Gordon and Breach; 1969.
  15. Koshelev AI. A priori estimates in Lp and generalized solutions of elliptic equations and systems. Amer Math Soc Transl. 1962;20(2):105–171.
  16. Steklov AV. The problem of cooling of an heterogeneous rigid rod. Communs Kharkov Math Soc Ser. 1896;2(5):136–181.
  17. de Rham G. Vari et es dierentiables. Herman, Paris; 1960.
  18. Marcinkiewicz J. Sur les multiplicateurs des s eries de Fourier. Studia Mathematica. 1939;8(1):78–91.
  19.  Temam R. Behavior at t=0 of the solutions of the semi-linear evoluation equations. J Di Equats. 1982;43(1):73–92.
Creative Commons Attribution License

©2018 Larkin, et al. This is an open access article distributed under the terms of the, which permits unrestricted use, distribution, and build upon your work non-commercially.