Finite Element Clifford Algebra: A New Toolkit for Evolution Problems
|
|
- Nathaniel Hutchinson
- 5 years ago
- Views:
Transcription
1 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 agillette/ Andrew Gillette - UCSD Finite ( ) Element Clifford Algebra SIAM PD11 - Nov / 18
2 Motivation Poisson s equation: Given f find u(x) such that { 0 = u + f in Ω R n u = 0 on Ω a Heat equation: Given f and g, find u(x, t) such that u t = u + f in Ω R n, for t > 0, u = 0 on Ω, for t > 0, u t=0 = g in Ω Finite element exterior calculus (FEEC) provides: abstract framework for analyzing numerical approximation of elliptic PDEs classification of stable finite element methods with optimal convergence rates How can the FEEC framework be expanded to classify stable finite element methods for evolutionary PDEs? Andrew Gillette - UCSD Finite ( ) Element Clifford Algebra SIAM PD11 - Nov / 18
3 Outline of approach Two possible methods for extending Finite Element Exterior Calculus: Semi-discrete: Finite element method in space, ODE in time t Ω domain: Ω [0, T ] R n R solution basis: φ h t=t0 : Ω R, for each t 0 [0, T ] error analysis: FEEC + Bochner space theory Fully discrete: Finite element method in space and time Ω domain: solution basis: error analysis: Ω [0, T ] R n R φ h : Ω [0, T ] R Finite Element Clifford Algebra This talk: Initial results on semi-discrete approach + a preview of FECA Andrew Gillette - UCSD Finite ( ) Element Clifford Algebra SIAM PD11 - Nov / 18
4 Table of Contents 1 Motivation 2 Background: FEEC, Bochner Spaces, Semi-Discrete methods 3 Results: New error estimates in Bochner norms 4 Preview: Why FECA is needed Andrew Gillette - UCSD Finite ( ) Element Clifford Algebra SIAM PD11 - Nov / 18
5 Outline 1 Motivation 2 Background: FEEC, Bochner Spaces, Semi-Discrete methods 3 Results: New error estimates in Bochner norms 4 Preview: Why FECA is needed Andrew Gillette - UCSD Finite ( ) Element Clifford Algebra SIAM PD11 - Nov / 18
6 Finite Element Exterior Calculus in context Consider a mixed method for Poisson s problem on a domain Ω R n : continuous u + f = 0, u H 2 mixed weak (div σ, φ) + (f, φ) = 0, φ L 2 = Λ n (σ, ω) + (u, div ω) = 0, ω H(div) = Λ n 1 mixed FEM (div σ h, φ h ) + (f, φ h ) = 0, φ h Λ n h L 2 (σ h, ω h ) + (u h, div ω h ) = 0, ω h Λ n 1 h H(div) Major Conclusions from FEEC The finite elements spaces Λ n 1 h and Λ n 1 h should be chosen from two classes of piecewise polynomial spaces, denoted P r Λ k h and Pr If this choice is made in a compatible manner implied by the exterior calculus structure, then optimal a priori error estimates are guaranteed Λ k h Andrew Gillette - UCSD Finite ( ) Element Clifford Algebra SIAM PD11 - Nov / 18
7 Finite Element Exterior Calculus in context Theorem [Arnold, Falk, Winther; Bulletin of AMS, 2010] Assume the elliptic regularity estimate u H s+2 + u H s+1 + σ H s+1 + div σ H s c f H s holds for 0 s s max. Choose finite element spaces Λ n 1 h =, Λn h = P r+1 Λ n 1 (T ) or P r+1 Λn 1 (T ) P r+1 Λn (T ) or P r Λ n (T ) Then for 0 s s max, the following error estimates hold { ch f L 2 if Λ n h = P 1 u u h L 2 Λn (T ), ch 2+s f H s otherwise, if s r 1 σ h σ L 2 ch s+1 f H s if { Λ n 1 h = P r+1 Λ n 1, s r + 1, Λ n 1 h = P r+1 Λn 1, s r, div (σ h σ) L 2 ch s f H s, if s r + 1. Andrew Gillette - UCSD Finite ( ) Element Clifford Algebra SIAM PD11 - Nov / 18
8 Semi-discrete Mixed Formulation Consider a mixed method for the heat equation on Ω R n for t I := [0, T ]. continuous u t u = f, u t=0 = g. mixed weak (u t, φ) (div σ, φ) = (f, φ), φ Λ n, t I, (σ, ω) + (u, div ω) = 0, ω Λ n 1, t I, u t=0 = g. mixed FEM (u h,t, φ h ) (div σ h, φ h ) = (f, φ h ), φ h Λ n h, t I, (σ h, ω h ) + (u h, div ω h ) = 0, ω h Λ n 1 h, t I, u h t=0 = g h. linear system AU t BΣ = F B T U + DΣ = 0 AU t + BD 1 B T U = F Andrew Gillette - UCSD Finite ( ) Element Clifford Algebra SIAM PD11 - Nov / 18
9 Semi-discrete Error Bounds Theorem [Thomée; Galerkin FEM for Parabolic Problems, 1997] Fix n = 2 and set Λ 2 h := discontinuous linear, Λ 1 h := Raviart-Thomas elements. Let g h be the solution to the elliptic problem with f = g. Then for t 0: t ) u h (t) u(t) L 2 ch ( u(t) 2 H 2 + u t H 2 ds, 0 ( t ) 1/2 ) σ h (t) σ(t) L 2 ch ( u(t) 2 H 3 + u t 2 H ds. 2 0 Homogeneous case (f = 0), g h as above, t 0: u h (t) u(t) L 2 ch 2 g H 2, if g Ḣ2, σ h (t) σ(t) L 2 ch 3 g H 3, if g Ḣ3. Homogeneous case (f = 0), g h := orthogonal projection of g on to Λ 2 h, t > 0: u h (t) u(t) L 2 ch 2 t 1 g L 2 σ h (t) σ(t) L 2 ch 2 t 3/2 g L 2 Note: These bounds are space-only and restricted to the case n = 2. Andrew Gillette - UCSD Finite ( ) Element Clifford Algebra SIAM PD11 - Nov / 18
10 Bochner spaces and norms Our new error bounds will employ the theory of Bochner spaces Definition Let X be a Banach space and I = (0, T ). Define Equip this space with the norm C(I, X) := {u : I X u bounded and continuous} u C(I,X) := sup u(t) X. t I The Bochner space L P (I, X) is defined to be the completion of C(I, X) with respect to the norm: ( ) 1/p u L p (I,X) := u(t) p X dt. We combine notations to get Bochner differential form spaces: I L 2 X k := L 2 (I, L 2 Λ k (Ω)) These are parametrized differential form spaces. Andrew Gillette - UCSD Finite ( ) Element Clifford Algebra SIAM PD11 - Nov / 18
11 Outline 1 Motivation 2 Background: FEEC, Bochner Spaces, Semi-Discrete methods 3 Results: New error estimates in Bochner norms 4 Preview: Why FECA is needed Andrew Gillette - UCSD Finite ( ) Element Clifford Algebra SIAM PD11 - Nov / 18
12 Bochner-FEEC Parabolic Error Estimates We combine the FEEC and parabolic error estimates to derive the following. Theorem [G, Holst, 2011] Let n 2 and fix I := [0, T ]. Suppose regularity estimate u(t) H s+2 + u(t) H s+1 + σ(t) H s+1 + div σ(t) H s c f (t) H s holds for 0 s s max and t I. Choose finite element spaces P r+1 Λ n 1 (T ) P Λ n 1 h = or P, r+1 Λn (T ) Λn h = or r+1 Λn 1 (T ) P r Λ n (T ) Then for 0 s s max and g h the solution to the elliptic problem we have ( ch f L 2 (I,L 2 ) + ) T f t L 1 (I,L 2 ) if Λ n h = P 1 Λn (T ) u h u L 2 X n ( ch 2+s f L 2 (I,H s ) + ) T f t L 1 (I,H s ) otherwise, if s r 1 and... Andrew Gillette - UCSD Finite ( ) Element Clifford Algebra SIAM PD11 - Nov / 18
13 Bochner-FEEC Paraoblic Error Estimates Theorem [G, Holst, 2011] Let n 2 and fix I := [0, T ]. Suppose regularity estimate u(t) H s+2 + u(t) H s+1 + σ(t) H s+1 + div σ(t) H s c f (t) H s holds for 0 s s max and t I. Choose finite element spaces P r+1 Λ n 1 (T ) P Λ n 1 h = or P, r+1 Λn (T ) Λn h = or r+1 Λn 1 (T ) P r Λ n (T ) Then for 0 s s max and g h the solution to the elliptic problem we have: Λ n 1 h = P 1 Λ n 1 (T ), s 1 If or Λ n 1 h = P and Λn h = P 1 Λn (T ), then 1 Λn 1 (T ), s = 0 ( ) σ h σ L 2 X c h 1+s f n 1 L 2 (I,H s ) + h T f t L 2 (I,L 2 ) For any other choice of spaces, if s r 1, ( ) σ h σ L 2 X c h 1+s f n 1 L 2 (I,H s ) + h2+s T f t L 2 (I,L 2 ) Andrew Gillette - UCSD Finite ( ) Element Clifford Algebra SIAM PD11 - Nov / 18
14 Proof and Significance Key idea of the proof: u(t) u h (t) L 2 }{{} error between weak and semi-discrete u(t) ũ h (t) L 2 }{{} error between weak and time-ignorant elliptic + ũ h (t) u h (t) L 2 }{{} error between time-ignorant elliptic and semi-discrete Significance of the error estimates These results give a priori estimates of convergence rates for the semi-discrete Galerkin FEM for the heat equation. By using the FEEC framework, we have classified choices of semi-discrete finite element spaces that guarantee optimal convergence rates. The results hold for arbitrary spatial dimension n, not just n = 2. For the homogeneous case (f = 0) with sufficiently regular g, we expect to find stronger error estimates akin to Thomée s. G, HOLST, Finite Element Exterior Calculus for Evolution Problems, in preparation. Andrew Gillette - UCSD Finite ( ) Element Clifford Algebra SIAM PD11 - Nov / 18
15 Outline 1 Motivation 2 Background: FEEC, Bochner Spaces, Semi-Discrete methods 3 Results: New error estimates in Bochner norms 4 Preview: Why FECA is needed Andrew Gillette - UCSD Finite ( ) Element Clifford Algebra SIAM PD11 - Nov / 18
16 The Bochner Complex FEEC theory studies discretizations of the L 2 derham complex: 0 HΛ 0 d Ω (grad) HΛ 1 d Ω d Ω (div) HΛ n d Ω 0 We can define a parametrized exterior derivative operator on Bochner spaces: d : HX k HX k+1 where (dµ)(t) := d Ω (µ(t)). This gives rise to a Bochner domain complex: 0 HX 0 d HX 1 d d HX n d 0 Ω For a fully discrete method, we need an exterior derivative operator on spacetime elements which can distinguish spacelike and timelike dimensions. Such an operator needs the Lorentzian signature of basis elements - a tool available in Clifford Algebra (or Geometric Calculus) but not exterior calculus. Andrew Gillette - UCSD Finite ( ) Element Clifford Algebra SIAM PD11 - Nov / 18
17 Beyond the derham Complex... The derivative operator in Clifford algebra is a formal sum of d and its adjoint: := d + δ [0, T ] Λ k+2 d The derham complex appears as diagonals in a full Clifford complex The Bochner complex appears as parametrizations of these diagonals [0, T ] Λ k d δ Λ k+1 Λ k 1 δ d Λ k δ Λ k 2 Finite Element Clifford Algebra will study discretizations of this larger complex. Andrew Gillette - UCSD Finite ( ) Element Clifford Algebra SIAM PD11 - Nov / 18
18 Questions? Slides and pre-prints available at agillette Andrew Gillette - UCSD Finite ( ) Element Clifford Algebra SIAM PD11 - Nov / 18
Geometric 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 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 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 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 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 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 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 informationSECOND ORDER TIME DISCONTINUOUS GALERKIN METHOD FOR NONLINEAR CONVECTION-DIFFUSION PROBLEMS
Proceedings of ALGORITMY 2009 pp. 1 10 SECOND ORDER TIME DISCONTINUOUS GALERKIN METHOD FOR NONLINEAR CONVECTION-DIFFUSION PROBLEMS MILOSLAV VLASÁK Abstract. We deal with a numerical solution of a scalar
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 informationA space-time Trefftz method for the second order wave equation
A space-time Trefftz method for the second order wave equation Lehel Banjai The Maxwell Institute for Mathematical Sciences Heriot-Watt University, Edinburgh & Department of Mathematics, University of
More informationFinite 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 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 informationFinite Differences for Differential Equations 28 PART II. Finite Difference Methods for Differential Equations
Finite Differences for Differential Equations 28 PART II Finite Difference Methods for Differential Equations Finite Differences for Differential Equations 29 BOUNDARY VALUE PROBLEMS (I) Solving a TWO
More informationFinite element approximation of the stochastic heat equation with additive noise
p. 1/32 Finite element approximation of the stochastic heat equation with additive noise Stig Larsson p. 2/32 Outline Stochastic heat equation with additive noise du u dt = dw, x D, t > u =, x D, t > u()
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 informationNumerical Analysis and Methods for PDE I
Numerical Analysis and Methods for PDE I A. J. Meir Department of Mathematics and Statistics Auburn University US-Africa Advanced Study Institute on Analysis, Dynamical Systems, and Mathematical Modeling
More informationUne méthode de pénalisation par face pour l approximation des équations de Navier-Stokes à nombre de Reynolds élevé
Une méthode de pénalisation par face pour l approximation des équations de Navier-Stokes à nombre de Reynolds élevé CMCS/IACS Ecole Polytechnique Federale de Lausanne Erik.Burman@epfl.ch Méthodes Numériques
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 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 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 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 informationA space-time Trefftz method for the second order wave equation
A space-time Trefftz method for the second order wave equation Lehel Banjai The Maxwell Institute for Mathematical Sciences Heriot-Watt University, Edinburgh Rome, 10th Apr 2017 Joint work with: Emmanuil
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 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 informationMotivations. Outline. Finite element exterior calculus and the geometrical basis of numerical stability. References. Douglas N.
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
More informationINTRODUCTION TO PDEs
INTRODUCTION TO PDEs In this course we are interested in the numerical approximation of PDEs using finite difference methods (FDM). We will use some simple prototype boundary value problems (BVP) and initial
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 informationPast, present and space-time
Past, present and space-time Arnold Reusken Chair for Numerical Mathematics RWTH Aachen Utrecht, 12.11.2015 Reusken (RWTH Aachen) Past, present and space-time Utrecht, 12.11.2015 1 / 20 Outline Past. Past
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 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 informationarxiv: v1 [math.na] 29 Feb 2016
EFFECTIVE IMPLEMENTATION OF THE WEAK GALERKIN FINITE ELEMENT METHODS FOR THE BIHARMONIC EQUATION LIN MU, JUNPING WANG, AND XIU YE Abstract. arxiv:1602.08817v1 [math.na] 29 Feb 2016 The weak Galerkin (WG)
More informationImplementation of Sparse Wavelet-Galerkin FEM for Stochastic PDEs
Implementation of Sparse Wavelet-Galerkin FEM for Stochastic PDEs Roman Andreev ETH ZÜRICH / 29 JAN 29 TOC of the Talk Motivation & Set-Up Model Problem Stochastic Galerkin FEM Conclusions & Outlook Motivation
More informationELLIPTIC RECONSTRUCTION AND A POSTERIORI ERROR ESTIMATES FOR PARABOLIC PROBLEMS
ELLIPTIC RECONSTRUCTION AND A POSTERIORI ERROR ESTIMATES FOR PARABOLIC PROBLEMS CHARALAMBOS MAKRIDAKIS AND RICARDO H. NOCHETTO Abstract. It is known that the energy technique for a posteriori error analysis
More informationSUPERCONVERGENCE PROPERTIES FOR OPTIMAL CONTROL PROBLEMS DISCRETIZED BY PIECEWISE LINEAR AND DISCONTINUOUS FUNCTIONS
SUPERCONVERGENCE PROPERTIES FOR OPTIMAL CONTROL PROBLEMS DISCRETIZED BY PIECEWISE LINEAR AND DISCONTINUOUS FUNCTIONS A. RÖSCH AND R. SIMON Abstract. An optimal control problem for an elliptic equation
More informationA local-structure-preserving local discontinuous Galerkin method for the Laplace equation
A local-structure-preserving local discontinuous Galerkin method for the Laplace equation Fengyan Li 1 and Chi-Wang Shu 2 Abstract In this paper, we present a local-structure-preserving local discontinuous
More informationA simplicial gauge theory on spacetime
A simplicial gauge theory on spacetime Tore G. Halvorsen, NTNU, Trondheim Norway. Joint work with Snorre H. Christiansen, UiO. 1 Abstract In this talk I will introduce the variational form of the SU(N)
More informationIterative methods for positive definite linear systems with a complex shift
Iterative methods for positive definite linear systems with a complex shift William McLean, University of New South Wales Vidar Thomée, Chalmers University November 4, 2011 Outline 1. Numerical solution
More informationTutorial 2. Introduction to numerical schemes
236861 Numerical Geometry of Images Tutorial 2 Introduction to numerical schemes c 2012 Classifying PDEs Looking at the PDE Au xx + 2Bu xy + Cu yy + Du x + Eu y + Fu +.. = 0, and its discriminant, B 2
More informationSpace-time sparse discretization of linear parabolic equations
Space-time sparse discretization of linear parabolic equations Roman Andreev August 2, 200 Seminar for Applied Mathematics, ETH Zürich, Switzerland Support by SNF Grant No. PDFMP2-27034/ Part of PhD thesis
More informationA Primal-Dual Weak Galerkin Finite Element Method for Second Order Elliptic Equations in Non-Divergence Form
A Primal-Dual Weak Galerkin Finite Element Method for Second Order Elliptic Equations in Non-Divergence Form Chunmei Wang Visiting Assistant Professor School of Mathematics Georgia Institute of Technology
More informationFINITE ELEMENT EXTERIOR CALCULUS FOR PARABOLIC EVOLUTION PROBLEMS ON RIEMANNIAN HYPERSURFACES
FINITE ELEMENT EXTERIOR CALCULUS FOR PARABOLIC EVOLUTION PROBLEMS ON RIEMANNIAN HYPERSURFACES MICHAEL HOLST AND CHRIS TIEE ABSTRACT. Over te last ten years, te Finite Element Exterior Calculus (FEEC) as
More informationNumerical Methods for PDEs
Numerical Methods for PDEs Partial Differential Equations (Lecture 1, Week 1) Markus Schmuck Department of Mathematics and Maxwell Institute for Mathematical Sciences Heriot-Watt University, Edinburgh
More informationJuan Vicente Gutiérrez Santacreu Rafael Rodríguez Galván. Departamento de Matemática Aplicada I Universidad de Sevilla
Doc-Course: Partial Differential Equations: Analysis, Numerics and Control Research Unit 3: Numerical Methods for PDEs Part I: Finite Element Method: Elliptic and Parabolic Equations Juan Vicente Gutiérrez
More informationYunrong Zhu: Curriculum Vitae
Yunrong Zhu: Curriculum Vitae Department of Mathematics Office : 5824 AP&M Building University of California, San Diego Phone : 1-(858) 534-5862 9500 Gilman Drive, Dept. 0112 Email : zhu@math.ucsd.edu
More informationSharp Sobolev Strichartz estimates for the free Schrödinger propagator
Sharp Sobolev Strichartz estimates for the free Schrödinger propagator Neal Bez, Chris Jeavons and Nikolaos Pattakos Abstract. We consider gaussian extremisability of sharp linear Sobolev Strichartz estimates
More informationFrom Completing the Squares and Orthogonal Projection to Finite Element Methods
From Completing the Squares and Orthogonal Projection to Finite Element Methods Mo MU Background In scientific computing, it is important to start with an appropriate model in order to design effective
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 informationA spacetime DPG method for acoustic waves
A spacetime DPG method for acoustic waves Rainer Schneckenleitner JKU Linz January 30, 2018 ( JKU Linz ) Spacetime DPG method January 30, 2018 1 / 41 Overview 1 The transient wave problem 2 The broken
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 informationAnalysis and Applications of Polygonal and Serendipity Finite Element Methods
Analysis and Applications of Polygonal and Serendipity Finite Element Methods Andrew Gillette Department of Mathematics University of California, San Diego http://ccom.ucsd.edu/ agillette/ Andrew Gillette
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 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 informationSimple Examples on Rectangular Domains
84 Chapter 5 Simple Examples on Rectangular Domains In this chapter we consider simple elliptic boundary value problems in rectangular domains in R 2 or R 3 ; our prototype example is the Poisson equation
More 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 information[2] (a) Develop and describe the piecewise linear Galerkin finite element approximation of,
269 C, Vese Practice problems [1] Write the differential equation u + u = f(x, y), (x, y) Ω u = 1 (x, y) Ω 1 n + u = x (x, y) Ω 2, Ω = {(x, y) x 2 + y 2 < 1}, Ω 1 = {(x, y) x 2 + y 2 = 1, x 0}, Ω 2 = {(x,
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 informationAPPLIED PARTIAL DIFFERENTIAL EQUATIONS
APPLIED PARTIAL DIFFERENTIAL EQUATIONS AN I N T R O D U C T I O N ALAN JEFFREY University of Newcastle-upon-Tyne ACADEMIC PRESS An imprint of Elsevier Science Amsterdam Boston London New York Oxford Paris
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 11 Partial Differential Equations Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign Copyright c 2002.
More informationOn Pressure Stabilization Method and Projection Method for Unsteady Navier-Stokes Equations 1
On Pressure Stabilization Method and Projection Method for Unsteady Navier-Stokes Equations 1 Jie Shen Department of Mathematics, Penn State University University Park, PA 1682 Abstract. We present some
More informationChapter Two: Numerical Methods for Elliptic PDEs. 1 Finite Difference Methods for Elliptic PDEs
Chapter Two: Numerical Methods for Elliptic PDEs Finite Difference Methods for Elliptic PDEs.. Finite difference scheme. We consider a simple example u := subject to Dirichlet boundary conditions ( ) u
More informationWe denote the space of distributions on Ω by D ( Ω) 2.
Sep. 1 0, 008 Distributions Distributions are generalized functions. Some familiarity with the theory of distributions helps understanding of various function spaces which play important roles in the study
More informationTowards parametric model order reduction for nonlinear PDE systems in networks
Towards parametric model order reduction for nonlinear PDE systems in networks MoRePas II 2012 Michael Hinze Martin Kunkel Ulrich Matthes Morten Vierling Andreas Steinbrecher Tatjana Stykel Fachbereich
More informationSwitching, sparse and averaged control
Switching, sparse and averaged control Enrique Zuazua Ikerbasque & BCAM Basque Center for Applied Mathematics Bilbao - Basque Country- Spain zuazua@bcamath.org http://www.bcamath.org/zuazua/ WG-BCAM, February
More informationEfficient Solvers for Stochastic Finite Element Saddle Point Problems
Efficient Solvers for Stochastic Finite Element Saddle Point Problems Catherine E. Powell c.powell@manchester.ac.uk School of Mathematics University of Manchester, UK Efficient Solvers for Stochastic Finite
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 informationOptimal control in fluid mechanics by finite elements with symmetric stabilization
Computational Sciences Center Optimal control in fluid mechanics by finite elements with symmetric stabilization Malte Braack Mathematisches Seminar Christian-Albrechts-Universität zu Kiel VMS Worshop
More informationGreedy Control. Enrique Zuazua 1 2
Greedy Control Enrique Zuazua 1 2 DeustoTech - Bilbao, Basque Country, Spain Universidad Autónoma de Madrid, Spain Visiting Fellow of LJLL-UPMC, Paris enrique.zuazua@deusto.es http://enzuazua.net X ENAMA,
More informationBasic Concepts of Adaptive Finite Element Methods for Elliptic Boundary Value Problems
Basic Concepts of Adaptive Finite lement Methods for lliptic Boundary Value Problems Ronald H.W. Hoppe 1,2 1 Department of Mathematics, University of Houston 2 Institute of Mathematics, University of Augsburg
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University The Residual and Error of Finite Element Solutions Mixed BVP of Poisson Equation
More informationAn introduction to Isogeometric Analysis
An introduction to Isogeometric Analysis A. Buffa 1, G. Sangalli 1,2, R. Vázquez 1 1 Istituto di Matematica Applicata e Tecnologie Informatiche - Pavia Consiglio Nazionale delle Ricerche 2 Dipartimento
More informationMULTIGRID PRECONDITIONING IN H(div) ON NON-CONVEX POLYGONS* Dedicated to Professor Jim Douglas, Jr. on the occasion of his seventieth birthday.
MULTIGRID PRECONDITIONING IN H(div) ON NON-CONVEX POLYGONS* DOUGLAS N ARNOLD, RICHARD S FALK, and RAGNAR WINTHER Dedicated to Professor Jim Douglas, Jr on the occasion of his seventieth birthday Abstract
More informationA Very General Hamilton Jacobi Theorem
University of Salerno 9th AIMS ICDSDEA Orlando (Florida) USA, July 1 5, 2012 Introduction: Generalized Hamilton-Jacobi Theorem Hamilton-Jacobi (HJ) Theorem: Consider a Hamiltonian system with Hamiltonian
More information(1) u (t) = f(t, u(t)), 0 t a.
I. Introduction 1. Ordinary Differential Equations. In most introductions to ordinary differential equations one learns a variety of methods for certain classes of equations, but the issues of existence
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 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 informationPartial Differential Equations
Part II Partial Differential Equations Year 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2015 Paper 4, Section II 29E Partial Differential Equations 72 (a) Show that the Cauchy problem for u(x,
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 information1. Differential Equations (ODE and PDE)
1. Differential Equations (ODE and PDE) 1.1. Ordinary Differential Equations (ODE) So far we have dealt with Ordinary Differential Equations (ODE): involve derivatives with respect to only one variable
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 # 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 informationNew DPG techniques for designing numerical schemes
New DPG techniques for designing numerical schemes Jay Gopalakrishnan University of Florida Collaborator: Leszek Demkowicz October 2009 Massachusetts Institute of Technology, Boston Thanks: NSF Jay Gopalakrishnan
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 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 informationMEASURE VALUED DIRECTIONAL SPARSITY FOR PARABOLIC OPTIMAL CONTROL PROBLEMS
MEASURE VALUED DIRECTIONAL SPARSITY FOR PARABOLIC OPTIMAL CONTROL PROBLEMS KARL KUNISCH, KONSTANTIN PIEPER, AND BORIS VEXLER Abstract. A directional sparsity framework allowing for measure valued controls
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 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 informationPreconditioning in H(div) and Applications
1 Preconditioning in H(div) and Applications Douglas N. Arnold 1, Ricard S. Falk 2 and Ragnar Winter 3 4 Abstract. Summarizing te work of [AFW97], we sow ow to construct preconditioners using domain decomposition
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 informationApplied Math Qualifying Exam 11 October Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems.
Printed Name: Signature: Applied Math Qualifying Exam 11 October 2014 Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems. 2 Part 1 (1) Let Ω be an open subset of R
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 informationProper Orthogonal Decomposition. POD for PDE Constrained Optimization. Stefan Volkwein
Proper Orthogonal Decomposition for PDE Constrained Optimization Institute of Mathematics and Statistics, University of Constance Joined work with F. Diwoky, M. Hinze, D. Hömberg, M. Kahlbacher, E. Kammann,
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 informationLocal flux mimetic finite difference methods
Local flux mimetic finite difference methods Konstantin Lipnikov Mikhail Shashkov Ivan Yotov November 4, 2005 Abstract We develop a local flux mimetic finite difference method for second order elliptic
More informationDS-GA 1002 Lecture notes 0 Fall Linear Algebra. These notes provide a review of basic concepts in linear algebra.
DS-GA 1002 Lecture notes 0 Fall 2016 Linear Algebra These notes provide a review of basic concepts in linear algebra. 1 Vector spaces You are no doubt familiar with vectors in R 2 or R 3, i.e. [ ] 1.1
More informationYunrong Zhu: Curriculum Vitae
Yunrong Zhu: Curriculum Vitae Appointments Education Department of Mathematics & Statistics Office : Physical Sciences 328B Idaho State University, Pocatello Phone : 1-(208) 282-3819 921 South 8th Ave.,
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 informationZdzislaw Brzeźniak. Department of Mathematics University of York. joint works with Mauro Mariani (Roma 1) and Gaurav Dhariwal (York)
Navier-Stokes equations with constrained L 2 energy of the solution Zdzislaw Brzeźniak Department of Mathematics University of York joint works with Mauro Mariani (Roma 1) and Gaurav Dhariwal (York) Stochastic
More informationSpace-time XFEM for two-phase mass transport
Space-time XFEM for two-phase mass transport Space-time XFEM for two-phase mass transport Christoph Lehrenfeld joint work with Arnold Reusken EFEF, Prague, June 5-6th 2015 Christoph Lehrenfeld EFEF, Prague,
More informationMultiscale Diffusion Modeling in Charged and Crowded Biological Environments
Multiscale Diffusion Modeling in Charged and Crowded Biological Environments Andrew Gillette Department of Mathematics University of Arizona joint work with Pete Kekenes-Huskey (U. Kentucky) and J. Andrew
More informationLaplace s Equation FEM Methods. Jacob White. Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy, Jaime Peraire and Tony Patera
Introduction to Simulation - Lecture 19 Laplace s Equation FEM Methods Jacob White Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy, Jaime Peraire and Tony Patera Outline for Poisson Equation
More information