arxiv: v1 [math-ph] 21 May 2011

Transcription

1 Exact, Rotational, Infinite Energy, Blowup Solutions to the 3-Dimensional Euler Equations arxiv: v1 [math-ph] 1 May 011 Manwai Yuen Department of Applied Mathematics, The Hong Kong Polytechnic University, Hung Hom, Kowloon, Hong Kong Revised 1-May-011 Abstract In this paper, we construct a new class of blowup solutions with elementary functions to the 3-dimensional compressible or incompressible Euler and Navier-Stokes equations. In detail, we obtain a class of global rotational exact solutions for the compressible fluids with γ > 1: [ ρ = max [ x +y +z xy +yz +xz) ] ȧt)x+y +z)+bt) ] 1,0} u 1 = at)+y z) 1) u = at)+ x+z) u 3 = at)+x y) where and at) = c 0 +c 1t ) bt) = 3c 0c 1t+ 3 c 1t +c 3) with, c 0, c 1 and c are arbitrary constants; And the corresponding blowup or global solutions for the incompressible Euler equations are also given. Our constructed solutions are similar to the famous Arnold-Beltrami-hildress AB) flow. The solutions with infinite energy can exhibit the interesting behaviors locally. Besides, the corresponding global solutions are also given for the compressible Euler equations. Furthermore, due to div u = 0 for the solutions, the solutions also work for the 3-dimnsional incompressible Euler and Navier-Stokes equations. MS: 35B40 35Q D05 Key Words: Euler Equations, Exact Solutions, Rotational, ompressible, Incompressible, Symmetry Reductions, Blowup, Global Solutions, Infinite Energy, Navier-Stokes Equations, Free Boundary 1 Introduction The N-dimensional Euler equations can be formulated as the follows: ρ t + ρ u) = 0 ρ[ u t + u ) u]+ Pρ) = 0 4) where x = x 1,x,...,x N ) R N, ρ = ρt, x) and u = ut, x) = u 1,u,...u n ) R N are the density and the velocity respectively. The γ-law could be applied to the pressure function, i.e. address: nevetsyuen@hotmail.com Pρ) = Kρ γ 5) 1

2 Manwai Yuen with K 0 and γ 1. The Euler equations 4) can be rewritten by the scalar form, ρ t + N k=1 u k ρ x k +ρ N u k k=1 ) ρ ui t + N k=1 u k ui x k x k =0, + P x i = 0, for i = 1,,...N. 6) The Euler equations 4) are the fundamental model in fluid mechanics [] and [5]. onstructing exact solutions is a very important part in mathematical physics to understand the nonlinear behaviors of the system. For the pressureless fluids K = 0, a class of exact blowup solutions were given in Yuen and Yeung s papers [14] and [11]: ρ = fx 1+d 1 a 1 t),x +d a t),...,x N +d N a N t) ), u N i = ȧit) Π i=1 ait) a x it) i +d i ) for i = 1,,...,N 7) a i t) = a i0 +a i1 t with an arbitrary function f 0 and f 1 ; and d i, a i0 > 0 and a i1 are constants. In particular, for a < 0, the solutions blow up in the finite time T = a 1 /a. There are other analytical solutions for the compressible Euler or Navier-Stokes equations in [9], [8], [6], [1], [13], [14], [11], [15] and [16]. In addition, there are two papers to investigate for a special structural form which generalizes the exact solutions for Burgers vortices in [7] and [3]. In this paper, we manipulate the elementary functions to construct some exact rotational solutions for the 3-dimensional compressible Euler equations 4): Theorem 1 For the 3-dimensional compressible Euler equations 4), there exists a class of rotational solutions: for γ > 1: ρ = max [ [ x +y +z xy +yz +xz) ] ȧt)x+y +z)+bt) ] 1, 0} u 1 = at)+y z) 8) u = at)+ x+z) u 3 = at)+x y) where and at) = c 0 +c 1 t 9) bt) = 3c 0 c 1 t+ 3 c 1t +c 10) with, c 0, c 1 and c are arbitrary constants; for γ = 1: ρ = e [x +y +z xy+yz+xz)] ȧt)x+y+z)+bt) u 1 = at)+y z) u = at)+ x+z) u 3 = at)+x y). The solutions 8) and 11) globally exist. 11) Remark In 1965, Arnold first introduced the famous Arnold-Beltrami-hildress AB) flow u 1 = Asinz +cosy u = Bsinx+Acosz 1) u 3 = siny +Bcosx with constants A, B, and a suitable pressure function P only for the incompressible Euler equations in [1]. We observe that our solutions 8) and 11) are similar to the AB flow.

3 Exact Solutions to the 3-dim Euler Equations 3 Remark 3 The solutions with infinite energy can exhibit the interesting behaviors locally. The exact solutions with infinite energy of the systems may be regionally applicable to understand the great complexity that exists in turbulent phenomena. Remark 4 We notice that the rational functional form with at) = 0 and = 1 in solutions 8) and 11) for the velocity u comes from Senba and Suzuki s book [10]. The velocities u in solutions 8) and 11) are not spherically symmetric. Remark 5 The exact rotational solutions 8) and 11) could be good examples for testing numerical methods for fluid dynamics. ompressible Rotational Fluids The main technique of this article is just to use the primary assumption about the velocities u: u 1 = at)+y z) u = at)+ x+z) 13) u 3 = at)+x y) to substitute the governing equations 4) to construct the density function ρ to balance the system. Therefore, the proof is simple to be checked by direct computation: Proof. For the vacuum solutions ρ = 0, it is the trivial solutions for the system 4). For γ > 1, with non-vacuum solutions, we have the first momentum equation 4),1 : with t u 1 +u 1 x u 1 +u y u 1 +u 3 z u 1 +ρ γ x ρ 14) = t u 1 +u y u 1 +u 3 z u 1 + γ 1 x ρ 15) γ 1 ρ = x +y +z xy +yz +xz) ) ȧt)x+y +z)+bt) 16) which we could require the condition for the density function: γ 1 [ ρ = max [ x +y +z xy +yz +xz) ] ȧt)x+y +z)+bt) ] } 1, 0, 17) to have = ȧt)+[at)+ x+z)] y [at)+y z)]+[at)+x y)] [at)+y z)] z + [ [ x +y +z xy +yz +xz) ] ȧt)x+y +z)+bt) ] 18) x = ȧt)+[at)+ x+z)] +[at)+x y) )] + [ x y z ȧt) ] 19) = ȧt) x+ z + y + x y z ȧt) 0) = 0. 1) Similarly by the nice symmetry of the functions 8), we could balance the second momentum equation 4), : t u +u 1 x u +u 3 z u + γ 1y ρ )

4 4 Manwai Yuen = ȧt)+[at)+y z)] x [at)+ x+z)]+[at)+x y)] z [at)+ x+z)] + [ [ x +y +z ) xy +yz +xz) ] ȧt)x+y +z)+bt) ] 3) y = ȧt)+[at)+y z)] )+[at)+x y)] 4) + y x z ȧt) For the last equation 4),3, we have t u 3 +u 1 x u 3 +u y u 3 + γ 1 = 0. 5) z ρ 6) = ȧt)+[at)+y z)] x [at)+x y)]+[at)+ x+z)] y [at)+x y)] + [ [ x +y +z ) xy +yz +xz) ] ȧt)x+y+z)+bt) ] 7) z = ȧt)+[at)+y z)] +[at)+ x+z)] )+ z y x ȧt) 8) For the mass equation 4) 1, we have with the density 17) to obtain: = 0. 9) ρ t + ρ u) 30) = [ γ 1 t + [ γ 1 x + [ γ 1 y + [ γ 1 x +y +z xy +yz +xz) ) ȧt)x+y +z)+bt) )] 1 z x +y +z xy +yz +xz) ) ȧt)x+y +z)+bt) )] 1 x +y +z xy +yz +xz) ) ȧt)x+y +z)+bt) )] 1 at)+y z)) x +y +z xy +yz +xz) ) ȧt)x+y +z)+bt) )] 1 at)+ x+z)) at)+x y)) 31) = 1 γ = 1 [ γ 1 x +y +z xy +yz +xz) ) ȧt)x+y +z)+bt) )] 1 1 [ ät)x+y +z)+ḃt) ] [ x y z ȧt) ) at)+y z)) ] [ y x z ȧt) ) at)+ x+z)) ] [ z x y ȧt) ) at)+x y)) ] [ γ 1 x +y +z xy +yz +xz) ) ȧt)x+y +z)+bt) )] 1 1 ät)x+y +z)+ḃt) + x y z ȧt) ) at)+ x y z ȧt) ) y z) + y x z ȧt) ) at)+ y x z ȧt) ) x+z) + z x y ȧt) ) at)+ z x y ȧt) ) x y) 3) 33)

5 Exact Solutions to the 3-dim Euler Equations 5 = 1 [ γ 1 x +y +z xy +yz +xz) ) ȧt)x+y +z)+bt) )] 1 1 } ät)x+y +z)+ḃt) 3ȧt)at) 34) = 0 35) by choosing and at) = c 0 +c 1 t 36) ḃt) = 3ȧt)at) 37) where c is the arbitrary constant. As for γ = 1, the proof is similar, we may skip the details here. It is clear to see that the solutions 8) and 11) globally exist. The proof is completed. Here the masses of the solutions 8) and 11) are infinite: 3 Incompressible Rotational Fluids ḃt) = 3c 1 c 0 +c 1 t) 38) bt) = 3c 0 c 1 t+ 3 c 1t +c 39) R 3 ρdx = +. 40) For the 3-dimensional incompressible Euler equations: div u = 0 ρ[ u t + u ) u]+ Pρ) = 0, 41) the solutions with elliptical symmetric velocity in [9], [8] and [16] for the compressible Euler equations 4), could not be applied to the incompressible Euler equations 41) because of the linear functional form: u 1 = c 1 t)x+d 1 t) u = c t)y +d t) 4) u 3 = c 3 t)z +d 3 t) with some functions c i t) and d i t) to have div u 0 43) for non-trivial solutions c 1 t) + c t) + c 3 t) 0. However, our solutions 8) and 11) could be applied to the incompressible Euler equations 41). The set of the corresponding solutions for the incompressible Euler equations 41) could be larger than the compressible ones 4): Theorem 6 For the 3-dimensional incompressible Euler equations 41), there exists a class of rotational solutions: [ [ Pρ) x ρ = +y +z xy +yz +xz) ] ] [ȧ 1 t)+a t)+a 3 t))]x [ȧ t)+a 1 t)+a 3 t))]y [ȧ 3 t)+a 1 t)+a t))]z +bt) u 1 = a 1 t)+y z) 44) u = a t)+ x+z) u 3 = a 3 t)+x y)

6 6 Manwai Yuen and Pρ) ρ = x +y +z +xy +yz +xz ) [ȧ 1 t)+a t)+a 3 t))]x [ȧ t)+a 1 t)+a 3 t))]y [ȧ 3 t)+a 1 t)+a t))]z +bt) u 1 = a 1 t)+y +z) u = a t)+x+z) u 3 = a 3 t)+x+y) where a i t) is an arbitrary local 1 function for i = 1, or 3, bt) is an arbitrary function, is an arbitrary constant. In particular, 1. if a i T) = or ȧ i T) = with the first finite positive constant T, the solutions 44) and 45) blow up in the finite time T;. For global 1 functions a i t), the solutions 44) and 45) globally exist. Proof. For the checking of the functions 44), the detail is similar to the proof of Theorem 1: For the checking of the functions 45), we have for the first momentum equation 41), : t u 1 +u 1 x u 1 +u y u 1 +u 3 z u 1 + ) P 46) x ρ = t u 1 +u y u 1 +u 3 z u 1 + ) P 47) x ρ by requiring the pressure function: to have P ρ = [ [ x +y +z xy +yz +xz) ] [ȧ 1 t)+a t)+a 3 t))]x [ȧ t)+a 1 t)+a 3 t))]y [ȧ 3 t)+a 1 t)+a t))]z +bt) = ȧ 1 t)+[a t)+ x+z)] y [a 1t)+y z)]+[a 3 t)+x y)] z [a 1t)+y z)] + [ [ x +y +z xy +yz +xz) ] ] [ȧ 1 t)+a t)+a 3 t))]x 49) x [ȧ t)+a 1 t)+a 3 t))]y [ȧ 3 t)+a 1 t)+a t))]z +bt) ] 45) 48) = ȧ 1 t)+[a t)+ x+z)] +[a 3 t)+x y) )] + [ x y z ȧt) a t)+a 3 t)) ] 50) = ȧ 1 t)+a t)+a 3 t)) x+ z+ y+ x y z ȧt) a t)+a 3 t)) 51) = 0. 5) Similarly by the nice symmetry of the functions 44), we could balance the second momentum equation 41), : t u +u 1 x u +u 3 z u + ) P 53) y ρ = ȧ t)+[a 1 t)+y z)] x [a t)+ x+z)]+[a 3 t)+x y)] z [a t)+ x+z)] + [ [ x +y +z xy +yz +xz) ] ] [ȧ 1 t)+a t)+a 3 t))]x 54) y [ȧ t)+a 1 t)+a 3 t))]y [ȧ 3 t)+a 1 t)+a t))]z +bt) = ȧ t)+[a 1 t)+y z)] )+[a 3 t)+x y)] 55) + y x z ȧ t) a 1 t)+a 3 t))

7 Exact Solutions to the 3-dim Euler Equations 7 = 0. 56) For the last equation 41),3, we have t u 3 +u 1 x u 3 +u y u 3 + ) P z ρ 57) = ȧ 3 t)+[a 1 t)+y z)] x [a 3t)+x y)]+[a t)+ x+z)] y [a 3t)+x y)] + [ [ x +y +z xy +yz +xz) ] ] [ȧ 1 t)+a t)+a 3 t))]x 58) z [ȧ t)+a 1 t)+a 3 t))]y [ȧ 3 t)+a 1 t)+a t))]z +bt) = ȧ 3 t)+[a 1 t)+y z)]+[a t)+ x+z)] )+ z y x ȧ 3 t) a 1 t)+a t)) 59) = 0. 60) Next, for the checking of the solutions 45), we have for the first momentum equation 41),1 : = t u 1 +u y u 1 +u 3 z u 1 + x ) P ρ 61) with the pressure function P ρ = [ x +y +z +xy +yz +xz ) [ȧ 1 t)+a t)+a 3 t))]x [ȧ t)+a 1 t)+a 3 t))]y [ȧ 3 t)+a 1 t)+a t))]z +bt) ] 6) to have = ȧ 1 t)+[a t)+x+z)] y [a 1t)+y +z)]+[a 3 t)+x+y)] z [a 1t)+y +z)] + [ x +y +z +xy +yz +xz ) ] [ȧ 1 t)+a t)+a 3 t))]x 63) x [ȧ t)+a 1 t)+a 3 t))]y [ȧ 3 t)+a 1 t)+a t))]z +bt) = ȧ 1 t)+[a t)+x+z)]+[a 3 t)+x+y)] x y z ȧt) a t)+a 3 t)) 64) = ȧ 1 t)+a t)+a 3 t))+ x+ z+ y x y z ȧ 1 t) a t)+a 3 t)) 65) = 0. 66) Similarly, we can balance the second momentum equation 41), with the symmetric form of the functions 45): = t u +u 1 x u +u 3 z u + ) P 67) y ρ = ȧ t)+[a 1 t)+y +z)] x [a t)+x+z)]+[a 3 t)+x+y)] z [a t)+x+z)]+ P y ρ 68) = ȧ t)+[a 1 t)+y +z)]+[a 3 t)+x+y)] y x z ȧ t) a 1 t)+a 3 t)) For the last equation 41),3, we could get = ȧ 3 t)+[a 1 t)+y +z)] = t u 3 +u 1 x u 3 +u y u 3 + z = 0. 69) ) P ρ 70) x [a 3t)+x+y)]+[a t)+x+z)] y [a 3t)+x+y)]+ ) P z ρ )

8 8 Manwai Yuen = ȧ 3 t)+[a 1 t)+y +z)]+[a t)+x+z)] z y x ȧ 3 t) a 1 t)+a t)) 71) = 0. 7) It is clear to see that 1. if a i T) = or ȧ i T) = for i = 1, or 3,with the first finite positive constant T, the solutions 44) and 45) blow up in the finite time T;. for all global 1 functions a i t), the solutions 44) and 45) globally exist. The proof is completed. Here the kinetic energy of the solutions 44) and 45) for the incompressible fluids are infinite: 1 u = +. 73) R 3 And the blowup solutions 44) and 45) are the other examples with infinite energy for the incompressible Euler equations 41), with respect to the blowup cylindrical ones in Gibbon, Moore and Stuart s work [4]. Additionally it is because u = 0 74) in the solutions 8) and 11) to be the corresponding solutions for the 3-dimensional compressible Navier-Stokes equations: ρ t + ρ u) = 0 75) ρ[ u t + u ) u]+ Pρ) = µ u with µ > 0. And the solutions 44) and 45) are the corresponding ones for the 3-dimensional incompressible Navier-Stokes equations: div u = 0 76) ρ[ u t + u ) u]+ Pρ) = µ u with µ > 0. 4 Discussion 1. We observe that the constructed density functions ρ 8) and 11), and the pressure functions 44) and 45) share some nice algebraic geometry, but we do not understand for the implication of this geometric structure for the fluids. ould the solutions with small perturbations of these solutions converge to the original ones?. ould we just cut a finite partial volume of the solutions to construct new week solutions? If we could, what happens for the corresponding weak solutions of the classical ones with infinite energy? 3. By taking consideration with Makino s blowup solutions form 4), we could guess that the more general functional form for the 3-dimensional velocities u in artesian coordinate: u 1 = a 1 t)+b 11 t)x+b 1 t)y +b 13 t)z u = a t)+b 1 t)x+b t)y +b 3 t)z 77) u 3 = a 3 t)+b 31 t)x+b 3 t)y +b 33 t)z with the local smooth time 1 functions a i t) and b ij t) for i, j = 1, and 3, could be the more general solutions structure for the systems. Then, it is more difficult to handle thelinearfunctionalstructureforthehigherdimensionalcasesr N withn 4. ouldweconstruct a blowup example with rotation for the incompressible fluids? In principle, we could adopt the guessing approach to construct the density functions. Some trials with luck are needed to obtain a right novel functional form to construct solutions for the systems, but it costs a lot of time for computation. Thus, we hope to arouse the readers interests to investigate more systematic

9 Exact Solutions to the 3-dim Euler Equations 9 approaches to calculate the possible particular solutions with the linear functional velocities u 77) in artesian coordinate or other coordinates. Finally, further researches are needed for understanding this class 77) of blowup or global solutions for the systems in the future. References [1] V. I. Arnold, Sur la Topologie des Ecoulements Stationnaires des Fluides Parfaits, French). R. Acad. Sci. Paris ), [] G.Q. hen and D.H. Wang, The auchy Problem for the Euler Equations for ompressible Fluids, Handbook of Mathematical Fluid Dynamics, Vol. I, , North-Holland, Amsterdam, 00. [3] P. onstantin, The Euler Equations and Nonlocal onservative Riccati Equations, Internat. Math. Res. Notices 000, [4] J. D. Gibbon, D. R. Moore and J. T. Stuart, Exact, Infinite Energy, Blow-up Solutions of the Three-dimensional Euler Equations, Nonlinearity ), [5] P.L. Lions, Mathematical Topics in Fluid Mechanics, Vols. 1,, 1998, Oxford: larendon Press. [6] T. Makino, Exact Solutions for the ompressible Eluer equation, Journal of Osaka Sangyo University Natural sciences ), [7] K. Ohkitani and J. D. Gibbon, Numerical Study of Singularity Formation in a lass of Euler and Navier-Stokes Flows, Phys. Fluids 1 000), [8] L. V. Ovsiannikov, Dokl. Akad. Nauk SSSR ), 47. [9] L. I. Sedov, Dokl. Akad. Nauk SSSR ), 753. [10] T. Senba and T. Suzuki, Applied Analysis, Mathematical Methods In Natural Science, Imperial ollege Press, London, 004. xiv+378 pp. [11] L.H. Yeung and M.W. Yuen, Note for Some Exact Blowup Solutions to the Pressureless Euler Equations in R N [ommun. Nonlinear Sci. Numer. Simul ), ], ommun. in Nonlinear Sci. and Numer. Simul., In Press, DOI: /j.cnsns [1] M.W. Yuen, Blowup Solutions for a lass of Fluid Dynamical Equations in R N, J. Math. Anal. Appl ), [13] M.W. Yuen, Analyitcal Solutions to the Navier-Stokes Equations, J. Math. Phys., ), 11310, 10pp. [14] M.W. Yuen, Some Exact Solutions to the Pressureless Euler Equations in R N, ommun. Nonlinear Sci. Numer. Simul ), [15] M.W. Yuen, Perturbational Blowup Solutions to the 1-dimensional ompressible Euler Equations, Pre-print, arxiv: [16] M.W. Yuen, Self-Similar Solutions with Elliptic Symmetry for the ompressible Euler and Navier-Stokes Equations in R N, Pre-print, arxiv: