arxiv: v2 [math.na] 5 Feb 2017

Transcription

1 THE EDDY CURRENT LLG EQUATIONS: FEM-BEM COUPLING AND A PRIORI ERROR ESTIMATES MICHAEL FEISCHL AND THANH TRAN arxiv: v2 [mat.na] 5 Feb 2017 Abstract. We analyze a numerical metod for te coupled system of te eddy current equations in R 3 wit te Landau-Lifsitz-Gilbert equation in a bounded domain. Te unbounded domain is discretized by means of finite-element/boundary-element coupling. Even toug te considered problem is strongly nonlinear, te numerical approac is constructed suc tat only two linear systems per time step ave to be solved. We prove unconditional weak convergence of a subsequence of te finite-element solutions towards a weak solution. We establis a priori error estimates if a sufficiently smoot strong solution exists. Numerical experiments underlining te teoretical results are presented. Key words. Landau Lifsitz Gilbert equation, eddy current, finite element, boundary element, coupling, a priori error estimates, ferromagnetism AMS subject classifications. Primary 35Q40, 35K55, 35R60, 60H15, 65L60, 65L20, 65C30; Secondary 82D45 1. Introduction. Tis paper deals wit te coupling of finite element and boundary element metods to solve te system of te eddy current equations in te wole 3D spatial space and te Landau-Lifsitz-Gilbert equation LLG, te so-called ELLG system or equations. Te system is also called te quasi-static Maxwell-LLG MLLG system. Te LLG is widely considered as a valid model of micromagnetic penomena occurring in, e.g., magnetic sensors, recording eads, and magneto-resistive storage device [21, 24, 31]. Classical results concerning existence and non-uniqueness of solutions can be found in [5, 33]. In a ferro-magnetic material, magnetization is created or affected by external electro-magnetic fields. It is terefore necessary to augment te Maxwell system wit te LLG, wic describes te influence of a ferromagnet; see e.g. [19, 23, 33]. Existence, regularity and local uniqueness for te MLLG equations are studied in [18]. Trougout te literature, tere are various works on numerical approximation metods for te LLG, ELLG, and MLLG equations [3, 4, 10, 11, 19, 25, 26] te list is not exausted, and even wit te full Maxwell system on bounded domains [7, 8], and in te wole R 3 [17]. Originating from te seminal work [3], te recent works [25, 26] consider a similar numeric integrator for a bounded domain. Wile te numerical integrator of [26] treated LLG and eddy current simultaneously per time step, [25] adapted an idea of [8] and decoupled te time-steps for LLG and te eddy current equation. Our approac follows [25]. Tis work studies te ELLG equations were we consider te electromagnetic field on te wole R 3 and do not need to introduce artificial boundaries. Differently from [17] were te Faedo-Galerkin metod is used to prove existence of weak solutions, we extend te analysis for te integrator used in [3, 25, 26] to a finiteelement/boundary-element FEM/BEM discretization of te eddy current part on Supported by te Australian Researc Council under grant numbers DP and DP Scool of Matematics and Statistics, Te University of New Sout Wales, Sydney 2052, Australia m.feiscl@unsw.edu.au. Scool of Matematics and Statistics, Te University of New Sout Wales, Sydney 2052, Australia tan.tran@unsw.edu.au. 1

2 2 M. FEISCHL AND T. TRAN R 3. Tis is inspired by te FEM/BEM coupling approac designed for te pure eddy current problem in [13], wic allows us to treat unbounded domains witout introducing artificial boundaries. Two approaces are proposed in [13]: te so-called magnetic or H-based approac wic eliminates te electric field, retaining only te magnetic field as te unknown in te eddy-current system, and te electric or E-based approac wic considers a primitive of te electric field as te only unknown. Te coupling of te eddy-current system wit te LLG dictates tat te first approac is more appropriate, because tis coupling involves te magnetic field in te LLG equation rater tan te electric field; see 1. Te main results of tis work are weak convergence of te discrete approximation towards a weak solution witout any condition on te space and time discretization as well as a priori error estimates under te condition tat te exact strong solution is sufficiently smoot. In particular, te first result implies te existence of weak solutions, wereas te latter sows tat te smoot strong solution is unique. To te best of our knowledge, no suc results for te tangent plane sceme ave been proved for te LLG equation. Terefore, we present te proof for tis equation in a separate section, before proving te result for te ELLG system. As in [1], te proof is facilitated by use of an idea of [9] for te armonic map eat analyzed for LLG in [1], wic avoids te normalization of te solution in eac time-step, and terefore allows us to use a linear update formula. Tis also enables us to consider general quasi-uniform triangulations for discretization and removes te requirement for very sape-regular elements all diedral angles smaller tan π/2 present in previous works on tis topic. Te remainder of tis work is organized as follows. Section 2 introduces te coupled problem and te notation, presents te numerical algoritm, and states te main results of tis paper. Section 3 is devoted to te proofs of tese main results. Numerical results are presented in Section 4. Te final section, te Appendix, contains te proofs of some rater elementary or well-known results. 2. Model Problem & Main Results Te problem. Consider a bounded Lipscitz domain D R 3 wit connected boundary Γ aving te outward normal vector n. We define D := R 3 \ D, D T := 0, T D, Γ T := 0, T Γ, D T := 0, T D, and R 3 T := 0, T R3 for T > 0. For simplicity, we assume tat D is simply connected. We start wit te quasi-static approximation of te full Maxwell-LLG system from [33] wic reads as 1a 1b 1c 1d 1e m t αm m t = m H eff in D T, σe H = 0 in R 3 T, µ 0 H t + E = µ 0 m t in R 3 T, divh + m = 0 in R 3 T, dive = 0 in DT, were m is te zero extension of m to R 3 and H eff is te effective field defined by H eff = C e m + H for some constant C e > 0. Here te parameter α > 0 and permeability µ 0 0 are constants, wereas te conductivity σ takes a constant positive value in D and te zero value in D. Equation 1d is understood in te distributional sense because tere is a jump of m across Γ. Note tat H eff contains only te ig order term for simplicity. A refined analysis migt also allow us to include lower order terms anisotropy, exterior applied field as done in [14].

3 THE EDDY CURRENT LLG EQUATIONS 3 It follows from 1a tat m is constant. We follow te usual practice to normalize m and tus te same condition is required for m 0. Te following conditions are imposed on te solutions of 1: 2a 2b 2c 2d 2e n m = 0 on Γ T, m = 1 in D T, m0, = m 0 in D, H0, = H 0 in R 3, Ht, x = O x 1 as x, were n denotes te normal derivative. Te initial data m 0 and H 0 satisfy m 0 = 1 in D and 3 divh 0 + m 0 = 0 in R 3. Te condition 2b togeter wit basic properties of te cross product leads to te following equivalent formulation of 1a: 4 αm t + m m t = H eff m H eff m in D T. Below, we focus on an H-based formulation of te problem. recover E once H and m are known; see 12 It is possible to 2.2. Function spaces and notations. Before introducing te concept of weak solutions to problem 1 2 we need te following definitions of function spaces. Let L 2 D := L 2 D; R 3 and Hcurl, D := { w L 2 D : w L 2 D }. We define H 1/2 Γ as te usual trace space of H 1 D and define its dual space H 1/2 Γ by extending te L 2 -inner product on Γ. For convenience we denote X := { ξ, ζ Hcurl, D H 1/2 Γ : n ξ Γ = n Γ ζ in te sense of traces }. Recall tat n ξ Γ is te tangential trace or twisted tangential trace of ξ, and Γ ζ is te surface gradient of ζ. Teir definitions and properties can be found in [15, 16]. Finally, if X is a normed vector space ten, for m 0 and p N { }, L 2 0, T ; X, H m 0, T ; X, and W m,p 0, T ; X denote te usual Lebesgue and Sobolev spaces of functions defined on 0, T and taking values in X. We finis tis subsection wit te clarification of te meaning of te cross product between different matematical objects. For any vector functions u, v, w we denote u v := u v, u v, u v, u v := x 1 x 2 x 3 3 i=1 u v x i x i and u v w := 3 i=1 u v w. x i x i 2.3. Weak solutions. A weak formulation for 1a is well-known, see e.g. [3, 26]. Indeed, by multiplying 4 by φ C D T ; R 3, using integration by parts, we deduce α m t, m φ DT + m m t, m φ DT +C e m, m φ DT = H, m φ DT.

4 4 M. FEISCHL AND T. TRAN To tackle te eddy current equations on R 3, we aim to employ FE/BE coupling metods. To tat end, we employ te magnetic approac from [13], wic eventually results in a variant of te Trifou-discretization of te eddy-current Maxwell equations. Multiplying 1c by ξ C D, R 3 satisfying ξ = 0 in D, integrating over R 3, and using integration by parts, we obtain for almost all t [0, T ] µ 0 H t t, ξ R 3 + Et, ξ R 3 = µ 0 m t t, ξ D. Using ξ = 0 in D and 1b we deduce µ 0 H t t, ξ R 3 + σ 1 Ht, ξ D = µ 0 m t t, ξ D. Since H = ξ = 0 in D and D is simply connected by definition a workaround for non-simply connected D is presented in [22], tere exists ϕ and ζ suc tat H = ϕ and ξ = ζ in D. Terefore, te above equation can be rewritten as µ 0 H t t, ξ D + µ 0 ϕ t t, ζ D + σ 1 Ht, ξ D = µ 0 m t t, ξ D. Since 1d implies divh = 0 in D, we ave ϕ = 0 in D, so tat formally ϕ t = 0 in D. Hence integration by parts yields 5 µ 0 H t t, ξ D µ 0 + n ϕ t t, ζ Γ + σ 1 Ht, ξ D = µ 0 m t t, ξ D, were n + is te exterior Neumann trace operator wit te limit taken from D. Te advantage of te above formulation is tat no integration over te unbounded domain D is required. Te exterior Neumann trace n + ϕ t can be computed from te exterior Diriclet trace λ of ϕ by using te Diriclet-to-Neumann operator S, wic is defined as follows. Let γ be te interior Diriclet trace operator and n be te interior normal derivative or Neumann trace operator. Te sign indicates tat te trace is taken from D. Recalling te fundamental solution of te Laplacian Gx, y := 1/4π x y, we introduce te following integral operators defined formally on Γ as Vλ := γ Vλ, Kλ := γ Kλ + 1 2, and Wλ := n Kλ, were, for x / Γ, Vλx := Γ Gx, yλy ds y and Kλx := Γ ny Gx, yλy ds y, see, e.g., [29] for furter details. Moreover, let K denote te adjoint operator of K wit respect to te extended L 2 -inner product. Ten te exterior Diriclet-to-Neumann map S: H 1/2 Γ H 1/2 Γ can be represented as 6 S = V 1 1/2 K. Anoter representation is 7 S = 1/2 K V 1 1/2 K W. Recall tat ϕ satisfies H = ϕ in D. We can coose ϕ satisfying ϕx = O x 1 as x. Now if λ = γ + ϕ ten λ t = γ + ϕ t. Since ϕ = ϕ t = 0 in D,

5 THE EDDY CURRENT LLG EQUATIONS 5 and since te exterior Laplace problem as a unique solution we ave Sλ = + n ϕ and Sλ t = + n ϕ t. Hence 5 can be rewritten as 8 H t t, ξ D Sλ t t, ζ Γ + µ 1 0 σ 1 Ht, ξ D = m t t, ξ D. We remark tat if Γ denotes te surface gradient operator on Γ ten it is wellknown tat Γ λ = ϕ Γ + n ϕn = H Γ + n ϕn; see e.g. [30, Section 3.4]. Hence n Γ λ = n H Γ. Te above analysis prompts us to define te following weak formulation. Definition 1. A triple m, H, λ satisfying m H 1 D T and m t ΓT L 2 0, T ; H 1/2 Γ, H L 2 0, T ; Hcurl, D H 1 0, T ; L 2 D, λ H 1 0, T ; H 1/2 Γ is called a weak solution to 1 2 if te following statements old 1. m = 1 almost everywere in D T ; 2. m0, = m 0, H0, = H 0, and λ0, = γ + ϕ 0 were ϕ 0 is a scalar function satisfies H 0 = ϕ 0 in D te assumption 3 ensures te existence of ϕ 0 ; 3. For all φ C D T ; R 3 9a α m t, m φ DT + m m t, m φ DT + C e m, m φ DT = H, m φ DT ; 4. Tere olds n Γ λ = n H Γ in te sense of traces; 5. For ξ C D; R 3 and ζ C Γ satisfying n ξ Γ = n Γ ζ in te sense of traces 9b H t, ξ DT Sλ t, ζ ΓT + σ 1 µ 1 0 H, ξ D T = m t, ξ DT ; 6. For almost all t [0, T ] mt 2 L 2 D + Ht 2 Hcurl,D + λt 2 H 1/2 Γ 10 + m t 2 L 2 D t + H t 2 L 2 D t + λ t 2 H 1/2 Γ t C, were te constant C > 0 is independent of t. A triple m, H, λ is called a strong solution of te ELLG system 1 2 if it is a weak solution and additionally it is sufficiently smoot suc tat 4 is satisfied in te strong sense. Remark 2. A refinement of te arguments in Teorem 5 would allow us to prove tat te weak solutions wic appear as limits of te approximations from Algoritm 2.5, are energy dissipative, i.e., C e 2 mt 2 L 2 D + Ht 2 L 2 D Sλt, λt Γ C e 2 m0 2 L 2 D + H0 2 L 2 D Sλ0, λ0 Γ for all t [0, T ]. Te proof works along te lines of [1, Teorem 24] or [14, Appendix A] and is terefore omitted.

6 6 M. FEISCHL AND T. TRAN Te reason we integrate over [0, T ] in 8 to ave 9b is to facilitate te passing to te limit in te proof of te main teorem. Te following lemma justifies te above definition. Lemma 3. Let m, H, E be a strong solution of 1 2. If ϕ H0, T ; H 1 D satisfies ϕ = H D T, and if λ := γ + ϕ, ten te triple m, H DT, λ is a weak solution in te sense of Definition 1. Conversely, given a weak solution m, H DT, λ in te sense of Definition 1, let ϕ be te solution of 11 ϕ = 0 in D, ϕ = λ on Γ, ϕx = O x 1 as x and define H D T := φ as well as E via E = σ 1 H DT in D T and outside of D T as te solution of 12a 12b 12c E = µ 0 H t in D T, dive = 0 in D T, n E D T = n E DT on Γ T. If m, H, and E are sufficiently smoot, m, H, E is a strong solution of 1 2. Proof. We follow [13]. Assume tat m, H, E satisfies 1 2. Ten, Item 1, Item 2 and Item 6 in Definition 1 old, noting 3. Item 3, Item 4 and Item 5 also old due to te analysis above Definition 1. Te converse is also true due to te well-posedness of 12 as stated in [13, Equation 15]. Remark 4. Te solution ϕ to 11 can be represented as ϕ = 1/2 + Kλ VSλ. Te next subsection defines te spaces and functions to be used in te approximation of te weak solution te sense of Definition Discrete spaces and functions. For time discretization, we use a uniform partition 0 t i T, i = 0,..., N wit t i := ik and k := T/N. Te spatial discretization is determined by a sape regular triangulation T of D into compact tetraedra τ T wit diameter τ /C C τ 1/3 for some uniform constant C > 0. Denoting by N te set of nodes of T, we define te following spaces S 1 T := { φ CD : φ τ P 1 τ for all τ T }, K φ := { ψ S 1 T 3 : ψ z φ z = 0 for all z N }, φ S 1 T 3, were P 1 τ is te space of polynomials of degree at most 1 on τ. For te discretization of 9b, we employ te space N D 1 T of first order Nédélec edge elements for H and and te space S 1 T Γ for λ. Here T Γ denotes te restriction of te triangulation to te boundary Γ. It follows from Item 4 in Definition 1 tat for eac t [0, T ], te pair Ht, λt X. We approximate te space X by X := { ξ, ζ N D 1 T S 1 T Γ : n Γ ζ = n ξ Γ }. To ensure te condition n Γ ζ = n ξ Γ, we observe te following. For any ζ S 1 T Γ, if e denotes an edge of T on Γ, ten e ξ τ ds = e ζ τ ds = ζz 0 ζz 1, were τ is te unit direction vector on e, and z 0, z 1 are te endpoints of e. Tus, taking as degrees of freedom all interior edges of T i.e. e i ξ τ ds as well as all nodes of T Γ i.e. ζz i, we fully determine a function pair ξ, ζ X. Due to te considerations above, it is clear tat te above space can be implemented directly witout use of Lagrange multipliers or oter extra equations.

7 THE EDDY CURRENT LLG EQUATIONS 7 Te density properties of te finite element spaces {X } >0 are sown in Subsection 3.1; see Lemma 10. Given functions w i : D Rd, d N, for all i = 0,..., N we define for all t [t i, t i+1 ] w k t := t i+1 t k Moreover, we define w i + t t i k wi+1, w k t := wi, w + k t := wi d t w i+1 := wi+1 w i k for all i = 0,..., N 1. Finally, we denote by Π S te usual interpolation operator on S 1 T.We are now ready to present te algoritm to compute approximate solutions to problem Numerical algoritm. In te sequel, wen tere is no confusion we use te same notation H for te restriction of H : R 3 T R3 to te domain D T. Input: Initial data m 0 S1 T 3, H 0, λ 0 X, and parameter θ [0, 1]. For i = 0,..., N 1 do: 1. Compute te unique function v i K m i satisfying for all φ K m i 14 α v i, φ D + m i v i, φ D + C e θk v i, φ D = C e m i, φ D + H i, φ D. 2. Define m i+1 S 1 T 3 nodewise by 15 m i+1 z = mi z + kv i z for all z N. 3. Compute te unique functions H i+1 X, λi+1 X satisfying for all ξ, ζ d t H i+1, ξ D d t S λ i+1, ζ Γ + σ 1 µ 1 0 Hi+1, ξ D 16 = v i, ξ D, were S : H 1/2 Γ S 1 T Γ is te discrete Diriclet-to-Neumann operator to be defined later. Output: Approximations m i, Hi, λ i for all i = 0,..., N. Te linear formula 15 was introduced in [9] for armonic map eat flow and adapted for LLG in [1]. As already observed in [25, 8] for bounded domains, we note tat te linear systems 14 and 16 are decoupled and can be solved successively. Equation 16 requires te computation of S λ for any λ H 1/2 Γ. Tis is done by use of te boundary element metod. Let µ H 1/2 Γ and µ P 0 T Γ be, respectively, te solution of 17 Vµ = K 1/2λ and Vµ, ν Γ = K 1/2λ, ν Γ ν P 0 T Γ, were P 0 T Γ is te space of piecewise-constant functions on T Γ.

8 8 M. FEISCHL AND T. TRAN If te representation 6 of S is used, ten Sλ = µ, and we can uniquely define S λ by solving 18 S λ, ζ Γ = µ, ζ Γ ζ S 1 T Γ. Tis is known as te Jonson-Nédélec coupling. If we use te representation 7 for Sλ ten Sλ = 1/2 K µ Wλ. In tis case we can uniquely define S λ by solving 19 S λ, ζ Γ = 1/2 K µ, ζ Γ Wλ, ζ Γ ζ S 1 T Γ. Tis approac yields an almost symmetric system and is called symmetric coupling. In practice, 16 only requires te computation of S λ, ζ Γ for any λ, ζ S 1 T Γ. So in te implementation, neiter 18 nor 19 as to be solved. It suffices to solve te second equation in 17 and compute te rigt-and side of eiter 18 or 19. It is proved in [6, Appendix A] tat symmetric coupling results in a discrete operator wic is uniformly elliptic and continuous: 20 S ζ, ζ Γ C 1 S ζ 2 H 1/2 Γ for all ζ S 1 T Γ, S ζ 2 H 1/2 Γ C S ζ 2 H 1/2 Γ for all ζ H 1/2 Γ, for some constant C S > 0 wic depends only on Γ. Even toug we are convinced tat te proposed algoritm works for bot approaces, we are not aware of te essential ellipticity result of te form 20 for te Jonson-Nédélec approac. Tus, ereafter, S is understood to be defined by te symmetric coupling Main results. Before stating te main results, we first state some general assumptions. Firstly, te weak convergence of approximate solutions requires te following conditions on and k, depending on te value of te parameter θ in 14: k = o 2 wen 0 θ < 1/2, 21 k = o wen θ = 1/2, no condition wen 1/2 < θ 1. Some supporting lemmas wic ave teir own interests do not require any condition wen θ = 1/2. For tose results, a sligtly different condition is required, namely { k = o 2 wen 0 θ < 1/2, 22 no condition wen 1/2 θ 1. Te initial data are assumed to satisfy m 0 H1 D + H 0 Hcurl,D + λ 0 H 1/2 Γ <, 23 sup >0 lim 0 m0 m 0 L 2 D = 0. Te following tree teorems state te main results of tis paper. Te first teorem proves existence of weak solutions. Te second teorem establises a priori error estimates for te pure LLG case of 1a, i.e., H eff = C e m and tere is no coupling wit te eddy current equations 1b 1d. Te tird teorem provides a priori error estimates for te ELLG system.

9 THE EDDY CURRENT LLG EQUATIONS 9 Teorem 5 Existence of solutions. Under te assumptions 21 and 23, te problem 1 2 as a solution m, H, λ in te sense of Definition 1. Teorem 6 Error estimates for LLG. Let m W 2, 0, T ; H 1 D W 1, 0, T ; W 1, D H 2 D denote a strong solution of 1a and 2a 2c wit H eff = C e m. Ten for θ > 1/2 were θ is te parameter in 14 and for all, k satisfying 0 <, k 1 and k α/2c e, te following statements old max mt i m i H 0 i N 1 D C conv m 0 m 0 H 1 D + + k 24 and 25 m m k L 2 0,T ;H 1 D C conv m 0 m 0 H 1 D + + k. Te constant C conv > 0 depends only on te regularity of m, te sape regularity of T, and te values of α and θ. Moreover, te strong solution m is unique and coincides wit te weak solution from Teorem 5. Teorem 7 Error estimates for for ELLG. Let m, H, λ be a strong solution of ELLG in te sense of Definition 1 wit te following properties m W 2, 0, T ; H 1 D W 1, 0, T ; W 1, D H 2 D, H W 2, 0, T ; H 2 D L D T, λ W 1, 0, T ; H 1 Γ suc tat Sλ t L 0, T ; H 1/2 pw Γ, were Hpw 1/2 Γ is defined piecewise on eac smoot part of Γ. Ten for θ > 1/2 and for all, k satisfying 0 <, k 1/2 and k α/2c e, tere old 26 and max 0 i N mt i m i 2 H 1 D + Ht i H i 2 L 2 D + λt i λ i 2 H 1/2 Γ + k Ht i H i 2 L 2 D C conv m 0 m 0 2 H 1 D + λ0 λ 0 2 H 1/2 Γ + k H 0 H 0 2 L 2 D k 2 m m k 2 L 2 0,T ;H 1 D + H H k, λ λ k 2 L 2 0,T ;X 27 C conv m 0 m 0 2 H 1 D + H0 H 0 2 L 2 D + k H 0 H 0, λ 0 λ 0 2 X k 2. Particularly, te strong solution m, H, λ is unique and coincides wit te weak solution from Teorem 5. Te constant C conv > 0 depends only on te smootness of m, H, λ, and on te sape regularity of T. Remark 8. It is possible to replace te assumption mt W 1, D in Teorems 6 7 by mt L 4 D and mt L D. Tis, owever, results in a reduced rate of convergence k instead of k; see Remark 22 for furter discussion.

10 10 M. FEISCHL AND T. TRAN Remark 9. Little is known about te regularity of te solutions of 1 in 3D; see [18]. However, for te 2D case of D being te flat torus R 2 /Z 2, Teorem 5.2 in [31] states te existence of arbitrarily smoot local solutions of te full MLLG system, given tat m 0 L2 D and te initial values of te Maxwell system are sufficiently small. Since te eddy-current equations are a particular simplification of Maxwell s equations, tis strongly endorses te assumption tat also for te 3D case of te ELLG equations tere exist arbitrarily smoot local solutions. 3. Proof of te main results Some lemmas. In tis subsection we prove all important lemmas wic are directly related to te proofs of te main results. Te first lemma proves density properties of te discrete spaces. Lemma 10. Provided tat te meses {T } >0 are regular, te union >0 X is dense in X. Tere exists an interpolation operator Π X := Π X,D, Π X,Γ : H 2 D H 2 Γ X X wic satisfies 28 1 Π X ξ, ζ Hcurl,D H 1/2 Γ C X ξ H2 D + 1/2 ζ H2 Γ, were C X > 0 depends only on D, Γ, and te sape regularity of T. Proof. Te interpolation operator Π X := Π X,D, Π X,Γ : H 2 D H 2 Γ X X is constructed as follows. Te interior degrees of freedom edges of Π X ξ, ζ are equal to te interior degrees of freedom of Π N D ξ N D 1 T, were Π N D is te usual interpolation operator onto N D 1 T. Te degrees of freedom of Π X ξ, ζ wic lie on Γ nodes are equal to Π S ζ. By te definition of X, tis fully determines Π X. Particularly, since n ξ Γ = n Γ ζ, tere olds Π N D ξ Γ = Π X,Γ ξ, ζ. Hence, te interpolation error can be bounded by 1 Π X ξ, ζ Hcurl,D H 1/2 Γ 1 Π N D ξ Hcurl,D + 1 Π S ζ H 1/2 Γ ξ H2 D + 1/2 ζ H2 Γ. Since H 2 D H 2 Γ X is dense in X, tis concludes te proof. Te following lemma gives an equivalent form to 9b and sows tat Algoritm 2.5 is well-defined. Lemma 11. Let a, : X X R, a, : X X R, and b, : Hcurl, D Hcurl, D R be bilinear forms defined by aa, B := ψ, ξ D Sη, ζ Γ, a A, B := ψ, ξ D S η, ζ Γ, bψ, ξ := σ 1 µ 1 0 ψ, ξ D, for all ψ, ξ Hcurl, D, A := ψ, η, B := ξ, ζ X, A = ψ, η, B = ξ, ζ X. Ten 1. Te bilinear forms satisfy, for all A = ψ, η X and A = ψ, η X, aa, A C ell ψ 2 L 2 D + η 2 H 1/2 Γ, 29 a A, A C ell ψ 2 L 2 D + η 2 H 1/2 Γ, bψ, ψ C ell ψ 2 L 2 D.

11 2. Equation 9b is equivalent to THE EDDY CURRENT LLG EQUATIONS T 0 T aa t t, B dt + bht, ξ dt = m t, ξ DT 0 for all B = ξ, ζ X, were A = H, λ. 3. Equation 16 is of te form 31 a d t A i+1, B + bh i+1, ξ = v i, ξ Γ were A i+1 := H i+1, λi+1 and B := ξ, ζ. 4. Algoritm 2.5 is well-defined in te sense tat 14 and 16 ave unique solutions. 5. Te norm 32 B 2 := ξ 2 L 2 D S ζ, ζ Γ B = ξ, ζ X, is equivalent to te grap norm of L 2 D H 1/2 Γ uniformly in. Proof. Te unique solvability of 16 follows immediately from te continuity and ellipticity of te bilinear forms a, and b,. Te unique solvability of 14 follows from te positive definiteness of te leftand side, te linearity of te rigt-and side, and te finite space dimension. Te ellipticity 20 sows te norm equivalence in te final statement. Te following lemma establises an energy bound for te discrete solutions. Lemma 12. Under te assumptions 22 and 23, tere olds for all k < 2α and j = 1,..., N 33 j 1 H i+1 H i 2 L 2 D + λi+1 λ i 2 H 1/2 Γ j 1 + k H i+1 2 L 2 D + Hj 2 Hcurl,D + λj 2 H 1/2 Γ + m j 2 L 2 D + max{2θ 1, j 1 j 1 0}k2 v i 2 L 2 D + k v i 2 L 2 D j 1 + k d t H i+1 2 L 2 D + d tλ i+1 j 1 + H i+1 2 H 1/2 Γ H i 2 L 2 D C ener. Proof. Coosing B = A i+1 we obtain in 31 and multiplying te resulting equation by k 34 a A i+1 A i, A i+1 +kbhi+1, Hi+1 = k vi, H i D k v i, H i+1 H i D. Following te lines of [25, Lemma 5.2] using te definition of a, and 20, we end

12 12 M. FEISCHL AND T. TRAN up wit a A i+1 A i, A i+1 + kbhi+1, Hi+1 + C e 2 m i+1 2 L 2 D mi 2 L 2 D + θ 1/2k 2 C e v i 2 L 2 D + α ɛ/2k vi 2 L 2 D k 2ɛ a A i+1 A i, A i+1 A i. Summing over i from 0 to j 1 and for te first term on te left-and side applying Abel s summation by parts formula we derive as in [25, Lemma 5.2] 35 j 1 H i+1 H i 2 L 2 D + λi+1 λ i 2 H 1/2 Γ j 1 + k H i+1 2 L 2 D + Hj 2 L 2 D + λ j 2 H 1/2 Γ + mj 2 L 2 D j 1 j 1 + 2θ 1k 2 v i 2 L 2 D + k v i 2 L 2 D C m 0 2 L 2 D + H0 2 L 2 D + λ0 2 H 1/2 Γ C, were in te last step we used 23. It remains to consider te last tree terms on te left-and side of 33. Again, we consider 31 and select B = d t A i+1 to obtain after multiplication by 2k 2ka d t A i+1, d ta i+1 + 2bHi+1, Hi+1 H i = 2k v i, d t H i+1 D k v i 2 L 2 D + k d th i+1 2 L 2 D, so tat, noting 35 and 29, 36 j 1 k d t H i+1 2 L 2 D + d tλ i+1 2 H 1/2 Γ j 1 k v i 2 L 2 D C. j bh i+1, Hi+1 H i Using Abel s summation by parts formula as in [25, Lemma 5.2] for te second sum on te left-and side, and noting te ellipticity of te bilinear form b, and 23,

13 we obtain togeter wit 35 j 1 H i+1 THE EDDY CURRENT LLG EQUATIONS 13 H i 2 L 2 D + λi+1 λ i 2 H 1/2 Γ j 1 + k H i+1 2 L 2 D + Hj 2 Hcurl,D + λj 2 H 1/2 Γ + mj 2 L 2 D 37 j 1 j 1 + 2θ 1k 2 v i 2 L 2 D + k v i 2 L 2 D j 1 + k d t H i+1 2 L 2 D + d tλ i+1 j 1 2 H 1/2 Γ + H i+1 H i 2 L 2 D C. Clearly, if 1/2 θ 1 ten 37 yields 33. If 0 θ < 1/2, we argue as in [25, Remark 6] to conclude te proof. Collecting te above results we obtain te following equations satisfied by te discrete functions defined from m i, Hi, λ i, and vi. Lemma 13. Let m k, A± k := H± k, λ± k, and v k be defined from mi, Hi, λ i, and v i as described in Subsection 2.4. Ten 38a α v k, φ k DT + m k v k, φ k DT + C e θk v k, φ k DT = C e m k, φ k DT + H k, φ k DT and wit t denoting time derivative 38b T 0 T a t A k t, B dt + bh + k t, ξ dt = v k, ξ DT 0 for all φ k and B := ξ, ζ satisfying φ k t, K m i for t [t i, t i+1 and B X. Proof. Te lemma is a direct consequence of 14 and 31. In te following lemma, we state some auxiliary results, already proved in [1]. Lemma 14. Tere olds m j 2 L 2 D m L 2 D CkC ener, as well as 40 m ± k L 2 0,T ;H 1 D T m ± k L 0,T ;H 1 D CC ener, were C > 0 depends only on te sape regularity of T and T. Te esti- Proof. Te estimate 39 follows analogously to [1, Proposition 9]. mate 40 ten follows from 23, 33, and 39. Te next lemma sows tat te functions defined in te above lemma form sequences wic ave convergent subsequences.

14 14 M. FEISCHL AND T. TRAN Lemma 15. Assume tat te assumptions 22 and 23 old. As, k 0, te following limits exist up to extraction of subsequences all limits old for te same subsequence 41a 41b 41c 41d 41e 41f 41g m k m in H 1 D T, m ± k m in L2 0, T ; H 1 D, m ± k m in L2 D T, H k, λ k H, λ in L 2 0, T ; X, H ± k, λ± k H, λ in L2 0, T ; X, H k, λ k H, λ in H 1 0, T ; L 2 D H 1/2 Γ, v k m t in L 2 D T, for certain functions m, H, and λ satisfying m H 1 D T, H H 1 0, T ; L 2 D, and H, λ L 2 0, T ; X. Here denotes te weak convergence and denotes te strong convergence in te relevant space. Moreover, if te assumption 23 olds ten tere olds additionally m = 1 almost everywere in D T. Proof. Te proof works analogously to [1] and is terefore omitted. We also need te following strong convergence property. Lemma 16. Under te assumptions 21 and 23 tere olds 42 m k m L 2 0,T ;H 1/2 D 0 as, k 0. Proof. It follows from te triangle inequality and te definitions of m k and m k tat m k m 2 L 2 0,T ;H 1/2 D m k m k 2 L 2 0,T ;H 1/2 D + m k m 2 L 2 0,T ;H 1/2 D N 1 N 1 k 3 v i 2 H 1/2 D + m k m 2 L 2 0,T ;H 1/2 D k 3 v i 2 H 1 D + m k m 2 L 2 0,T ;H 1/2 D. Te second term on te rigt-and side converges to zero due to 41a and te compact embedding of H 1 D T {v v L 2 0, T ; H 1 D, v t L 2 0, T ; L 2 D} into L 2 0, T ; H 1/2 D; see [27, Teorem 5.1]. For te first term on te rigt-and side, wen θ > 1/2, 33 implies N 1 k3 v i 2 H 1 D k 0. Wen 0 θ 1/2, a standard inverse inequality, 33 and 21 yield N 1 N 1 k 3 v i 2 H 1 D completing te proof of te lemma. 2 k 3 v i 2 L 2 D 2 k 2 0,

15 THE EDDY CURRENT LLG EQUATIONS 15 Te following lemma involving te L 2 -norm of te cross product of two vectorvalued functions will be used wen passing to te limit of equation 38a. Lemma 17. Tere exists a constant C sob > 0 wic depends only on D suc tat 43a 43b w 1/2 w 1 L2 D C sob w 1/2 H 1/2 D w 1 H1 D w 0 w 1/2 H 1 D C sob w 0 L 2 D w 1/2 H 1/2 D, for all w 0 L 2 D, w 1/2 H 1/2 D, and all w 1 H 1 D. Proof. It is sown in [2, Teorem 5.4, Part I] tat te embedding ι: H 1 D L 6 D is continuous. Obviously, te identity ι: L 2 D L 2 D is continuous. By real interpolation, we find tat ι: [L 2 D, H 1 D] 1/2 [L 2 D, L 6 D] 1/2 is continuous. Well-known results in interpolation teory sow [L 2 D, H 1 D] 1/2 = H 1/2 D and [L 2 D, L 6 D] 1/2 = L 3 D wit equivalent norms; see e.g. [12, Teorem 5.2.1]. By using Hölder s inequality, we deduce w 0 w 1 L 2 D w 0 L 3 D w 1 L 6 D w 0 H 1/2 D w 1 H 1 D proving 43a. For te second statement, tere olds wit te well-known identity 44 a b c = b c a = c a b a, b, c R 3 tat w 0 w 1/2 H 1 D = sup w 0 w 1/2, w 1 D w 1 H 1 D\{0} w 1 H1 D sup w 1 H 1 D\{0} Te estimate 43a concludes te proof. w 0 L2 D w 1/2 w 1 L2 D. w 1 H 1 D Finally, to pass to te limit in equation 38b we need te following result. Lemma 18. For any sequence {λ } H 1/2 Γ and any function λ H 1/2 Γ, if 45 lim 0 λ, ν Γ = λ, ν Γ ν H 1/2 Γ ten 46 lim 0 S λ, ζ Γ = Sλ, ζ Γ ζ H 1/2 Γ. Proof. Let µ and µ be defined by 17 wit λ in te second equation replaced by λ. Ten recalling tat Costabel s symmetric coupling is used Sλ and S λ are defined via µ and µ by 7 and 19, respectively, namely, Sλ = 1/2 K µ Wλ and S λ, ζ Γ = 1/2 K µ, ζ Γ Wλ, ζ Γ for all ζ S 1 T Γ. For any ζ H 1/2 Γ, let {ζ } be a sequence in S 1 T Γ satisfying lim 0 ζ ζ H 1/2 Γ = 0. By using te triangle inequality and te above representations of Sλ and S λ we

16 16 M. FEISCHL AND T. TRAN deduce S λ, ζ Sλ, ζ Γ S λ Sλ, ζ Γ + S λ Sλ, ζ ζ Γ K µ µ, ζ Γ + Wλ λ, ζ Γ + S λ Sλ, ζ ζ Γ 1 2 K µ µ, ζ Γ + Wλ λ, ζ Γ + Wλ λ, ζ ζ Γ + S λ Sλ, ζ ζ Γ. Te second term on te rigt-and side of 47 goes to zero as 0 due to 45 and te self-adjointness of W. Te tird term converges to zero due to te strong convergence ζ ζ in H 1/2 Γ and te boundedness of {λ } in H 1/2 Γ, wic is a consequence of 45 and te Banac-Steinaus Teorem. Te last term tends to zero due to te convergence of {ζ } and te boundedness of {S λ }; see 20. Hence 46 is proved if we prove 48 lim 0 1/2 K µ µ, ζ Γ = 0. Tis, owever, follows from standard convergence arguments in boundary element metods and concludes te proof Furter lemmas for te proofs of Teorems 6 and 7. In tis section, we prove all necessary results for te proof of te a priori error estimates. Lemma 19. Recall te operators S and S defined in 7 and 19. Given λ suc tat Sλ H 1/2 Γ, tere olds 49 Sλ S λ, ξ Γ sup C S Sλ 1/2 ξ S 1 T Γ\{0} ξ H H 1/2 pw Γ. Γ Proof. By definition of µ in 17 as te Galerkin approximation of µ = Sλ, tere olds by standard arguments µ µ H 1/2 Γ µ H 1/2 pw Γ. Hence, tere olds wit te mapping properties of K Sλ S λ, ξ Γ = 1/2 K µ µ, ξ Γ µ H 1/2 pw Γ ξ H 1/2 Γ. Tis concludes te proof. Te next lemma proves tat te time derivative of te exact solution, namely m t t i, can be approximated in te discrete tangent space K m i in suc a way tat te error can be controlled by te error in te approximation of mt i by m i. Lemma 20. Assume te following regularity of te strong solution m of 1: C reg := m W 1, 0,T ;H 2 D + m W 1, 0,T ;W 1, D <. For any i = 1,..., N let P i : H1 D K m i K m i. Ten denote te ortogonal projection onto m t t i P i m t t i H 1 D C P + mti m i H 1 D, were C P > 0 depends only on C reg and te sape regularity of T.

17 THE EDDY CURRENT LLG EQUATIONS 17 Proof. We fix i {1,..., N}. Due to te well-known result m t t i P i m t t i H 1 D = inf w K m i m t t i w H 1 D and te estimate recalling te definition of te interpolation Π S in Subsection 2.4 m t t i w H1 D m t t i Π S m t t i H1 D + Π S m t t i w H1 D C reg + Π S m t t i w H1 D w K m i, it suffices to prove tat 50 inf w K m i Π S m t t i w H1 D + mt i m i H1 D, were te constant depends only on C reg and te sape regularity of T. To tis end we first note tat te assumption on te regularity of te exact solution m implies tat after a modification on a set of measure zero, we can assume tat m is continuous in D T. Hence mt i, x = 1 and m t t i, x mt i, x = 0 x D. Tus by using te elementary identity it can be easily sown tat a b c = a cb a bc a, b, c R 3 51 Π S m t t i, z = m t t i, z = R i zmt i, z were, for any z N, te mapping R i z : R 3 R 3 is defined by R i za = a mt i, z m t t i, z a R 3. We note tat tis mapping as te following properties 52 R i za a = 0 and R i za m t t i, z a a R 3. Next, prompted by 51 and 52 in order to prove 50 we define w S 1 T by 53 wz = R i zm i z z N. Ten w K m i and we can estimate Π S m t t i w 2 H 1 D by Π S m t t i w 2 H 1 D = T T Π S m t t i w 2 H 1 T T T Π S m t t i w w T 2 H 1 T + were w T is polynomial of degree 1 on T defined by 54 w T z = R i z mti, z T m i z T z N T. Here, z T is te node in T satisfying 55 mt i, z T m i z T mt i, z m i z z N T. T T w T 2 H 1 T

18 18 M. FEISCHL AND T. TRAN Tus, 50 is proved if we prove 56 Π S m t t i w w T 2 H 1 T + mt i m i H 1 D T T and 57 T T w T 2 H 1 T mt i m i H 1 D. To prove 56 we denote A := T T Π S m t t i w w T 2 H 1 T standard inverse estimate and te equivalence [26, Lemma 3.2] to ave A 2 Π S m t t i w w T 2 L 2 T T T ΠS m t t i, z wz w T z 2. T T z N T and use a Tis togeter wit and te regularity assumption of m yields A T T z N T T T z N T = T T z N T Rz i mt i, z m i z mt i, z T + m i 2 z T mti, z m i z mt i, z T + m i z T 2 Π S mt i, z m i z Π S mt i, z T m i z T 2. Since Π S m m i is polynomial of degree 1 on T, Lemma 28 in te Appendix gives A T T ΠS mt i m i 2 L 2 T = Π S mt i m i 2 L 2 D Π S mt i m i 2 H 1 D. By using te triangle inequality and te approximation property of te interpolation operator Π S, noting tat m L 0, T ; H 2 D, we obtain 56. It remains to prove 57. Denoting ϕz = mt i, z m i z for z N T, it follows successively from [26, Lemma 3.2], Lemma 28, 54 and 52 tat 58 w T 2 H 1 T 3 z N T = 3 3 z N T z N T w T z 2 + z N T R i zϕz T 2 + ϕz T 2 + z N T z N T w T z w T z T 2 R i z R i z T ϕzt 2 R i z R i z T ϕzt 2. For te term in te last sum on te rigt-and side, we use te triangle inequality

19 THE EDDY CURRENT LLG EQUATIONS 19 and te regularity of m to obtain R i z R i z T ϕzt 2 ϕz T 2 mt i, z m t t i, z mt i, z T m t t i, z T 2 ϕz T 2 mt i, z 2 m t t i, z m t t i, z T 2 + ϕz T 2 m t t i, z T 2 mt i, z mt i, z T 2 ϕz T 2 m t t i, z m t t i, z T 2 + mt i, z mt i, z T 2 2 ϕz T 2, were in te last step we used also Taylor s Teorem and z z T, noting tat m L D T and m t L D T. Terefore, 58, 55 and [26, Lemma 3.2] imply w T 2 H 1 T 3 ϕz T 2 3 ϕz 2 mt i m i 2 L 2 T. z N T z N T Summing over T T we obtain 57, completing te proof of te lemma. Te following lemma sows tat for te pure LLG equation, m t t i solves te same equation as v i, up to an error term. Lemma 21. Let m denote a strong solution of 1a, 2a 2c wic satisfies m L D T and m t L 0, T ; H 2 D. Ten, for i = 0,..., N tere olds 59 were α m t t i, φ D + mt i m t t i, φ D + C e θk m t t i, φ D = C e mt i, φ D + Rφ φ K m i, 60 Rφ C R + θk + mti m i L2 D φ L2 D. Here, te constant C R > 0 depends only on te regularity assumptions on m and te sape regularity of T. Proof. Note tat m = 1 implies m m = m 2. Recalling tat 1a is equivalent to 4, tis identity and 4 give, for all φ K m i, 61 α m t t i, φ D + mt i m t t i, φ D = C e mt i, φ D + C e mt i 2 mt i, φ D = C e mt i, φ D + C e mt i 2 mt i, φ D were in te last step we used 2a and integration by parts. Hence 59 olds wit 62 Rφ := C e θk m t t i, φ D + C e mt i 2 mt i, φ D. It remains to sow 60. Te condition 2a implies n m t = 0 on Γ T, and tus te first term on te rigt-and side of 62 can be estimated as 63 C e θk m t t i, φ D θk m t t i L2 D φ L2 D θk φ L2 D, were in te last step we used te fact tat m t L 0, T ; H 2 D. Te second term on te rigt-and side of 62 can be estimated as mti 2 mt i, φ D m 2 L D T mt i φ L1 D 64 mt i m i φ L 1 D + m i φ L 1 D. mt i m i L 2 D φ L 2 D + m i φ L 1 D.

20 20 M. FEISCHL AND T. TRAN Since m i φ is a quadratic function on eac T T wic vanises at every node in T, te semi-norm D 2 m i φ L1 T is a norm, were D 2 is te partial derivative operator of order 2. If T is te reference element, ten a scaling argument and norm equivalence on finite dimensional spaces give m i φ L1 D = m i φ L1 T 3 m i φ L1 T T T T T 3 D 2 m i φ L 1 T 65 2 D 2 m i φ L 1 T. T T T T Let i, i = 1, 2, 3, denote te directional derivatives in R 3. Since m i and φ are polynomials of degree 1 on T, tere olds i j m i φ = i j m i φ + m i j φ = j m i i φ + i m i j φ 2 m i φ, implying D 2 m i φ m i φ. Tis and 65 yield m i φ L1 D 2 66 T T m i φ L1 T = 2 m i φ L1 D 2 m i L2 D φ L2 D φ L2 D, were in te last step we used te energy bound 33, and a standard inverse estimate. Estimate 60 now follows from and 66, completing te proof of te lemma. Remark 22. It is noted tat te assumption m L D T can be replaced by m L 0, T ; L 4 D to obtain a weaker bound Rφ θk φ L2 D + + mt i m i H1 D φ H1 D. Tis can be done by use of te continuous embedding H 1 D L 4 D in 64 and obvious modifications in te remainder of te proof. Wit straigtforward modifications, te proof of Lemma 23 is still valid to prove a weaker estimate wit k instead of k 2 on te rigt-and side of 67. Tis eventually results in a reduced rate of convergence k in Teorem 6. Analogous arguments old true for te corresponding results in Teorem 7. As a consequence of te above lemma, we can estimate te approximation of m t t i by v i as follows. Lemma 23. Under te assumptions of Lemmas 20 and 21 tere olds, wit 1/2 < θ 1, α m t t i v i C 2 L 2 D + k m tt i v i 2 L 2 D e mt i m i, m t t i v i Proof. Subtracting 14 from 59 we obtain C m 2 + k 2 + mt i m i 2 H 1 D α m t t i v i, φ D + C e θk m t t i v i, φ D +C e mt i m i, φ D = m i v i mt i m t t i, φ D + Rφ, φ K m i..

21 THE EDDY CURRENT LLG EQUATIONS 21 Using te above equation, writing φ := m t t i v i side of 67 we try to estimate and referring to te left-and 68 α φ 2 L 2 D + C eθk φ 2 L 2 D + C e mt i m i, φ = α φ 2 L 2 D φ, φ D + Ce θk φ 2 L 2 D φ, φ D + C e mt i m i, φ φ + m i v i mt i m t t i, φ D + Rφ =: T T 5. Coosing φ = P i m tt i v i K m i from tat lemma tat were Pi is defined in Lemma 20, we deduce T 1 φ L 2 D φ φ L 2 D φ L 2 D + mti m i H 1 D, 69 T 2 k φ L 2 D φ φ L 2 D k φ L2 D + mti m i H 1 D, T 3 mt i m i L2 D φ φ L2 D mt i m i L2 D + mti m i H 1 D 2 + mt i m i 2 H 1 D. For te term T 4 since m i φ, φ D = 0 implies m i vi, φ D = m i m tt i, φ D we deduce, wit te elp of Lemma 20 again, T 4 m i v i mt i m t t i, φ φ D + m i v i mt i m t t i, φ D = m i v i mt i m t t i, φ φ D + m i mt i m t t i, φ D m i v i mt i m t t i H 1 D φ φ H 1 D + m i mt i L2 D φ L2 D m i v i mt i m t t i D H 1 + mti m i H 1 D + m i mt i L2 D φ L2 D. Lemma 17 implies m i v i mt i m t t i H 1 D m i v i m t t i H 1 D + mi mt i m t t i H 1 D m i H1 D m t t i v i L 2 D + m t t i L 2 D mt i m i H 1 D φ L 2 D + mt i m i H 1 D, were in te last step we used te regularity of m t and te bound 39 and 33. Terefore T 4 φ L2 D + mt i m i H1 D + mti m i H1 D 70 φ L 2 D + mti m i H 1 D mt i m i 2 H 1 D.

22 22 M. FEISCHL AND T. TRAN Finally, for T 5, Lemma 21, te triangle inequality, and Lemma 20 give 71 T 5 + k + mt i m i L2 D φ L 2 D + k + mt i m i L D 2 φ L 2 D + + k + mt i m i L D 2 + mti m i H 1 D + k + mt i m i L2 D φ L2 D k 2 + mt i m i 2 H 1 D. Altogeter, yield, for any ɛ > 0, α φ 2 L C 2 D + θk φ 2 L 2 D + mt i m i, φ e φ L2 D + k φ L2 D + k + mti m i H1 D k 2 + mt i m i 2 H 1 D ɛ φ 2 L 2 D + ɛk2 φ 2 L 2 D ɛ k 2 + mt i m i 2 H 1 D. Te required estimate 67 is obtained for ɛ = min{α/2c e, θ 1/2}. Tree more lemmas are required for te proof of Teorem 7. Analogously to Lemma 21 we sow tat for te ELLG system, m t t i solves te same equation as v i, up to an error term. Lemma 24. Let m, H, λ denote a strong solution of ELLG wic satisfies Ten, for i = 0,..., N tere olds m L 0, T ; W 1, D, m t L 0, T ; H 2 D, H L D T. α m t t i, φ D + mt i m t t i, φ D + C e θk m t t i, φ D 72 were = C e mt i, φ D + Ht i, φ D + Rφ φ K m i, 73 Rφ C R + θk + mti m i L2 D φ L2 D. Here, te constant C R > 0 depends only on te regularity of m and H, and te sape regularity of T. Proof. Te proof is similar to tat of Lemma 21. Instead of 61 we now ave and tus α m t t i, φ D + mt i m t t i, φ D = C e mt i, φ D + C e mt i 2 mt i, φ D + Ht i, φ D mt i Ht i mt i, φ D, Rφ := Rφ mt i Ht i mt i, φ D, were Rφ is given in 62. Terefore, it suffices to estimate te term mt i Ht i mt i, φ D. Since mt i Ht i mt i, φ D m L D T H L D T mt i φ L 1 D, te proof follows exactly te same way as tat of Lemma 21; cf. 64. prove 73. Tus we

23 THE EDDY CURRENT LLG EQUATIONS 23 We will establis a recurrence estimate for te eddy current part of te solution. We first introduce 74 e i := Π N D Ht i H i, f i := Π S λt i λ i, and E i := e i, f i, were Π N D denotes te usual interpolation operator onto N D 1 T. Lemma 25. Let m, H, λ be a strong solution of ELLG wic satisfies m W 1, 0, T ; L 2 D W 1, D T, H L 0, T ; H 2 D W 1, 0, T ; H 1 D W 2, 0, T ; L 2 D, λ W 1, 0, T ; H 1 Γ suc tat Sλ t L 0, T ; H 1/2 pw Γ. Ten for all k 0, 1/2] and i {1,..., N 1} tere olds 75 E i k 2 e i+1 2 L 2 D 1 + 2k E i 2 + C H k m t t i v i 2 L 2 D k 2 were te constant C H > 0 depends only on te smootness of H, λ, and on te sape-regularity of T. Proof. Recalling equation 31 and in view of 74, we will establis a similar equation for Π N D Ht i. Wit A := H, λ and A := Π N D H, Π S λ X, it follows from 9b tat a t A t, B + bπ N D Ht, ξ = m t t, ξ D a t A A t, B + b1 Π N D Ht, ξ + S S λ t t, ζ Γ =: m t t, ξ D + R 0 At, B B = ξ, ζ X. Te continuity of a, and b,, te approximation properties of Π N D and Π S, and 20, 49 give R 0 At, B 1 Π N D H t t L2 D ξ L2 D + 1 Π S λ t t H 1/2 Γ ζ H 1/2 Γ + 1 Π N D Ht Hcurl,D ξ L 2 D + Sλ t t 1/2 H pw Γ ζ H 1/2 Γ H t t H1 D + Ht H2 D + Sλ t t 1/2 B X H. pw Γ Te regularity assumptions on H and λ yield 76 R 0 At, B B X. Recalling definition 13 and using Taylor s Teorem, we ave d t A t i+1 = t A t i + r i, were r i is te remainder in te integral form of te Taylor expansion wic satisfies r i L 2 D H 1/2 Γ k H tt L 0,T ;L 2 D + λ tt L 0,T ;H 1/2 Γ. Terefore, Proof 16 and 76 imply for all B = ξ, ζ X 77 a d t A t i+1, B + bπ N D Ht i, ξ = m t t i, ξ D + R 1 At i, B,

24 24 M. FEISCHL AND T. TRAN were R 1 At i, B := R 0 At i, B a r i, B so tat 78 R 1 At i, B + k H tt L 0,T ;L 2 D + λ tt L 0,T ;H 1/2 Γ B X + k B X, were we used te uniform continuity of a, in obtained from 20. Subtracting 31 from 77 and setting B = E i+1 yield a d t E i+1, E i+1 + be i+1, e i+1 = v i m t t i, e i+1 D + R 1 At i, E i+1. Multiplying te above equation by k and using te ellipticity properties of a,, b,, and te Caucy-Scwarz inequality, and recalling definition 32 of te -norm, we deduce E i k e i+1 2 L 2 D E i E i+1 + k v i m t t i L 2 D e i+1 L 2 D + Ck + k E i+1 X, for some constant C > 0 wic does not depend on or k. Wit Young s inequality, tis implies 79 E i k e i+1 2 L 2 D 1 2 E i E i k C S 2 C S 1 vi m t t i 2 L 2 D + k C S 1 e i+1 2 L 2 C 2 D + k E i+1 2 X + kc S S 2C S 2 C2 + k 2. Note tat k CS 1 e i+1 2 L 2 C 2 D + 1 E i+1 2 X S C S = k ei+1 2 L 2 2 D + 1 e i+1 2 L C 2 D + 1 f i 2 H S C 1/2 Γ S k ei+1 2 L 2 2 D + 1 e i+1 2 L C 2 D S f i, f i Γ S k 2 E i k e i+1 2 L 2C 2 D. S Hence 79 yields after multiplying by 2 1 k E i k2 1 C S e i+1 2 L 2 D E i 2 + kc S C S 1 vi m t t i 2 L 2 D + kc SC 2 + k 2. Dividing by 1 k, using te fact tat 1 1/1 k 1 + 2k 2 since 0 < k 1/2, and noting tat C S 1, we obtain te desired estimate 75, concluding te proof. Similarly to Lemma 23 we now prove te following lemma for te ELLG system. Lemma 26. Let m, H, λ be a strong solution of ELLG wic satisfies m W 2, 0, T ; H 1 D W 1, 0, T ; W 1, D H 2 D, H L 0, T ; H 2 D L D W 2, 0, T ; L 2 D.

25 Ten, for 1/2 < θ 1, we ave THE EDDY CURRENT LLG EQUATIONS 25 α C e m t t i v i 2 L 2 D + k m tt i v i 2 L 2 D + 2 mt i m i, m t t i v i C H 2 + k 2 + mt i m i 2 H 1 D + Ht i H i 2 L 2 D. Proof. Te proof follows tat of Lemma 23. Subtracting 14 from 72 and putting φ := m t t i v i we obtain for φ := P i m tt i v i K m i α φ, φ D + C e θk φ, φ D + C e mt i m i, φ D = m i v i mt i m t t i, φ D + Ht i H i, φ D + Rφ. Similarly to 68 we now ave α φ 2 L 2 D + C eθk φ 2 L 2 D + C e mt i m i, φ D = T T 4 + T 5 + T 6, were T 1,..., T 4 are defined as in 68 wereas T 5 := Ht i H i, φ D and T6 := Rφ. Estimates for T 1,..., T 4 ave been carried out in te proof of Lemma 23. For T 5 we ave T 5 Ht i H i L 2 D φ L 2 D Ht i H i L 2 D φ L 2 D + Ht i H i L D 2 + mti m i L 2 D Ht i H i L2 D φ L2 D Ht i H i 2 L 2 D + mt i m i 2 L 2 D, were we used te triangle inequality and invoked Lemma 20 to estimate φ φ L 2 D. Finally, for T 6 we use 73, te triangle inequality, and Lemma 20 to obtain T 6 + k + mt i m i L2 D φ L 2 D + k + mt i m i L D 2 φ L 2 D + + k + mt i m i L D 2 + mti m i H 1 D 2 + k 2 + mt i m i 2 H 1 D + + k + mt i m i L2 D φ L2 D. Te proof finises in exactly te same manner as tat of Lemma Proof of Teorem 5. We are now ready to prove tat te problem 1 2 as a weak solution. Proof. We recall from 41a 41g tat m H 1 D T, H, λ L 2 0, T ; X and H H 1 0, T ; L 2 D. By virtue of Lemma 11 it suffices to prove tat m, H, λ satisfies 9a and 30. Let φ C D T and B := ξ, ζ L 2 0, T ; X. On te one and, we define te test function φ k := Π S m k φ as te usual interpolant of m k φ into S1 T 3.

26 26 M. FEISCHL AND T. TRAN By definition, φ k t, K m j for all t [t j, t j+1. On te oter and, it follows from Lemma 10 tat tere exists B := ξ, ζ X converging to B X. Equation 38 old wit tese test functions. Te main idea of te proof is to pass to te limit in 38a and 38b to obtain 9a and 30, respectively. In order to prove tat 38a implies 9a we need to prove tat as, k 0 80a 80b 80c 80d 80e v k, φ k DT m t, m φ DT, m k v k, φ k DT m m t, m φ DT, k v k, φ k DT 0, m k, φ k DT m, m φ DT, H k, φ k DT H, m φ DT. Te proof as been carried out in [1, 3, 25] and is terefore omitted. Next, recalling tat B B in X we prove tat 38b implies 30 by proving 81a 81b 81c 81d t H k, ξ DT H t, ξ DT, S t λ k, ζ ΓT Sλ t, ζ ΓT, H + k, ξ DT H, ξ DT, v k, ξ DT v, ξ DT. Te proof is similar to tat of 80 were we use Lemma 18 for te proof of 81b and is terefore omitted. Passing to te limit in 38a 38b and using properties prove Items 3 and 5 of Definition 1. Finally, we obtain m0, = m 0, H0, = H 0, and λ0, = λ 0 from te weak convergence and te continuity of te trace operator. Tis and m = 1 yield Statements 1 2 of Definition 1. To obtain 4, note tat Γ : H 1/2 Γ H 1/2 Γ and n n : Hcurl, D H 1/2 Γ are bounded linear operators; see [16, Section 4.2] for exact definition of te spaces and te result. Weak convergence ten proves Item 4 of Definition 1. Estimate 10 follows by weak lower-semicontinuity and te energy bound 33. Tis completes te proof of te teorem Proof of Teorem 6. In tis subsection we invoke Lemmas 20, 21, and 23 to prove a priori error estimates for te pure LLG equation. Proof. It follows from Taylor s Teorem, 15, and Young s inequality tat mt i+1 m i+1 2 H 1 D 1 + k mt i + km t t i m i + kv i 2 H 1 D k 1 k4 4 m tt 2 L 0,T ;H 1 D 1 + k mt i m i 2 H 1 D kk2 m t t i v i 2 H 1 D + 2k1 + k mt i m i, m t t i v i D + 2k1 + k mt i m i, m t t i v i D + k 3 m 2 W 2, 0,T ;H 1 D,

27 THE EDDY CURRENT LLG EQUATIONS 27 recalling tat 0 < k 1. Te tird term on te rigt-and side is estimated as 2k1 + k mti m i, m t t i v i D δ 1 k1 + k mt i m i 2 L 2 D for any δ > 0, so tat 82 becomes mt i+1 m i+1 2 H 1 D + δk1 + k m t t i v i 2 L 2 D, k1 + δ 1 k mt i m i 2 H 1 D + k1 + k k + δ m t t i v i 2 L 2 D + k m tt i v i 2 L 2 D + 2 mt i m i, m t t i v i D + k 3 m 2 W 2, 0,T ;H 1 D. Due to te assumption k < α/2c e we can coose δ = α/2c e suc tat k+δ α/c e and use 67 to deduce mt i+1 m i+1 2 H 1 D 1 + k1 + δ 1 k mt i m i 2 H 1 D + 2kC m 2 + k 2 + mt i m i 2 H 1 D + k 3 m 2 W 2, 0,T ;H 1 D = δ 1 + 2C m k + δ 1 k 2 mt i m i 2 H 1 D + 2kC m 2 + k 2 + k 3 m 2 W 2, 0,T ;H 1 D. Applying Lemma 29 in te Appendix below wit a i := mt i m i 2 H 1 D, b i := 1 + δ 1 + 2C m k + δ 1 k 2, and c i := k 2C m 2 + k 2 + k 2 m 2 W 2, 0,T ;H 1 D, we deduce mt j m j 2 H 1 D e m0 tj m 0 2 H 1 D + t j 2Cm 2 + k 2 + k 2 m 2 W 2, 0,T ;H 1 D m 0 m 0 2 H 1 D k 2, proving 24. To prove 25 we first note tat 84 m m k 2 L 2 0,T ;H 1 D N 1 ti+1 = t i+1 t mt m i t i k + t t i mt m i+1 2 k H D dt 1 N 1 ti+1 mt m i 2 H 1 D t i + mt mi+1 2 H 1 D max 0 i N mt i m i 2 H 1 D + k2 m t 2 L 0,T ;H 1 D, were in te last step we used Taylor s Teorem. Te uniqueness of te strong solution m follows from 25 and te fact tat te m k are uniquely determined by Algoritm 2.5. Wit te weak convergence proved in Teorem 5, we obtain tat tis weak solution coincides wit m. Tis concludes te proof. dt

28 28 M. FEISCHL AND T. TRAN Remark 27. Since m m k 2 H 1 D T = T mt mk t 2 H 1 D + m tt t m k t 2 L 2 D dt 0 N 1 ti+1 t i + m t t v i 2 L 2 D mt m i 2 H 1 D dt + mt mi+1 2 H 1 D max 0 i N mt i m j 2 H 1 D + max 0 i N m tt i v i 2 L 2 D + k 2 m t 2 L 0,T ;H 1 D + k2 m tt 2 L 0,T ;L 2 D, by using 67 for te second term on te rigt-and side and using Young s inequality 2ab ka 2 + k 1 b 2 for te inner product in 67, we obtain a weaker convergence in te H 1 D T -norm, namely m m k H1 D T C conv k 1/2 m 0 m 0 H1 D + + k, provided tat k 1/2 0 wen, k Proof of Teorem 7. Tis section bootstraps te results of te previous section to include te full ELLG system into te analysis. Proof. Similarly to te proof of Teorem 6, we derive 83 wit δ = α/4c e. Multiplying 75 by β = α/4c e C H and adding te resulting equation to 83 yields mt i+1 m i+1 2 H 1 D + β E i β k 2 e i+1 2 L 2 D 1 + k1 + δ 1 k mt i m i 2 H 1 D + k1 + k k + δ + βc H m t t i v i 2 L 2 D + k m tt i v i 2 L 2 D + 2 mt i m i, m t t i v i D + k 3 m 2 W 2, 0,T ;H 1 D kβ E i 2 + βc H k 2 + k 2. Te assumption k α/2c e, see Teorem 7, implies k + δ + βc H α/c e. By invoking Lemma 26 we infer mt i+1 m i+1 2 H 1 D + β E i β k 2 e i+1 2 L 2 D 1 + k1 + δ 1 k mt i m i 2 H 1 D + k1 + kc H 2 + k 2 + mt i m i 2 H 1 D + Ht i H i 2 L 2 D + k 3 m 2 W 2, 0,T ;H 1 D kβ E i 2 + βc H k 2 + k 2 = δ 1 + C H k + δ 1 + C H k 2 mt i m i 2 H 1 D 85 + k1 + kc H Ht i H i 2 L 2 D kβ E i 2 + k 3 m 2 W 2, 0,T ;H 1 D + kc H1 + k + β 2 + k 2. Te approximation properties of Π N D and te regularity of H imply Ht i H i 2 L 2 D E i 2 L 2 D + 2,

29 THE EDDY CURRENT LLG EQUATIONS 29 were te idden constant depends only on te sape regularity of T and on te regularity of H. Hence, we obtain from 85 mt i+1 m i+1 2 H 1 D + β E i β k 2 e i+1 2 L 2 D 1 + C comb k mt i m i 2 H 1 D Ckβ E i 2 + k 3 m 2 W 2, 0,T ;H 1 D + kc2 + k 2, were C comb := 1 + 2δ 1 + 2C H and for some constant C > 0 wic is independent of k, and i. Hence, we find a constant C comb > 0 suc tat 86 mt i+1 m i+1 2 H 1 D + β E i β k 2 e i+1 2 L 2 D 1 + C comb k mt i m i 2 H 1 D + β E i 2 + k Ccomb 2 + k 2. Applying Lemma 29 in te Appendix below wit a i := mt i+1 m i+1 2 H 1 D + β E i β k 2 e i+1 2 L 2 D, b i = C comb k, and c i = k C comb 2 + k 2 we deduce, for all i = 0,..., N, 87 Since 88 mt i+1 m i+1 2 H 1 D + E i k e i+1 2 L 2 D m 0 m 0 2 H 1 D + E k e 0 2 L 2 D + C comb 2 + k 2. Hti H i 2 L 2 D e i 2 L 2 D 2, Ht i H i, λt i λ i 2 X E i 2 X 2, wic is a result of te approximation properties of Π N D and Π S and te regularity assumptions on H and λ estimate 26 follows immediately. To prove 27 it suffices to estimate te term wit k factor on te left-and side of tat inequality because te oter terms can be estimated in exactly te same manner as in te proof of Teorem 6. By using Taylor s Teorem and 88 we deduce 89 H H k 2 L 2 D T N i=1 k Ht i+1 H i+1 2 L 2 D + C Hk 3 N k e i+1 2 L 2 D k 2, i=1 were C H := H 2 W 1, 0,T ;Hcurl,D. On te oter and, it follows from 86 tat k e i+1 2 L 2 D mt i m i 2 H 1 D mt i+1 m i+1 2 H 1 D + k mt i m i 2 H 1 D + β E i 2 E i k E i 2 + k 2 + k 2,

30 30 M. FEISCHL AND T. TRAN wic ten implies by using telescoping series, 87, and 88 N k e i+1 2 L 2 D m0 0 m 0 2 H 1 D + E 0 2 i=1 + max mt i m i 2 H 0 i N 1 D + E i k 2 m 0 m 0 2 H 1 D + H0 H 0 2 L 2 D + λ0 λ 0 2 H 1/2 Γ + max mt i m i 2 H 0 i N 1 D + Ht i H i 2 L 2 D + λt i λ i 2 H 1/2 Γ k 2. Te required result now follows from 89 and 26. Uniqueness is also obtained as in te proof of Teorem 6, completing te proof te teorem. 4. Numerical experiments. Te following numerical experiments are carried out by use of te FEM toolbox FEniCS [28] fenicsproject.org and te BEM toolbox BEM++ [32] bempp.org. We use GMRES to solve te linear systems and blockwise diagonal scaling as preconditioners. Te values of te constants in tese examples are taken from te standard problem #1 proposed by te Micromagnetic Modelling Activity Group at te National Institute of Standards and Tecnology [20]. As domain serves te unit cube D = [0, 1] 3 wit initial conditions { m 0 0, 0, 1 for dx 1/4, x 1, x 2, x 3 := 2Ax 1, 2Ax 2, A 2 dx/a 2 + dx for dx < 1/4, were dx := x x and A := 1 2 dx 4 /4 and H 0 = { 0, 0, 3 in D, 0, 0, 3 m 0 in D. We coose te constants { 1 in D, α = 0.5, σ = 0 in D, µ 0 = , C e = µ Example 1. For time and space discretisation of D T := [0, 5] D, we apply a uniform partition in space = 0.1 and time k = Figure 1 plots te corresponding energies over time. Figure 2 sows a series of magnetisations mt i at certain times t i [0, 5]. Figure 3 sows tat same for te magnetic field Ht i Example 2. We use uniform time and space discretisation of te domain D T := [0, 0.1] D to partition [0, 0.1] into 1/k time intervals for k {0.001, 0.002, 0.004, 0.008, 0.016} and D into ON 3 tetraedra for N {5, 10, 15, 20, 25}. Figure 4 sows convergence rates wit respect to te space discretisation and Figure 5 wit respect to te time discretization. Since te exact solution is unknown, we use te finest computed approximation as a reference solution. Te convergence plots reveal tat te space discretization error dominates te time discretization error by far. Te expected convergence rate Ok can be observed in Figure 5 wic underlines te

31 THE EDDY CURRENT LLG EQUATIONS energy menergy energy sum time Fig. 1. Plot of m k t L 2 D and H kt Hcurl,D over te time. teoretical results of Teorem 7. It is less clear in Figure 4 if tere is a convergence of order O. Preconditioners, a topic of furter study, are required for implementation wit larger values of N. Appendix. Below, we state some well-known results. Lemma 28. Given T, tere exists a constant C norm > 0 wic depends solely on te sape regularity of T suc tat Cnorm w 1 L2 D T T z N T wz wz T 2 1/2 Cnorm w L2 D for all w S 1 T and some arbitrary coice of nodes z T N T for all T T. Proof. Te proof follows from scaling arguments. Lemma 29. If {a i }, {b i }, {c i } are sequences of non-negative numbers satisfying ten for all j N 0 tere olds a i b i a i + c i for all i N 0 a j exp j 1 j 1 b i a0 + c i. Proof. Te lemma can be easily sown by induction. REFERENCES [1] C. Abert, G. Hrkac, M. Page, D. Praetorius, M. Ruggeri, and D. Süss, Spin-polarized transport in ferromagnetic multilayers: An unconditionally convergent FEM integrator, Comput. Mat. Appl., , pp [2] R. A. Adams, Sobolev Spaces, Academic Press [A subsidiary of Harcourt Brace Jovanovic, Publisers], New York-London, Pure and Applied Matematics, Vol. 65. [3] F. Alouges, A new finite element sceme for Landau-Lifcitz equations, Discrete Contin. Dyn. Syst. Ser. S, , pp , doi: /dcdss

32 32 M. FEISCHL AND T. TRAN Fig. 2. Slice of te magnetisation m k t i at [0, 1] 2 {1/2} for i = 0,..., 10 wit t i = 0.2i. Te colour of te vectors represents te magnitude m k. We observe tat te magnetisation aligns itself wit te initial magnetic field H 0 by performing a damped precession. [4] F. Alouges, E. Kritsikis, J. Steiner, and J.-C. Toussaint, A convergent and precise finite element sceme for Landau-Lifscitz-Gilbert equation, Numer. Mat., , pp , doi: /s [5] F. Alouges and A. Soyeur, On global weak solutions for Landau-Lifsitz equations: existence and nonuniqueness, Nonlinear Anal., , pp , doi: / x L. [6] M. Aurada, M. Feiscl, and D. Praetorius, Convergence of some adaptive FEM-BEM coupling for elliptic but possibly nonlinear interface problems, ESAIM Mat. Model. Numer. Anal., , pp , doi: /m2an/ [7] L. Baňas, S. Bartels, and A. Prol, A convergent implicit finite element discretization of te Maxwell Landau Lifsitz Gilbert equation, SIAM J. Numer. Anal., , pp [8] L. Baňas, M. Page, and D. Praetorius, A convergent linear finite element sceme for te Maxwell-Landau-Lifsitz-Gilbert equations, Electron. Trans. Numer. Anal., , pp [9] S. Bartels, Projection-free approximation of geometrically constrained partial differential equations, Mat. Comp., 2015, doi: /mcom/3008. [10] S. Bartels, J. Ko, and A. Prol, Numerical analysis of an explicit approximation sceme for te Landau-Lifsitz-Gilbert equation, Mat. Comp., , pp , doi: /S [11] S. Bartels and A. Prol, Convergence of an implicit finite element metod for te Landau- Lifsitz-Gilbert equation, SIAM J. Numer. Anal., , pp electronic, doi: / [12] J. Berg and J. Löfström, Interpolation spaces. An introduction, Springer-Verlag, Berlin- New York, Grundleren der Matematiscen Wissenscaften, No [13] A. Bossavit, Two dual formulations of te 3-D eddy-currents problem, COMPEL, , pp , doi: /eb [14] F. Bruckner, D. Süss, M. Feiscl, T. Fürer, P. Goldenits, M. Page, D. Praeto-

33 33 THE EDDY CURRENT LLG EQUATIONS Fig. 3. Slice of te magnetic field H k ti at [0, 1]2 {1/2} for i = 0,..., 10 wit ti = 0.2i. Te colour of te vectors represents te magnitude H k. We observe only a sligt movement in te middle of te cube combined wit an overall reduction of field strengt. mnconv 10 1 HNconv error error N N Fig. 4. Te convergence rate of km? mk kl 0,T ;H 1 D left and kh? H k kl2 0,T ;Hcurl,D rigt, were m? = mk and H? = H k for = 1/25 and k = Te dased line indicates O. rius, and M. Ruggeri, Multiscale modeling in micromagnetics: existence of solutions and numerical integration, Mat. Models Metods Appl. Sci., , pp , doi: /s [15] A. Buffa and P. Ciarlet, Jr., On traces for functional spaces related to Maxwell s equations. I. An integration by parts formula in Lipscitz polyedra, Mat. Metods Appl. Sci., , pp. 9 30, doi: / :19::aid-mma191i3.0.co;2-2. [16] A. Buffa and P. Ciarlet, Jr., On traces for functional spaces related to Maxwell s equations. II. Hodge decompositions on te boundary of Lipscitz polyedra and applications, Mat. Metods Appl. Sci., , pp , doi: / :19::