arxiv: v2 [nlin.cd] 9 Mar 2012

Transcription

1 Conditional regularity of solutions of the three diensional Navier-Stokes equations and iplications for interittency J. D. Gibbon Departent of Matheatics, Iperial College London SW7 2AZ, UK arxiv: v2 [nlin.cd] 9 Mar 212 Two unusual tie-integral conditional regularity proofs are presented for the three-diensional Navier-Stokes equations. The ideas are based on L 2 -nors of the vorticity, denoted by Ω t, and particularly on = Ω α, where α = 2/4 3 for 1. The first result, ore appropriate for the unforced case, can be stated siply: if there exists an 1 < and an < ε < 2 such that Dε +1dτ c ε, Dε dτ with c ε, > 1 then no singularity can occur on [, t]. Secondly, for the forced case, by iposing a critical lower bound on Ddτ, no singularity can occur in Dt for large initial data. Moveent across this critical lower bound shows how solutions can behave interittently, in analogy with a relaxation oscillator. Potential singularities that drive Ddτ over this critical value can be ruled out whereas others cannot. I. INTRODUCTION The twin related thees of this paper are: i the regularity proble for solutions of the three-diensional incopressible Navier-Stokes equations; and ii the interittent behaviour of these solutions. Traditionally the first proble, which still reains tantalizingly open, has lain in the doain of the analyst, whereas the phenoenon of interittency has tended to be ore of interest to the physics and engineering fluid dynaics counities. It will be deonstrated in this paper that these two issues are intiately related & require siultaneous study. A. History Forally, a weak solution ux, t of the three-diensional Navier-Stokes equations t +u u = ν u p+fx, 1 with divu =, is called regular if the H 1 -nor is continuous [1]. What is coonly referred to as conditional regularity can be achieved if it is found necessary to ipose assuptions on certain syste variables such as the velocity field. The early work of Prodi [3], Serrin [4] and Ladyzhenskaya [5] can be suarized thus [6, 7]: every Leray-Hopf solution of the incopressible three-diensional Navier-Stokes equations with u L r, T;L s is regular on, T] provided 2/r+3/s = 1, with s 3, ], or if u L, T;L p with p > 3. The long-standing case s = 3 was finally settled by von Wahl [8] and Giga [9] who proved regularity in the space C, T];L 3 : see also Kozono and Sohr [1] and Escauriaza, Seregin and Sverák [11]. In suary, the s = 3 case sees tantalizingly close to the bounded case s = 2, but not quite close enough. More recent results exist that ipose conditions on either the pressure or on one derivative of the velocity field: see references in Kukavica and Ziane [12, 13], Zhou [14], Cao & Titi [15, 16], Cao [17], Cao, Qin and Titi for channel flows [18], Chen and Gala [19] and the review by Doering [2]. Results on the direction of vorticity can be found in Constantin and Fefferan [21] and asseur [22], and those on the use of Besov spaces in Cheskidov and Shvydkoy [23]. Recent work by Biswas and Foias [24] considers analyticity properties of Navier-Stokes solutions in which they study the axial space analyticity radius associated with a regular solution involving Gevrey-class nors. B. Motivation and notation In an entirely different thread of intellectual endeavour, the seinal experiental paper of Batchelor and Townsend [25] investigated the phenoenon of interittency in wind tunnel turbulence by observing that the flatness of their signals the ratio of the fourth order oent to the square of the second order oent took uch higher values than Electronic address: j.d.gibbon@ic.ac.uk; URL:

2 flatreg3.tex 9th/3/12 For the Constantin 6th birthday volue 2 those expected for rando Gaussian behaviour. They concluded that the vorticity is neither spatially nor teporally distributed in an even fashion but undergoes local clustering or spottiness, which is consistent with spiky behaviour in signals interspersed by longer quiescent periods. This is considered to be a classic characteristic of interittency. These ideas have been confired, developed and extended in any subsequent experients and coputations: see the papers by Kuo and Corrsin [26], Sreenivasan [27], Meneveau and Sreenivasan [28] and the books by Frisch [29] and Davidson [3] for further references. Most of these discussions have been based around Kologorov s statistical theory with the widespread use of velocity structure functions to study interittent behaviour. However, structure functions are not easily translatable into results in Navier-Stokes analysis see soe of the arguents in Kuksin [31]. The ain difficulty lies in translating the special conditions in I A into sensible physics while conversely aking sense of the experiental observations in ters of Navier-Stokes variables. The two threads can be erged if the spiky nature of the vorticity field is considered in the context of L p -nors of the vorticity ω = curlu with p = 2 Ω t = L 3 ω 2 d 1/2 +, 1, 2 on a periodic box [, L] 3. The additive frequency = νl 2 is present for technical reasons. Clearly the Ω t are ordered for all t such that Ω 1 t Ω 2 t... Ω t Ω +1 t..., 3 where L 3/2 Ω 1 t is the H 1 -nor. Control fro above on any one of the Ω will also control the H 1 -nor fro above which is the ultiate key to regularity. The Navier-Stokes equations have a well-known invariance property under the transforations x = ǫx, t = ǫ 2 t, u = ǫu and p = ǫ 2 p. Under these Ω scales as Ω α = ǫω α, α = Thus it is natural to define t = [Ω t] α. 5 While the nors Ω are ordered with increasing as in 3, the α decrease with. Thus there is no natural ordering aong the. Note also that D 1 is the square of the H 1 -nor. For the forced case, the diensionless Grashof nuber Gr is based on the bounded-ness of the root-ean-square frs 2 = L 3 f 2 2 of the divergence-free forcing fx and is defined as Gr = L3 f rs ν 2. 6 C. Two fundaental results There are two results that for the basis of those given in later sections. The first is a theore on tie integrals or averages of [32] which uses a result of Foias, Guillopé and Tea [33]. This proof will not be repeated: Theore 1 For 1 and α defined as α = weak solutions obey dτ c tgr 2 +η 1, 7 where η 1 = Lν 3 E and E is the initial energy. In the unforced case the right hand side is just cη 1. In [32] this result was converted into a set of length scales. Let the tie average up to tie T be defined by T = lisup D 1 T τdτ 8 T

3 flatreg3.tex 9th/3/12 For the Constantin 6th birthday volue 3 in which case 7 can be re-expressed as T cgr 2 +O T 1. 9 Then, otivated by the definition of the Kologorov length for = 1, a set of length scales can be defined thus: Lλ 1 2α := D T. 1 In [32] the bounds in 9 were expressed as Re 3 instead of Gr 2 based on a device of Doering and Foias [34] who used the definition of the velocity U 2 = L 3 u 2 2 T within the Reynolds nuber Re = U Lν 1 ; this led to the inequality Gr cre 2. However, the preference in this paper is to reain with the Grashof nuber Gr in this paper; as a rough guide Gr 2 can be replaced by Re 3. In ters of Re equation 1 becoes Lλ 1 cre 3 2α. 11 When = 1, α 1 = 2 and thus Lλ 1 1 cre 3/4, which is consistent with Kolorgorov s statistical theory [29]. The second result is a differential inequality for the. Atie-differentiation of the vorticityfield createsprobles because only weak solutions exist. Circuvention of this difficulty requires a contradiction strategy coonly used in geoetric analysis: assue that there is a axial interval of existence and uniqueness [, T, which iplies that H 1 T =. In any subsequent calculation, if the H 1 -nor turned out to be bounded in the liit t T, then a contradiction would result and so the interval [, T could not be axial. Moreover, it cannot be zero, so T would have to be infinite. Define three frequencies 1, = α c 1 1, 2, = α c 2, 3, = α c 3,. 12 where the constants c n, n = 1,2,3 are algebraically increasing with. With the dot representing differentiation with respect to tie, we have: Theore 2 For 1 < on [, T, the t forally satisfy the set of inequalities { ρ Ḋ D 3 D+1 1, + 2,}, 13 where ρ = In the forced case there is an additive ter 3,Gr. Proof: The proof requires soe variations on a previous result [35] and is relegated to Appendix A. The rest of the paper is based on two integrations of the differential inequality for given in Theore 2. The first integration in II, ore appropriate for the unforced case, produces a conditional regularity result based on the relative sizes of dτ and +1dτ and ay serve as an alternative to the Beale-Kato-Majda theore [36]. The result in III, ore appropriate for the forced case, is based on an assued lower bound on +1dτ in ters of the Grashof nuber Gr. Both cases take advantage of the bounded tie integrals of fro Theore 1. II. A CONDITIONAL REGULARITY RESULT FOR UNFORCED NAIER-STOKES A. First integration of the inequality An iediate conclusion that can be drawn fro Theore 2 is that solutions coe under control pointwise in t provided +1 t c ρ t, 14 where c ε,ρ = [c 1, c 2, ] ε/ρ. However this can be bettered through a foral solution of 13 which is displayed in the Theore below. Theore 3 For any value of 1 < and ε unifor in the range < ε < 2, if the integral condition is satisfied then t obeys t on the interval [, t]. D+1dτ ε c ε,ρ Ddτ ε 15

4 flatreg3.tex 9th/3/12 For the Constantin 6th birthday volue 4 Reark: If 15 holds then clearly the Navier-Stokes equations are regular. While ε being allowed arbitrarily close to zero ay suggest the use of logariths, the case ε = is analyzed in III. Proof: Divide 13 by D 3 ε and integrate to obtain [ t] ε 2 [ ] ε 2 ε 1, D+1 ρ D ε dτ ε 2, where ε n, = 2 ε n,. Noting that ρ 1/3, a Hölder inequality then easily shows that [ ρ ] ε D+1 ε dτ = D+1 D ε ρ [D ε ] ρ ε ρ dτ 16 can then be re-written as t [] ε 2 + ε 1, D+1 ρ ε Ddτ ε ρ t ρ ε Ddτ ε ρ t ρ/ε Dε +1 dτ t ρ ε/ε Dε dτ ε 2, Ddτ ε D ε dτ, ε Clearly no singularity can for if a zero cannot appear in the denoinator, the condition for which is 15. If this be the case then t. It follows that Ω t is bounded above and thus so is the H 1 -nor Ω 1 so the contradiction arguent explained above can be invoked. B. Range of scales The conclusion that can be drawn fro Theore 3 is that it needs to hold for only one value of for the Navier- Stokes equations to be regular. For ε = 1 and converting the integrals of and +1 into the length scales λ defined in equation 1 turns the condition for the regular regie into [32] Lλ 1 2α+1 +1 cρ Lλ 1 2α 19 which can be re-written as c ρ λ +1 λ α λ α 1 +1 L where c ρ = cρ 1/2α+1. Given the sall size of the right hand side and the positivity of the exponent in 2 iplies that λ +1 is squeezed down in size copared to λ. If the Theore holds for a string of values of then the corresponding λ would represent a decreasing cascade of length scales representing a regular flow. Because Theore 3 needs to be satisfied for only one value of to obtain control of the Navier-Stokes equations this eans that the regularity proble is only open if, for all values of 1 <, the reverse of 2 is true. The final average on the right hand side of the reverse of 19 is Lλ 1 2α... c1 Lλ c Gr 2 c Re 3 21 which siply eans that Lλ 1 c Re 3/2α, which is consistentwith the results oftheore1. This reversedirection produces no conclusion: no cascade is iplied as the λ can ix in sizes because the α decrease with increasing. The lack of a regularity proof for this regie eans that singularities are still potentially possible. III. SECOND INTEGRATION OF THE INEQUALITY Inclusion of the forcing in 13 odifies it to { Ḋ D 3 1 1, A. A lower bound on Ddτ D+1 ρ + 2,} + 3, Gr, 22

5 flatreg3.tex 9th/3/12 For the Constantin 6th birthday volue 5 where ρ = To proceed, divide by D3 the case ε = of Theore 3 to write 22 as where 1 d 2dt D 2 X D 2 2, 23 A lower bound for X dτ can be estiated in a anner siilar to?? ρ D+1 X = 1, D 2 3, Gr. 24 and so +1 dτ = [ D+1 D+1 ρ ] 1 D 2 ρ D ρ 2 ρ dτ ρ 1 Ddτ 2 ρ t ρ 2 ρ dτ t 1/ρ 25 X dτ 1, t 1 D ρ +1dτ dτ ρ 2 3,tGr integrates to [ t] 2 { exp 2 } t X dτ [ ] 2 2 2, exp{ 2 τ X dτ } dτ. 27 Let us recall that ρ = and let us also define then γ = α , 28 Theore 4 On the interval [, t] if there exists a value of lying in the range 1 < <, with initial data [ ] 2 < C Gr, for which the integral lies on or above the critical value where η 2 η 1 Gr 2δ 1 and where 1 4 c tgr 2δ t +η 2 dτ 29 = 4{δ 1+ρ γ ρ γ } then t decays exponentially on [, t]. with 1+4ρ γ 41+ρ γ < δ < 1, 3 Reark: δ ց 11/2 for large so enough slack lies between the upper and lower bounds on dτ. Proof: It is not difficult to prove that Ω 2 Ω Ω 1 for > 1, fro which it is easily found that +1dτ dτ dτ D 1dτ γ. 31

6 flatreg3.tex 9th/3/12 For the Constantin 6th birthday volue 6 Thus 26 can be re-written as t ργ +2 X dτ 1, t 1 dτ t ργ 3, tgr 1dτ c t 1, Gr 3, Gr 32 having used the assued lower bound in the theore and the upper bound of D 1dτ. Moreover, to have the dissipation greater than forcing requires > 1 so δ ust lie in the range as in 3 because 1 < 4. For large Gr the negative Gr-ter in 27 is dropped so the integral in the denoinator of 27 is estiated as τ exp 2 X dτ dτ [2 c 1, ] 1 Gr 1 exp [ 2 1, c tgr ], 33 and so the denoinator of 27 satisfies Denoinator [ ] 2 c 2, c 1, 2 c 1 Gr 1 exp [ 2 1, c tgr ]. 34 This can never go negative if [ ] 2 > c 1, c 2, 2 c 1 Gr, which eans < C Gr 1 2. B. A echanis for interittency A ajor feature of interittent flows lies in the strong, spiky, excursions of the vorticity away fro averages with periods of relative inactivity between the spikes. How do these aperiodic cycles appear in solutions and does the critical lower bound iposed as an assuption in Theore 4 lead to this? Using the average notation t, equation 27 shows that if t lies above critical then t collapses exponentially. Experientally, signals go through cycles of growth and collapse so it is not realistic to expect this critical lower bound to hold for all tie.. t.. Gr 2 upper bound of D t... critical lower-bound of D Gr 2δ... t potential singularity... t t t 1 t 2 Figure 1: A cartoon of t versus t illustrating the four phases of interittency. To understand interittency we turn to the figure and draw the horizontal line at Gr 2δ as the critical lower bound on t. Within this allowed range, t will decay exponentially fast. Because integrals ust take account of history, there will be a delay before t decreases below the value above which a zero in the denoinator of 27 can be prevented at t 1 : at this point all constraints are reoved and the pointwise solution t is free to grow rapidly again in the interval t 1 t t 2. Curiously, and perhaps counter-intuitively, if the value of this integral drops below critical then it is in this interval that the occurrence of a singular event ust still forally be considered. Then the growth in will be such that, after another delay, it will force t above critical and the syste, with a re-set of initial conditions at t 2, is free to go through another cycle as in a relaxation oscillator.

7 flatreg3.tex 9th/3/12 For the Constantin 6th birthday volue 7 I. CONCLUSION The results of Theores 3 and 4 are new. Taking the in reverse order, Theore 4 can be suarized thus: if dτ c Gr 2δ +η 2 then t collapses exponentially. The question that stands out is the eaning of the regie where 35 dτ < c Gr 2δ +η It ight be natural to suppose, on intuitive grounds, that potentially singular behaviour would be less likely to occur when the upper bound is suppressed, yet this behaviour cannot be wholly ruled out. What can be ruled out are potentially singular spikes that substantially contribute to the tie integral of t because, in driving up the value of the integral, they would push it over its critical value thereby forcing exponential collapse in t itself. There still reains the possibility of needle-like potentially singular spikes that contribute little or nothing to the tie integral of t. More appropriate for the unforced case, Theore 3 can be suarized thus: for < ε < 2 if D+1dτ ε c ε,ρ Ddτ ε 37 then t. The question that stands out is the eaning of the regie where D ε +1 dτ < c ε,ρ D ε dτ 38 Is the regular regie physically reachable? Fro fairly general initial conditions, not all flows fro would obey the consequence of Theore 3 that t. Nuerical experients ay decide which of the two regies are the ost anifest. Acknowledgents Thanks are due to Robert Kerr Warwick and Darryl Hol Iperial for discussions on this proble and to Mihaela Ifri UC Davis for a critical reading of the anuscript. Appendix A: Proof of Theore 2 The following proof has soe iportant variations on that given in [35]. Consider J = ω 2 d such that 1 J 2 = ω 2 1 ω {ν ω +ω u+curlf} d. A1 Bounds on the three constituent parts of A1 are dealt with in turn, culinating in a differential inequality for J. a The Laplacian ter: Let φ = ω 2 = ω ω. Then ω 2 1 ω ωd = } 1 φ { ω 1 2 φ 2 1 φ 1 2 φ d. Using the fact that φ = { 1φ 2 φ 2 +φ 1 φ } we obtain ω 2 1 ω ωd φ 2 φ 2 d = ω 2 d, d φ d A2 A3

8 flatreg3.tex 9th/3/12 For the Constantin 6th birthday volue 8 having used the Divergence Theore. Thus we have { ω 2 1 ω ωd ω 2 ]d = 1, 2 c 1, A 2 d 2. A4 where A = ω and c 1, = 2 / 1 with equality at = 1. The negativity of the right hand side of A4 is iportant. Both A 2 and A 2 will appear later in the proof. b The nonlinear ter in A1: After a Hölder inequality, the second ter in A1 becoes ω 2 u d c 1 1 u d ω d ω 2 2 d ω 2+1 d = c J 1/2 +1 J1/2, ω 2 2 d where the inequality u p c p ω p for p 1, has been used, which is based on a Riesz transfor: note the exclusion of the case =. Together with A2 this akes A1 into 1 J 2 ν ω 2 d +c J 1/2 +1 c J1/2 + ω 2 1 ω curlf d. A6 1, c The forcing ter in A1: Now we use the sallest scale in the forcing l with l = L/2π to estiate the last ter in A6 ω 2 1 ω curlf d ω f 2 A7 However, by going up to at least N 3 derivatives in a Sobolev inequality and using our restriction of single-scale forcing at k l 1 with l = L/2π it can easily be shown that f 2 c f 2 l and so ω 2 1 ω curlf d cω L Gr. A8 d A differential inequality for J : Recalling that A = ω allows us to re-write J +1 as A5 J +1 = A 2+1/ 2+1/. A9 A Gagliardo-Nirenberg inequality yields A 2+1/ c A 3/2+1 2 A 2 1/2+1 2 A1 which eans that With the definition J +1 c 3/2 ω 2 d J 2 1/2. A11 β = A12 the factor of 4 3 is different fro that in [35], A11 can be used to for Ω +1 1/2+1 Ω +1 = L 3 J c L 1 1/β ω 2 d + 2 [ L 3 J 1/2 + ] 2 1/2+1 A13

9 flatreg3.tex 9th/3/12 For the Constantin 6th birthday volue 9 which converts to c L 1 ω 2 d + 2 Ω+1 Ω β Ω 2. A14 This otivates us to re-write A6 as 1 L 3 J 2 c 1, L 1 ω 2 d + c 2, L 3 J +1 1/2 L 3 J 1/2 + c 3, 2 Ω2 1 Gr. A15 Converting the J into Ω and using the fact that Ω Ω Ω { c 4, Ω+1 Ω 4+1/3 +c 5, Ω+1 Ω +1 Ω +c 6, Gr} A16 Using a Hölder inequality on the central ter on the right hand side, with the definition of β given in A12, A16 finally becoes { Ω Ω 1 } β Ω+1 1 +c 2, c 1, Ω Ω 2α +c3, Gr. A17 Given the definition = 1 Ω α it is found that having used the fact that Ω+1 Ω β = D+1 ρ D β = 2 α +1 α A18 A19 and ρ = β /α +1. This converts A17 to 13 of the theore with = excluded. The constants c n, for n = 1,2,3 grow algebraically with. [1] J. Leray, Acta Math., 63, [2] C. Cao and E. S. Titi, Ann. Math., 166, [3] G. Prodi, Ann. Mat. Pura Appl., 48, [4] J. Serrin, The initial value proble for the Navier-Stokes equations, Nonlinear Probles, Proc. Sypos., Madison, Wis. pp , edited R. E. Langer, Univ. of Wisconsin Press, Madison, Wis., [5] O. A. Ladyzhenskaya, Matheatical Theory of iscous Incopressible Flow, Gordon and Breach, New York, English translation, 2nd ed., [6] P. Constantin and C. Foias, The Navier-Stokes equations, Chicago University Press, Chicago, [7] C. Foias, O. Manley, R. Rosa and R. Tea, Navier-Stokes equations and turbulence, Cabridge University Press, Cabridge, 21. [8] W. von Wahl, Proc. Syp. Pure Math., 45, [9] Y. Giga, J. Diff. Eqn., 62, [1] H. Kozono and H. Sohr, Adv. Diff. Equns, 2, [11] L. Escauriaza, G. Seregin, and. Sverák, Russ. Math. Surveys, 58, [12] I. Kukavica and M. Ziane, Nonlinearity, 19, [13] I. Kukavica and M. Ziane, J. Math. Phys., 48, [14] Y. Zhou, Methods Appl. Anal., 9, [15] C. Cao and E. S. Titi, Indiana Univ. Math. J., 57, [16] C. Cao and E. S. Titi, Global regularity criterion for the 3D Navier-Stokes equations involving one entry of the velocity gradient tensor, arxiv: v1, 25th May, 21. [17] C. Cao, Discrete and Continuous Dynaical Systes, 26, [18] C. Cao, J. Qin and E. S. Titi, Co. PDEs, 33,

10 flatreg3.tex 9th/3/12 For the Constantin 6th birthday volue 1 [19] Wenying Chen, S. Gala, Elec. J. Diff. Equns., 211, No. 6, [2] C. R. Doering, Annu. Rev. Fluid Mech., 41, [21] P. Constantin and C. Fefferan, Indiana Univ. Math. J., 42, [22] A. asseur, Applications of Matheatics, 54, No. 1, [23] A. Cheskidov and R. Shvydkoy, Arch. Ration. Mech. Anal., 195, [24] A. Biswas and C. Foias, On the axial spatial analyticity radius for the 3D Navier-Stoke equations and turbulence, subitted 211. [25] G. K. Batchelor and A. A. Townsend, Proc R. Soc. Lond. A, 199, [26] A. Y-S. Kuo and S. Corrsin, J. Fluid Mech., 5, [27] K. Sreenivasan, J. Fluid Mech., 151, [28] C. Meneveau and K. Sreenivasan, J. Fluid Mech., 224, [29] U. Frisch, Turbulence: the legacy of A. N. Kologorov, Cabridge University Press, Cabridge, [3] P. A. Davidson, Turbulence, Oxford University Press, Oxford, 24. [31] S. B. Kuksin, Geo. Funct. Anal., 9, [32] J. D. Gibbon, Co Math. Sci., 1, [33] C. Foias, C. Guillopé and R. Tea, Co. PDEs. 6, [34] C. R. Doering and C. Foias, J. Fluid Mech, 467, [35] J. D. Gibbon, Proc. Royal Soc A, 466, [36] J. T. Beale, T. Kato and A. J. Majda, 94,