Mixed Finite Elements Method
|
|
- Marylou Ramsey
- 6 years ago
- Views:
Transcription
1 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 Poisson equation 2 3 Stokes equation 4 3. First weak formulation Second weak formulation Time Harmonic Maxwell problem 5 5 Time Domain Maxwell problem 6 5. First formulation Second formulation Exercises 7
2 Introduction. Notations Scalar-valued test functions will be denoted by ϕ, while ϕ i will denote a scalar basis function, after reordering all the basis functions of a given discrete space. Vector-valued functions will be written in bold, like u, v. Ψ will denote a vector-valued test function, whileψ i will denote a vector-valued basis function (after reordering the basis functions). Even if most of what follows holds for both the 2D and 3D cases, we will restrict our studies to the 2D one. We recall that in 2D, there are two curl operators, one acting on scalars φ = ( y φ, x φ ) and one acting on vectors Ψ = x Ψ y y Ψ x. Differential operators that return a vector (grad, curl) will be written in bold (, ). We also recall the Green formula for the divergence and curl/rotational operators ( F) G = F G + (F n) G, F H(div, ), G H () () ( G) F = G F (G n) F, F H(curl, ), G H () (2) If R d, u = (u, u 2,..., u d ) and Ψ = (Ψ, Ψ 2,..., Ψ d ), we recall the notation 2 Poisson equation u : Ψ = d i= u i Ψ i (3) Let s start with the Poisson problem { 2 φ = f in, φ = 0 on, (4) As we know, the associated variational formulation can be written as a minimization problem over H (). Let us now, introduce an auxiliary variable u = φ. The Poisson problem is then equivalent to the following system u = f in, u = φ in, φ = 0 on, In order to have a weak formulation, we need to multiply the first equation by a scalar-valued test function φ, while the second one will be multiplied by a vector-valued test function Ψ. Therefor, ( u) ϕ = f ϕ, u Ψ = (6) φ Ψ, (5) 2
3 Now let us use the Green s formula for the divergence. Note that the integration by parts will only be done on the second equation. This leads to the following weak formulation ( u) ϕ + f ϕ = 0, u Ψ + φ Ψ = 0, (7) Now, let us consider the bilinear form B (( ) ( )) ( ) u Ψ, = ( u) ϕ + u Ψ + φ Ψ φ ϕ ( ) u Where we consider the trial vector-valued function and the test vector-valued function φ For the sake of clarity, we rewrite the blinear form as the following B (( ) ( )) ( ) u Ψ, = u Ψ + ( u ) ϕ + φ Ψ φ ϕ (8) ( ) Ψ. ϕ Therefor, it is easy to see that our blinear form is symmetric. Up to now, we didn t mention what are the correct spaces for this formulation. Because, there is no differential operator over the unknown φ and the scalar-valued test function ϕ, these functions can be placed in L 2 (). On the other hand, both u, Ψ, u and Ψ must be in L 2 (). Therefor, we need the Hilbert space H(div, ) defined as (9) H(div, ) = { Ψ L 2 (), Ψ L 2 () } (0) Therefor, the mixed weak formulation for the Poisson equation is Find (u, φ) H(div, ) L 2 () such that (Ψ, ϕ) H(div, ) L 2 () u Ψ + φ Ψ = 0, ( u) ϕ = f ϕ, () Now let us introduce the bilinear form a(u, Ψ) = u Ψ, u, Ψ H(div, ) and b(φ, Ψ) = φ Ψ, (φ, Ψ) L2 () H(div, ). The later problem can be written as Find (u, φ) H(div, ) L 2 () such that (Ψ, ϕ) H(div, ) L 2 () a(u, Ψ) + b(φ, Ψ) = r (Ψ), b(ϕ, u) = r 2 (ϕ), (2) 3
4 3 Stokes equation Let u be the velocity of a fluid and p the pressure in the fluid. We denote by f the body forces. The Stokes problem writes 2 u + p = f in, u = g in, u = 0 on, When g = 0, the fluid is said to be incompressible. Because u = u n = 0, the source term g must fullfill the compatitbility condition g = 0. Therefor, the function g will be in the space L0 2 (), defined as } L0 {ϕ 2 () = L 2 (); ϕ = 0 (4) Solving the Stokes equation consists of looking for u H () and p L 2 0 (), where f H () and g L 2 0 (). 3. First weak formulation The first formulation consists on putting the constraint u = 0 directly into the continuous (and discrete) space. We start by introducing the space Find u V such that Ψ H 0 (div, ) (3) H 0 (div, ) = { Ψ H 0 (), Ψ = 0} (5) u : Ψ = f Ψ (6) This bilinear form can be shown to be coercive and one can directly apply the Lax-Milgram theorem. However, this formulation has a major drawback from the practical point of view: the pressure is not solved directly. Moreover, on the discrete level, we need to construct an increasing sequence (V h ) h of subspaces of V. This can be done thanks to the discrete DeRham sequence. 3.2 Second weak formulation Now let us keep the space for u as big as possible (which means the space H0 ()) and then add the constraint u = 0, as a Lagrange multiplier. In this case, the problem writes Find (u, p) H0 () L2 0 () such that (Ψ, ϕ) H 0 () L2 0 () u : Ψ p Ψ = f Ψ (7) ϕ u = 0 4
5 Now let us introduce the bilinear form a(u, Ψ) = u : Ψ, u, Ψ H(div, ) and b(p, Ψ) = p Ψ, (p, Ψ) L2 () H(div, ). The later problem can be written as Find (u, p) H 0 () L2 0 () such that (Ψ, ϕ) H 0 () L2 0 () a(u, Ψ) + b(p, Ψ) = r (Ψ), b(ϕ, u) = r 2 (ϕ), (8) 4 Time Harmonic Maxwell problem The Time Harmonic Maxwell problem writes Find E 0 such that { E = ω 2 E in, E n = 0 on, (9) Two cases occure: either ω is provided and then we only look for E. The other case is known as an eigenvalue problem, and the problem must be solved for both E and ω. If we multiply by a vector-valued function Ψ, and integrating by parts using the curl/rotational Green formula, we get E Ψ = ω 2 E Ψ (20) Again, the space where E lives was not mentioned; let us find what this space can be. We need both E, Ψ, E and Ψ to be in L 2 (), while satisfying the boundary condition E n = 0 on the boundary. Therefor, we need the Hilbert space H 0 (curl, ) defined as where H 0 (curl, ) = {Ψ H(curl, ), Ψ n = 0 on } (2) H(curl, ) = { Ψ L 2 (), Ψ L 2 () } (22) Therefor, in weak form, the problem writes Find E H 0 (curl, ), E 0 such that, Ψ H 0 (curl, ) E Ψ = ω 2 E Ψ (23) 5
6 5 Time Domain Maxwell problem The Maxwell equations write D = ρ B = 0 t B = E t D = H J In the case of linear, isotropic and non-dispersive materials, we have two additional relations: B = µh and D = ɛe. Furthermore, we assume that J = σe. Pluging these relations in (Eq. 24), we get (ɛe) = ρ (µh) = 0 (25) µ t H = E ɛ t E = H J Now, let us assume that both ɛ and µ are constants and equal to. We also, restrict our study to the 2D case and we consider the Transverse Electric (TE) mode. Therefor, the problem is only involving the variables E x, E y and H z. In this case, the Time Domain Maxwell problem writes t H = E t E = H J (26) E = ρ where we consider here that E = (E x, E y ) and H = H z. 5. First formulation The first variational formulation, in the case of perfectly conductin boundary conditions, can be derived using the curl/rotational Green formula in Ampere s law. Therefor, the TE Maxwell problem reads Find (E, H) H 0 (curl, ) L 2 () such that (Ψ, ϕ) H 0 (curl, ) L 2 () d dt E Ψ H Ψ = J Ψ d dt Hϕ + (27) ( E)ϕ = 0, 5.2 Second formulation Using the divergence Green formula in Faraday s law yields the second variational formulation. Find (E, H) H(div, ) H () such that (Ψ, ϕ) H(div, ) H () d dt E Ψ ( H) Ψ = J Ψ d dt Hϕ + (28) E ( ϕ) = 0 (24) 6
7 6 Exercises Exercise 6.. Let us consider the mixed formulation (Eq. 29) of the 2D Poisson problem. Let V h H(div, ) and W h L 2 () be finite dimensional subspaces. The Galerkin method applied to (Eq. 29) leads to the following discrete problem Find (u h, φ h ) V h W h such that (Ψ, ϕ) V h W h u h Ψ + φ h Ψ = 0, We assume that the computational domain is the unit square. ( u h) ϕ = f ϕ, (29) Figure : The computational domain and its triangulation.. Let us assume that V h = span {Ψ, Ψ 2,..., Ψ n } and W h = span {ϕ, ϕ 2,..., ϕ m }. Write the global matrices and the corresponding right hand sides. 2. The latter linear system can be written as ( A B T B 0 ) ( ) u = φ ( ) f g (30) (a) We assume that A is invertible. Show that the system 30 has a unique solution if the matrix BA B T is invertible. (This method is called the Schur Complement) (b) If A is positive definite, show that BA B T is also positive definite. (c) Use the symmetry of A to show that BA B T is also a symmetric matrix. 7
8 3. We consider a triangulation T h constructed as the following: we take a uniform mesh of n n subsquares, where every subsquare is devided into two triangles as in the figure (Fig. ). We consider V h to be the space of piecewise linear vector fields and W h the space of continuous linear scalar fields, (a) For T T h, a given triangle, construct a function v W h, such that T v = 0. Show how to construct a piecewise linear function that is orthogonal to all piecewise constants. (b) Show that v ϕ = 0, ϕ V h, (c) Show that the stiffness matrix is singular. (Hint: Show that the function (0, v) belongs to the kernel of the stiffness matrix), Exercise 6.2 (The Checkerboard instability). In this exercise, we study the Stokes equation with the second weak formulation, when the inf-sup condition is not fullfilled. The compuational domain will be the unit square, with a uniform mesh of n n subsquares Q h = n i, j= Q i, j, as described in figure (Fig. 2). Let us recall that the variational problem can be written as Figure 2: The computational domain and the its mesh. Find (u, p) H 0 () L2 0 () such that (Ψ, ϕ) H 0 () L2 0 () u : Ψ p Ψ = f Ψ ϕ u = 0 (3) We consider V h to be the space of continuous linear vector fields and W h the space of piecewise 8
9 constants scalar fields, i.e. V h = { u h [C 0 ( )] 2 ; T T h, u h [P ] 2, u h = 0 } { } W h = φ h L 2 () C 0 ( ); T T h, φ h P, φ h = 0 For better clarity, we define the function ϕ W h, which is a piecewise constant, by its value in the center of our meshes, i.e. ϕ := ϕ i+ 2, j+, while for any function Ψ V h is defined by its 2 values on the vertices of our mesh, i.e. Ψ i, j = ( Ψi, x j, Ψy i, j). The aim of this exercise is to construct a non zero scalar-valued function q such that q Ψ = 0, Ψ V h.. Show that for all (Ψ, ϕ) V h W h, we have ϕ Ψ = ( Ψ x ( n 2 i, j δ x [ϕ] ) i, j + ( Ψy i, j δ y [ϕ] ) i, j) 2 i, j n (32) where the operators δ x and δ y are defined as ( δ x [ϕ] ) i, j := n { } ϕ 2 i+ 2, j+ ϕ 2 i 2, j+ + ϕ 2 i+ 2, j ϕ 2 i 2, j 2 ( δ y [ϕ] ) i, j := n { } ϕ 2 i+ 2, j+ + ϕ 2 i 2, j+ ϕ 2 i+ 2, j ϕ 2 i 2, j 2 and Ψ = (Ψ x, Ψ y ) 2. Prove that ϕ Ψ = 0, Ψ V h if and only if for any i, j n we have ϕ i+ 2, j+ 2 ϕ i 2, j+ 2 = ϕ i 2, j 2 = ϕ i+ 2, j 2 3. Construct such a function by using the values and +. Such oscillatory function is known as a spurious mode. 9
Mixed 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 informationDiscontinuous Petrov-Galerkin Methods
Discontinuous Petrov-Galerkin Methods Friederike Hellwig 1st CENTRAL School on Analysis and Numerics for Partial Differential Equations, November 12, 2015 Motivation discontinuous Petrov-Galerkin (dpg)
More informationA primer on Numerical methods for elasticity
A primer on Numerical methods for elasticity Douglas N. Arnold, University of Minnesota Complex materials: Mathematical models and numerical methods Oslo, June 10 12, 2015 One has to resort to the indignity
More informationA Finite Element Method for the Surface Stokes Problem
J A N U A R Y 2 0 1 8 P R E P R I N T 4 7 5 A Finite Element Method for the Surface Stokes Problem Maxim A. Olshanskii *, Annalisa Quaini, Arnold Reusken and Vladimir Yushutin Institut für Geometrie und
More informationA NEW APPROXIMATION TECHNIQUE FOR DIV-CURL SYSTEMS
MATHEMATICS OF COMPUTATION Volume 73, Number 248, Pages 1739 1762 S 0025-5718(03)01616-8 Article electronically published on August 26, 2003 A NEW APPROXIMATION TECHNIQUE FOR DIV-CURL SYSTEMS JAMES H.
More informationFinite Element Exterior Calculus and Applications Part V
Finite Element Exterior Calculus and Applications Part V Douglas N. Arnold, University of Minnesota Peking University/BICMR August 15 18, 2015 De Rham complex 0 H 1 (Ω) grad H(curl, Ω) curl H(div, Ω) div
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 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 informationSpline Element Method for Partial Differential Equations
for Partial Differential Equations Department of Mathematical Sciences Northern Illinois University 2009 Multivariate Splines Summer School, Summer 2009 Outline 1 Why multivariate splines for PDEs? Motivation
More informationRadiation Integrals and Auxiliary Potential Functions
Radiation Integrals and Auxiliary Potential Functions Ranga Rodrigo June 23, 2010 Lecture notes are fully based on Balanis [?]. Some diagrams and text are directly from the books. Contents 1 The Vector
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 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 informationWeierstraß-Institut. für Angewandte Analysis und Stochastik. Leibniz-Institut im Forschungsverbund Berlin e. V. Preprint ISSN
Weierstraß-Institut für Angewandte Analysis und Stochastik Leibniz-Institut im Forschungsverbund Berlin e. V. Preprint ISSN 2198-5855 On the divergence constraint in mixed finite element methods for incompressible
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 informationGlowinski Pironneau method for the 3D ω-ψ equations
280 GUERMOND AND QUARTAPELLE Glowinski Pironneau method for the 3D ω-ψ equations Jean-Luc Guermond and Luigi Quartapelle 1 LIMSI CNRS, Orsay, France, and Dipartimento di Fisica, Politecnico di Milano,
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 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 informationA Generalization for Stable Mixed Finite Elements
A Generalization for Stable Mixed Finite Elements Andrew Gillette joint work with Chandrajit Bajaj Department of Mathematics Institute of Computational Engineering and Sciences University of Texas at Austin,
More informationApproximation in Banach Spaces by Galerkin Methods
2 Approximation in Banach Spaces by Galerkin Methods In this chapter, we consider an abstract linear problem which serves as a generic model for engineering applications. Our first goal is to specify the
More informationA SADDLE POINT LEAST SQUARES METHOD FOR SYSTEMS OF LINEAR PDES. Klajdi Qirko
A SADDLE POINT LEAST SQUARES METHOD FOR SYSTEMS OF LINEAR PDES by Klajdi Qirko A dissertation submitted to the Faculty of the University of Delaware in partial fulfillment of the requirements for the degree
More informationWeighted 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 informationINF-SUP CONDITION FOR OPERATOR EQUATIONS
INF-SUP CONDITION FOR OPERATOR EQUATIONS LONG CHEN We study the well-posedness of the operator equation (1) T u = f. where T is a linear and bounded operator between two linear vector spaces. We give equivalent
More information7.4 The Saddle Point Stokes Problem
346 CHAPTER 7. APPLIED FOURIER ANALYSIS 7.4 The Saddle Point Stokes Problem So far the matrix C has been diagonal no trouble to invert. This section jumps to a fluid flow problem that is still linear (simpler
More informationChapter 12. Partial di erential equations Di erential operators in R n. The gradient and Jacobian. Divergence and rotation
Chapter 12 Partial di erential equations 12.1 Di erential operators in R n The gradient and Jacobian We recall the definition of the gradient of a scalar function f : R n! R, as @f grad f = rf =,..., @f
More informationNitsche-type Mortaring for Maxwell s Equations
Nitsche-type Mortaring for Maxwell s Equations Joachim Schöberl Karl Hollaus, Daniel Feldengut Computational Mathematics in Engineering at the Institute for Analysis and Scientific Computing Vienna University
More informationPhysic-based Preconditioning and B-Splines finite elements method for Tokamak MHD
Physic-based Preconditioning and B-Splines finite elements method for Tokamak MHD E. Franck 1, M. Gaja 2, M. Mazza 2, A. Ratnani 2, S. Serra Capizzano 3, E. Sonnendrücker 2 ECCOMAS Congress 2016, 5-10
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 informationINNOVATIVE FINITE ELEMENT METHODS FOR PLATES* DOUGLAS N. ARNOLD
INNOVATIVE FINITE ELEMENT METHODS FOR PLATES* DOUGLAS N. ARNOLD Abstract. Finite element methods for the Reissner Mindlin plate theory are discussed. Methods in which both the tranverse displacement and
More informationSchur Complements on Hilbert Spaces and Saddle Point Systems
Schur Complements on Hilbert Spaces and Saddle Point Systems Constantin Bacuta Mathematical Sciences, University of Delaware, 5 Ewing Hall 976 Abstract For any continuous bilinear form defined on a pair
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 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 informationMaxwell s Equations. In the previous chapters we saw the four fundamental equations governging electrostatics and magnetostatics. They are.
Maxwell s Equations Introduction In the previous chapters we saw the four fundamental equations governging electrostatics and magnetostatics. They are D = ρ () E = 0 (2) B = 0 (3) H = J (4) In the integral
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 informationNumerical Analysis of Nonlinear Multiharmonic Eddy Current Problems
Numerical Analysis of Nonlinear Multiharmonic Eddy Current Problems F. Bachinger U. Langer J. Schöberl April 2004 Abstract This work provides a complete analysis of eddy current problems, ranging from
More informationFinite Element Exterior Calculus. Douglas N. Arnold, University of Minnesota The 41st Woudschoten Conference 5 7 October 2016
Finite Element Exterior Calculus Douglas N. Arnold, University of Minnesota The 41st Woudschoten Conference 5 7 October 2016 The Hilbert complex framework Discretization of Hilbert complexes Finite element
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 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 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 informationAn inverse problem for eddy current equations
An inverse problem for eddy current equations Alberto Valli Department of Mathematics, University of Trento, Italy In collaboration with: Ana Alonso Rodríguez, University of Trento, Italy Jessika Camaño,
More informationLECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES. Sergey Korotov,
LECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES Sergey Korotov, Institute of Mathematics Helsinki University of Technology, Finland Academy of Finland 1 Main Problem in Mathematical
More informationThe Dirichlet-to-Neumann operator
Lecture 8 The Dirichlet-to-Neumann operator The Dirichlet-to-Neumann operator plays an important role in the theory of inverse problems. In fact, from measurements of electrical currents at the surface
More information2.3 Variational form of boundary value problems
2.3. VARIATIONAL FORM OF BOUNDARY VALUE PROBLEMS 21 2.3 Variational form of boundary value problems Let X be a separable Hilbert space with an inner product (, ) and norm. We identify X with its dual X.
More informationFinite Element methods for hyperbolic systems
Finite Element methods for hyperbolic systems Eric Sonnendrücker Max-Planck-Institut für Plasmaphysik und Zentrum Mathematik, TU München Lecture notes Wintersemester 14-15 January 19, 15 Contents 1 Introduction
More informationChapter 2 Finite Element Spaces for Linear Saddle Point Problems
Chapter 2 Finite Element Spaces for Linear Saddle Point Problems Remark 2.1. Motivation. This chapter deals with the first difficulty inherent to the incompressible Navier Stokes equations, see Remark
More informationThe Divergence Theorem Stokes Theorem Applications of Vector Calculus. Calculus. Vector Calculus (III)
Calculus Vector Calculus (III) Outline 1 The Divergence Theorem 2 Stokes Theorem 3 Applications of Vector Calculus The Divergence Theorem (I) Recall that at the end of section 12.5, we had rewritten Green
More informationYongdeok Kim and Seki Kim
J. Korean Math. Soc. 39 (00), No. 3, pp. 363 376 STABLE LOW ORDER NONCONFORMING QUADRILATERAL FINITE ELEMENTS FOR THE STOKES PROBLEM Yongdeok Kim and Seki Kim Abstract. Stability result is obtained for
More informationIntroduction. J.M. Burgers Center Graduate Course CFD I January Least-Squares Spectral Element Methods
Introduction In this workshop we will introduce you to the least-squares spectral element method. As you can see from the lecture notes, this method is a combination of the weak formulation derived from
More informationChapter 2. General concepts. 2.1 The Navier-Stokes equations
Chapter 2 General concepts 2.1 The Navier-Stokes equations The Navier-Stokes equations model the fluid mechanics. This set of differential equations describes the motion of a fluid. In the present work
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 informationApproximation of fluid-structure interaction problems with Lagrange multiplier
Approximation of fluid-structure interaction problems with Lagrange multiplier Daniele Boffi Dipartimento di Matematica F. Casorati, Università di Pavia http://www-dimat.unipv.it/boffi May 30, 2016 Outline
More informationWELL POSEDNESS OF PROBLEMS I
Finite Element Method 85 WELL POSEDNESS OF PROBLEMS I Consider the following generic problem Lu = f, where L : X Y, u X, f Y and X, Y are two Banach spaces We say that the above problem is well-posed (according
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 informationDiscontinuous Galerkin Methods
Discontinuous Galerkin Methods Joachim Schöberl May 20, 206 Discontinuous Galerkin (DG) methods approximate the solution with piecewise functions (polynomials), which are discontinuous across element interfaces.
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 informationChapter 1 Foundations of Elliptic Boundary Value Problems 1.1 Euler equations of variational problems
Chapter 1 Foundations of Elliptic Boundary Value Problems 1.1 Euler equations of variational problems Elliptic boundary value problems often occur as the Euler equations of variational problems the latter
More informationFinite Element Methods for Optical Device Design
DFG Research Center Matheon mathematics for key technologies and Zuse Institute Berlin Finite Element Methods for Optical Device Design Frank Schmidt Sven Burger, Roland Klose, Achim Schädle, Lin Zschiedrich
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 informationAn introduction to the mathematical theory of finite elements
Master in Seismic Engineering E.T.S.I. Industriales (U.P.M.) Discretization Methods in Engineering An introduction to the mathematical theory of finite elements Ignacio Romero ignacio.romero@upm.es October
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 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 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 informationPhysics 6303 Lecture 3 August 27, 2018
Physics 6303 Lecture 3 August 27, 208 LAST TIME: Vector operators, divergence, curl, examples of line integrals and surface integrals, divergence theorem, Stokes theorem, index notation, Kronecker delta,
More informationDot Products. K. Behrend. April 3, Abstract A short review of some basic facts on the dot product. Projections. The spectral theorem.
Dot Products K. Behrend April 3, 008 Abstract A short review of some basic facts on the dot product. Projections. The spectral theorem. Contents The dot product 3. Length of a vector........................
More informationAMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends
AMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends Lecture 3: Finite Elements in 2-D Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Finite Element Methods 1 / 18 Outline 1 Boundary
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 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 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 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 informationConvergence of the MAC scheme for incompressible flows
Convergence of the MAC scheme for incompressible flows T. Galloue t?, R. Herbin?, J.-C. Latche??, K. Mallem????????? Aix-Marseille Universite I.R.S.N. Cadarache Universite d Annaba Calanque de Cortiou
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 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 informationRobust Domain Decomposition Preconditioners for Abstract Symmetric Positive Definite Bilinear Forms
www.oeaw.ac.at Robust Domain Decomposition Preconditioners for Abstract Symmetric Positive Definite Bilinear Forms Y. Efendiev, J. Galvis, R. Lazarov, J. Willems RICAM-Report 2011-05 www.ricam.oeaw.ac.at
More informationON LEAST-SQUARES FINITE ELEMENT METHODS FOR THE POISSON EQUATION AND THEIR CONNECTION TO THE DIRICHLET AND KELVIN PRINCIPLES
ON LEAST-SQUARES FINITE ELEMENT METHODS FOR THE POISSON EQUATION AND THEIR CONNECTION TO THE DIRICHLET AND KELVIN PRINCIPLES PAVEL BOCHEV 1,3 AND MAX GUNZBURGER 2 Abstract. Least-squares finite element
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 informationSolutions of Selected Problems
1 Solutions of Selected Problems October 16, 2015 Chapter I 1.9 Consider the potential equation in the disk := {(x, y) R 2 ; x 2 +y 2 < 1}, with the boundary condition u(x) = g(x) r for x on the derivative
More informationFast Iterative Solution of Saddle Point Problems
Michele Benzi Department of Mathematics and Computer Science Emory University Atlanta, GA Acknowledgments NSF (Computational Mathematics) Maxim Olshanskii (Mech-Math, Moscow State U.) Zhen Wang (PhD student,
More informationInverse Eddy Current Problems
Inverse Eddy Current Problems Bastian von Harrach bastian.harrach@uni-wuerzburg.de (joint work with Lilian Arnold) Institut für Mathematik - IX, Universität Würzburg Oberwolfach workshop on Inverse Problems
More informationCONVERGENCE THEORY. G. ALLAIRE CMAP, Ecole Polytechnique. 1. Maximum principle. 2. Oscillating test function. 3. Two-scale convergence
1 CONVERGENCE THEOR G. ALLAIRE CMAP, Ecole Polytechnique 1. Maximum principle 2. Oscillating test function 3. Two-scale convergence 4. Application to homogenization 5. General theory H-convergence) 6.
More informationKonstantinos Chrysafinos 1 and L. Steven Hou Introduction
Mathematical Modelling and Numerical Analysis Modélisation Mathématique et Analyse Numérique Will be set by the publisher ANALYSIS AND APPROXIMATIONS OF THE EVOLUTIONARY STOKES EQUATIONS WITH INHOMOGENEOUS
More informationLinear Algebra Massoud Malek
CSUEB Linear Algebra Massoud Malek Inner Product and Normed Space In all that follows, the n n identity matrix is denoted by I n, the n n zero matrix by Z n, and the zero vector by θ n An inner product
More informationMathematical Notes for E&M Gradient, Divergence, and Curl
Mathematical Notes for E&M Gradient, Divergence, and Curl In these notes I explain the differential operators gradient, divergence, and curl (also known as rotor), the relations between them, the integral
More informationFEM for Stokes Equations
FEM for Stokes Equations Kanglin Chen 15. Dezember 2009 Outline 1 Saddle point problem 2 Mixed FEM for Stokes equations 3 Numerical Results 2 / 21 Stokes equations Given (f, g). Find (u, p) s.t. u + p
More informationMIXED FINITE ELEMENT APPROXIMATION OF THE VECTOR LAPLACIAN WITH DIRICHLET BOUNDARY CONDITIONS
MIXED FINITE ELEMENT APPROXIMATION OF THE VECTOR LAPLACIAN WITH DIRICHLET BOUNDARY CONDITIONS DOUGLAS N. ARNOLD, RICHARD S. FALK, AND JAY GOPALAKRISHNAN Abstract. We consider the finite element solution
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 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 informationKevin James. MTHSC 3110 Section 2.2 Inverses of Matrices
MTHSC 3110 Section 2.2 Inverses of Matrices Definition Suppose that T : R n R m is linear. We will say that T is invertible if for every b R m there is exactly one x R n so that T ( x) = b. Note If T is
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Sobolev Embedding Theorems Embedding Operators and the Sobolev Embedding Theorem
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 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 informationSpectral Properties of Elliptic Operators In previous work we have replaced the strong version of an elliptic boundary value problem
Spectral Properties of Elliptic Operators In previous work we have replaced the strong version of an elliptic boundary value problem L u x f x BC u x g x with the weak problem find u V such that B u,v
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 informationReductions of Operator Pencils
Reductions of Operator Pencils Olivier Verdier Department of Mathematical Sciences, NTNU, 7491 Trondheim, Norway arxiv:1208.3154v2 [math.na] 23 Feb 2014 2018-10-30 We study problems associated with an
More informationON THE EXISTENCE OF TRANSMISSION EIGENVALUES. Andreas Kirsch1
Manuscript submitted to AIMS Journals Volume 3, Number 2, May 29 Website: http://aimsciences.org pp. 1 XX ON THE EXISTENCE OF TRANSMISSION EIGENVALUES Andreas Kirsch1 University of Karlsruhe epartment
More informationarxiv: v3 [math.na] 8 Sep 2015
A Recovery-Based A Posteriori Error Estimator for H(curl) Interface Problems arxiv:504.00898v3 [math.na] 8 Sep 205 Zhiqiang Cai Shuhao Cao Abstract This paper introduces a new recovery-based a posteriori
More informationMaxwell s equations in Carnot groups
Maxwell s equations in Carnot groups B. Franchi (U. Bologna) INDAM Meeting on Geometric Control and sub-riemannian Geometry Cortona, May 21-25, 2012 in honor of Andrey Agrachev s 60th birthday Researches
More informationScientific Computing I
Scientific Computing I Module 8: An Introduction to Finite Element Methods Tobias Neckel Winter 2013/2014 Module 8: An Introduction to Finite Element Methods, Winter 2013/2014 1 Part I: Introduction to
More informationProblem of Second grade fluids in convex polyhedrons
Problem of Second grade fluids in convex polyhedrons J. M. Bernard* Abstract This article studies the solutions of a three-dimensional grade-two fluid model with a tangential boundary condition, in a polyhedron.
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 informationAnalysis of Hybrid Discontinuous Galerkin Methods for Incompressible Flow Problems
Analysis of Hybrid Discontinuous Galerkin Methods for Incompressible Flow Problems Christian Waluga 1 advised by Prof. Herbert Egger 2 Prof. Wolfgang Dahmen 3 1 Aachen Institute for Advanced Study in Computational
More informationNew Model Stability Criteria for Mixed Finite Elements
New Model Stability Criteria for Mixed Finite Elements Andrew Gillette Department of Mathematics Institute of Computational Engineering and Sciences University of Texas at Austin http://www.math.utexas.edu/users/agillette
More information