Institut de Recherche MAthématique de Rennes
|
|
- Anna Skinner
- 5 years ago
- Views:
Transcription
1 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 DAUGE Institut de Recherche MAthématique de Rennes
2 References Theory and analysis of the h -version: M. Costabel, M. Dauge: Weighted regularization of Maxwell equations in polyhedral domains. Numer. Math. Online publication DOI.7/s88 () Numerical tests, including hp version, benchmark problems: Numerical results in D: M. Costabel, M. Dauge, D. Martin, G. Vial: Weighted regularization of Maxwell equations computations in curvilinear polygons. Proceedings of ENUMATH () Our Finite Element package: D. Martin: Mélina. Online documentation: Analysis of the hp version: M. Costabel, M. Dauge, Ch. Schwab: Exponential convergence of the hp -FEM for the weighted regularization of Maxwell equations in polygonal domains. In preparation
3 Convergence of Eigenvalues:
4 The Ideal... Valeur propre Neumann n sur le Maillage - Interp Interp Interp Interp Interp Interp Interp 7 Interp st Neumann eigenvalue on a curved rectangle Nodal elements Two elements p version Interp 9 - Relative error vs # DOF
5 ... and the reality Valeur propre Maxwell n pénalisée au bord par λ = ; Maillage Interp Interp Interp st Maxwell Interp Interp eigenvalue on a curved L shape Nodal elements Regularization without weight - Relative error vs # DOF
6 The standard nodal approximation of the first singularity The exact solution ( grad(r sin θ )). Computation with Q elements, regularization without weight. Non- H singularities are not approximated.
7 The Ideal Recovered - - Geo. deg. Geo. deg st Maxwell eigenvalue on a curved L shape Geo. deg. 7 Nodal elements - Geo. deg. Geo. deg. Geo. deg Weighted regularization -8 Relative error vs # DOF
8 Outline The Maxwell source problem and eigenproblem and their variational formulation with weighted regularization. Convergence analysis of the weighted regularization method: The h version The hp version Numerical results in D: Straight L-shaped domain Curved L-shaped domain Evidence for exponential convergence Numerical results in D: The thick L Fichera s corner. 7
9 Maxwell source problem Permittivity ε and permeability µ set to. Ω polyhedron in R (or polygon in R ). Given ω R, J L (Ω) with div J =, find (E, H) L (Ω) such that curl E iωh = & curl H + iωe = J in Ω, (Maxwell SP) E n = & H n = on Ω Maxwell eigenvalue problem Find non-zero ω such that (E, H) curl E iωh = & curl H + iωe = (Maxwell EVP) E n = & H n = on Ω in Ω, In both cases, div E = div B = is implied if ω. 8
10 Variational formulation:the WRM Bilinear form for s >, < α, a s (u, v) = Ω curl u curl v + s Ω rα div u div v, where r(x) dist(x, { nonconvex edges and corners }). Find u X such that v X : a s (u, v) ω Ω u v = Ω f v (SP): ω, f given, u unknown; (EVP): f = ; u, ω unknown. Energy space X = X γ N := H (curl) {u L (Ω), r γ div u L (Ω)}. with γ = α/. Recall: σ(p [s] rγ ) = σ(m) s σ(ldir r ). γ 9
11 Galerkin approximation: Everything standard here With the test and trial subspace X h X, the Galerkin solution u h X h is defined by v h X h : a s (u h, v h ) ω Ω u h v h = Ω f v h Since a s (, ) is positive definite and symmetric, we have the usual stability and error estimates, in particular: Céa s lemma: For the source problem and ω = : u u h X γ N inf{ u v h X γ N v h X h } Convergence of the eigenvalues: For the eigenvalue problem, f = : Suppose {X h } h> is a family of subspaces approximating X. If ω σ(p [s] r γ discrete problem such that ω ω h = O((energy norm error) )). The rest is approximation theory... is a simple eigenvalue, then there is a simple eigenvalue ω h of the
12 Birman-Solomyak decompositions Basic idea for the estimate of the approximation error:. Decomposition of the solution u = w + grad ϕ, with w regular and ϕ carrying the main corner singularities Energy norm: u X γ N w H + grad ϕ L + rγ ϕ L. Approximation of u = w + grad ϕ by u h = w h + grad ϕ h, u u h X γ N w w h H + ϕ ϕ h V γ, where w h X h, and ϕ h belongs to an auxiliary space Φ h satisfying Note: grad Φ h X h. Since X h consists of globally continuous functions, Φ h must consist of globally C functions. The space Φ h is not used in the computations, but only in the estimation of the approximation error.
13 B-S decomposition suitable for the h version Suitable hypotheses on the weight, smooth right hand side (or EV problem), D,... Theorem: The solution can be decomposed as u = w + grad ϕ, with w H +τ (Ω) and ϕ m= K m β Here (D, one reentrant corner of opening ω ) H +τ is the regularity of the Neumann problem for the Laplace operator in Ω, i.e. τ < π/ω, and β is determined by the regularity of the Dirichlet problem, β > π/ω. The weighted Sobolev spaces K m β ϕ K m β are defined by the seminorms = α =m r β+ α α ϕ L Result: Energy norm estimate of order O(h τ ) for any mesh that allows C elements.
14 B-S decomposition suitable for the hp version At present: Ω has to be a polygon in R. Weighted spaces of analytic functions, compare BABUŠKA-GUO: A β (Ω) m= K m β (Ω), C : m : ϕ K m β Cm+ m! One assumes that the right hand side is analytic (true for EVP), for example in A (Ω). Theorem: For γ (γ Ω, ] there is a δ γ > such that for all < δ < δ γ, the solution u can be decomposed as u = w + grad ϕ, with w A δ (Ω) and ϕ H A γ δ(ω) Result: For suitably (geometrically) refined meshes, one has an energy norm error estimate of order O(e bp ) for some b >.
15 For the analysis of the hp version, types of exponentially refined meshes: In CYAN, the patch to consider for the construction of C interpolants.
16 Rectangles with hanging nodes
17 Q -quadrilaterals
18 Triangles
19 Problems with Q -quadrilaterals Polynomials on the reference square K are not polynomials on the physical element K, in general. The Jacobian J P ( K) is not constant, therefore the gradient (in physical space) does not map Q p ( K) into Q p ( K). The (auxiliary) approximation FE space Φ h for the potentials has to be chosen not as Q p ( K), but as J m Q p m ( K) (depends on K!), for some m. There is a geometrical condition on the patches (the jacobians J on both sides have to be proportional on the interface) which is not always satisfied, but satisfied in our case of the L and two axiparallel sides. a a b b a b = a b 8
20 Computations in D The straight L: Ω is the L -shaped domain : [., ] [., ] \ [.7, ] [.7, ]. Computation of the first eigenvalues Λ[s] for s = [ : ], α =,,, q =,,, The computed eigenvalues are sorted depending on whether curl u s r α/ div u is ρ or ρ with a fixed ρ. 9
21 Legend For the next 8 figures. In Cyan, the Laplace-Neumann eigenvalues, which coincide for D domains with Maxwell eigenvalues. In Blue, the circles represent the computed Λ[s] with curl -dominant eigenvectors. In Red, the lines join the computed Λ[s] with div -dominant eigenvectors. Abscissa : s / Ordinate : Λ[s].
22 q =, α = 7 8 9
23 q =, α = 7 8 9
24 q =, α = 7 8 9
25 q =, α = 7 8 9
26 q =, α = 7 8 9
27 q =, α = 7 8 9
28 q =, α =
29 q =, α =
30 Geometrical mesh refinement, on straight L-shape domain (sl) Ratio layers. layers. 9
31 Geometrical mesh refinement, on straight L-shape domain (sl) Ratio layers. layers.
32 First Maxwell ev. on L-domain, with α = and ref. ratio On domain (sl) with number of layers m & funct. interp. p. - - Nb of layers Nb of layers Nb of layers Nb of layers Relative error versus the number of DOF. Computed with MÉLINA FEM library Nb of layers 7 developped by - Nb of layers 8 D. MARTIN Nb of layers 9 (Rennes, ENSTA). -8 Nb of layers
33 First Maxwell ev. on L-domain, with α = and ref. ratio On domain (sl) with number of layers m & funct. interp. p Nb of layers Nb of layers Nb of layers Nb of layers Nb of layers Nb of layers Nb of layers 9 9 Nb of layers Relative error versus the number of DOF. Computed with MÉLINA FEM library developped by D. MARTIN (Rennes, ENSTA).
34 Second Maxwell ev. on L-domain, with α = and ref. ratio On domain (sl) with number of layers m & funct. interp. p Nb of layers Nb of layers Nb of layers Nb of layers Nb of layers 7 Nb of layers 8 Nb of layers 9 Nb of layers Relative error versus the number of DOF. Computed with MÉLINA FEM library developped by D. MARTIN (Rennes, ENSTA).
35 Second Maxwell ev. on L-domain, with α = and ref. ratio On domain (sl) with number of layers m & funct. interp. p Nb of layers Nb of layers Nb of layers Nb of layers Nb of layers 7 Nb of layers 8 Nb of layers 9 Nb of layers Relative error versus the number of DOF. Computed with MÉLINA FEM library developped by D. MARTIN (Rennes, ENSTA).
36 Curved L-shaped domain..... The curved L-shaped domain (cl) with a π/ corner.
37 Meshes on the approximate curved L-domain Ω q,m, q =, m =,
38 First Maxwell ev. on curved L, with weight & refinement On domain (cl) with ratio, number of layers m =, and α = Geo. deg. Geo. deg. Geo. deg. Geo. deg. Geo. deg. Geo. deg Geom. interp. n & Funct. interp. p. Relative error versus the number of DOF. Computed with MÉLINA
39 Second Maxwell ev. on curved L, with weight & refinement On domain (cl) with ratio, number of layers m = and α = Geo. deg Geo. deg. Geo. deg. Geo. deg. Geo. deg. Geo. deg. Geom. interp. n & Funct. interp. p. Relative error versus the number of DOF. Computed with MÉLINA. Something weird going on here... 8
40 First Maxwell ev. on curved L, with weight & refinement On domain (cl) with ratio, number of layers m =, and α =. - - Geo. deg. Geo. deg Nodes = Gauss-Lobatto points. Geom. interp. n & Funct. interp. Geo. deg. 7 p. - Geo. deg. Geo. deg Relative error Geo. deg. 9 9 versus the # of DOF. -8 9
41 Second Maxwell ev. on curved L, with weight & refinement On domain (cl) with ratio, number of layers m = and α = Geo. deg. Geo. deg. Geo. deg Geo. deg. 7 Geo. deg. Geo. deg Nodes = Gauss-Lobatto points. Geom. interp. n & Funct. interp. p. Relative error versus the # of DOF....this is much better!
42 Exponential convergence Domain (sl) with ratio, α =, lowest eigenvalues eigenvalue No eigenvalue No eigenvalue No eigenvalue No h p version: eigenvalue No p = #layers Relative error vs #DOF
43 Exponential convergence Domain (sl) with ratio, α =, lowest eigenvalues h p version: p = #layers.. eigenvalue No eigenvalue No eigenvalue No log(-log(error)) vs eigenvalue No eigenvalue No log(#dof)
44 Exponential convergence Domain (sl) with ratio, α =, lowest eigenvalues eigenvalue No eigenvalue No eigenvalue No eigenvalue No eigenvalue No h p version: p = #layers log(error) vs (#DOF) /
45 D hp version: Computations are working nicely, Theory: not yet finished... In the following graphs, we show results for two D domains: The thick L L [, ], where L = [, ] \ [, ] is the D straight L-shaped domain. Fichera s corner (7/8th cube) [, ] \ [, ]
46 The Thick L: Legend The Maxwell eigenvalues of the thick L L [, ] are known: They are sums of Dirichlet and Neumann eigenvalues of the D domain L (which can easily be computed to high precision) and of the interval [, ] (which are of the form π n, n N ). The following 9 graphs show the error as a function of the number of degrees of freedom for the first 8 Maxwell eigenvalues and computed convergence rates for the th eigenvalue. The computations were done for Q p elements on a Q -hexahedral mesh, obtained by using a geometrically refined mesh on L with refinement ratio 8 and a variable number of layers indicated by color, and two layers in the third dimension. The degree p varies from to. For the weight sr α, we used s =, α =, and r the distance to the singular edge {(, )} [, ].
47 The Thick L: Maxwell Eigenvalue vp Maxwell (D) nº -- Ratio de maillage Nb of layers Nb of layers Nb of layers 7 7 rel. error vs #dof - - Nb of layers Nb of layers Nb of layers Nb of layers 7 p #layers 7 - -
48 The Thick L: Maxwell Eigenvalue vp Maxwell (D) nº -- Ratio de maillage Nb of layers Nb of layers Nb of layers Nb of layers Nb of layers Nb of layers Nb of layers rel. error vs #dof p #layers
49 The Thick L: Maxwell Eigenvalue vp Maxwell (D) nº -- Ratio de maillage Nb of layers Nb of layers rel. error vs #dof - - Nb of layers Nb of layers Nb of layers Nb of layers Nb of layers 7 7 p #layers
50 The Thick L: Maxwell Eigenvalue vp Maxwell (D) nº -- Ratio de maillage Nb of layers Nb of layers Nb of layers Nb of layers Nb of layers Nb of layers Nb of layers rel. error vs #dof p #layers
51 The Thick L: Maxwell Eigenvalue vp Maxwell (D) nº -- Ratio de maillage convergence rates p #layers
52 The Thick L: Maxwell Eigenvalue vp Maxwell (D) nº -- Ratio de maillage 8 - Nb of layers Nb of layers rel. error vs #dof - - Nb of layers Nb of layers Nb of layers Nb of layers Nb of layers p #layers 7 - -
53 The Thick L: Maxwell Eigenvalue vp Maxwell (D) nº -- Ratio de maillage Nb of layers Nb of layers Nb of layers Nb of layers Nb of layers Nb of layers Nb of layers rel. error vs #dof p #layers 7-8 -
54 The Thick L: Maxwell Eigenvalue 7 vp Maxwell (D) nº 7 -- Ratio de maillage Nb of layers Nb of layers Nb of layers Nb of layers Nb of layers Nb of layers Nb of layers rel. error vs #dof p #layers 7-8 -
55 The Thick L: Maxwell Eigenvalue 8 vp Maxwell (D) nº 8 -- Ratio de maillage Nb of layers Nb of layers Nb of layers 7 7 rel. error vs #dof - - Nb of layers Nb of layers Nb of layers Nb of layers 7 p #layers 7 - -
56 Fichera Corner The eigenvalues for this domain are not known. Therefore we show the validity of our computations with the weighted regularization method in two ways: Dependence of the eigenvalues on the parameter s in the weight sr α : We compute the first Maxwell eigenvalues for s ranging from to. The curl dominant eigenvalues are shown in blue, the div dominant eigenvalues in red, and the indifferent ones in green (factor ρ =. ). We see that for α = the eigenvalues show the required behavior, i.e. they lie on two families of straight lines. Without weight ( α = ), the behavior is different. Approximation of the cube. We deform the 7/8th cube so that it approximates the unit cube [, ]. The Maxwell eigenvalues for the latter are known ( π times sums of squares of integers). In our graphs, we show the eigenvalues divided by π to make this approximation more evident.
57 Fichera Corner cont d We show computations for domains: F, the 7/8th cube [, ] \ [, ]. The eigenvalues have high multiplicity due to symmetries. F = F ([.8, ] [., ] [, ]) Here the symmetries are broken. From F to F and from F to F, the thickness of the protruding hexahedra is reduced by a factor of, so that F is % away from the unit cube [, ] in Hausdorff distance. We use Q hexahedral elements on a tensor product mesh generated from the -dimensional mesh {, /, /, } [, ]. For readability reasons, in the following picture the refinement ratio / is replaced by / and the geometric approximation factor / by /.
58 Fichera Corner F [a,b,c] = [,,] Fichera Corner F [a,b,c] = [.,.8,] The domains F... F Fichera Corner F [a,b,c] = [.,.,.] Fichera Corner F [a,b,c] = [.,.,.]
59 Ratio de tri. et α =. Domaine F F, no weight: α =
60 Ratio de tri. et α =. Domaine F F, α =
61 Ratio de tri. et α =. Domaine F F, α =...
62 7 Ratio de tri. et α =. Domaine F F, α =
63 8 Ratio de tri. et α =. Domaine F F, α = 7
64 Conclusion The Weighted Regularizaton Method requires only a very simple modification of the standard regularized variational formulation of time-harmonic Maxwell source and eigenvalue problems. It allows the use of standard nodal finite elements even in the case of domains with reentrant corners and solutions with strong non- H singularities where the standard regularization gives incorrect results. In combination with the hp version of the finite element method, an efficient and robust algorithm is obtained that allows highly precise approximations even for D problems.
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 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 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 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 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 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 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 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 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 informationMesh Grading towards Singular Points Seminar : Elliptic Problems on Non-smooth Domain
Mesh Grading towards Singular Points Seminar : Elliptic Problems on Non-smooth Domain Stephen Edward Moore Johann Radon Institute for Computational and Applied Mathematics Austrian Academy of Sciences,
More informationTrefftz-Discontinuous Galerkin Methods for Maxwell s Equations
Trefftz-Discontinuous Galerkin Methods for Maxwell s Equations Ilaria Perugia Faculty of Mathematics, University of Vienna Vienna P D E Ralf Hiptmair (ETH Zürich), Andrea Moiola (University of Reading)
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 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 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 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 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 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 informationThe Plane Stress Problem
The Plane Stress Problem Martin Kronbichler Applied Scientific Computing (Tillämpad beräkningsvetenskap) February 2, 2010 Martin Kronbichler (TDB) The Plane Stress Problem February 2, 2010 1 / 24 Outline
More informationTrefftz-discontinuous Galerkin methods for time-harmonic wave problems
Trefftz-discontinuous Galerkin methods for time-harmonic wave problems Ilaria Perugia Dipartimento di Matematica - Università di Pavia (Italy) http://www-dimat.unipv.it/perugia Joint work with Ralf Hiptmair,
More information3. Numerical integration
3. Numerical integration... 3. One-dimensional quadratures... 3. Two- and three-dimensional quadratures... 3.3 Exact Integrals for Straight Sided Triangles... 5 3.4 Reduced and Selected Integration...
More informationMaximum norm estimates for energy-corrected finite element method
Maximum norm estimates for energy-corrected finite element method Piotr Swierczynski 1 and Barbara Wohlmuth 1 Technical University of Munich, Institute for Numerical Mathematics, piotr.swierczynski@ma.tum.de,
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 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 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 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 informationNumerical Solutions of Laplacian Problems over L-Shaped Domains and Calculations of the Generalized Stress Intensity Factors
WCCM V Fifth World Congress on Computational Mechanics July 7-2, 2002, Vienna, Austria Eds.: H.A. Mang, F.G. Rammerstorfer, J. Eberhardsteiner Numerical Solutions of Laplacian Problems over L-Shaped Domains
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 informationSTOKES PROBLEM WITH SEVERAL TYPES OF BOUNDARY CONDITIONS IN AN EXTERIOR DOMAIN
Electronic Journal of Differential Equations, Vol. 2013 2013, No. 196, pp. 1 28. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu STOKES PROBLEM
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 informationPreconditioned space-time boundary element methods for the heat equation
W I S S E N T E C H N I K L E I D E N S C H A F T Preconditioned space-time boundary element methods for the heat equation S. Dohr and O. Steinbach Institut für Numerische Mathematik Space-Time Methods
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 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 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 informationFinal Ph.D. Progress Report. Integration of hp-adaptivity with a Two Grid Solver: Applications to Electromagnetics. David Pardo
Final Ph.D. Progress Report Integration of hp-adaptivity with a Two Grid Solver: Applications to Electromagnetics. David Pardo Dissertation Committee: I. Babuska, L. Demkowicz, C. Torres-Verdin, R. Van
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 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 informationMath 660-Lecture 15: Finite element spaces (I)
Math 660-Lecture 15: Finite element spaces (I) (Chapter 3, 4.2, 4.3) Before we introduce the concrete spaces, let s first of all introduce the following important lemma. Theorem 1. Let V h consists of
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 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 informationhp-dgfem FOR SECOND-ORDER MIXED ELLIPTIC PROBLEMS IN POLYHEDRA
MATHEMATICS OF COMPUTATION Volume 00, Number 0, Pages 000 000 S 0025-5718(XX)0000-0 hp-dgfem FOR SECOND-ORDER MIXED ELLIPTIC PROBLEMS IN POLYHEDRA DOMINIK SCHÖTZAU, CHRISTOPH SCHWAB, AND THOMAS P. WIHLER
More informationThomas Apel 1, Ariel L. Lombardi 2 and Max Winkler 1
Mathematical Modelling and Numerical Analysis Modélisation Mathématique et Analyse Numérique Will be set by the publisher ANISOTROPIC MESH REFINEMENT IN POLYHEDRAL DOMAINS: ERROR ESTIMATES WITH DATA IN
More informationApproximation by plane and circular waves
LMS EPSRC DURHAM SYMPOSIUM, 8 16TH JULY 2014 Building Bridges: Connections and Challenges in Modern Approaches to Numerical PDEs Approximation by plane and circular waves Andrea Moiola DEPARTMENT OF MATHEMATICS
More informationTrefftz-Discontinuous Galerkin Methods for the Time-Harmonic Maxwell Equations
Trefftz-Discontinuous Galerkin Methods for the Time-Harmonic Maxwell Equations Ilaria Perugia Dipartimento di Matematica - Università di Pavia (Italy) http://www-dimat.unipv.it/perugia Joint work with
More informationA non-standard Finite Element Method based on boundary integral operators
A non-standard Finite Element Method based on boundary integral operators Clemens Hofreither Ulrich Langer Clemens Pechstein June 30, 2010 supported by Outline 1 Method description Motivation Variational
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 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 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 informationGeoPDEs. An Octave/Matlab software for research on IGA. R. Vázquez. IMATI Enrico Magenes, Pavia Consiglio Nazionale della Ricerca
GeoPDEs An Octave/Matlab software for research on IGA R. Vázquez IMATI Enrico Magenes, Pavia Consiglio Nazionale della Ricerca Joint work with C. de Falco and A. Reali Supported by the ERC Starting Grant:
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 informationIsogeometric Analysis:
Isogeometric Analysis: some approximation estimates for NURBS L. Beirao da Veiga, A. Buffa, Judith Rivas, G. Sangalli Euskadi-Kyushu 2011 Workshop on Applied Mathematics BCAM, March t0th, 2011 Outline
More informationApproximation of the LBB constant on corner domains
Approximation of the LBB constant on corner domains Martin Costabel Collaboration with Monique Dauge, Michel Crouzeix, Christine Bernardi, Vivette Girault, Yvon Lafranche IRMAR, Université de Rennes 1
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 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 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 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 informationOn an Approximation Result for Piecewise Polynomial Functions. O. Karakashian
BULLETIN OF THE GREE MATHEMATICAL SOCIETY Volume 57, 010 (1 7) On an Approximation Result for Piecewise Polynomial Functions O. arakashian Abstract We provide a new approach for proving approximation results
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 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 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 information1. Introduction. We consider the model problem that seeks an unknown function u = u(x) satisfying
A SIMPLE FINITE ELEMENT METHOD FOR LINEAR HYPERBOLIC PROBLEMS LIN MU AND XIU YE Abstract. In this paper, we introduce a simple finite element method for solving first order hyperbolic equations with easy
More informationCoupling of eddy-current and circuit problems
Coupling of eddy-current and circuit problems Ana Alonso Rodríguez*, ALBERTO VALLI*, Rafael Vázquez Hernández** * Department of Mathematics, University of Trento ** Department of Applied Mathematics, University
More informationPartial Differential Equations and the Finite Element Method
Partial Differential Equations and the Finite Element Method Pavel Solin The University of Texas at El Paso Academy of Sciences ofthe Czech Republic iwiley- INTERSCIENCE A JOHN WILEY & SONS, INC, PUBLICATION
More informationWEAK GALERKIN FINITE ELEMENT METHODS ON POLYTOPAL MESHES
INERNAIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 12, Number 1, Pages 31 53 c 2015 Institute for Scientific Computing and Information WEAK GALERKIN FINIE ELEMEN MEHODS ON POLYOPAL MESHES LIN
More informationTime domain boundary elements for dynamic contact problems
Time domain boundary elements for dynamic contact problems Heiko Gimperlein (joint with F. Meyer 3, C. Özdemir 4, D. Stark, E. P. Stephan 4 ) : Heriot Watt University, Edinburgh, UK 2: Universität Paderborn,
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 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 informationFinite Element Error Estimates in Non-Energy Norms for the Two-Dimensional Scalar Signorini Problem
Journal manuscript No. (will be inserted by the editor Finite Element Error Estimates in Non-Energy Norms for the Two-Dimensional Scalar Signorini Problem Constantin Christof Christof Haubner Received:
More informationAdaptive Finite Element Methods Lecture Notes Winter Term 2017/18. R. Verfürth. Fakultät für Mathematik, Ruhr-Universität Bochum
Adaptive Finite Element Methods Lecture Notes Winter Term 2017/18 R. Verfürth Fakultät für Mathematik, Ruhr-Universität Bochum Contents Chapter I. Introduction 7 I.1. Motivation 7 I.2. Sobolev and finite
More informationANALYSIS OF A MODIFIED SCHRÖDINGER OPERATOR IN 2D: REGULARITY, INDEX, AND FEM
ANALYSIS OF A MODIFIED SCHRÖDINGER OPERATOR IN 2D: REGULARITY, INDEX, AND FEM HENGGUANG LI AND VICTOR NISTOR Version: 2.; Revised: April 2 by V.; Run: January 4, 2008 Abstract. Let r = (x 2 +x2 2 )/2 be
More informationDiscontinuous Galerkin Methods: An Overview and Some Applications. Daya Reddy UNIVERSITY OF CAPE TOWN
Discontinuous Galerkin Methods: An Overview and Some Applications Daya Reddy UNIVRSITY OF CAP TOWN Joint work with Beverley Grieshaber and Andrew McBride SANUM, Stellenbosch, 22 24 March 2016 Structure
More informationRemarks on the analysis of finite element methods on a Shishkin mesh: are Scott-Zhang interpolants applicable?
Remarks on the analysis of finite element methods on a Shishkin mesh: are Scott-Zhang interpolants applicable? Thomas Apel, Hans-G. Roos 22.7.2008 Abstract In the first part of the paper we discuss minimal
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 informationProjected Surface Finite Elements for Elliptic Equations
Available at http://pvamu.edu/aam Appl. Appl. Math. IN: 1932-9466 Vol. 8, Issue 1 (June 2013), pp. 16 33 Applications and Applied Mathematics: An International Journal (AAM) Projected urface Finite Elements
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 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 informationME FINITE ELEMENT ANALYSIS FORMULAS
ME 2353 - FINITE ELEMENT ANALYSIS FORMULAS UNIT I FINITE ELEMENT FORMULATION OF BOUNDARY VALUE PROBLEMS 01. Global Equation for Force Vector, {F} = [K] {u} {F} = Global Force Vector [K] = Global Stiffness
More informationMixed Finite Element Methods. Douglas N. Arnold, University of Minnesota The 41st Woudschoten Conference 5 October 2016
Mixed Finite Element Methods Douglas N. Arnold, University of Minnesota The 41st Woudschoten Conference 5 October 2016 Linear elasticity Given the load f : Ω R n, find the displacement u : Ω R n and the
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 informationAxioms of Adaptivity (AoA) in Lecture 3 (sufficient for optimal convergence rates)
Axioms of Adaptivity (AoA) in Lecture 3 (sufficient for optimal convergence rates) Carsten Carstensen Humboldt-Universität zu Berlin 2018 International Graduate Summer School on Frontiers of Applied and
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 informationMaximum-norm a posteriori estimates for discontinuous Galerkin methods
Maximum-norm a posteriori estimates for discontinuous Galerkin methods Emmanuil Georgoulis Department of Mathematics, University of Leicester, UK Based on joint work with Alan Demlow (Kentucky, USA) DG
More informationSpace-time Finite Element Methods for Parabolic Evolution Problems
Space-time Finite Element Methods for Parabolic Evolution Problems with Variable Coefficients Ulrich Langer, Martin Neumüller, Andreas Schafelner Johannes Kepler University, Linz Doctoral Program Computational
More informationA posteriori error estimates in FEEC for the de Rham complex
A posteriori error estimates in FEEC for the de Rham complex Alan Demlow Texas A&M University joint work with Anil Hirani University of Illinois Urbana-Champaign Partially supported by NSF DMS-1016094
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 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 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 informationSolving Differential Equations on 2-D Geometries with Matlab
Solving Differential Equations on 2-D Geometries with Matlab Joshua Wall Drexel University Philadelphia, PA 19104 (Dated: April 28, 2014) I. INTRODUCTION Here we introduce the reader to solving partial
More informationOn angle conditions in the finite element method. Institute of Mathematics, Academy of Sciences Prague, Czech Republic
On angle conditions in the finite element method Michal Křížek Institute of Mathematics, Academy of Sciences Prague, Czech Republic Joint work with Jan Brandts (University of Amsterdam), Antti Hannukainen
More information1 Discretizing BVP with Finite Element Methods.
1 Discretizing BVP with Finite Element Methods In this section, we will discuss a process for solving boundary value problems numerically, the Finite Element Method (FEM) We note that such method is a
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 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 informationENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED FINITE ELEMENT METHODS
MATHEMATICS OF COMPUTATION Volume 75, Number 256, October 2006, Pages 1659 1674 S 0025-57180601872-2 Article electronically published on June 26, 2006 ENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED
More informationThe Discontinuous Galerkin Finite Element Method
The Discontinuous Galerkin Finite Element Method Michael A. Saum msaum@math.utk.edu Department of Mathematics University of Tennessee, Knoxville The Discontinuous Galerkin Finite Element Method p.1/41
More informationETNA Kent State University
Electronic Transactions on Numerical Analysis. Volume 11, pp. 1-24, 2000. Copyright 2000,. ISSN 1068-9613. ETNA NEUMANN NEUMANN METHODS FOR VECTOR FIELD PROBLEMS ANDREA TOSELLI Abstract. In this paper,
More informationTrefftz-Discontinuous Galerkin Methods for Acoustic Scattering - Recent Advances
Trefftz-Discontinuous Galerkin Methods for Acoustic Scattering - Recent Advances Ilaria Perugia Dipartimento di Matematica Università di Pavia (Italy) In collaboration with Ralf Hiptmair, Christoph Schwab,
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 informationFEM techniques for nonlinear fluids
FEM techniques for nonlinear fluids From non-isothermal, pressure and shear dependent viscosity models to viscoelastic flow A. Ouazzi, H. Damanik, S. Turek Institute of Applied Mathematics, LS III, TU
More informationFEM Convergence for PDEs with Point Sources in 2-D and 3-D
FEM Convergence for PDEs with Point Sources in -D and 3-D Kourosh M. Kalayeh 1, Jonathan S. Graf, and Matthias K. Gobbert 1 Department of Mechanical Engineering, University of Maryland, Baltimore County
More informationNew class of finite element methods: weak Galerkin methods
New class of finite element methods: weak Galerkin methods Xiu Ye University of Arkansas at Little Rock Second order elliptic equation Consider second order elliptic problem: a u = f, in Ω (1) u = 0, on
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 information