REGULARITY OF WEAK SOLUTIONS TO THE LANDAU-LIFSHITZ SYSTEM IN BOUNDED REGULAR DOMAINS

Transcription

1 Electronic Journal of Differential Equations, Vol , No. 141, pp ISSN: URL: or ftp ejde.math.txstate.edu login: ftp REGULARITY OF WEAK SOLUTIONS TO THE LANDAU-LIFSHITZ SYSTEM IN BOUNDED REGULAR DOMAINS KÉVIN SANTUGINI-REPIQUET Abstract. In this paper, we study the regularity, on the boundary, of wea solutions to the Landau-Lifshitz system in the framewor of the micromagnetic model in the quasi-static approximation. We establish the existence of global wea solutions to the Landau-Lifshitz system whose tangential space gradient on the boundary is square integrable. 1. Introduction The Landau-Lifshitz equation models the behavior of ferromagnetic materials, and it is used in crystallography. This nonlinear systems of partial differential equations is m = m h αm m h. 1.1 with h magnetic excitation depending on m magnetization, and α R + dampening parameter. Existence of global wea solutions to the Landau-Lifshitz equation have been proved in [14, 1, 9]. However, these solutions are only required to belong to L R + ; H 1 O and to H 1 O 0, T for all time T > 0. Besides, they are not unique. There are existence results for of more regular solutions in the case when h = m, see [10, 11]. For 3D and a full characterization of the unique good solution in 2D has been obtained by Harpes[8]. These results are not easily generalized to more complicated forms of h as they often rely on harmonic analysis. Some authors have studied the regularity of stationary maps when the excitation h has a more complicated form: Critical points of the energy are regular away from a set of zero two-dimensional Hausdorff measure; see Carbou [3] and Hardt and Kinderlehrer [7]. Using standard analysis, one can easily prove that any wea solution m belongs to C[0, + ; H 1 wo, and to the Niol sii space L 2 0, T ; B1, O 2 3, and satisfy m L 2 0, T ; L 1 O, see [13, 6.2]. In this paper, we establish that for any initial condition, at least one wea solution to the Landau-Lifshitz system has a trace that belongs to L 2 0, T ; H 1 O Mathematics Subject Classification. 35D10, 35K55, 35K60. Key words and phrases. Micromagnetism; Landau-Lifshitz equation; wea solutions; boundary regularity. c 2007 Texas State University - San Marcos. Submitted July 18, Published October 24, Wor done at Institut Galilée, Université Paris 13. France. 1

2 2 K. SANTUGINI-REPIQUET EJDE-2007/141 This result is not limited to the case h = A m and can be generalized to a wide range of forms of h. In section 1, we introduce a common notation used throughout this paper. Then, in section 2, we recall briefly the micromagnetic model. We recall nown results concerning the existence of wea solutions in section 3. We state and prove our main theorem in section 4. Notation. Given an open set O, we denote by L p O the set of all measurable functions u over O such that O u p dx < +. This is a Banach space for the norm 1/p. u Lp O = u dx p For any integer m 0, we denote by H m O, the space of all measurable functions u over O such that for any multi-indices α, α < m, D α u belongs to L 2 O. This is an Hilbert space for the norm 1/2. u H m O = D α u 2 L O 2 α m We set H m O = H m O 3 and L p O = L p O 3. By O, we denote the Lebesgue measure of set O. Given a smooth surface O, ν represents the unit outward normal vector to the surface, u ν the normal trace of u on O, and Tu the tangential gradient of u on O. In all this paper, is a bounded open set of R 3 with a smooth boundary. We also define = 0, T. 2. The micromagnetic model In this section, we recall briefly the micromagnetic model. We begin by introducing the more common energies and excitations that model completely the static behavior, then we introduce the nonlinear PDE that models the evolution problem. From now on, represents the domain filled with a ferromagnetic material The static problem. The magnetic state of a ferromagnetic material is represented by two vector fields: the magnetization m and the magnetic excitation h. In the micromagnetic model the magnetization must verify a non convex constraint: m = 1 in and be null outside. The excitation h depends on m. To each interaction p is associated an energy E p m and an operator H p related by: DE p m v = H p m v dx, E p 0 = 0. The excitation contributed by p is h p = H p m and the total excitation is h = p h p. In this paper we consider three interactions, see Brown [2] for details: Exchange: The exchange energy and its associated excitation are given by: E e m = A m 2 dx, H e m = A m, 2 where A > 0. Anisotropy: E a m = 1 2 O m Km dx, H a m = Km,

3 EJDE-2007/141 REGULARITY OF WEAK SOLUTIONS 3 where K is a symmetric positive quadratic form. Demagnetization field: E d m = 1 H d m 2 dx, H d m = h d, 2 R 3 where h d is the only solution in L 2 R 3 to the magnetostatic system in the sense of distributions. divh d + m = 0, curlh d = 0, The demagnetization field operator has been extensively studied by Friedman [4, 5, 6]. This operator is symmetric negative, is linear continuous from L p R 3 to L p R 3, 1 < p < +, and satisfies H d m m dx = H d m 2 dx. R 3 The LL 2 R 3, L 2 R 3 norm of H d is less than 1. We also define H = H e + H a + H d, H d,a = H d + H a, H e,a = H e + H a, E = E e + E a + E d. The solutions to the static problem are the local minimizers of the total energy E satisfying the constraint m = 1 a.e. in The evolution system. The Landau-Lifshitz system models the evolution of ferromagnetic materials. Let α > 0 be a dampening parameter. The Landau- Lifshitz system comprises the nonlinear system the initial condition the non convex constraint m = m Hm αm m Hm in R +, 2.1a and the Neumann boundary conditions m, 0 = m 0 in, 2.1b m = 1 in R +, 2.1c m ν = 0 on R+. 2.1d 3. Known results We give the definition of wea solutions to the Landau-Lifshitz system 2.1. Definition 3.1. Given m 0 in H 1, m 0 = 1 a.e. in, we call m a wea solution to the Landau-Lifshitz system 2.1 if 1 For all T > 0, m belongs to H 1 0, T, and m = 1 a.e. in R +.

4 4 K. SANTUGINI-REPIQUET EJDE-2007/141 2 For all φ in H 1 0, T, m φ α T 3 = 1 + α 2 A i=1 m m φ m m φ 1 + α 2 m H d,a m φ. 3 m, 0 = m 0 in the sense of traces. 4 For all T > 0, EmT + α m α 2 dt dx Em0, where 3.1a 3.1b Eu = A 2 u 2 L Ku 2 L H dm 2 L 2. It was proved by Alouges and Soyeur [1] that the Landau-Lifshitz system 2.1 has at least one wea solution when h = A m. This result was generalized to the full h in [9]: Theorem 3.2. Let m 0 belongs to H 1, such that m 0 = 1 a.e. in. Then, there exists at least one solution to the Landau-Lifshitz system m in the sense of definition 3.1. Proof. The proof is based on the study of a penalized system whose solution converges to a wea solution to the Landau-Lifshitz system. The penalized system was: α m + m m = 1 + α 2 Hm m 2 1m, 3.2a m = 0, ν 3.2b m, 0 = m c See [1, 9] for details. See also [12] for a generalization with surface energies. We will prove a slightly stronger result. However, we need a uniform L R + bound of m. But, the non locality of the H d operator prevents us from proving this result. For this reason, we introduce another penalized, less simple, system 4.1 to prove the existence of solutions more regular on the boundary. 4. Existence of a solution with H 1 regularity in space on the boundary We state our main result as follows. Theorem 4.1. Let m 0 belong to H 1, such that m 0 = 1 a.e. in. Then, there exists at least one solution to the Landau-Lifshitz system m in the sense of definition 3.1 such that γm L 2 0, T ; H 1.

5 EJDE-2007/141 REGULARITY OF WEAK SOLUTIONS 5 This theorem is a consequence of propositions 4.3 and 4.6. The rest of this section is dedicated to the proof of this theorem. We first introduce a penalized system: = 1 + α 2 H e,a m + 1 {m 0}H d m α m + m m H d m m m m m m 2 1m in R +, 4.1a m ν = 0 on R+, 4.1b m, 0 = m 0 in. 4.1c We first show that this penalized system has a wea solution that converges to a wea solution to the Landau-Lifshitz system as tends to +. We then show that the L 2 0, T ; H 1 norm of m is bounded independently of. This requires a uniform L R + bound of m first: something we could not establish for system 3.2, hence the modification of the penalized system Properties of solutions to the penalized system 4.1. One can easily adapt the proof of existence of solution to the original penalized system 3.2 found in [9, 1] to the modified penalized system 4.1: Proposition 4.2. Let m 0 in H 1, m 0 = 1 a.e. in. Then, there exists a solution m to system 4.1 in L R + ; H 1 H 1 0, T, i.e. satisfying: α m + m m ψ = 1 + α 2 A m ψ α 2 1 {m 0}H d,a m ψ 1 + α 2 Hd m m m m 2 ψ 1 + α 2 m 2 1m ψ, for all ψ in H 1 0, T, and for all time T > 0, all η > 0: Em α, T α 2 η T + Em 0 exp T 2η a m 2 dt + L 2 4 m, T L 2 4.2b Proof. The proof is the same as the one in [1, 9] with a minor complication arising from the supplementary term. Let w 1,..., w n,... be the orthonormal hilbertian basis of L 2 comprising the eigenfunctions of the Laplace operator with homogeneous Neumann boundary conditions. The w i belongs to C. Let V n be the subspace of L 2 spanned by functions w 1,..., w n. Let P n be the orthogonal projector on V n in L 2. The basis w 1,..., w n,... is also an hilbertian orthogonal basis of H and P n

6 6 K. SANTUGINI-REPIQUET EJDE-2007/141 is also an orthogonal projector in H 1. We search for m n with the form m n = n i=1 φ itw i x, where the φ i are in C R + ; R 3 such that α m n + P n m n m n = 1 + α2 P n Hm n H dm n m nm n m n 2 + n α 2 P n m n 2 1m n m n, 0 = P n m 0. Let Φ n = φ 1,..., φ n. Equation 4.3a is equivalent to dφ n dt AΦ n t dφ n dt = F Φ n t, 4.3a 4.3b where F is of class C thans to the n 1 in the denominator which serves no other purpose and Φ n t AΦ n t is linear continuous, thus smooth. Moreover, AΦ is an antisymmetric matrix for all Φ. So the matrix I AΦ is nonsingular and the function Φ I AΦ 1 is of class C. By Cauchy-Lipshitz, Φ i t exists locally in time. Equation 4.3a can be expressed as α m n + m n m n ψ = 1 + α 2 A 1 + α α 2 m n ψ α 2 Hd m n m nm n m n 2 + n 1 m n 2 1m nψ, ψ H d,a m n ψ 4.4 for all ψ in V n C R + ; R 3. In 4.4, we tae ψ = m n, we obtain for all T > 0: T Em n, T + α m n α 2 0 dt + L 2 4 m n, T L 2 Em P nm L 2 Hd m n m nm n m n 2 + n 1 m n. Therefore, for all η > 0, all time T > 0: Em α n, T α 2 η T 0 Em P nm L η By Gronwall, for all η > 0, all time T > 0: T m n 2 dt + L 2 4 m n, T L 2 T m n Em α n, T α 2 η m n 2 0 dt + L 2 4 m n, T L 2 + Em P nm T L 2 exp 2η Thus, m n exists in global time. Since, P n m 0 tends to m 0 in H 1, 4 P nm L 2 tends to 0 as n tends to +. Therefore, m n is bounded in H 1 and in L 0, T ; H 1, independently of n. There exists m in L 0, T ; H 1 and in

7 EJDE-2007/141 REGULARITY OF WEAK SOLUTIONS 7 H 1 0, T and v in L 2, such that, modulo a subsequence, as n tends to + : m n m strongly in L 2 0, T, 4.6a m n m wealy in H 1 0, T, 4.6b m n m wealy in L 0, T ; H 1, 4.6c H d m 1 n m nm n {m =0} m n 2 + n 1 v wealy in L 2 0, T. 4.6d Also, by Aubin s lemma, for all 1 < p < +, 1 < q < 6, T > 0, m n m strongly in L p 0, T ; L q. 4.6e The limit m has the required properties: P n m 0 converges to m 0 in H 1. Computing the limit of 4.5, yields 4.2b. Since on the set {x, t m = 0}, m = 0 and m = 0 a.e., computing the limit of 4.4 yields first v = 1 {m =0}H d m, then 4.2a for all ψ in n=1 V n C [0, T ]; R 3. By density, 4.2a holds for all ψ in H 1 0, T. Now, we can prove that m converges to a wea solution to the Landau-Lifshitz system. Proposition 4.3. Let m 0 be in H 1, m 0 = 1 a.e. in. Let m be a wea solution to the penalized system 4.1. Then, there exists a subsequence of m, that converges to a wea solution m to the Landau-Lifshitz system 2.1 wealy in H 1 0, T. Proof. By 4.2b, m is bounded in L 0, T ; H 1 and in H 1 0, T. Besides, m L 0,T ;L 2 is bounded. There exists m in H1 0, T, such that, up to a subsequence, m m strongly in L 2 0, T, 4.7a m m wealy in H 1 0, T, 4.7b m m wealy in L 0, T ; H 1, 4.7c m strongly in L 2 0, T, 4.7d Also, by Aubin s lemma, for all 1 < p < +, 1 < q < 6, T > 0, m m strongly in L p 0, T ; L q. 4.7e Obviously, m, 0 = m 0. By 4.7d, m = 1. We compute the limit of 4.2b as tends to + as in [1, 9]. Then, we have η tend to 0: energy inequality 3.1b is satisfied. It only remains to prove that m satisfy 3.1a: in 4.2a, we tae ψ = m φ, with φ in C 0, T ; R 3 and tae the limit as tends to +. Since the supplementary term containing H d m m m m disappears, we can then conclude m as in [1, 9] and obtain 3.1a for all φ in C 0, T. By density, 3.1a also holds for all φ in H 1 0, T. We establish an important proposition that we were unable to prove for the original penalized system 3.2 due to the non locality of H d.

8 8 K. SANTUGINI-REPIQUET EJDE-2007/141 Proposition 4.4. Let m 0 be in H 1, m 0 = 1 a.e. in. Let m be a wea solution to system 4.1. Then, m satisfy m 1 a.e. in R +. Proof. We follow Alouges-Soyeur [1], and introduce the map g : R R defined by 0 if x < 0, gx = x if 0 x 1, 1 if x 1, We set ψx, t = g m x, t 2 1m x, t. The function ψ belongs to H 1 0, T and ψ = 2g m 2 1 m m m + g m 2 1 m. Reporting ψ in 4.2a yields α g m T x, t 2 1m x, t m 3 = 1 + α 2 A 2g m 2 1 m m 2 i=1 1 + α 2 g m 2 1 m 2 T 1 + α 2 g m x, t 2 1m Km T 1 + α 2 m 2 1 m 2 g m x, t 2 1. Should we have used system 3.2 instead, then the term containing the global operator would have been very difficult, if not outright impossible, to estimate. Therefore, g m x, t 2 1m x, t m 0. Thus, for all T > 0, G m, T 2 1 dx G m dx = 0, where Gx = x 0 gs ds. Since G 0, G m, T 2 1 = 0 a.e. in for all T > 0. Therefore, m 1 a.e. in R +. As a corollary to proposition 4.4, we have Corollary 4.5. Any wea solution m to system 4.1 belongs to the space H 2,1 0, T Uniform H 1 bound of the penalized solution. To prove Theorem 4.1, we only need to prove the following proposition. Proposition 4.6. Let m 0 be in H 1, m 0 = 1. Let m be a wea solution to the penalized system 4.1. Then, for all time T > 0 the quantity T γm 2 L 2 0,T is bounded uniformly in. Proof. Let φ i be in C 0, T for any integer i, 1 i 3. By corollary 4.5 and proposition 4.4, we can multiply 4.1a by m φ i and integrate over any open

9 EJDE-2007/141 REGULARITY OF WEAK SOLUTIONS 9 set O with a smooth boundary: m α = 1 + α α 2 m φ i + m m m φ i 1 + α α 2 A 1 {m 0}H d,a m m φ i H d m m m m m m φ i m 2 1 m m φ i } {{ } I 3 m m 1 m 2 j=1 x j x j 2 x j φ i }{{} II However, in the above equality, I = 4 m φ i + m νi φ i dσx dt, 4 and II = + 3 j=1 m m m φi + 1 x j x j 2 m φ i dσx dt 1 ν 2 m 2 φ i m 2 ν i φ i dσx dt.

10 10 K. SANTUGINI-REPIQUET EJDE-2007/141 Therefore, m α 1 + α α α 2 A 3 j=1 1 + α 2 A α 2 4 = 1 + α α 2 A 1 + α 2 A 2 m φ i + m m m φ i 1 {m 0}H d,a m m φ i H d m m m m m m m φ i m x j φi x j m 2 φ i m φ i m νi φ i dσx dt m m φ i dσx dt ν m 2 ν i φ i dσx dt, 4.8 for all φ i in C 0, T. In 4.8, we choose φ i independent of the time t such that φ i = ν i on O and sum over i. We obtain, denoting by φ the vector valued function φ 1, φ 2, φ 3, m α 1 + α α α 2 A φ m + m m φ m 3 i,j=1 1 + α 2 A α 2 4 = 1 + α α 2 A 1 {m 0}H d,a m φ m H d m m m m m m φ m m 2 div φ m x j φi x j m div φ m dσx dt m 2 dσx dt ν

11 EJDE-2007/141 REGULARITY OF WEAK SOLUTIONS α 2 A 2 m 2 dσx dt. The left hand-side is bounded uniformly in. Therefore, A 2 Tγm 2 L 2 A m 2 ν + 4 γm L 2 CO. 2 L 2 Since m ν = 0 on, taing O = yields the wanted result. We derive from the previous proof the following corollary. Corollary 4.7. Let m 0 be in H 1, m 0 = 1 a.e. in. Let m be a wea solution to the penalized system 4.1. Then, the quantity γm L 2 0,T is bounded uniformly in. Conclusion. In this paper, we have proved the existence of wea solutions to the Landau-Lifshitz system with an H 1 regularity in space on the boundary of the domain. This result holds for very general form of h and is not limited to the case h = A m. These inds of results are important because there is currently no perfect concept of what is a good wea solution to the Landau-Lifshitz system. Any criteria allowing to discriminate among those wea solutions is always welcome. It is natural to prefer among wea solutions those that are more regular. This result is also interesting as it opens the possibility to use first order transmission conditions between adjacent domains. Acnowledgments. The author wansts to than S. Labbé and L. Halpern for their support. References [1] F. Alouges and A. Soyeur. On global wea solutions for Landau-Lifshitz equations : existence and nonuniqueness. Nonlinear Analysis. Theory, Methods & Applications, 1811: , [2] W. F. Brown. Micromagnetics. Interscience Publishers, [3] G. Carbou. Regularity for critical points of a non local energy. Calc. Var. PDE, 5: , [4] M. J. Friedman. Mathematical study of the nonlinear singular integral magnetic field equation I. SIAM J. Appl. Math., 391:14 20, Août [5] M. J. Friedman. Mathematical study of the nonlinear singular integral magnetic field equation II. SIAM J. Numer. Math., 184: , August [6] M. J. Friedman. Mathematical study of the nonlinear singular integral magnetic field equation III. SIAM J. Math. Anal., 124: , July [7] R. Hardt and D. Kinderlehrer. Some regularity results in micromagnetism. Communications in Partial Differential Equations, 257 8: , [8] P. Harpes. Uniqueness and bubbling of the 2-dimensional Landau-Lifshitz flow. Calc. Var. Partial Differential Equations, 202: , [9] S. Labbé. Simulation numérique du comportement hyperfréquence des matériaux ferromagnétiques. PhD thesis, Université Paris 13, Décembre [10] C. Melcher. Existence of partially regular solutions for Landau-Lifshitz equations. preprint, Institute for mathematics and its applications, University of Minnesota, November [11] R. Moser. Partial regularity for the Landau-Lifshitz equation in small dimensions. Technical Report 26, Max-Planc-Institut fur Mathemati in den Naturwissenschaften, Leipzig, 2002.

12 12 K. SANTUGINI-REPIQUET EJDE-2007/141 [12] K. Santugini. Solutions to the Landau-Lifshitz system with nonhomogenous boundary conditions. preprint submitted, March [13] K. Santugini-Repiquet. Matériaux ferromagnetiques: influence d un espaceur mince non magnétique, et homogénéisation d agencements multicouches, en présence de couplage sur la frontière. Thése de doctorat, Université Paris 13, Villetaneuse, dec [14] A. Visintin. On Landau-Lifshitz equations for ferromagnetism. Japan J. Appl. Math., 21:69 84, Kévin Santugini-Repiquet Institut Mathématique de Bordeaux, Université Bordeaux 1, 351 cours de la libération, Talence Cedex, France address: Kevin.Santugin@math.u-bordeaux1.fr URL: santugin