LOCAL WELL-POSEDNESS FOR THE HALL-MHD EQUATIONS WITH FRACTIONAL MAGNETIC DIFFUSION DONGHO CHAE 1, RENHUI WAN 2 AND JIAHONG WU 3 arxiv:144.486v2 [math.ap] 3 Jan 215 Abstract. The Hall-magnetohydrodynamics (Hall-MHD) equations, rigorously derived from kinetic models, are useful in describing many physical phenomena in geophysics and astrophysics. This paper studies the local well-posedness of classical solutions to the Hall-MHD equations with the magnetic diffusion given by a fractional Laplacian operator, ( ) α. Due to the presence of the Hall term in the Hall-MHD equations, standard energy estimates appear to indicate that we need α 1 in order to obtain the local well-posedness. This paper breaks the barrier and shows that the fractionalhall-mhdequationsarelocallywell-posedforanyα > 1 2. The approachhere fully exploits the smoothing effects of the dissipation and establishes the local bounds for the Sobolev norms through the Besov space techniques. The method presented here may be applicable to similar situations involving other partial differential equations. 1. Introduction This paper focuses on the Hall-magnetohydrodynamics (Hall-MHD) equations with fractional magnetic diffusion, t u+u u+ p = B B, t B +u B + (( B) B)+( ) α B = B u, u =, B =, u(x,) = u (x), B(x,) = B (x), (1.1) where x R d with d 2, u = u(x,t) and B = B(x,t) are vector fields representing the velocity and the magnetic field, respectively, p = p(x,t) denotes the pressure, α > is a parameter and the fractional Laplacian ( ) α is defined through the Fourier transform, ( ) α f(ξ) = ξ 2α f(ξ). For notational convenience, we also use Λ for ( ) 1 2. The Hall-MHD equations with the usual Lapalcian dissipation were derived in [1] from kinetic models. The Hall- MHD equations differ from the standard incompressible MHD equations in the Hall term (( B) B), which is important in the study of magnetic reconnection (see, e.g., [8, 13]). The Hall-MHD equations have been mathematically investigated in several works ([1, 4, 5, 6, 7]). Global weak solutions of (1.1) with both u and B and local classical solutions of (1.1) with B (with or without u) were obtained in [4]. In addition, a blowup criterion and the global existence of small classical solutions were also established in [4]. These results were later sharpened by [5]. 21 Mathematics Subject Classification. 35Q35, 35B65, 35Q85, 76W5. Key words and phrases. Hall-MHD equations, fractional magnetic diffusion, local well-posedness. 1
2 D. CHAE, R. WAN AND J. WU We examine the issue of whether or not (1.1) is locally well-posedness when the fractional power α < 1. Previously local solutions of (1.1) were obtained for α = 1 ([4, 5]). Standard energy estimates appear to indicate that α 1 is necessary in order to obtain local bounds for the solutions in Sobolev spaces. This requirement comes from the estimates of the regularity-demanding Hall term (( B) B). To understand more precisely the issue at hand, we perform a short energy estimate on the essential part of the equation for B, The global L 2 -bound t B + (( B) B)+( ) α B =. follows from the simple fact (( B) B) B = t B(t) 2 L +2 Λ α B(τ) 2 2 L 2dτ = B 2 L (1.2) 2 (( B) B) ( B) =. To obtain the H 1 -bound, we invoke the equation for B 2 L, 2 1 d d 2dt B 2 L 2 + Λα B 2 L = 2 i (( B) B) i B. i=1 We can indeed shift one-derivative, namely i (( B) B) i B = Hölder s inequality allows us to conclude that (( B) i B) i B. 1 d 2dt B 2 L 2 + Λα B 2 L B 2 L 2 B L B L 2. Therefore, it appears that we need α 1 in order to bound the term B L 2 on the right-hand side. More generally, the energy inequality involving the H σ -norm d dt B 2 H σ + Λα B 2 H σ C B H σ B L B H σ. also appears to demand that α 1 in order to bound B H σ. This paper obtains the local existence and uniqueness of solutions to (1.1) with any α > 1. More precisely, we prove the following theorem. 2 Theorem 1.1. Consider (1.1) with α > 1 2. Assume (u,b ) H σ (R d ) with σ > 1+ d 2, and u = B =. Then there exist T = T ( (u,b ) H σ) > and a unique solution (u,b) of (1.1) on [,T ] such that In addition, for any σ < σ, (u,b) L ([,T ];H σ (R d )). (u,b) C([,T ];H σ (R d )) and (u(t),b(t)) H σ is continuous from the right on [,T ).
HALL-MHD EQUATIONS 3 The essential idea of proving Theorem 1.1 is to fully exploit the dissipation in the equation for B and estimate the Sobolev norm (u,b) H σ via Besov space techniques. We identify H σ with the Besov space B2,2 σ and suitably shift the derivatives in the nonlinear term. The definition of Besov spaces and related facts used in this paper are provided in the appendix. The rest of this paper is divided into two sections followed by an appendix. Section 2 states and proves the result for the local a priori bound. Section 3 presents the complete proof of Theorem 1.1. The appendix supplies the definitions of the Littlewood-Paley decomposition and Besov spaces. 2. Local a priori bound This section establishes a local a priori bound for smooth solutions of (1.1), which is the key component in the proof of Theorem 1.1. The result for the local a priori bound can be stated as follows. Proposition 2.1. Consider (1.1) with α > 1. Assume the initial data (u 2,B ) H σ (R d ) with σ > 1 + d. Let (u,b) be the corresponding solution. Then, there exists 2 T = T ( (u,b ) H σ) > such that, for t [,T ], and T (u(t),b(t)) H σ C(α,T, (u,b ) H σ) Λ α B(s) 2 H σ ds C(α,T, (u,b ) H σ). Proof of Proposition 2.1. TheproofidentifiestheSobolevspaceH σ withthebesovspace B σ 2,2 and resorts to Besov space techniques. Let l 1 be an integer and let l denote the homogeneous frequency localized operator. Applying l to (1.1) yields t l u+ l (u u)+ l p = l (B B), t l B + l (u B)+ l (( B) B)+( ) α l B = l (B u). Taking the inner product with ( l u, l B) and integrating by parts, we have where 1 d ( l u 2 L 2dt + ) lb 2 2 L 2 +C 2 2αl l B 2 L = K 2 1 +K 2 +K 3 +K 4 +K 5, (2.1) K 1 = K 3 = K 5 = [ l,u ]u l u, K 2 = [ l,b ]B l u, K 4 = l (( B) B) l B. Note that we have used the standard commutator notation, [ l,u ]B l B, [ l,b ]u l B, [ l,u ]u = l (u u) u ( l u)
4 D. CHAE, R. WAN AND J. WU and applied the lower bound, for a constant C >, l B ( ) α l B C 2 2αl l B 2 L 2. Using the notion of paraproducts, we write K 1 = K 11 +K 12 +K 13, where K 11 = ( l (S k 1 u k u) S k 1 u l k u) l u, k l 2 K 12 = k l 2 K 13 = ( l ( k u S k 1 u) k u l S k 1 u) l u, ( l ( k u k u) k u l k u) l u with k = k 1 + k + k+1. By Hölder s inequality and a standard commutator estimate, K 11 C S l 1 u L l u L 2 k u L 2 C u L l u L 2 k l 2 k l 2 k u L 2. Since the summation over k for fixed l above consists of only a finite number of terms and, as we shall later in the proof, the norm generated by each term is a multiple of that generated by the typical term, it suffices to keep the typical term with k = l and ignore the summation. This would help keep our presentation concise. We will invoke this practice throughout the rest of the paper. By Hölder s inequality, K 12 is bounded by K 12 C u L l u 2 L 2. By Hölder s inequality and Bernstein s inequality, K 13 C l u L 2 u L 2 l k k u L 2. Therefore, K 1 C l u L 2 u L Similarly, K 2, K 3 and K 4 are bounded by ( l u L 2 + 2 l k k u L 2 K 2 C u L l B 2 L +C B 2 L lu L 2 l B L 2 + C u L l B L 2 2 l k k B L 2, ).
K 3 C B L l u L 2 HALL-MHD EQUATIONS 5 ( l B L 2 + ) 2 l k k B L 2, Using the simple fact and the vector identity K 4 C B L l u L 2 l B L 2 +C u L l B 2 L 2 +C B L l u L 2 2 l k k B L 2. (B ( l B)) l B = B ( B) = 1 (B B) (B )B, 2 we can rewrite K 5 as K 5 = ( l (B ( B)) B ( l B)) l B = [ l,b ]B l B (2.2) ( ) ) 12 + ( l (B B) ( l B) B l B. (2.3) The term in (2.2) can be estimated in a similar way as K 3. To estimate the term in (2.3), we use the the notion of paraproducts to write ( ) ) 12 ( l (B B) ( l B) B l B = K 51 +K 52 +K 53, where, K 51 = k l 2 K 52 = k l 2 K 53 = By Hölder s inequality, ( l (( S k 1 B) k B) ( l S k 1 B) k B) l B, ( l (S k 1 B ( k B)) (S k 1 B) ( l k B)) l B, ( l ( ( )) ) 1 2 kb k B ( l k B) k B l B. K 51 l (( S k 1 B) k B) ( l S k 1 B) k B L 2 l B L 2 C2 l S l 1 B L l B 2 L 2 C2l B L l B 2 L 2. By Hölder s inequality and a standard commutator estimate, K 52 C2 l B L l B 2 L 2. By Hölder s inequality and Bernstein s inequality, K 53 C2 l B L l B L 2 2 l k k B L 2.
6 D. CHAE, R. WAN AND J. WU Therefore, K 5 C2 l B L l B L 2 ( l B L 2 + 2 l k k B L 2 Inserting the estimates above in (2.1), we obtain d ( l u 2 L dt + ) lb 2 2 L 2 +C 2 2αl l B 2 L 2 C ( u, B) L ( l u 2 L + lb 2 2 L 2) ( 2 ( ) 2 +C ( u, B) L 2 l k k u L 2) + 2 l k k B L 2 + C2 l B L l B 2 L 2 +C2l B L l B L 2 ) 2 l k k B L 2. Multiplying the inequality above by 2 2σl and summing over l 1, invoking the global bound for the L 2 -norm of (u,b) and the equivalence of the norms f 2 H 2 2σl σ l f 2 L 2, we have u(t) 2 H σ + B(t) 2 H σ + C t B(τ) 2 Hσ+α dτ t u 2 H + B 2 σ H +C ( u, B) σ L ( u(τ) 2 H σ + B(τ) 2 H σ)dτ +C t 2 (2σ+1)l B L l B 2 L 2dτ +C ( 2 t 2 (2σ+1)l B L 2 l k k B L 2) dτ. (2.4) To derive the inequality above, we have used Young s inequality for series convolution ( ) 2 2 2σl 2 l k k u L 2 = ( ) 2 2 (σ+1)(l k) 2 σk k u L 2 C 2 2σl l u 2 L 2 C u 2 H σ. We further bound the last two terms in (2.4), L 1 C t 2 (2σ+1)l B L l B 2 L dτ, 2 L 2 C ( 2 t 2 (2σ+1)l B L 2 l k k B L 2) dτ..
HALL-MHD EQUATIONS 7 Set θ = 1 1. For α > 1, θ (,1). By Hölder s inequality, 2α 2 t ( ) L 1 =C B L 2 2σl l B 2 θ ( ) L 2 2(σ+α)l 2 l B 2 (1 θ) L dτ 2 C C t t ( ) θ ( ) (1 θ) B L 2 2σl l B 2 L 2 2(σ+α)l 2 l B 2 L dτ 2 B 1 θ L B 2 H σ dτ + C 4 t B(τ) 2 H σ+αdτ. By Young s inequality for series convolution and an interpolation inequality, ( 2 t L 2 =C B L 2 (l k)(σ+1 2 ) 2 (σ+1 2 )k k B L 2) dτ C C t t B L B 2 dτ C H σ+1 2 B 1 θ L B 2 H σ dτ + C 4 t t B L B 2θ H σ B 2(1 θ) H σ+α dτ B(τ) 2 H σ+αdτ. Inserting the estimates above in (2.4) and invoking the embedding inequalities we have B L C B H σ for σ > 1+ d 2, u(t) 2 H σ + B(t) 2 H σ + C t B(τ) 2 Hσ+α dτ t u 2 H + B 2 σ H +C ( ) u(t) 2 σ H σ + B(t) 2 γ H dτ, (2.5) σ for a constant γ > 1. This inequality implies a local bound for u(t) 2 H σ + B(t) 2 H σ, namely for some T = T ( (u,b ) H σ) > such that, for t [,T ], and u(t) H σ + B(t) H σ C(u,B,α,T ) T This completes the proof of Proposition 2.1. B(τ) 2 Hσ+αdτ <. (2.6) This section proves Theorem 1.1. 3. Local existence and uniqueness Proof of Theorem 1.1. The local existence and uniqueness can be obtained through an approximation procedure. Here we use the Friedrichs method, a smoothing approach through filtering the high frequencies. For each positive integer n, we define Ĵ n f(ξ) = χ Bn (ξ) f(ξ),
8 D. CHAE, R. WAN AND J. WU where B n denotes the closed ball of radius n centered at and χ Bn denotes the characteristic functions on B n. Denote } Hn σ {f H σ (R d ), supp f B n. We seek a solution (u,b) Hn σ satisfying t u+j n P(J n Pu J n Pu) = J n P(J n PB J n PB), t B +J n P(J n Pu J n PB)+J n P( (( J n PB) J n PB)) +( ) α J n PB = J n P(J n PB J n Pu), u(x,) = (J n u )(x), B(x,) = (J n B )(x), (3.1) where P denotes the projection onto divergence-free vector fields. For each fixed n 1, it is not very hard, although tedious, to verify that the righthand side of (3.1) satisfies the Lipschitz condition in Hn σ and, by Picard s theorem, (3.1) has a unique global (in time) solution. The uniqueness implies that J n Pu = u, J n PB = B and ensures the divergence-free conditions u = and B =. Then, (3.1) is simplified to { t u+j n P(u u) = J n P(B B), t B +J n P(u B)+J n P( (( B) B))+( ) α B = J n P(B u). We denote this solution by (u n,b n ). As in the proof of Proposition 2.1, we can show that (u n,b n ) satisfies (u n,b n ) 2 H σ (un,bn ) 2 H σ +C t (u n (s),b n (s)) 2γ Hσ ds (3.2) for some γ > 1. Due to (u n,bn ) H σ (u,b ) H σ, this inequality is uniform in n. This allows us to obtain a uniform local bound sup (u n (t),b n (t)) H σ M(α,T, (u,b ) H σ). (3.3) t [,T ] As in (2.6), we also have the uniform local bound for the time integral T (Λ α B n )(s) 2 H σ ds M(α,T, (u,b ) H σ). Furthermore, these uniform bounds allow us to show that (u n,b n ) (u m,b m ) L 2 as n,m. (3.4) This is shown through standard energy estimates for (u n,b n ) (u m,b m ) L 2. The process involves many terms, but most of them can be handled in a standard fashion (see, e.g., [9, p.17]). We provide the detailed energy estimate for the term that is special here, namely the Hall term (( B) B). In the process of the energy estimates,
HALL-MHD EQUATIONS 9 we need to bound the term ( (( B n ) B n ) (( B m ) B m )) (B n B m )dx = ( (( (B n B m )) B n )) (B n B m )dx + ( (( B m ) (B n B m )) (B n B m )dx. (3.5) The first term on the right-hand side of (3.5) is zero, ( (( (B n B m )) B n )) (B n B m )dx = (( (B n B m )) B n ) ( (B n B m ))dx =. For the second term on the right of (3.5), by the simple vector identity we have (( B m ) (B n B m )) where we have used = (B n B m ) ( B m ) ( B m ) (B n B m ), ( (( B m ) (B n B m )) (B n B m )dx ( B m ) L d α B n B m L 2d d 2α Bn B m L 2 C Λ α B m 2 H σ Bn B m 2 L 2 + 1 8 Λα (B n B m ) 2 L 2, 2 f L d α C Λ α f H σ, f L 2d d 2α C Λα f L 2. Putting together the estimates for all the terms, we obtain d dt Bn B m 2 L 2 C Λα B m 2 H σ Bn B m 2 L 2 +C ( 1 n + 1 ). m Noticing that Λ α B m 2 Hσ is time integrable, Gronwall s inequality yields the desired convergence (3.4). Let (u,b) be the limit. Due to the uniform bound (3.3), (u,b) H σ for t [,T ]. By the interpolation inequality, for any < σ < σ, f H σ C σ f 1 σ σ L 2 we further obtain the strong convergence (u n,b n ) (u,b) H σ f σ σ H σ, as n and consequently, (u,b) C([,T ];H σ ). This strong convergence makes it easy to check that (u, B) satisfies the Hall-MHD equation in (1.1). In addition, the time continuity in (u,b) C([,T ];H σ ) allows to show the weak time continuity (u,b) C W ([,T ];H σ ) or t (u(x,t),b(x,t)) φ(x)dx is continuous
1 D. CHAE, R. WAN AND J. WU for any φ H σ. To show the right (in time) continuity of (u(t),b(t)) H σ, we make use of the energy inequality, for any t > t, t (u(t),b(t)) 2 H (u( t),b( t)) 2 σ H +C (u(s),b(s)) γ σ H ds, 2σ This inequality can be obtained in a similar fashion as (3.2). Then, By the weak continuity in time, lim (u(t),b(t)) H σ (u( t),b( t)) H σ. t t+ (u( t),b( t)) H σ lim (u(t),b(t)) H σ. t t+ The desired right (in time) continuity of (u(t),b(t)) H σ then follows. This completes the proof of Theorem 1.1. t Appendix A. Besov spaces This appendix provides the definitions of some of the functional spaces and related facts used in the previous sections. Materials presented in this appendix can be found in several books and many papers (see, e.g., [2, 3, 1, 11, 12]). We start with several notations. S denotes the usual Schwartz class and S its dual, the space of tempered distributions. S denotes a subspace of S defined by { } S = φ S : φ(x)x γ dx =, γ =,1,2, R d and S denotes its dual. S can be identified as S = S /S = S /P where P denotes the space of multinomials. To introduce the Littlewood-Paley decomposition, we write for each j Z A j = { ξ R d : 2 j 1 ξ < 2 j+1}. The Littlewood-Paley decomposition asserts the existence of a sequence of functions {Φ j } j Z S such that and supp Φ j A j, j= Φj (ξ) = Φ (2 j ξ) or Φ j (x) = 2 jd Φ (2 j x), Φ j (ξ) = { 1, if ξ R d \{},, if ξ =. Therefore, for a general function ψ S, we have Φ j (ξ) ψ(ξ) = ψ(ξ) for ξ R d \{}. j=
In addition, if ψ S, then That is, for ψ S, and hence j= HALL-MHD EQUATIONS 11 Φ j (ξ) ψ(ξ) = ψ(ξ) for any ξ R d. j= j= Φ j ψ = ψ Φ j f = f, f S in the sense of weak- topology of S. For notational convenience, we define j f = Φ j f, j Z. (A.1) Definition A.1. For s R and 1 p,q, the homogeneous Besov space B s p,q consists of f S satisfying f Bs p,q 2js j f L p l q <. We now choose Ψ S such that Then, for any ψ S, and hence Ψ(ξ) = 1 Ψ ψ + Ψ f + Φ j (ξ), ξ R d. j= Φ j ψ = ψ j= Φ j f = f in S for any f S. To define the inhomogeneous Besov space, we set, if j 2, j f = Ψ f, if j = 1, Φ j f, if j =,1,2,. j= (A.2) Definition A.2. The inhomogeneous Besov space B s p,q with 1 p,q and s R consists of functions f S satisfying f B s p,q 2 js j f L p l q <. The Besov spaces B p,q s and Bp,q s with s (,1) and 1 p,q can be equivalently defined by the norms (Rd ) ( f(x+t) f(x) L p) q 1/q f Bs p,q = dt, t d+sq
12 D. CHAE, R. WAN AND J. WU ) (Rd ( f(x+t) f(x) L p) q 1/q f B s p,q = f L p + dt. t d+sq When q =, the expressions are interpreted in the normal way. Many frequently used function spaces are special cases of Besov spaces. The following proposition lists some useful equivalence and embedding relations. Proposition A.3. For any s R, For any s R and 1 < q <, H s B s 2,2, Hs B s 2,2. B s q,min{q,2} W s q B s q,max{q,2}. In particular, B q,min{q,2} Lq B q,max{q,2}. For notational convenience, we write j for j. There will be no confusion if we keep inmindthat j sassociatedwiththehomogeneousbesovspacesisdefinedin(a.1)while those associated with the inhomogeneous Besov spaces are defined in (A.2). Besides the Fourier localization operators j, the partial sum S j is also a useful notation. For an integer j, S j j 1 k= 1 where k is given by (A.2). For any f S, the Fourier transform of S j f is supported on the ball of radius 2 j. Bernstein s inequalities are useful tools in dealing with Fourier localized functions and these inequalities trade integrability for derivatives. The following proposition provides Bernstein type inequalities for fractional derivatives. k, Proposition A.4. Let α. Let 1 p q. 1) If f satisfies supp f {ξ R d : ξ K2 j }, for some integer j and a constant K >, then 2) If f satisfies ( ) α f L q (R d ) C 1 2 2αj+jd(1 p 1 q ) f L p (R d ). supp f {ξ R d : K 1 2 j ξ K 2 2 j } for some integer j and constants < K 1 K 2, then C 1 2 2αj f L q (R d ) ( ) α f L q (R d ) C 2 2 2αj+jd(1 p 1 q ) f L p (R d ), where C 1 and C 2 are constants depending on α,p and q only.
HALL-MHD EQUATIONS 13 Acknowledgements Chae was partially supported by NRF grant No.26-93854 and No.29-83521. Wu was partially supported by NSF grant DMS129153 and the AT&T Foundation at Oklahoma State University. References [1] M. Acheritogaray, P. Degond, A. Frouvelle and J.-G. Liu, Kinetic formulation and global existence for the Hall-Magneto-hydrodynamics system, Kinet. Relat. Models 4 (211), 91-918. [2] H. Bahouri, J.-Y. Chemin and R. Danchin, Fourier Analysis and Nonlinear Partial Differential Equations, Springer, 211. [3] J. Bergh and J. Löfström, Interpolation Spaces, An Introduction, Springer-Verlag, Berlin- Heidelberg-New York, 1976. [4] D. Chae, P. Degond and J.-G. Liu, Well-posedness for Hall-magnetohydrodynamics, Ann. Inst. H. Poincaré Anal. Non Linéaire 31 (214), 555-565. [5] D. Chae and J. Lee, On the blow-up criterion and small data global existence for the Hallmagnetohydrodynamics, J. Differential Equations 256 (214), 3835-3858. [6] D. Chae and M. Schonbek, On the temporal decay for the Hall-magnetohydrodynamic equations, J. Differential Equations 255 (213), 3971-3982. [7] D. Chae and S. Weng, Singularity formation for the incompressible Hall-MHD equations without resistivity, arxiv:1312.5519 [math.ap] 19 Dec 213. [8] H. Homann and R. Grauer, Bifurcation analysis of magnetic reconnection in Hall-MHD systems, Phys. D 28 (25), 59-72. [9] A.J. Majda and A.L. Bertozzi, Vorticity and Incompressible Flow, Cambridge University Press, Cambridge, UK, 21. [1] C. Miao, J. Wu and Z. Zhang, Littlewood-Paley Theory and its Applications in Partial Differential Equations of Fluid Dynamics, Science Press, Beijing, China, 212 (in Chinese). [11] T. Runst and W. Sickel, Sobolev Spaces of fractional order, Nemytskij operators and Nonlinear Partial Differential Equations, Walter de Gruyter, Berlin, New York, 1996. [12] H. Triebel, Theory of Function Spaces II, Birkhauser Verlag, 1992. [13] M. Wardle, Star formation and Hall effect, Astrophysics and Space Science 292 (24), 317-323. 1 Department of Mathematics, College of Natural Science, Chung-Ang University, Seoul 156-756, Republic of Korea E-mail address: dchae@cau.ac.kr 2 Department of Mathematics, Zhejiang University, Hanzhou 3127, China E-mail address: rhwanmath@163.com,212352@zju.edu.cn 3 Department of Mathematics, Oklahoma State University, Stillwater, OK 7478, USA; and Department of Mathematics, Chung-Ang University, Seoul 156-756, Republic of Korea E-mail address: jiahong@math.okstate.edu