FROM SECOND GRADE FUIDS TO TE NAVIER-STOKES EQUATIONS A. V. BUSUIOC arxiv:607.06689v [math.ap] Jul 06 Abstract. We consider the limit α 0 for a second grade fluid on a bounded domain with Dirichlet boundary conditions. We show convergence towards a solution of the Navier-Stokes equations under two different types of hypothesis on the initial velocity u 0. Iftheproduct u 0 u 0 issufficientlysmallweproveglobal-in-timeconvergence. If there is no smallness assumption we obtain local-in-time convergence up to the time C/ u 0 4.. Introduction We consider in this paper the incompressible second grade fluid equations: () t (u α u) ν u+u (u α u)+ (u α u) u = p, divu = 0, where α and ν are some non-negative constants. The fluid is assumed to be enclosed in a bounded smooth region of R 3 and the homogeneous Dirichlet boundary conditions are imposed () u(t, ) = 0 for all t 0. The initial value problem is considered and we denote by u 0 the initial velocity: (3) u(0,x) = u 0 (x). The equations () were deduced in [8] from physical principles. et us ust mention here that the second grade fluids are characterized by the following fact: the stress tensor is a polynomial of degree two in the first two Rivlin-Ericksen tensors which are the deformation tensor D and the tensor ( t +u )D. The vanishing viscosity case ν = 0 is also known under the name α Euler or Euler α equations and was later obtained via an averaging procedure performed on the classical incompressible Euler equations. Two main boundary conditions were used for () in the mathematical literature: the no-slip boundary conditions and the frictionless slip Navier boundary conditions where the fluid is allowed to slip on the boundary without friction. The second boundary condition is more complex but allows for better mathematical results; also it has less physical relevance. The classical well-posedness results for () are the following (see [6, 7,, ] for the Dirichlet boundary conditions and [5] for the Navier boundary conditions): In dimension two there exists a unique global 3 solution if u 0 3. In dimension three there exists a unique local 3 solution if u 0 3. The solution is global if u 0 is small in 3. We call 3 solution a divergence free vector field verifying the boundary conditions and the PDE () and who is bounded in time (up to time t = 0) with values in 3 (). et us also mention the paper [] where solutions in W,p, p > 3, are constructed. et us observe that when α = 0 relation () becomes the Navier-Stokes equations (4) t u ν u+u u = p
and when α = ν = 0 it becomes the Euler equations. It is interesting to know if the solutions of () converge to the solutions of the limit equation when α 0 and ν > 0 is fixed or when α,ν 0. This was already studied in several papers as we shall see below. et us first mention that in the absence of boundaries one can obtain 3 estimates uniform in α and ν in both dimensions two and three and pass to the limit. This was performed in [4], see also [3] for a simpler proof. But such a result cannot hold true on domains with boundaries. Indeed, if the solutions of ()-() are bounded in 3 uniformly in α then one can easily pass to the limit α 0 and obtain at the limit a solution of the Navier-Stokes equations which must also be bounded in 3. For such a solution to exist, the initial data must verify a compatibility condition. Indeed, one can apply the eray proector to (4) to obtain that t u νp u + P(u u) = 0. Since u vanishes at the boundary, so does t u. We infer that νp u + P(u u) = 0 at the boundary. Observe that u being in 3 implies that these two terms are in so the trace at the boundary makes sense. By time continuity we infer that the initial data must verify the compatibility condition νp u 0 = P(u 0 u 0 ) at the boundary. This is of course in general not verified if we only assume that u 0 3 is divergence free and vanishing on the boundary. We review now the results available for domains with boundary. Concerning the limit α,ν 0, we proved in [3] the expected convergence in D for weak solutions in the case of the Navier boundary conditions. We also proved the convergence in 3D but under the additional hypothesis that the solutions exist on a time interval independent of α and ν. We proved in [4] that the hypothesis of existence of a uniform time existence is verified if ν = 0 and α 0. In the case of the Dirichlet boundary conditions, there is only the paper [5] which shows convergence in D. The case of Dirichlet boundary conditions in 3D is open. Concerning the limit α 0 and ν > 0 fixed, it was proved in [3] the expected convergence for weak solutions. That result is stated for = R n, n =,3, but the proof relies only on energy estimates and standard compactness arguments so it goes through to bounded domains without difficulty (Dirichlet and Navier boundary conditions likewise), see also [3, Remarque 4]. There is however a maor drawback to the result of [3]: the author assumes that the sequence of weak solutions exist on a time interval independent of α and proves convergence on any such time interval. This raises the question of proving the existence of such a uniform time interval. et us also mention the paper [] where the author considers the D case with Navier boundary conditions and shows stronger convergence of solutions together with some estimates for the rate of convergence. Observe however that in dimension two the solutions are global in time so the hypothesis assumed in [3] that the solutions exist on a uniform time interval is automatically satisfied. In this paper we aim to prove that the hypothesis of [3] about the existence of solutions on a uniform time interval is verified in various situations for Dirichlet boundary conditions. Since global existence of solutions holds true in dimension two for Dirichlet boundary conditions, in D the hypothesis of [3] is automatically satisfied and the problem is settled. So we restrict ourselves to the 3D case. Recall that the Navier-Stokes equations in dimension three are locally well-posed for large data and globally well-posed for small data. ikewise, we will prove two results on the uniform time of existence: a local result for large data and a global result for small data. Surprisingly, we find that if the norm of the initial data is sufficiently small and if α is sufficiently small too, then the solutions of () are global. et us emphasize that the smallness of the initial data is measured only in the norm (in fact in a weaker space, see Theorem in the next section) and not in the 3 norm as required by the classical global well-posedness result for (). It was expected to find that, if the norm
of the initial data is small then the maximal time of existence of () goes to infinity as α 0. Indeed the limit equations, i.e. the Navier-Stokes equations, are globally wellposed for small initial data and so we expect convergence for the times of existence of solutions. But we did not expect the maximal time of existence to actually be infinite if α is sufficiently small. For large data, we prove in particular that if α is sufficiently small then the time existence of the solution has a lower bound that depends only on the norm of the initial velocity. This should be compared with previous results on local existence of solutions where the time of existence depends on the 3 norm or the W,p, p > 3, norm of the initial velocity. When α 0 we obtain convergence of solutions of () towards solutions of (4) up to the time C/ u 0 4. The plan of the paper is the following. In the next section we state and comment our results. In Section 3 we prove the global result for small data. In the last section, we prove our uniform local existence result.. Statement of the results Since we are interested in the limit α 0 with ν fixed, we will assume throughout this paper that α and ν > 0. et us first state the result of [3] which motivates the present work. Theorem (see [3]). et ν > 0, T > 0 and u α some solutions of () defined up to the time T such that u α (0) Cα and u α (0) converges weakly in to some u 0. Then there exists a weak eray solution u of the Navier-Stokes equations with initial velocity u 0 such that, after extraction of a sub-sequence, u α u in loc ([0,T); ) weak* and in loc ([0,T); ) weak. Our result about global existence of solutions reads as follows. Theorem. There exists some small constant ε = ε() depending only on such that if the initial velocity u 0 belongs to 3 (), is divergence free, vanishes at the boundary and verifies the following smallness conditions (5) and (6) u 0 u 0 ε ν u 0 u 0 ε ν α, u0 u 0 ε ν α and u 0 3 ενα 5 4, then the solution u of ()-(3) is global: u (R + ; 3 ()). The important thing to observe here is that the smallness conditions in (6) involve only negative powers of α, so they disappear when α 0. Only the smallness condition (5) subsists, and this is in accordance with the classical global well-posedness results for small data for the Navier-Stokes equations. Note that condition (5) is scaling invariant and that it implies the smallness of the norm of u 0 by the interpolation inequality u 0 u 0 u 0. This in turn gives the existence of a unique global solution of ()-(4) by the celebrated result of Fuita and Kato [0]. We have the following immediate corollary. Corollary 3. Suppose that u 0 3 () is divergence free and vanishes on the boundary. There exists some small constant ε = ε() depending only on and some constant α 0 > 0 such that if u 0 u 0 ε ν and α α 0 then the 3 solution of ()-(3) is globally defined. 3
Combining this corollary with Theorem we immediately obtain the global convergence of solutions of () towards solutions of the Navier-Stokes equations if we assume that the initial velocity u 0 does not depend on α, belongs to 3 and is small in. This is ust an example of convergence result, more general results can be obtained by allowing the initial data to depend on α and combining Theorems and. Finally, let us conclude our results on global solutions with the observation that the smallness conditions (5) and (6) are implied by the simpler but less general condition: u 0 εν, u 0 ενα and u 0 3 ενα 5 4. Indeed, if the above relation holds true one can readily obtain (5) and (6) simply by estimating u 0 u 0. We state now our result on the uniform local existence of solutions. Theorem 4. There exist two constants ε = ε() and K = K() depending only on such that if u 0 3 () is divergence free and vanishes on the boundary and (7) u 0 ενα 4 and u 0 3 ενα 5 4 then there exists a 3 solution u (0,T; 3 ()) of ()-(3) which is defined at least up to the time T = K( u 0 + α u 0 ) 4. et us emphasize again that in the above theorem the time of existence of the solution depends on the norm of the initial data, but not on the 3 norm or the W,p, p > 3, norm as in the previous results on local existence of solutions. Another important thing to note is that both the smallness conditions for u 0 and u 0 3 involve negative powers of α. So they disappear when taking the limit α 0. More precisely, we have the following corollary. Corollary 5. Suppose that u 0 3 () is divergence free and vanishes on the boundary. et T α be the maximal time of existence of the 3 solution u of ()-(3), i.e. Then u loc ([0,T α); 3 ())\ (0,T α ; 3 ()). lim inf α 0 T α K u 0 4. Combining this corollary with Theorem yields, under the assumptions of the corollary, the convergence of the solutions of () towards the solutions of the Navier-Stokes equations up to the time T = ν3 K u 0. 4 3. Global existence of solutions InthissectionweshowTheorem. Inthesequel, wewilldenotebyc agenericconstant which depends only on and whose value can change from one line to another. We will use the standard notation for the m norms ( u m = β N 3, β m ) β u(x) dx. We denote by P the eray proector, i.e. the orthogonal proector on the subspace of divergence free and tangent to the boundary vector fields. 4
We will prove that there exists a sufficiently small constant ε = ε () depending only on such that if (8) u ν and u 3 ε ν α on some time interval [0,T], then the two quantities above are even smaller: (9) u ν and u 3 ε ν 4α on the same time interval [0,T]. By time continuity, this implies that (8) never breaks down if it is verified at the initial time and all the estimates that follow hold true globally intime. etusnowobserve that(8)holdstrueattheinitialtimeifεandε aresufficiently small independently of α and ν. We use the Gagliardo-Nirenberg inequality and (5) to write u 0 3 C u 0 u 0 Cεν < ε ν provided that Cε < ε. This is a condition that we assume in what follows. Next, we use an interpolation inequality and relation (6) to estimate u 0 u 0 u 0 ενα 4. From the Gagliardo-Nirenberg inequality and (6) we deduce that u 0 C u 0 4 u 0 3 4 C ε α < ν α provided that ε < which is another condition that we will assume in the sequel. C We conclude from the above estimates that (8) holds true at the initial time with strict inequality provided that Cε < ε and ε <. C We assume in the following that (8) holds true. Because we need a precise dependence on α of the constants, we have to split the estimates for u in three different parts: estimates, estimates and 3 estimates. The estimates below should be viewed as a priori estimates. They can be turned into rigorous estimates via a standard Galerkin approximation procedure. estimates. We multiply () by u and integrate in space to obtain t( u +α u )+ν u = 0 so, after an integration in time, t (0) u(t) +α u(t) +ν u(s) ds = u 0 +α u 0. 0 estimates. To perform estimates on () we need to write it under a different form. First we remove the pressure by applying the eray proector P. We obtain the following equivalent equation: t (u αp u) νp u+p [ u (u α u) ] +P [ ] (u α u) u = 0. We observe first that Next, we have that P ( u q + u u ) = P ( u ) = 0. q u = (u q) 5
is a gradient so ( P u q + q u ) = 0 for any q. Recalling that P u u is a gradient, we infer that () can be written under the following equivalent form: t (u αp u) νp u+p [ u (u αp u) ] αp [ ] (P u) u = 0. We multiply the above equation by P u and integrate in space. We recall that P is a self-adoint proector. We have that P = P and Pu = u because u is divergence free and tangent to the boundary. We infer that t( u +α P u )+ν P u = P [ u (u αp u) ] P u α = Using the cancellation we infer that u (u αp u) P u α u P u P u = 0 () t( u +α P u )+ν P u = We bound () α and P [ (P u) u ] P u (P u) u P u. u u P u α (P u) u P u α P u u ν P u u u P u u 3 u 6 P u C u 3 u P u C u 3 P u Cε ν P u ν 4 P u (P u) u P u. provided that ε /4C which is the only condition we will impose on ε. We used above the Sobolev embedding 6 and the classical regularity result for the stationary Stokes operator which claims that P u u. We infer from the previous relations that so, after an integration in time, (3) u(t) +α P u(t) + ν t ( u +α P u )+ ν P u 0 t 0 P u(s) ds u 0 +α P u 0. 6
3 estimates. et us introduce the notation ω = curlu, ω α = ω α ω. To performthe 3 estimates, we applythecurl operatorto(). Weobtainthefollowing PDE: t ω α ν ω +u ω α ω α u = 0. We multiply by ω α and integrate in space to obtain t ω α ν ω ω α = ω α u ω α ω α u ν α ω α. Next, since ω α = ω α ω we can write ν ω α = ω ν (ω ω α ) ω α = 3ν α 4α ω α + ν ω α α ω ν α ω. We infer that t ω α + 3ν 4α ω α + ν ω α α ω ν α ω ν α ω α so t ω α + ν α ω α ν α ω = ν α u. The Gronwall inequality implies that (4) ω α (t) ω α(0) e t ν α + ν α t 0 e (s t) ν α u(s) ds ω α (0) ν ν t e t α + α sup u(s) e (s t) ν α ds s [0,t] 0 ω α (0) +4 sup u(s) s [0,t] ω α (0) +4 u 0 +4α P u 0 where we used the bound for u given in (3). Adding this relation to (3) implies that (5) F(t) CF(0) where The following lemma holds true: emma 6. We have that with constants depending only on. F(t) = u(t) + α P u(t) + ω α (t). F(t) u(t) +α u(t) 3 Proof. Because P is bounded in and ω α = curlu α curlu we immediately get the bound F C( u + α u +α u 3). From the interpolation inequality u u u we deduce that α u 3 u +α u 3 so F C( u +α u 3). 7
To prove the reverse bound we observe first by the Poincaré inequality that u u. Therefore, it suffices to show that F Cα u 3. Observing that u = ω we can write ω α = ω α ω α ω ω = α curl u u. Next, we recall that u P u is a gradient so curl u = curlp u. We infer from the previous relations that F α P u +α curlp u α( P u + curlp u ) where we used that α. Since P u is divergence free and tangent to the boundary, one can apply [9, Proposition.4] to deduce that P u + curlp u C P u. Finally, the classical regularity result for the stationary Stokes operator says in particular that P u C u 3. We infer that F Cα u 3 and this completes the proof of the lemma. We go back to the proof of Theorem. From the Gagliardo-Nirenberg estimate u C u 4 u 3 4 3, from relation (5) and from emma 6 we infer that α u(t) Cα 4 F(t) Cα 4 F(0) C(α 4 u0 +α 5 4 u0 3). From the Gagliardo-Nirenberg and Poincaré inequalities and from relations (0) and (3) we can bound u 3 C u u C u u C( u 0 +α u 0 ) 4 ( u0 +α P u 0 ) 4. We conclude that for (9) to be verified it suffices to assume the following conditions on the initial data: (6) and C(α 4 u0 +α 5 4 u0 3) ν 4 (7) C( u 0 +α u 0 ) 4 ( u0 +α P u 0 ε ν ) 4. We prove now that if our smallness assumptions (5) and (6) hold true for a sufficiently small constant ε, then (6) and (7) are verified. Observe first that (6) implies (6). Indeed, this follows immediately from the interpolation inequality u 0 u 0 u 0 To prove (7), we notice that it is equivalent to u 0 u 0 +α u0 +α u0 u 0 +α u 0 u 0 ε ν C for some constant C. The necessary bounds for the first, the third and the fourth term on the left-hand side are included in the hypothesis (5) and (6) for ε small enough. The bound for the second term shows up in (6), so it is already proved. We conclude that if (5) and (6) hold true for a sufficiently small constant ε, then (8) never fails to be true and the solution u exists globally. This completes the proof of Theorem. 8.
4. Uniform local existence of solutions We show in this section Theorem 4. To prove the uniform local existence of solutions, we will adapt the estimates from the previous section in the following manner. We search for a constant M and a time T such that if we assume that (8) u M and u ν α on some sub-interval [0,T ] [0,T] then (9) u M and u ν 4α on the same sub-interval [0,T ]. By time continuity it follows that relation (9) holds true on [0,T] if it holds true at the initial time. As in the previous section one easily checks that the second condition from (8) is verified at time t = 0 if ε is sufficiently small. Indeed, the condition on the ipschitz norm of u 0 is the same as in (8) and in Section 3 we proved it using bounds for u 0 and u 0 3 which are the same as those from (7). For the first bound in (8) to hold true at time t = 0 we impose the following condition on the constant M (recall that M will be chosed later): u 0 < M. So let us assume that (8) holds true and let us prove relation (9). We adapt the estimates from the previous section. The estimates are not needed anymore. The estimates must be modified. Relations () and () remain valid and imply that (0) t ( u +α P u )+ν P u u u P u. The estimate of the right-hand side must be modified in the following manner u u P u u u P u C u 3 u P u C u 3 P u 3 ν P u + C u 6. for some constant C. We used the Poincaré and Young inequalities, the Gagliardo- Nirenberg estimate u C u u and the relation u C P u. Using this estimate in (0) yields t ( u +α P u )+ ν P u C u 6. Recalling that we assumed u M and integrating in time yields u(t) +α P u(t) u 0 +α P u 0 + C tm 6 If we make the assumption that () u 0 +α P u 0 + C TM 6 then we get that M 4 () u(t) +α P u(t) M 4. 9
In particular, we have that u(t) M which implies the first half of (9). The 3 estimates remain valid up to the relation (4) that we recall now ω α (t) ω α(0) +4 sup u(s). s [0,t] We can further bound u(s) M to obtain that Adding this estimate to () yields ω α (t) ω α(0) +4M. F(t) CF(0)+CM. We observed in the previous section that u Cα 3 4F so we deduce that u(t) Cα 3 4 (F(0)+M) C (α 3 4 u0 +α 4 u0 3 +α 3 4 M) for some constant C, where we used emma 6. We conclude that if we assume () together with (3) C (α 3 4 u0 +α 4 u0 3 +α 3 4 M) ν 4α then (9) holds true. To satisfy the condition () it suffices to assume that (4) u 0 M 3 and to choose T such that that is and C TM 6 (5) T = α P u0 M = M C M 4. 3 On the other hand, to ensure that (3) holds true it suffices to assume that ν ν ν (6) u 0, u C α 0 3 and M. 4 C α 5 4 C α 4 Clearly, to be able to find a constant M such that (4) and (6) are satisfied we need to assume the following conditions on the initial data: (7) u 0 ενα 4, u0 ενα 3 4 and u 0 3 ενα 5 4 for a sufficiently small ε. In fact, the condition on the norm is not necessary since the interpolation inequality shows that it follows from the other two 3 conditions. Once (7) is satisfied, we can choose M = C 3 ( u 0 + α u 0 ) for some constant C 3. According to (5), this gives a time of existence of the solution of the form T = C C3( u 4 0 + α u 0 ) 4. This completes the proof of Theorem 4. 0
Acknowledgments This work was supported by the ABEX MIYON (ANR-0-ABX-0070) of Université de yon, within the program Investissements d Avenir (ANR--IDEX-0007) operated by the French National Research Agency (ANR). References [] N. Arada. On the convergence of the two-dimensional second grade fluid model to the Navier Stokes equation. Journal of Differential Equations, 60(3):557 586, 06. [] D. Bresch and J. emoine. On the existence of solutions for non-stationary second-grade fluids. In Navier-Stokes equations and related nonlinear problems (Palanga, 997), pages 5 30. VSP, Utrecht, 998. [3] A. V. Busuioc, D. Iftimie, M. C. opes Filho and. J. Nussenzveig opes. Incompressible Euler as a limit of complex fluid models with Navier boundary conditions. Journal of Differential Equations, 5():64 640, 0. [4] A. V. Busuioc, D. Iftimie, M. C. opes Filho and. J. Nussenzveig opes. Uniform time of existence for the alpha Euler equations. Journal of Functional Analysis, 7(5):34 375, 06. [5] A. V. Busuioc and T. S. Ratiu. The second grade fluid and averaged Euler equations with Navier-slip boundary conditions. Nonlinearity, 6(3):9 49, 003. [6] D. Cioranescu and V. Girault. Weak and classical solutions of a family of second grade fluids. International Journal of Non-inear Mechanics, 3():37 335, 997. [7] D. Cioranescu and E.. Ouazar. Existence and uniqueness for fluids of second grade. In Nonlinear partial differential equations and their applications. Collège de France seminar, Vol. VI (Paris, 98/983), pages 78 97. Pitman, Boston, MA, 984. [8] J. E. Dunn and R.. Fosdick. Thermodynamics, stability, and boundedness of fluids of complexity and fluids of second grade. Arch. Rational Mech. Anal., 56:9 5, 974. [9] C. Foias and R. Temam. Remarques sur les équations de Navier-Stokes stationnaires et les phénomènes successifs de bifurcation. Ann. Scuola Norm. Sup. Pisa, Ser. IV, 5():9 63, 978. [0]. Fuita and T. Kato. On the Navier-Stokes initial value problem. I. Archive for Rational Mechanics and Analysis, 6(4):69 35, 964. [] G. P. Galdi, M. Grobbelaar-Van Dalsen and N. Sauer. Existence and uniqueness of classical solutions of the equations of motion for second-grade fluids. Archive for Rational Mechanics and Analysis, 4(3): 37, 993. [] G. P. Galdi and A. Sequeira. Further existence results for classical solutions of the equations of a second-grade fluid. Archive for Rational Mechanics and Analysis, 8(4):97 3, 994. [3] D. Iftimie. Remarques sur la limite α 0 pour les fluides de grade. In Nonlinear partial differential equations and their applications. Collège de France Seminar, Vol. XIV(Paris, 997/998), volume 3 of Stud. Math. Appl., pages 457 468. North-olland, Amsterdam, 00. [4] J. S. inshiz and E. S. Titi. On the convergence rate of the Euler-α, an inviscid second-grade complex fluid, model to the Euler equations. Journal of Statistical Physics, 38(-3):305 33, 00. [5] M. C. opes Filho,. J. Nussenzveig opes, E. S. Titi and A. Zang. Convergence of the D Eulerα to Euler equations in the Dirichlet case: Indifference to boundary layers. Physica D: Nonlinear Phenomena, 9 93:5 6, 05. (A.V. Busuioc) Université de yon, Université de Saint-Etienne CNRS UMR 508 Institut Camille Jordan Faculté des Sciences 3 rue Docteur Paul Michelon 403 Saint-Etienne Cedex, France E-mail address: valentina.busuioc@univ-st-etienne.fr