An Effective Integrator for the Landau-Lifshitz-Gilbert Equation
|
|
- June Jones
- 5 years ago
- Views:
Transcription
1 ASC Report No. 02/2012 An Effective Integrator for te Landau-Lifsitz-Gilbert Equation P. Goldenits, G. Hrkac, D. Praetorius, D. Suess Institute for Analysis and Scientific Computing Vienna University of Tecnology TU Wien ISBN
2 Most recent ASC Reports 01/2012 A. Arnold, L. Neumann, W. Hocauser Stability of glued and embedded Glass Panes: Dunkerley straigt Line as a conservative Estimate of superimposed buckling Coefficients 43/2011 G. Hastermann, P. Lima, L. Morgado, E. Weinmüller Numerical Solution of te Density Profile Equation wit p-laplacians 42/2011 M. Aurada, J.M. Melenk, D. Praetorius Mixed Conforming elements for te large-body limit in micromagnetics: a finite element approac 41/2011 P. Amodio, T. Levitina, G. Settanni, E.B. Weinmüller On te Calculation of te Finite Hankel Transform Eigenfunctions 40/2011 D.P. Hewett, S. Langdon, J.M. Melenk A ig frequency p boundary element metod for scattering by convex polygons 39/2011 A. Jüngel Semiconductor Device Problems 38/2011 L. Neumann, A. Arnold, W. Hocauser Zur Stabilität von geklebten und geklotzten Glassceiben: Beurteilung der Dunkerley scen Geraden zur Beulwertbestimmung 37/2011 A. Dick, O. Koc, R. März, E. Weinmüller Collocation Scemes for Nonlinear Index 1 DAEs wit a Singular Point 36/2011 J.A. Carrillo, S. Hittmeir, A. Jüngel Cross diffusion and nonlinear diffusion preventing blow up in te Keller-Segel model 35/2011 P. Fucs, A. Jüngel, M. von Renesse On te Lagrangian structure of quantum fluid models Institute for Analysis and Scientific Computing Vienna University of Tecnology Wiedner Hauptstraße Wien, Austria admin@asc.tuwien.ac.at WWW: ttp:// FAX: ISBN c Alle Recte vorbealten. Nacdruck nur mit Genemigung des Autors. ASC TU WIEN
3 An Effective Integrator for te Landau-Lifsitz-Gilbert Equation P. Goldenits G. Hrkac D. Praetorius D. Suess Institute for Analysis and Scientific Computing, Vienna University of Tecnology, 1040 Wien, Wiedner Hauptstraße 8-10, Austria, ( Institute for Analysis and Scientific Computing, Vienna University of Tecnology, 1040 Wien, Wiedner Hauptstraße 8-10, Austria and Department of Materials Science and Engineering, University of Seffield, S1 3JD Seffield, Mappin Street, UK, ( Institute for Analysis and Scientific Computing, Vienna University of Tecnology, 1040 Wien, Wiedner Hauptstraße 8-10, Austria, ( Institute of Solid State Pysics, Vienna University of Tecnology, 1040 Wien, Wiedner Hauptstraße 8-10, Austria, ( Abstract: We consider a lowest-order finite element sceme for te Landau-Lifsitz-Gilbert equation (LLG) wic describes te dynamics of micromagnetism. In contrast to previous works from te matematics literature, we examine LLG including te total magnetic field induced by pysical penomena described in terms of excange energy, anisotropy energy, magnetostatic energy, as well as Zeeman energy. Besides a strong non-linearity and a non-convex side constraint, te non-local dependence of te demagnetization field from te magnetization represents a callenging task for te numerical integrator. In our numerical sceme, only te igest order term, namely te excange contribution, is treated implicitly, wereas te remaining contributions are computed explicitly. Tis is, in particular, advantageous for te computation of te demagnetization field by means of te popular approac of Fredkin et al. (1990). Furtermore, our sceme requires to solve only one linear system per time-step and allows a simplified computation of te arising system matrices by mass-lumping. Finally, te proposed integrator is matematically reliable in te sense tat we prove unconditional convergence for te approximation of a weak solution. Keywords: finite elements, micromagnetism, Landau-Lifsitz-Gilbert equation, time-integration, non-convex, non-local, non-linear, demagnetization field 1. INTRODUCTION Te understanding of te dynamic beaviour of a micromagnetic body under te influence of certain micromagnetic penomena is essential and of utter relevance for te development of magnetic materials. In te literature, te Landau-Lifsitz-Gilbert equation (LLG) is well-accepted to model te dynamics of micromagnetism. A variety of applications suc as for example te development of magneto-resistive storage devices as well as te amount of numerical issues makes LLG of interest for bot, pysicists and matematicians. In our contribution we consider a polyedral bounded Lipscitz domain R 3 and a fixed time interval (0,τ). Let m 0 : S 2 = {x R 3 x = 1} be some given initial state. Ten, te non-dimensional formulation of LLG reads m t = α 1+α 2m (m eff(m,f)) 1 1+α 2m eff(m,f), (1) m(0) = m 0 in, (2) ν m = 0 on (0,τ), (3) were te (unknown) magnetization is denoted by a vector-valued function m : (0,τ) S 2 wic satisfies te non-convex side constraint m = 1 a.e. in (0,τ). Here,α > 0referstoteGilbertdampingparameterwic depends only on te material. Moreover, m t is te time derivative of m, and eff (m,f) denotes te total magnetic field and is given by te negative variation of te Gibbs Free energy eff (m,f) = δe(m) δm. (4) In tis work, te bulk energy e( ) consists of excange energy, anisotropy energy, magnetostatic energy, as well
4 as Zeeman energy and tus reads e(m) = C ex m 2 +C ani Φ(m) (5) m u f m. 2 Here, Φ refers to te anisotropy density, f denotes an applied external field, and u is te demagnetization field. Te latter is obtained from te magnetostatic Maxwell s equations wit u being te solution u = (u int,u ext ) of te full space transmission problem u int = divm in, u ext = 0 in R 3 \, [u] = 0 on Γ, [ ν u] = m ν on Γ, u ext (x) = O ( 1/ x ) as x. Here, [u] and [ ν u] denote te jumps of u and its normal derivative across te boundary Γ of. Te contribution of te magnetostatic potential, tus involves certain integral operators. Terefore, te computation of te demagnetization field is te most time and memory consuming part in numerical simulations and as to be realized effectively. In Goldenits et al. (2012), we discuss several approaces from te literature to solve (6) numerically. In te present work, we restrict ourselves to te ybrid FEM-BEM approac proposed in Fredkin et al. (1990), wic is mostly used in te pysics literature: Let u 1 be te (up to an additive constant) unique solution of te Neumann problem u 1 = divm in, ν u 1 = m ν on Γ, and extend u 1 by zero to te entire space R 3. Wit u te magnetostatic potential determined by (6), te remainder u 2 = u u 1 satisfies u 2 = 0 in R 3 \Γ, [u 2 ] = u 1 Γ on Γ, [ ν u 2 ] = 0 on Γ, u 2 (x) = O(1/ x ) as x. Here and in te following, u 1 Γ denotes te interior trace u 1 Γ of u 1. As is known from potential teory, te unique solution of (8) is te double-layer potential u 2 = Ku 1 Γ were (Ku 1 Γ )(x) := 1 (y x) ν(y) 4π Γ y x 3 u 1 (y) (9) for all x R 3 \Γ and wit ν te exterior unit normal vector on Γ. According to te jump of K across Γ, one can sow tat u 2 on is caracterized by te inomogeneous Diriclet problem u 2 = 0 in, u 2 = (K 1/2)u 1 Γ on Γ, and we ave u = u 1 + u 2 in. (6) (7) (8) (10) Remark. We stress tat te factor c = 1/2 for te trace jump (K c)u 1 Γ of te double-layer potential Ku 1 Γ in olds only almost everywere on Γ, were Γ is flat. At corners or on edges of Γ, te factor c depends on te interior angle of. In order to provide a numerical sceme, we obtain tat supplemented by te same initial and boundary conditions (2) and (3), te classical formulation of LLG, cf. (1), can equivalently be stated as αm t +m m t = eff (m,f) (m eff (m,f))m, (11) cf. e.g. Goldenits (2012) for a detailed proof. We empasis, tat tis alternative formulation still is non-linear in consideration of te magnetization m but is linear in consideration of its time derivative m t. Formulation (11) will serve as te basis for our finite element (FE) sceme to solve LLG numerically,wereweapproximatem (t, ) m(t, ) and v (t, ) m t (t, ) for all times t (0,τ). Note tat, due to te non-convex constraint m = 1 a.e., te time derivative m t belongs to te tangential space of m, i.e. m m t = 0 a.e. in τ. 2. NUMERICAL ALGORITHM Let T denote a quasi-uniform and regular triangulation of te domain into tetraedra and N be te set of its nodes. To discretize te magnetization m in te spatial variable, we use te vector-valued Courant FE space V = S 1 (T ) 3 of piecewise linear and globally continuous functions. To discretize te time interval, we consider a uniform partition 0 = τ 0 < τ 1 < < τ J = τ wit timestep size k = k j := τ j+1 τ j for j = 0,...,J 1. For eac (discrete) function ϕ, ϕ j denotes te evaluation ϕ(τ j ) at time τ j. Now, let M = { } φ V φ (z) = 1 for all z N (12) be te restricted finite element set, were our solution m j m(τ j) is sougt due to te non-convex constraint m = 1. Furtermore, for φ M, let K φ = { } ψ V ψ (z) φ (z) = 0 for all z N (13) be te discrete tangential space associated wit φ M, were te discrete time derivative v j m t(τ j ) is sougt. To obtain a numerical integrator for LLG, we follow te idea of Alouges (2008), were te small particle limit eff (m) = m is considered: We proceed by setting v = m t in (11) and by discretizing te weak form of it according to our framework.we treat te term of igest order implicitly, namely te excange contribution, wereas te remaining tree terms of te effective magnetic field eff (m,f) = C ex m+ low (m,f) wit low (m,f) := C ani DΦ(m) u+f (14) are computed explicitly. For te numerical integrator, we set low (m j,fj ) := C anidφ(m j ) uj +fj. (15) Here, u j := uj 1 + uj 2 S1 (T ) denotes an FE solution of te superposition ansatz of Fredkin et al. (1990), were we proceed as follows: First, let u 1 S 1 (T ) wit e.g. u 1 = 0 be te unique FE solution of u 1 v = m v (16) for all v S 1 (T ) wit v = 0. Second, let u 2 S 1 (T ) be te unique solution of te Diriclet problem u 2 v = 0 (17)
5 for all v S 1 (T ) wit v Γ = 0, wic additionally satisfies te inomogeneous discrete Diriclet condition u 2 Γ = S (K 1/2)u 1 Γ. Here, S maps te continuous boundary data (K 1/2)u 1 Γ onto some discrete boundary data in S 1 (T Γ ), were T Γ denotes te induced triangulation of te boundary Γ into flat surface triangles. A possible coice for te approximation operator S is given by te Clément-Operator J v := z N Γ v z β z wit v z := 1 γ z Here, γ z denotes te surface node-patc γ z := {T Γ : T T wit z T} Γ γ z v. of a boundary node z N Γ. Moreover, β z S 1 (T ) is te corresponding at function. Finally, te discrete demagnetization field is defined by u = u 1 + u 2. Considering te last term in (15), we stress tat f j f(τ j) approximates te given applied field f. If f is continuous in time, a valid coice is te evaluation f j := f(τ j) of f at time τ j. If f is continuous in space and time, f j V can be cosen to be te nodal interpolant of f(τ j ) wic furter simplifies te implementation. Finally, wit te nodal interpolation operator I : C() S 1 (T ), we include te so-called mass-lumping of te L 2 -scalar-product to compute te arising massmatrices only approximately. Tis results in corresponding matrix blocks wic are diagonal instead of sparse only. Te proposed time-splitting sceme now reads as follows: Algoritm. Input: Discretized initial data m 0 M, damping parameter α > 0, parameter 0 < θ 1, counter j = 0. (i) Compute v j K m j by solving te (regular) linear system α I (v j ψ )+ I ((m j vj ) ψ ) = C ex (m j +θkvj ) ψ (18) + low (m j,fj ) ψ for all test functions ψ K m j tangential space. (ii) Define m j+1 M by setting from te discrete m j+1 (z) = mj (z)+kvj (z) m j (z)+kvj (z) (19) for all nodes z N and go to (i). Output: Sequence of functions v j K m j as well as m j+1 M for j 0. Some remarks are in order to comment on te wellposedness and te advantages of te proposed algoritm: According to te Lemma of Lax-Milgram,(18) admits a unique solution v j K m j. Asaconsequenceofteortogonalityrelationm j (z) v j (z) = 0, one as mj (z) + kvj (z) 2 = m j (z) 2 + k 2 v j (z) 2 1. Terefore, te discretized magnetization m j+1 M defined in step (ii) of our Algoritm is well-defined. We stress tat only one (sparse) linear system (18) as to be solved per time-step and te non-convex side constraint m = 1 is fulfilled node-wise. Te assembly of tis system is te topic of te subsequent Section 3. Altoug, one may also drop te nodal interpolation I in te linear system (18), we put empasis on te fact tat te use of mass-lumping for te cross product contribution of (18) is implementationally very attractive. An explicit treatment of te non-local contribution stemming from te magnetostatic potential u is included, i.e. te computation of m j+1 only requires te approximate field u j from te previous timestep. Put differently, te approac by Fredkin et al. (1990) is only used once per time-step. Tis results in te solution of two additional (sparse) linear systems per time-step. Formally, te Crank-Nicolson type sceme θ = 1/2 is of second order in time, wereas θ = 1 corresponds to an implicit Euler sceme. Instead of te approac of Fredkin et al. (1990) to approximate te demagnetization field u, also oter approaces can be used. Te same applies for te assumptions on and te discretization of te exterior field f. We refer to te remarksin Section 4 and Goldenits et al. (2012) for furter details. 3. IMPLEMENTATION In tis section, we focus on te computation of v j K m j in step (i) of our algoritm. Empasis is put on te assembly and structure of te matrices of te linear system(18),wicisposedontesubspacek m j oftefe space V = S 1 (T ) 3. Let β i S 1 (T ) denote te canonical at function associated wit te node z i N. We define β i+(l 1)n := β i e l wit te l-t unit vector e l R 3 and note tat β 1,...,β 3N is te canonical basis of V. In a first step, we rewrite te variational form (18) as follows: Find v j K m j suc tat a(v j,ψ )+b j (v j,ψ ) = L j (ψ ) (20) for all ψ K m j, were we use te abbreviate notation
6 a(φ,ψ ) = α I (φ ψ ) +θkc ex φ ψ, b j (φ,ψ ) = (21) I ((m j φ ) ψ ), (22) L j (ψ ) = C ex m j ψ C ani DΦ(m j ) ψ u j ψ + f j ψ. (23) Due to te coice of te basis functions β i+(l 1)n = β i e l, te bilinear form a(, ) from (21) corresponds to a symmetric block diagonal matrix A = A0 A 0 R 3N 3N A 0 sym (24) wit symmetric blocks A 0 R N N sym, were A 0 ii = α I (β i β i )+θkc ex β i β i = α ω i δ ii +θkc ex β i β i, (25) wit ω i = {T T : z i T} te volume patc associated wit te node z i N. Here, δ ii denotes Kronecker s delta wit δ ii = 1 for i = i and δ ii = 0 oterwise. Te occurring integrals can be computed by closed formulae. Note tat A 0 is te (positively weigted) sum of te standard stiffness and an approximated diagonal mass matrix. Consequently, A 0 and ence A are positive definite sparse matrices and do not depend on te time step τ j. For te bilinear form b j (, ) from (22), we use te identity a (b c) = (a b) c to see b j (φ,ψ ) = I ((m j φ ) ψ ) = I ((φ ψ ) m j ). To derive te corresponding matrix B R 3N 3N, note tat β i e l β i e l = β i β i e l e l and e 1 e 2 = e 3, e 2 e 3 = e 1, e 3 e 1 = e 2, and e l e l = 0. By coice of te basis functions β i+(l 1)n = β i e l, B terefore as some block structure of te type B = 0 +B3 B 2 B 3 0 +B 1 R 3N 3N (26) +B 2 B 1 0 wit diagonal blocks B l R N N, were B l ii = I (β i β i m j e l) = ω i m j (z i) e l δ ii. Clearly, B is a sparse matrix wic is skew-symmetric and positive semidefinite, since b j (φ,ψ ) = b j (ψ,φ ) and b j (φ,φ ) = 0. Altogeter, te system matrix M = A + B R 3N 3N is positive definite and ence regular. Terefore te variational formulation (20) as even a unique solution if posed on te entire space V = S 1 (T ) 3. However, v j is determined by solving (20) in te subspace K m j. We realize te linear constraints m j (z) vj (z) = 0 for all z N (27) by a Lagrange multiplier ansatz, i.e. we obtain te unknown coefficients ν R 3N of 3N v j = ( )( ) ( M Λ ν m β m by solution of T ν b =. Λ 0 λ 0) m=1 Here, λ R N is te Lagrange multiplier and Λν = 0 realizes te constraints (27). Consequently, te Lagrange matrix reads Λ = ( Λ 1 Λ 2 Λ 3) R N 3N wit Λ l R N N, Λ l ii = k2 m j (z i) e l δ ii, i.e. te matricesλ l arediagonaland scaledby k 2 in order to stabilize te sceme. 4. CONVERGENCE RESULT Te definition of a weak solution to LLG is based on te idea of Alouges et al. (1992) and reads as follows: Definition. Let m 0 H 1 (;S 2 ) be a given initial magnetization. Ten, m is called a weak solution to LLG, if tere olds for all times τ > 0: (i) m H 1 ( τ ;S 2 ) wit m(x,0) = m 0 (x) in te sense of traces; (ii) for all φ C0 ( τ ;R 3 ), tere olds m t φ α (m m t ) φ τ τ (28) = ( eff (m,f) m) φ; τ (iii) for almost all t (0,τ), tere olds 1 m(t) 2 +C 1 m t 2 2 τ 1 (29) m 0 2 +C 2, 2 wit positive constants C 1,C 2 > 0 wic depends only on f, and α > 0. We interpolate te discrete solution m j for j = 1,...,J of te numerical integrator from Section 2 as a continuous, piecewise affine function in time: For all x and all times t (0,τ) wit j = {0,...,J 1} suc tat t [jk,(j +1)k), we define m k (t,x) := t jk m j+1 (j +1)k t (x)+ m j (x). (30) k k Te following convergence teorem generalizes te result of Alouges (2008), yields reliability of te proposed algoritm, and even proves existence of global weak solutions. We stress tat no coupling of te time-step size k and te space-mes size is imposed.
7 Convergence Teorem. Let θ (1/2,1] be a fixed parameter and let T be a family of sape-regular triangulations of te magnetic domain wit mes-size 0. Let te at functions {β i } associated wit T fulfill β i β i 0 for all i i. (31) For te discrete initial magnetization m 0 M, we assume m 0 m 0 in H 1 () as 0. Ten, tere olds weak convergence in H 1 ( τ ) of a subsequence of te in time interpolated discrete magnetization m k from (30) to a weak solution m of LLG as (,k) 0. Remark. Due to Bartels (2005), condition (31) implies te following energy estimate for te excange energy in step (ii) of te algoritm: I ( φ φ ) 2 φ 2 (32) for all φ V wit φ (z) 1 and for all nodes z N. Remark. Under certain assumptions one may approximate te demagnetization field u as well as te applied external field f differently witout canging te result of te convergence teorem: We assume tat te approximation P m j related to te demagnetization field Pm = u as its discrete counterpart, fulfills te following properties P m j L 2 () C P m j L 2 () (33) as well as Pm P m L2 ( τ) 0 as 0 (34) wit a positive constant C P > 0. We note tat te approximate stray-field operator P m j = u,wicisgivenbyteapproacoffredkin et al. (1990) and discussed above, fulfills property (33) as well as (34). A detailed proof is given in Goldenits et al. (2012). We consider te approximation f k in space and time of an applied external field f, wic is given by f k (t,x) = f j (x) for all t j t t j+1,x. Wenotetatanydiscretizationf k off wicsatisfies f f k L2 ( τ) 0 as (,k) 0 (35) will provide an admissible coice in te sense tat te convergence teorem olds. In particular, it is also possible to deal wit discontinuous applied fields. Te discretization of continuous f by nodal interpolation in time or in space-time, as proposed in Section 2, guarantees (35). 5. SKETCH OF PROOF In tis section we briefly comment on te essential arguments to prove te convergence teorem. A rigorously elaborated proof is given in Goldenits (2012). In a first step, we aim for convergence properties of te output v j and mj for j > 0 of our algoritm. To tis end, we also interpret te output function v j as a piecewise constant function in time: For all x and all times t (0,τ) wit j = {0,...,J} suc tat t [jk,(j +1)k), we define v k (t,x) := v j (x). On te one and, te definition of m k in (30) as well as te non-convex property m j (x) 1 for j 0 and all x, wic olds due to step (ii) of our time-splitting sceme, implies m k (t,x) 1. On te oter and and by use of te test function ψ = v j in (18), one may deduce some stability estimate J 1 1 m J 2 2 +kc v j 2 j=0 1 m (36) + k J 1 low (m j C,fj ) vj. ex j=0 Tis yields uniform boundedness of m k in H 1 ( τ ) as well as of v k in L 2 ( τ ). Terefore, we may extract a subsequence of v k and m k, respectively, suc tat m k m weakly in H 1 ( τ ) as (,k) 0, (37) m k m strongly in L 2 ( τ ) as (,k) 0, (38) v k v weakly in L 2 ( τ ) as (,k) 0, (39) wit certain v L 2 ( τ ;R 3 ) and m H 1 ( τ ;R 3 ). Te remaining part of te proof is dedicated to te stepwise verification tat te limit m is a weak solution to LLG, see Definition 4. We remark tat tis also involves te proof of v = m t to link te limits (37) (39). ACKNOWLEDGEMENTS Te autors acknowledge support of te Viennese Science and Tecnology Fund (WWTF) troug te researc project MA Micromagnetic Simulation and Computational Design of Future Devices. REFERENCES F. Alouges. A new finite element sceme for Landau- Lifcitz equations. Discrete and Continuous Dynamical Systems Series S, Vol. 1, 2008, F. Alouges, and A. Soyeur. On global weak solutions for Landau-Lifsitz equations: existence and nonuniqueness. Nonlinear Analysis, 18: , S. Bartels. Stability and convergence of finite-element approximation scemes for armonic maps. SIAM J. Numer. Anal., 43: , D. Fredkin, and T. Koeler. Hybrid metod for computing demagnetizing fields. IEEE Transactions on Magnetics, 26: , P. Goldenits, G. Hrkac, M. Page, D. Praetorius, T. Screfl, and D. Süss. Effective Simulation of te Dynamics of Ferromagnetism. Work in progress, P. Goldenits. Pd-Tesis. Institute for Analysis and Scientific Computing, Vienna University of Tecnology.
Multiscale Modeling in Micromagnetics: Well-Posedness and Numerical Integration
ASC Report No. 34/212 Multiscale Modeling in Micromagnetics: Well-Posedness and Numerical Integration F. Bruckner, M. Feiscl, T. Fürer, P. Goldenits, M. Page, D. Praetorius and D. Süss Institute for Analysis
More informationCONVERGENCE OF AN IMPLICIT FINITE ELEMENT METHOD FOR THE LANDAU-LIFSHITZ-GILBERT EQUATION
CONVERGENCE OF AN IMPLICIT FINITE ELEMENT METHOD FOR THE LANDAU-LIFSHITZ-GILBERT EQUATION SÖREN BARTELS AND ANDREAS PROHL Abstract. Te Landau-Lifsitz-Gilbert equation describes dynamics of ferromagnetism,
More informationA finite element approximation for the quasi-static Maxwell Landau Lifshitz Gilbert equations
ANZIAM J. 54 (CTAC2012) pp.c681 C698, 2013 C681 A finite element approximation for te quasi-static Maxwell Landau Lifsitz Gilbert equations Kim-Ngan Le 1 T. Tran 2 (Received 31 October 2012; revised 29
More informationAn extended midpoint scheme for the Landau-Lifschitz-Gilbert equation in computational micromagnetics. Bernhard Stiftner
AANMPDE(JS) Strobl, July 7, 2016 An extended midpoint sceme for te Landau-Lifscitz-Gilbert equation in computational micromagnetics Bernard Stiftner joint work wit Dirk Praetorius, Micele Ruggeri TU Wien
More informationarxiv: v2 [math.na] 8 Nov 2017
CONVERGENCE OF AN IMPLICIT-EXPLICIT MIDPOINT SCHEME FOR COMPUTATIONAL MICROMAGNETICS DIRK PRAETORIUS, MICHELE RUGGERI, AND BERNHARD STIFTNER arxiv:1611.2465v2 [mat.na] 8 Nov 217 Abstract. Based on lowest-order
More informationOn the Landau-Lifshitz-Gilbert equations with magnetostriction
ASC Report No. 8/213 On te Landau-Lifsitz-Gilbert equations wit magnetostriction L. Banas, M. Page, D. Praetorius, and J. Rocat Institute for Analysis and Scientific Computing Vienna University of Tecnology
More informationNUMERICAL APPROXIMATION OF THE LANDAU-LIFSHITZ-GILBERT EQUATION AND FINITE TIME BLOW-UP OF WEAK SOLUTIONS
NUMERICAL APPROXIMATION OF THE LANDAU-LIFSHITZ-GILBERT EQUATION AND FINITE TIME BLOW-UP OF WEAK SOLUTIONS SÖREN BARTELS, JOY KO, AND ANDREAS PROHL Abstract. Te Landau-Lifsitz-Gilbert equation describes
More informationarxiv: v1 [math.na] 29 Nov 2017
LINEAR SECOND-ORDER IMEX-TYPE INTEGRATOR FOR THE (EDDY CURRENT) LANDAU LIFSHITZ GILBERT EQUATION GIOVANNI DI FRATTA, CARL-MARTIN PFEILER, DIRK PRAETORIUS, MICHELE RUGGERI, BERNHARD STIFTNER arxiv:1711.1715v1
More informationPoisson Equation in Sobolev Spaces
Poisson Equation in Sobolev Spaces OcMountain Dayligt Time. 6, 011 Today we discuss te Poisson equation in Sobolev spaces. It s existence, uniqueness, and regularity. Weak Solution. u = f in, u = g on
More informationarxiv: v2 [math.na] 5 Feb 2017
THE EDDY CURRENT LLG EQUATIONS: FEM-BEM COUPLING AND A PRIORI ERROR ESTIMATES MICHAEL FEISCHL AND THANH TRAN arxiv:1602.00745v2 [mat.na] 5 Feb 2017 Abstract. We analyze a numerical metod for te coupled
More informationAn Implicit Finite Element Method for the Landau-Lifshitz-Gilbert Equation with Exchange and Magnetostriction
SMA An Implicit Finite Element Metod for te Landau-Lifsitz-Gilbert Equation wit Excange and Magnetostriction Jonatan ROCHAT under te supervision of L. Banas and A. Abdulle 2 ACKNOWLEDGEMENTS I would like
More information1. Introduction. We consider the model problem: seeking an unknown function u satisfying
A DISCONTINUOUS LEAST-SQUARES FINITE ELEMENT METHOD FOR SECOND ORDER ELLIPTIC EQUATIONS XIU YE AND SHANGYOU ZHANG Abstract In tis paper, a discontinuous least-squares (DLS) finite element metod is introduced
More informationNumerical quadratic energy minimization bound to convex constraints in thin-film micromagnetics
ASC Report No. 32/2011 Numerical quadratic energy minimization bound to convex constraints in tin-film micromagnetics S. Ferraz-Leite, JM. Melenk, D. Praetorius Institute for Analysis and Scientific Computing
More informationPROJECTION-FREE APPROXIMATION OF GEOMETRICALLY CONSTRAINED PARTIAL DIFFERENTIAL EQUATIONS
PROJECTION-FREE APPROXIMATION OF GEOMETRICALLY CONSTRAINED PARTIAL DIFFERENTIAL EQUATIONS Abstract. We devise algoritms for te numerical approximation of partial differential equations involving a nonlinear,
More informationParameter Fitted Scheme for Singularly Perturbed Delay Differential Equations
International Journal of Applied Science and Engineering 2013. 11, 4: 361-373 Parameter Fitted Sceme for Singularly Perturbed Delay Differential Equations Awoke Andargiea* and Y. N. Reddyb a b Department
More informationA Hybrid Mixed Discontinuous Galerkin Finite Element Method for Convection-Diffusion Problems
A Hybrid Mixed Discontinuous Galerkin Finite Element Metod for Convection-Diffusion Problems Herbert Egger Joacim Scöberl We propose and analyse a new finite element metod for convection diffusion problems
More informationPreconditioning in H(div) and Applications
1 Preconditioning in H(div) and Applications Douglas N. Arnold 1, Ricard S. Falk 2 and Ragnar Winter 3 4 Abstract. Summarizing te work of [AFW97], we sow ow to construct preconditioners using domain decomposition
More informationMIXED DISCONTINUOUS GALERKIN APPROXIMATION OF THE MAXWELL OPERATOR. SIAM J. Numer. Anal., Vol. 42 (2004), pp
MIXED DISCONTINUOUS GALERIN APPROXIMATION OF THE MAXWELL OPERATOR PAUL HOUSTON, ILARIA PERUGIA, AND DOMINI SCHÖTZAU SIAM J. Numer. Anal., Vol. 4 (004), pp. 434 459 Abstract. We introduce and analyze a
More information1. Introduction. Consider a semilinear parabolic equation in the form
A POSTERIORI ERROR ESTIMATION FOR PARABOLIC PROBLEMS USING ELLIPTIC RECONSTRUCTIONS. I: BACKWARD-EULER AND CRANK-NICOLSON METHODS NATALIA KOPTEVA AND TORSTEN LINSS Abstract. A semilinear second-order parabolic
More informationDifferent Approaches to a Posteriori Error Analysis of the Discontinuous Galerkin Method
WDS'10 Proceedings of Contributed Papers, Part I, 151 156, 2010. ISBN 978-80-7378-139-2 MATFYZPRESS Different Approaces to a Posteriori Error Analysis of te Discontinuous Galerkin Metod I. Šebestová Carles
More informationRate optimal adaptive FEM with inexact solver for strongly monotone operators
ASC Report No. 19/2016 Rate optimal adaptive FEM with inexact solver for strongly monotone operators G. Gantner, A. Haberl, D. Praetorius, B. Stiftner Institute for Analysis and Scientific Computing Vienna
More informationSome Error Estimates for the Finite Volume Element Method for a Parabolic Problem
Computational Metods in Applied Matematics Vol. 13 (213), No. 3, pp. 251 279 c 213 Institute of Matematics, NAS of Belarus Doi: 1.1515/cmam-212-6 Some Error Estimates for te Finite Volume Element Metod
More informationA SHORT INTRODUCTION TO BANACH LATTICES AND
CHAPTER A SHORT INTRODUCTION TO BANACH LATTICES AND POSITIVE OPERATORS In tis capter we give a brief introduction to Banac lattices and positive operators. Most results of tis capter can be found, e.g.,
More informationarxiv: v2 [math.na] 5 Jul 2017
Trace Finite Element Metods for PDEs on Surfaces Maxim A. Olsanskii and Arnold Reusken arxiv:1612.00054v2 [mat.na] 5 Jul 2017 Abstract In tis paper we consider a class of unfitted finite element metods
More informationNumerische Mathematik
Numer. Mat. (2012) 122:101 131 DOI 10.1007/s00211-012-0454-z Numerisce Matematik Numerical quadratic energy minimization bound to convex constraints in tin-film micromagnetics Samuel Ferraz-Leite Jens
More informationStability and convergence of finite element approximation schemes for harmonic maps
Stability and convergence of finite element approximation scemes for armonic maps Sören Bartels Department of Matematics, Humboldt-Universität zu Berlin, Unter den Linden 6, D-199 Berlin, Germany E-Mail:
More informationNew Streamfunction Approach for Magnetohydrodynamics
New Streamfunction Approac for Magnetoydrodynamics Kab Seo Kang Brooaven National Laboratory, Computational Science Center, Building 63, Room, Upton NY 973, USA. sang@bnl.gov Summary. We apply te finite
More informationStability properties of a family of chock capturing methods for hyperbolic conservation laws
Proceedings of te 3rd IASME/WSEAS Int. Conf. on FLUID DYNAMICS & AERODYNAMICS, Corfu, Greece, August 0-, 005 (pp48-5) Stability properties of a family of cock capturing metods for yperbolic conservation
More informationA Mixed-Hybrid-Discontinuous Galerkin Finite Element Method for Convection-Diffusion Problems
A Mixed-Hybrid-Discontinuous Galerkin Finite Element Metod for Convection-Diffusion Problems Herbert Egger Joacim Scöberl We propose and analyse a new finite element metod for convection diffusion problems
More informationQuasi-optimal a priori estimates for fluxes in mixed finite element methods and applications to the Stokes-Darcy coupling
ASC Report No. 05/2012 Quasi-optimal a priori estimates for fluxes in mixed finite element metods and applications to te Stokes-Darcy coupling J.M. Melenk, H. Rezaijafari, B. Wolmut Institute for Analysis
More informationERROR BOUNDS FOR THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BRADLEY J. LUCIER*
EO BOUNDS FO THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BADLEY J. LUCIE* Abstract. Te expected error in L ) attimet for Glimm s sceme wen applied to a scalar conservation law is bounded by + 2 ) ) /2 T
More informationFEM solution of the ψ-ω equations with explicit viscous diffusion 1
FEM solution of te ψ-ω equations wit explicit viscous diffusion J.-L. Guermond and L. Quartapelle 3 Abstract. Tis paper describes a variational formulation for solving te D time-dependent incompressible
More informationDownloaded 06/08/17 to Redistribution subject to SIAM license or copyright; see
SIAM J. NUMER. ANAL. Vol. 55, No. 3, pp. 357 386 c 7 Society for Industrial and Applied Matematics Downloaded 6/8/7 to 67.96.45.. Redistribution subject to SIAM license or copyrigt; see ttp://www.siam.org/journals/ojsa.pp
More informationarxiv: v2 [math.na] 11 Dec 2016
Noname manuscript No. will be inserted by te editor Sallow water equations: Split-form, entropy stable, well-balanced, and positivity preserving numerical metods Hendrik Ranoca arxiv:609.009v [mat.na]
More informationWeierstraß-Institut. im Forschungsverbund Berlin e.v. Preprint ISSN
Weierstraß-Institut für Angewandte Analysis und Stocastik im Forscungsverbund Berlin e.v. Preprint ISSN 0946 8633 Stability of infinite dimensional control problems wit pointwise state constraints Micael
More informationLEAST-SQUARES FINITE ELEMENT APPROXIMATIONS TO SOLUTIONS OF INTERFACE PROBLEMS
SIAM J. NUMER. ANAL. c 998 Society for Industrial Applied Matematics Vol. 35, No., pp. 393 405, February 998 00 LEAST-SQUARES FINITE ELEMENT APPROXIMATIONS TO SOLUTIONS OF INTERFACE PROBLEMS YANZHAO CAO
More informationDifferentiation in higher dimensions
Capter 2 Differentiation in iger dimensions 2.1 Te Total Derivative Recall tat if f : R R is a 1-variable function, and a R, we say tat f is differentiable at x = a if and only if te ratio f(a+) f(a) tends
More informationDiscrete Monogenic Functions in Image processing
Discrete Function Teory Discrete Monogenic Functions in Image processing U. Käler CIDMA and Departamento de Matemática Universidade de Aveiro ukaeler@ua.pt MOIMA Scloss Herrenausen Hannover, June 20-24,
More informationOn the best approximation of function classes from values on a uniform grid in the real line
Proceedings of te 5t WSEAS Int Conf on System Science and Simulation in Engineering, Tenerife, Canary Islands, Spain, December 16-18, 006 459 On te best approximation of function classes from values on
More informationMANY scientific and engineering problems can be
A Domain Decomposition Metod using Elliptical Arc Artificial Boundary for Exterior Problems Yajun Cen, and Qikui Du Abstract In tis paper, a Diriclet-Neumann alternating metod using elliptical arc artificial
More information5 Ordinary Differential Equations: Finite Difference Methods for Boundary Problems
5 Ordinary Differential Equations: Finite Difference Metods for Boundary Problems Read sections 10.1, 10.2, 10.4 Review questions 10.1 10.4, 10.8 10.9, 10.13 5.1 Introduction In te previous capters we
More informationAPPROXIMATION OF CRYSTALLINE DENDRITE GROWTH IN TWO SPACE DIMENSIONS. Introduction
Acta Mat. Univ. Comenianae Vol. LXVII, 1(1998), pp. 57 68 57 APPROXIMATION OF CRYSTALLINE DENDRITE GROWTH IN TWO SPACE DIMENSIONS A. SCHMIDT Abstract. Te pase transition between solid and liquid in an
More informationMATH745 Fall MATH745 Fall
MATH745 Fall 5 MATH745 Fall 5 INTRODUCTION WELCOME TO MATH 745 TOPICS IN NUMERICAL ANALYSIS Instructor: Dr Bartosz Protas Department of Matematics & Statistics Email: bprotas@mcmasterca Office HH 36, Ext
More informationarxiv: v1 [math.na] 20 Jul 2009
STABILITY OF LAGRANGE ELEMENTS FOR THE MIXED LAPLACIAN DOUGLAS N. ARNOLD AND MARIE E. ROGNES arxiv:0907.3438v1 [mat.na] 20 Jul 2009 Abstract. Te stability properties of simple element coices for te mixed
More informationChemnitz Scientific Computing Preprints
G. Of G. J. Rodin O. Steinbac M. Taus Coupling Metods for Interior Penalty Discontinuous Galerkin Finite Element Metods and Boundary Element Metods CSC/11-02 Cemnitz Scientific Computing Preprints Impressum:
More informationAN ANALYSIS OF THE EMBEDDED DISCONTINUOUS GALERKIN METHOD FOR SECOND ORDER ELLIPTIC PROBLEMS
AN ANALYSIS OF THE EMBEDDED DISCONTINUOUS GALERKIN METHOD FOR SECOND ORDER ELLIPTIC PROBLEMS BERNARDO COCKBURN, JOHNNY GUZMÁN, SEE-CHEW SOON, AND HENRYK K. STOLARSKI Abstract. Te embedded discontinuous
More informationImplicit-explicit variational integration of highly oscillatory problems
Implicit-explicit variational integration of igly oscillatory problems Ari Stern Structured Integrators Worksop April 9, 9 Stern, A., and E. Grinspun. Multiscale Model. Simul., to appear. arxiv:88.39 [mat.na].
More informationNumerical Experiments Using MATLAB: Superconvergence of Nonconforming Finite Element Approximation for Second-Order Elliptic Problems
Applied Matematics, 06, 7, 74-8 ttp://wwwscirporg/journal/am ISSN Online: 5-7393 ISSN Print: 5-7385 Numerical Experiments Using MATLAB: Superconvergence of Nonconforming Finite Element Approximation for
More informationNotes on Multigrid Methods
Notes on Multigrid Metods Qingai Zang April, 17 Motivation of multigrids. Te convergence rates of classical iterative metod depend on te grid spacing, or problem size. In contrast, convergence rates of
More informationChapter 5 FINITE DIFFERENCE METHOD (FDM)
MEE7 Computer Modeling Tecniques in Engineering Capter 5 FINITE DIFFERENCE METHOD (FDM) 5. Introduction to FDM Te finite difference tecniques are based upon approximations wic permit replacing differential
More informationAnalysis of A Continuous Finite Element Method for H(curl, div)-elliptic Interface Problem
Analysis of A Continuous inite Element Metod for Hcurl, div)-elliptic Interface Problem Huoyuan Duan, Ping Lin, and Roger C. E. Tan Abstract In tis paper, we develop a continuous finite element metod for
More informationRecent Advances in Time-Domain Maxwell s Equations in Metamaterials
Recent Advances in Time-Domain Maxwell s Equations in Metamaterials Yunqing Huang 1, and Jicun Li, 1 Hunan Key Laboratory for Computation and Simulation in Science and Engineering, Xiangtan University,
More informationAdaptive Finite Element Method
39 Capter 3 Adaptive inite Element Metod 31 Introduction As already pointed out in capter 2, singularities occur in interface problems Wen discretizing te problem (221) wit inite Elements, te singularities
More informationVariational Localizations of the Dual Weighted Residual Estimator
Publised in Journal for Computational and Applied Matematics, pp. 192-208, 2015 Variational Localizations of te Dual Weigted Residual Estimator Tomas Ricter Tomas Wick Te dual weigted residual metod (DWR)
More information3 Parabolic Differential Equations
3 Parabolic Differential Equations 3.1 Classical solutions We consider existence and uniqueness results for initial-boundary value problems for te linear eat equation in Q := Ω (, T ), were Ω is a bounded
More informationMass Lumping for Constant Density Acoustics
Lumping 1 Mass Lumping for Constant Density Acoustics William W. Symes ABSTRACT Mass lumping provides an avenue for efficient time-stepping of time-dependent problems wit conforming finite element spatial
More informationFinite Difference Methods Assignments
Finite Difference Metods Assignments Anders Söberg and Aay Saxena, Micael Tuné, and Maria Westermarck Revised: Jarmo Rantakokko June 6, 1999 Teknisk databeandling Assignment 1: A one-dimensional eat equation
More informationComputing eigenvalues and eigenfunctions of Schrödinger equations using a model reduction approach
Computing eigenvalues and eigenfunctions of Scrödinger equations using a model reduction approac Suangping Li 1, Ziwen Zang 2 1 Program in Applied and Computational Matematics, Princeton University, New
More informationA h u h = f h. 4.1 The CoarseGrid SystemandtheResidual Equation
Capter Grid Transfer Remark. Contents of tis capter. Consider a grid wit grid size and te corresponding linear system of equations A u = f. Te summary given in Section 3. leads to te idea tat tere migt
More informationarxiv: v3 [math.na] 25 Jun 2017
ANALYSIS OF A HIGH ORDER UNFITTED FINITE ELEMENT METHOD FOR ELLIPTIC INTERFACE PROBLEMS CHRISTOPH LEHRENFELD AND ARNOLD REUSKEN arxiv:1602.02970v3 [mat.na] 25 Jun 2017 Abstract. In te context of unfitted
More informationSIMULATION OF Q-TENSOR FIELDS WITH CONSTANT ORIENTATIONAL ORDER PARAMETER IN THE THEORY OF UNIAXIAL NEMATIC LIQUID CRYSTALS
SIMULATION OF Q-TENSOR FIELDS WITH CONSTANT ORIENTATIONAL ORDER PARAMETER IN THE THEORY OF UNIAXIAL NEMATIC LIQUID CRYSTALS SÖREN BARTELS AND ALEXANDER RAISCH Abstract. We propose a practical finite element
More informationWYSE Academic Challenge 2004 Sectional Mathematics Solution Set
WYSE Academic Callenge 00 Sectional Matematics Solution Set. Answer: B. Since te equation can be written in te form x + y, we ave a major 5 semi-axis of lengt 5 and minor semi-axis of lengt. Tis means
More informationFinite Element Methods for Linear Elasticity
Finite Element Metods for Linear Elasticity Ricard S. Falk Department of Matematics - Hill Center Rutgers, Te State University of New Jersey 110 Frelinguysen Rd., Piscataway, NJ 08854-8019 falk@mat.rutgers.edu
More informationOrder of Accuracy. ũ h u Ch p, (1)
Order of Accuracy 1 Terminology We consider a numerical approximation of an exact value u. Te approximation depends on a small parameter, wic can be for instance te grid size or time step in a numerical
More informationFourier Type Super Convergence Study on DDGIC and Symmetric DDG Methods
DOI 0.007/s095-07-048- Fourier Type Super Convergence Study on DDGIC and Symmetric DDG Metods Mengping Zang Jue Yan Received: 7 December 06 / Revised: 7 April 07 / Accepted: April 07 Springer Science+Business
More informationNON STANDARD FITTED FINITE DIFFERENCE METHOD FOR SINGULAR PERTURBATION PROBLEMS USING CUBIC SPLINE
Global and Stocastic Analysis Vol. 4 No. 1, January 2017, 1-10 NON STANDARD FITTED FINITE DIFFERENCE METHOD FOR SINGULAR PERTURBATION PROBLEMS USING CUBIC SPLINE K. PHANEENDRA AND E. SIVA PRASAD Abstract.
More informationError estimates for a semi-implicit fully discrete finite element scheme for the mean curvature flow of graphs
Interfaces and Free Boundaries 2, 2000 34 359 Error estimates for a semi-implicit fully discrete finite element sceme for te mean curvature flow of graps KLAUS DECKELNICK Scool of Matematical Sciences,
More informationNumerical analysis of a free piston problem
MATHEMATICAL COMMUNICATIONS 573 Mat. Commun., Vol. 15, No. 2, pp. 573-585 (2010) Numerical analysis of a free piston problem Boris Mua 1 and Zvonimir Tutek 1, 1 Department of Matematics, University of
More informationNumerical Solution to Parabolic PDE Using Implicit Finite Difference Approach
Numerical Solution to arabolic DE Using Implicit Finite Difference Approac Jon Amoa-Mensa, Francis Oene Boateng, Kwame Bonsu Department of Matematics and Statistics, Sunyani Tecnical University, Sunyani,
More informationDual-Primal Isogeometric Tearing and Interconnecting Solvers for large-scale systems of multipatch continuous Galerkin IgA equations
www.oeaw.ac.at Dual-Primal Isogeometric Tearing and Interconnecting Solvers for large-scale systems of multipatc continuous Galerkin IgA equations C. Hofer, U. Langer RICAM-Report 2015-44 www.ricam.oeaw.ac.at
More informationNumerical Differentiation
Numerical Differentiation Finite Difference Formulas for te first derivative (Using Taylor Expansion tecnique) (section 8.3.) Suppose tat f() = g() is a function of te variable, and tat as 0 te function
More informationMathematics 5 Worksheet 11 Geometry, Tangency, and the Derivative
Matematics 5 Workseet 11 Geometry, Tangency, and te Derivative Problem 1. Find te equation of a line wit slope m tat intersects te point (3, 9). Solution. Te equation for a line passing troug a point (x
More informationThe derivative function
Roberto s Notes on Differential Calculus Capter : Definition of derivative Section Te derivative function Wat you need to know already: f is at a point on its grap and ow to compute it. Wat te derivative
More informationA method of Lagrange Galerkin of second order in time. Une méthode de Lagrange Galerkin d ordre deux en temps
A metod of Lagrange Galerkin of second order in time Une métode de Lagrange Galerkin d ordre deux en temps Jocelyn Étienne a a DAMTP, University of Cambridge, Wilberforce Road, Cambridge CB3 0WA, Great-Britain.
More informationON THE CONVERGENCE OF A DUAL-PRIMAL SUBSTRUCTURING METHOD. January 2000 Revised April 2000
ON THE CONVERGENCE OF A DUAL-PRIMAL SUBSTRUCTURING METHOD JAN MANDEL AND RADEK TEZAUR January 2000 Revised April 2000 Abstract In te Dual-Primal FETI metod, introduced by Farat et al [5], te domain is
More informationA First-Order System Approach for Diffusion Equation. I. Second-Order Residual-Distribution Schemes
A First-Order System Approac for Diffusion Equation. I. Second-Order Residual-Distribution Scemes Hiroaki Nisikawa W. M. Keck Foundation Laboratory for Computational Fluid Dynamics, Department of Aerospace
More informationA UNIFORM INF SUP CONDITION WITH APPLICATIONS TO PRECONDITIONING
A UNIFORM INF SUP CONDIION WIH APPLICAIONS O PRECONDIIONING KEN ANDRE MARDAL, JOACHIM SCHÖBERL, AND RAGNAR WINHER Abstract. A uniform inf sup condition related to a parameter dependent Stokes problem is
More informationarxiv: v1 [math.na] 17 Jul 2014
Div First-Order System LL* FOSLL* for Second-Order Elliptic Partial Differential Equations Ziqiang Cai Rob Falgout Sun Zang arxiv:1407.4558v1 [mat.na] 17 Jul 2014 February 13, 2018 Abstract. Te first-order
More informationTHE DISCRETE PLATEAU PROBLEM: CONVERGENCE RESULTS
MATHEMATICS OF COMPUTATION Volume 00, Number 0, Pages 000 000 S 0025-5718XX0000-0 THE DISCRETE PLATEAU PROBLEM: CONVERGENCE RESULTS GERHARD DZIUK AND JOHN E. HUTCHINSON Abstract. We solve te problem of
More informationA SYMMETRIC NODAL CONSERVATIVE FINITE ELEMENT METHOD FOR THE DARCY EQUATION
A SYMMETRIC NODAL CONSERVATIVE FINITE ELEMENT METHOD FOR THE DARCY EQUATION GABRIEL R. BARRENECHEA, LEOPOLDO P. FRANCA 1 2, AND FRÉDÉRIC VALENTIN Abstract. Tis work introduces and analyzes novel stable
More informationUniform estimate of the constant in the strengthened CBS inequality for anisotropic non-conforming FEM systems
Uniform estimate of te constant in te strengtened CBS inequality for anisotropic non-conforming FEM systems R. Blaeta S. Margenov M. Neytceva Version of November 0, 00 Abstract Preconditioners based on
More informationNew Fourth Order Quartic Spline Method for Solving Second Order Boundary Value Problems
MATEMATIKA, 2015, Volume 31, Number 2, 149 157 c UTM Centre for Industrial Applied Matematics New Fourt Order Quartic Spline Metod for Solving Second Order Boundary Value Problems 1 Osama Ala yed, 2 Te
More information2009 Elsevier Science. Reprinted with permission from Elsevier.
P. Hansbo and M. Juntunen. 009. Weakly imposed Diriclet boundary conditions for te Brinkman model of porous media flow. Applied Numerical Matematics, volume 59, number 6, pages 174 189. doi:10.1016/j.apnum.008.07.003.
More informationRobotic manipulation project
Robotic manipulation project Bin Nguyen December 5, 2006 Abstract Tis is te draft report for Robotic Manipulation s class project. Te cosen project aims to understand and implement Kevin Egan s non-convex
More informationRunge-Kutta methods. With orders of Taylor methods yet without derivatives of f (t, y(t))
Runge-Kutta metods Wit orders of Taylor metods yet witout derivatives of f (t, y(t)) First order Taylor expansion in two variables Teorem: Suppose tat f (t, y) and all its partial derivatives are continuous
More informationFINITE ELEMENT APPROXIMATIONS AND THE DIRICHLET PROBLEM FOR SURFACES OF PRESCRIBED MEAN CURVATURE
FINITE ELEMENT APPROXIMATIONS AND THE DIRICHLET PROBLEM FOR SURFACES OF PRESCRIBED MEAN CURVATURE GERHARD DZIUK AND JOHN E. HUTCHINSON Abstract. We give a finite element procedure for te Diriclet Problem
More informationarxiv: v1 [math.na] 29 Jan 2019
Iterated total variation regularization wit finite element metods for reconstruction te source term in elliptic systems Micael Hinze a Tran Nan Tam Quyen b a Department of Matematics, University of Hamburg,
More informationComputers and Mathematics with Applications. A nonlinear weighted least-squares finite element method for Stokes equations
Computers Matematics wit Applications 59 () 5 4 Contents lists available at ScienceDirect Computers Matematics wit Applications journal omepage: www.elsevier.com/locate/camwa A nonlinear weigted least-squares
More informationVolume 29, Issue 3. Existence of competitive equilibrium in economies with multi-member households
Volume 29, Issue 3 Existence of competitive equilibrium in economies wit multi-member ouseolds Noriisa Sato Graduate Scool of Economics, Waseda University Abstract Tis paper focuses on te existence of
More informationA Demonstration of the Advantage of Asymptotic Preserving Schemes over Standard Finite Volume Schemes
A Demonstration of te Advantage of Asymptotic Preserving Scemes over Standard Finite Volume Scemes Jocen Scütz Berict Nr. 366 Juni 213 Key words: conservation laws, asymptotic metods, finite volume metods,
More information3. THE EXCHANGE ECONOMY
Essential Microeconomics -1-3. THE EXCHNGE ECONOMY Pareto efficient allocations 2 Edgewort box analysis 5 Market clearing prices 13 Walrasian Equilibrium 16 Equilibrium and Efficiency 22 First welfare
More informationDepartment of Mathematical Sciences University of South Carolina Aiken Aiken, SC 29801
RESEARCH SUMMARY AND PERSPECTIVES KOFFI B. FADIMBA Department of Matematical Sciences University of Sout Carolina Aiken Aiken, SC 29801 Email: KoffiF@usca.edu 1. Introduction My researc program as focused
More informationMATH1151 Calculus Test S1 v2a
MATH5 Calculus Test 8 S va January 8, 5 Tese solutions were written and typed up by Brendan Trin Please be etical wit tis resource It is for te use of MatSOC members, so do not repost it on oter forums
More informationAn approximation method using approximate approximations
Applicable Analysis: An International Journal Vol. 00, No. 00, September 2005, 1 13 An approximation metod using approximate approximations FRANK MÜLLER and WERNER VARNHORN, University of Kassel, Germany,
More informationN igerian Journal of M athematics and Applications V olume 23, (2014), 1 13
N igerian Journal of M atematics and Applications V olume 23, (24), 3 c N ig. J. M at. Appl. ttp : //www.kwsman.com CONSTRUCTION OF POLYNOMIAL BASIS AND ITS APPLICATION TO ORDINARY DIFFERENTIAL EQUATIONS
More informationAN ANALYSIS OF NEW FINITE ELEMENT SPACES FOR MAXWELL S EQUATIONS
Journal of Matematical Sciences: Advances and Applications Volume 5 8 Pages -9 Available at ttp://scientificadvances.co.in DOI: ttp://d.doi.org/.864/jmsaa_7975 AN ANALYSIS OF NEW FINITE ELEMENT SPACES
More informationContinuity and Differentiability of the Trigonometric Functions
[Te basis for te following work will be te definition of te trigonometric functions as ratios of te sides of a triangle inscribed in a circle; in particular, te sine of an angle will be defined to be te
More informationarxiv: v1 [math.na] 28 Apr 2017
THE SCOTT-VOGELIUS FINITE ELEMENTS REVISITED JOHNNY GUZMÁN AND L RIDGWAY SCOTT arxiv:170500020v1 [matna] 28 Apr 2017 Abstract We prove tat te Scott-Vogelius finite elements are inf-sup stable on sape-regular
More informationReflection of electromagnetic waves from magnetic having the ferromagnetic spiral
Reflection of electromagnetic waves from magnetic aving te ferromagnetic spiral Igor V. Bycov 1a Dmitry A. Kuzmin 1b and Vladimir G. Savrov 3 1 Celyabins State University 51 Celyabins Br. Kasiriny Street
More informationThe Laplace equation, cylindrically or spherically symmetric case
Numerisce Metoden II, 7 4, und Übungen, 7 5 Course Notes, Summer Term 7 Some material and exercises Te Laplace equation, cylindrically or sperically symmetric case Electric and gravitational potential,
More information