arxiv: v1 [math.ap] 12 Apr 2016

Transcription

1 arxiv: v [math.ap] Apr 06 Confinement of a hot temperature patch in the modified SQG model. Roberto Garra, Department of Statistical Sciences, Sapienza University of Rome, Piazzale Aldo Moro 5, 0085 Rome, Italy. Abstract In this paper we study the time evolution of a temperature patch in R according to the modified Surface Quasi-Geostrophic Equation (SQG) patch equation. In particular we give a temporal estimate on the growth of the support, providing a rigorous proof of the confinement of a hot patch of temperature in absence of external forcing, under the quasi-geostrophic approximation. Keywords: Modified Surface Quasi-Geostrophic equation, Confinement of temperature patch, Active scalar flow MSC 00: 35Q86, 76E0 Introduction In the present paper we study the time evolution of a patch of temperature in R according to the modified surface quasi-geostrophic (SQG) equation, considered in recent years by different authors (see e.g. [] and the references therein). We recall the basic equations of the α-patch model: let θ(x,t), x R be the solution of the equation t θ+u θ = 0, (.) where ( ψ u = (u,u ) =, ψ ), x x (.) θ(x,t) = Λ α ψ(x,t), (.3) and Λ =. We denote this system as active scalar flow, where the scalar field θ(x,t) plays the role of a temperature field. In this paper we consider the evolution of a positive address: roberto.garra@uniroma.it

2 compactly supported initial temperature θ 0 L (R ) according to (.) for α (0,). Here we assume that the support of θ 0 is contained in the ball centred at the origin of bounded diameter R 0. The main aim of this paper is to give an estimate on the growth rate of the diameter of the support of a hot patch of temperature evolving according to (.). From the physical point of view the problem here considered is the confinement in time of a hot patch of temperature in absence of an external forcing, under the quasi-geostrophic approximation. We briefly recall the main results obtained about the growth rate of the diameter of the support of a vortex patch in the case of the Euler equation in D (that is (.3) for α = 0). The first result, obtained by Marchioro in [6], gives a bound of O(t /3 ) on the diameter of the support. An improvement of this result was given by Iftimie et al. in [9], where a bound of O ( (tlnt) /4) was proved. For the sake of completness we also recall that a similar estimate was independently obtained at the same time by Serfati in [0]. In the paper by Iftimie et al. [9], the main result was obtained in two different ways. Indeed a first proof of Proposition. in [9] was based on the application of sharp estimates on the mass of vorticity near the boundary by using a mollifier (similarly to the proof given in [6]). The second proof in the Appendix of the same paper, was given by Gamblin and is based on estimates of higher order momenta m n (t) = x 4n θ(x,t)dx. (.4) R In this paper we will generalize the last result, for the estimate of the growth of the support of a temperature patch according to the α-patch model. We will show that the proof is based on non-trivial adaptation of the tecniques used in [9], due to the different Green function of the problem. The physical motivation: confinement of an hot patch of temperature in the SQG model Here we give an outline of the physical problem considered in this paper, firstly recalling the derivation of the model equations and their meaning. The behavior of the atmosphere in the mid-latitudes is described by the equations of rotating fluids and in many cases geophysicists adopt the so-called quasi-geostrophic approximation, that is a perturbative expansion around a uniformly rotating state (see [7]). Under this approximation, the long-scale dynamics of the atmosphere in mid-latitudes is governed by the balance between the Coriolis force and the pressure gradient, i.e. fu = p, x fu = p (.), x where we denoted with u (u,u ) the velocity field, p(x,y,t) the pressure field and f is the Coriolis parameter (assumed constant). We observe that, by means of (.), we can consider

3 the pressure field as the stream function in this kind of problems in fluid mechanics. In order to obtain the evolution equations for the pressure field, a perturbative expansion with respect to the zeroth-order geostrophic field should be done. At the first order of this expansion, the equation of conservation of the potential vorticity arise, that is ( ) t +(u ) q = 0, (.) where q = ψ x + ψ x ( ( + ) ) f ψ. (.3) x 3 N x 3 In equation (.3), q( ) is the potential vorticity and N the buoyancy frequency. The further assumption in the mathematical derivation of the surface quasi-geostrophic equation (SQG) is to take (f/n) = and q = const. (for simplicity null), therefore obtaining from (.3) ( x + x + x 3 ) ψ = 0. (.4) Is moreover usual in geophysics to apply the Boussinesq approximation, for which the vertical derivative of the pressure field is proportional to the temperature field θ, that is θ p x 3 = ψ x 3. (.5) Moreover, from the energy balance law, we have that also the temperature field is advected by the fluid ( t +(u ))θ = 0. (.6) By using equations (.4) and (.5) we finally have where θ = x 3 ψ = ( ) / ψ, (.7) = x + x. We recall that the operator appearing in (.7) is a pseudo-differential operator, defined by the Fourier symbol as ( ˆ ) / ψ(k) = k ˆψ(k) (.8) and it is the operator mapping Neumann to Dirichlet boundary value problems (for a discussion about this point see e.g. []). Assuming that the fluid is incompressible, i.e. that the velocity field is divergence free, we finally obtain the basic equations of the SQG model { t θ+(u )θ = 0, (.9) u = 0, 3

4 where θ = ( ) / ψ. Starting from the seminal paper of Constantin, Majda and Tabak [5] about the singular front formation in a model for a quasi-geostrophic flow, a great number of works about the SQG equation and its generalizations with dissipation appeared in literature (see e.g. [3],[4], [3],[4],[9] and the references therein). The modified SQG equation was recently studied in different papers as a family of active scalar flows interpolating the D Euler and SQG equations. Indeed, we observe that, equations (.)-(.)-(.3), for α = 0 describe the flow of an ideal incompressible fluid according to the Euler equation in two dimensions in which θ has the meaning of vorticity, while for α = the Surface Quasi-Geostrophic equation (SQG) is recovered in which θ has the meaning of temperature. This last system has relevance in geophysics as just shown (see for instance [7]). The modified SQG equation was firstly studied by Pierrhaumbert et.al [8], in relation to the spectrum of turbulence in -d, beginning the research on α turbulence (see also [8]). On the other hand, from a mathematical point of view, this dynamical system was deeply studied (see [5],[4]) in relation to the problem of formation of singular fronts. The authors observed a formal analogy with the 3D Euler equation and studied the singular behavior of the solutions in the more suitable framework of the D quasi-geostrophic equation. The patch problem for SQG-type equations has been considered by different authors in recent papers and leads to nontrivial problems. Local existence results were proved by Gancedo in [7] for the SQG and modified SQG patches with boundaries which are H 3 closed curves. Computational studies have been also performed in [6] and [5], suggesting a finite time singularity with two patches touching each other. Finally, in the recent paper [], Kiselev et al. have studied the modified SQG patch dynamics in domains with boundaries. Some results in this direction have been also proved by Cavallaro et al. in [] where the dynamics of N concentrated temperature fields has been considered. As far as we know, a detailed analysis of the growth in time of the radius of a single patch of temperature evolving according to the modified SQG equation in the whole plane is still missing in the literature, this will be the main object of our paper. We also briefly add some considerations on the modified SQG that will be useful in the next sections. It is possible to introduce a weak form of eq.s (.)-(.4): d dt θ[f] = θ[u f]+θ[ tf], (.0) where f(x,t) is a bounded smooth function and θ[f] = dx θ(x,t) f(x,t). (.) In (.0) the velocity field u is given by u(x,t) = K(x y)θ(y, t)dy, (.) where K(x y) = G(x y), being = ( x, x ) and G(x) the Green function of ( ) ( α/) with vanishing boundary condition at infinity. 4

5 Moreover the divergence-free of the velocity field and the fact that θ is transported by the flow show that the Lebesgue measure and the maximum of θ are conserved duringthe motion. We remember that the Green function of ( ) ( α/) in R with vanishing boundary conditions at infinity is (see for instance []): G(r) = ψ(α)r α, ψ(α) = Γ(α/) π ( α) Γ( α ), r = x +x, (.3) where Γ( ) denotes the Euler Gamma function. Finally, we recall, for the utility of the reader, that several quantities are conserved in the time evolution for the α-patch model (see e.g. []), that is the total mass θ(x,t)dx = θ 0 (x)dx = m 0. the maximum norm θ(x,t) L = θ 0 (x) L = M 0. (.4) the center of mass m 0 x θ(x,t)dx = m 0 x θ 0 (x)dx = c 0. the moment of inertia x θ(x,t)dx = x θ 0 (x)dx = i 0. These conservation laws play a central role for the proof of the main theorem in the next section. Theaim ofthis paperistorigorouslystudytheconfinementof anhotpatch oftemperature in absence of an external forcing. We are therefore interested in proving some rigorous results on the concentration of a temperature field under the quasi-geostrophic approximation, in the more flexible formulation of the modified SQG equation in which the parameter α (that appears in the fractional Laplacian of equation (.3)) describe the difference in the behavior of the air mass between the description given by the ideal incompressible fluid (according to the Euler equation) and the one provided by the SQG equation. In particular in the next section we generalize some rigorous results previously proved by Iftimie et al. in [9] for the Euler equation in the plane, in order to estimate the growth in time of a hot patch of temperature according to the α-patch model. 5

6 3 Estimates on the growth of the support of a temperature patch We now give an estimate on the growth of the support of a temperature patch, assuming that the support of the initial temperature field θ 0 is contained in the ball centered at the origin of radius R 0. Theorem 3. Let θ(x,t) be a solution of the active scalar flow equation with a positive compactly supported initial temperature θ 0 L (R ). There exists a constant C 0 C 0 (i 0,R 0,m 0,M 0 ) such that, for every t 0 the support of θ(x,t) is contained in the ball x < 4R 0 + C 0 [tln(+t)] 4+α, for any α (0,). Proof The proof of the theorem is strongly based on the use of the conserved quantities and in particular of the moment of inertia i 0 and the total mass m 0. The main difference and technical difficulty with respect to the case α = 0 is due to the different Green function given by (.3). Without loss of generality, we will take the center of mass located at the origin. We should prove that the radial component of the velocity field satisfies the following estimate x x u(x,t) C 0 x3+α, (3.) for all x 4d 0 +C 0 [tln(+t)] 4+α. The radial component of the velocity field in the case of the active scalar flow is given by x x u = x x x G(x,y)θ(y)dy = x (x y) x y +αθ(y)dy, (3.) where (x,x ) = ( x,x ). Let us divide the last integral in two regions: A := x y < x A := x y x. The integral over the region A is simply bounded by C A θ(y) x y α+dy. (3.3) By using the fact that x (x y) = x y and the conservation of the center of the mass (that is located at the origin), we have x A x (x y) x y A α+θ(y)dy = x y x x y +αθ(y)dy x y = A ( x x y +α ) x y x +α θ(y)dy + A x 3+αθ(y)dy (3.4) 6

7 The first integral can be bounded as follows x y A ( x x y α+ ) x α+ θ(y) dy const i 0, (3.5) x 3+α where we have used the conservation of the moment of inertia. Regarding the integral over the region A, being x y < x 3 x, then we have that y and therefore x y θ(y)dy A x 3+α const. θ(y) x y α+dy. (3.6) At this stage, we have obtained the following estimate for the radial component of the velocity field x x u(x) const. i 0 θ(y) +const. x 3+α x y α+dy (3.7) We now introduce the following Lemma Lemma 3. Let S R and h : S R + a function belonging to L (S) L (S). Then β h(y) x y βdy C h(y) L h(y) β/ L, (3.8) with β (0,). S The proof is essentially the same of Lemma. in [9], details of the proof are given in the Appendix. By using the fact that {y : x y < x /} {y : y > x /} and the Lemma 3. we have the following estimate on the radial component of the velocity ( x x u C α+ +C M x 3+α 0 A A θ(y)dy y > x / ) α. (3.9) In the next proposition we prove that also the last integral is of order O( x 3 α ) for x large, so that inequality (3.) holds and the proof of the theorem is achieved. Proposition 3. There exists a constant C 0 C 0 (i 0,R 0,m 0,M 0 ) such that, for any k > 0 θ(y,t)dy C 0 x k, (3.0) y > x / for all x > 4R 0 +C 0 (t ln(+t)) 4+α. Proof. The proof of this proposition is based on the following technical Lemma: 7

8 Lemma 3.3 Let be m n,α (t) = x (4+α)n θ(x,t)dx, (3.) then, there exists a constant C 0 such that, for any n m n,α (t) m 0 ( R 4+α 0 +C 0 i 0 n +α t ) n. (3.) We postpone the proof of this Lemma and use it to prove Proposition 3.. Let us suppose that r 4+α [R 4+α 0 +C 0 i 0 ktln(+t)], (3.3) taking n such that then we have from Lemma 4. that x r kln(+t) < n α+ kln(+t), (3.4) θ(x,t)dx m n,α(t) r 4(n+α) m 0 r k when r 4+α [R 4+α 0 +C 0 i 0 ktln(+t)], as claimed. ( R 4+α 0 +C 0 i 0 n α+ t ) n r 4n k+nα const. m 0 r k, Proof of Lemma 4.. By using (.) and (.) and integrating by parts, we have that m n,α (t) = x (4+α)n t θ(x,t)dx = x (4+α)n u r θ(x,t)dx (3.5) = (4+α)n x (x y) π x y +α x 4n+αn θ(x,t)θ(y,t)dxdy. We define so that, we can write in a compact way m n,α(t) = (4+α)n π We consider the following partition A = {(x,y) : y A = {(x,y) : K(x,y) = x (x y), (3.6) x y +α ( n A 3 = {(x,y) : x K(x,y) x 4n+αn θ(x,t)θ(y,t)dxdy. (3.7) ( n ) x } ) x y ) y }. ( n ( ) x } n 8

9 Then, we have that d dt m n,α(t) = 3 i= (4+α)n π A i K(x,y) x 4n+αn θ(x,t)θ(y,t)dxdy. (3.8) We first consider the region A, where x y x /n and we have the following bound such that K(x,y) x nα+ const. x y α+ x α, (3.9) (4+α)n K(x,y) x 4n+αn θ(x,t)θ(y,t)dxdy const. n +α x (4+α)n α θ(x,t)θ(y,t)dxdy π A A n 4+α const. ( n) x (4+α)(n ) y θ(x,t)θ(y,t)dxdy const. n +α i 0 m n,α (t), A (3.0) where we used also the fact that y / x ( n ) in A and the conservation of the moment of inertia. In the set A 3, we have that x y y /n and we can obtain the following estimate ( ) x n α+ y K(x,y) const. nα+ x y α+ y α const. n x α+, (3.) where we also used the fact that in A 3, ( n ) y x. It follows that, also in this case, (4+α)n K(x,y) x 4n+αn θ(x,t)θ(y,t)dxdy const.n +α i 0 m n,α (t). (3.) π A Finally, the contribution related to the set A can be bounded in a similar way by a term proportional to m n,α (t). Thus we have the following estimate We now rearrange the last term x (4+α)(n ) θ(x,t)dx = and applying the Holder s inequality, we finally obtain and integrating we obtain the claimed result. d dt m n,α(t) C 0 i 0 n +α m n,α (t). (3.3) x (4+α)(n ) θ /n (x,t)θ /n (x,t)dx, (3.4) d dt m n,α(t) C 0 i 0 n +α m /n 0 m /n n,α (t), (3.5) 9

10 Remark 3. As a final remark, we observe that uniqueness result for the global solution of the modified SQG equation with an initially positive compactly supported datum is still missing and therefore we are forced to state our main result in a weak form, that is a rigorous estimate of the growth of the support of a solution of the governing equation. Acknowledgments We wish to thank Carlo Marchioro (Rome) for fruitful discussions. 4 Appendix: Proof of Lemma 3. For the sake of completness, we here provide the detailed proof of Lemma 3.. In particular, we here recall a more general Lemma (see [0]) that includes Lemma 3. as a special case Lemma 4. Let β (0,), S R and h : S R + a function belonging to L (S) L p (S), with p > β. Then with β > 0. Proof S β /p β h(y) x y βdy C h(y) /p /p L (S) h(y) L p (S), (4.) Let k > 0, by using the Hölder inequality, we have that h(y) S x y β = h(y) S { x y >k} x y βdy + h L (S) k β + h L p (S) x β p L = h L (S) k β +C h L p (S)k β /p. p ( x k) S { x y <k} h(y) x y βdy ( ) and by taking k = h L (S) h /p L p (S), we conclude the proof. Finally observe that the Lemma 3. is a particular case, for h L (S) L (S). References [] Caffarelli, L., Silvestre, L.: An extension problem related to the fractional Laplacian. Communications in partial differential equations, 3.8: 45 60, (007). [] Cavallaro, G., Garra, R. and Marchioro, C.: Localization and stability of active scalar flows, Riv. Mat. Univ. Parma, 4, 75 96, (03) 0

11 [3] Chae, D., Constantin, P., Wu, J.: Inviscid models generalizing the two-dimensional Euler and the Surface Quasi-geostrophic equations, Arch. Rational Mech. Anal., 0, 35 6, (0) [4] Constantin, P., Cordoba, D., Wu, J.: On the Critical Dissipative Quasi-geostrophic Equation, Indiana University mathematics journal, 50, 97 06, (00) [5] Constantin, P., Majda, A.J., Tabak, E.: Formation of strong fronts in the -d quasigeostrophic thermal active scalar, Nonlinearity, 7, , (994) [6] Cordoba, D., Fontelos, M. A., Mancho, A.M., Rodrigo, J. L.: Evidence of singularities for a family of contour dynamics equations. Proceedings of the National Academy of Sciences, 0(7), , (005) [7] Gancedo, F.: Existence for the alpha-patch model and the QG sharp front in Sobolev spaces, Advances in Mathematics,7, , (008) [8] Held, I.M., Pierrehumbert, R.T., Garner, S.T., Swanson, K.L.: Surface quasi-geostrophic dynamics, J. Fluid Mech., 8, 0, (995) [9] Iftimie, D., Sideris, T.C., Gamblin, P.: On the evolution of compactly supported planar vorticity, Comm. Partial Differential Equations, 4(9-0), ,(999), [0] Iftimie, D.: Evolution de tourbillon à support compact, Actes du Colloque de Saint-Jeande-Monts, (999) [] Iwayama, T., Watanabe, T.: Green s function for a generalized two-dimensional fluid, Physical Review E, 8, , (00) [] Kiselev, A., Ryzhik, L., Yao, Y., Zlatos, A.: Finite time singularity formation for the modified SQG patch equation, arxiv: , (05) [3] Li, D.: Existence theorems for the D quasi-geostrophic equation with plane wave initial conditions, Nonlinearity,, , (009) [4] Majda, A.J., Tabak, E.: A two-dimensional model for quasi-geostrophic flow: comparison with the two-dimensional Euler flow, Physica D, 98, 55 5, (996) [5] Mancho, A.: Numerical studies on the self-similar collapse of the alpha-patches problem, Communications in Nonlinear Science and Numerical Simulation, 6., 5 66, (05) [6] Marchioro, C.: Bounds on the growth of the support of a vortex patch, Comm. Math. Phys., 64, , (994). [7] Pedlosky, J.: Geophysical Fluid Dynamics, Springer-Verlag, (987) [8] Pierrehumbert, R.T., Held, I.M., Swanson, K.L.: Spectra of local and nonlocal twodimensional turbulence, Chaos, Solitons and Fractals, 4, 6, (994)

12 [9] Rodrigo, J.L.: On the evolution of sharp Fronts for the Quasi-Geostrophic equation, Commun. Pure and Appl. Math., 58, 8 866, (005) [0] Serfati, P.: Borne en temps des caractéristiques de l équation d Euler D á tourbillon positif et localisation pour le modéle point-vortex, Preprint, (998)