arxiv: v1 [math.ap] 8 Jan 2018

Transcription

1 GLOBAL STRONG SOLUTIONS FOR THE THREE-DIMENSIONAL HASEGAWA-MIMA MODEL WITH PARTIAL DISSIPATION arxiv:181.37v1 [math.ap] 8 Jan 18 CHONGSHENG CAO, YANQIU GUO, AND EDRISS S. TITI Abstract. We study the three-dimensional Hasegawa-Mima model of turbulent magnetized plasma with horizontal viscous terms and a weak vertical dissipative term. In particular, we establish the global existence and uniqueness of strong solutions for this model. 1. Introduction 1.1. Literature. In 1977, Hasegawa and Mima introduced a system in [1, 11] to elucidate the drift wave turbulence in Tokamak, the most advanced magnetic confinement device. The three-dimensional inviscid Hasegawa-Mima equations can be written as cf. [, 3, 1, 11, 17, ]: w φ +Jφ,w+ =, t z 1.1 t hφ φ+jφ, h φ+γ φ y w =, z 1. where Jf,g = f g f g is the Jacobian and x y y x h = + is the horizontal x y Laplacian. System describes the coupling of the drift modes to the ionacoustic waves that propagate along the magnetic field. Here, φ is the electrostatic potential, and simultaneously is the stream function for the horizontal flow in the xy-plane. Moreover, w represents the normalized ion velocity in the z-direction, and γ is a constant which is proportional to the density gradient. Like the three-dimensional Euler equations of inviscid incompressible fluid, the only conserved quantity for the 3D Hasegawa-Mima equations is the kinetic energy, and the global regularity problem is open. Nevertheless, by adding the full viscosity to , Zhang and Guo [] proved the global regularity and the existence of global attractors for a viscous and forced 3D Hasegawa-Mima model using standard tools from the theory of Navier-Stokes equations. On the other hand, Cao, Farhat and Titi [3] proposed and studied an inviscid three-dimensional modified Date: January 6, 18. Key words and phrases. Hasegawa-Mima equations, plasma turbulence, global strong solutions, partial dissipation, fluid dynamics. 1

2 C. CAO, Y. GUO, AND E. S. TITI version of , the pseudo-hasegawa-mima equations: w t +u hw U L ω =, 1.3 z ω t +u hω U w =, 1.4 L z with h u =, for some constant U, where u = u,v tr is the horizontal component of the velocity vector field u,v,w tr, and ω = h u is the vorticity. The operator h =, x y tr is the horizontal gradient. In particular, the global well-posedness of the weak solutions to was established in [3]. Observe that ω in plays the role of the term h φ φ in Therefore, system is a modified version of the Hasegawa-Mima equations , with the essential difference that the term φ ω is replace by. Nevertheless, model is simpler z z than in the sense that it has a nice mathematical structure. Indeed, adding and subtracting 1.3 and 1.4 yield a three-dimensional coupled transport system with collinear transport velocities in opposite directions leading to an intensified shear in the vertical direction, which results in exponential growth in the relevant estimates for in [3]. It is worth mentioning other interesting models describing plasma turbulence. For instance, Hasegawa and Wakatani proposed equations for a two-fluid model which describe the resistive drift wave turbulence in Tokamak cf. [1, 13]. The existence and uniqueness of strong solutions to the Hasegawa-Wakatani equations have been established by Kondo and Tani [14]. In the context of geophysical fluid dynamics, there are certain models resemble the structure of Hasegawa-Mima equations In particular, Charney [5] and Obukhov [18] derived the following two-dimensional shallow water model from the Euler equations with free surface under a quasi-geostrophic velocity field assumption: t hφ Fφ +Jφ, h φ +Jφ,φ B +βy =. 1.5 Here φ x,y is the amplitude of the surface perturbation at the lowest order in the Rossby number, and the equation z = φ B x,y describes the given bottom topography. F is the Froude number. One may refer to [] for a derivation of model 1.5. For the simple case when φ B is a constant representing a flat bottom, 1.5 reduces to the Hasegawa-Mima-Charney-Obukhov equation: t hφ Fφ +Jφ, h φ +β φ =. 1.6 x Since 1.6 bears a close resemblance to the two-dimensional Euler equations, the standard tools for handling the D Euler equations can be adopted to analyze 1.6. Indeed, Guo and Han [7] proved the global existence and uniqueness of solutions for 1.6. For other results concerning 1.6 see, e.g., Paumond [19], and Gao and Zhu [6].

3 THE HASEGAWA-MIMA MODEL WITH PARTIAL DISSIPATION 3 It is worth mentioning that one may refer to the monographs [16, ] as well as the papers [8, 9] for other relevant geophysical models. 1.. The model. Motivated by the Hasegawa-Mima equations and the Charney- Obukhov equations mentioned in subsection 1.1, we introduce and study in this paper the following three-dimensional Hasegawa-Mima model with horizontal viscous terms and a weak vertical dissipative term: w t +u hw ψ z = 1 Re hw, 1.7 ω t +u hω w z = 1 Re hω +ǫ ψ z, 1.8 h u =. 1.9 The velocity vector field u,v,w tr defined in = [,L] [,1] satisfies the periodic boundary condition with the horizontal velocity u = u,v tr. The stream function ψ for the horizontal flow is defined as ψ = h 1 ω with ψdxdy =, and [,L] ω = h u. We denote h =, x y tr and h = x + y. The constant Re is the Reynolds number. System bears a resemblance as the three-dimensional Hasegawa-Mima equations with the difference that Hasegawa-Mima equations are inviscid, whereas model is regularized by horizontal viscosity and a partial vertical dissipation. The purpose of introducing and investigating is to shed light on the analysis of the inviscid Hasegawa-Mima equations Mathematically, the difficulty of establishing the global regularity for system lies in the following aspects: i The physical domain is three-dimensional. ii The regularizing viscosity acts only on the horizontal variables. iii The system contains the troublesome term ψ. z Since the lack of the viscosity in the vertical direction provides great challenge for establishing the global regularity, we impose a weak dissipative term ǫ ψ in the z equation 1.8. Since ψ = h 1 ω, we remark that, as a dissipation, ψ is weaker z than the vertical viscosity ω. In a priori estimates conducted in section, the z dissipative term ǫ ψ plays a vital role in controlling the terms w and z ψ with z z the help of an anisotropic Ladyzhenskaya type inequality see Lemma Preliminaries. In this subsection, we introduce some preliminaries that will be used later in our analysis. Recall the three-dimensional periodic space domain = [,L] [,1]. Throughout, the norm for the L p space, for p [1, ], is denoted by f p. The inner product of f and g in the L space is denoted by

4 4 C. CAO, Y. GUO, AND E. S. TITI f,g = fgdxdydz. As usual, the Sobolev space H1 = {f L : f L }. In addition, we define the following Hilbert space: H 1 h = { f L : h f L }, 1.1 that features the inner product f,g H 1 h = f,g+ h f, h g. For sufficiently smooth functions f, g and u, with h u =, integration by parts yields which immediately implies that u h f,g = u h g,f, 1.11 u h f,f =. 1.1 Recall that the horizontal velocity u, the vertical vorticity ω, and the stream function ψ for the horizontal flow have the following relations: ω = h u = v x u y, ω = h ψ, u = ψ y, ψ x tr, 1.13 where [,L] ψdxdy =. It follows that, if ω L, then ω,ψ = u In addition, for sufficiently smooth functions f, u and ψ such that u = ψ y, ψ x tr, observe that u h ψ =, then apply 1.11 to deduce u h f,ψ = u h ψ,f = Main result. Before we state the main result of the paper, we give a definition of a strong solution for system Definition 1.1. We call u,w tr = u,v,w tr a strong solution on [,T] for system if i u,w tr has the following regularity: u, w L,T;H 1 C[,T];L ; h u, h w, ω z, h w z, ψ zz L,T; 1.16 u t, w t L,T; ii the equations below hold in the following sense: w t +u hw ψ z = 1 Re hw, in L,T; ω t +u hω w z = 1 Re hω +ǫ ψ z, in L,T;Hh 1, with h u =, where ω = h u, ψ = h 1 ω with ψdxdy =, and [,L] Hh 1 is the dual of the space Hh 1, defined in 1.1.

5 THE HASEGAWA-MIMA MODEL WITH PARTIAL DISSIPATION 5 Now we are ready to state the main result of the paper: the global existence, uniqueness, and continuous dependence on initial data of strong solutions for our model Theorem 1.. Let T >. Assume u,w tr H 1 3, then system admits a unique strong solution u,w tr on [,T] in the sense of Definition 1.1 satisfying the initial condition u,w tr = u,w tr. Moreover, the energy equality is valid for every t [,T]: 1 t wt + ut + = 1 w + u [ ] 1 h w + h u +ǫ ψ z ds Re In addition, the H 1 norm of the solution u,w tr has a uniform bound independent of T. That is, sup ut H 1 + wt H 1 K, t T where K is independent of T, but depends only on Re, ǫ, L, u H 1 and w H 1. Furthermore, if {u n,wn tr } is a bounded sequence of initial data in H 1 such that u n,wn tr u,w tr in L, then the corresponding strong solutions u n,w n tr and u,w tr satisfy u n,w n tr u,w tr in C[,T];L.. A priori estimates In this section, we assume that system holds for smooth functions and we establish the following formal a priori estimates. However, as we will show in section 3, these formal estimates can be justified rigorously by establishing them first for the Galerkin approximation system and then passing to the limit using the appropriate Aubin compactness theorem..1. Estimate for w + u. Taking the L inner product of the system with w,ψ tr yields 1 d w dt + u 1 + h w Re + hu +ǫ ψ z =,.1 where we have used identities 1.1, 1.14 and Integrating.1 over the interval [, t] yields t wt + ut + h w Re + hu +ǫ ψ z ds = w + u..

6 6 C. CAO, Y. GUO, AND E. S. TITI.. Estimate for ω. Taking the inner product of 1.8 with ω yields 1 d dt ω + 1 Re hω +ǫ u z = w z,ω,.3 where 1.1 and 1.14 have been used. Thanks to 1.13, we have w z,ω = w z h ψdxdydz = h w h ψ z dxdydz h w h ψ z = h w u z ǫ u z + 1 ǫ hw..4 Combining.3 and.4 implies d dt ω + Re hω +ǫ u z 1 ǫ hw..5 By integrating.5 over the interval [, t], we obtain t ωt + Re hω +ǫ u z ds ω + 1ǫ t h w ds ω + Re ω ǫ + u,.6 where the last inequality is due to...3. An anisotropic Ladyzhenskaya type inequality. We state here the following anisotropic Ladyzhenskaya type inequality which will be useful in subsequent a priori estimates. It is worth mentioning that similar inequalities can be found in [4]. However, for the sake of completeness we present the proof of this technical lemma in the appendix. Lemma.1. Let f H 1, g Hh 1 and h L. Then fgh dxdydz C f + h f 1 f + f z 1 1 g g + h g 1 h..4. Estimate for h w. Taking the inner product of 1.7 with h w yields 1 d dt hw + 1 Re hw u h w h w dxdydz+ ψ z h w C ω 1/ u + u z 1/ h w 1/ h w 3/ + ψ z h w, where we have used Lemma.1 and the Poincaré inequality since [,L] udxdy = [,L] ψ y, ψ x tr dxdy = and [,L] h wdxdy =.

7 THE HASEGAWA-MIMA MODEL WITH PARTIAL DISSIPATION 7 By employing the Young s inequality, we obtain d dt hw + 1 Re hw C ω u + u z h w +C ψ z. Thanks to the Gronwall s inequality, we have h wt + 1 Re t h w ds C h w + t ψ z ds e t C ω u + uz ds C w, h w, ω..7 The uniform bound.7 is due to estimates. and Estimate for w z + u z. We take the L inner product of with w zz, ψ zz tr. After conducting integration by parts, one has 1 d wz + u z 1 dt + h w z + ω z +ǫ ψ zz Re u z h ww z dxdydz + u z h ψ z ω dxdydz+ u h ψ z ω z dxdydz..8 Next, we estimate each term on the right-hand side of.8. By Lemma.1 with f = h w, g = u z and h = w z, and along with the Poincaré inequality, we obtain u z h ww z dxdydz C h w 1/ h w + h w z 1/ u z 1/ ω z 1/ w z 1 h w z 6Re + ω z +C h w + u z wz Also using Lemma.1 with f = ω, g = u z and h = h ψ z gives us u z h ψ z ω dxdydz C h ω 1/ ω + ω z 1/ u z 3/ ω z 1/ 1 6Re ω z +C h ω + u z +1 u z..1 In addition, due to Lemma.1 with f = u, g = h ψ z and h = ω z, one has u h ψ z ω z dxdydz C ω 1/ u + u z 1/ u z 1/ ω z 3/ 1 6Re ω z +C ω u + u z uz..11

8 8 C. CAO, Y. GUO, AND E. S. TITI Apply estimates to the inequality.8 yields d wz + u z 1 dt + h w z + ω z +ǫ ψ zz Re C h w + u z wz +1 +C h ω + u z + ω u + ω u z uz. Thanks to Gronwall s inequality, we obtain t [ ] 1 w z t + u zt + h w z Re + ω z +ǫ ψ zz ds t z w + zu +C h w + u z + ω u u z ds { t } exp C h w + u z + hω + ω u z ds C w H 1, u H 1..1 The uniform bound.1 is due to.,.6 and Rigorous justification of the a priori estimates and the existence of strong solutions This section is devoted to proving the existence of global strong solutions for the model by assuming the initial data u,w tr H 1 3. We employ the standard Galerkin method and use the analogue of the a priori estimates that were established in section. Let e j = expπi[j 1 x+j y/l+j 3 z] for j = j 1,j,j 3 tr. For m N, let P m L be a subspace of L spanned by {e j } j m. Also, for any L function f = j Z 3 α j e j, with α j = f,e j, we write P m f = j m α je j. Let us consider the Galerkin approximation for our model : t w m +P m u m h w m z ψ m = 1 Re hw m, 3.1 t ω m +P m u m h ω m z w m = 1 Re hω m +ǫ zz ψ m, 3. h u m =, 3.3 u m = P m u, w m = P m w, 3.4 whereu m,w m P m L andω m = h u m,ψ m = h 1 ω m with [,L] ψ m dxdy =. For each m 1, the Galerkin approximation corresponds to a first order system of ordinary differential equations with quadratic nonlinearity. Therefore, by the theory of ordinary differential equations, there exists some T m > such that

9 THE HASEGAWA-MIMA MODEL WITH PARTIAL DISSIPATION 9 system admits a unique solution u m,w m tr on [,T m ]. Since u m and w m have finitely many modes, they are smooth functions, and therefore all of the a priori estimates established in section are valid for the Galerkin approximate solution u m,w m tr. In particular, the H 1 norm of u m,w m tr is uniformly bounded for all time. Hence, the Galerkin approximate solution u m,w m tr exists globally in time, in particular, over [,T], for every T >. Furthermore, by the a priori estimates in section, one has the following uniform bounds for the sequence of the Galerkin approximate solutions. u m, w m are uniformly bounded in L,T;H 1 ; 3.5 h ω m, h w m, z ω m, h z w m, zz ψ m are uniformly bounded in L,T. 3.6 Therefore, there exist a subsequence, denoted also by u m, w m, ω m, ψ m, and corresponding limits, u, w, ω, and ψ, respectively, such that u m u, w m w, weakly in L,T;H 1 ; 3.7 h ω m h ω, h w m h w, weakly in L,T; 3.8 z ω m z ω, h z w m h w z, zz ψ m zz ψ, weakly in L,T. 3.9 Moreover, due to the a priori estimates in section, we find that sup t T u m t H 1 + w mt H 1 K 3.1 where K is independent of T, but depends only on parameters Re, ǫ, L as well as the H 1 -norm, u H 1 and w H 1 of the initial data. Also thanks to the weak- convergence stated in 3.7, one has u L,T;H 1 liminf m u m L,T;H 1 and w L,T;H 1 liminf m w m L,T;H 1. Therefore, we obtain from 3.1 that sup ut H 1 + wt H 1 K. t T In order to obtain the strong convergence of the approximate solutions, we shall derive uniform bounds for t w m and t u m. First, we claim that the sequence t w m is uniformly bounded in L,T. Indeed, for any function ϕ L 4/3,T;L,

10 1 C. CAO, Y. GUO, AND E. S. TITI we use Lemma.1 to estimate T u m h w m ϕ dxdydzdt C T C sup t [,T] ω m 1/ u m + z u m 1/ h w m 1/ h w m 1/ ω m 1/ u m + z u m 1/ h w m 1/ T 1/4 T h w m dt 3/4 ϕ 4/3 dt ϕ dt C u H 1, w H 1 ϕ L 4/3,T;L, 3.11 where the last inequality is due to the a priori estimates.,.6,.7 and.1. Consequently, the sequence u m h w m is uniformly bounded in L 4,T;L. 3.1 As a result, from and 3.1, we obtain from 3.1 that the sequence t w m is uniformly bounded in L,T Next, we show that t u m is uniformly bounded in L,T. Recall the Hilbert space Hh 1 = {f L : h f L } associated with the norm f = Hh 1 f + h f. For any function φ L,T;Hh 1, we apply Lemma.1 in order to estimate T u m h ω m φ dxdydzdt C T C sup t [,T] ω m 1/ u m + z u m 1/ h ω m φ 1/ T ω m 1/ u m + z u m 1/ φ + h φ 1/ dt 1/ T 1/ h ω m dt φ + h φ dt C u H 1, w H 1 φ L,T;Hh 1, 3.14 where we have used the a priori estimates.,.6 and.1. Therefore, the sequence u m h ω m is uniformly bounded in L,T;H 1 h, 3.15 where Hh 1 is the dual space of Hh 1. Consequently, according to 3.6 and 3.15, we obtain from 3. that the sequence t ω m is uniformly bounded in L,T;H 1 h, 3.16

11 and thus THE HASEGAWA-MIMA MODEL WITH PARTIAL DISSIPATION 11 t u m is uniformly bounded in L,T Then, we infer from 3.13 and 3.17 that there is a subsequence such that t w m t w, t u m t u weakly in L,T By 3.5, 3.13, 3.17, and thanks to the Aubin s compactness theorem, we have, for a subsequence, the following strong convergence holds: u m u, w m w in L,T Next, we show the convergence of the nonlinear terms in Let η be a trigonometric polynomial with continuous coefficients. For m larger than the degree of η we have T P m u m h ω m ηdxdydz T = u h ω m ηdxdydz + T u m u h ω m ηdxdydz. 3. Since h ω m h ω weakly in L,T, u m u in L,T, and h ω m is uniformly bounded in L,T, we can pass to the limit in 3.: T T P m u m h ω m ηdxdydz = u h ωηdxdydz. 3.1 lim m An analogous argument yields T P m u m h w m ηdxdydz = lim m T u h wηdxdydz. 3. Therefore, due to , 3.18, 3.1 and 3., we pass to the limit for the Galerkin approximate equations It follows that T t w+u h w z ψ 1 Re hw ηdxdydzdt =, 3.3 T T t ω +u h ω z w 1 Re hω ǫ zz ψ ηdxdydzdt =, 3.4 h uηdxdydzdt =, 3.5 for any trigonometric polynomial η with continuous coefficients. By applying Lemma.1 as the arguments in 3.11, we can deduce that u h w L 4,T;L. Then, since t w, zz ψ and h w L,T, one has t w+u h w z ψ 1 Re hw L,T. 3.6

12 1 C. CAO, Y. GUO, AND E. S. TITI Also, using Lemma.1 as the estimates in 3.14, one may derive that u h ω L,T;H 1 h. Since z w, h ω, zz ψ L,T and t ω L,T;H 1 h, we obtain t ω +u h ω z w 1 Re hω ǫ zz ψ L,T;H 1 h. 3.7 On account of 3.6 and 3.7, we obtain from that t w+u h w z ψ = 1 Re hw, in L,T; 3.8 t ω +u h ω z w = 1 Re hω +ǫ zz ψ, in L,T;H 1 h, 3.9 with h u =. It follows from that u,w C,T;L. 3.3 Dueto3.3and3.19, onehas, foreveryt [,T], u m t utandw m t wt in L. In particular, u m u and w m w. On the other hand, by 3.4, we find that u m u and w m w. As a result, u,w tr satisfies the desired initial condition: u = u and w = w. Finally, due to the regularity of solutions, we can multiply by w,ψ tr and integrate the result over,t for t [,T]. Then the energy identity 1.17 follows. 4. Uniqueness of strong solutions This section is devoted to proving that strong solutions for the system are unique and depend continuously on the initial data. Assume there are two strong solutions u 1,w 1 tr and u,w tr on [,T] in the sense of Definition 1.1. Set u = u 1 u and w = w 1 w. Therefore, w t +u hw 1 +u h w ψ z = 1 Re hw, in L,T; 4.1 ω t +u hω 1 +u h ω w z = 1 Re hω +ǫ ψ z, in L,T;Hh 1, 4. with h u =. Since u and w satisfy the regularity 1.16, we can multiply by w,ψ tr and integrate over. By using 1.11, 1.1, 1.14 and 1.15, we obtain, for a.e. t [,T], 1 d w dt + u 1 + h w Re + hu +ǫ ψ z u h ww 1 dxdydz + u h ψω dxdydz. 4.3

13 THE HASEGAWA-MIMA MODEL WITH PARTIAL DISSIPATION 13 Next we estimate the two integrals on the right-hand side of 4.3. Using Lemma 1 with f = w 1, g = u and h = h w, we obtain u h ww 1 dxdydz C w 1 + h w 1 1/ w 1 + z w 1 1/ u 1/ h u 1/ h w h w + h u 1 4Re +C w1 + h w 1 w1 + z w 1 u. 4.4 Also, using Lemma 1 with f = u, g = h ψ, h = ω, we have u h ψω dxdydz C ω 1/ u + z u 1/ u 1/ h u 3/ 1 4Re hu +C ω u + z u u. 4.5 Now, we combine the estimates to deduce, for a.e. t [,T], d w dt + u 1 + h w Re + h u +ǫ ψ z [ w1 C + hw 1 w1 + ] zw 1 + ω u + zu u. By Gronwall s inequality, it follows that wt + ut w + u e C t w 1 + hw 1 w 1 + zw 1 + ω u + zu ds w + u e tc w 1 H 1, u H 1, 4.6 for any t [,T]. In particular, if u,w tr =, i.e., the initial values of the two solutions u 1,w 1 tr and u,w tr coincide, then 4.6 implies wt + ut = for all t [,T]. This completes the proof for the uniqueness of strong solutions. To see the continuous dependence on the initial data, we let ũ, w tr H 1 3 andtakeaboundedsequence{u n,wn tr }ofinitialdatainh 1 suchthatu n,wn tr ũ, w tr in L, and u n H 1, wn H 1, ũ H 1, w H 1 M for some M >. Denotethecorresponding strongsolutionsby u n,w n tr andũ, w tr, respectively. Then, on account of 4.6, we have, for all t [,T], w w n + ũ un w w n + ũ u n e tc w H 1, u n H 1 w w n + ũ u n e T CM. It follows that u n,w n tr ũ, w tr in C[,T];L. This completes the proof for the continuous dependence on the initial data with respect to the L -norm for the strong solutions.

14 14 C. CAO, Y. GUO, AND E. S. TITI 5. Appendix We prove the anisotropic Ladyzhenskaya type inequality stated in Lemma.1. Proof. It suffices to prove the inequality in Lemma.1 for smooth periodic functions, and then pass to the limit using a standard density argument. Recall = [,L] [,1]. For a fixed x,y [,L] and for every z, σ [,1], we have z f 4 d x,y,z = f 4 x,y,ξ dξ +f 4 x,y,σ dξ = 4 4 σ z σ 1 f 3 x,y,ξf ξ x,y,ξdξ +f 4 x,y,σ Integrating 5.1 with respect to σ over [,1], we obtain f 4 x,y,z 4 1 fx,y,ξ 3 f ξ x,y,ξ dξ+f 4 x,y,σ. 5.1 fx,y,ξ 3 f ξ x,y,ξ dξ + and by Cauchy-Schwarz inequality, we have Here we denote Now, if we denote f L z written as 1 f 4 x,y,σdσ, f 4 x,y,z 4 f 3 L 6 z f z L z + f 4 L 4 z /p f L p z = fx,y,z dz p. = sup z [,1] fx,y,z, then the inequality 5. can be f L z C f 3/4 L 6 z f z 1/4 L z + f L 4 z. 5.3 Thanks to the Hölder s inequality and 5.3, we have f gh dxdydz f L z g L z h L z dxdy [,L] C f 3/4 L f z 1/4 6 [,L] z L + f L 4 g z z L z h L z dxdy [ ] 1/4 1/4 C f 3 L f z 6 [,L] z L z dxdy + f 4 g 4 L dxdy h [,L] z C f 3/4 6 f z 1/4 + f 4 [,L] g 4 L z dxdy 1/4 h. 5.4

15 THE HASEGAWA-MIMA MODEL WITH PARTIAL DISSIPATION 15 Recall the Ladyzhenskaya inequality see, e.g., [15] in the three-dimensional periodic domain : 6 p p ϕ p C p ϕ ϕ + ϕ x p p ϕ + ϕ y p p ϕ + ϕ z p p, 5.5 for ϕ H 1 and p [,6]. By 5.5, one has and f 6 C f + h f /3 f + f z 1/3, 5.6 f 4 C f 1/4 f + h f 1/ f + f z 1/ By virtue of 5.6 and 5.7, we have f 3/4 6 f z 1/4 + f 4 C f + h f 1/ f + f z 1/. 5.8 Recall the Agmon s inequality see, e.g., [1] in one dimension: φ L [,L] C φ 1/ L [,L] φ 1/ H 1 [,L]. 5.9 By using 5.9, we obtain 1 1 g 4 L dxdy = g dz g dz dxdy [,L] z [,L] [ L 1 L 1 C g L ] dx g +gx 1 dx dzdy [ L 1 L 1 g L ] dy g +gy 1 dy dzdx C g g + h g. 5.1 By combining 5.4, 5.8 and 5.1, we conclude the proof of Lemma.1. Acknowledgment. The work of E.S.T. was supported in part by the ONR grant N References 1. S. Agmon, Lectures on elliptic boundary value problems, AMS Chelsea Publishing, Providence, RI, 1.. H. O. Akerstedt, J. Nycander, V. P. Pavlenko, Three-dimensional stability of drift vortices, Phys. Plasmas , no. 1, C. Cao, A. Farhat, E. S. Titi, Global well-posedness of an inviscid three-dimensional pseudo- Hasegawa-Mima model, Comm. Math. Phys , no. 1,

16 16 C. CAO, Y. GUO, AND E. S. TITI 4. C. Cao, E. S. Titi, Global regularity criterion for the 3D Navier-Stokes equations involving one entry of the velocity gradient tensor, Arch. Ration. Mech. Anal. 11, no. 3, J. G. Charney, On the scale of atmospheric motions, Geofys. Publ. Norske Vid.-Akad. Oslo 17, no., 17 pp. 6. H. Gao, A. Zhu, The global strong solutions of Hasegawa-Mima-Charney-Obukhov equation, J. Math. Phys. 46 5, no. 8, 83517, 6 pp. 7. B. Guo, Y. Han, Existence and uniqueness of global solution of the Hasegawa-Mima equation, J. Math. Phys. 45 4, no. 4, K. Julien, E. Knobloch, R. Milliff, J. Werne, Generalized quasi-geostrophy for spatially anisotropic rotationally constrained flows, J. Fluid Mech , K. Julien, E. Knobloch, Reduced models for fluid flows with strong constraints, J. Math. Phys. 48 7, no. 6, 6545, 34 pp. 1. A. Hasegawa, K. Mima, Pseudo-three-dimensional turbulence in magnetized nonuniform plasma, Phys. Fluids , A. Hasegawa, K. Mima, Stationary spectrum of strong turbulence in magnetized nonuniform plasma, Phys. Rev. Lett , A. Hasegawa, M. Wakatani, Plasma edge turbulence, Phys. Rev. Lett , A. Hasegawa, M. Wakatani, A collisional drift wave description of plasma edge turbulence, Phys. Fluids, , S. Kondo, A. Tani, Initial boundary value problem for model equations of resistive drift wave turbulence, SIAM J. Math. Anal , no., O. A. Ladyzhenskaya, The boundary value problems of mathematical physics, Springer-Verlag, New York, A. Majda, Introduction to PDEs and waves for the atmosphere and ocean, Courant Lecture Notes in Mathematics, J. D. Meiss, W. Horton, Solitary drift waves in the presence of magnetic shear, Phys. Fluids , A. M. Obukhov, Concerning the geostrophic wind, Russian Izvestiya Akad. Nauk SSSR. Ser. Geograf. Geofiz , L. Paumond, Some remarks on a Hasegawa-Mima-Charney-Obukhov equation, Phys. D 195 4, no. 3-4, J. Pedlosky, Geophysical fluid dynamics, Springer-Verlag, M. Sprague, K. Julien, E. Knobloch, J. Werne, Numerical simulation of an asymptotically reduced system for rotationally constrained convection, J. Fluid Mech , R. Zhang, B. Guo, Dynamical behavior for the three-dimensional generalized Hasegawa-Mima equations, J. Math. Phys. 48 7, no. 1, 173, 11 pp. Department of Mathematics & Statistics, Florida International University, Miami, Florida 33199, USA address: caoc@fiu.edu Department of Mathematics & Statistics, Florida International University, Miami, Florida 33199, USA address: yanguo@fiu.edu Department of Mathematics, Texas A&M University, College Station, TX 77843, USA AND Department of Computer Science and Applied Mathematics, Weizmann Institute of Science, Rehovot 7611 Israel address: titi@math.tamu.edu, edriss.titi@weizmann.ac.il