Nedelec elements for computational electromagnetics
|
|
- Peter Farmer
- 6 years ago
- Views:
Transcription
1 Nedelec elements for computational electromagnetics Per Jacobsson, June 5, 2007 Maxwell s equations E = jωµh () H = jωεe + J (2) (εe) = ρ (3) (µh) = 0 (4) J = jωρ (5) Only three of the equations are needed, (), (2) and (3) or (), (2) and (5). Combining () and (2) gives the wave equation (µ E) ω 2 εe = jωj (6) for E and similarly for H. E,J H(curl;Ω) corresponds to finite energy solutions, where H(curl;Ω) = {u [L 2 (Ω)] 3 : u [L 2 (Ω)] 3 (7) in three dimensions. The norm associated with H(curl;Ω) is defined as Some other vector spaces used: where ν is the outward normal to Ω. v H(curl;Ω) = ( v [L2 (Ω)] 3 + v [L 2 (Ω)] 3)/2 (8) H (Ω) = {v L 2 (Ω) : v [L 2 (Ω)] 3 (9) H 0(Ω) = {v H (Ω) : v = 0 on Ω (0) H 0 (curl;ω) = {v H(curl;Ω) : ν v = 0 () H(div;Ω,ε) = {v [L 2 (Ω)] 3 : (εv) L 2 (Ω) (2) H 0 (div;ω,ε) = {v [L 2 (Ω)] 3 : (εv) = 0 (3) V = H 0 (curl;ω) H(div;Ω,ε) (4) 2 Problems when using nodal elements When using nodal elements spurious solutions can occur and we can also get convergence to solutions of a problem different from our original problem. As an example we look at a simple cavity eigenvalue problem with perfectly conducting walls (µ E) = ω 2 εe in Ω (5) (εe) = 0 in Ω (6) ν E = 0 on Ω (7)
2 The weak form of this eigenvalue problem is to find ω R, E H 0 (curl;ω) H 0 (div;ω,ε), E 0 such that (µ E, v) = ω 2 (εe,v) v H 0 (curl;ω) H 0 (div;ω,ε). (8) We see that ω 2 = 0 is not an eigenvalue to (8) since that would imply E = 0. This means that we cannot have a physical mode with non-zero energy and ω = 0, the only physical modes are the ones with ω > 0. For this case, (5) implies (6) and we can write the problem on weak form as Find ω R, E H 0 (curl;ω), E 0 s.t. (µ E, v) = ω 2 (εe,v) v H 0 (curl;ω). (9) This adds an infinite-dimensional space of non-physical eigensolutions that coincides with H 0(Ω) with zero eigenvalues. It is important to recognize these non-physical solutions so that they can be discarded. When using FE we seek solutions in a finite-dimensional subspace V h of H 0 (curl;ω): Find ω h R, E h V h, E h 0 s.t. (µ E h, v) = ω h (εe h,v) v V h. (20) For nodal elements the approximation space is V h = [P k ] 3 H 0 (curl;ω). When using nodal based FE on an unstructured grid, we obtain many spurious eigenvalues approaching zero and it can be hard to recognize these so that they can be discarded. For certain structured grids, the zero eigenvalues can be well approximated but other spurious solutions appear in the spectrum. Since these spurious eigenvalues can be distributed all over the spectrum, they are hard to distinguish from the physical eigenvalues. For non-convex domains, it has been proved that the discrete eigenvalues computed using nodal elements converge to the eigenvalues of a problem different than the original problem. One possible explanation for the failure of nodal based FE in the case of non-convex domains in two dimensions is that the nodal elements are constrained in [H (Ω)] 2, which is a proper and closed subset of V. This means that there are eigenfunctions in V that cannot be approximated by nodal elements in [H (Ω)] 2. Another important drawback of nodal-based FE when solving Maxwell s equations is the inability to model field singularities at conducting corners and tips. The boundary condition for the electric field at a conducting surface is ν E = 0, but enforcing this boundary condition at each node also enforces normal continuity, which is not desired. 3 Nedelec s edge elements We start by introducing a subspace S k of homogeneous vector polynomials of degree k in d dimensions by S k = {p [ P k ] d : x p = 0, (2) where P k is the space of homogeneous polynomials of order k. We can now introduce the space R k which is needed for Nedelec s vector elements, R k = [P k ] d S k. (22) The dimension of the space R k in three dimensions is dim(r k ) = 2 (k + 3)(k + 2)k, so for k = we have dim(r ) = 6. In two dimensions, a polynomial p R has dimension 3 and can be represented as p = (c + c 2 y,c 3 c 2 x). (23) For k =, the local degrees of freedom in two dimensions are given by N i (p) = p τ i ds, (24) e i 2
3 where τ i is the tangential unit vector of edge e i. The local shape functions can be expressed in terms of barycentric coordinates. In two dimensions we can define the local shape function θ associated with edge oriented from node to node 2 as θ = λ λ 2 λ 2 λ, (25) where the barycentric coordinates for node and 2 can be written as λ = a + b x + c y and λ 2 = a 2 + b 2 x + c 2 y respectively, where a i and b i depend on the nodal coordinates of the triangle. Now θ = λ λ 2 λ 2 λ = (a + b x + c y)(b 2,c 2 ) (a 2 + b 2 x + c 2 y)(b,c ) = (d + d 2 y,d 3 d 2 x) R. (26) The local shape functions for triangular elements in two dimensions are shown in Fig.. Figure : The local shape functions for k = in two dimensions on a triangle. 3. Affine maps When transforming to and from the reference element we need to transform in a special way since we are working in H(curl; ˆ) with vectors in R k. For the affine map F ˆx = B ˆx + b and û R k, the transformation is given by u F = (B T ) û (27) and the curl transforms as u = det(b ) B ˆ û. (28) Lemma. The space R k is invariant under the transformation (27) Proof. Since for û R k we can write û = ˆp + ˆp 2, where ˆp [P k ] 3 and ˆp 2 S k. Now, u(x) = [(B T ) ˆp + (B T ) ˆp 2 ](F (x)) = [(B T ) ˆp ](F (x)) + (BT ) ˆp 2 (B x B b ). (29) 3
4 Since ˆp 2 [ P k ] 3, we have that ˆp 2 (B x B b ) = ˆp 2 (B x) + ˆp 3(x), where ˆp 3 [P k ] 3. This means u(x) = [(B) T ˆp (F (x)) + ˆp 3(x)] + (B) T ˆp 2 (B x). (30) We have that (B T ) ˆp F + ˆp 3 [P k ] 3 and also that (B T ) ˆp 2 (B x) x = ˆp 2(B x) (B x) = 0, since ˆp 2 S k. Thus, u R k. Lemma 2. If τ h is a regular mesh and s 0, then we have for v transformed by (27) to give ˆv and ˆv [Hs ( ˆ)] 3 Ch s 2 v [H s ()] 3 (3) ˆ ˆv [H s ( ˆ)] Ch s+ 3 2 v [H s ()] 3. (32) 3.2 Degrees of freedom For the reference element ˆ, which is a tetrahedron defined by the vertices â,...,â 4 given by â = (0,0,0) T, â 2 = (,0,0) T, â 3 = (0,,0) T and â 4 = (0,0,) T, we define the curl-conforming finite element ( ˆ, P ˆ, N ˆ) by ˆ is the reference element P ˆ = R k Three types of degrees of freedom associated with edges ê of ˆ, faces ˆf of ˆ and ˆ itself. The unit vector along edge ê is denoted ˆτ. Mê(û) = û ˆτ ˆqdŝ, ˆq P k (ê), (33) ê { M ˆf(û) = area( ˆf) û ˆqdÂ, ˆq [P k 2( ˆf)] 3 and ˆq ˆν = 0, (34) ˆf M ˆ(û) = û ˆqdˆV, ˆq [P k 3 (. ˆ)] 3 (35) ˆ The total set of degrees of freedom is then N ˆ = Mê(û) M ˆf(û) M ˆ(û). By using the transformation (27) we can define the finite element (, P, N ) for a general tetrahedron as is a tetrahedron P = R k Three types of degrees of freedom associated with edges e of, faces f of and itself. The unit vector along edge e is denoted τ. M e (u) = u τqds, q P k (e), (36) e { M f (u) = u qda, q = B ˆq, ˆq [P k 2 ( area(f) ˆf)] 3 and ˆq ˆν = 0, (37) f M (u) = u qdv, q mapped by q F = det(b ) B ˆq, ˆq [P k 3 ()] 3. The total set of degrees of freedom for a general tetrahedron is then N = M e (u) M f (u) M (u). Lemma 3. If det(b ) 0, τ = B ˆτ B ˆτ and the transformation (27) is used, the degrees of freedom (36) (38) on are identical to those on ˆ. (38) 4
5 3.3 Unisolvence Theorem. If u R k is such that all degrees of freedom (36) (38)vanish, then u = 0. To prove Theorem we need two Lemmas. Lemma 4. If u R k is such that the degrees of freedom (37) vanish on f and the degrees of freedom (36) vanish on all edges of f, then u ν = 0 on f. Lemma 5. If u R k is such that u = 0, then u = p for p P k. Proof of Theorem. Since the degrees of freedom on the general tetrahedron and the reference tetrahedron are identical, we prove the theorem on the reference element. Because of Lemma 4 we know that û ˆν = 0 on ˆ. By integration by parts and using that the degrees of freedom (38) are zero, we have ˆ û ˆqdˆV = û ˆ ˆqdˆV = 0, ˆq [P k 2 ] 3. (39) ˆ ˆ Using Stokes theorem on each face ˆf of ˆ and that the degrees of freedom (37) are zero, we have ˆ ˆf û T ˆqd = û T ˆ ˆf ˆqd = 0, ˆq P k ( ˆf), (40) ˆf ˆf where û T = (ˆν û) ˆν on ˆ. The surface vector curl operator ˆ ˆf is defined as ˆ ˆf q = ˆν ˆ ˆfq. This means that ˆ ˆf û T = (ˆ û) ˆν = 0 on ˆf and hence on ˆ. Because of this and the fact that ˆ û [P k ] 3, we can write ˆ û = (ˆx Ψ, ˆx 2,Ψ 2, ˆx 3 Ψ 3 ) T, (4) where Ψ = (Ψ,Ψ 2,Ψ 3 ) T [P k 2 ] 3. Picking ˆq = Ψ in (39) shows that ˆ û = 0 in ˆ. Using Lemma 5 this means that we can write û = ˆ ˆp for some ˆp P k. The fact that û ˆν = 0 on ˆ implies that we can take ˆp = 0 on ˆ and write ˆp as ˆp = ˆx ˆx 2ˆx 3ˆr, for some ˆr P k 3. (42) The vanishing of the degrees of freedom (38) implies that ˆp = 0 and thus û = The interpolant We can define the global finite element space on a mesh τ h as V h = {u H(curl;Ω) : u R k τ h. (43) The local interpolant r u R k for τ h is characterized by the vanishing of the degrees of freedom on u r u, The global interpolant is then defined by M e (u r u) = M f (u r u) = M (u r u) = 0. (44) r h u = r u, τ h. (45) Lemma 6. If u [L p ()] 3, u [L p ()] 3 and u ν [L p ( )] 3 for p > 2, then r h u is well-defined and bounded. Theorem 2. Let τ h be a regular mesh on Ω. If u [H s (Ω)] 3 and u [H s (Ω)] 3 for s k, then u r h u [L 2 (Ω)] 3 + (u r hu) [L 2 (Ω)] 3 Chs ( u [H s (Ω)] 3 + u [H s (Ω)] 3). (46) 5
6 The error estimates for the interpolant holds for (u r h u) and ( (u r h u) separately: u r h u [L2 (Ω)] 3 Chs ( u [Hs (Ω)] 3 + u [H s (Ω)] 3) (47) u r h u [L 2 (Ω)] 3 Chs u [H s (Ω)] 3 (for s > ) (48) (u r h u) [L2 (Ω)] 3 Chs u [Hs (Ω)] 3. (49) Theorem 3. Let τ h be a regular mesh on Ω. If u [H 2 +δ ()] 3, 0 < δ 2 and u D k = [P k ] 3 P k x, then u r h u [L 2 ()] 3 C(h 2 +δ u [H /2+δ ()] + h u 3 [L 2 ()] 3). (50) The error estimates show that the convergence is u r h u [L 2 (Ω)] 3 = O(hk ), according to the Bramble-Hilbert lemma since [P k ] 3 R k, but [P k ] 3 R k. For triangular elements in two dimensions it has been shown that superconvergence is obtained, so that the order of convergence is similar to nodal based elements, which have O(h k+ ) convergence for degree k. To obtain higher order convergence, a second family of vector elements was introduced by Nedelec: is a tetrahedron P = [P k ] 3 Three types of degrees of freedom associated with edges e of, faces f of and itself. The unit vector along edge e is denoted τ. M e (u) = u τqds, q P k (e), (5) e { M f (u) = u T qda, q D k (f), (52) area(f) f M (u) = u qdv, q D k (). (53) The total set of degrees of freedom is then N = M e (u) M f (u) M (u). The dimension of the second family of elements is 2 (k + )(k + 2)(k + 3). References [] C. Amrouce, C. Bernardi, M. Dauge, and V. Girault. Vector potentials in three-dimensional non-smooth domains. Mathematical Methods in the Applied Sciences, 2: , 998. [2] D. Boffi, P. Fernandes, L. Gastaldi, and I. Perugia. Computational models of electromagnetic resonators: Analysis of edge element approximation. SIAM J. Numer. Anal., 36: , 999. [3] Daniele Boffi, Marco Farina, and Lucia Gastaldi. On the approximation of Maxwell s eigenproblem in general 2D domains. Computers and Structures, 79: , 200. [4] Daniele Boffi and Lucia Gastaldi. Interpolation estimates for edge finite elements and application to band gap computation. Applied Numerical Mathematics, 56: , [5] Martin Costabel and Monique Dauge. Maxwell and lame eigenvalues on polyhedra. Mathematical Methods in the Applied Sciences, 22: , 999. [6] R. Hiptmair. Finite elements in computational electromagnetics. Acta Numerica, : ,
7 [7] Jianming Jin. The Finite Element Method in Electromagnetics. John Wiley & Sons, 2 edition, [8] Daniel R. Lynch and eith D. Paulsen. Origin of vector parasites in numerical Maxwell solutions. IEEE Transactions on Microwave Theory and Techniques, 39, 99. [9] Peter Monk. Finite Element Methods for Maxwell s Equations. Oxford Science Publications, [0] J. C. Nedelec. Mixed finite elements in R 3. Numerische Mathematik, 35(3):35 34, September 980. [] J.C. Nedelec. A new family of mixed finite elements in R 3. Numerische Mathematik, 50:57 8,
Weighted Regularization of Maxwell Equations Computations in Curvilinear Polygons
Weighted Regularization of Maxwell Equations Computations in Curvilinear Polygons Martin Costabel, Monique Dauge, Daniel Martin and Gregory Vial IRMAR, Université de Rennes, Campus de Beaulieu, Rennes,
More informationPAijpam.eu NEW H 1 (Ω) CONFORMING FINITE ELEMENTS ON HEXAHEDRA
International Journal of Pure and Applied Mathematics Volume 109 No. 3 2016, 609-617 ISSN: 1311-8080 (printed version); ISSN: 1314-3395 (on-line version) url: http://www.ijpam.eu doi: 10.12732/ijpam.v109i3.10
More informationWe consider the Maxwell s equation with periodic coefficients
1 Summary We consider the Maxwell s equation with periodic coefficients as it is usually done for the modeling of photonic crystals. Using Bloch/Floquet theory, the problem reduces in a standard way to
More informationFinite Element Methods for Maxwell Equations
CHAPTER 8 Finite Element Methods for Maxwell Equations The Maxwell equations comprise four first-order partial differential equations linking the fundamental electromagnetic quantities, the electric field
More informationHEXAHEDRAL H(DIV) AND H(CURL) FINITE ELEMENTS
HEXAHEDRAL H(DIV) AND H(CURL) FINITE ELEMENTS RICHARD S. FAL, PAOLO GATTO, AND PETER MON Abstract. We study the approximation properties of some finite element subspaces of H(div; Ω) and H(curl ; Ω) defined
More informationNumerical approximation of output functionals for Maxwell equations
Numerical approximation of output functionals for Maxwell equations Ferenc Izsák ELTE, Budapest University of Twente, Enschede 11 September 2004 MAXWELL EQUATIONS Assumption: electric field ( electromagnetic
More informationMultigrid Methods for Maxwell s Equations
Multigrid Methods for Maxwell s Equations Jintao Cui Institute for Mathematics and Its Applications University of Minnesota Outline Nonconforming Finite Element Methods for a Two Dimensional Curl-Curl
More informationStandard Finite Elements and Weighted Regularization
Standard Finite Elements and Weighted Regularization A Rehabilitation Martin COSTABEL & Monique DAUGE Institut de Recherche MAthématique de Rennes http://www.maths.univ-rennes1.fr/~dauge Slides of the
More informationANALYSIS OF A FINITE ELEMENT PML APPROXIMATION FOR THE THREE DIMENSIONAL TIME-HARMONIC MAXWELL PROBLEM
MATHEMATICS OF COMPUTATION Volume 77, Number 261, January 2008, Pages 1 10 S 0025-5718(07)02037-6 Article electronically published on September 18, 2007 ANALYSIS OF A FINITE ELEMENT PML APPROXIMATION FOR
More informationA Multigrid Method for Two Dimensional Maxwell Interface Problems
A Multigrid Method for Two Dimensional Maxwell Interface Problems Susanne C. Brenner Department of Mathematics and Center for Computation & Technology Louisiana State University USA JSA 2013 Outline A
More informationA High-Order Nodal Finite Element Formulation for Microwave Engineering
A High-Order Nodal Finite Element Formulation for Microwave Engineering Rubén Otín CIMNE - International Center For Numerical Methods in Engineering Parque Tecnológico del Mediterráneo, Edificio C4 - Despacho
More informationINTRODUCTION TO FINITE ELEMENT METHODS
INTRODUCTION TO FINITE ELEMENT METHODS LONG CHEN Finite element methods are based on the variational formulation of partial differential equations which only need to compute the gradient of a function.
More informationMixed Finite Element Method for 2D Vector Maxwell s Eigenvalue Problem in Anisotropic Media
Progress In Electromagnetics Research, Vol. 48, 59 7, 24 Mixed Finite Element Method for 2D Vector Maxwell s Eigenvalue Problem in Anisotropic Media Wei Jiang, Na Liu, Yifa Tang 2, and Qing Huo Liu 3,
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Nonconformity and the Consistency Error First Strang Lemma Abstract Error Estimate
More informationAn interpolation operator for H 1 functions on general quadrilateral and hexahedral meshes with hanging nodes
An interpolation operator for H 1 functions on general quadrilateral and hexahedral meshes with hanging nodes Vincent Heuveline Friedhelm Schieweck Abstract We propose a Scott-Zhang type interpolation
More informationChapter 5 A priori error estimates for nonconforming finite element approximations 5.1 Strang s first lemma
Chapter 5 A priori error estimates for nonconforming finite element approximations 51 Strang s first lemma We consider the variational equation (51 a(u, v = l(v, v V H 1 (Ω, and assume that the conditions
More informationLOW ORDER CROUZEIX-RAVIART TYPE NONCONFORMING FINITE ELEMENT METHODS FOR APPROXIMATING MAXWELL S EQUATIONS
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 5, Number 3, Pages 373 385 c 28 Institute for Scientific Computing and Information LOW ORDER CROUZEIX-RAVIART TYPE NONCONFORMING FINITE ELEMENT
More informationQUADRILATERAL H(DIV) FINITE ELEMENTS
QUADRILATERAL H(DIV) FINITE ELEMENTS DOUGLAS N. ARNOLD, DANIELE BOFFI, AND RICHARD S. FALK Abstract. We consider the approximation properties of quadrilateral finite element spaces of vector fields defined
More informationPOLYNOMIAL APPROXIMATIONS OF REGULAR AND SINGULAR VECTOR FIELDS WITH APPLICATIONS TO PROBLEMS OF ELECTROMAGNETICS
POLYNOMIAL APPROXIMATIONS OF REGULAR AND SINGULAR VECTOR FIELDS WITH APPLICATIONS TO PROBLEMS OF ELECTROMAGNETICS Alex Besalov School of Mathematics The University of Manchester Collaborators: Norbert
More informationA Sufficient Condition of Having Independent TE and TM Modes in a Waveguide Filled with Homogenous Anisotropic Lossless Medium
1 A Sufficient Condition of Having Independent TE and TM Modes in a Waveguide Filled with Homogenous Anisotropic Lossless Medium Wei Jiang, Jie Liu, and Qing Huo Liu, Fellow, IEEE arxiv:1603.06439v1 [math.na]
More informationInstitut de Recherche MAthématique de Rennes
LMS Durham Symposium: Computational methods for wave propagation in direct scattering. - July, Durham, UK The hp version of the Weighted Regularization Method for Maxwell Equations Martin COSTABEL & Monique
More informationA NONCONFORMING PENALTY METHOD FOR A TWO DIMENSIONAL CURL-CURL PROBLEM
A NONCONFORMING PENALTY METHOD FOR A TWO DIMENSIONAL CURL-CURL PROBLEM SUSANNE C. BRENNER, FENGYAN LI, AND LI-YENG SUNG Abstract. A nonconforming penalty method for a two-dimensional curl-curl problem
More informationCorner Singularities
Corner Singularities Martin Costabel IRMAR, Université de Rennes 1 Analysis and Numerics of Acoustic and Electromagnetic Problems Linz, 10 15 October 2016 Martin Costabel (Rennes) Corner Singularities
More informationA P4 BUBBLE ENRICHED P3 DIVERGENCE-FREE FINITE ELEMENT ON TRIANGULAR GRIDS
A P4 BUBBLE ENRICHED P3 DIVERGENCE-FREE FINITE ELEMENT ON TRIANGULAR GRIDS SHANGYOU ZHANG DEDICATED TO PROFESSOR PETER MONK ON THE OCCASION OF HIS 6TH BIRTHDAY Abstract. On triangular grids, the continuous
More informationLocal discontinuous Galerkin methods for elliptic problems
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2002; 18:69 75 [Version: 2000/03/22 v1.0] Local discontinuous Galerkin methods for elliptic problems P. Castillo 1 B. Cockburn
More informationTHE APPROXIMATION OF THE MAXWELL EIGENVALUE PROBLEM USING A LEAST-SQUARES METHOD
MATHEMATICS OF COMPUTATION Volume 00, Number 0, Pages 000 000 S 0025-5718(XX)0000-0 THE APPROXIMATION OF THE MAXWELL EIGENVALUE PROBLEM USING A LEAST-SQUARES METHOD JAMES H. BRAMBLE, TZANIO V. KOLEV, AND
More informationMatrix construction: Singular integral contributions
Matrix construction: Singular integral contributions Seminar Boundary Element Methods for Wave Scattering Sophie Haug ETH Zurich November 2010 Outline 1 General concepts in singular integral computation
More informationENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED FINITE ELEMENT METHODS
ENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED FINITE ELEMENT METHODS CARLO LOVADINA AND ROLF STENBERG Abstract The paper deals with the a-posteriori error analysis of mixed finite element methods
More informationRichard S. Falk 1, Paolo Gatto 2 and Peter Monk 3
ESAIM: M2AN 45 (2011) 115 143 DOI: 10.1051/m2an/2010034 ESAIM: Mathematical Modelling and Numerical Analysis www.esaim-m2an.org HEXAHEDRAL H(DIV) AND H(CURL) FINITE ELEMENTS Richard S. Falk 1, Paolo Gatto
More informationProblems of Corner Singularities
Stuttgart 98 Problems of Corner Singularities Monique DAUGE Institut de Recherche MAthématique de Rennes Problems of Corner Singularities 1 Vertex and edge singularities For a polyhedral domain Ω and a
More informationA Mixed Method for Axisymmetric Div-Curl Systems
Portland State University PDXScholar Mathematics and Statistics Faculty Publications and Presentations Fariborz Maseeh Department of Mathematics and Statistics 2008 A Mixed Method for Axisymmetric Div-Curl
More informationDaniele Boffi 1 and Lucia Gastaldi 2
Mathematical Modelling and Numerical Analysis ESAIM: M2AN Modélisation Mathématique et Analyse Numérique M2AN, Vol. 36, N o 2, 2002, pp. 293 305 DOI: 10.1051/m2an:2002013 EDGE FINITE ELEMENTS FOR THE APPROXIMATION
More informationplasmas Lise-Marie Imbert-Gérard, Bruno Després. July 27th 2011 Laboratoire J.-L. LIONS, Université Pierre et Marie Curie, Paris.
Lise-Marie Imbert-Gérard, Bruno Després. Laboratoire J.-L. LIONS, Université Pierre et Marie Curie, Paris. July 27th 2011 1 Physical and mathematical motivations 2 approximation of the solutions 3 4 Plan
More informationA priori error analysis of the BEM with graded meshes for the electric eld integral equation on polyhedral surfaces
A priori error analysis of the BEM with graded meshes for the electric eld integral equation on polyhedral surfaces A. Bespalov S. Nicaise Abstract The Galerkin boundary element discretisations of the
More informationInterpolation in h-version finite element spaces
Interpolation in h-version finite element spaces Thomas Apel Institut für Mathematik und Bauinformatik Fakultät für Bauingenieur- und Vermessungswesen Universität der Bundeswehr München Chemnitzer Seminar
More informationFinite element approximation on quadrilateral meshes
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2001; 17:805 812 (DOI: 10.1002/cnm.450) Finite element approximation on quadrilateral meshes Douglas N. Arnold 1;, Daniele
More informationarxiv: v1 [math.na] 29 Feb 2016
EFFECTIVE IMPLEMENTATION OF THE WEAK GALERKIN FINITE ELEMENT METHODS FOR THE BIHARMONIC EQUATION LIN MU, JUNPING WANG, AND XIU YE Abstract. arxiv:1602.08817v1 [math.na] 29 Feb 2016 The weak Galerkin (WG)
More informationAn introduction to Isogeometric Analysis
An introduction to Isogeometric Analysis A. Buffa 1, G. Sangalli 1,2, R. Vázquez 1 1 Istituto di Matematica Applicata e Tecnologie Informatiche - Pavia Consiglio Nazionale delle Ricerche 2 Dipartimento
More informationICES REPORT Direct Serendipity Finite Elements on Convex Quadrilaterals
ICES REPORT 17-8 October 017 Direct Serendipity Finite Elements on Convex Quadrilaterals by Todd Arbogast and Zhen Tao The Institute for Computational Engineering and Sciences The University of Texas at
More informationSemi-discrete finite element approximation applied to Maxwell s equations in nonlinear media
Semi-discrete finite element approximation applied to Maxwell s equations in nonlinear media Lutz Angermann arxiv:9.365v math.na Jan 9 January 9, 9 In this paper the semi-discrete finite element approximation
More informationOverlapping Schwarz Preconditioners for Spectral. Problem in H(curl)
Overlapping Schwarz Preconditioners for Spectral Nédélec Elements for a Model Problem in H(curl) Technical Report TR2002-83 November 22, 2002 Department of Computer Science Courant Institute of Mathematical
More informationAN OPTIMAL DOMAIN DECOMPOSITION PRECONDITIONER FOR LOW-FREQUENCY TIME-HARMONIC MAXWELL EQUATIONS
MATHEMATICS OF COMPUTATION Volume 68, Number 226, April 1999, Pages 607 631 S 0025-5718(99)01013-3 AN OPTIMAL DOMAIN DECOMPOSITION PRECONDITIONER FOR LOW-FREQUENCY TIME-HARMONIC MAXWELL EQUATIONS ANA ALONSO
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Variational Problems of the Dirichlet BVP of the Poisson Equation 1 For the homogeneous
More informationFinite Elements for Elastic Shell Models in
Elastic s in Advisor: Matthias Heinkenschloss Computational and Applied Mathematics Rice University 13 April 2007 Outline Elasticity in Differential Geometry of Shell Geometry and Equations The Plate Model
More informationEXISTENCE AND REGULARITY OF SOLUTIONS FOR STOKES SYSTEMS WITH NON-SMOOTH BOUNDARY DATA IN A POLYHEDRON
Electronic Journal of Differential Equations, Vol. 2017 (2017), No. 147, pp. 1 10. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu EXISTENCE AND REGULARITY OF SOLUTIONS FOR
More informationFrom the Boundary Element Domain Decomposition Methods to Local Trefftz Finite Element Methods on Polyhedral Meshes
From the Boundary Element Domain Decomposition Methods to Local Trefftz Finite Element Methods on Polyhedral Meshes Dylan Copeland 1, Ulrich Langer 2, and David Pusch 3 1 Institute of Computational Mathematics,
More informationAdaptive approximation of eigenproblems: multiple eigenvalues and clusters
Adaptive approximation of eigenproblems: multiple eigenvalues and clusters Francesca Gardini Dipartimento di Matematica F. Casorati, Università di Pavia http://www-dimat.unipv.it/gardini Banff, July 1-6,
More informationA MIXED METHOD FOR AXISYMMETRIC DIV-CURL SYSTEMS
A MIXED METHOD FOR AXISYMMETRIC DIV-CURL SYSTEMS DYLAN M. COPELAND, JAYADEEP GOPALAKRISHNAN, AND JOSEPH E. PASCIAK Abstract. We present a mixed method for a three-dimensional axisymmetric div-curl system
More informationMACRO STOKES ELEMENTS ON QUADRILATERALS
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 15, Number 4-5, Pages 729 745 c 2018 Institute for Scientific Computing and Information MACRO STOES ELEMENTS ON QUADRILATERALS MICHAEL NEILAN
More informationSUPERCONVERGENCE PROPERTIES FOR OPTIMAL CONTROL PROBLEMS DISCRETIZED BY PIECEWISE LINEAR AND DISCONTINUOUS FUNCTIONS
SUPERCONVERGENCE PROPERTIES FOR OPTIMAL CONTROL PROBLEMS DISCRETIZED BY PIECEWISE LINEAR AND DISCONTINUOUS FUNCTIONS A. RÖSCH AND R. SIMON Abstract. An optimal control problem for an elliptic equation
More informationA local-structure-preserving local discontinuous Galerkin method for the Laplace equation
A local-structure-preserving local discontinuous Galerkin method for the Laplace equation Fengyan Li 1 and Chi-Wang Shu 2 Abstract In this paper, we present a local-structure-preserving local discontinuous
More informationScientific Computing WS 2018/2019. Lecture 15. Jürgen Fuhrmann Lecture 15 Slide 1
Scientific Computing WS 2018/2019 Lecture 15 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 15 Slide 1 Lecture 15 Slide 2 Problems with strong formulation Writing the PDE with divergence and gradient
More informationPREPRINT 2010:23. A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO
PREPRINT 2010:23 A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO Department of Mathematical Sciences Division of Mathematics CHALMERS UNIVERSITY OF TECHNOLOGY UNIVERSITY OF GOTHENBURG
More informationIntroduction to finite element exterior calculus
Introduction to finite element exterior calculus Ragnar Winther CMA, University of Oslo Norway Why finite element exterior calculus? Recall the de Rham complex on the form: R H 1 (Ω) grad H(curl, Ω) curl
More informationChapter 3 Conforming Finite Element Methods 3.1 Foundations Ritz-Galerkin Method
Chapter 3 Conforming Finite Element Methods 3.1 Foundations 3.1.1 Ritz-Galerkin Method Let V be a Hilbert space, a(, ) : V V lr a bounded, V-elliptic bilinear form and l : V lr a bounded linear functional.
More informationA posteriori error estimates for a Maxwell type problem
Russ. J. Numer. Anal. Math. Modelling, Vol. 24, No. 5, pp. 395 408 (2009) DOI 0.55/ RJNAMM.2009.025 c de Gruyter 2009 A posteriori error estimates for a Maxwell type problem I. ANJAM, O. MALI, A. MUZALEVSKY,
More informationMULTIGRID PRECONDITIONING IN H(div) ON NON-CONVEX POLYGONS* Dedicated to Professor Jim Douglas, Jr. on the occasion of his seventieth birthday.
MULTIGRID PRECONDITIONING IN H(div) ON NON-CONVEX POLYGONS* DOUGLAS N ARNOLD, RICHARD S FALK, and RAGNAR WINTHER Dedicated to Professor Jim Douglas, Jr on the occasion of his seventieth birthday Abstract
More informationConstruction of H(div)-Conforming Mixed Finite Elements on Cuboidal Hexahedra
Numerische Mathematik manuscript No. (will be inserted by the editor) Construction of H(div)-Conforming Mixed Finite Elements on Cuboidal Hexahedra Todd Arbogast Zhen Tao Received: July 9, 27 / Accepted:
More informationA posteriori error estimates for non conforming approximation of eigenvalue problems
A posteriori error estimates for non conforming approximation of eigenvalue problems E. Dari a, R. G. Durán b and C. Padra c, a Centro Atómico Bariloche, Comisión Nacional de Energía Atómica and CONICE,
More informationEXACT DE RHAM SEQUENCES OF SPACES DEFINED ON MACRO-ELEMENTS IN TWO AND THREE SPATIAL DIMENSIONS
EXACT DE RHAM SEQUENCES OF SPACES DEFINED ON MACRO-ELEMENTS IN TWO AND THREE SPATIAL DIMENSIONS JOSEPH E. PASCIAK AND PANAYOT S. VASSILEVSKI Abstract. This paper proposes new finite element spaces that
More informationweak Galerkin, finite element methods, interior estimates, second-order elliptic
INERIOR ENERGY ERROR ESIMAES FOR HE WEAK GALERKIN FINIE ELEMEN MEHOD HENGGUANG LI, LIN MU, AND XIU YE Abstract Consider the Poisson equation in a polytopal domain Ω R d (d = 2, 3) as the model problem
More informationarxiv: v1 [math.na] 27 Jan 2016
Virtual Element Method for fourth order problems: L 2 estimates Claudia Chinosi a, L. Donatella Marini b arxiv:1601.07484v1 [math.na] 27 Jan 2016 a Dipartimento di Scienze e Innovazione Tecnologica, Università
More informationA posteriori error estimation for elliptic problems
A posteriori error estimation for elliptic problems Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in
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 informationHölder regularity for Maxwell's equations under minimal assumptions on the coefficients
Hölder regularity for Maxwell's equations under minimal assumptions on the coefficients G. S. Alberti Research Report No. 2016-22 April 2016 Seminar für Angewte Mathematik Eidgenössische Technische Hochschule
More informationBasic Principles of Weak Galerkin Finite Element Methods for PDEs
Basic Principles of Weak Galerkin Finite Element Methods for PDEs Junping Wang Computational Mathematics Division of Mathematical Sciences National Science Foundation Arlington, VA 22230 Polytopal Element
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University The Residual and Error of Finite Element Solutions Mixed BVP of Poisson Equation
More informationMultigrid Methods for Elliptic Obstacle Problems on 2D Bisection Grids
Multigrid Methods for Elliptic Obstacle Problems on 2D Bisection Grids Long Chen 1, Ricardo H. Nochetto 2, and Chen-Song Zhang 3 1 Department of Mathematics, University of California at Irvine. chenlong@math.uci.edu
More informationA LOWEST-ORDER COMPOSITE FINITE ELEMENT EXACT SEQUENCE ON PYRAMIDS arxiv: v1 [math.na] 28 Apr 2017
A LOWEST-ORDER COMPOSITE FINITE ELEMENT EXACT SEQUENCE ON PYRAMIDS arxiv:1705.00064v1 [math.na] 8 Apr 017 MARK AINSWORTH AND GUOSHENG FU Abstract. Composite basis functions for pyramidal elements on the
More informationA Remark on the Regularity of Solutions of Maxwell s Equations on Lipschitz Domains
A Remark on the Regularity of Solutions of Maxwell s Equations on Lipschitz Domains Martin Costabel Abstract Let u be a vector field on a bounded Lipschitz domain in R 3, and let u together with its divergence
More informationDivergence-free or curl-free finite elements for solving the curl-div system
Divergence-free or curl-free finite elements for solving the curl-div system Alberto Valli Dipartimento di Matematica, Università di Trento, Italy Joint papers with: Ana Alonso Rodríguez Dipartimento di
More informationNumerical analysis of problems in electromagnetism
Numerical analysis of problems in electromagnetism ALBERTO VALLI Department of Mathematics, University of Trento Numerical analysis of problems in electromagnetism p.1/195 Finite elements A finite element
More informationA posteriori error estimates for Maxwell Equations
www.oeaw.ac.at A posteriori error estimates for Maxwell Equations J. Schöberl RICAM-Report 2005-10 www.ricam.oeaw.ac.at A POSTERIORI ERROR ESTIMATES FOR MAXWELL EQUATIONS JOACHIM SCHÖBERL Abstract. Maxwell
More informationHigh order, finite volume method, flux conservation, finite element method
FLUX-CONSERVING FINITE ELEMENT METHODS SHANGYOU ZHANG, ZHIMIN ZHANG, AND QINGSONG ZOU Abstract. We analyze the flux conservation property of the finite element method. It is shown that the finite element
More informationANDREA TOSELLI. Abstract. Two-level overlapping Schwarz methods are considered for nite element problems
OVERLAPPING SCHWARZ METHODS FOR MAXWELL'S EQUATIONS IN THREE DIMENSIONS ANDREA TOSELLI Abstract. Two-level overlapping Schwarz methods are considered for nite element problems of 3D Maxwell's equations.
More informationFrom the Boundary Element DDM to local Trefftz Finite Element Methods on Polyhedral Meshes
www.oeaw.ac.at From the Boundary Element DDM to local Trefftz Finite Element Methods on Polyhedral Meshes D. Copeland, U. Langer, D. Pusch RICAM-Report 2008-10 www.ricam.oeaw.ac.at From the Boundary Element
More informationHigh Order Differential Form-Based Elements for the Computation of Electromagnetic Field
1472 IEEE TRANSACTIONS ON MAGNETICS, VOL 36, NO 4, JULY 2000 High Order Differential Form-Based Elements for the Computation of Electromagnetic Field Z Ren, Senior Member, IEEE, and N Ida, Senior Member,
More informationAMS subject classifications. Primary, 65N15, 65N30, 76D07; Secondary, 35B45, 35J50
A SIMPLE FINITE ELEMENT METHOD FOR THE STOKES EQUATIONS LIN MU AND XIU YE Abstract. The goal of this paper is to introduce a simple finite element method to solve the Stokes equations. This method is in
More informationA Least-Squares Finite Element Approximation for the Compressible Stokes Equations
A Least-Squares Finite Element Approximation for the Compressible Stokes Equations Zhiqiang Cai, 1 Xiu Ye 1 Department of Mathematics, Purdue University, 1395 Mathematical Science Building, West Lafayette,
More informationRegularized Maxwell Equations and Nodal Finite Elements for Electromagnetic Field Computations
Regularized Maxwell Equations and Nodal Finite Elements for Electromagnetic Field Computations Ruben Otin CIMNE - International Center For Numerical Methods in Engineering Parque Tecnológico del Mediterráneo,
More informationLecture Notes: African Institute of Mathematics Senegal, January Topic Title: A short introduction to numerical methods for elliptic PDEs
Lecture Notes: African Institute of Mathematics Senegal, January 26 opic itle: A short introduction to numerical methods for elliptic PDEs Authors and Lecturers: Gerard Awanou (University of Illinois-Chicago)
More informationAPPROXIMATION OF THE EIGENVALUE PROBLEM FOR THE TIME HARMONIC MAXWELL SYSTEM BY CONTINUOUS LAGRANGE FINITE ELEMENTS
APPROXIMATION OF THE EIGENVALUE PROBLEM FOR THE TIME HARMONIC MAXWELL SYSTEM BY CONTINUOUS LAGRANGE FINITE ELEMENTS ANDREA BONITO, AND JEAN-LUC GUERMOND,3 Abstract. We propose and analyze an approximation
More informationICES REPORT Analysis of the DPG Method for the Poisson Equation
ICES REPORT 10-37 September 2010 Analysis of the DPG Method for the Poisson Equation by L. Demkowicz and J. Gopalakrishnan The Institute for Computational Engineering and Sciences The University of Texas
More informationFinite Elements. Colin Cotter. February 22, Colin Cotter FEM
Finite Elements February 22, 2019 In the previous sections, we introduced the concept of finite element spaces, which contain certain functions defined on a domain. Finite element spaces are examples of
More informationHigher order hierarchical H(curl) Legendre basis functions applied in the finite element method: emphasis on microwave circuits
Higher order hierarchical H(curl) Legendre basis functions applied in the finite element method: emphasis on microwave circuits Johannesson, Peter 20 Link to publication Citation for published version
More informationThe Pole Condition: A Padé Approximation of the Dirichlet to Neumann Operator
The Pole Condition: A Padé Approximation of the Dirichlet to Neumann Operator Martin J. Gander and Achim Schädle Mathematics Section, University of Geneva, CH-, Geneva, Switzerland, Martin.gander@unige.ch
More informationTHE PATCH RECOVERY FOR FINITE ELEMENT APPROXIMATION OF ELASTICITY PROBLEMS UNDER QUADRILATERAL MESHES. Zhong-Ci Shi and Xuejun Xu.
DISCRETE AND CONTINUOUS Website: http://aimsciences.org DYNAMICAL SYSTEMS SERIES B Volume 9, Number, January 2008 pp. 63 82 THE PATCH RECOVERY FOR FINITE ELEMENT APPROXIMATION OF ELASTICITY PROBLEMS UNDER
More informationSolving the curl-div system using divergence-free or curl-free finite elements
Solving the curl-div system using divergence-free or curl-free finite elements Alberto Valli Dipartimento di Matematica, Università di Trento, Italy or: Why I say to my students that divergence-free finite
More informationb i (x) u + c(x)u = f in Ω,
SIAM J. NUMER. ANAL. Vol. 39, No. 6, pp. 1938 1953 c 2002 Society for Industrial and Applied Mathematics SUBOPTIMAL AND OPTIMAL CONVERGENCE IN MIXED FINITE ELEMENT METHODS ALAN DEMLOW Abstract. An elliptic
More informationDefect-based a-posteriori error estimation for implicit ODEs and DAEs
1 / 24 Defect-based a-posteriori error estimation for implicit ODEs and DAEs W. Auzinger Institute for Analysis and Scientific Computing Vienna University of Technology Workshop on Innovative Integrators
More informationarxiv: v2 [math.na] 23 Apr 2016
Improved ZZ A Posteriori Error Estimators for Diffusion Problems: Conforming Linear Elements arxiv:508.009v2 [math.na] 23 Apr 206 Zhiqiang Cai Cuiyu He Shun Zhang May 2, 208 Abstract. In [8], we introduced
More informationThe Nullspace free eigenvalue problem and the inexact Shift and invert Lanczos method. V. Simoncini. Dipartimento di Matematica, Università di Bologna
The Nullspace free eigenvalue problem and the inexact Shift and invert Lanczos method V. Simoncini Dipartimento di Matematica, Università di Bologna and CIRSA, Ravenna, Italy valeria@dm.unibo.it 1 The
More informationLecture Note III: Least-Squares Method
Lecture Note III: Least-Squares Method Zhiqiang Cai October 4, 004 In this chapter, we shall present least-squares methods for second-order scalar partial differential equations, elastic equations of solids,
More informationWHITNEY/NÉDÉLEC ELEMENTS METHOD APPROACH APPLIED IN THE MAXWELL S EQUATIONS
WHITNEY/NÉDÉLEC ELEMENTS METHOD APPROACH APPLIED IN THE MAXWELL S EQUATIONS Jean Eduardo Sebold jeansebold@ufpr.br/jean.sebold@ifc-araquari.edu.br 1,2 Saulo Pomponet Oliveira saulopo@ufpr.br 2 Luiz Alkimin
More informationPOLYNOMIAL PRESERVING GRADIENT RECOVERY AND A POSTERIORI ESTIMATE FOR BILINEAR ELEMENT ON IRREGULAR QUADRILATERALS
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume, Number, Pages 24 c 2004 Institute for Scientific Computing and Information POLYNOMIAL PRESERVING GRADIENT RECOVERY AND A POSTERIORI ESTIMATE
More informationOn Friedrichs inequality, Helmholtz decomposition, vector potentials, and the div-curl lemma. Ben Schweizer 1
On Friedrichs inequality, Helmholtz decomposition, vector potentials, and the div-curl lemma Ben Schweizer 1 January 16, 2017 Abstract: We study connections between four different types of results that
More informationOverlapping Schwarz preconditioners for Fekete spectral elements
Overlapping Schwarz preconditioners for Fekete spectral elements R. Pasquetti 1, L. F. Pavarino 2, F. Rapetti 1, and E. Zampieri 2 1 Laboratoire J.-A. Dieudonné, CNRS & Université de Nice et Sophia-Antipolis,
More informationLocal flux mimetic finite difference methods
Local flux mimetic finite difference methods Konstantin Lipnikov Mikhail Shashkov Ivan Yotov November 4, 2005 Abstract We develop a local flux mimetic finite difference method for second order elliptic
More informationSuperconvergence and H(div) Projection for Discontinuous Galerkin Methods
INTRNATIONAL JOURNAL FOR NUMRICAL MTHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2000; 00:1 6 [Version: 2000/07/27 v1.0] Superconvergence and H(div) Projection for Discontinuous Galerkin Methods Peter Bastian
More informationKey words. Incompressible magnetohydrodynamics, mixed finite element methods, discontinuous Galerkin methods
A MIXED DG METHOD FOR LINEARIZED INCOMPRESSIBLE MAGNETOHYDRODYNAMICS PAUL HOUSTON, DOMINIK SCHÖTZAU, AND XIAOXI WEI Journal of Scientific Computing, vol. 40, pp. 8 34, 009 Abstract. We introduce and analyze
More information