Hybridized Discontinuous Galerkin Methods
|
|
- Joshua Watts
- 6 years ago
- Views:
Transcription
1 Hybridized Discontinuous Galerkin Methods Theory and Christian Waluga Joint work with Herbert Egger (Uni Graz) 1st DUNE User Meeting, Stuttgart Christian Waluga (AICES) HDG Methods October 6-8, / 26
2 Outline Outline 1 Hybridization: From DG to HDG Preliminaries A simple DG method Hybridization: HDG Example: Oseen problem 2 Hybridized Mortar Methods Relation between HDG and Hybridized Mortar Example: Stokes problem 3 Conclusion Christian Waluga (AICES) HDG Methods October 6-8, / 26
3 Outline Hybridization: From DG to HDG 1 Hybridization: From DG to HDG Preliminaries A simple DG method Hybridization: HDG Example: Oseen problem 2 Hybridized Mortar Methods Relation between HDG and Hybridized Mortar Example: Stokes problem 3 Conclusion Christian Waluga (AICES) HDG Methods October 6-8, / 26
4 Model problem Hybridization: From DG to HDG Preliminaries Ω Model problem: Poisson s equation Given f L 2 (Ω), find u H 1 (Ω), such that Ω u = f, u = 0, on Ω on Ω Assume Ω R d is bounded, sufficiently regular domain. Notation: (u, v) D := D u v dx D Ω v 0,D := (v, v) D, v 1,D := ( v, v) D We will abbreviate H 1 (Ω) := [ H 1 (Ω) ] d etc. Christian Waluga (AICES) HDG Methods October 6-8, / 26
5 Hybridization: From DG to HDG Meshes and basic notation Preliminaries Meshes: (hanging nodes also possible) T i T e ij j T h := { T i }, where Ω = T i T h := { T i } E h := { e ij : e ij = T i T j, i > j } Broken Sobolev spaces: H s (T h ) := { v L 2 (Ω) : v T H s (T ) for all T T h } (u, v) Th := (u, v) L 2 (T ), T T h v Th := (v, v) 1/2 T h, etc. Jump and average: v := v Ti v Tj, {v} := 1 2 (v T i + v Tj ) Note: functions defined on E h can be interpreted as functions on T h. Christian Waluga (AICES) HDG Methods October 6-8, / 26
6 A simple DG method Hybridization: From DG to HDG A simple DG method Example: Symmetric Interior Penalty Galerkin (SIPG 1 ) Find u H s (T h ) H 1 0 (Ω), s > 3/2, such that a SIP G h (u, v) = (f, v) Th, v H s (T h ), where we define a symmetric bilinear form a SIP h G (u, v) :=( u, v) Th ({ n u}, v ) Eh ({ n v}, u ) Eh + τ( u, v ) Eh. Here, τ := α p2 h is a penalty parameter, with α R+. Natural energy norm for the discrete analysis: v 1,h := ( v 2 T h + τ 1/2 v 2 E h ) 1/2. 1 e.g. [Arnold 1982], [Rivière 2008] Christian Waluga (AICES) HDG Methods October 6-8, / 26
7 Hybridization: From DG to HDG Hybridized DG method Hybridization: HDG Example: Hybridized Symmetric Interior Penalty Galerkin Find (u, û) (H s (T h ) H0 1 (Ω)) L 2 (E h ), s > 3/2, such that a h (u, û; v, ˆv) = (f, v) Th, (v, ˆv) H s (T h ) L 2 (E h ), where û := {u} and ˆv := {v}. We define a symmetric bilinear form a h (u, û; v, ˆv) :=( u, v) Th ( n u, v ˆv) Th ( n v, u û) Th + τ(u û, v ˆv) Th. Here, τ := α p2 h is a penalty parameter, with α R+. Natural energy norm for the discrete analysis: (v, ˆv) 1,h := ( v 2 T h + τ 1/2 (v ˆv) 2 T h ) 1/2. Christian Waluga (AICES) HDG Methods October 6-8, / 26
8 Hybridization: From DG to HDG Hybridization: HDG Hybridization: From DG to HDG (4/4) Given a triangulation T h, define finite dimensional spaces V h and V h : V h := { v H 1 (T h ) : v T P p (T ), T T h } ˆV h := { ˆv L 2 (E h ) : v E P p (E), E E h, ˆv = 0 on Ω } Discrete problem: Find (u h, û h ) V h V h such that a h (u h, û h ; v h, ˆv h ) = (f, v h ) Ω, (v h, ˆv h ) V h V h. Stability analysis: For any α > α 0, there holds 2 a h (v h, ˆv h ; v h, ˆv h ) 1 2 (v h, ˆv h ) 2 1,h, (v h, ˆv h ) V h V h, where α 0 depends on the element shape. One can choose α large enough (e.g. α = 10) or explicitly compute α 0 T on d-simplices and d-hypercubes. 2 e.g. [Egger 2008] Christian Waluga (AICES) HDG Methods October 6-8, / 26
9 Some remarks Hybridization: From DG to HDG Hybridization: HDG The hybridized method is consistent by construction! Standard error analysis yields optimal error estimates. Stability is independent of the particular choice of V h. More general boundary conditions are possible. Other problems were also investigated, e.g. Convection-diffusion [Egger and Schöberl 2009] Stokes problem [Cockburn et. al. 2010], [Egger and W. 2010b] Oseen problem (in preparation) Assembly in an element-wise fashion. Static condensation on element level. Upwinding can be easily incorporated. Locally varying polynomial degrees and nonconforming h-refinements possible. Christian Waluga (AICES) HDG Methods October 6-8, / 26
10 Hybridization: From DG to HDG HDG methods: Implementation (1/2) Implementation uses dune-pdelab and the core modules. Approximations in the interior of elements by monomials (using MonomLocalFiniteElementMap) Approximations on the element borders require some extra work: We use the IntersectionIndexSet provided in PDELab and a self-written FaceMonomLocalFiniteElementMap to define a grid function space on the faces of the mesh. Problem: The orientation of intersections may differ in two intersecting elements. This causes problems when mapping from intersection coordinates to coordinates inside the element. Christian Waluga (AICES) HDG Methods October 6-8, / 26
11 HDG methods (2/2) Hybridization: From DG to HDG The only solution we know so far is a little helper that finds the corresponding intersection in the outside element if it has a lower index or a higher level than the current intersection. We then map the coordinates as follows: if(wrongintersection) return rightintersection->geometryinoutside().global(x); else return rightintersection->geometryininside().global(x); The big disadvantage here is that we need to instantiate a quadratic number of intersections (performance issues?!) Christian Waluga (AICES) HDG Methods October 6-8, / 26
12 Hybridization: From DG to HDG Example: Oseen problem Numerical example: Kovaznay (1/3) Example: Oseen problem; Ω = ( 0.5, 2) ( 0.5, 1.5) [Kovasznay 1947] } u + w u + p = 0 on Ω, u = u div u = 0 exact on Ω Exact solution: ( 1 exp(λx) cos(2πy) u(x, y) = λ 2π exp(λx) sin(2πy) p(x, y) = 1 2 exp(2λx) + p, with parameters λ and p given as 8π 2 λ := ν 1 + ν π 2 ), and p = exp(2λx) dx. 2 Figure: Velocity vectors (ν = 1/10) Christian Waluga (AICES) HDG Methods October 6-8, / 26
13 Hybridization: From DG to HDG Example: Oseen problem Numerical example: Kovaznay (2/3) level P 1 P P 2 P P 3 P dimension L 2 error energy error K S error rate error rate K S error rate error rate K S error rate error rate Table: Kovsznay flow: Errors of the numerical solution for different inf-sup stable finite element approximations and a sequence of uniformly refined meshes. Christian Waluga (AICES) HDG Methods October 6-8, / 26
14 Hybridization: From DG to HDG Example: Oseen problem Numerical example: Kovaznay (3/3) y x Figure: Kovasznay flow (ν = 1): Streamlines and nonconforming mesh after 4 subsequent h-refinements. Christian Waluga (AICES) HDG Methods October 6-8, / 26
15 Outline Hybridized Mortar Methods 1 Hybridization: From DG to HDG Preliminaries A simple DG method Hybridization: HDG Example: Oseen problem 2 Hybridized Mortar Methods Relation between HDG and Hybridized Mortar Example: Stokes problem 3 Conclusion Christian Waluga (AICES) HDG Methods October 6-8, / 26
16 Hybridized Mortar Methods Hybridized mortar method (1/2) Relation between HDG and Hybridized Mortar Similarly as for the HDG method, we can define a hybridized mortar method 3 : V h := { v H 1 0 (Ω h ) : v T P p (T ), T T h } ˆV h := { ˆv L 2 (Γ h ) : v E P p (E), E T h (Γ h ) } Γ 13 Ω 1 Ω 2 Ω 3 Ω 4 Γ 34 partition: Ω h := {Ω 1, Ω 2,..., Ω N } interfaces: Γ ij := Ω i Ω j Γ h := {Γ ij } skeleton: Γ := Γ ij 3 cf. [Egger 2008] Christian Waluga (AICES) HDG Methods October 6-8, / 26
17 Hybridized Mortar Methods Hybridized mortar method (2/2) Relation between HDG and Hybridized Mortar Discrete problem: Find (u h, û h ) V h V h such that a h (u h, û h ; v h, ˆv h ) = (f, v h ) Ω, (v h, ˆv h ) V h V h. where we define Important results: a h (u, û, v, ˆv) :=( u, v) Ωh ( n u, v ˆv) Ωh ( n v, u û) Ωh + τ(u û, v ˆv) Ωh. No direct coupling between subdomain solutions. Elimination on subdomains leads to a Schur complement system for the hybrid variables only ( domain decomposition methods). Space for the hybrid variable can be chosen with great flexibility (no inf-sup-condition necessary for multiplier). Hybridized mortar methods for other problems were also analyzed, e.g. Maxwell [Hollaus et. al. 2010], Stokes [Egger and W. 2010a]. For the finest partition Ω h = T h, we recover the hybridized DG method. Christian Waluga (AICES) HDG Methods October 6-8, / 26
18 Hybridized Mortar Methods Hybrid mortar: Implementation In many applications, there exists a parametrization for the interface (e.g. planar, cylindrical, spherical,...) Idea: Lagrange multipliers live on d 1 dimensional structured meshes that are transformed to the physical space. Overlaps between multiplier mesh and subdomain meshes are computed using dune-grid-glue by C. Engwer and O. Sander. Christian Waluga (AICES) HDG Methods October 6-8, / 26
19 Hybridized Mortar Methods Hybrid mortar: Implementation In many applications, there exists a parametrization for the interface (e.g. planar, cylindrical, spherical,...) Idea: Lagrange multipliers live on d 1 dimensional structured meshes that are transformed to the physical space. Overlaps between multiplier mesh and subdomain meshes are computed using dune-grid-glue by C. Engwer and O. Sander. Christian Waluga (AICES) HDG Methods October 6-8, / 26
20 Hybridized Mortar Methods Hybrid mortar: Implementation In many applications, there exists a parametrization for the interface (e.g. planar, cylindrical, spherical,...) Idea: Lagrange multipliers live on d 1 dimensional structured meshes that are transformed to the physical space. Overlaps between multiplier mesh and subdomain meshes are computed using dune-grid-glue by C. Engwer and O. Sander. Christian Waluga (AICES) HDG Methods October 6-8, / 26
21 Hybridized Mortar Methods Hybrid mortar: Stokes problem Example: Stokes problem Stokes problem [Egger and W. 2010a]: Partitioning and triangulations: The subdomain interfaces are extracted using a Codim1Extractor The entire OneDGrid is extracted with a Codim0Extractor A transformation is given to place the interface grids on the coupling boundaries. We use the PSurface backend to generate remote intersections. Christian Waluga (AICES) HDG Methods October 6-8, / 26
22 Numerical example Hybridized Mortar Methods Example: Stokes problem Example: Stokes problem: Colliding flow, Ω = ( 1, 1) 2 u + p = 0 div u = 0 } on Ω, u = u exact on Ω u exact = ( 20xy 3, 5x 4 5y 4), p exact = 60x 2 y 20y 3 Plot of analytic solution; velocity vectors and pressure field Christian Waluga (AICES) HDG Methods October 6-8, / 26
23 Hybridized Mortar Methods Mortar: Numerical results (1/2) Example: Stokes problem Figure: Numerical results (p=2): h: 1.00; L 2 -error: 1.418; energy error: Christian Waluga (AICES) HDG Methods October 6-8, / 26
24 Hybridized Mortar Methods Mortar: Numerical results (1/2) Example: Stokes problem Figure: Numerical results (p=2): h: 0.50; L 2 -error: 0.171; energy error: Christian Waluga (AICES) HDG Methods October 6-8, / 26
25 Hybridized Mortar Methods Mortar: Numerical results (1/2) Example: Stokes problem Figure: Numerical results (p=2): h: 0.25; L 2 -error: 0.021; energy error: Christian Waluga (AICES) HDG Methods October 6-8, / 26
26 Hybridized Mortar Methods Mortar: Numerical results (2/2) Example: Stokes problem Figure: Simple example for a diffusion problem with non-matching meshes in 3D Christian Waluga (AICES) HDG Methods October 6-8, / 26
27 Outline Conclusion 1 Hybridization: From DG to HDG Preliminaries A simple DG method Hybridization: HDG Example: Oseen problem 2 Hybridized Mortar Methods Relation between HDG and Hybridized Mortar Example: Stokes problem 3 Conclusion Christian Waluga (AICES) HDG Methods October 6-8, / 26
28 Conclusion Conclusion Summary Presentation of hybridized DG methods. Implementation of HDG and hybridized mortar methods in DUNE. Possible applications include different interface problems (e.g. propellers). Known issues HDG methods in DUNE not (yet) naturally implementable? Future work Implementation of parallel codes for high performance computing. Analysis and implementation for time-dependent problems. Use of efficient domain decomposition solvers. Christian Waluga (AICES) HDG Methods October 6-8, / 26
29 Conclusion Conclusion Summary Presentation of hybridized DG methods. Implementation of HDG and hybridized mortar methods in DUNE. Possible applications include different interface problems (e.g. propellers). Known issues HDG methods in DUNE not (yet) naturally implementable? Future work Implementation of parallel codes for high performance computing. Analysis and implementation for time-dependent problems. Use of efficient domain decomposition solvers. Thanks for your attention! Christian Waluga (AICES) HDG Methods October 6-8, / 26
30 Acknowledgement Acknowledgement Financial support from the German Research Association (DFG) through grant GSC 111 is gratefully acknowledged. Christian Waluga (AICES) HDG Methods October 6-8, / 26
31 Literature Literature D. N. Arnold, F. Brezzi, B. Cockburn, and L. D. Marini. Unified Analysis of Discontinuous Galerkin Methods for Elliptic Problems. SIAM Journal for Numerical Analysis, 39: , R. Becker, P. Hansbo, and R. Stenberg. A finite element method for domain decomposition with non-matching grids. ESAIM: Mathematical Modeling and Numerical Analysis, 37(2): , B. Cockburn, J. Gopalakrishnan, and R. Lazarov. Unified hybridization of discontinuous Galerkin, mixed and continuous Galerkin methods for second order elliptic problems. SIAM Journal for Numerical Analysis, 47(2): , H. Egger. A class of hybrid mortar finite element methods for interface problems with non-matching meshes. Technical report, AICES , RWTH-Aachen University, H. Egger, and C. Waluga. A hybridized discontinuous Galerkin method for incompressible flow. Part I: Stokes problem. in preparation, H. Egger, and C. Waluga. A hybrid mortar finite element method for the Stokes problem. submitted, Christian Waluga (AICES) HDG Methods October 6-8, / 26
Analysis 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 informationAn Implementation of Hybrid Discontinuous Galerkin Methods in DUNE
Aachen Institute for Advanced Study in Computational Engineering Science Preprint: AICES-2010/12-1 15/December/2010 An Implementation of Hybrid Discontinuous Galerkin Methods in DUNE C. Waluga and H. Egger
More informationA Hybrid Discontinuous Galerkin Method 2 for Darcy-Stokes Problems 3 UNCORRECTED PROOF
1 A Hybrid Discontinuous Galerkin Method 2 for Darcy-Stokes Problems 3 Herbert Egger 1 and Christian Waluga 2 4 1 Center of Mathematics, Technische Universität München, Boltzmannstraße 3, 85748 5 Garching
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 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 informationA note on discontinuous Galerkin divergence-free solutions of the Navier-Stokes equations
A note on discontinuous Galerkin divergence-free solutions of the Navier-Stokes equations Bernardo Cockburn Guido anschat Dominik Schötzau June 1, 2007 Journal of Scientific Computing, Vol. 31, 2007, pp.
More informationA Two-grid Method for Coupled Free Flow with Porous Media Flow
A Two-grid Method for Coupled Free Flow with Porous Media Flow Prince Chidyagwai a and Béatrice Rivière a, a Department of Computational and Applied Mathematics, Rice University, 600 Main Street, Houston,
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 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 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 informationNon-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions
Non-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions Wayne McGee and Padmanabhan Seshaiyer Texas Tech University, Mathematics and Statistics (padhu@math.ttu.edu) Summary. In
More informationHybrid Discontinuous Galerkin methods for incompressible flow problems
Hybrid Discontinuous Galerkin methods incompressible flow problems Christoph Lehrenfeld, Joachim Schöberl MathCCES Computational Mathematics in Engineering Workshop Linz, May 31 - June 1, 2010 Contents
More informationAn Equal-order DG Method for the Incompressible Navier-Stokes Equations
An Equal-order DG Method for the Incompressible Navier-Stokes Equations Bernardo Cockburn Guido anschat Dominik Schötzau Journal of Scientific Computing, vol. 40, pp. 188 10, 009 Abstract We introduce
More informationA LAGRANGE MULTIPLIER METHOD FOR ELLIPTIC INTERFACE PROBLEMS USING NON-MATCHING MESHES
A LAGRANGE MULTIPLIER METHOD FOR ELLIPTIC INTERFACE PROBLEMS USING NON-MATCHING MESHES P. HANSBO Department of Applied Mechanics, Chalmers University of Technology, S-4 96 Göteborg, Sweden E-mail: hansbo@solid.chalmers.se
More informationA note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations
A note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations S. Hussain, F. Schieweck, S. Turek Abstract In this note, we extend our recent work for
More informationAn A Posteriori Error Estimate for Discontinuous Galerkin Methods
An A Posteriori Error Estimate for Discontinuous Galerkin Methods Mats G Larson mgl@math.chalmers.se Chalmers Finite Element Center Mats G Larson Chalmers Finite Element Center p.1 Outline We present an
More informationOn an Approximation Result for Piecewise Polynomial Functions. O. Karakashian
BULLETIN OF THE GREE MATHEMATICAL SOCIETY Volume 57, 010 (1 7) On an Approximation Result for Piecewise Polynomial Functions O. arakashian Abstract We provide a new approach for proving approximation results
More informationHp-Adaptivity on Anisotropic Meshes for Hybridized Discontinuous Galerkin Scheme
Hp-Adaptivity on Anisotropic Meshes for Hybridized Discontinuous Galerkin Scheme Aravind Balan, Michael Woopen and Georg May AICES Graduate School, RWTH Aachen University, Germany 22nd AIAA Computational
More informationSTABILIZED DISCONTINUOUS FINITE ELEMENT APPROXIMATIONS FOR STOKES EQUATIONS
STABILIZED DISCONTINUOUS FINITE ELEMENT APPROXIMATIONS FOR STOKES EQUATIONS RAYTCHO LAZAROV AND XIU YE Abstract. In this paper, we derive two stabilized discontinuous finite element formulations, symmetric
More informationWEAK GALERKIN FINITE ELEMENT METHODS ON POLYTOPAL MESHES
INERNAIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 12, Number 1, Pages 31 53 c 2015 Institute for Scientific Computing and Information WEAK GALERKIN FINIE ELEMEN MEHODS ON POLYOPAL MESHES LIN
More informationASM-BDDC Preconditioners with variable polynomial degree for CG- and DG-SEM
ASM-BDDC Preconditioners with variable polynomial degree for CG- and DG-SEM C. Canuto 1, L. F. Pavarino 2, and A. B. Pieri 3 1 Introduction Discontinuous Galerkin (DG) methods for partial differential
More informationANALYSIS OF AN INTERFACE STABILISED FINITE ELEMENT METHOD: THE ADVECTION-DIFFUSION-REACTION EQUATION
ANALYSIS OF AN INTERFACE STABILISED FINITE ELEMENT METHOD: THE ADVECTION-DIFFUSION-REACTION EQUATION GARTH N. WELLS Abstract. Analysis of an interface stabilised finite element method for the scalar advectiondiffusion-reaction
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 informationThe Discontinuous Galerkin Finite Element Method
The Discontinuous Galerkin Finite Element Method Michael A. Saum msaum@math.utk.edu Department of Mathematics University of Tennessee, Knoxville The Discontinuous Galerkin Finite Element Method p.1/41
More 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 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 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 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 informationEnergy norm a-posteriori error estimation for divergence-free discontinuous Galerkin approximations of the Navier-Stokes equations
INTRNATIONAL JOURNAL FOR NUMRICAL MTHODS IN FLUIDS Int. J. Numer. Meth. Fluids 19007; 1:1 [Version: 00/09/18 v1.01] nergy norm a-posteriori error estimation for divergence-free discontinuous Galerkin approximations
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 informationFourier Type Error Analysis of the Direct Discontinuous Galerkin Method and Its Variations for Diffusion Equations
J Sci Comput (0) 5:68 655 DOI 0.007/s095-0-9564-5 Fourier Type Error Analysis of the Direct Discontinuous Galerkin Method and Its Variations for Diffusion Equations Mengping Zhang Jue Yan Received: 8 September
More informationOn a Discontinuous Galerkin Method for Surface PDEs
On a Discontinuous Galerkin Method for Surface PDEs Pravin Madhavan (joint work with Andreas Dedner and Bjo rn Stinner) Mathematics and Statistics Centre for Doctoral Training University of Warwick Applied
More informationFEM: Domain Decomposition and Homogenization for Maxwell s Equations of Large Scale Problems
FEM: and Homogenization for Maxwell s Equations of Large Scale Problems Karl Hollaus Vienna University of Technology, Austria Institute for Analysis and Scientific Computing February 13, 2012 Outline 1
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 informationA High-Order Discontinuous Galerkin Method for the Unsteady Incompressible Navier-Stokes Equations
A High-Order Discontinuous Galerkin Method for the Unsteady Incompressible Navier-Stokes Equations Khosro Shahbazi 1, Paul F. Fischer 2 and C. Ross Ethier 1 1 University of Toronto and 2 Argonne National
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 information1. Introduction. We consider the model problem: seeking an unknown function u satisfying
A DISCONTINUOUS LEAST-SQUARES FINITE ELEMENT METHOD FOR SECOND ORDER ELLIPTIC EQUATIONS XIU YE AND SHANGYOU ZHANG Abstract In tis paper, a discontinuous least-squares (DLS) finite element metod is introduced
More informationTrefftz-discontinuous Galerkin methods for time-harmonic wave problems
Trefftz-discontinuous Galerkin methods for time-harmonic wave problems Ilaria Perugia Dipartimento di Matematica - Università di Pavia (Italy) http://www-dimat.unipv.it/perugia Joint work with Ralf Hiptmair,
More informationDiscontinuous Galerkin Methods: Theory, Computation and Applications
Discontinuous Galerkin Methods: Theory, Computation and Applications Paola. Antonietti MOX, Dipartimento di Matematica Politecnico di Milano.MO. X MODELLISTICA E CALCOLO SCIENTIICO. MODELING AND SCIENTIIC
More informationC e n t r u m v o o r W i s k u n d e e n I n f o r m a t i c a
C e n t r u m v o o r W i s k u n d e e n I n f o r m a t i c a Modelling, Analysis and Simulation Modelling, Analysis and Simulation Bilinear forms for the recovery-based discontinuous Galerkin method
More informationALMOST OPTIMAL INTERIOR PENALTY DISCONTINUOUS APPROXIMATIONS OF SYMMETRIC ELLIPTIC PROBLEMS ON NON-MATCHING GRIDS
ALMOST OPTIMAL INTERIOR PENALTY DISCONTINUOUS APPROXIMATIONS OF SYMMETRIC ELLIPTIC PROBLEMS ON NON-MATCHING GRIDS R.D. LAZAROV, J.E. PASCIAK, J. SCHÖBERL, AND P.S. VASSILEVSKI Abstract. We consider an
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 informationTwo-scale Dirichlet-Neumann preconditioners for boundary refinements
Two-scale Dirichlet-Neumann preconditioners for boundary refinements Patrice Hauret 1 and Patrick Le Tallec 2 1 Graduate Aeronautical Laboratories, MS 25-45, California Institute of Technology Pasadena,
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 informationA PROJECTION-BASED ERROR ANALYSIS OF HDG METHODS BERNARDO COCKBURN, JAYADEEP GOPALAKRISHNAN, AND FRANCISCO-JAVIER SAYAS
A PROJECTION-BASED ERROR ANALYSIS OF HDG METHODS BERNARDO COCBURN, JAYADEEP GOPALARISHNAN, AND FRANCISCO-JAVIER SAYAS Abstract. We introduce a new technique for the error analysis of hybridizable discontinuous
More informationICES REPORT Analysis of the DPG Method for the Poisson Equation
ICES REPORT 10-37 September 2010 Analysis of the DPG Method for the Poisson Equation by L. Demkowicz and J. Gopalakrishnan The Institute for Computational Engineering and Sciences The University of Texas
More informationThe CG1-DG2 method for conservation laws
for conservation laws Melanie Bittl 1, Dmitri Kuzmin 1, Roland Becker 2 MoST 2014, Germany 1 Dortmund University of Technology, Germany, 2 University of Pau, France CG1-DG2 Method - Motivation hp-adaptivity
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 informationDiscontinuous Galerkin Methods: Theory, Computation and Applications
Discontinuous Galerkin Methods: Theory, Computation and Applications Paola. Antonietti MOX, Dipartimento di Matematica Politecnico di Milano.MO. X MODELLISTICA E CALCOLO SCIENTIICO. MODELING AND SCIENTIIC
More informationDiscontinuous Galerkin Method for interface problem of coupling different order elliptic equations
Discontinuous Galerkin Method for interface problem of coupling different order elliptic equations Igor Mozolevski, Endre Süli Federal University of Santa Catarina, Brazil Oxford University Computing Laboratory,
More information1. Introduction. The Stokes problem seeks unknown functions u and p satisfying
A DISCRETE DIVERGENCE FREE WEAK GALERKIN FINITE ELEMENT METHOD FOR THE STOKES EQUATIONS LIN MU, JUNPING WANG, AND XIU YE Abstract. A discrete divergence free weak Galerkin finite element method is developed
More informationAN EQUILIBRATED A POSTERIORI ERROR ESTIMATOR FOR THE INTERIOR PENALTY DISCONTINUOUS GALERKIN METHOD
AN EQUILIBRATED A POSTERIORI ERROR ESTIMATOR FOR THE INTERIOR PENALTY DISCONTINUOUS GALERIN METHOD D. BRAESS, T. FRAUNHOLZ, AND R. H. W. HOPPE Abstract. Interior Penalty Discontinuous Galerkin (IPDG) methods
More informationVirtual Element Methods for general second order elliptic problems
Virtual Element Methods for general second order elliptic problems Alessandro Russo Department of Mathematics and Applications University of Milano-Bicocca, Milan, Italy and IMATI-CNR, Pavia, Italy Workshop
More informationAn Extended Finite Element Method for a Two-Phase Stokes problem
XFEM project An Extended Finite Element Method for a Two-Phase Stokes problem P. Lederer, C. Pfeiler, C. Wintersteiger Advisor: Dr. C. Lehrenfeld August 5, 2015 Contents 1 Problem description 2 1.1 Physics.........................................
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 informationPREPRINT 2010:25. Fictitious domain finite element methods using cut elements: II. A stabilized Nitsche method ERIK BURMAN PETER HANSBO
PREPRINT 2010:25 Fictitious domain finite element methods using cut elements: II. A stabilized Nitsche method ERIK BURMAN PETER HANSBO Department of Mathematical Sciences Division of Mathematics CHALMERS
More informationThe Mortar Boundary Element Method
The Mortar Boundary Element Method A Thesis submitted for the degree of Doctor of Philosophy by Martin Healey School of Information Systems, Computing and Mathematics Brunel University March 2010 Abstract
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 class of finite element methods: weak Galerkin methods
New class of finite element methods: weak Galerkin methods Xiu Ye University of Arkansas at Little Rock Second order elliptic equation Consider second order elliptic problem: a u = f, in Ω (1) u = 0, on
More informationarxiv: v1 [math.na] 15 Nov 2017
Noname manuscript No. (will be inserted by the editor) An HDG method with orthogonal projections in facet integrals Issei Oikawa arxiv:1711.05396v1 [math.na] 15 Nov 2017 Received: date / Accepted: date
More informationEfficient Augmented Lagrangian-type Preconditioning for the Oseen Problem using Grad-Div Stabilization
Efficient Augmented Lagrangian-type Preconditioning for the Oseen Problem using Grad-Div Stabilization Timo Heister, Texas A&M University 2013-02-28 SIAM CSE 2 Setting Stationary, incompressible flow problems
More informationBUBBLE STABILIZED DISCONTINUOUS GALERKIN METHOD FOR PARABOLIC AND ELLIPTIC PROBLEMS
BUBBLE STABILIZED DISCONTINUOUS GALERKIN METHOD FOR PARABOLIC AND ELLIPTIC PROBLEMS ERIK BURMAN AND BENJAMIN STAMM Abstract. In this paper we give an analysis of a bubble stabilized discontinuous Galerkin
More informationBUBBLE STABILIZED DISCONTINUOUS GALERKIN METHOD FOR STOKES PROBLEM
BUBBLE STABILIZED DISCONTINUOUS GALERKIN METHOD FOR STOKES PROBLEM ERIK BURMAN AND BENJAMIN STAMM Abstract. We propose a low order discontinuous Galerkin method for incompressible flows. Stability of the
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 informationAnalysis of a DG XFEM Discretization for a Class of Two Phase Mass Transport Problems
Analysis of a DG XFEM Discretization for a Class of Two Phase Mass Transport Problems Christoph Lehrenfeld and Arnold Reusken Bericht Nr. 340 April 2012 Key words: transport problem, Nitsche method, XFEM,
More informationFROM GODUNOV TO A UNIFIED HYBRIDIZED DISCONTINUOUS GALERKIN FRAMEWORK FOR PARTIAL DIFFERENTIAL EQUATIONS
FROM GODUNOV TO A UNIFIED HYBRIDIZED DISCONTINUOUS GALERKIN FRAMEWORK FOR PARTIAL DIFFERENTIAL EQUATIONS TAN BUI-THANH Abstract. By revisiting the basic Godunov approach for linear system of hyperbolic
More informationPARTITION OF UNITY FOR THE STOKES PROBLEM ON NONMATCHING GRIDS
PARTITION OF UNITY FOR THE STOES PROBLEM ON NONMATCHING GRIDS CONSTANTIN BACUTA AND JINCHAO XU Abstract. We consider the Stokes Problem on a plane polygonal domain Ω R 2. We propose a finite element method
More informationNumerical Simulation of Flows in Highly Heterogeneous Porous Media
Numerical Simulation of Flows in Highly Heterogeneous Porous Media R. Lazarov, Y. Efendiev, J. Galvis, K. Shi, J. Willems The Second International Conference on Engineering and Computational Mathematics
More informationElectronic Transactions on Numerical Analysis Volume 32, 2008
Electronic Transactions on Numerical Analysis Volume 32, 2008 Contents 1 On the role of boundary conditions for CIP stabilization of higher order finite elements. Friedhelm Schieweck. We investigate the
More informationThe Virtual Element Method: an introduction with focus on fluid flows
The Virtual Element Method: an introduction with focus on fluid flows L. Beirão da Veiga Dipartimento di Matematica e Applicazioni Università di Milano-Bicocca VEM people (or group representatives): P.
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 informationICES REPORT Global properties of DPG test spaces for convection-diffusion problems
ICES REPORT 13-05 February 2013 Global properties of DPG test spaces for convection-diffusion problems by Jesse Chan, Jay Gopalakrishnan, and Leszek Demkowicz The Institute for Computational Engineering
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 informationLOCAL MASS-CONSERVATION FOR STABILIZED LINEAR FINITE ELEMENT INTERPOLATIONS OF INCOMPRESSIBLE FLOW
LOCAL MASS-CONSERVATION FOR STABILIZED LINEAR FINITE ELEMENT INTERPOLATIONS OF INCOMPRESSIBLE FLOW B. GMEINER, C. WALUGA, AND B. WOHLMUTH (a) (b) Abstract. In this work, we discuss a class of stabilized
More informationNUMERICAL STUDIES OF VARIATIONAL-TYPE TIME-DISCRETIZATION TECHNIQUES FOR TRANSIENT OSEEN PROBLEM
Proceedings of ALGORITMY 212 pp. 44 415 NUMERICAL STUDIES OF VARIATIONAL-TYPE TIME-DISCRETIZATION TECHNIQUES FOR TRANSIENT OSEEN PROBLEM NAVEED AHMED AND GUNAR MATTHIES Abstract. In this paper, we combine
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 informationMaximum norm estimates for energy-corrected finite element method
Maximum norm estimates for energy-corrected finite element method Piotr Swierczynski 1 and Barbara Wohlmuth 1 Technical University of Munich, Institute for Numerical Mathematics, piotr.swierczynski@ma.tum.de,
More 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 information20. A Dual-Primal FETI Method for solving Stokes/Navier-Stokes Equations
Fourteenth International Conference on Domain Decomposition Methods Editors: Ismael Herrera, David E. Keyes, Olof B. Widlund, Robert Yates c 23 DDM.org 2. A Dual-Primal FEI Method for solving Stokes/Navier-Stokes
More informationOn some numerical convergence studies of mixed finite element methods for flow in porous media
On some numerical convergence studies of mixed finite element methods for flow in porous media Gergina Pencheva Abstract We consider an expanded mixed finite element method for solving second-order elliptic
More informationBasics and some applications of the mortar element method
GAMM-Mitt. 28, No. 2, 97 123 (2005) Basics and some applications of the mortar element method Christine Bernardi 1, Yvon Maday 1, and Francesca Rapetti 2 1 Laboratoire Jacques-Louis Lions, C.N.R.S. & université
More informationINTRODUCTION TO FINITE ELEMENT METHODS
INTRODUCTION TO FINITE ELEMENT METHODS LONG CHEN Finite element methods are based on the variational formulation of partial differential equations which only need to compute the gradient of a function.
More informationA posteriori error estimator based on gradient recovery by averaging for discontinuous Galerkin methods
A posteriori error estimator based on gradient recovery by averaging for discontinuous Galerkin methods Emmanuel Creusé, Serge Nicaise October 7, 29 Abstract We consider some (anisotropic and piecewise
More informationADAPTIVE HYBRIDIZED INTERIOR PENALTY DISCONTINUOUS GALERKIN METHODS FOR H(CURL)-ELLIPTIC PROBLEMS
ADAPTIVE HYBRIDIZED INTERIOR PENALTY DISCONTINUOUS GALERKIN METHODS OR HCURL-ELLIPTIC PROBLEMS C. CARSTENSEN, R. H. W. HOPPE, N. SHARMA, AND T. WARBURTON Abstract. We develop and analyze an adaptive hybridized
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 super convergent hybridisable discontinuous Galerkin method for linear elasticity
A super convergent hybridisable discontinuous Galerkin method for linear elasticity arxiv:82.223v [math.na] 6 Feb 28 Ruben Sevilla Zienkiewicz Centre for Computational Engineering, College of Engineering,
More informationNitsche XFEM with Streamline Diffusion Stabilization for a Two Phase Mass Transport Problem
Nitsche XFEM with Streamline Diffusion Stabilization for a Two Phase Mass Transport Problem Christoph Lehrenfeld and Arnold Reusken Bericht Nr. 333 November 2011 Key words: transport problem, Nitsche method,
More informationFEniCS Course. Lecture 0: Introduction to FEM. Contributors Anders Logg, Kent-Andre Mardal
FEniCS Course Lecture 0: Introduction to FEM Contributors Anders Logg, Kent-Andre Mardal 1 / 46 What is FEM? The finite element method is a framework and a recipe for discretization of mathematical problems
More informationA Nonoverlapping Subdomain Algorithm with Lagrange Multipliers and its Object Oriented Implementation for Interface Problems
Contemporary Mathematics Volume 8, 998 B 0-88-0988--03030- A Nonoverlapping Subdomain Algorithm with Lagrange Multipliers and its Object Oriented Implementation for Interface Problems Daoqi Yang. Introduction
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 informationc 2008 Society for Industrial and Applied Mathematics
SIAM J. NUMER. ANAL. Vol. 46, No. 3, pp. 640 65 c 2008 Society for Industrial and Applied Mathematics A WEIGHTED H(div) LEAST-SQUARES METHOD FOR SECOND-ORDER ELLIPTIC PROBLEMS Z. CAI AND C. R. WESTPHAL
More informationAn Iterative Substructuring Method for Mortar Nonconforming Discretization of a Fourth-Order Elliptic Problem in two dimensions
An Iterative Substructuring Method for Mortar Nonconforming Discretization of a Fourth-Order Elliptic Problem in two dimensions Leszek Marcinkowski Department of Mathematics, Warsaw University, Banacha
More informationHigh Order Unfitted Finite Element Methods for Interface Problems and PDEs on Surfaces
N O V E M B E R 2 0 1 6 P R E P R I N T 4 5 9 High Order Unfitted Finite Element Methods for Interface Problems and PDEs on Surfaces Christoph Lehrenfeld and Arnold Reusken Institut für Geometrie und Praktische
More informationA NEW FAMILY OF STABLE MIXED FINITE ELEMENTS FOR THE 3D STOKES EQUATIONS
MATHEMATICS OF COMPUTATION Volume 74, Number 250, Pages 543 554 S 0025-5718(04)01711-9 Article electronically published on August 31, 2004 A NEW FAMILY OF STABLE MIXED FINITE ELEMENTS FOR THE 3D STOKES
More information1. Introduction. We consider the model problem that seeks an unknown function u = u(x) satisfying
A SIMPLE FINITE ELEMENT METHOD FOR LINEAR HYPERBOLIC PROBLEMS LIN MU AND XIU YE Abstract. In this paper, we introduce a simple finite element method for solving first order hyperbolic equations with easy
More informationCoupling of HDG with a double-layer potential BEM
Coupling of HDG with a double-layer potential BEM Zhixing Fu, Norbert Heuer, & Francisco Javier Sayas June 29, 2018 arxiv:1212.2128v1 [math.na] 10 Dec 2012 Abstract In this paper we propose and analyze
More informationAdditive Average Schwarz Method for a Crouzeix-Raviart Finite Volume Element Discretization of Elliptic Problems
Additive Average Schwarz Method for a Crouzeix-Raviart Finite Volume Element Discretization of Elliptic Problems Atle Loneland 1, Leszek Marcinkowski 2, and Talal Rahman 3 1 Introduction In this paper
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 informationDiscontinuous Galerkin Methods: An Overview and Some Applications. Daya Reddy UNIVERSITY OF CAPE TOWN
Discontinuous Galerkin Methods: An Overview and Some Applications Daya Reddy UNIVRSITY OF CAP TOWN Joint work with Beverley Grieshaber and Andrew McBride SANUM, Stellenbosch, 22 24 March 2016 Structure
More informationETNA Kent State University
Electronic Transactions on Numerical Analysis. Volume 37, pp. 166-172, 2010. Copyright 2010,. ISSN 1068-9613. ETNA A GRADIENT RECOVERY OPERATOR BASED ON AN OBLIQUE PROJECTION BISHNU P. LAMICHHANE Abstract.
More information