Unified Hybridization Of Discontinuous Galerkin, Mixed, And Continuous Galerkin Methods For Second Order Elliptic Problems.
|
|
- Melvin Cain
- 5 years ago
- Views:
Transcription
1 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 of January, 2011 Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
2 Inhalt 1 introduction (an example) Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
3 Inhalt 1 introduction (an example) 2 the framework Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
4 Inhalt 1 introduction (an example) 2 the framework 3 examples for hybridizable methods Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
5 Inhalt 1 introduction (an example) 2 the framework 3 examples for hybridizable methods 4 Other novel methods hybridizable methods well suited for adaptivity The RT-method on meshes with hanging nodes Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
6 Inhalt literature D.N. Arnold, F. Brezzi, B. Cockburn, and L.D. Marini. Unified analysis of discontinuous Galerkin methods for elliptic problems. SIAM journal on numerical analysis, 39(5): , B. Cockburn and J. Gopalakrishnan. A characterization of hybridized mixed methods for second order elliptic problems. SIAM Journal on Numerical Analysis, 42(1): , B. Cockburn, J. Gopalakrishnan, and R. Lazarov. Unified hybridization of discontinuous Galerkin, mixed and continuous Galerkin methods for second order elliptic problems. SIAM J. Numer. Anal, 47(2): , R.M. Kirby, S.J. Sherwin, and B. Cockburn. To CG or to HDG: A Comparative Study. M. Vohraĺık, J. Maryska, and O. Severýn. Mixed and nonconforming finite element methods on a system of polygons. Applied Numerical Mathematics, 57(2): , OC Zienkiewicz. Displacement and equilibrium models in the finite element method by B. Fraeijs de Veubeke, Chapter 9, Pages of Stress Analysis, Edited by OC Zienkiewicz and GS Holister, Published by John Wiley & Sons, International Journal for Numerical Methods in Engineering, 52(3): , Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
7 introduction (an example) 1 introduction (an example) 2 the framework 3 examples for hybridizable methods 4 Other novel methods hybridizable methods well suited for adaptivity The RT-method on meshes with hanging nodes Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
8 introduction (an example) Second order elliptic boundary value problem Search solution u for (a u)+du = f in Ω Ê n, u = g on Ω, where a(x) is a bounded symmetric positive definite matrix-valued function, Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
9 introduction (an example) Second order elliptic boundary value problem Search solution u for (a u)+du = f in Ω Ê n, u = g on Ω, where a(x) is a bounded symmetric positive definite matrix-valued function, d(x) is a bounded nonnegative function. Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
10 introduction (an example) Second order elliptic boundary value problem Search solution u for (a u)+du = f in Ω Ê n, u = g on Ω, where a(x) is a bounded symmetric positive definite matrix-valued function, d(x) is a bounded nonnegative function. Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
11 introduction (an example) Second order elliptic boundary value problem Search solution u for where (a u)+du = f in Ω Ê n, u = g on Ω, a(x) is a bounded symmetric positive definite matrix-valued function, d(x) is a bounded nonnegative function. This can be rewritten mixed form Search solution (q, u) for q+a u = 0 in Ω, q+du = f in Ω, u = g on Ω. Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
12 introduction (an example) We re choosing a triangulation of Ω, T h T h Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
13 introduction (an example) We re choosing a triangulation of Ω, T h and using Raviart-Thomas elements P k (K) n +xp k (K) T h RT 0 RT 1 RT 2 Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
14 introduction (an example) Notation E h interior edges, E h boundary faces of Ω, Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
15 introduction (an example) Notation Eh interior edges, E h boundary faces of Ω, E h = Eh E h, Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
16 introduction (an example) Notation Eh interior edges, E h boundary faces of Ω, E h = Eh E h, (u,v) D := D uvdx for u,v L2 (D), D Ê n, Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
17 introduction (an example) Notation Eh interior edges, E h boundary faces of Ω, E h = Eh E h, (u,v) D := D uvdx for u,v L2 (D), D Ê n, u,v D := D uvdx for u,v L2 (D), D Ê n 1, Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
18 introduction (an example) Notation Eh interior edges, E h boundary faces of Ω, E h = Eh E h, (u,v) D := D uvdx for u,v L2 (D), D Ê n, u,v D := D uvdx for u,v L2 (D), D Ê n 1, (u,v) Th = K T h (v,w) K for u,v on Ω, Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
19 introduction (an example) Notation Eh interior edges, E h boundary faces of Ω, E h = Eh E h, (u,v) D := D uvdx for u,v L2 (D), D Ê n, u,v D := D uvdx for u,v L2 (D), D Ê n 1, (u,v) Th = K T h (v,w) K for u,v on Ω, µ,λ E = e E µ,λ e for µ,λ on E E h. Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
20 introduction (an example) Notation Eh interior edges, E h boundary faces of Ω, E h = Eh E h, (u,v) D := D uvdx for u,v L2 (D), D Ê n, u,v D := D uvdx for u,v L2 (D), D Ê n 1, (u,v) Th = K T h (v,w) K for u,v on Ω, µ,λ E = e E µ,λ e for µ,λ on E E h. H(div,Ω) := {v L 2 (Ω) n v L 2 (Ω)} Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
21 introduction (an example) RT finite element approximation The approximation (q h,u h ) is sought in the finite element space V RT h Wh RT given by V RT h = {v H(div,Ω) v K P k (K) n +xp k (K) K T h }, W RT h = {w L 2 (Ω) w K P k (K) K T h }. and is defined by requiring that, Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
22 introduction (an example) RT finite element approximation The approximation (q h,u h ) is sought in the finite element space V RT h Wh RT given by V RT h = {v H(div,Ω) v K P k (K) n +xp k (K) K T h }, W RT h = {w L 2 (Ω) w K P k (K) K T h }. and is defined by requiring that, (cq h,v) Ω (u h, v) Ω = g,v n Ω (w, q h ) Ω (du h,w) Ω = (f,w) Ω v V RT h, w W RT h, with c = a 1 for each element K T h Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
23 introduction (an example) This yields to a matrix equation of the form: ( t )( ) ( É = Í À ), where É and Í are vectors of coefficients of q h and u h with respect to their corresponding finite element basis. Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
24 introduction (an example) problem System is not positive definite. Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
25 introduction (an example) problem System is not positive definite. solutions using expensive solver, Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
26 introduction (an example) problem System is not positive definite. solutions using expensive solver, using a positive definite system by elimination of É: solve equation ( 1 t + )Í = + 1 for Í. Required inversion 1 is difficult to compute and a full matrix. Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
27 introduction (an example) problem System is not positive definite. solutions using expensive solver, using a positive definite system by elimination of É: solve equation ( 1 t + )Í = + 1 for Í. Required inversion 1 is difficult to compute and a full matrix. hybridization: Introduce new unkowns λ h (Lagrangian multipliers) and relax the continuity constraints between element interfaces for q h. Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
28 introduction (an example) Definition (hybridized RT finite element approximation) The approximation (q h,u h,λ h ) is sought in the finite element space V h RT Ŵ h RT M h RT given by V h RT = {v L 2 (Ω) n v K P k (K) n +xp k (K) K T h }, Ŵ h RT = {w L 2 (Ω) w K P k (K) K T h }, M hrt = {µ L 2 (E h ) µ e P k (e) e E h } Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
29 introduction (an example) Definition (hybridized RT finite element approximation) The approximation (q h,u h,λ h ) is sought in the finite element space V h RT Ŵ h RT M h RT given by RT V h = {v L 2 (Ω) n v K P k (K) n +xp k (K) K T h }, RT Ŵ h = {w L 2 (Ω) w K P k (K) K T h }, M hrt = {µ L 2 (Eh ) µ e P k (e) e Eh } and is defined by requiring that, (cq h,v) Ω (u h, v) Ω + λ h, v E h = g,v n Ω v V RT h, (w, q h ) Ω (du h,w) Ω = (f,w) Ω w ŴhRT, µ, q h E h = 0 µ M hrt. Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
30 introduction (an example) This yields to a matrix of the form A B t C t Q G B D 0 U = F, C 0 0 Λ 0 where Λ is the vector of dofs associated to the multiplier λ h. Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
31 introduction (an example) This yields to a matrix of the form A B t C t Q G B D 0 U = F, C 0 0 Λ 0 where Λ is the vector of dofs associated to the multiplier λ h. New vectors of dofs Q and U define the same approximation (q h,u h ) as the original method. Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
32 introduction (an example) This yields to a matrix of the form A B t C t Q G B D 0 U = F, C 0 0 Λ 0 where Λ is the vector of dofs associated to the multiplier λ h. New vectors of dofs Q and U define the same approximation (q h,u h ) as the original method. Both Q and U can be easily eliminated Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
33 introduction (an example) This yields to a matrix of the form A B t C t Q G B D 0 U = F, C 0 0 Λ 0 where Λ is the vector of dofs associated to the multiplier λ h. New vectors of dofs Q and U define the same approximation (q h,u h ) as the original method. Both Q and U can be easily eliminated Λ = À, =CA 1 (A B t (BA 1 B t +D) 1 B)A 1 C t, À = CA 1 (A B t (BA 1 B t +D) 1 B)A 1 G CA 1 B t (BA 1 B t +D) 1 F. Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
34 introduction (an example) advantages of hybridization A is now a block diagonal matrix, Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
35 introduction (an example) advantages of hybridization A is now a block diagonal matrix, is symmetric positive definite, Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
36 introduction (an example) advantages of hybridization A is now a block diagonal matrix, is symmetric positive definite, number of dofs is remarkably smaller, Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
37 introduction (an example) advantages of hybridization A is now a block diagonal matrix, is symmetric positive definite, number of dofs is remarkably smaller, once Λ has been obtained, both Q and U can be computed efficiently element by element, Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
38 introduction (an example) advantages of hybridization A is now a block diagonal matrix, is symmetric positive definite, number of dofs is remarkably smaller, once Λ has been obtained, both Q and U can be computed efficiently element by element, the multiplier λ h can be used to improve the approximation to u by means of local post-processing. Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
39 introduction (an example) history of hybridization year event 1965 first hybridization of finite elements for solving equations of linear elasticity (often: static condensation)/ implementation trick [6] 1985 proof: hybrid variable contains extra information about the exact solution ( local post-processing)[5] 2004 hybridized RT and BDM methods of arbitrary order[2] 2005 extended to finite element methods for stationary Stokes equation (DG, mixed) 2009 unifying framework[3] 2010 comparison CG-H vs. LDG-H[4] Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
40 the framework 1 introduction (an example) 2 the framework 3 examples for hybridizable methods 4 Other novel methods hybridizable methods well suited for adaptivity The RT-method on meshes with hanging nodes Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
41 the framework Unified framework provides approximations for 1 (q,u) in the interior of the elements K T h, (q h,u h ), 2 u on the interior border of the elements, λ h. Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
42 the framework Definition (local solvers) For any single valued function m L 2 ( K), the functions (Qm,Um) are the solutions of cqm+ Um = 0 on K, Qm+dUm = 0 on K, Um = m on K Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
43 the framework Definition (local solvers) For any single valued function m L 2 ( K), the functions (Qm,Um) are the solutions of cqm+ Um = 0 on K, Qm+dUm = 0 on K, Um = m on K and for any single valued function f L 2 (K), the functions (Qf,Uf) are the solutions of cqf + Uf = 0 on K, Qf +duf = f on K, Uf = 0 on K. Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
44 the framework With u = λ + g and the linearity of the problem we have that (q,u) = (Qλ+Qg +Qf,Uλ+Ug +Uf). Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
45 the framework With u = λ + g and the linearity of the problem we have that (q,u) = (Qλ+Qg +Qf,Uλ+Ug +Uf). The above property only holds iff Definition (transmission condition) This completely characterizes λ. Qλ+Qg +Qf = 0 Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
46 the framework Definition (discrete local solvers) (Qm,Um) and (Qf,Uf) are the discrete versions of the local solvers. Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
47 the framework Definition (discrete local solvers) (Qm,Um) and (Qf,Uf) are the discrete versions of the local solvers. Then the discrete solution can be written as (q h,u h ) = (Qλ h +Qg h +Qf,Uλ h +Ug h +Uf), where λ h M h and g h M h are approximations to the values of u. u h λ h M h g h M h Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
48 the framework Definition (discrete version of transmission condition (weak)) a h (λ h,µ) = b h (µ) µ M h is the discrete version of the transmission condition. Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
49 the framework Notation V(K) polynomial space in which q is approximated, Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
50 the framework Notation V(K) polynomial space in which q is approximated, W(K) polynomial space in which u is approximated, Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
51 the framework Notation V(K) polynomial space in which q is approximated, W(K) polynomial space in which u is approximated, V h = {v v K V(K)}, Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
52 the framework Notation V(K) polynomial space in which q is approximated, W(K) polynomial space in which u is approximated, V h = {v v K V(K)}, W h = {w w K W(K)}. Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
53 the framework Notation V(K) polynomial space in which q is approximated, W(K) polynomial space in which u is approximated, V h = {v v K V(K)}, W h = {w w K W(K)}. In general: A function h on K is double-valued, each branch is denoted by h K + or h K. Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
54 the framework To define a hybridizable method, we had to define the discrete local solvers (Qm,Um) and (Qf,Uf), Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
55 the framework To define a hybridizable method, we had to define the discrete local solvers (Qm,Um) and (Qf,Uf), the transmission condition determining λ h, a h and b h, Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
56 the framework To define a hybridizable method, we had to define the discrete local solvers (Qm,Um) and (Qf,Uf), the transmission condition determining λ h, a h and b h, the finite element spaces Mh and M h, Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
57 the framework To define a hybridizable method, we had to define the discrete local solvers (Qm,Um) and (Qf,Uf), the transmission condition determining λ h, a h and b h, the finite element spaces Mh and M h, the finite element spaces V(K) and W(K), Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
58 the framework To define a hybridizable method, we had to define the discrete local solvers (Qm,Um) and (Qf,Uf), the transmission condition determining λ h, a h and b h, the finite element spaces Mh and M h, the finite element spaces V(K) and W(K), the trace spaces ˆQm and ˆQf (used in discrete local solvers). Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
59 the framework To define a hybridizable method, we had to define the discrete local solvers (Qm,Um) and (Qf,Uf), the transmission condition determining λ h, a h and b h, the finite element spaces Mh and M h, the finite element spaces V(K) and W(K), the trace spaces ˆQm and ˆQf (used in discrete local solvers). Finally, think about existence and uniqueness of λ h. Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
60 the framework Definition (trace space M h ) M h := {µ M h µ = 0 on Ω}. Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
61 the framework Definition (discrete local solver) 1 maps m M h to the function (Qm,Um): (cqm,v) K (Um, v) K = m,v n K ( w,qm) K + w, ˆQm n K +(dum,w) K = 0 v V(K), w W(K). Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
62 the framework Definition (discrete local solver) 1 maps m M h to the function (Qm,Um): (cqm,v) K (Um, v) K = m,v n K ( w,qm) K + w, ˆQm n K +(dum,w) K = 0 v V(K), w W(K). 2 maps f L 2 (Ω) to the pair (Qf,Uf): (cqf,v) K (Uf, v) K = 0 ( w,qf) K + w, ˆQf n K +(duf,w) K = (f,w) K v V(K), w W(K). Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
63 the framework Definition (characterization of λ h ) a h (η,µ) = (cqη,qµ) Th +(duη,uµ) Th + 1, (Uµ µ)(ˆqη Qη) Eh, b h (µ) = g h, ˆQµ Eh +(f,uµ) Th 1, (Uµ µ)(ˆqf Qf) Eh + 1, Uf(ˆQµ Qµ) Eh 1, (Uµ µ)(ˆqg h Qg h ) Eh + 1, (Ug h g)(ˆqµ Qµ) Eh. Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
64 the framework Theorem (existence and uniqueness of λ h ) If the three assumptions are fulfilled, then there is a unique solution λ h of the weak formulation. Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
65 the framework Theorem (existence and uniqueness of λ h ) If the three assumptions 1 existence and uniqueness of the local solvers, 2 3 are fulfilled, then there is a unique solution λ h of the weak formulation. Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
66 the framework Theorem (existence and uniqueness of λ h ) If the three assumptions 1 existence and uniqueness of the local solvers, 2 positive semidefiniteness of the local solvers and 3 are fulfilled, then there is a unique solution λ h of the weak formulation. Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
67 the framework Theorem (existence and uniqueness of λ h ) If the three assumptions 1 existence and uniqueness of the local solvers, 2 positive semidefiniteness of the local solvers and 3 the gluing condition (If µ Mh, then for every interior face there exists a branch K with P K µ = µ.) are fulfilled, then there is a unique solution λ h of the weak formulation. Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
68 examples for hybridizable methods 1 introduction (an example) 2 the framework 3 examples for hybridizable methods 4 Other novel methods hybridizable methods well suited for adaptivity The RT-method on meshes with hanging nodes Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
69 examples for hybridizable methods general assumptions same local solvers in every element K, Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
70 examples for hybridizable methods general assumptions same local solvers in every element K, conforming simplicial triangulation. Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
71 examples for hybridizable methods general assumptions same local solvers in every element K, conforming simplicial triangulation. M c h,k := {µ C(E h) µ e P k (e) for all faces e E h } M h,k := {µ L 2 (E h ) µ e P k (e) for all faces e E h } Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
72 examples for hybridizable methods RT-H (Raviart Thomas): V(K) W(K) P k (K) n +xp k (K) P k (K) M h M h,k ˆQm Qm ˆQf Qf a h (η,µ) (cqη,qµ) Th +(duη,uµ) Th b h (µ) g h,qµ n Ω +(f,uµ) Th Conservativity strong Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
73 examples for hybridizable methods BDM-H (Brezzi Douglas Marini): V(K) W(K) P k (K) n P k 1 (K) M h M h,k ˆQm Qm ˆQf Qf a h (η,µ) (cqη,qµ) Th +(duη,uµ) Th b h (µ) g h,qµ n Ω +(f,uµ) Th Conservativity strong Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
74 examples for hybridizable methods LDG-H (Local Discontinuous Galerkin): V(K) W(K) 1.) P k (K) n P k 1 (K) 2.) P k (K) n P k (K) 3.) P k 1 (K) n P k (K) M h M h,k ˆQm Qm+τ K (Um m)n ˆQf Qf +τ K (Uf)n a h (η,µ) (cqη,qµ) Th +(duη,uµ) Th + 1, (Uµ µ)(τ K (Uη η)n) Eh b h (µ) g h,qµ n+τ K Uµ Ω +(f,uµ) Th Conservativity strong The stabilization parameter τ K is nonnegative constant on each face in E h, double-valued on Eh. Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
75 examples for hybridizable methods CG-H (Continuous Galerkin): V(K) W(K) P k 1 (K) n P k (K) M h M c h,k ˆQm a new unkown variable ˆQf a new unkown variable a h (η,µ) (a Uη, Uµ) Th +(duη,uµ) Th b h (µ) g h, ˆQµ Eh +(f,uµ) Th Conservativity weak Assume that a(x) is a constant on each element. CG-H is an LDG-H method with τ. Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
76 examples for hybridizable methods Are all methods hybridizable? Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
77 examples for hybridizable methods Are all methods hybridizable? No. Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
78 examples for hybridizable methods Are all methods hybridizable? No. Example (DG-methods) Remember the lecture Scientific Computing or see [1]. A hybrizable method has to be single valued for û h (= λ h ). Method û h Bassi-Rebay {u h } Brezzi et al. {u h } LDG {u h } β u h IP {u h } Bassi et al. {u h } Baumann-Oden {u h }+n K u h NIPG {u h }+n K u h Babuska-Zlamal (u h K ) K Brezzi et al. (u h K ) K Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
79 Other novel methods 1 introduction (an example) 2 the framework 3 examples for hybridizable methods 4 Other novel methods hybridizable methods well suited for adaptivity The RT-method on meshes with hanging nodes Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
80 Other novel methods hybridizable methods well suited for adaptivity 1 introduction (an example) 2 the framework 3 examples for hybridizable methods 4 Other novel methods hybridizable methods well suited for adaptivity The RT-method on meshes with hanging nodes Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
81 Other novel methods hybridizable methods well suited for adaptivity hybridizable methods well suited for adaptivity: V(K) W(K) 1.) P k(k) (K) n +xp k(k) P k(k) (K) 2.) P k(k) (K) n P k(k) 1 (K) 3.) P k(k) (K) n P k(k) (K) 4.) P k(k) 1 (K) n P k(k) (K) Mh {µ µ e M h,k(e) e Eh } {µ µ K C({x K τ K (x) = })} ˆQm Qm+τ K (Um m)n ˆQf Qf +τ K (Uf)n a h (η,µ) depending on V(K),W(K) b h (µ) depending on V(K),W(K) with τ K 0 for 3.) (on at least one face) and 4.) (for alle faces). Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
82 Other novel methods hybridizable methods well suited for adaptivity For e = K + K, we set max{k(k + ),k(k )}, if τ + < and τ < k(e) := k(k ± ), if τ ± = and τ < min{k(k + ),k(k )}, if τ + = and τ =. Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
83 Other novel methods hybridizable methods well suited for adaptivity assumption K τ K (0, ) interior face E h [0, ] τ K face of T h corresponding to LDG-H depending on the selected method Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
84 Other novel methods hybridizable methods well suited for adaptivity Main features of this class of methods: 1 variable degree approximation spaces on conforming meshes (k(k)), Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
85 Other novel methods hybridizable methods well suited for adaptivity Main features of this class of methods: 1 variable degree approximation spaces on conforming meshes (k(k)), 2 automatical coupling of different methods on conforming meshes (τ K determines method), Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
86 Other novel methods hybridizable methods well suited for adaptivity Main features of this class of methods: 1 variable degree approximation spaces on conforming meshes (k(k)), 2 automatical coupling of different methods on conforming meshes (τ K determines method), 3 mortaring capabilities for nonconforming meshes (choice of τ for hanging nodes). Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
87 Other novel methods The RT-method on meshes with hanging nodes 1 introduction (an example) 2 the framework 3 examples for hybridizable methods 4 Other novel methods hybridizable methods well suited for adaptivity The RT-method on meshes with hanging nodes Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
88 Other novel methods The RT-method on meshes with hanging nodes Considere the case of variable degree RT-H method with τ 0 everywhere (assumption in last method is violated). Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
89 Other novel methods The RT-method on meshes with hanging nodes Considere the case of variable degree RT-H method with τ 0 everywhere (assumption in last method is violated). Mesh is locally refined into four congruent triangles. Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
90 Other novel methods The RT-method on meshes with hanging nodes Considere the case of variable degree RT-H method with τ 0 everywhere (assumption in last method is violated). Mesh is locally refined into four congruent triangles. Idea: Impose special conditions on the meshes and link the definition of k(k) to the structure of the mesh. Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
91 Other novel methods The RT-method on meshes with hanging nodes k(k 2 ) " " max{k(k 2 ),k(k 4 )} k(k 4 ) arbitrarily max{k(k 3 ),k(k 4 )} k(k 3 ) " " k(k 1 ) Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
92 Other novel methods The RT-method on meshes with hanging nodes Thanks for your attention! Stefan Girke (WWU Münster Institut für Numerische und Angewandte Hybridization Mathematik) 10th of January, / 43
Hybridized 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 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 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 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 informationAnalysis of Hybrid Discontinuous Galerkin Methods for Incompressible Flow Problems
Analysis of Hybrid Discontinuous Galerkin Methods for Incompressible Flow Problems Christian Waluga 1 advised by Prof. Herbert Egger 2 Prof. Wolfgang Dahmen 3 1 Aachen Institute for Advanced Study in Computational
More 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 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 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 informationStatic condensation, hybridization, and the devising of the HDG methods
Static condensation, hybridization, and the devising of the HDG methods Bernardo Cockburn Abstract In this paper, we review and refine the main ideas for devising the so-called hybridizable discontinuous
More informationarxiv: v1 [math.na] 19 Dec 2018
THE WEAK GALERKIN METHODS ARE REWRITINGS OF THE HYBRIDIZABLE DISCONTINUOUS GALERKIN METHODS arxiv:1812.08146v1 [math.na] 19 Dec 2018 BERNARDO COCKBURN Abstract. We establish that the Weak Galerkin methods
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 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 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 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 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 informationTo CG or to HDG: A Comparative Study
To CG or to HDG: A Comparative Study Robert M. Kirby a,1, Spencer J. Sherwin b,3 Bernardo Cockburn c,2, a School of Computing, Univ. of Utah, Salt Lake City, UT, USA b Department of Aeronautics, Imperial
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 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 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 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 informationMultilevel Preconditioning of Graph-Laplacians: Polynomial Approximation of the Pivot Blocks Inverses
Multilevel Preconditioning of Graph-Laplacians: Polynomial Approximation of the Pivot Blocks Inverses P. Boyanova 1, I. Georgiev 34, S. Margenov, L. Zikatanov 5 1 Uppsala University, Box 337, 751 05 Uppsala,
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 informationPerformance comparison between hybridizable DG and classical DG methods for elastic waves simulation in harmonic domain
March 4-5, 2015 Performance comparison between hybridizable DG and classical DG methods for elastic waves simulation in harmonic domain M. Bonnasse-Gahot 1,2, H. Calandra 3, J. Diaz 1 and S. Lanteri 2
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 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 informationarxiv: v1 [math.na] 14 Apr 2011
A WEAK GALERKIN FINIE ELEMEN MEHOD FOR SECOND-ORDER ELLIPIC PROBLEMS JUNPING WANG AND XIU YE arxiv:1104.2897v1 [math.na] 14 Apr 2011 Abstract. In this paper, authors shall introduce a finite element method
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 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 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 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 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 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 informationMIXED FINITE ELEMENT METHODS FOR PROBLEMS WITH ROBIN BOUNDARY CONDITIONS
MIXED FINITE ELEMENT METHODS FOR PROBLEMS WITH ROBIN BOUNDARY CONDITIONS JUHO KÖNNÖ, DOMINIK SCHÖTZAU, AND ROLF STENBERG Abstract. We derive new a-priori and a-posteriori error estimates for mixed nite
More 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 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 informationDiscontinuous Galerkin method for a class of elliptic multi-scale problems
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 000; 00: 6 [Version: 00/09/8 v.0] Discontinuous Galerkin method for a class of elliptic multi-scale problems Ling Yuan
More informationSuperconvergence and H(div) Projection for Discontinuous Galerkin Methods
INTRNATIONAL JOURNAL FOR NUMRICAL MTHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2000; 00:1 6 [Version: 2000/07/27 v1.0] Superconvergence and H(div) Projection for Discontinuous Galerkin Methods Peter Bastian
More informationMATH 676. Finite element methods in scientific computing
MATH 676 Finite element methods in scientific computing Wolfgang Bangerth, Texas A&M University Lecture 33.25: Which element to use Part 2: Saddle point problems Consider the stationary Stokes equations:
More 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 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 informationDiscontinuous Galerkin and Finite Difference Methods for the Acoustic Equations with Smooth Coefficients. Mario Bencomo TRIP Review Meeting 2013
About Me Mario Bencomo Currently 2 nd year graduate student in CAAM department at Rice University. B.S. in Physics and Applied Mathematics (Dec. 2010). Undergraduate University: University of Texas at
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 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 informationarxiv: v2 [math.na] 23 Apr 2016
Improved ZZ A Posteriori Error Estimators for Diffusion Problems: Conforming Linear Elements arxiv:508.009v2 [math.na] 23 Apr 206 Zhiqiang Cai Cuiyu He Shun Zhang May 2, 208 Abstract. In [8], we introduced
More informationTHE MULTISCALE DISCONTINUOUS GALERKIN METHOD FOR SOLVING A CLASS OF SECOND ORDER ELLIPTIC PROBLEMS WITH ROUGH COEFFICIENTS
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 8, Number, Pages 28 47 c 20 Institute for Scientific Computing and Information THE MULTISCALE DISCONTINUOUS GALERKIN METHOD FOR SOLVING A
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 informationA SIMPLE TUTORIAL ON DISCONTINUOUS GALERKIN METHODS. Jennifer Proft CERMICS, ENPC. J. Proft CERMICS, ENPC
A SIMPLE TUTORIAL ON DISCONTINUOUS GALERKIN METHODS Jennifer Proft Outline Definitions A simple example Issues Historical development elliptic equations hyperbolic equations Discontinuous vs. continuous
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 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 informationBubble stabilized discontinuous Galerkin methods on conforming and non-conforming meshes
Calcolo 2011) 48: 189 209 DOI 10.1007/s10092-010-0036-3 Bubble stabilized discontinuous Galerkin methods on conforming and non-conforming meshes Erik Burman Benjamin Stamm Received: 4 January 2010 / Accepted:
More informationUNIFIED A POSTERIORI ERROR ESTIMATOR FOR FINITE ELEMENT METHODS FOR THE STOKES EQUATIONS
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 10, Number 3, Pages 551 570 c 013 Institute for Scientific Computing and Information UNIFIED A POSTERIORI ERROR ESTIMATOR FOR FINITE ELEMENT
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 informationAN ANALYSIS OF THE EMBEDDED DISCONTINUOUS GALERKIN METHOD FOR SECOND ORDER ELLIPTIC PROBLEMS
AN ANALYSIS OF THE EMBEDDED DISCONTINUOUS GALERKIN METHOD FOR SECOND ORDER ELLIPTIC PROBLEMS BERNARDO COCKBURN, JOHNNY GUZMÁN, SEE-CHEW SOON, AND HENRYK K. STOLARSKI Abstract. Te embedded discontinuous
More informationError analysis of piecewise constant approximations of Darcy s law
Error analysis of piecewise constant approximations of Darcy s law F. Brezzi 1,2, M. Fortin 3, L.D. Marini 1,2 January 17, 2005 Keywords Darcy, Finite Elements, Mixed formulations, Constant pressures,
More informationSTABILIZATION IN RELATION TO WAVENUMBER IN HDG METHODS
STABILIZATION IN RELATION TO WAVENUMBER IN HDG METHODS J. GOPALAKRISHNAN, S. LANTERI, N. OLIVARES, AND R. PERRUSSEL Abstract. Simulation of wave propagation through complex media relies on proper understanding
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 informationUnified A Posteriori Error Control for all Nonstandard Finite Elements 1
Unified A Posteriori Error Control for all Nonstandard Finite Elements 1 Martin Eigel C. Carstensen, C. Löbhard, R.H.W. Hoppe Humboldt-Universität zu Berlin 19.05.2010 1 we know of Guidelines for Applicants
More 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 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 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 informationTHE COMPACT DISCONTINUOUS GALERKIN (CDG) METHOD FOR ELLIPTIC PROBLEMS
SIAM J. SCI. COMPUT. Vol. 30, No. 4, pp. 1806 1824 c 2008 Society for Industrial and Applied Mathematics THE COMPACT DISCONTINUOUS GALERKIN (CDG) METHOD FOR ELLIPTIC PROBLEMS J. PERAIRE AND P.-O. PERSSON
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 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 informationIMA Preprint Series # 2146
AN ANALYSIS OF THE MINIMAL DISSIPATION LOCAL DISCONTINUOUS GALERKIN METHOD FOR CONVECTION-DIFFUSION PROBLEMS By Bernardo Cockburn and Bo Dong IMA Preprint Series # 2146 ( November 2006 ) INSTITUTE FOR
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 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 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 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 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 informationFINITE ELEMENT APPROXIMATION OF STOKES-LIKE SYSTEMS WITH IMPLICIT CONSTITUTIVE RELATION
Proceedings of ALGORITMY pp. 9 3 FINITE ELEMENT APPROXIMATION OF STOKES-LIKE SYSTEMS WITH IMPLICIT CONSTITUTIVE RELATION JAN STEBEL Abstract. The paper deals with the numerical simulations of steady flows
More informationA SECOND ELASTICITY ELEMENT USING THE MATRIX BUBBLE WITH TIGHTENED STRESS SYMMETRY
A SECOND ELASTICITY ELEMENT USING THE MATRIX BUBBLE WITH TIGHTENED STRESS SYMMETRY J. GOPALAKRISHNAN AND J. GUZMÁN Abstract. We presented a family of finite elements that use a polynomial space augmented
More 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 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 informationMixed Discontinuous Galerkin Methods for Darcy Flow
Journal of Scientific Computing, Volumes 22 and 23, June 2005 ( 2005) DOI: 10.1007/s10915-004-4150-8 Mixed Discontinuous Galerkin Methods for Darcy Flow F. Brezzi, 1,2 T. J. R. Hughes, 3 L. D. Marini,
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 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 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 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 informationGeneral linear methods for time-dependent PDEs
General linear methods for time-dependent PDEs Alexander Jaust and Jochen Schütz UHasselt Computational Mathematics Preprint Nr. UP-17-01 January 19, 2017 General linear methods for time-dependent PDEs
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 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 informationHybrid (DG) Methods for the Helmholtz Equation
Hybrid (DG) Methods for the Helmholtz Equation Joachim Schöberl Computational Mathematics in Engineering Institute for Analysis and Scientific Computing Vienna University of Technology Contributions by
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 gradient recovery method based on an oblique projection and boundary modification
ANZIAM J. 58 (CTAC2016) pp.c34 C45, 2017 C34 A gradient recovery method based on an oblique projection and boundary modification M. Ilyas 1 B. P. Lamichhane 2 M. H. Meylan 3 (Received 24 January 2017;
More informationPRECONDITIONING OF DISCONTINUOUS GALERKIN METHODS FOR SECOND ORDER ELLIPTIC PROBLEMS. A Dissertation VESELIN ASENOV DOBREV
PRECONDITIONING OF DISCONTINUOUS GALERKIN METHODS FOR SECOND ORDER ELLIPTIC PROBLEMS A Dissertation by VESELIN ASENOV DOBREV Submitted to the Office of Graduate Studies of Texas A&M University in partial
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 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 informationStochastic multiscale modeling of subsurface and surface flows. Part III: Multiscale mortar finite elements for coupled Stokes-Darcy flows
Stochastic multiscale modeling of subsurface and surface flows. Part III: Multiscale mortar finite elements for coupled Stokes-Darcy flows Ivan otov Department of Mathematics, University of Pittsburgh
More informationInterior superconvergence in mortar and non-mortar mixed finite element methods on non-matching grids
Interior superconvergence in mortar and non-mortar mixed finite element methods on non-matching grids Gergina Pencheva a,, Ivan Yotov b a Center for Subsurface Modeling, Institute for Computational Engineering
More informationNONCONFORMING MIXED ELEMENTS FOR ELASTICITY
Mathematical Models and Methods in Applied Sciences Vol. 13, No. 3 (2003) 295 307 c World Scientific Publishing Company NONCONFORMING MIXED ELEMENTS FOR ELASTICITY DOUGLAS N. ARNOLD Institute for Mathematics
More 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 informationEQUADIFF 6. Jean-Claude Nédélec Mixed finite element in 3D in H(div) and H(curl) Terms of use:
EQUADIFF 6 Jean-Claude Nédélec Mixed finite element in 3D in H(div) and H(curl) In: Jaromír Vosmanský and Miloš Zlámal (eds.): Equadiff 6, Proceedings of the International Conference on Differential Equations
More 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 informationFourier analysis for discontinuous Galerkin and related methods. Abstract
Fourier analysis for discontinuous Galerkin and related methods Mengping Zhang and Chi-Wang Shu Abstract In this paper we review a series of recent work on using a Fourier analysis technique to study the
More informationAdaptive Finite Elements with Thin Obstacles
International Mathematical Forum, Vol. 8, 2013, no. 39, 1921-1931 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/imf.2013.310194 Adaptive Finite Elements with Thin Obstacles J. Frohne and F.T.
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 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 informationFIRST-ORDER SYSTEM LEAST SQUARES (FOSLS) FOR PLANAR LINEAR ELASTICITY: PURE TRACTION PROBLEM
SIAM J. NUMER. ANAL. c 1998 Society for Industrial Applied Mathematics Vol. 35, No. 1, pp. 320 335, February 1998 016 FIRST-ORDER SYSTEM LEAST SQUARES (FOSLS) FOR PLANAR LINEAR ELASTICITY: PURE TRACTION
More information