Discontinuous Galerkin Methods: An Overview and Some Applications. Daya Reddy UNIVERSITY OF CAPE TOWN
|
|
- Audrey Holland
- 6 years ago
- Views:
Transcription
1 Discontinuous Galerkin Methods: An Overview and Some Applications Daya Reddy UNIVRSITY OF CAP TOWN Joint work with Beverley Grieshaber and Andrew McBride SANUM, Stellenbosch, March 2016
2 Structure of talk A model elliptic problem: weak or variational formulations The Galerkin finite element method: analysis and approximations Discontinuous Galerkin (DG) formulations Near-incompressibility in elasticity
3 Model problem: deformation of a membrane Minimization of an energy min J(v) J(v) = 1 v 2 rv 2 dx fvdx f The solution u satisfies the weak problem ru rvdx= fvdx for all functions v which satisfy v = 0 on the boundary Sufficiently smooth u satisfies the Poisson equation and boundary condition u = f on u = 0 on u 2 2
4 The membrane problem, continued min v2v J(v) := 1 2 ru rv dx= rv 2 dx fv dx fv dx Define the bilinear form a(, ) and linear functional `( ) a : V V! R, a(u, v) = ru rv dx ` : V! R, `(v) = fv dx Thus the above problem is min v2v 1 2a(v, v) `(v) or equivalently a(u, v) =`(v) 8v 2 V
5 Interlude: the Sobolev spaces H m () Built from the Lebesgue space of square-integrable functions: L 2 () = v : v 2 dx := kvk 2 0 < 1 Define, for integer m 0, H m () = Seminorm v : D v 2 L 2 (), apple m v 2 m = X D v 2 dx =m D v = for problem in R 2 Hilbert space with induced norm kvk 2 m = X e.g. kvk 2 1 = applem v 2 m h v i dx We will also need H 1 0 () ={v 2 H 1 () : v = 0 on }
6 Well-posedness of the variational problem min w2w 1 2a(w, w) `(w) This problem has a unique solution if: W is a closed subspace of a Hilbert space H a is coercive or W-elliptic: a(w, w) kwk 2 H a is continuous: a(w, z) apple Mkwk H kzk H ` is continuous: `(z) apple Ckzk H The model problem has a unique solution in H 1 0 ()
7 Finite element approximations Aim: to pose the variational problem on a finite-dimensional subspace V h V 1. Partition the domain into subdomains or finite elements 2. Construct a basis {' i } N i=1 for V h comprising continuous functions that are polynomials on each element portion of hip replacement: physical object and finite element model
8 The Galerkin finite element method 3. The piecewise-polynomial approximations can be written u h = X i ' i (x)u i 'u ' i v h = X i ' i (x)v i 'v 4. Substitute in the weak formulation a(u h,v h )=`(v h ) to obtain X X v i [a(' i, ' j )u i `(' j )] = 0 a(' i, ' j ) u i = `(' j ) {z } {z } i,j i K ji F j Ku = F
9 Convergence of finite element approximations Construct V h V and seek u h 2 V h such that for all v h 2 V h a(u h,v h )=`(v h ) for all v h 2 V h Ku = F h T = diameter of T mesh size h = max T 2T h T h T Define the error by sense that u u h : under what conditions do we have convergence in the lim u h = u? h!0 Orthogonality of the error: a(u u h,v h ) = a(u, v h ) a(u h,v h ) = `(v h ) `(v h ) = 0
10 An a priori estimate ku u h k 2 V apple a(u u h,u u h ) = a(u u h,u v h )+a(u u h,v h u h ) (V- ellipticity) apple a(u u h,u v h ) apple Mku u h k V ku v h k V (continuity) (orthogonality of error) ku u h k 1 apple C inf v h 2V h ku v h k 1 Céa s lemma Strategy for obtaining error bound: a) choose v h to be the interpolate of uin V h b) use interpolation error estimate to bound actual error
11 Finite element interpolation theory Ciarlet and Raviart Arch. Rat. Mech. Anal h T = diameter of T h T T =sup{diameter of B; B a ball contained in T } T T = h T / T Let T be a triangulation of a bounded domain with polygonal boundary: = [ T 2T T Define the mesh size h = max T 2T h T A family of triangulations is regular as h! 0 if there exists > 0 such that T apple for all T 2 T h
12 Finite element interpolation theory (Local estimate) For a regular triangulation with v 2 H k+1 (T ), k +1 m and the interpolation operator which maps functions to polynomials of degree apple k, v v m,t apple Ch k+1 m v k+1,t (Global estimate) Let T domain. Define be a uniformly regular triangulation of a polygonal V h = {v h 2 C( ) :v h T 2 P k } \ V h : H 2 ()! V h h = max T h T (global interpolator) h v v h T kv h vk m, apple Ch k+1 m v k+1, m =0, 1
13 Convergence of finite element approximations ku u h k V apple C inf vh 2V h ku v h k V apple Cku h uk V apple Ch min(k,r 1) u r So for the simplest approximation, by piecewise-linear simplices, ku u h k 1 apple Ch u 2 for the second-order elliptic equations Could use piecewise-quadratic simplices (k = 2) in which case ku u h k 1 = O(h 2 ) -- provided that the solution is smooth enough to belong to H 3 ()!
14 Discontinuous Galerkin (DG) methods
15 Discontinuous Galerkin (DG) methods Drop the condition of continuity So V h 62 V A member v h 2 V h is now a polynomial on each element, but not continuous across element boundaries An immediate consequence: much greater number of degrees of freedom u { u} N 1 K 1 K 2 N 2 e 12
16 What do solutions look like? 4.3 Implementation of the predictor corrector solution strategy (a) Conforming approximation (b) Discontinuous Galerkin approximation (c) Discontinuous Galerkin approximation with a 101
17 Why bother with DG? Can accommodate hanging nodes Useful in adaptive mesh refinement log(energy error) log(interface jump)
18 fficiency and accuracy = = = 10 Conforming DG Displacement rror (m) 1 ten yck and Lew 2010 (nonlinear elas7city) Total DOF Figure 6. Plot of the L 2 -norm over B 0 of the error in the displacement field as a function of the total number of degrees of freedom for the test in Section 6.1. Curves are shown for the discontinuous Galerkin method as well as for the conforming one. Remarkably, the two curves overlap, indicating that both methods provide the same accuracy for the same computational cost, i.e. for this example they are equally efficient.
19 Our objective For the problem of deformations of elastic bodies DG methods show good behaviour for near-incompressibility with the use of low-order triangles in two dimensions, and tetrahedra in three DG with quadrilaterals and hexahedra less straightforward: poor behaviour, including locking, for low-order approximations We show why this is so, and propose some remedies
20 First, review derivation of heat equation Heat equation: Balance of energy + div q = s q = s = heat flux vector heat source Fourier s law q = kr# give the (steady) heat equation k # = s
21 Governing equations for elasticity displacement vector u x u(x) stress tensor or matrix σ lasticity: quilibrium + div = f give Navier s equation Hooke s law = div u +(ru +[ru] T ) u = [ u +(1+ )rdiv u] =f
22 quivalent to minimization problem min u J(u) = 1 2 r s u 2 +(1+ )(div u) 2 dv f u dv Define a(u, v) = [r s u : r s v +(1+ )div u div v] dx `(v) = f v dx r s u = 1 2 [ru +(ru)t ] (r s u) ij = j i Then the minimization problem is or equivalently min v2v 1 2a(v, v) `(v) a(u, v) =`(v) which has a unique solution in V := [H 1 0 ()] d Furthermore (Brenner and Sung 1992) the solution satisfies u 2 [H 2 ()] 2 kuk H 2 + k div uk H 1 apple Ckfk L 2
23 Locking The use of low-order elements leads to a non-physical solution in the incompressible limit a(u, v) = [r s u : r s v +(1+ )div u div v] dx!1 div u! 0 Q 1 0,8 0,7 0,6 0,5 0,4 Case III 0,3 Case V (AS) 0,2 MHW Case 1 0,1 HW Case I or Q1 0-0,6-0,4-0,2 0-0,1 0,2 0,4 0,6-0,2 We explore the use of DG methods as a remedy for locking Q 1
24 DG formulation: some definitions Jumps and averages: on an interior edge of element T, n + T + T n [v] = v + n + + v n, [[v]] = v + n + + v n, {v} = 1 2 (v+ + v ) [[ ]] = + n + + n, { } = 1 2 ( + + ) For example [[ ]] =( + )n + On an edge that forms part of the boundary, [[v]] = v n { } = n
25 Setting it up take the inner product with a test function integrate, and integrate by parts div v dv = f v dv u = [ u +(1+ {z )rdiv u] } div V = div u +(ru +[ru] T ) = f Consider the left-hand side: div v dv = X div T T v dv div v dv = X n v ds X T T : r s v dv
26 Thus weak equation is now X n v ds + X T T : r s v dv = f v Next, we need the magic formula T X v n ds = X [[v]] : { } ds + X int int {v} [[ ]] ds
27 This gives X [[v]] : { } ds X int int {v} [[ ]] ds + : r s v dv = f v dv Since the exact solution is smooth ( u 2 [H 2 ()] d ) and the stress is continuous we can assume that [[ ]] = 0 which leaves X [[v]] : { } ds + : r s v dv = f v dv
28 xploiting again the smoothness of the exact solution we can add a term to symmetrize the problem: X [[v]] : { (u)} ds X [[u]] : { (v)} ds + X [[u]] = 0 : r s v dv = f v dv Finally, we may need to stabilize the problem: this gives us the symmetric interior penalty (SIP) formulation (Douglas and Dupont 1976) X [[v]] : { (u)} ds X : { (v)} ds + k [[u]] X + X : r s v dv = 1 h [[u]][[v]] f v dv
29 The DG formulation V h = v h : v h 2 [L 2 ()] d, v h T 2 P 1 (T ) or Q 1 (T ) P 1 (T ): v h = a 0 + a 1 x + a 2 y Q 1 (T ): v h = a 0 + a 1 x + a 2 y + a 3 xy
30 V h = v h : v h 2 [L 2 ()] d, v h T 2 P 1 (T ) or Q 1 (T ) a h (u h, v h )=`(v h ) X T T (u) :r s v dv + X a h (u, v) [[u]] : { (v)} ds X +k X [[v]] : { (u)} ds 1 h [[u]][[v]] = f v dv = 8 >< >: Nonsymmetric Interior Penalty Galerkin (NIPG) RIVIÈR & WHLR 1999, 2000; ODN, BABUSKA & BAUMANN 1998 Symmetric Interior Penalty Galerkin (SIPG) DOUGLAS & DUPONT 1976, ARNOLD 1982, HANSBO & LARSON 2002 Incomplete Interior Penalty Galerkin (IIPG) DAWSON, SUN AND WHLR 2004
31 DG with triangles is uniformly convergent WIHLR 2002 y 62 Chapter 3. Locking-Freeh DGFM for lasticity Problems DGFM on graded mesh λ=1 λ=100 λ=500 λ=1000 λ= ements O x absolute L 2 error Figure 3.2: Graded mesh with refinement towards the origin (γ = (1/2, 0, 0, 0, 0)) number of degrees of freedom Figure 3.3: Uniform mesh (i.e. graded mesh with γ = (0, 0, 0, 0, 0)). Figure 3.5: Performance of the DGFM on the L-sh mesh) and with γ = 0 (uniform mesh). in order to obtain the optimal convergence rate, a graded mesh with refinement towards the origin must be used for the numerical simulations. The first picture of Figure 3.4 shows the errors of the DGFM for in the energy norm u 2 h = K T λ {1, 100, 500, 1000, 5000} (µ = 1) ϵ(u) 2 L 2 (K ) + 1 m elast e int,d e 1 e [u] 2 ds on a graded mesh with grading vector γ = (1/2, 0, 0, 0, 0) (cf. Figure 3.2). Ob-
32 Some computational results with quads appear to show locking-free behaviour OBB SIPG NIPG IIPG Fig. 8. Tensile stress contours ðr y > 0:15 MPaÞ of DG solutions of beam problem with Poisson s ratio (a) OBB; (b) SIPG; d p ¼ 30; (c) NIPG; d p ¼ 10; and (d) IIPG; d p ¼ 50. LIU, WHLR AND DAWSON 2009 These authors report good results using all methods
33 but not so in other cases Grieshaber, McBride and R. (2015) ALL ALL
34 We know that DG works for simplicial elements Some authors report good results for quads We find locking for a simple example What s going on?
35 Let s see why DG with triangles works a h (u, v) := X T 2T h +k µ µ X 2 T (u) :r s v + X 2 1 h [[u]] : [[v]] + k [[u]] : { (v)} X 1 h 2 X 2 [[u]] : [[v]] { (u)} :[[v]] := 2µ + " X " T 2T h + X T 2T h T r s u : r s v dx + X 2 T (div u)(div v)+ X 2 [[u]] : {r s v} (note the use of two stabilization terms) X 2 (I) (II) (III) 2 [[u]] : {(div v)i} X {ru} :[[v]] + k µ X 2 {(div u)i} :[[v]] ds +k X 2 1 h # 1 [[u]] : [[v]] 2h [[u]] : [[v]] (IV) #
36 rror analysis for linear triangles: NIPG Approximation error e = u u h u = interpolant of u = u u {z } interpolation error + u u h {z } Use Crouzeix-Raviart interpolant: linear on triangles, continuous at midpoints of edges Properties: (u u) n ds =0 T div(u u) dv =0
37 e = + = u u = u u h 2 V h ) div, ( ) const. λ-terms in error bound can be handled as follows: for example, (I) = X T to give eventually T (div )(div ) ds = X T div div ds =0 kek 2 DG apple Ch 2 kuk 2 H 2 () + 2 k div uk 2 H 1 () T Use estimate kuk H 2 + k div uk H 1 apple Ckfk L 2 to get uniform convergence
38 Convergence analysis for quadrilaterals First need to construct a suitable interpolant u 2 V h As before e = u u h = u u {z } + u u h {z } and u must satisfy (u u) n ds [[u u]] n ds =0 (u u) n ds =0
39 The interpolant Inspired by DOUGLAS, SANTOS, SHN & Y 1999; CIA, DOUGLAS, Y 1999: Construct interpolant with span 1,x,y, (x)+apple(y) 1,x,y,apple(x)+ (y) (x) =3x 2 10x 4 +7x 6 apple(x) = 3x 2 +5x 4 Define u as orthogonal projection on element
40 rror bound ventually get kũ u h k 2 DG apple C X T h 2 T kuk 2 H 2 (T ) + 2 kuk 2 H 2 (T ) + 2 k div uk 2 H 1 (T ) so not possible to bound λ-dependent term
41 xample: square plate BRNNR SINUM 1993 IIPG SIPG NIPG ALL =(0, 1) 2 µ =1 u 1 (x, y) = sin2 y 1 + cos 2 x u 2 (x, y) = sin2 x 1 cos 2 y sin x sin y sin x sin y
42 A remedy: selective reduced integration (SRI) Recall classical SRI for elasticity problem: to overcome volumetric locking one replaces rs u h : r s v h + (div u)(div v) dv {z } volumetric term by or rs u h : r s v h dv + rs u h : r s v h dv + X w T div u(x T )divv(x T ) T {z } underintegration of volumetric term X T T ( 0 div u)( 0 div v) dv 0 div v = 1 T T div v dv
43 Underintegration applied to (II) gives [[ ]]{div } ds ' 0 [[ ]] 0 {div } ds = 0 {div } 0 [[ ]] ds = 0 {div } [[ ]] ds =0 Replace (IV) with k X 0 [[ ]] 0 [[ ]] = k X = k X [[ 0 ]] [[ ]] 0 [[ ]] [[ ]] = 0 But we now have to check coercivity and consistency of the modified bilinear form!
44 Underintegration applied to (II) gives [[ ]]{div } ds ' = 0 {div } = 0 {div } ventually get 0 [[ ]] 0 {div } ds 0 [[ ]] ds [[ ]] ds =0 a h (, ) apple Ck k DG P T h2 T kuk 2 H kdiv uk 2 H 1 {z } appleckfk 2 1/2 But we now have to check coercivity and consistency of the modified bilinear form!
45 Coercivity NIPG: OK for SRI on terms (II) and (III) SIPG: OK for SRI on terms (II), (III), (IV) IIPG: OK for SRI on terms (III) and (IV) Term (III) when underintegrated leads to a consistency error RI (u, ) = X 0 {div u} tr 0 [[ ]] {div u} tr[[ ]] ds which can be controlled
46 Locking-free (uniformly convergent) behaviour P 1 Q 1 Q 1 +SRI SIP NIP IIP
47 Square plate: DG with triangles quads with and without SRI Standard finite elements / DG quads without SRI DG triangles / quads with SRI
48 NIPG IIPG DG DG DG+SRI DG+SRI SIPG DG DG+SRI
49 Cube with prescribed body force u 1 = (cos 2 x 1)(sin 2 y sin z sin y sin 2 z) u 2 = (cos 2 y 1)(sin 2 z sin x sin z sin 2 x) u 3 = (cos 2 z 1)(sin 2 x sin y sin x sin 2 y) sin x sin y sin z, sin x sin y sin z, sin x sin y sin z.
50 10 2 IIPG 10 1 Q 1 full SG, P 1 full 10 2 NIPG log(h1 error) 10 0 P 1 SRI Q 1 SRI log(h1 error) SG P 1 full Q 1 full Q 1 SRI log(h) Q1 full Q1 UI P1 full SG log(h1 error) log(h) SIPG SG, P 1 full, Q 1 full P 1 SRI Q 1 SRI Q1 full Q1 UI P1 full P1 UI SG Q1 full Q1 UI P1 full P1 UI SG
51 Back to the T-bar example LIU. WHLR, DAWSON 2009 appeared to show good behaviour for IIPG 3D problem but effectively plane stress (Neumann or flux-free boundary condition) in out-of-plane direction, along ABCD Here treated as plane strain (Dirichlet or constrained displacement) along ABCD A C B C Without SRI With SRI
52 Concluding remarks For near-incompressibility DG with quadrilateral elements is not straightforward A remedy, viz. selective under-integration λ-dependent edge terms, has been proposed, analysed, and shown to converge uniformly at the optimal rate Arbitrary quadrilaterals: numerical experiments indicate behaviour similar to that for rectangles. Analysis would be quite complex Current work: nonlinear problems; other parameter-dependent problems References Grieshaber BJ, McBride AT and Reddy BD, Uniformly convergent interior penalty approximations using multilinear approximations for problems in elasticity. SIAM J. Numer. Anal. 53 (2015) Grieshaber BJ, McBride AT and Reddy BD, Computational aspects of Discontinuous Galerkin approximations using quadrilateral and hexahedral elements. In review.
53 Alternative approach: P 1 approximation on quads Much of the manipulations carry over, all the bounds are as before, but this time = u u h 2 V h P 1 u h 2 P 1 Problematic terms: (II) = } ds = [[ ]]{div X ]{div } [[ ]] ds = 0 given properties of interpolant But (IV) = k k X 1 [[ [[ ]] [[w]] [[ ]] ds 6= 6= 0 h so remains a problem for SIPG and IIPG, but is absent for NIPG
54 Square plate, DG with quads and P 1 NIPG IIPG SIPG
PREPRINT 2010:23. A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO
PREPRINT 2010:23 A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO Department of Mathematical Sciences Division of Mathematics CHALMERS UNIVERSITY OF TECHNOLOGY UNIVERSITY OF GOTHENBURG
More 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 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 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 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 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 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 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 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 informationANALYSIS OF THE FEM AND DGM FOR AN ELLIPTIC PROBLEM WITH A NONLINEAR NEWTON BOUNDARY CONDITION
Proceedings of EQUADIFF 2017 pp. 127 136 ANALYSIS OF THE FEM AND DGM FOR AN ELLIPTIC PROBLEM WITH A NONLINEAR NEWTON BOUNDARY CONDITION MILOSLAV FEISTAUER, ONDŘEJ BARTOŠ, FILIP ROSKOVEC, AND ANNA-MARGARETE
More informationLECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES. Sergey Korotov,
LECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES Sergey Korotov, Institute of Mathematics Helsinki University of Technology, Finland Academy of Finland 1 Main Problem in Mathematical
More 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 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 informationNumerical Analysis of Higher Order Discontinuous Galerkin Finite Element Methods
Numerical Analysis of Higher Order Discontinuous Galerkin Finite Element Methods Contents Ralf Hartmann Institute of Aerodynamics and Flow Technology DLR (German Aerospace Center) Lilienthalplatz 7, 3808
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 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 informationGraded Mesh Refinement and Error Estimates of Higher Order for DGFE-solutions of Elliptic Boundary Value Problems in Polygons
Graded Mesh Refinement and Error Estimates of Higher Order for DGFE-solutions of Elliptic Boundary Value Problems in Polygons Anna-Margarete Sändig, Miloslav Feistauer University Stuttgart, IANS Journées
More informationScientific Computing WS 2017/2018. Lecture 18. Jürgen Fuhrmann Lecture 18 Slide 1
Scientific Computing WS 2017/2018 Lecture 18 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 18 Slide 1 Lecture 18 Slide 2 Weak formulation of homogeneous Dirichlet problem Search u H0 1 (Ω) (here,
More 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 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 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 informationIt is known that Morley element is not C 0 element and it is divergent for Poisson equation (see [6]). When Morley element is applied to solve problem
Modied Morley Element Method for a ourth Order Elliptic Singular Perturbation Problem Λ Wang Ming LMAM, School of Mathematical Science, Peking University Jinchao u School of Mathematical Science, Peking
More informationTHE PATCH RECOVERY FOR FINITE ELEMENT APPROXIMATION OF ELASTICITY PROBLEMS UNDER QUADRILATERAL MESHES. Zhong-Ci Shi and Xuejun Xu.
DISCRETE AND CONTINUOUS Website: http://aimsciences.org DYNAMICAL SYSTEMS SERIES B Volume 9, Number, January 2008 pp. 63 82 THE PATCH RECOVERY FOR FINITE ELEMENT APPROXIMATION OF ELASTICITY PROBLEMS UNDER
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 informationHybridized Discontinuous Galerkin Methods
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, 2010 1
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 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 informationA Subspace Correction Method for Discontinuous Galerkin Discretizations of Linear Elasticity Equations
www.oeaw.ac.at A Subspace Correction Method for Discontinuous Galerkin Discretizations of Linear lasticity quations B. Ayuso de Dios, I. Georgiev, J. Kraus, L. Zikatanov RICAM-Report 2011-29 www.ricam.oeaw.ac.at
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 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 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 informationMIXED FINITE ELEMENTS FOR PLATES. Ricardo G. Durán Universidad de Buenos Aires
MIXED FINITE ELEMENTS FOR PLATES Ricardo G. Durán Universidad de Buenos Aires - Necessity of 2D models. - Reissner-Mindlin Equations. - Finite Element Approximations. - Locking. - Mixed interpolation or
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 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 informationPDEs, part 1: Introduction and elliptic PDEs
PDEs, part 1: Introduction and elliptic PDEs Anna-Karin Tornberg Mathematical Models, Analysis and Simulation Fall semester, 2013 Partial di erential equations The solution depends on several variables,
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 informationMultigrid Methods for Maxwell s Equations
Multigrid Methods for Maxwell s Equations Jintao Cui Institute for Mathematics and Its Applications University of Minnesota Outline Nonconforming Finite Element Methods for a Two Dimensional Curl-Curl
More informationA posteriori error estimation for elliptic problems
A posteriori error estimation for elliptic problems Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in
More informationarxiv: v1 [math.na] 27 Jan 2016
Virtual Element Method for fourth order problems: L 2 estimates Claudia Chinosi a, L. Donatella Marini b arxiv:1601.07484v1 [math.na] 27 Jan 2016 a Dipartimento di Scienze e Innovazione Tecnologica, Università
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 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 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 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 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 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 informationDISCRETE MAXIMUM PRINCIPLES in THE FINITE ELEMENT SIMULATIONS
DISCRETE MAXIMUM PRINCIPLES in THE FINITE ELEMENT SIMULATIONS Sergey Korotov BCAM Basque Center for Applied Mathematics http://www.bcamath.org 1 The presentation is based on my collaboration with several
More informationInstitut de Recherche MAthématique de Rennes
LMS Durham Symposium: Computational methods for wave propagation in direct scattering. - July, Durham, UK The hp version of the Weighted Regularization Method for Maxwell Equations Martin COSTABEL & Monique
More informationFinite Element Methods for Fourth Order Variational Inequalities
Louisiana State University LSU Digital Commons LSU Doctoral Dissertations Graduate School 2013 Finite Element Methods for Fourth Order Variational Inequalities Yi Zhang Louisiana State University and Agricultural
More information1. Let a(x) > 0, and assume that u and u h are the solutions of the Dirichlet problem:
Mathematics Chalmers & GU TMA37/MMG800: Partial Differential Equations, 011 08 4; kl 8.30-13.30. Telephone: Ida Säfström: 0703-088304 Calculators, formula notes and other subject related material are not
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 angle conditions in the finite element method. Institute of Mathematics, Academy of Sciences Prague, Czech Republic
On angle conditions in the finite element method Michal Křížek Institute of Mathematics, Academy of Sciences Prague, Czech Republic Joint work with Jan Brandts (University of Amsterdam), Antti Hannukainen
More informationA very short introduction to the Finite Element Method
A very short introduction to the Finite Element Method Till Mathis Wagner Technical University of Munich JASS 2004, St Petersburg May 4, 2004 1 Introduction This is a short introduction to the finite element
More 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 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 informationDifferent Approaches to a Posteriori Error Analysis of the Discontinuous Galerkin Method
WDS'10 Proceedings of Contributed Papers, Part I, 151 156, 2010. ISBN 978-80-7378-139-2 MATFYZPRESS Different Approaces to a Posteriori Error Analysis of te Discontinuous Galerkin Metod I. Šebestová Carles
More informationDiscontinuous Galerkin methods for nonlinear elasticity
Discontinuous Galerkin methods for nonlinear elasticity Preprint submitted to lsevier Science 8 January 2008 The goal of this paper is to introduce Discontinuous Galerkin (DG) methods for nonlinear elasticity
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 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 informationMaximum-norm a posteriori estimates for discontinuous Galerkin methods
Maximum-norm a posteriori estimates for discontinuous Galerkin methods Emmanuil Georgoulis Department of Mathematics, University of Leicester, UK Based on joint work with Alan Demlow (Kentucky, USA) DG
More 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 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 informationChapter 1 Foundations of Elliptic Boundary Value Problems 1.1 Euler equations of variational problems
Chapter 1 Foundations of Elliptic Boundary Value Problems 1.1 Euler equations of variational problems Elliptic boundary value problems often occur as the Euler equations of variational problems the latter
More informationMULTIGRID METHODS FOR MAXWELL S EQUATIONS
MULTIGRID METHODS FOR MAXWELL S EQUATIONS A Dissertation Submitted to the Graduate Faculty of the Louisiana State University and Agricultural and Mechanical College in partial fulfillment of the requirements
More informationVariational Formulations
Chapter 2 Variational Formulations In this chapter we will derive a variational (or weak) formulation of the elliptic boundary value problem (1.4). We will discuss all fundamental theoretical results that
More informationDISCONTINUOUS GALERKIN FINITE ELEMENT DISCRETIZATION FOR STEADY STOKES FLOWS WITH THRESHOLD SLIP BOUNDARY CONDITION
DISCONTINUOUS GALRKIN FINIT LMNT DISCRTIZATION FOR STADY STOKS FLOWS WITH THRSHOLD SLIP BOUNDARY CONDITION J.K. Djoko Department of Mathematics and Applied Mathematics, University of Pretoria, Private
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 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 informationA Least-Squares Finite Element Approximation for the Compressible Stokes Equations
A Least-Squares Finite Element Approximation for the Compressible Stokes Equations Zhiqiang Cai, 1 Xiu Ye 1 Department of Mathematics, Purdue University, 1395 Mathematical Science Building, West Lafayette,
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 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 informationA BIVARIATE SPLINE METHOD FOR SECOND ORDER ELLIPTIC EQUATIONS IN NON-DIVERGENCE FORM
A BIVARIAE SPLINE MEHOD FOR SECOND ORDER ELLIPIC EQUAIONS IN NON-DIVERGENCE FORM MING-JUN LAI AND CHUNMEI WANG Abstract. A bivariate spline method is developed to numerically solve second order elliptic
More informationPDE & FEM TERMINOLOGY. BASIC PRINCIPLES OF FEM.
PDE & FEM TERMINOLOGY. BASIC PRINCIPLES OF FEM. Sergey Korotov Basque Center for Applied Mathematics / IKERBASQUE http://www.bcamath.org & http://www.ikerbasque.net 1 Introduction The analytical solution
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 informationChapter 5 A priori error estimates for nonconforming finite element approximations 5.1 Strang s first lemma
Chapter 5 A priori error estimates for nonconforming finite element approximations 51 Strang s first lemma We consider the variational equation (51 a(u, v = l(v, v V H 1 (Ω, and assume that the conditions
More informationChapter 3 Conforming Finite Element Methods 3.1 Foundations Ritz-Galerkin Method
Chapter 3 Conforming Finite Element Methods 3.1 Foundations 3.1.1 Ritz-Galerkin Method Let V be a Hilbert space, a(, ) : V V lr a bounded, V-elliptic bilinear form and l : V lr a bounded linear functional.
More informationChapter 1: The Finite Element Method
Chapter 1: The Finite Element Method Michael Hanke Read: Strang, p 428 436 A Model Problem Mathematical Models, Analysis and Simulation, Part Applications: u = fx), < x < 1 u) = u1) = D) axial deformation
More informationA Finite Element Method Using Singular Functions for Poisson Equations: Mixed Boundary Conditions
A Finite Element Method Using Singular Functions for Poisson Equations: Mixed Boundary Conditions Zhiqiang Cai Seokchan Kim Sangdong Kim Sooryun Kong Abstract In [7], we proposed a new finite element method
More informationXIAO-CHUAN CAI AND MAKSYMILIAN DRYJA. strongly elliptic equations discretized by the nite element methods.
Contemporary Mathematics Volume 00, 0000 Domain Decomposition Methods for Monotone Nonlinear Elliptic Problems XIAO-CHUAN CAI AND MAKSYMILIAN DRYJA Abstract. In this paper, we study several overlapping
More informationNotes for Math 663 Spring Alan Demlow
Notes for Math 663 Spring 2016 Alan Demlow 2 Chapter 1 Introduction and preliminaries 1.1 Introduction In this chapter we will introduce the main topics of the course. We first briefly define the finite
More informationFinite Elements. Colin Cotter. February 22, Colin Cotter FEM
Finite Elements February 22, 2019 In the previous sections, we introduced the concept of finite element spaces, which contain certain functions defined on a domain. Finite element spaces are examples of
More 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 informationHigh order, finite volume method, flux conservation, finite element method
FLUX-CONSERVING FINITE ELEMENT METHODS SHANGYOU ZHANG, ZHIMIN ZHANG, AND QINGSONG ZOU Abstract. We analyze the flux conservation property of the finite element method. It is shown that the finite element
More informationChapter 6 A posteriori error estimates for finite element approximations 6.1 Introduction
Chapter 6 A posteriori error estimates for finite element approximations 6.1 Introduction The a posteriori error estimation of finite element approximations of elliptic boundary value problems has reached
More informationSome applications of discontinuous Galerkin methods in solid mechanics
Some applications of discontinuous Galerkin methods in solid mechanics Adrian Lew, Alex Ten Eyck, and Ramsharan Rangarajan Mechanical Engineering, Stanford University, USA. Email: {lewa,alabute,rram}@stanford.edu
More informationA posteriori error estimation of approximate boundary fluxes
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2008; 24:421 434 Published online 24 May 2007 in Wiley InterScience (www.interscience.wiley.com)..1014 A posteriori error estimation
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 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 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 informationWELL POSEDNESS OF PROBLEMS I
Finite Element Method 85 WELL POSEDNESS OF PROBLEMS I Consider the following generic problem Lu = f, where L : X Y, u X, f Y and X, Y are two Banach spaces We say that the above problem is well-posed (according
More informationReview of A Priori Error Estimation for Discontinuous Galerkin Methods
TICAM REPORT 00-27 October 7, 2000 Review of A Priori Error Estimation for Discontinuous Galerkin Methods S. Prudhomme Λ, F. Pascal y, J.T. Oden Λ, and A. Romkes Λ Λ Texas Institute for Computational and
More informationSolutions of Selected Problems
1 Solutions of Selected Problems October 16, 2015 Chapter I 1.9 Consider the potential equation in the disk := {(x, y) R 2 ; x 2 +y 2 < 1}, with the boundary condition u(x) = g(x) r for x on the derivative
More informationNumerical Methods for the Navier-Stokes equations
Arnold Reusken Numerical Methods for the Navier-Stokes equations January 6, 212 Chair for Numerical Mathematics RWTH Aachen Contents 1 The Navier-Stokes equations.............................................
More informationAdaptive Finite Element Methods Lecture Notes Winter Term 2017/18. R. Verfürth. Fakultät für Mathematik, Ruhr-Universität Bochum
Adaptive Finite Element Methods Lecture Notes Winter Term 2017/18 R. Verfürth Fakultät für Mathematik, Ruhr-Universität Bochum Contents Chapter I. Introduction 7 I.1. Motivation 7 I.2. Sobolev and finite
More 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 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 informationA Finite Element Method for an Ill-Posed Problem. Martin-Luther-Universitat, Fachbereich Mathematik/Informatik,Postfach 8, D Halle, Abstract
A Finite Element Method for an Ill-Posed Problem W. Lucht Martin-Luther-Universitat, Fachbereich Mathematik/Informatik,Postfach 8, D-699 Halle, Germany Abstract For an ill-posed problem which has its origin
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 informationMixed Hybrid Finite Element Method: an introduction
Mixed Hybrid Finite Element Method: an introduction First lecture Annamaria Mazzia Dipartimento di Metodi e Modelli Matematici per le Scienze Applicate Università di Padova mazzia@dmsa.unipd.it Scuola
More informationThe Plane Stress Problem
The Plane Stress Problem Martin Kronbichler Applied Scientific Computing (Tillämpad beräkningsvetenskap) February 2, 2010 Martin Kronbichler (TDB) The Plane Stress Problem February 2, 2010 1 / 24 Outline
More information