Weighted Regularization of Maxwell Equations Computations in Curvilinear Polygons
|
|
- Julie Clark
- 5 years ago
- Views:
Transcription
1 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, France URL: Maxwell eigenfrequency problem The problem consists in finding the non-zero frequencies ω> such that there exists an electromagnetic field (E, H) solving the equation curl εe iω µh = and curl µh + iω εe = in Ω, () in the cavity Ω. We assume perfect conductor boundary conditions, that is E n = and H n =on Ω, where n denotes the outer unit normal on Ω. Concentrating on the problems posed by the singularities of the boundary of Ω, we assume that the dielectric material filling Ω is homogeneous and isotropic, that is, after a possible change of unknowns ε = µ =. We want to solve () by a Galerkin procedure. The first step is to propose a variational formulation for (). Eliminating H from the equations, we find the formulation E H (curl; Ω), E H (curl; Ω), curl E curl E dx = ω E E dx, () Ω Ω where H (curl; Ω) is the space of L (Ω) fields E with curl E L (Ω) and tangential component E n =on Ω. Testing against gradients of functions ϕ H (Ω) we find that if ω, then div E =. Conversely, any such gradient is a solution of () with ω =. There are two classes of strategies to overcome the difficulties generated by the presence of this big kernel: (i) Discretize () directly by special families of finite elements possessing a discrete derham complex property: these are the edge elements first introduced by Nedelec [, ]. (ii) Regularize () by a term div, div and use more standard elements. Strategy (i) is widely spread but, still, it is worthwhile to develop an alternative, which will have different advantages.
2 Martin Costabel, Monique Dauge, Daniel Martin and Gregory Vial At first glance, strategy (ii) is a very bad idea if Ω has non-convex corners or edges. The reason for this is now well known: The variational space for the continuous problem is X N = { E H (curl; Ω), div E L (Ω) } () whereas any finite element field which is conforming for the bilinear form curl, curl L + div, div L is continuous across the faces of the elements, therefore belong to the subspace H N X N of H (Ω) fields. But H N is closed in X N for the norm curl L + div L and finally H N X N as soon as Ω has reentrant corners in D and reentrant edges in D. Thus we are in the somewhat unusual situation where the convergence of the Galerkin method does not mean that the limit is the solution of the continuous problem. Weighted regularization In order to overcome this difficulty due to the lack of density, a weighted regularization was proposed in [], that is replacing the regularization term div, div L by the scalar product div, div Y in another functional space Y containing L (Ω) and contained in H (Ω). Such a space can be conveniently realized as a weighted L space, where the weight is a power γ, γ, of a distance function d : x d(x) to the reentrant corners or edges of Ω. The variational space is, instead of () X N [γ] = { E H (curl; Ω), d(x) γ div E L (Ω) } () and the corresponding variational formulation is (with α := γ and s > chosen ad libitum) E X N [γ], E X N [γ], (curl E curl E + sd(x) α div E div E ) dx = ω Ω Ω E E dx. () It is clear that the eigenpairs (E,ω) with ω of () are eigenpairs of (). Any grad ϕ with ϕ H (Ω) is an eigenvector of () with ω =. For problem () also the spurious eigenvectors are gradients grad ϕ with ϕ H (Ω) and d γ ϕ L (Ω): ϕ is such that d γ ϕ is an eigenvector of the self-adjoint operator A[γ] :=d γ d γ with Dirichlet conditions. One can prove Theorem. The eigenpairs (E,ω) of problem () can be split into two types (i) div E =and ω : (E,ω) is then a Maxwell eigenpair. (ii) curl E =and ω = sν with ν and eigenvalue of the operator A[γ]. It is important to know that the spectrum of A[γ] is separated from :
3 Weighted Regularization of Maxwell Equations Theorem. The eigenvalues of A[γ] decrease when γ increases. For any γ, A[γ] is positive definite. For any γ<, A[γ] has a compact resolvent. From Theorems and we know that the continuous formulation () can be used to determine the Maxwell eigenvalues for any α. Now, we have to investigate whether there exists α so that smooth fields are dense in X N [γ]. The following statement is deduced from [] Theorem. Let Ω be non-convex polygon or a non-convex polyhedron. Then there exists γ Ω < such that for any γ, γ Ω <γ<, the space H N of H vector fields satisfying the electric boundary condition is dense in X N [γ]. If Ω is a polygon, let ω Ω be its largest corner. Then γ Ω = π/ω Ω. There holds a similar formula for polyhedra. Moreover, it is proved in [] that for the source problem in a polygon, if nodal elements of degree p are used, we have a convergence rate in energy norm for the h-version of finite elements in h τ ε for all ε>, where τ = min { γ + π ω Ω, + with ω conv the largest convex angle of Ω. π ω conv, π }, ω Ω Experiments in an L-shape domain The condition p is not necessary for practical convergence of the method (it was simply comfortable for the proof). On the other hand, owing to the classical structure of the singularities (though unbounded) we can guess that a geometrical refinement of the mesh near singular points of the boundary will improve the performance of the method layers layers. Fig.. Meshes on the L-shape domain Figure presents our first domain of computation Ω with the first two levels of refinement of the mesh. Here we use a ratio.
4 Martin Costabel, Monique Dauge, Daniel Martin and Gregory Vial We have computed the first Maxwell eigenvalue by a Galerkin discretization of formulation () with α = in eight nested meshes where the number of layers increases from to and with nodal elements based on tensor Q p polynomials, p ranging from to. Note that, since the edges of the boundary Ω are straight, the essential boundary condition E n = is implemented exactly. The computations are done with the Finite Element Library Mélina []. - - layers layers layers layers layers layers 9 9 layers 9 - layers Fig.. Relative errors for st eigenvalue In Figure relative errors for the first Maxwell eigenvalue are represented: The comparison value ω =. is obtained by solving the Neumann problem on the same domain with the finest mesh and p = the D Maxwell computations work as a benchmark problem since the Maxwell eigenpairs are the (curl ϕ, ν) with ϕ eigenvector of the Laplace Neumann problem with eigenvalue ν. The abscissae represent the number of degrees of freedom, and the ordinates the relative error. Each line connects the relative errors for a fixed number of layers and as p goes from to. The numbers above the lines indicate the values of the degree p. In Figure we represent the convergence rates of the previous relative errors with respect to the number N of degrees of freedom (N is proportional to p on each mesh) and as p increases. The hp effect is clearly visible.
5 Weighted Regularization of Maxwell Equations layers layers layers layers layers layers 9 layers layers - Fig.. Convergence rates of relative errors Experiments in a curved L-domain The domain Ω with curved boundaries (arcs of circle A, A and A of radii, and respectively) for which we did the computations is represented opposite. The discretization of its boundary depends now on the mesh we choose and on the degree q of the geometrical approximation Fig.. The curved L-domain The meshes are similar to those constructed in the straight L-domain of Section and refined by adding successive layers of elements around the reentrant corner. Let m be the number of layers of the mesh. Whatever the value of m we have elements abutting at the exterior arc A, and at the interior arc A. Along A there are m elements. We use a geometrical approximation of degree q ranging from to, independently of the functional
6 Martin Costabel, Monique Dauge, Daniel Martin and Gregory Vial degree p. The values of q and m determine different approximations Ω q,m of Ω, which consist in defining each element by a grid of (q +) of Gauss- Lobatto points. For example, when q =, the discretized domains Ω,m are polygonal approximations of Ω, cf Figure Fig.. Meshes on the approximate curved L-domains Ω q,m, q =,m =, If we represent the relative errors for the first eigenvalue with the eight best geometrical approximations Ω,m with m =,...,, as p increases, we obtain a graph very similar to Figure. We prefer to investigate the influence of the geometrical degree q on the quality of the approximation. We represent in Figures and the relative errors obtained by the computation of the first and second Maxwell eigenvalues with the variational formulation () on the meshes with layers and geometrical degree q. The comparison is done with the Neumann eigenvalues computed on the same mesh with q =. We have to note that our approximation is now non conforming for two reasons:. The domains of computation Ω q,m do not coincide exactly with Ω: the boundary is approximated by piecewise curves of degree q. Only at the grid points a i the boundary Ω q,m of the approximate domain coincides with Ω.. The essential boundary condition E n = is implemented only at the functional (Gauss-Lobatto) nodes b j. If and only if q =, these pointwise conditions imply that E n =on Ω q,m. We can see the relative influence of the position of the geometrical nodes a i and the functional nodes b j by modifying their position. If the a i are defined from an equidistant grid and b j unchanged, the results are very similar. The geometrical degree q is the degree of the local maps which transform the reference element ˆK intothe elements K. The functional degree p is the degree of the interpolating polynomials in ˆK for the test and trial functions. The case p = q corresponds to isoparametric elements.
7 Weighted Regularization of Maxwell Equations Geo. deg. Geo. deg. Geo. deg. - Geo. deg. Geo. deg. Geo. deg Fig.. Relative errors for the st eigenvalue Geo. deg. Geo. deg. Geo. deg. Geo. deg. Geo. deg. Geo. deg Fig.. Relative errors for the nd eigenvalue
8 Martin Costabel, Monique Dauge, Daniel Martin and Gregory Vial If, on the other hand, the a i are Gauss-Lobatto points and the b j are equidistant, we obtain for the second eigenvalue the numerical results contained in Figure Geo. deg. Geo. deg. Geo. deg. 9 - Geo. deg. Geo. deg. - Geo. deg. Fig.. Relative errors for the nd eigenvalue with equidistant functional nodes We see that the choice of equidistant nodes for the interpolation of the boundary condition E n = is a very serious variational crime that prevents further reduction of the error for degrees p. Note that this phenomenon is only observed for vector boundary conditions and curved boundaries. It does not appear in scalar problems nor on domains with piecewise straight boundaries. References. M. Costabel, M. Dauge. Weighted regularization of maxwell equations in polyhedral domains. Preprint -, Université de Rennes (). Toappear in Numer. Math.. D. Martin. Mélina. On line documentation: J.-C. Nédélec. Mixed finite elements in R. Numer. Math. (9).. J.-C. Nédélec. A new family of mixed finite elements in R. Numer. Math. () (9).
Standard 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 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 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 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 informationNedelec elements for computational electromagnetics
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,
More informationTransformation of corner singularities in presence of small or large parameters
Transformation of corner singularities in presence of small or large parameters Monique Dauge IRMAR, Université de Rennes 1, FRANCE Analysis and Numerics of Acoustic and Electromagnetic Problems October
More informationInequalities of Babuška-Aziz and Friedrichs-Velte for differential forms
Inequalities of Babuška-Aziz and Friedrichs-Velte for differential forms Martin Costabel Abstract. For sufficiently smooth bounded plane domains, the equivalence between the inequalities of Babuška Aziz
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 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 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 information1. Fast Solvers and Schwarz Preconditioners for Spectral Nédélec Elements for a Model Problem in H(curl)
DDM Preprint Editors: editor1, editor2, editor3, editor4 c DDM.org 1. Fast Solvers and Schwarz Preconditioners for Spectral Nédélec Elements for a Model Problem in H(curl) Bernhard Hientzsch 1 1. Introduction.
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 informationDomain Decomposition Preconditioners for Spectral Nédélec Elements in Two and Three Dimensions
Domain Decomposition Preconditioners for Spectral Nédélec Elements in Two and Three Dimensions Bernhard Hientzsch Courant Institute of Mathematical Sciences, New York University, 51 Mercer Street, New
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 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 informationSCHRIFTENREIHE DER FAKULTÄT FÜR MATHEMATIK. On Maxwell s and Poincaré s Constants. Dirk Pauly SM-UDE
SCHIFTENEIHE DE FAKULTÄT FÜ MATHEMATIK On Maxwell s and Poincaré s Constants by Dirk Pauly SM-UDE-772 2013 On Maxwell s and Poincaré s Constants Dirk Pauly November 11, 2013 Dedicated to Sergey Igorevich
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 informationInvertibility of the biharmonic single layer potential operator
Invertibility of the biharmonic single layer potential operator Martin COSTABEL & Monique DAUGE Abstract. The 2 2 system of integral equations corresponding to the biharmonic single layer potential in
More informationNumerische Mathematik
Numer. Math. (2012) 122:61 99 DOI 10.1007/s00211-012-0456-x Numerische Mathematik C 0 elements for generalized indefinite Maxwell equations Huoyuan Duan Ping Lin Roger C. E. Tan Received: 31 July 2010
More informationOn the Maxwell Constants in 3D
On the Maxwell Constants in 3D Dirk Pauly Fakultät für Mathematik Universität Duisburg-Essen, Campus Essen, Germany JSA 3 & MC 65 Journées Singulières Augmentés 03: Conférence en l honneur de Martin Costabel
More informationChapter 1 Mathematical Foundations
Computational Electromagnetics; Chapter 1 1 Chapter 1 Mathematical Foundations 1.1 Maxwell s Equations Electromagnetic phenomena can be described by the electric field E, the electric induction D, the
More informationSuperconductivity in domains with corners
Superconductivity in domains with corners Virginie BONNAILLIE-NOËL IRMAR, Université Rennes 1 and ENS Cachan Bretagne Each China Normal University, Shanghai 2007.5.15 Outline 1. Introduction 2. Linear
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 informationMixed Finite Elements Method
Mixed Finite Elements Method A. Ratnani 34, E. Sonnendrücker 34 3 Max-Planck Institut für Plasmaphysik, Garching, Germany 4 Technische Universität München, Garching, Germany Contents Introduction 2. Notations.....................................
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 informationBasics and some applications of the mortar element method
GAMM-Mitt. 28, No. 2, 97 123 (2005) Basics and some applications of the mortar element method Christine Bernardi 1, Yvon Maday 1, and Francesca Rapetti 2 1 Laboratoire Jacques-Louis Lions, C.N.R.S. & université
More informationScattering of electromagnetic waves by thin high contrast dielectrics II: asymptotics of the electric field and a method for inversion.
Scattering of electromagnetic waves by thin high contrast dielectrics II: asymptotics of the electric field and a method for inversion. David M. Ambrose Jay Gopalakrishnan Shari Moskow Scott Rome June
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 informationA very short introduction to the Finite Element Method
A very short introduction to the Finite Element Method Till Mathis Wagner Technical University of Munich JASS 2004, St Petersburg May 4, 2004 1 Introduction This is a short introduction to the finite element
More informationA Finite Element Method Using Singular Functions for Poisson Equations: Mixed Boundary Conditions
A Finite Element Method Using Singular Functions for Poisson Equations: Mixed Boundary Conditions Zhiqiang Cai Seokchan Kim Sangdong Kim Sooryun Kong Abstract In [7], we proposed a new finite element method
More informationA FV Scheme for Maxwell s equations
A FV Scheme for Maxwell s equations Convergence Analysis on unstructured meshes Stephanie Lohrengel * Malika Remaki ** *Laboratoire J.A. Dieudonné (UMR CNRS 6621), Université de Nice Sophia Antipolis,
More informationHIGHER ORDER FINITE-ELEMENT MODELLING FOR CURVED INTERFACE PROBLEMS BASED ON STRAIGHT- EDGED N-SIMPLEXES AND POLYNOMIAL ALGEBRA
6th European Conference on Computational Mechanics (ECCM 6) 7th European Conference on Computational Fluid Dynamics (ECFD 7) 1115 June 2018, Glasgow, UK HIGHER ORDER FINITE-ELEMENT MODELLING FOR CURVED
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 informationLecture 2. Introduction to FEM. What it is? What we are solving? Potential formulation Why? Boundary conditions
Introduction to FEM What it is? What we are solving? Potential formulation Why? Boundary conditions Lecture 2 Notation Typical notation on the course: Bolded quantities = matrices (A) and vectors (a) Unit
More informationA Balancing Algorithm for Mortar Methods
A Balancing Algorithm for Mortar Methods Dan Stefanica Baruch College, City University of New York, NY 11, USA Dan Stefanica@baruch.cuny.edu Summary. The balancing methods are hybrid nonoverlapping Schwarz
More informationVolume and surface integral equations for electromagnetic scattering by a dielectric body
Volume and surface integral equations for electromagnetic scattering by a dielectric body M. Costabel, E. Darrigrand, and E. H. Koné IRMAR, Université de Rennes 1,Campus de Beaulieu, 35042 Rennes, FRANCE
More informationScientific Computing WS 2017/2018. Lecture 18. Jürgen Fuhrmann Lecture 18 Slide 1
Scientific Computing WS 2017/2018 Lecture 18 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 18 Slide 1 Lecture 18 Slide 2 Weak formulation of homogeneous Dirichlet problem Search u H0 1 (Ω) (here,
More informationA note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations
A note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations S. Hussain, F. Schieweck, S. Turek Abstract In this note, we extend our recent work for
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 informationMimetic Finite Difference methods
Mimetic Finite Difference methods An introduction Andrea Cangiani Università di Roma La Sapienza Seminario di Modellistica Differenziale Numerica 2 dicembre 2008 Andrea Cangiani (IAC CNR) mimetic finite
More informationc 2008 Society for Industrial and Applied Mathematics
SIAM J. NUMER. ANAL. Vol. 46, No. 3, pp. 640 65 c 2008 Society for Industrial and Applied Mathematics A WEIGHTED H(div) LEAST-SQUARES METHOD FOR SECOND-ORDER ELLIPTIC PROBLEMS Z. CAI AND C. R. WESTPHAL
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 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 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 informationA Balancing Algorithm for Mortar Methods
A Balancing Algorithm for Mortar Methods Dan Stefanica Baruch College, City University of New York, NY, USA. Dan_Stefanica@baruch.cuny.edu Summary. The balancing methods are hybrid nonoverlapping Schwarz
More informationH(div) and H(curl)-conforming Virtual Element Methods
Noname manuscript No. (will be inserted by the editor) H(div) and H(curl)-conforming Virtual Element Methods L. Beirão da Veiga F. Brezzi L. D. Marini A. Russo the date of receipt and acceptance should
More informationAnisotropic mesh refinement in polyhedral domains: error estimates with data in L 2 (Ω)
Anisotropic mesh refinement in polyhedral domains: error estimates with data in L 2 (Ω) Thomas Apel Ariel L. Lombardi Max Winkler February 6, 2014 arxiv:1303.2960v1 [math.na] 12 Mar 2013 Abstract. The
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 informationExplicit kernel-split panel-based Nyström schemes for planar or axisymmetric Helmholtz problems
z Explicit kernel-split panel-based Nyström schemes for planar or axisymmetric Helmholtz problems Johan Helsing Lund University Talk at Integral equation methods: fast algorithms and applications, Banff,
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 informationMULTIGRID METHODS FOR MAXWELL S EQUATIONS
MULTIGRID METHODS FOR MAXWELL S EQUATIONS A Dissertation Submitted to the Graduate Faculty of the Louisiana State University and Agricultural and Mechanical College in partial fulfillment of the requirements
More informationGraded Mesh Refinement and Error Estimates of Higher Order for DGFE-solutions of Elliptic Boundary Value Problems in Polygons
Graded Mesh Refinement and Error Estimates of Higher Order for DGFE-solutions of Elliptic Boundary Value Problems in Polygons Anna-Margarete Sändig, Miloslav Feistauer University Stuttgart, IANS Journées
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 informationA u + b u + cu = f in Ω, (1.1)
A WEIGHTED H(div) LEAST-SQUARES METHOD FOR SECOND-ORDER ELLIPTIC PROBLEMS Z. CAI AND C. R. WESTPHAL Abstract. This paper presents analysis of a weighted-norm least squares finite element method for elliptic
More informationR T (u H )v + (2.1) J S (u H )v v V, T (2.2) (2.3) H S J S (u H ) 2 L 2 (S). S T
2 R.H. NOCHETTO 2. Lecture 2. Adaptivity I: Design and Convergence of AFEM tarting with a conforming mesh T H, the adaptive procedure AFEM consists of loops of the form OLVE ETIMATE MARK REFINE to produce
More informationfield using second order edge elements in 3D
The current issue and full text archive of this journal is available at http://www.emerald-library.com using second order edge elements in 3D Z. Ren Laboratoire de GeÂnie Electrique de Paris, UniversiteÂs
More informationFinite Element Method (FEM)
Finite Element Method (FEM) The finite element method (FEM) is the oldest numerical technique applied to engineering problems. FEM itself is not rigorous, but when combined with integral equation techniques
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 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 informationMATHEMATICS. Course Syllabus. Section A: Linear Algebra. Subject Code: MA. Course Structure. Ordinary Differential Equations
MATHEMATICS Subject Code: MA Course Structure Sections/Units Section A Section B Section C Linear Algebra Complex Analysis Real Analysis Topics Section D Section E Section F Section G Section H Section
More informationDivergence-conforming multigrid methods for incompressible flow problems
Divergence-conforming multigrid methods for incompressible flow problems Guido Kanschat IWR, Universität Heidelberg Prague-Heidelberg-Workshop April 28th, 2015 G. Kanschat (IWR, Uni HD) Hdiv-DG Práha,
More informationCorner singularities of Maxwell interface and eddy current problems
Corner singularities of Maxwell interface and eddy current problems Martin Costabel, Monique Dauge, and Serge Nicaise Abstract. The aim of this work is to provide a description of the corner asymptotics
More informationAdaptive C1 Macroelements for Fourth Order and Divergence-Free Problems
Adaptive C1 Macroelements for Fourth Order and Divergence-Free Problems Roy Stogner Computational Fluid Dynamics Lab Institute for Computational Engineering and Sciences University of Texas at Austin March
More informationScientific Computing: Numerical Integration
Scientific Computing: Numerical Integration Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 Course MATH-GA.2043 or CSCI-GA.2112, Fall 2015 Nov 5th, 2015 A. Donev (Courant Institute) Lecture
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 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 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 informationBasic Concepts of Adaptive Finite Element Methods for Elliptic Boundary Value Problems
Basic Concepts of Adaptive Finite lement Methods for lliptic Boundary Value Problems Ronald H.W. Hoppe 1,2 1 Department of Mathematics, University of Houston 2 Institute of Mathematics, University of Augsburg
More informationComputations of the first eigenpairs for the Schrödinger operator with magnetic field
Computations of the first eigenpairs for the Schrödinger operator with magnetic field V. Bonnaillie-Noël, M. Dauge, D. Martin and G. Vial Abstract This paper is devoted to computations of eigenvalues and
More informationStatic Maxwell Type Problems: Functional A Posteriori Error Estimates and Estimates for the Maxwell Constant in 3D. Dirk Pauly
Static Maxwell Type Problems: Functional A Posteriori Error Estimates and Estimates for the Maxwell Constant in 3D Dirk Pauly Fakultät für Mathematik Universität Duisburg-Essen, Campus Essen, Germany partially
More informationNEUMANN AND MIXED PROBLEMS ON CURVILINEAR POLYHEDRA
NEUMANN AND MIXED PROBLEMS ON CURVILINEAR POLYHEDRA Monique Dauge We prove regularity results in L p Sobolev spaces. On one hand, we state some abstract results by L p functional techniques : exponentially
More informationSingularly Perturbed Partial Differential Equations
WDS'9 Proceedings of Contributed Papers, Part I, 54 59, 29. ISN 978-8-7378--9 MTFYZPRESS Singularly Perturbed Partial Differential Equations J. Lamač Charles University, Faculty of Mathematics and Physics,
More informationThe Imaging of Anisotropic Media in Inverse Electromagnetic Scattering
The Imaging of Anisotropic Media in Inverse Electromagnetic Scattering Fioralba Cakoni Department of Mathematical Sciences University of Delaware Newark, DE 19716, USA email: cakoni@math.udel.edu Research
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationarxiv: v1 [math.na] 8 Feb 2018
arxiv:180.094v1 [math.na] 8 Feb 018 The nonconforming virtual element method for eigenvalue problems F. Gardini a, G. Manzini b, and G. Vacca c a Dipartimento di Matematica F. Casorati, Università di Pavia,
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 informationSimple Examples on Rectangular Domains
84 Chapter 5 Simple Examples on Rectangular Domains In this chapter we consider simple elliptic boundary value problems in rectangular domains in R 2 or R 3 ; our prototype example is the Poisson equation
More informationElectromagnetic wave propagation. ELEC 041-Modeling and design of electromagnetic systems
Electromagnetic wave propagation ELEC 041-Modeling and design of electromagnetic systems EM wave propagation In general, open problems with a computation domain extending (in theory) to infinity not bounded
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 informationElectromagnetic Forces on Parallel Current-
Page 1 of 5 Tutorial Models : Electromagnetic Forces on Parallel Current-Carrying Wires Electromagnetic Forces on Parallel Current- Carrying Wires Introduction One ampere is defined as the constant current
More information1 Fundamentals. 1.1 Overview. 1.2 Units: Physics 704 Spring 2018
Physics 704 Spring 2018 1 Fundamentals 1.1 Overview The objective of this course is: to determine and fields in various physical systems and the forces and/or torques resulting from them. The domain of
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 informationSome New Elements for the Reissner Mindlin Plate Model
Boundary Value Problems for Partial Differential Equations and Applications, J.-L. Lions and C. Baiocchi, eds., Masson, 1993, pp. 287 292. Some New Elements for the Reissner Mindlin Plate Model Douglas
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 informationIsogeometric modeling of Lorentz detuning in linear particle accelerator cavities
Isogeometric modeling of Lorentz detuning in linear particle accelerator cavities Mauro Bonafini 1,2, Marcella Bonazzoli 3, Elena Gaburro 1,2, Chiara Venturini 1,3, instructor: Carlo de Falco 4 1 Department
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 informationTechnische Universität Graz
Technische Universität Graz A note on the stable coupling of finite and boundary elements O. Steinbach Berichte aus dem Institut für Numerische Mathematik Bericht 2009/4 Technische Universität Graz A
More informationA DELTA-REGULARIZATION FINITE ELEMENT METHOD FOR A DOUBLE CURL PROBLEM WITH DIVERGENCE-FREE CONSTRAINT
A DELTA-REGULARIZATION FINITE ELEMENT METHOD FOR A DOUBLE CURL PROBLEM WITH DIVERGENCE-FREE CONSTRAINT HUOYUAN DUAN, SHA LI, ROGER C. E. TAN, AND WEIYING ZHENG Abstract. To deal with the divergence-free
More informationMultilevel Preconditioning of Graph-Laplacians: Polynomial Approximation of the Pivot Blocks Inverses
Multilevel Preconditioning of Graph-Laplacians: Polynomial Approximation of the Pivot Blocks Inverses P. Boyanova 1, I. Georgiev 34, S. Margenov, L. Zikatanov 5 1 Uppsala University, Box 337, 751 05 Uppsala,
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 informationTechnische Universität Graz
Technische Universität Graz Stability of the Laplace single layer boundary integral operator in Sobolev spaces O. Steinbach Berichte aus dem Institut für Numerische Mathematik Bericht 2016/2 Technische
More informationIntroduction. Finite and Spectral Element Methods Using MATLAB. Second Edition. C. Pozrikidis. University of Massachusetts Amherst, USA
Introduction to Finite and Spectral Element Methods Using MATLAB Second Edition C. Pozrikidis University of Massachusetts Amherst, USA (g) CRC Press Taylor & Francis Group Boca Raton London New York CRC
More informationarxiv: v1 [math.na] 11 Jul 2011
Multigrid Preconditioner for Nonconforming Discretization of Elliptic Problems with Jump Coefficients arxiv:07.260v [math.na] Jul 20 Blanca Ayuso De Dios, Michael Holst 2, Yunrong Zhu 2, and Ludmil Zikatanov
More informationLocking in the incompressible limit: pseudo-divergence-free element free Galerkin
Locking in the incompressible limit: pseudodivergencefree element free Galerkin Yolanda Vidal Pierre Villon and Antonio Huerta Departament de Matemàtica Aplicada III, ES de Ingenieros de Caminos, Canales
More informationA Framework for Analyzing and Constructing Hierarchical-Type A Posteriori Error Estimators
A Framework for Analyzing and Constructing Hierarchical-Type A Posteriori Error Estimators Jeff Ovall University of Kentucky Mathematics www.math.uky.edu/ jovall jovall@ms.uky.edu Kentucky Applied and
More informationDivergence Boundary Conditions for Vector Helmholtz Equations with Divergence Constraints
Divergence Boundary Conditions for Vector Helmholtz Equations with Divergence Constraints Urve Kangro Department of Mathematical Sciences Carnegie Mellon University Pittsburgh, PA 15213 Roy Nicolaides
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 informationNumerical methods for PDEs FEM convergence, error estimates, piecewise polynomials
Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen Numerical methods for PDEs FEM convergence, error estimates, piecewise polynomials Dr. Noemi Friedman Contents of the course Fundamentals
More information8 A pseudo-spectral solution to the Stokes Problem
8 A pseudo-spectral solution to the Stokes Problem 8.1 The Method 8.1.1 Generalities We are interested in setting up a pseudo-spectral method for the following Stokes Problem u σu p = f in Ω u = 0 in Ω,
More information