A finite element approximation for the quasi-static Maxwell Landau Lifshitz Gilbert equations
|
|
- Deirdre Parker
- 5 years ago
- Views:
Transcription
1 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 July 2013) Abstract Te quasi-static Maxwell Landau Lifsitz Gilbert equations wic describe te electromagnetic beaviour of a ferromagnetic material are igly nonlinear. Sopisticated numerical scemes are required to solve te equations, given teir nonlinearity and te constraint tat te solution stays on a spere. We propose an implicit finite element solution to te problem. Te resulting system of algebraic equations is linear wic facilitates te solution process compared to nonlinear metods. We present numerical results to sow te efficacy of te proposed metod. ttp://journal.austms.org.au/ojs/index.pp/anziamj/article/view/6318 gives tis article, c Austral. Matematical Soc Publised December 24, 2013, as part of te Proceedings of te 16t Biennial Computational Tecniques and Applications Conference. issn (Print two pages per seet of paper.) Copies of tis article must not be made oterwise available on te internet; instead link directly to tis url for tis article.
2 Contents C682 Contents 1 Introduction C682 2 A variational formulation of te MLLG equations C685 3 Te finite element sceme C A basis for W (j) C A basis for Y C691 4 Numerical experiments C692 References C695 1 Introduction Te Maxwell Landau Lifsitz Gilbert (mllg) equations describe te electromagnetic beaviour of a ferromagnetic material. For simplicity, we assume tat tere is a bounded cavity D R 3 (wit perfectly conducting outer surface D) into wic a ferromagnet D R 3 is embeded. We furter assume tat D\ D is a vacuum. Over time period (0, T) we let D T := (0, T) D and D T := (0, T) D, and let S 2 be te unit spere. We denote te unit vector of te magnetisation by m(t, x) : D T S 2, and te magnetic field by H(t, x) : DT R 3 over time t and space x. Te quasi-static mllg system is m t = λ 1 m H eff λ 2 m (m H eff ) in D T, (1a) µ 0 m t = µ 0 H t + σ ( H) in D T, (1b) in wic te subscript t indicates a partial derivative wit respect to time, λ 1 0, λ 2 > 0, σ 0 and µ 0 > 0 are constants and H eff is te effective magnetic field wic is dependent on bot H and m. Here m : DT R 3 is
3 1 Introduction C683 te zero extension of m onto D T, tat is { m(t, x) (t, x) D T, m(t, x) = 0 (t, x) D T \D T. Te system (1a) (1b) is supplemented wit initial conditions and boundary conditions m(0,.) = m 0 in D and H(0,.) = H 0 in D, (2) m n = 0 on D T and ( H) n = 0 on D T, (3) were n is te outward normal vector to te relevant surface. Equation (1a) is te first dynamical model for te precessional motion of te magnetisation, suggested by Landau and Lifsitz [8] in In tis model, te time derivative of te magnetisation m is a combination of te precessional movement m H eff and te dissipative movement m (m H eff ) ; see Figure 1. Cimrák [4] sowed existence and uniqueness of a local strong solution of (1a) (3). He also proposed a finite element metod to approximate tis local solution and provided an error estimation [3]. Baňas, Bartels and Prol [2] proposed an implicit nonlinear sceme using te finite element metod, and proved tat te approximate solution converges to a weak global solution. Teir metod required te condition k = O( 2 ) on te time step k and space step for te convergence of te nonlinear system of equations resulting from te discretisation. We propose an implicit linear finite element sceme to find a weak global solution to (1a) (3). Tis approac was initially developed by Alouges and Jaison [1] for te single Landau Lifsitz equation (1a). We extend teir approac to te system (1a) (1b). Te advantage of tis approac is tat tere is no condition imposed on te time step and te space step. For simplicity we coose te effective field H eff = m + H. We focus on implementation issues of te metod. In particular, we sow ow te
4 1 Introduction C684 Figure 1: Te combination of te precessional movement m H eff and te dissipative movement m (m H eff ). finite element spaces and teir bases are constructed. In anoter article we conducted a convergence analysis [9]. In Section 2 we rewrite (1a) in a form wic is more suitable for our approac; see (4). We ten introduce a variational formulation of te mllg system. Section 3 is devoted to te presentation of te implicit linear finite element sceme. Numerical experiments are presented in te last section.
5 2 A variational formulation of te MLLG equations C685 2 A variational formulation of te MLLG equations Before presenting a variational formulation for te mllg equations, we sow ow to rewrite (1a) in te form introduced by Gilbert [6] in wic µ = λ λ2 2. λ 1 m t + λ 2 m m t = µm H eff, (4) Lemma 1. Equation (1a) and equation (4) are equivalent. Proof: By using te elementary identity a (b c) = (a c)b (a b)c, for all a, b, c R 3, and te property m = 1 we obtain m [m (m H eff )] = m H eff. (5) Assume tat m is a solution to (1a). Multiplying bot side of (1a) by λ 2 m using te vector product and using (5) we obtain λ 2 m m t = λ 1 λ 2 m (m H eff ) + λ 2 2m H eff. Multiplying bot sides of (1a) by λ 1 and adding te resulting equation to te above equation, we deduce tat m satisfies (4). Now assume tat m is a solution to (4). On multiplying bot sides of tis equation by λ 2 m using te vector product and noting m (m m t ) = (m m t )m m 2 m t = m t, (6) we deduce λ 1 λ 2 m m t λ 2 2m t = λ 2 µm (m H eff ).
6 3 Te finite element sceme C686 Subtracting tis equation from equation (4) times λ 1 and dividing bot sides of te resulting equation by µ, we obtain (1a). Tis proves te lemma. Before presenting te variational form of tis problem, it is necessary to introduce te function spaces { } H 1 (D, R 3 ) = u L 2 (D, R 3 u ) L 2 (D, R 3 ) for i = 1, 2, 3, x i H(curl; D) = { u L 2 ( D, R 3 ) u L 2 ( D, R 3 ) }. Here L 2 (Ω, R 3 ) is te usual space of Lebesgue integrable functions defined on Ω and taking values in R 3. Following Lemma 1, instead of solving (1a) (3) we solve (1b) (4). A variational form of tis problem is as follows. For all φ C (D T, R 3 ) and ζ C ( D T, R 3 ), find m H 1 (D T ) and H L 2 ( D T ) suc tat H t L 2 ( D T ) and H L 2 ( D T ) to satisfy te Landau Lifsitz Gilbert (llg) equation λ 1 m t φ dx dt + λ 2 (m m t ) φ dx dt D T D T = µ m (m φ) dx dt + µ (m H) φ dx dt, (7) D T D T and Maxwell s equation µ 0 H t ζ dx dt + σ H ζ dx dt = µ 0 m t ζ dx dt. (8) D T D T D T In te following section we introduce a finite element sceme to approximate te solution (m, H) of (7) (8). 3 Te finite element sceme Let T be a regular tetraedrisation of te domain D into a tetraedra of maximal mes size, and let T D be te restriction of T to te domain D.
7 3 Te finite element sceme C687 Te set of N vertices is N := {x 1,..., x N } and te set of M edges is M := {e 1,..., e M }. To discretise te llg equation (7) we introduce te finite element space V of all continuous piecewise linear functions on T D, wic is a subspace of H 1 (D, R 3 ). A basis for V is cosen to be (φ n ) 1 n N, were φ n (x m ) = δ n,m, and δ n,m is te Kronecker delta. Te interpolation operator from C 0 (D, R 3 ) onto V is N I V (v) = v(x n )φ n for all v C 0 (D, R 3 ). n=1 To discretise Maxwell s equation (8), we use te space Y of lowest order edge elements of Nedelec s first family [10] wic is a subspace of H(curl; D). For a function u wic is Lebesgue integrable on all edges in M, we define [10] te interpolation I Y onto Y as I Y (u) = M u q ψ q for all u C 0 ( D, R 3 ), q=1 were u q = u τ q ds, e q in wic τ q is te unit vector in te direction of edge e q. Fixing a positive integer J, we coose te time step k = T/J, and define t j = jk for j = 0,, J. Te functions m(t j, ) and H(t j, ) are approximated by m (j) V and H (j) Y, respectively, for j = 1, 2,..., J. Since m t (t j, ) m(t j+1, ) m(t j, ) k m(j+1) k m (j),
8 3 Te finite element sceme C688 we evaluate m (j+1) from m (j) using were v (j) is an approximation of m t (t j, ). However, to maintain te con- = 1, we normalise te rigt and side of (9) and terefore belonging to V by dition m (j+1) define m (j+1) m (j+1) := m (j) + kv(j), (9) ( m (j+1) m (j) = I + kv(j) V m (j) + kv(j) ) = N m (j) n=1 m (j) Hence it suffices to propose a sceme to compute v (j). Motivated by m t m = 0, we find v (j) Given m (j) (x n) + kv (j) (x n) (x n) + kv (j) (x φ n. n) in te space W(j) defined by }. (10) W (j) := {w V w(x n ) m j (x n) = 0, n = 1,..., N V, if (7) is used to compute v (j), an approximation to m t(t j, ), ten a different test space from V is required because te test function φ in (7) is not perpendicular to m, unlike m t. To circumvent tis difficulty, we use (6) to rewrite (7) as λ 2 m t w dx dt λ 1 (m m t ) w dx dt D T D T = µ m w dx dt + µ H w dx dt (11) D T D T were w = m φ. Now bot m t and w are perpendicular to m for all (t, x) D T. Hence, given m (j) V and H (j) Y, we compute te approximations v (j) and H (j+1) of m t (t j, ) and H(t j+1, ), respectively, as follows. Find v (j) W(j) and H(j+1) Y satisfying, for all w (j) W(j) and
9 3 Te finite element sceme C689 ζ Y, λ 2 and Here, v (j) D = µ D µ 0 w(j) dx λ 1 ( m (j) D + θkv(j) (m (j) ) v(j) ) w(j) dx w (j) dx + µ H (j+1/2) D w (j) dx dt, (12) d t H (j+1) ξ dx dt + σ H (j+1/2) ξ dx dt D D = µ 0 H (j+1/2) v (j) ξ dx dt. (13) D := H(j+1) 2 + H (j) and d t H (j+1) := k 1 (H (j+1) H (j) ). Te parameter θ is arbitrarily cosen to be in [0, 1]. Te metod is explicit wen θ = 0 and fully implicit wen θ = 1. Te algoritm for te numerical approximation of te mllg system is summerised in Algoritm 1. By te Lax Milgram teorem, for eac j > 0 tere exists a unique solution (v (j), Hj+1 ) W(j) Y of equations (12) (13). Since m (0) (x n) = 1 and v (j) (x n) m (j) (x n) = 0 for all n = 1,..., N, and j = 0,..., J, by induction, m (j) (x n) + kv (j) (x n) 1 and m (j) (x n) = 1. Terefore, te algoritm is well defined. We now comment on te construction of basis functions for W (j) wic are necessary in solving (12) (13) in Step 5 of Algoritm 1. and Y
10 3 Te finite element sceme C690 Algoritm 1: Numerical approximation of te mllg system 1 begin 2 Set j = 0 ; 3 Coose m 0 = I V m 0 and H 0 = I Y H 0 ; 4 for j = 0, 1, 2,..., J 1 do 5 Solve (12) and (13) to obtain (v (j), H(j+1) ) W (j) Y ; 6 Define 7 end 8 end m (j+1) (x) := N m (j) n=1 m (j) (x n) + kv (j) (x n) (x n) + kv (j) (x φ n (x) ; n) 3.1 A basis for W (j) From (10), te basis functions of W (j) at eac iteration depend on te solution m (j) computed in te previous iteration. Terefore tey must be computed again for eac iteration of j in Algoritm 1. For eac w W (j), let α n = w(x n ) n = 1,, N, and let α (1) n and α (2) n be two basis vectors of te plane tangential to te vector m (j) (x n) = (u 1, u 2, u 3 ) T R 3. It follows from (10) tat α n = β n α (1) n + γ n α (2) n, for some real numbers β n and γ n. In our computation we take α (1) n = Am (j) (x n) and α (2) n = m (j) (x n) α (1) n,
11 3 Te finite element sceme C691 were u 2 u A = 1 u 2 u u 2 u 3 A basis for W (j) is defined from bases of te 2D planes wic are perpendicular to m (j) (x n) for all n = 1,..., N. Terefore dim(w (j) ) = 2N and w is expressed in terms of α (1) n and α (2) n by It can be sown tat space W (j). w(x) = N n=1 ( ) βn α (1) n + γ n α (2) n φn (x). (14) { } (α (1) n φ n, α (2) n φ n ) 1 n N is a basis for te vector 3.2 A basis for Y A basis {ψ 1,..., ψ M } of Y is defined as follows [10, Section 5.5.1]. Consider an edge e q, q = 1,..., M, and let K be te tetraedron aving e q as an edge. Let λ (1) q and λ (2) q be te barycentric coordinate functions corresponding to te endpoints of e q. We define Te relation ψ q K := λ (1) q λ (2) q λ (2) q λ (1) q. ( ψ q K ) = 2 λ (1) q λ (2) q, is useful in te computation of ( ψ q K ).
12 4 Numerical experiments C692 4 Numerical experiments In order to carry out pysically relevant experiments, te initial conditions of te mllg equations sould be cosen to satisfy te divergence-free constraint [7] div(h 0 + χ D m 0 ) = 0 in D, were χ D is te caracteristic function of D. Tis is acieved by taking H 0 = H 0 χ D m 0, (15) were H 0 is some function defined on D satisfying div H 0 = 0. In our experiment, for simplicity, we coose H 0 to be a constant. We solve te standard problem #4 proposed by te Micromagnetic Modeling Activity Group at te National Institute of Standards and Tecnology [5]. In tis model, te initial conditions m 0 and H 0, and te effective field H eff are m 0 = (1, 0, 0) in D, H 0 = (0.01, 0.01, 0.01) in D, H eff = 2A µ 0 M m+h. s Te parameters λ 1 = γ and λ 1 + α 2 2 = γα 1 + α, 2 were te positive pysical constants are te damping parameter α, te gyromagnetic ratio γ, te vacuum permeability µ 0, te excange constant A, and te magnitude of magnetisation M s. Te values of te pysical constants are α = 1, σ = S 1 m s 1, µ 0 = H m 1 s 1, A = J m 1, µ 0 = H m 1, γ = m A 1 s 1, M s = A m 1. Te domains D and D are cosen to be D = (0, 0.5) (0, 0.125) (0, 0.003),
13 4 Numerical experiments C Figure 2: Te magnetisation domain D (in red) is a tin film. Te magnetic domain D is in blue. and D = ( 0.583, 1.083) ( 0.146, 0.271) ( 0.767, 0.770), wit distances measured in µm; see Figure 2. Te domain D is uniformly partitioned into cubes of dimensions µm 3, were eac cube consists of six tetraedra. We generate a nonuniform mes for te magnetic domain D in suc a way tat it is identical to te mes for D in te region near D, and te mes size gradually increases away from D. A cross section of te mes at x 3 = 0 is displayed in Figure 3. At eac iteration of Algoritm 1, te system to be solved is linear and of
14 4 Numerical experiments C694 Table 1: A comparision between Baňas, Bartels and Prol s (bbp) metod and our metod. bbp Our metod Discrete system nonlinear linear (uses fixed-point iteration) (solved directly) Degrees of freedom 3N + M 2N + M Required condition k = O( 2 ) None Basis functions same different of solution space for all iterations in eac iteration size 2N + M. A comparision of our metod and te metod proposed by Baňas, Bartels and Prol [2] is presented in Table 1. For j = 0, 1, 2,..., let E (j) T E (j) T := 2σA M s D := E (j) ex + E (j) H + E(j) E, be te total energy at time t j = jk defined by H (j) 2 dx + λ 2 H (j) D µ 2 dx D m (j) 2 dx + 2σµ 0 were E (j) ex, E (j) H and E(j) E are te excange energy, magnetic field energy and electric field energy, respectively. Our computation sows tat te total energy decreases, tat is E (j+1) T E (j) T for all j 0 ; (16) see Figure 4. Te decrease of te discrete energy E (j) T, j = 1,..., J, suggests te gradient stability of te mllg solutions. In Figure 4 we plot different scaled versions of te different energies versus log t. Te cange in E (j) T dominated by te cange in E (j) H. Figure 4 also sows tat te sequence { m (j) } j 0 is bounded in L 2 (D), and te sequences {H (j) } j 0 and { H (j) } j 0 are bounded in L 2 ( D). In a fortcoming paper, we will sow tat our numerical solution converges to is
15 References C695 Figure 3: Mes for te domain D at x 3 = 0. a weak solution of te problem (1b) (4) by proving tat te inequality (16) olds. Acknowledgements Te autors acknowledge financial support troug te Australian ARC project DP References [1] F. Alouges. A new finite element sceme for Landau Lifcitz equations. Discrete Contin. Dyn. Syst. Ser. S, 1: (2008). doi: /dcdss C683 [2] L. Baňas, S. Bartels, and A. Prol. A convergent implicit finite element discretization of te Maxwell Landau Lifsitz Gilbert equation. SIAM
16 References C E energy H Energy excange energy total energy 1 Energy (J) log(t) (ns) Figure 4: Evolution of te energies, log(t) E E (t)/600, E H (t)/10 6, 10 4 E ex (t), E T (t)/10 6.
17 References C697 J. Numer. Anal., 46: (2008). doi: / C683, C694 [3] I. Cimrák. Error analysis of a numerical sceme for 3D Maxwell Landau Lifsitz system. Mat. Metods Appl. Sci., 30: (2008). doi: /mma.863 C683 [4] I. Cimrák. Existence, regularity and local uniqueness of te solutions to te Maxwell Landau Lifsitz system in tree dimensions. J. Mat. Anal. Appl., 329: (2007). doi: /j.jmaa C683 [5] Micromagnetic Modeling Activity Group, Center for Teoretical and Computational Materials Science, National Institute of Standards and Tecnology, USA. ttp:// C692 [6] T. L. Gilbert. A Lagrangian formulation of te gyromagnetic equation of te magnetic field. Pys Rev, 100: (1955). Abstract only; full report, Armor Researc Foundation Project No. A059, Supplementary Report, May 1, 1956 (unpublised). C685 [7] B. Guo and S. Ding. Landau Lifsitz Equations, volume 1 of Frontiers of Researc wit te Cinese Academy of Sciences. World Scientific, Hackensack NJ (2008). doi: / C692 [8] L. Landau and E. Lifscitz. On te teory of te dispersion of magnetic permeability in ferromagnetic bodies. Pys Z Sowjetunion, 8: (1935); Perspectives in Teoretical Pysics, pp , Pergamon, Amsterdam (1992). doi: /b C683 [9] K.-N. Le and T. Tran. A convergent finite element approximation for te quasi-static Maxwell Landau Lifsitz Gilbert equations. Researc Report, arxiv: , 2012, submitted to Computers & Matematics wit Applications. ttp://arxiv.org/abs/ C684
18 References C698 [10] P. Monk. Finite Element Metods for Maxwell s equations. Numerical Matematics and Scientific Computation. Oxford University Press, New York (2003). C687, C691 Autor addresses 1. Kim-Ngan Le, Scool of Matematics and Statistics, Te University of New Sout Wales, Sydney 2052, Australia. mailto:n.le-kim@student.unsw.edu.au 2. T. Tran, Scool of Matematics and Statistics, Te University of New Sout Wales, Sydney 2052, Australia. mailto:tan.tran@unsw.edu.au
CONVERGENCE 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 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] 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 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 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 informationWEAK MARTINGALE SOLUTIONS TO THE STOCHASTIC LANDAU LIFSHITZ GILBERT EQUATION WITH MULTI-DIMENSIONAL NOISE VIA A CONVERGENT FINITE-ELEMENT SCHEME
WEAK MARTINGALE SOLUTIONS TO THE STOCHASTIC LANDAU LIFSHITZ GILBERT EQUATION WITH MULTI-DIMENSIONAL NOISE VIA A CONVERGENT FINITE-ELEMENT SCHEME BENIAMIN GOLDYS, JOSEPH GROTOWSKI, AND KIM-NGAN LE Abstract.
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 informationDedicated to the 70th birthday of Professor Lin Qun
Journal of Computational Matematics, Vol.4, No.3, 6, 4 44. ACCELERATION METHODS OF NONLINEAR ITERATION FOR NONLINEAR PARABOLIC EQUATIONS Guang-wei Yuan Xu-deng Hang Laboratory of Computational Pysics,
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 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 information= 0 and states ''hence there is a stationary point'' All aspects of the proof dx must be correct (c)
Paper 1: Pure Matematics 1 Mark Sceme 1(a) (i) (ii) d d y 3 1x 4x x M1 A1 d y dx 1.1b 1.1b 36x 48x A1ft 1.1b Substitutes x = into teir dx (3) 3 1 4 Sows d y 0 and states ''ence tere is a stationary point''
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 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 informationch (for some fixed positive number c) reaching c
GSTF Journal of Matematics Statistics and Operations Researc (JMSOR) Vol. No. September 05 DOI 0.60/s4086-05-000-z Nonlinear Piecewise-defined Difference Equations wit Reciprocal and Cubic Terms Ramadan
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 informationThe entransy dissipation minimization principle under given heat duty and heat transfer area conditions
Article Engineering Termopysics July 2011 Vol.56 No.19: 2071 2076 doi: 10.1007/s11434-010-4189-x SPECIAL TOPICS: Te entransy dissipation minimization principle under given eat duty and eat transfer area
More informationCONVERGENCE ANALYSIS OF YEE SCHEMES FOR MAXWELL S EQUATIONS IN DEBYE AND LORENTZ DISPERSIVE MEDIA
INTRNATIONAL JOURNAL OF NUMRICAL ANALYSIS AND MODLING Volume XX Number 0 ages 45 c 03 Institute for Scientific Computing and Information CONVRGNC ANALYSIS OF Y SCHMS FOR MAXWLL S QUATIONS IN DBY AND LORNTZ
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 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 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 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 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 informationERROR BOUNDS FOR FINITE-DIFFERENCE METHODS FOR RUDIN OSHER FATEMI IMAGE SMOOTHING
ERROR BOUNDS FOR FINITE-DIFFERENCE METHODS FOR RUDIN OSHER FATEMI IMAGE SMOOTHING JINGYUE WANG AND BRADLEY J. LUCIER Abstract. We bound te difference between te solution to te continuous Rudin Oser Fatemi
More informationA = h w (1) Error Analysis Physics 141
Introduction In all brances of pysical science and engineering one deals constantly wit numbers wic results more or less directly from experimental observations. Experimental observations always ave inaccuracies.
More informationCONVERGENCE ANALYSIS OF YEE SCHEMES FOR MAXWELL S EQUATIONS IN DEBYE AND LORENTZ DISPERSIVE MEDIA
INTRNATIONAL JOURNAL OF NUMRICAL ANALYSIS AND MODLING Volume Number 4 ages 657 687 c 04 Institute for Scientific Computing and Information CONVRGNC ANALYSIS OF Y SCHMS FOR MAXWLL S QUATIONS IN DBY AND
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 informationFunction Composition and Chain Rules
Function Composition and s James K. Peterson Department of Biological Sciences and Department of Matematical Sciences Clemson University Marc 8, 2017 Outline 1 Function Composition and Continuity 2 Function
More informationSECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY
(Section 3.2: Derivative Functions and Differentiability) 3.2.1 SECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY LEARNING OBJECTIVES Know, understand, and apply te Limit Definition of te Derivative
More informationLIMITS AND DERIVATIVES CONDITIONS FOR THE EXISTENCE OF A LIMIT
LIMITS AND DERIVATIVES Te limit of a function is defined as te value of y tat te curve approaces, as x approaces a particular value. Te limit of f (x) as x approaces a is written as f (x) approaces, as
More informationlecture 26: Richardson extrapolation
43 lecture 26: Ricardson extrapolation 35 Ricardson extrapolation, Romberg integration Trougout numerical analysis, one encounters procedures tat apply some simple approximation (eg, linear interpolation)
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 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 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 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 informationSOLUTIONS OF FOURTH-ORDER PARTIAL DIFFERENTIAL EQUATIONS IN A NOISE REMOVAL MODEL
Electronic Journal of Differential Equations, Vol. 7(7, No., pp.. ISSN: 7-669. URL: ttp://ejde.mat.txstate.edu or ttp://ejde.mat.unt.edu ftp ejde.mat.txstate.edu (login: ftp SOLUTIONS OF FOURTH-ORDER PARTIAL
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 informationDownloaded 11/15/17 to Redistribution subject to SIAM license or copyright; see
SIAM J. NUMER. ANAL. Vol. 55, No. 6, pp. 2787 2810 c 2017 Society for Industrial and Applied Matematics EDGE ELEMENT METHOD FOR OPTIMAL CONTROL OF STATIONARY MAXWELL SYSTEM WITH GAUSS LAW IRWIN YOUSEPT
More informationHow to Find the Derivative of a Function: Calculus 1
Introduction How to Find te Derivative of a Function: Calculus 1 Calculus is not an easy matematics course Te fact tat you ave enrolled in suc a difficult subject indicates tat you are interested in te
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 informationarxiv: v1 [math.dg] 4 Feb 2015
CENTROID OF TRIANGLES ASSOCIATED WITH A CURVE arxiv:1502.01205v1 [mat.dg] 4 Feb 2015 Dong-Soo Kim and Dong Seo Kim Abstract. Arcimedes sowed tat te area between a parabola and any cord AB on te parabola
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 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 informationSuperconvergence of energy-conserving discontinuous Galerkin methods for. linear hyperbolic equations. Abstract
Superconvergence of energy-conserving discontinuous Galerkin metods for linear yperbolic equations Yong Liu, Ci-Wang Su and Mengping Zang 3 Abstract In tis paper, we study superconvergence properties of
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 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 informationDiscontinuous Galerkin Methods for Relativistic Vlasov-Maxwell System
Discontinuous Galerkin Metods for Relativistic Vlasov-Maxwell System He Yang and Fengyan Li December 1, 16 Abstract e relativistic Vlasov-Maxwell (RVM) system is a kinetic model tat describes te dynamics
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 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 informationNUMERICAL DIFFERENTIATION. James T. Smith San Francisco State University. In calculus classes, you compute derivatives algebraically: for example,
NUMERICAL DIFFERENTIATION James T Smit San Francisco State University In calculus classes, you compute derivatives algebraically: for example, f( x) = x + x f ( x) = x x Tis tecnique requires your knowing
More information1 Calculus. 1.1 Gradients and the Derivative. Q f(x+h) f(x)
Calculus. Gradients and te Derivative Q f(x+) δy P T δx R f(x) 0 x x+ Let P (x, f(x)) and Q(x+, f(x+)) denote two points on te curve of te function y = f(x) and let R denote te point of intersection of
More informationMath 31A Discussion Notes Week 4 October 20 and October 22, 2015
Mat 3A Discussion Notes Week 4 October 20 and October 22, 205 To prepare for te first midterm, we ll spend tis week working eamples resembling te various problems you ve seen so far tis term. In tese notes
More informationHigh-Order Energy and Linear Momentum Conserving Methods for the Klein-Gordon Equation
matematics Article Hig-Order Energy and Linear Momentum Conserving Metods for te Klein-Gordon Equation He Yang Department of Matematics, Augusta University, Augusta, GA 39, USA; yang@augusta.edu; Tel.:
More informationLinearized Primal-Dual Methods for Linear Inverse Problems with Total Variation Regularization and Finite Element Discretization
Linearized Primal-Dual Metods for Linear Inverse Problems wit Total Variation Regularization and Finite Element Discretization WENYI TIAN XIAOMING YUAN September 2, 26 Abstract. Linear inverse problems
More informationParametric Spline Method for Solving Bratu s Problem
ISSN 749-3889 print, 749-3897 online International Journal of Nonlinear Science Vol4202 No,pp3-0 Parametric Spline Metod for Solving Bratu s Problem M Zarebnia, Z Sarvari 2,2 Department of Matematics,
More informationSmoothness of solutions with respect to multi-strip integral boundary conditions for nth order ordinary differential equations
396 Nonlinear Analysis: Modelling and Control, 2014, Vol. 19, No. 3, 396 412 ttp://dx.doi.org/10.15388/na.2014.3.6 Smootness of solutions wit respect to multi-strip integral boundary conditions for nt
More informationMath 161 (33) - Final exam
Name: Id #: Mat 161 (33) - Final exam Fall Quarter 2015 Wednesday December 9, 2015-10:30am to 12:30am Instructions: Prob. Points Score possible 1 25 2 25 3 25 4 25 TOTAL 75 (BEST 3) Read eac problem carefully.
More informationLecture 15. Interpolation II. 2 Piecewise polynomial interpolation Hermite splines
Lecture 5 Interpolation II Introduction In te previous lecture we focused primarily on polynomial interpolation of a set of n points. A difficulty we observed is tat wen n is large, our polynomial as to
More informationIntegral Calculus, dealing with areas and volumes, and approximate areas under and between curves.
Calculus can be divided into two ke areas: Differential Calculus dealing wit its, rates of cange, tangents and normals to curves, curve sketcing, and applications to maima and minima problems Integral
More informationDerivation Of The Schwarzschild Radius Without General Relativity
Derivation Of Te Scwarzscild Radius Witout General Relativity In tis paper I present an alternative metod of deriving te Scwarzscild radius of a black ole. Te metod uses tree of te Planck units formulas:
More informationDefinition of the Derivative
Te Limit Definition of te Derivative Tis Handout will: Define te limit grapically and algebraically Discuss, in detail, specific features of te definition of te derivative Provide a general strategy of
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 informationOverlapping domain decomposition methods for elliptic quasi-variational inequalities related to impulse control problem with mixed boundary conditions
Proc. Indian Acad. Sci. (Mat. Sci.) Vol. 121, No. 4, November 2011, pp. 481 493. c Indian Academy of Sciences Overlapping domain decomposition metods for elliptic quasi-variational inequalities related
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 informationAn Effective Integrator for the Landau-Lifshitz-Gilbert Equation
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
More informationEXISTENCE OF SOLUTIONS FOR A CLASS OF VARIATIONAL INEQUALITIES
Journal of Matematics and Statistics 9 (4: 35-39, 3 ISSN: 549-3644 3 doi:.3844/jmssp.3.35.39 Publised Online 9 (4 3 (ttp://www.tescipub.com/jmss.toc EXISTENCE OF SOLUTIONS FOR A CLASS OF ARIATIONAL INEQUALITIES
More informationA MONTE CARLO ANALYSIS OF THE EFFECTS OF COVARIANCE ON PROPAGATED UNCERTAINTIES
A MONTE CARLO ANALYSIS OF THE EFFECTS OF COVARIANCE ON PROPAGATED UNCERTAINTIES Ronald Ainswort Hart Scientific, American Fork UT, USA ABSTRACT Reports of calibration typically provide total combined uncertainties
More informationc 2011 Society for Industrial and Applied Mathematics Key words. total variation, variational methods, finite-difference methods, error bounds
SIAM J. NUMER. ANAL. Vol. 49, No. 2, pp. 845 868 c 211 Society for Industrial and Applied Matematics ERROR BOUNDS FOR FINITE-DIFFERENCE METHODS FOR RUDIN OSHER FATEMI IMAGE SMOOTHING JINGYUE WANG AND BRADLEY
More informationSymmetry Labeling of Molecular Energies
Capter 7. Symmetry Labeling of Molecular Energies Notes: Most of te material presented in tis capter is taken from Bunker and Jensen 1998, Cap. 6, and Bunker and Jensen 2005, Cap. 7. 7.1 Hamiltonian Symmetry
More informationParallel algorithm for convection diffusion system based on least squares procedure
Zang et al. SpringerPlus (26 5:69 DOI.86/s464-6-3333-8 RESEARCH Parallel algoritm for convection diffusion system based on least squares procedure Open Access Jiansong Zang *, Hui Guo 2, Hongfei Fu 3 and
More informationPre-Calculus Review Preemptive Strike
Pre-Calculus Review Preemptive Strike Attaced are some notes and one assignment wit tree parts. Tese are due on te day tat we start te pre-calculus review. I strongly suggest reading troug te notes torougly
More informationDepartment of Statistics & Operations Research, Aligarh Muslim University, Aligarh, India
Open Journal of Optimization, 04, 3, 68-78 Publised Online December 04 in SciRes. ttp://www.scirp.org/ournal/oop ttp://dx.doi.org/0.436/oop.04.34007 Compromise Allocation for Combined Ratio Estimates of
More informationINTRODUCTION TO CALCULUS LIMITS
Calculus can be divided into two ke areas: INTRODUCTION TO CALCULUS Differential Calculus dealing wit its, rates of cange, tangents and normals to curves, curve sketcing, and applications to maima and
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 informationChapter 2. Limits and Continuity 16( ) 16( 9) = = 001. Section 2.1 Rates of Change and Limits (pp ) Quick Review 2.1
Capter Limits and Continuity Section. Rates of Cange and Limits (pp. 969) Quick Review..... f ( ) ( ) ( ) 0 ( ) f ( ) f ( ) sin π sin π 0 f ( ). < < < 6. < c c < < c 7. < < < < < 8. 9. 0. c < d d < c
More information1 The concept of limits (p.217 p.229, p.242 p.249, p.255 p.256) 1.1 Limits Consider the function determined by the formula 3. x since at this point
MA00 Capter 6 Calculus and Basic Linear Algebra I Limits, Continuity and Differentiability Te concept of its (p.7 p.9, p.4 p.49, p.55 p.56). Limits Consider te function determined by te formula f Note
More informationSolutions to the Multivariable Calculus and Linear Algebra problems on the Comprehensive Examination of January 31, 2014
Solutions to te Multivariable Calculus and Linear Algebra problems on te Compreensive Examination of January 3, 24 Tere are 9 problems ( points eac, totaling 9 points) on tis portion of te examination.
More informationSection 15.6 Directional Derivatives and the Gradient Vector
Section 15.6 Directional Derivatives and te Gradient Vector Finding rates of cange in different directions Recall tat wen we first started considering derivatives of functions of more tan one variable,
More informationLecture 21. Numerical differentiation. f ( x+h) f ( x) h h
Lecture Numerical differentiation Introduction We can analytically calculate te derivative of any elementary function, so tere migt seem to be no motivation for calculating derivatives numerically. However
More informationReflection Symmetries of q-bernoulli Polynomials
Journal of Nonlinear Matematical Pysics Volume 1, Supplement 1 005, 41 4 Birtday Issue Reflection Symmetries of q-bernoulli Polynomials Boris A KUPERSHMIDT Te University of Tennessee Space Institute Tullaoma,
More informationThe Verlet Algorithm for Molecular Dynamics Simulations
Cemistry 380.37 Fall 2015 Dr. Jean M. Standard November 9, 2015 Te Verlet Algoritm for Molecular Dynamics Simulations Equations of motion For a many-body system consisting of N particles, Newton's classical
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 informationarxiv: v1 [math.ap] 4 Aug 2017
New Two Step Laplace Adam-Basfort Metod for Integer an Non integer Order Partial Differential Equations arxiv:178.1417v1 [mat.ap] 4 Aug 17 Abstract Rodrigue Gnitcogna*, Abdon Atangana** *Department of
More informationArbitrary order exactly divergence-free central discontinuous Galerkin methods for ideal MHD equations
Arbitrary order exactly divergence-free central discontinuous Galerkin metods for ideal MHD equations Fengyan Li, Liwei Xu Department of Matematical Sciences, Rensselaer Polytecnic Institute, Troy, NY
More informationMath 312 Lecture Notes Modeling
Mat 3 Lecture Notes Modeling Warren Weckesser Department of Matematics Colgate University 5 7 January 006 Classifying Matematical Models An Example We consider te following scenario. During a storm, a
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 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 informationMath Spring 2013 Solutions to Assignment # 3 Completion Date: Wednesday May 15, (1/z) 2 (1/z 1) 2 = lim
Mat 311 - Spring 013 Solutions to Assignment # 3 Completion Date: Wednesday May 15, 013 Question 1. [p 56, #10 (a)] 4z Use te teorem of Sec. 17 to sow tat z (z 1) = 4. We ave z 4z (z 1) = z 0 4 (1/z) (1/z
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 informationMore on generalized inverses of partitioned matrices with Banachiewicz-Schur forms
More on generalized inverses of partitioned matrices wit anaciewicz-scur forms Yongge Tian a,, Yosio Takane b a Cina Economics and Management cademy, Central University of Finance and Economics, eijing,
More informationKey words. Finite element method; convection-diffusion-reaction; nonnegativity; boundedness
PRESERVING NONNEGATIVITY OF AN AFFINE FINITE ELEMENT APPROXIMATION FOR A CONVECTION-DIFFUSION-REACTION PROBLEM JAVIER RUIZ-RAMÍREZ Abstract. An affine finite element sceme approximation of a time dependent
More informationQuasiperiodic phenomena in the Van der Pol - Mathieu equation
Quasiperiodic penomena in te Van der Pol - Matieu equation F. Veerman and F. Verulst Department of Matematics, Utrect University P.O. Box 80.010, 3508 TA Utrect Te Neterlands April 8, 009 Abstract Te Van
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 informationTaylor Series and the Mean Value Theorem of Derivatives
1 - Taylor Series and te Mean Value Teorem o Derivatives Te numerical solution o engineering and scientiic problems described by matematical models oten requires solving dierential equations. Dierential
More informationA SPLITTING LEAST-SQUARES MIXED FINITE ELEMENT METHOD FOR ELLIPTIC OPTIMAL CONTROL PROBLEMS
INTERNATIONAL JOURNAL OF NUMERICAL ANALSIS AND MODELING Volume 3, Number 4, Pages 6 626 c 26 Institute for Scientific Computing and Information A SPLITTING LEAST-SQUARES MIED FINITE ELEMENT METHOD FOR
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 informationPre-lab Quiz/PHYS 224 Earth s Magnetic Field. Your name Lab section
Pre-lab Quiz/PHYS 4 Eart s Magnetic Field Your name Lab section 1. Wat do you investigate in tis lab?. For a pair of Helmoltz coils described in tis manual and sown in Figure, r=.15 m, N=13, I =.4 A, wat
More informationOn convergence of the immersed boundary method for elliptic interface problems
On convergence of te immersed boundary metod for elliptic interface problems Zilin Li January 26, 2012 Abstract Peskin s Immersed Boundary (IB) metod is one of te most popular numerical metods for many
More informationDifferentiation. Area of study Unit 2 Calculus
Differentiation 8VCE VCEco Area of stud Unit Calculus coverage In tis ca 8A 8B 8C 8D 8E 8F capter Introduction to limits Limits of discontinuous, rational and brid functions Differentiation using first
More informationGlobal Existence of Classical Solutions for a Class Nonlinear Parabolic Equations
Global Journal of Science Frontier Researc Matematics and Decision Sciences Volume 12 Issue 8 Version 1.0 Type : Double Blind Peer Reviewed International Researc Journal Publiser: Global Journals Inc.
More information