Virtual Element Methods for general second order elliptic problems
|
|
- Jordan Blair
- 6 years ago
- Views:
Transcription
1 Virtual Element Methods for general second order elliptic problems Alessandro Russo Department of Mathematics and Applications University of Milano-Bicocca, Milan, Italy and IMATI-CNR, Pavia, Italy Workshop on Polytopal Element Methods in Mathematics and Engineering October 25 Georgia Tech, Atlanta A. Russo (Milan) VEM October 26, 25 / 27
2 Outline of the presentation Outline of the presentation Virtual Element Spaces Projectors VEM approximation of general elliptic equations Primal form Mixed form Numerical Experiments Conclusions Joint work with L. Beirão da Veiga, F. Brezzi, D. Marini A. Russo (Milan) VEM October 26, 25 2 / 27
3 Elliptic problem in primal form Elliptic problem in primal form We want to approximate a general second-order elliptic equation in 3D with an arbitrary polyhedral mesh with a conforming finite element method of order k. { div (κ u) + β u + αu = f in Ω R 3 u = g on Ω Beirão da Veiga, Brezzi, Marini, R.: Virtual Element Methods for general second order elliptic problems on polygonal meshes, to appear in Mathematical Models and Methods in Applied Sciences (also arxiv: ) We start with the two-dimensional case. A. Russo (Milan) VEM October 26, 25 3 / 27
4 Virtual Element spaces The local finite element space V k (E) Let E be a polygon. We would like to define a finite element space V k (E) on E such that: V k (E) contains the space P k (E) of polynomials of degree less than or equal to k (plus other non-polynomial functions); a function in V k (E) restricted to an edge e is in P k (e), so that if two polygons E and E have an edge in common, the two spaces V k (E) and V k (E ) must glue in C (E E ); given a function v h V k (E), I can compute its L 2 projection Π k v h in P k (E), and this is all the information of v h that I need in order to approximate my equation. A. Russo (Milan) VEM October 26, 25 4 / 27
5 Virtual Element spaces Meaning of I can compute In what follows the precise meaning of the statement is: I can compute Π k v h given the array dof i (v h ), I can directly compute Π k v h The same applies for all other quantities which are computable from v h. A. Russo (Milan) VEM October 26, 25 5 / 27
6 Virtual Element spaces The space V k,k (E) Given a polygon E, consider the following space (already mentioned in Franco s talk): V k,k (E) := {v h C (E) such that: v h e P k (e), v h P k (E)} If the polygon E has N E edges (and vertices), it is clear that dim V k,k (E) = k N E + dim P k (E) = k N E + Is the space V k,k (E) OK for us? (k + )(k + 2) 2 A. Russo (Milan) VEM October 26, 25 6 / 27
7 Virtual Element spaces Degrees of freedom in V k,k (E) Let (x E, y E ) be the centroid of E and h E its diameter. If α = (α, α 2 ) is a multiindex we define the scaled monomials of degree α = α + α 2 : ( ) x α ( ) xe y α2 ye m α (x, y) :=. The set {m α, with α k} is a basis for P k (E). h E As degrees of freedom in V k,k (E), we choose: the value of v h at the vertices and at k equally spaced points on each edge; the (scaled) moments v h m α for α k. E E h E A. Russo (Milan) VEM October 26, 25 7 / 27
8 Virtual Element spaces Degrees of freedom in V k,k (E) It can be easily shown that: the degrees of freedom above are unisolvent in V k,k (E). Short proof: suppose that v h V k,k (E) is zero on the boundary and v h m α = for α k We show that v h. For, v h 2 = E Hence v h. E E v h v h = since v h P k (E) A. Russo (Milan) VEM October 26, 25 8 / 27
9 Virtual Element spaces Degrees of freedom in V k,k (E) The choice of the scaled moments v h m α among the degrees of E E freedom implies that, starting from the degrees of freedom of v h, I can trivially compute Π k v h := L 2 projection of v h onto P k (E). It is clear that if I could compute directly the bilinear form on the space V k,k (E), the resulting finite element method would converge with the right optimal rates. For the time being, assume that the L 2 projection onto P k (E) is enough (we have a general way of approximating the bilinear form just using the L 2 projection). A. Russo (Milan) VEM October 26, 25 9 / 27
10 Virtual Element spaces Basis functions in V k,k (E) For i =,..., N dof we define ϕ i as the function in V k,k (E) such that { if i = j dof j (ϕ i ) = j-th degree of freedom of ϕ i = if i j We have the usual Lagrange-type expansion v h = dof i (v h ) ϕ i. N dof i= A. Russo (Milan) VEM October 26, 25 / 27
11 Virtual Element spaces Is the space V k,k (E) OK for the VEM?... yes, but IT HAS TOO MANY DEGREES OF FREEDOM! We can satisfy the same properties with a (proper) subspace of V k,k (E). Original VEM: define a projector Π k : V k,k(e) P k (E) which uses only the moments up to degree k 2; identify a subspace V k (E) V k,k (E) requiring that v h m α = Π k v h m α for α = k, k E E the dimension of this subspace is k N E + dim P k 2 (E) and the internal degrees of freedom are the (scaled) moments up to k 2; in this subspace we can compute Π k v h. A. Russo (Milan) VEM October 26, 25 / 27
12 Virtual Element spaces Is the space V k,k (E) OK for the VEM? Serendipity VEM: define a projector Π S k : V k,k(e) P k (E) which uses only the moments up to degree k N E ; identify a subspace V k (E) V k,k (E) requiring that v h m α = Π S k v h m α for α = k N E +,..., k E E the dimension of this subspace is k N E + dim P k NE (E) and the internal degrees of freedom are the (scaled) moments up to k N E ; in this subspace we can compute Π k v h. A. Russo (Milan) VEM October 26, 25 2 / 27
13 Virtual Element spaces The local finite element space V k (E) Let E be a polygon. We have defined a finite element space V k (E) on E such that: V k (E) contains the space P k (E) of polynomials of degree less than or equal to k (plus other non-polynomial functions); a function in V k (E) restricted to an edge e is in P k (e), so that if two polygons E and E have an edge in common, the two spaces V k (E) and V k (E ) must glue in C (E E ); given a function v h V k (E), I can compute its L 2 projection Π k v h in P k (E), and this is all the information of v h that I need in order to approximate my equation. A. Russo (Milan) VEM October 26, 25 3 / 27
14 Virtual Element spaces Computing another projector in V k (E): Π k v h We show now that we can easily compute Π k v h. To compute Π k v h we need to know the moments of v h up to order k : v h x m m β β = v h E x + v h m β nx, β k E E and both terms are computable directly from the degrees of freedom of v h. A. Russo (Milan) VEM October 26, 25 4 / 27
15 VEM approximation of general elliptic equations VEM approximation of general elliptic equations We consider now a general second order elliptic operator with variable coefficients: div (κ u) + β u + αu = f and we approximate the various local consistency terms as: κ u h v h κ [ Π k u ] [ ] h Π k v h E E ( ) ( [ ]) β uh vh β Π k u h Π k v h E E α u h v h and for the right-hand-side: f v h E E E E α [ Π k u ] [ ] h Π k v h f Π k v h ( Π k 2 v h is enough for k 2 ) A. Russo (Milan) VEM October 26, 25 5 / 27
16 VEM approximation of general elliptic equations VEM approximation of general elliptic equations The approximations above produce, as usual, rank-deficient matrices that must be stabilized. For the stabilization we can take the following term: with S(ϕ i, ϕ j ) = δ ij. S((I Π k )u h, (I Π k )v h) This means that we expand Π k ϕ i in terms of the ϕ j s and then we replace S(ϕ i, ϕ j ) with the identity matrix. A. Russo (Milan) VEM October 26, 25 6 / 27
17 US mesh 49 polygons, 5686 vertices A. Russo (Milan) VEM October 26, 25 7 / 27
18 US mesh vertices GEORGIA A. Russo (Milan) VEM October 26, 25 8 / 27
19 US mesh VEM solution on skeleton exact solution = 4 th degree polynomial...2 L 2 error k = 2:.6 3 k = 3:.43 4 k = 4: (patch test) A. Russo (Milan) VEM October 26, 25 9 / 27
20 Mixed VEM The continuous problem in mixed form Ω R 2 (polygonal) computational domain, κ, γ, b smooth. Assume that the problem { L p := div( κ(x) p + b(x)p) + γ(x) p = f (x) in Ω p = on Γ is solvable for any f H (Ω), and p,ω C f,ω, p 2,Ω C f,ω Existence and uniqueness as well for the adjoint operator L p := div( κ(x) p) b(x) p + γ(x) p In particular, f L 2 (Ω) ϕ H 2 (Ω) H (Ω): L ϕ = f, ϕ 2,Ω C f,ω A. Russo (Milan) VEM October 26, 25 2 / 27
21 Mixed VEM The continuous problem in mixed form Setting: ν = κ, β = κ b and νu := p + βp, we have div u + γ p = f in Ω, p = on Γ V = H(div, Ω), Q = L 2 (Ω) Find (u, p) V Q such that (νu, v) (p, div v) (β v, p) = v V (div u, q) + (γp, q) = (f, q) q Q A. Russo (Milan) VEM October 26, 25 2 / 27
22 Mixed VEM VEM approximation: the discrete spaces The spaces: We define, for k integer, Vh k (E) := {v H(div; E) H(rot; E) : v n e P k (e) e E, div v P k (E), and rot v P k (E)}. Vh k := {v H(div; Ω) such that v E Vh k (E) E T h} Qh k := {q L2 (Ω) such that: q E P k (E) element E int h } degrees of freedom in Vh k : v n p k ds e, p k P k (e) e v g k dx E, g k G k (E), E v g k dx E, g k G k (E) E (G k = P k, G k (E) = L2 (E) orthogonal of G k (E) in (P k (E)) 2 ) A. Russo (Milan) VEM October 26, / 27
23 Mixed VEM The bilinear forms a E h (v, w) := (νπ k v, Π k w),e + S E (v Π k v, w Π k w) Set α a E (v, v) S E (v, v) α a E (v, v) a h (v, w) := E ah E (v, w). v V k h Lemma The bilinear form a h (, ) is continuous and elliptic in (L 2 ) 2 : M > such that a h (v, w) M v w v, w Vh k, α > such that a h (v, v) α v 2 v Vh k, with M and α depending on ν but independent of h. A. Russo (Milan) VEM October 26, / 27
24 Mixed VEM The discrete problem Find (u h, p h ) Vh k Qk h such that ( ) a h (u h, v h ) (p h, div v h ) (β Π k v h, p h ) = (div u h, q h ) + (γp h, q h ) = (f, q h ) q h Q h. v h V k h Theorem Assume that the continuous problem has a unique solution p. Then, for h sufficiently small, ( ) has a unique solution (u h, p h ) V k h Qk h, and p p h Ch k+( u k+ + p k+ ), u u h Ch k+( u k+ + p k+ ), div(u u h ) Ch k+( f k+ + p k+ ) with C a constant depending on ν, β, and γ but independent of h. A. Russo (Milan) VEM October 26, / 27
25 Mixed VEM Superconvergence results Theorem Let p h be the discrete solution, and let p I Qh k Then, for h sufficiently small, be the interpolant of p. p I p h C h k+2( u k+ + p k+ + f k+ ), where C is a constant depending on ν, β, and γ but independent of h. A. Russo (Milan) VEM October 26, / 27
26 Numerical Experiments Joining meshes joined 832 polygons, Ndof=6849, h max =.77e, h mean = 4.8e 2 VEM solution A. Russo (Milan) VEM July 4, / 43
27 Numerical Experiments Joining meshes error A. Russo (Milan) VEM July 4, / 43
28 Numerical Experiments Joining meshes 2 joined 852 polygons, Ndof=733, h max =.77e, h mean = 4.3e 2 VEM solution A. Russo (Milan) VEM July 4, / 43
29 Numerical Experiments Joining meshes VEM solution A. Russo (Milan) VEM July 4, / 43
30 Numerical Experiments Joining meshes 3 joined 573 polygons, Ndof=389, h max = 3.73e, h mean = 3.46e 2 A. Russo (Milan) VEM July 4, / 43
31 Numerical Experiments Joining meshes VEM solution A. Russo (Milan) VEM July 4, / 43
32 Numerical results div( κ p + bp) + γp = f in [, ] [, ] κ(x, y) = ( x 2 + y 2 xy ) xy x 2 + y 2 b(x, y) = [x, y] γ(x, y) = x 2 + y 2 Exact solution: p(x, y) = x 2 + y 2 + cos(7y) sin(3x) + 2, u = κ p + bp Donatella Marini (Pavia) VEM variable Ferrara 25 / 23
33 exact solution x 2 +y 2 +cos(y 7.) sin(x 3.) y x.6.8
34 Errors: uniform triangulation, k = k = 2 k= u u h,ω p p I,Ω p I p h,ω Donatella Marini (Pavia) VEM variable Ferrara 25 3 / 23
35 Errors: uniform triangulation, k = 3 3 k = 3 3 k=3 u u h,ω p p I,Ω p I p h,ω Donatella Marini (Pavia) VEM variable Ferrara 25 4 / 23
36 Voronoi mesh, k = lloyd Donatella Marini (Pavia) VEM variable Ferrara 25 5 / 23
37 Voronoi mesh, k = k = k= u u h,ω p p I,Ω p I p h,ω Donatella Marini (Pavia) VEM variable Ferrara 25 5 / 23
38 Voronoi mesh, k = 3 2 k = 3 2 k=3 u u h,ω p p I,Ω p I p h,ω Donatella Marini (Pavia) VEM variable Ferrara 25 6 / 23
39 Poisson problem; load= Dirac mass at center Finite elements of degree 4: mesh and solution Generated by conn2mesh on 2 Mar 25 8:45: Donatella Marini (Pavia) VEM variable Ferrara 25 7 / 23
40 Poisson problem; load= Dirac mass at center Virtual elements of degree 4: mesh and solution joined Donatella Marini (Pavia) VEM variable Ferrara 25 8 / 23
41 Advection dominated case ε p + β p = in [, ] [, ] p = ( cos(2πy)/2 on x =, p = on x = p = on y =, y = n Donatella Marini (Pavia) VEM variable Ferrara 25 9 / 23
42 Advection dominated case ε p + β p = in [, ] [, ] p = ( cos(2πy)/2 on x =, p = on x = p = on y =, y = n β = (, /3). Test case: ε = 2, Virtual elements of degree : mesh and solution Donatella Marini (Pavia) VEM variable Ferrara 25 9 / 23
43 Advection dominated case ε p + β p = in [, ] [, ] p = ( cos(2πy)/2 on x =, p = on x = p = on y =, y = n β = (, /3). Test case: ε = 2, Virtual elements of degree : mesh and solution joined VEM solution k= Donatella Marini (Pavia) VEM variable Ferrara 25 9 / 23
44 Advection dominated case joined Donatella Marini (Pavia) VEM variable Ferrara 25 2 / 23
45 Advection dominated case VEM solution k= Donatella Marini (Pavia) VEM variable Ferrara 25 2 / 23
46 Advection dominated case joined Donatella Marini (Pavia) VEM variable Ferrara 25 2 / 23
47 Advection dominated case VEM solution k= Donatella Marini (Pavia) VEM variable Ferrara 25 2 / 23
48 Advection dominated case joined Donatella Marini (Pavia) VEM variable Ferrara / 23
49 Advection dominated case VEM solution k= Donatella Marini (Pavia) VEM variable Ferrara / 23
50 The VEM paradigm The VEM paradigm We believe that the Virtual Element Method has a wide range of applicability. The keypoints of VEM are: The definition of the local finite element spaces and of the associated degrees of freedom. These spaces contain polynomials plus other functions which are not computable. The costruction of variuos projectors onto polynomial spaces. The definition of a consistent and stable bilinear form using these projectors which is the VEM approximation of the exact bilinear form. A general theorem that guarantees convergence for a consistent and stable approximate bilinear form. A. Russo (Milan) VEM October 26, / 27
51 Basic references Basic References Virtual Element Methods for general second order elliptic problems on polygonal meshes, Beirão da Veiga, Brezzi, Marini, R.: to appear in Math. Models and Methods in Applied Sciences (also arxiv: [math.na]) Mixed Virtual Element Methods for general second order elliptic problems on polygonal meshes, Beirão da Veiga, Brezzi, Marini, R.: to appear in ESAIM: Math. Modelling and Numerical Analysis (also arxiv: [math.na]) The Hitchhiker s Guide to the Virtual Element Method Beirão da Veiga, Brezzi, Marini, R., Math. Models and Methods in Applied Sciences Vol. 24, No. 8 (24), pp Thanks for your attention! alessandro.russo@unimib.it A. Russo (Milan) VEM October 26, / 27
The 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 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 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 informationH(div) and H(curl)-conforming Virtual Element Methods
Noname manuscript No. (will be inserted by the editor) H(div) and H(curl)-conforming Virtual Element Methods L. Beirão da Veiga F. Brezzi L. D. Marini A. Russo the date of receipt and acceptance should
More informationA C 1 virtual element method for the Cahn-Hilliard equation with polygonal meshes. Marco Verani
A C 1 virtual element method for the Cahn-Hilliard equation with polygonal meshes Marco Verani MOX, Department of Mathematics, Politecnico di Milano Joint work with: P. F. Antonietti (MOX Politecnico di
More informationFebruary 18, :41 WSPC/INSTRUCTION FILE Brezzi-Corrections VIRTUAL ELEMENT METHOD FOR GENERAL SECOND ORDER ELLIPTIC PROBLEMS ON POLYGONAL MESHES
February 8, 206 7:4 WSPC/INSTRUCTION FIL Brezzi-Corrections Mathematical Models and Methods in Applied Sciences c World Scientific Publishing Company VIRTUAL LMNT MTHOD FOR GNRAL SCOND ORDR LLIPTIC PROBLMS
More informationVEM s for the numerical solution of PDE s. F. Brezzi
VEM s for the numerical solution of PDE s F. Brezzi IUSS-Istituto Universitario di Studi Superiori, Pavia, Italy IMATI-C.N.R., Pavia, Italy From Joint works with L. Beirão da Veiga, L.D. Marini, and A.
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 informationDivergence free Virtual Elements for the Stokes problem on polygonal meshes
Divergence free Virtual Elements for the Stokes problem on polygonal meshes L. Beirão da Veiga, C. Lovadina and G. Vacca Abstract arxiv:1510.01655v1 [math.na] 6 Oct 2015 In the present paper we develop
More informationarxiv: v1 [math.na] 11 Sep 2017
Dual virtual element method in presence of an inclusion Alessio Fumagalli arxiv:1709.03519v1 [math.na] 11 Sep 2017 Abstract We consider a Darcy problem for saturated porous media written in dual formulation
More informationTHE NONCONFORMING VIRTUAL ELEMENT METHOD
THE NONCONFORMING VIRTUAL ELEMENT METHOD BLANCA AYUSO DE DIOS, ONSTANTIN LIPNIOV, AND GIANMARCO MANZINI Abstract. We introduce the nonconforming Virtual Element Method (VEM) for the approximation of second
More informationINSTITUTE OF MATHEMATICS THE CZECH ACADEMY OF SCIENCES
INSTITUTE OF MATHEMATICS THE CZECH ACADEMY OF SCIENCES On the approximation of a virtual coarse space for domain decomposition methods in two dimensions Juan Gabriel Calvo Preprint No. 31-2017 PRAHA 2017
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 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 informationMimetic Finite Difference methods
Mimetic Finite Difference methods An introduction Andrea Cangiani Università di Roma La Sapienza Seminario di Modellistica Differenziale Numerica 2 dicembre 2008 Andrea Cangiani (IAC CNR) mimetic finite
More informationFinite Element Multigrid Framework for Mimetic Finite Difference Discretizations
Finite Element Multigrid Framework for Mimetic Finite ifference iscretizations Xiaozhe Hu Tufts University Polytopal Element Methods in Mathematics and Engineering, October 26-28, 2015 Joint work with:
More informationVirtual Elements for plates and shells
Virtual Elements for plates and shells Claudia Chinosi Dipartimento di Scienze e Innovazione Tecnologica, Università del Piemonte Orientale, Alessandria, Italia E-mail: claudia.chinosi@unipmn.it eywords:
More informationH(div) and H(curl)-conforming VEM
H(div) and H(curl)-conforming VM L. Beirão da Veiga, F. Brezzi, L. D. Marini, and A. Russo arxiv:1407.6822v1 [math.na] 25 Jul 2014 Abstract In the present paper we construct Virtual lement Spaces that
More informationQUADRILATERAL H(DIV) FINITE ELEMENTS
QUADRILATERAL H(DIV) FINITE ELEMENTS DOUGLAS N. ARNOLD, DANIELE BOFFI, AND RICHARD S. FALK Abstract. We consider the approximation properties of quadrilateral finite element spaces of vector fields defined
More informationICES REPORT A Generalized Mimetic Finite Difference Method and Two-Point Flux Schemes over Voronoi Diagrams
ICS RPORT 15-17 July 2015 A Generalized Mimetic Finite Difference Method and Two-Point Flux Schemes over Voronoi Diagrams by Omar Al-Hinai, Mary F. Wheeler, Ivan Yotov The Institute for Computational ngineering
More informationAn introduction to Isogeometric Analysis
An introduction to Isogeometric Analysis A. Buffa 1, G. Sangalli 1,2, R. Vázquez 1 1 Istituto di Matematica Applicata e Tecnologie Informatiche - Pavia Consiglio Nazionale delle Ricerche 2 Dipartimento
More informationSerendipity Finite Element Methods for Cardiac Electrophysiology Models
Serendipity Finite Element Methods for Cardiac Electrophysiology Models Andrew Gillette Department of Mathematics University of California, San Diego http://ccom.ucsd.edu/ agillette/ Andrew Gillette -
More 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 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 informationarxiv: v1 [math.na] 8 Feb 2018
arxiv:180.094v1 [math.na] 8 Feb 018 The nonconforming virtual element method for eigenvalue problems F. Gardini a, G. Manzini b, and G. Vacca c a Dipartimento di Matematica F. Casorati, Università di Pavia,
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 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 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 informationBASIC PRINCIPLES OF VIRTUAL ELEMENT METHODS
Mathematical Models and Methods in Applied Sciences Vol. 23, No. 1 (2013) 199 214 c World Scientific Publishing Company DOI: 10.1142/S0218202512500492 BASIC PRINCIPLES OF VIRTUAL ELEMENT METHODS L. BEIRÃO
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 informationSome New Elements for the Reissner Mindlin Plate Model
Boundary Value Problems for Partial Differential Equations and Applications, J.-L. Lions and C. Baiocchi, eds., Masson, 1993, pp. 287 292. Some New Elements for the Reissner Mindlin Plate Model Douglas
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University The Residual and Error of Finite Element Solutions Mixed BVP of Poisson Equation
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Nonconformity and the Consistency Error First Strang Lemma Abstract Error Estimate
More informationarxiv: v1 [math.na] 3 Jun 2018
SUPG stabilization for the nonconforming virtual element method for advection-diffusion-reaction equations S Berrone a, A Borio a, G Manzini b,c arxiv:186879v1 [mathna] 3 Jun 18 Abstract a Dipartimento
More informationUC Irvine UC Irvine Previously Published Works
UC Irvine UC Irvine Previously Published Works Title Some error analysis on virtual element methods Permalink https://escholarship.org/uc/item/1ph610x7 Journal CALCOLO, 55(1) ISSN 0008-0624 Authors Chen,
More informationDiscontinuous Petrov-Galerkin Methods
Discontinuous Petrov-Galerkin Methods Friederike Hellwig 1st CENTRAL School on Analysis and Numerics for Partial Differential Equations, November 12, 2015 Motivation discontinuous Petrov-Galerkin (dpg)
More informationHamburger Beiträge zur Angewandten Mathematik
Hamburger Beiträge zur Angewandten Mathematik Numerical analysis of a control and state constrained elliptic control problem with piecewise constant control approximations Klaus Deckelnick and Michael
More informationarxiv: v1 [math.na] 3 Jun 2016
SRNDIPITY FAC AND DG VM SPACS L. BIRÃO DA VIGA[1], F. BRZZI [2], L.D. MARINI [3], A. RUSSO [1] arxiv:1606.01048v1 [math.na] 3 Jun 2016 Abstract. We extend the basic idea of Serendipity Virtual lements
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 informationProjected Surface Finite Elements for Elliptic Equations
Available at http://pvamu.edu/aam Appl. Appl. Math. IN: 1932-9466 Vol. 8, Issue 1 (June 2013), pp. 16 33 Applications and Applied Mathematics: An International Journal (AAM) Projected urface Finite Elements
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 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 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 informationA Locking-Free MHM Method for Elasticity
Trabalho apresentado no CNMAC, Gramado - RS, 2016. Proceeding Series of the Brazilian Society of Computational and Applied Mathematics A Locking-Free MHM Method for Elasticity Weslley S. Pereira 1 Frédéric
More informationMIXED FINITE ELEMENT METHODS FOR 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 informationA new approach for Kirchhoff-Love plates and shells
A new approach for Kirchhoff-Love plates and shells Walter Zulehner Institute of Computational Mathematics JKU Linz, Austria AANMPDE 10 October 2-6, 2017, Paleochora, Crete, Greece Walter Zulehner (JKU
More informationINSTITUTE OF MATHEMATICS THE CZECH ACADEMY OF SCIENCES. A virtual overlapping Schwarz method for scalar elliptic problems in two dimensions
INSTITUTE OF MATHEMATICS THE CZECH ACADEMY OF SCIENCES A virtual overlapping Schwarz method for scalar elliptic problems in two dimensions Juan Gabriel Calvo Preprint No. 25-2017 PRAHA 2017 A VIRTUAL
More informationIntroduction. J.M. Burgers Center Graduate Course CFD I January Least-Squares Spectral Element Methods
Introduction In this workshop we will introduce you to the least-squares spectral element method. As you can see from the lecture notes, this method is a combination of the weak formulation derived from
More informationarxiv: v1 [math.na] 8 Jan 2019
arallel solvers for virtual element discretizations of elliptic equations in mixed form F. Dassi a, S. Scacchi b a Dipartimento di Matematica, Università degli Studi di Milano-Bicocca, Via Cozzi, 2133
More informationIsogeometric Analysis:
Isogeometric Analysis: some approximation estimates for NURBS L. Beirao da Veiga, A. Buffa, Judith Rivas, G. Sangalli Euskadi-Kyushu 2011 Workshop on Applied Mathematics BCAM, March t0th, 2011 Outline
More informationUniform inf-sup condition for the Brinkman problem in highly heterogeneous media
Uniform inf-sup condition for the Brinkman problem in highly heterogeneous media Raytcho Lazarov & Aziz Takhirov Texas A&M May 3-4, 2016 R. Lazarov & A.T. (Texas A&M) Brinkman May 3-4, 2016 1 / 30 Outline
More informationABHELSINKI UNIVERSITY OF TECHNOLOGY
ABHELSINKI UNIVERSITY OF TECHNOLOGY TECHNISCHE UNIVERSITÄT HELSINKI UNIVERSITE DE TECHNOLOGIE D HELSINKI A posteriori error analysis for Kirchhoff plate elements Jarkko Niiranen Laboratory of Structural
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 informationFinite Elements. Colin Cotter. January 18, Colin Cotter FEM
Finite Elements January 18, 2019 The finite element Given a triangulation T of a domain Ω, finite element spaces are defined according to 1. the form the functions take (usually polynomial) when restricted
More informationStandard Finite Elements and Weighted Regularization
Standard Finite Elements and Weighted Regularization A Rehabilitation Martin COSTABEL & Monique DAUGE Institut de Recherche MAthématique de Rennes http://www.maths.univ-rennes1.fr/~dauge Slides of the
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 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 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 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 informationNumerical methods for PDEs FEM convergence, error estimates, piecewise polynomials
Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen Numerical methods for PDEs FEM convergence, error estimates, piecewise polynomials Dr. Noemi Friedman Contents of the course Fundamentals
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 informationAdaptive methods for control problems with finite-dimensional control space
Adaptive methods for control problems with finite-dimensional control space Saheed Akindeinde and Daniel Wachsmuth Johann Radon Institute for Computational and Applied Mathematics (RICAM) Austrian Academy
More informationMIXED FINITE ELEMENT APPROXIMATION OF THE VECTOR LAPLACIAN WITH DIRICHLET BOUNDARY CONDITIONS
MIXED FINITE ELEMENT APPROXIMATION OF THE VECTOR LAPLACIAN WITH DIRICHLET BOUNDARY CONDITIONS DOUGLAS N. ARNOLD, RICHARD S. FALK, AND JAY GOPALAKRISHNAN Abstract. We consider the finite element solution
More informationABHELSINKI UNIVERSITY OF TECHNOLOGY
ABHELSINKI UNIVERSITY OF TECHNOLOGY TECHNISCHE UNIVERSITÄT HELSINKI UNIVERSITE DE TECHNOLOGIE D HELSINKI A posteriori error analysis for the Morley plate element Jarkko Niiranen Department of Structural
More informationarxiv: v1 [math.na] 21 Dec 2018
The p- and hp-versions of the virtual element method for elliptic eigenvalue problems O. Čertík, F. Gardini, G. Manzini, L. Mascotto, G. Vacca arxiv:8.090v [math.na] Dec 08 Abstract We discuss the p- and
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 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 informationKey words. Incompressible magnetohydrodynamics, mixed finite element methods, discontinuous Galerkin methods
A MIXED DG METHOD FOR LINEARIZED INCOMPRESSIBLE MAGNETOHYDRODYNAMICS PAUL HOUSTON, DOMINIK SCHÖTZAU, AND XIAOXI WEI Journal of Scientific Computing, vol. 40, pp. 8 34, 009 Abstract. We introduce and analyze
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 informationNumerical methods for PDEs FEM convergence, error estimates, piecewise polynomials
Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen Numerical methods for PDEs FEM convergence, error estimates, piecewise polynomials Dr. Noemi Friedman Contents of the course Fundamentals
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 informationarxiv: v1 [math.na] 5 Jun 2018
PRIMAL-DUAL WEAK GALERKIN FINIE ELEMEN MEHODS FOR ELLIPIC CAUCHY PROBLEMS CHUNMEI WANG AND JUNPING WANG arxiv:1806.01583v1 [math.na] 5 Jun 2018 Abstract. he authors propose and analyze a well-posed numerical
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 informationIsogeometric mortaring
Isogeometric mortaring E. Brivadis, A. Buffa, B. Wohlmuth, L. Wunderlich IMATI E. Magenes - Pavia Technical University of Munich A. Buffa (IMATI-CNR Italy) IGA mortaring 1 / 29 1 Introduction Splines Approximation
More informationHEXAHEDRAL H(DIV) AND H(CURL) FINITE ELEMENTS
HEXAHEDRAL H(DIV) AND H(CURL) FINITE ELEMENTS RICHARD S. FAL, PAOLO GATTO, AND PETER MON Abstract. We study the approximation properties of some finite element subspaces of H(div; Ω) and H(curl ; Ω) defined
More informationSpace-time Finite Element Methods for Parabolic Evolution Problems
Space-time Finite Element Methods for Parabolic Evolution Problems with Variable Coefficients Ulrich Langer, Martin Neumüller, Andreas Schafelner Johannes Kepler University, Linz Doctoral Program Computational
More 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 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 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 informationRobust Domain Decomposition Preconditioners for Abstract Symmetric Positive Definite Bilinear Forms
www.oeaw.ac.at Robust Domain Decomposition Preconditioners for Abstract Symmetric Positive Definite Bilinear Forms Y. Efendiev, J. Galvis, R. Lazarov, J. Willems RICAM-Report 2011-05 www.ricam.oeaw.ac.at
More informationFinite Elements. Colin Cotter. January 15, Colin Cotter FEM
Finite Elements January 15, 2018 Why Can solve PDEs on complicated domains. Have flexibility to increase order of accuracy and match the numerics to the physics. has an elegant mathematical formulation
More informationSUPERCONVERGENCE PROPERTIES FOR OPTIMAL CONTROL PROBLEMS DISCRETIZED BY PIECEWISE LINEAR AND DISCONTINUOUS FUNCTIONS
SUPERCONVERGENCE PROPERTIES FOR OPTIMAL CONTROL PROBLEMS DISCRETIZED BY PIECEWISE LINEAR AND DISCONTINUOUS FUNCTIONS A. RÖSCH AND R. SIMON Abstract. An optimal control problem for an elliptic equation
More informationA 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 informationOmar Al-Hinai 1, Mary F. Wheeler 1 and Ivan Yotov 2
SAIM: M2AN 51 (2017) 679 706 DOI: 10.1051/m2an/2016033 SAIM: Mathematical Modelling and Numerical Analysis www.esaim-m2an.org A GNRALIZD MIMTIC FINIT DIFFRNC MTHOD AND TWO-POINT FLUX SCHMS OVR VORONOI
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 informationTrefftz-discontinuous Galerkin methods for time-harmonic wave problems
Trefftz-discontinuous Galerkin methods for time-harmonic wave problems Ilaria Perugia Dipartimento di Matematica - Università di Pavia (Italy) http://www-dimat.unipv.it/perugia Joint work with Ralf Hiptmair,
More informationarxiv: v1 [math.na] 5 Oct 2017
Lowest order Virtual Element approximation o magnetostatic problems L. Beirão da Veiga a,b, F. Brezzi b,, F. Dassi a, L.D. Marini c,b, A. Russo a,b a Dipartimento di Matematica e Applicazioni, Università
More informationLocal flux mimetic finite difference methods
Local flux mimetic finite difference methods Konstantin Lipnikov Mikhail Shashkov Ivan Yotov November 4, 2005 Abstract We develop a local flux mimetic finite difference method for second order elliptic
More 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 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 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 informationA Primal-Dual Weak Galerkin Finite Element Method for Second Order Elliptic Equations in Non-Divergence Form
A Primal-Dual Weak Galerkin Finite Element Method for Second Order Elliptic Equations in Non-Divergence Form Chunmei Wang Visiting Assistant Professor School of Mathematics Georgia Institute of Technology
More informationA WEAK GALERKIN MIXED FINITE ELEMENT METHOD FOR BIHARMONIC EQUATIONS
A WEAK GALERKIN MIXED FINITE ELEMENT METHOD FOR BIHARMONIC EQUATIONS LIN MU, JUNPING WANG, YANQIU WANG, AND XIU YE Abstract. This article introduces and analyzes a weak Galerkin mixed finite element method
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 informationGeometric Decomposition of Serendipity Finite Element Spaces
Geometric Decomposition of Serendipity Finite Element Spaces Andrew Gillette Department of Mathematics University of California, San Diego http://ccom.ucsd.edu/ agillette/ Joint Mathematics Meetings 2013
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 34: Improving the Condition Number of the Interpolation Matrix Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu
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 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 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 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 information