Motivations. Outline. Finite element exterior calculus and the geometrical basis of numerical stability. References. Douglas N.
|
|
- Peregrine Norman
- 6 years ago
- Views:
Transcription
1 Outline Finite element exterior calculus and the geometrical basis of numerical stability Douglas N. Arnold joint with R. Falk and R. Winther School of Mathematics University of Minnesota April 009 Motivations Exterior calculus and PDE 3 Discretization 4 Finite element differential forms 5 Application to elasticity / 40 / 40 References Finite element exterior calculus, homological techniques, and applications, Acta Numerica 006, p. 55 Any young (or not so young) mathematician who spends the time to master this paper will have tools that will be useful for his or her entire career. Math Reviews Mixed finite element methods for linear elasticity with weakly imposed symmetry, Math. Comp. 007 Motivations Why do we need more theory for finite elements? Geometric decompositions and local bases for spaces finite element differential forms, to appear in CMAME everything is at arnold 3 / 40
2 Why do we need FEEC? The finite element method is incredibly successful. FEM often amenable to mathematical analysis, allowing validation and comparison of methods. But plenty of challenges remain, for algorithms and analysis! Approximability, consistency, and stability = convergence Stability, like well-posedness, can be extremely subtle Well-posedness + approximability + consistency stability Exterior calculus, de Rham cohomology, Hodge theory,..., are geometric tools to get at well-posedness. FEEC adapts these tools to the discrete level to get at stability. 5 / 40 Steady heat conduction problem: finite elements in H div C grad u = f strong C grad u grad v dx = f v dx v weak ( ) C grad u grad u f u u dx minimum variational grad u dx < u H () H : u L (), grad u L (; R n ) The right FE space: Lagrange elements { v H () : v T Pr (T ) T Th } 6 / 40 Lagrange finite elements First order (mixed) formulation Shape fns: P DOFs: vertex vals. & edge averages A σ = grad u, divσ = f strong A σ τ dx = div τ u dx τ, div σ v dx = f v dx v ( A σ σ + div σ u + f u) dx σ, u stationary pt. weak variational P P P3 σ H(div, ), u L () Lagrange elements? Unstable! 7 / 40 8 / 40
3 Thermal problem in D σ = u, σ = f on (, ) (σ + σ σ, u u + f u) dx stationary point H L Babuška Narasimhan Thermal problem in D σ = grad u, div σ = f ( σ σ, u + div σ u + f u) dx stationary point H(div) L P P0 P Λ P0 P-P (0 elts) P-P (40 elts) P-P0 (40 elts) 9 / 40 0 / 40 Raviart Thomas elements (P r Λ ) Shape functions P Λ : span [ ( 0), ( 0 ), ( x y) ] DOFs: A mixed FEM for nd order elliptic problems, Proc. conf. Math l Aspects of the FEM, Rome 975. Springer Lect. Notes in Math #606, 977. Generalizes to all degrees, and all dimensions (n = 3: Nédélec 80) Math & CS SIAM J. Numerical Analysis Numerische Mathematik Mathematics of Computation RAIRO M AN Num. Methods for PDEs Eng. & Apps CMAME Computational Geosciences J. Computational Physics IJNME COMPEL Maxwell eigenvalue problem, unstructured mesh µ curl E curl Ẽ = ω ɛe Ẽ Ẽ Right space is H(curl) λ = m + n = 0,,,, 4, 4, 5, 5, 8, (Lag.P) P Λ / 40 / 40
4 Maxwell eigenvalue problem, regular mesh Vector Laplacian curl curl u grad div u = f in Ω λ = m + n =,,, 4, 4, 5, 5, 8,... u n = 0, Z Ω rot u = 0 on Ω u ( curl u + div u ) f u minimum Lagrange finite elements converge nicely but not to the solution! (same problem with any conforming FE) A mixed formulation based on appropriate finite elements works fine Z Boffi-Brezzi-Gastaldi 99 3 / 40 Ω σ,u ( σ curl σ u div u f u) dx stationary point H(curl) H(div) 4 / 40 EM calculations based on the generalized RT elements Scho berl, Zaglmayr 006, NGSolve Exterior calculus and PDE The continuous problem K tets, P6 Λ Also: White EMSolve, Demkowicz 3Dhp90, Durufle Montjoie,... 6K tets, P3 Λ 5 / 40
5 Differential forms For R n, Λ k () consists of functions Alt k R n so if ω Λ k (), ωx(v,..., vk) R, x, vi R n Λ 0 (): real-valued functions on Λ (): covector fields, ω = n i= fidxi, fi functions (dxi(ej) = δij) Λ (): ω = i<j fijdxi dxj, (dxi dxj := dxi dxj dxj dxi) De Rham complex and cohomology 0 Λ 0 () d0 Λ () d dn Λ n () 0 # of components, i = 0, Zk := ker(d k ) # of holes, i =, Bk := range(d k dim Zk/Bk = ) # of voids, i =, vector proxies in R 3 : fidxi (f, f, f3), fijdxidxj (f3, f3, f) 0 C () grad C (; R 3 ) curl C (; R 3 ) div C () 0 d k : Λ k () Λ k+ () d(f dxj dxk) = i ω R ω Λ k, dim f = k f f xi dxi dxj dxk Physical vector quantities may be divided into two classes, in one of which the quantity is dened with reference to a line, while in the other the quantity is dened with reference to an area. James Clerk Maxwell, Treatise on Electricity & Magnetism, 89 7 / 40 8 / 40 Hodge theory Making use of the inner product: Hodge star: : Λ k () = Λ n k () formal adjoint: δ = ± d : Λ k () Λ k () Λ k () d δ Λ k () d Λ k+ () δ Hodge Laplacian: dδ + δd : Λ k Λ k harmonic forms: H k := { ζ Z k ζ B k } = Z k /B k Hodge decomposition: L Λ k () = B k H k (Z k ) Poincaré s inequality: ω L c dω L, ω (Z k ) Sobolev spaces: HΛ k () = { ω L Λ k () dω L Λ k+ () } 0 HΛ 0 () 0 H () d HΛ () grad H(curl, ) d curl H(div, ) d HΛ n () 0 div L () 0 Some applications Physical quantities: 0-forms: temperature; electric field potential -forms: temperature gradient; electric field -forms: heat flux; magnetic flux stress is a covector-valued -form 3-forms: heat density; charge density; mass density PDEs: div grad u = f (curl curl grad div)u = f curl curl u = f, div u = 0 div u = f, curl u = 0 Maxwell s equations elasticity dynamic problems, eigenvalue problems, lower order-terms variable coefficients, nonlinearities... 9 / 40 0 / 40
6 Hodge Laplace problem Well-posedness of the Hodge Laplacian Given f Λ k (0 k n), find u Λ k with (dδ + δd)u = f (plus BC) Harmonic functions determine well-posedness: u f H k, u is determined only mod H k This mixed formulation is always well-posed: Given f L Λ k (), find σ HΛ k, u HΛ k, p H k : σ, τ dτ, u = 0 τ HΛ k dσ, v + du, dv + p, v = f, v v HΛ k u, q = 0 q H k Equivalently σ, σ du, du dσ, u u, p + f, u saddle point / 40 σ HΛ k, u HΛ k, p H k : σ, τ dτ, u = 0 dσ, v + du, dv + p, v = f, v u, q = 0 τ HΛ k v HΛ k q H k Need to control σ HΛ + u HΛ + p by a bounded choice of τ, v, and q. τ = σ controls σ, v = dσ controls dσ, v = p controls p v = u controls du, How to control u?? Hodge decomp.: u = dη + s + z, η HΛ k, s H k, z (Z k ) τ = η controls dη and q = s controls s. To bound z we use Poincaré s inequality: z c dz = c du (which is under control) / 40 Abstract setting Λ k d k Λ k Λ k+ d k Discretization We have a well-posed variational PDE problem. How do we discretize it stably? Λ k h d k Λ k h d k Λ k+ h Complex of Hilbert spaces with d k bounded and closed range. For discretization, construct a finite dimensional subcomplex. Define H k h = (Bk h ) Z k h. Discrete Hodge decomp. follows: Λ k h = Bk h Hk h (Zk h ) Galerkin s method: Λ k, Λ k, H k Λ k h, Λ k h, Hk h When is it stable? 4 / 40
7 Bounded cochain projections Key property: The finite dimensional subcomplex admits a bounded cochain projection. Λ k d k Λ k π k h Λ k h π k h d k Λ k h π k h bounded π k h a projection πh kd k = d k π k h limh 0 πh k v = v, v Λk Theorem The induced map on cohomology is an isomorphism for h small. gap ( H k, H k h) 0 The discrete Poincaré inequality holds uniformly in h. Galerkin s method is stable and convergent. Proof of discrete Poincaré inequality Thm. There is a positive constant c, independent of h, such that ω c dω, ω Z k h. Proof. Given ω Z k h, define η Zk HΛ k () by dη = dω. By the Poincaré inequality, η c dω, so it is enough to show that ω c η. Now, ω πhη Λ k h and d(ω πhη) = 0, so ω πhη Z k h. Therefore ω = ω, πhη + ω, ω πhη = ω, πhη ω πhη, whence ω πhη, and the result follows. 5 / 40 6 / 40 Finite element differential forms Finite element differential forms How do we construct finite element spaces that fit together in de Rham subcomplexes with bounded cochain projections? Let T = Th be a triangulation of R n. We wish to construct finite element spaces Λ k (T ) HΛ k () which form a finite dimensional subcomplex with bounded cochain projections. We will construct them as usual for finite elements: On each simplex T T we specify a space of polynomials shape functions degrees of freedom, each associated to a face of the simplex It turns out that for each form degree k and polynomial degree r, there are just two natural finite element subspaces of HΛ k (): Pr Λ k (Th) and P r Λ k (Th) Pr Λ k (Th), k = 0, Pr Λ k (Th) = Pr Λ k (Th), k = n, strictly between, 0 < k < n 8 / 40
8 The Koszul complex Key tool: the Koszul differential κ : Λ k Λ k : Definition of P r Λ k (κω)x(v,..., v k ) = ωx(x, v,..., v k ), X = x x0 0 Pr Λ 0 κ Pr Λ κ κ Pr nλ n 0 C.f., the polynomial de Rham complex 0 Pr Λ 0 d Pr Λ d d Pr nλ n 0 For R 3 0 Pr () X Pr (; R 3 ) Key relation: (dκ + κd)ω = (r + k)ω X Pr (; R 3 ) X Pr 3() 0 ω Hr Λ k (homogeneous polys) Pr Λ k = Pr Λ k + Hr Λ k = Pr Λ k + κhr Λ k+ + dhr+λ k P r Λ k := Pr Λ k + κhr Λ k+ God made P rλ k and P r Λ k, all the rest is the work of man. Hr Λ k = dhr+λ k κhr Λ k+ 9 / / 40 Degrees of freedom Finite element differential forms and classical mixed FEM The other ingredient of a finite element space are the degrees of freedom, i.e., a decomposition of the dual spaces (Pr Λ k (T )) and (P r Λ k (T )), into subspaces associated to subsimplices f of T. DOF for Pr Λ k (T ): to a subsimplex f of dim. d k we associate ω Trf ω η, η P r+k d Λd k (f ) f DOF for Pr Λ k (T ): ω Trf ω η, η Pr+k d Λ d k (f ) Hiptmair f The resulting FE spaces have exactly the continuity required by HΛ k : Theorem. Pr Λ k (T ) = { ω HΛ k () : ω T Pr Λ k (T ) T T }. Similarly for P r. 3 / 40 Pr Λ 0 (T ) = Pr Λ 0 (T ) H Lagrange elts Pr Λ n (T ) = Pr Λ n (T ) L discontinuous elts n = : Pr Λ (T ) H(curl) Raviart Thomas elts n = : P r Λ (T ) H(curl) Brezzi Douglas Marini elts n = 3: Pr Λ (T ) H(curl) Nedelec st kind edge elts n = 3: P r Λ (T ) H(curl) Nedelec nd kind edge elts n = 3: Pr Λ (T ) H(div) Nedelec st kind face elts n = 3: P r Λ (T ) H(div) Nedelec nd kind face elts 3 / 40
9 Finite element de Rham subcomplexes From these spaces we want to build discrete de Rham complexes with bounded projections. It turns out that there are lots of ways to do this ( n for each r). Extreme cases are: 0 P r Λ 0 (T ) d P r Λ (T ) d d Pr Λ n (T ) 0 0 grad Whitney 957, Bossavit 988 curl div 0 Application to Elasticity What else can you do with FEEC? 0 Pr Λ 0 (T ) d Pr Λ (T ) d d Pr nλ n (T ) 0 0 grad curl div 0 33 / 40 Stress displacement mixed finite elements for elasticity Find stress σ : R 3 3 sym, displacement u : R 3 such that Aσ = ɛ(u), div σ = f ( ) σ,u A σ :σ + div σ u + f u dx stationary point H(div;S) L (R n ) Recent progress coming from the FEEC perspective First stable elements based on polynomials, D (Arnold Winther 00), all degrees r : Search for stable finite elements dates back to the 60s, very limited success.... to derive elements that exhibit complete continuity of the appropriate components along interfaces... was achieved by Raviart and Thomas in the case of the heat conduction problem.... Extension to the full stress problem is difficult and as yet such elements have not been successfully noted. Zienkiewicz, Taylor, Zhu The Finite Element Method: Its Basis & Fundamentals, 6th ed., vol., 005 3D stable elements, all degrees r (Arnold Awanou Winther 007): for r = stress space has 6 degrees of freedom (7 per component on average) 35 / / 40
10 A computation using the new elements From Eberhard, Hueber, Jiang, Wohlmuth 006 Mixed formulation with weak symmetry Idea goes back to Fraeijs de Veubeke 975, Amara Thomas 979 In the classical Hellinger Reissner principle, symmetry of the stress tensor (balance of angular momentum) is assumed to hold exactly. Instead we impose it weakly with a Lagrange multiplier (the rotation). ( ) σ,u A σ :σ + div σ u + f u dx stationary point H(div;S) L (R n ) ( ) σ,u,p A σ :σ + div σ u + σ :p + f u dx S.P. H(div;M) L (R n ) L (K) FEEC has led to very simple stable elements 37 / 40 σ u p 38 / 40 Features of the new mixed elements Based on HR formulation with weak symmetry; very natural Lowest degree element is very simple: full P for stress, P0 for displacement and rotation Works for every polynomial degree Works the same in and 3 (or more) dimensions Robust to material constraints like incompressibility Provably stable and convergent 39 / 40 Conclusions Exterior calculus clarifies the nature of physical quantities and the structure of the PDEs involving them. Capturing the right structure on the discrete level can be essential to get stable methods. FEEC provides a very natural framework for the design and understanding of subtle stability issues that arise in the discretization of a wide variety of PDE systems. It brings to bear tools from geometry, topology, and algebra to develop discretizations which are compatible with the geometric, topological, and algebraic structure of the PDE system, and so obtain stability. FEEC has been used to unify, clarify, and refine many known finite element methods. It is a mathematically rigorous theory. The Pr Λ k and P r Λ k spaces are the natural finite element discretizations for differential forms and the de Rham complex. Through FEEC we believe we have completed the long search for just the right mixed finite elements for elasticity. 40 / 40
Finite element exterior calculus: A new approach to the stability of finite elements
Finite element exterior calculus: A new approach to the stability of finite elements Douglas N. Arnold joint with R. Falk and R. Winther Institute for Mathematics and its Applications University of Minnesota
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 informationOutline. 1 Motivating examples. 2 Hilbert complexes and their discretization. 3 Finite element differential forms. 4 The elasticity complex 1 / 25
Outline 1 Motivating examples Hoge theory, Hilbert complexes, an finite element ifferential forms Douglas N. Arnol 2 Hilbert complexes an their iscretization 3 Finite element ifferential forms 4 The elasticity
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 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 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 informationSmoothed projections in finite element exterior calculus
Smooted projections in finite element exterior calculus Ragnar Winter CMA, University of Oslo Norway based on joint work wit: Douglas N. Arnold, Minnesota, Ricard S. Falk, Rutgers, and Snorre H. Cristiansen,
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 informationAposteriorierrorestimatesinFEEC for the de Rham complex
AposteriorierrorestimatesinFEEC 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 and a
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 informationAn abstract Hodge Dirac operator and its stable discretization
An abstract Hodge Dirac operator and its stable discretization Paul Leopardi Mathematical Sciences Institute, Australian National University. For presentation at ICCA, Tartu,. Joint work with Ari Stern,
More informationPublished in Computer Methods in Applied Mechanics and Engineering 198 (2009)
Published in Computer Methods in Applied Mechanics and Engineering 198 (2009) 1660 1672 GEOMETRIC DECOMPOSITIONS AND LOCAL BASES FOR SPACES OF FINITE ELEMENT DIFFERENTIAL FORMS DOUGLAS N. ARNOLD, RICHARD
More informationSpaces of Finite Element Differential Forms
Spaces of Finite Element Differential Forms Douglas N. Arnold Abstract We discuss the construction of finite element spaces of differential forms which satisfy the crucial assumptions of the finite element
More informationFinite element exterior calculus
Finite element exterior calculus Douglas N. Arnold School of Mathematics, University of Minnesota NSF/CBMS Conference, ICERM, June 11 15, 2012 collaborators: R. Falk, R. Winther, G. Awanou, F. Bonizzoni,
More informationA POSTERIORI ERROR ESTIMATES FOR FINITE ELEMENT EXTERIOR CALCULUS: THE DE RHAM COMPLEX
A POSTERIORI ERROR ESTIMATES FOR FINITE ELEMENT EXTERIOR CALCULUS: THE DE RHAM COMPLEX ALAN DEMLOW 1 AND ANIL N. HIRANI 2 Abstract. Finite element exterior calculus (FEEC) has been developed over the past
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 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 informationA Mixed Nonconforming Finite Element for Linear Elasticity
A Mixed Nonconforming Finite Element for Linear Elasticity Zhiqiang Cai, 1 Xiu Ye 2 1 Department of Mathematics, Purdue University, West Lafayette, Indiana 47907-1395 2 Department of Mathematics and Statistics,
More informationNONCONFORMING MIXED ELEMENTS FOR ELASTICITY
Mathematical Models and Methods in Applied Sciences Vol. 13, No. 3 (2003) 295 307 c World Scientific Publishing Company NONCONFORMING MIXED ELEMENTS FOR ELASTICITY DOUGLAS N. ARNOLD Institute for Mathematics
More informationFinite Element Clifford Algebra: A New Toolkit for Evolution Problems
Finite Element Clifford Algebra: A New Toolkit for Evolution Problems Andrew Gillette joint work with Michael Holst Department of Mathematics University of California, San Diego http://ccom.ucsd.edu/ agillette/
More informationarxiv: v2 [math.na] 8 Sep 2015
OMPLEXES OF DISRETE DISTRIBUTIONAL DIFFERENTIAL FORMS AND THEIR HOMOLOGY THEORY MARTIN WERNER LIHT arxiv:1410.6354v2 [math.na] 8 Sep 2015 Abstract. omplexes of discrete distributional differential forms
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 informationFinite Element Methods for Linear Elasticity
Finite Element Metods for Linear Elasticity Ricard S. Falk Department of Matematics - Hill Center Rutgers, Te State University of New Jersey 110 Frelinguysen Rd., Piscataway, NJ 08854-8019 falk@mat.rutgers.edu
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 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 informationSome Elliptic and Subelliptic Complexes from Geometry
Geometric and Nonlinear Partial Differential Equations, Xi An Jiaotong University p. 1/14 Some Elliptic and Subelliptic Complexes from Geometry Michael Eastwood [ based on joint work with Robert Bryant,
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 informationDouble Complexes and Local Cochain Projections
Double Complexes and Local Cochain Projections Richard S. Falk, 1 Ragnar Winther 2 1 Department o Mathematics, Rutgers University, Piscataway, New Jersey 08854 2 Centre o Mathematics or Applications and
More informationA SECOND ELASTICITY ELEMENT USING THE MATRIX BUBBLE WITH TIGHTENED STRESS SYMMETRY
A SECOND ELASTICITY ELEMENT USING THE MATRIX BUBBLE WITH TIGHTENED STRESS SYMMETRY J. GOPALAKRISHNAN AND J. GUZMÁN Abstract. We presented a family of finite elements that use a polynomial space augmented
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 informationarxiv: v2 [math.na] 23 Apr 2016
Improved ZZ A Posteriori Error Estimators for Diffusion Problems: Conforming Linear Elements arxiv:508.009v2 [math.na] 23 Apr 206 Zhiqiang Cai Cuiyu He Shun Zhang May 2, 208 Abstract. In [8], we introduced
More 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 informationIs a dual mesh really necessary?
Is a dual mesh really necessary? Paul Leopardi Mathematical Sciences Institute, Australian National University. For presentation at ICIAM, Vancouver, 2011. Joint work with Ari Stern, UCSD. 19 July 2011
More informationDISCRETE EXTENSION OPERATORS FOR MIXED FINITE ELEMENT SPACES ON LOCALLY REFINED MESHES
DISCRETE EXTENSION OPERATORS FOR MIXED FINITE ELEMENT SPACES ON LOCALLY REFINED MESHES MAR AINSWORTH, JOHNNY GUZMÁN, AND FRANCISCO JAVIER SAYAS Abstract. The existence of uniformly bounded discrete extension
More informationModeling, Structure and Discretization of Mixed-dimensional Partial Differential Equations
Modeling, Structure and Discretization of Mixed-dimensional Partial Differential Equations J. M. Nordbotten and W. M. Boon Department of Mathematics University of Bergen Norway Abstract Mixed-dimensional
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 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 informationFinite Element Multigrid Framework for Mimetic Finite Difference Discretizations
Finite Element Multigrid Framework for Mimetic Finite ifference iscretizations Xiaozhe Hu Tufts University Polytopal Element Methods in Mathematics and Engineering, October 26-28, 2015 Joint work with:
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 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 informationWeak Galerkin Finite Element Scheme and Its Applications
Weak Galerkin Finite Element Scheme and Its Applications Ran Zhang Department of Mathematics Jilin University, China IMS, Singapore February 6, 2015 Talk Outline Motivation WG FEMs: Weak Operators + Stabilizer
More informationMixed finite elements for elasticity on quadrilateral meshes
Adv Comput Math (2015) 41:553 572 DOI 10.1007/s10444-014-9376-x Mixed finite elements for elasticity on quadrilateral meshes Douglas N. Arnold Gerard Awanou Weifeng Qiu Received: 28 June 2013 / Accepted:
More informationSymmetric Matrix Fields in the Finite Element Method
Symmetry 2010, 2, 1375-1389; doi:10.3390/sym2031375 OPEN ACCESS symmetry ISSN 2073-8994 www.mdpi.com/journal/symmetry Review Symmetric Matrix Fields in the Finite Element Method Gerard Awanou Northern
More informationA MULTIGRID METHOD FOR THE PSEUDOSTRESS FORMULATION OF STOKES PROBLEMS
SIAM J. SCI. COMPUT. Vol. 29, No. 5, pp. 2078 2095 c 2007 Society for Industrial and Applied Mathematics A MULTIGRID METHOD FOR THE PSEUDOSTRESS FORMULATION OF STOKES PROBLEMS ZHIQIANG CAI AND YANQIU WANG
More informationA Second Elasticity Element Using the Matrix Bubble
Portland State University PDXScholar Mathematics and Statistics Faculty Publications and Presentations Fariborz Maseeh Department of Mathematics and Statistics 2011 A Second Elasticity Element Using the
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 informationFinite element exterior calculus framework for geophysical fluid dynamics
Finite element exterior calculus framework for geophysical fluid dynamics Colin Cotter Department of Aeronautics Imperial College London Part of ongoing work on UK Gung-Ho Dynamical Core Project funded
More informationUnified Hybridization Of Discontinuous Galerkin, Mixed, And Continuous Galerkin Methods For Second Order Elliptic Problems.
Unified Hybridization Of Discontinuous Galerkin, Mixed, And Continuous Galerkin Methods For Second Order Elliptic Problems. Stefan Girke WWU Münster Institut für Numerische und Angewandte Mathematik 10th
More informationInvariant differential operators on the sphere
Third CoE Talk at the University of Tokyo p. 1/21 Invariant differential operators on the sphere Michael Eastwood Australian National University Third CoE Talk at the University of Tokyo p. 2/21 The round
More informationStability, Consistency, and Convergence: A 21st Century Viewpoint. Convergence, consistency, and stability of discretizations
Stability, Consistency, and Convergence: A 1st Century Viepoint Douglas N. Arnold Scool of Matematics, University of Minnesota Society for Industrial and Applied Matematics Feng Kang Distinguised Lecture
More informationA posteriori error estimates for Maxwell Equations
www.oeaw.ac.at A posteriori error estimates for Maxwell Equations J. Schöberl RICAM-Report 2005-10 www.ricam.oeaw.ac.at A POSTERIORI ERROR ESTIMATES FOR MAXWELL EQUATIONS JOACHIM SCHÖBERL Abstract. Maxwell
More 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 informationTowards Discrete Exterior Calculus and Discrete Mechanics for Numerical Relativity
Towards Discrete Exterior Calculus and Discrete Mechanics for Numerical Relativity Melvin Leok Mathematics, University of Michigan, Ann Arbor. Joint work with Mathieu Desbrun, Anil Hirani, and Jerrold
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 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 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 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 informationHybridized DG methods
Hybridized DG methods University of Florida (Banff International Research Station, November 2007.) Collaborators: Bernardo Cockburn University of Minnesota Raytcho Lazarov Texas A&M University Thanks:
More informationSpace-Time Finite-Element Exterior Calculus and Variational Discretizations of Gauge Field Theories
21st International Symposium on Mathematical Theory of Networks and Systems July 7-11, 2014. Space-Time Finite-Element Exterior Calculus and Variational Discretizations of Gauge Field Theories Joe Salamon
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 informationEQUADIFF 6. Jean-Claude Nédélec Mixed finite element in 3D in H(div) and H(curl) Terms of use:
EQUADIFF 6 Jean-Claude Nédélec Mixed finite element in 3D in H(div) and H(curl) In: Jaromír Vosmanský and Miloš Zlámal (eds.): Equadiff 6, Proceedings of the International Conference on Differential Equations
More informationMATH 676. Finite element methods in scientific computing
MATH 676 Finite element methods in scientific computing Wolfgang Bangerth, Texas A&M University Lecture 33.25: Which element to use Part 2: Saddle point problems Consider the stationary Stokes equations:
More informationModel Reduction via Discretization : Elasticity and Boltzmann
Model Reduction via Discretization : Elasticity and Boltzmann Joachim Schöberl Computational Mathematics in Engineering Institute for Analysis and Scientific Computing Vienna University of Technology IMA
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 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 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 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 informationPhD dissertation defense
Isogeometric mortar methods with applications in contact mechanics PhD dissertation defense Brivadis Ericka Supervisor: Annalisa Buffa Doctor of Philosophy in Computational Mechanics and Advanced Materials,
More informationGeometric Decomposition of Serendipity Finite Element Spaces
Geometric Decomposition of Serendipity Finite Element Spaces Andrew Gillette Department of Mathematics University of California, San Diego http://ccom.ucsd.edu/ agillette/ Joint Mathematics Meetings 2013
More informationCS 468. Differential Geometry for Computer Science. Lecture 13 Tensors and Exterior Calculus
CS 468 Differential Geometry for Computer Science Lecture 13 Tensors and Exterior Calculus Outline Linear and multilinear algebra with an inner product Tensor bundles over a surface Symmetric and alternating
More informationAdaptive methods for control problems with finite-dimensional control space
Adaptive methods for control problems with finite-dimensional control space Saheed Akindeinde and Daniel Wachsmuth Johann Radon Institute for Computational and Applied Mathematics (RICAM) Austrian Academy
More informationQUADRILATERAL H(DIV) FINITE ELEMENTS
QUADRILATERAL H(DIV) FINITE ELEMENTS DOUGLAS N. ARNOLD, DANIELE BOFFI, AND RICHARD S. FALK Abstract. We consider the approximation properties of quadrilateral finite element spaces of vector fields defined
More 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 informationFrom point cloud data to the continuum model of geometry
From point cloud data to the continuum model of geometry J. Harrison University of California, Berkeley July 22, 2007 Continuum model of geometry By the continuum model of geometry we refer to smooth manifolds
More informationFinite Elements for Magnetohydrodynamics and its Optimal Control
Finite Elements for Magnetohydrodynamics and its Karl Kunisch Marco Discacciati (RICAM) FEM Symposium Chemnitz September 25 27, 2006 Overview 1 2 3 What is Magnetohydrodynamics? Magnetohydrodynamics (MHD)
More informationarxiv: v1 [math.na] 19 Dec 2017
Arnold-Winther Mixed Finite Elements for Stokes Eigenvalue Problems Joscha Gedicke Arbaz Khan arxiv:72.0686v [math.na] 9 Dec 207 Abstract This paper is devoted to study the Arnold-Winther mixed finite
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 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 informationMIXED FINITE ELEMENT METHODS FOR PROBLEMS WITH ROBIN BOUNDARY CONDITIONS
MIXED FINITE ELEMENT METHODS FOR PROBLEMS WITH ROBIN BOUNDARY CONDITIONS JUHO KÖNNÖ, DOMINIK SCHÖTZAU, AND ROLF STENBERG Abstract. We derive new a-priori and a-posteriori error estimates for mixed nite
More informationAn a posteriori error estimate and a Comparison Theorem for the nonconforming P 1 element
Calcolo manuscript No. (will be inserted by the editor) An a posteriori error estimate and a Comparison Theorem for the nonconforming P 1 element Dietrich Braess Faculty of Mathematics, Ruhr-University
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 informationMORTAR MULTISCALE FINITE ELEMENT METHODS FOR STOKES-DARCY FLOWS
MORTAR MULTISCALE FINITE ELEMENT METHODS FOR STOKES-DARCY FLOWS VIVETTE GIRAULT, DANAIL VASSILEV, AND IVAN YOTOV Abstract. We investigate mortar multiscale numerical methods for coupled Stokes and Darcy
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 informationLecture 4: Harmonic forms
Lecture 4: Harmonic forms Jonathan Evans 29th September 2010 Jonathan Evans () Lecture 4: Harmonic forms 29th September 2010 1 / 15 Jonathan Evans () Lecture 4: Harmonic forms 29th September 2010 2 / 15
More informationA Posteriori Error Estimation Techniques for Finite Element Methods. Zhiqiang Cai Purdue University
A Posteriori Error Estimation Techniques for Finite Element Methods Zhiqiang Cai Purdue University Department of Mathematics, Purdue University Slide 1, March 16, 2017 Books Ainsworth & Oden, A posteriori
More informationHybrid Discontinuous Galerkin Methods for Fluid Dynamics and Solid Mechanics
Hybrid Discontinuous Galerkin Methods for Fluid Dynamics and Solid Mechanics Joachim Schöberl Computational Mathematics in Engineering Institute for Analysis and Scientific Computing Vienna University
More informationLeast-Squares Finite Element Methods
Pavel В. Bochev Max D. Gunzburger Least-Squares Finite Element Methods Spri ringer Contents Part I Survey of Variational Principles and Associated Finite Element Methods 1 Classical Variational Methods
More informationAMS subject classifications. Primary, 65N15, 65N30, 76D07; Secondary, 35B45, 35J50
A SIMPLE FINITE ELEMENT METHOD FOR THE STOKES EQUATIONS LIN MU AND XIU YE Abstract. The goal of this paper is to introduce a simple finite element method to solve the Stokes equations. This method is in
More informationarxiv: v1 [math-ph] 25 Apr 2013
Higher-order compatible discretization on hexahedrals Jasper Kreeft and Marc Gerritsma arxiv:1304.7018v1 [math-ph] 5 Apr 013 Abstract We derive a compatible discretization method that relies heavily on
More informationSerendipity Finite Element Methods for Cardiac Electrophysiology Models
Serendipity Finite Element Methods for Cardiac Electrophysiology Models Andrew Gillette Department of Mathematics University of California, San Diego http://ccom.ucsd.edu/ agillette/ Andrew Gillette -
More informationFINITE ELEMENT APPROXIMATION OF STOKES-LIKE SYSTEMS WITH IMPLICIT CONSTITUTIVE RELATION
Proceedings of ALGORITMY pp. 9 3 FINITE ELEMENT APPROXIMATION OF STOKES-LIKE SYSTEMS WITH IMPLICIT CONSTITUTIVE RELATION JAN STEBEL Abstract. The paper deals with the numerical simulations of steady flows
More informationTHE LOCAL THEORY OF ELLIPTIC OPERATORS AND THE HODGE THEOREM
THE LOCAL THEORY OF ELLIPTIC OPERATORS AND THE HODGE THEOREM BEN LOWE Abstract. In this paper, we develop the local theory of elliptic operators with a mind to proving the Hodge Decomposition Theorem.
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 informationNumerische Mathematik
Numer. Math. (2002) 92: 401 419 Digital Object Identifier (DOI) 10.1007/s002110100348 Numerische Mathematik Mixed finite elements for elasticity Douglas N. Arnold 1, Ragnar Winther 2 1 Institute for Mathematics
More informationEXACT DE RHAM SEQUENCES OF SPACES DEFINED ON MACRO-ELEMENTS IN TWO AND THREE SPATIAL DIMENSIONS
EXACT DE RHAM SEQUENCES OF SPACES DEFINED ON MACRO-ELEMENTS IN TWO AND THREE SPATIAL DIMENSIONS JOSEPH E. PASCIAK AND PANAYOT S. VASSILEVSKI Abstract. This paper proposes new finite element spaces that
More informationComplex manifolds, Kahler metrics, differential and harmonic forms
Complex manifolds, Kahler metrics, differential and harmonic forms Cattani June 16, 2010 1 Lecture 1 Definition 1.1 (Complex Manifold). A complex manifold is a manifold with coordinates holomorphic on
More informationMultigrid Methods for Saddle Point Problems
Multigrid Methods for Saddle Point Problems Susanne C. Brenner Department of Mathematics and Center for Computation & Technology Louisiana State University Advances in Mathematics of Finite Elements (In
More informationUnified A Posteriori Error Control for all Nonstandard Finite Elements 1
Unified A Posteriori Error Control for all Nonstandard Finite Elements 1 Martin Eigel C. Carstensen, C. Löbhard, R.H.W. Hoppe Humboldt-Universität zu Berlin 19.05.2010 1 we know of Guidelines for Applicants
More informationA Locking-Free MHM Method for Elasticity
Trabalho apresentado no CNMAC, Gramado - RS, 2016. Proceeding Series of the Brazilian Society of Computational and Applied Mathematics A Locking-Free MHM Method for Elasticity Weslley S. Pereira 1 Frédéric
More informationMIXED FINITE ELEMENT METHODS FOR LINEAR VISCOELASTICITY USING WEAK SYMMETRY
MIXED FINITE ELEMENT METHODS FOR LINEAR VISCOELASTICITY USING WEAK SYMMETRY MARIE E. ROGNES AND RAGNAR WINTHER Abstract. Small deformations of a viscoelastic body are considered through the linear Maxwell
More information