Mimetic Finite Difference methods
|
|
- Ellen Holland
- 5 years ago
- Views:
Transcription
1 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 difference methods 1 / 30
2 Some cronology... RECENTLY ( ) Shashkov-Steinberg, JCP 1995 (support-operator methods) (Hyman)-Shashkov-Steinberg, JCP 1996-(7) (rough diffusion) Hyman-Shashkov, SINUM 1999 (Maxwell) Campbell-Shashkov, JCP 2001 (gas dynamics) MORE RECENTLY ( ) Kuznetsov-Lipnikov-Shashkov, Comp. Geos (polygons) Brezzi-Lipnikov-Shashkov, SINUM 2005 (error analysis) Brezzi-Lipnikov-Simoncini, M 3 AS 2005 (a new family of MFD) Brezzi-Lipnikov-Shashkov, M 3 AS 2006 (curved faces) Andrea Cangiani (IAC CNR) mimetic finite difference methods 2 / 30
3 Some cronology (continued) MORE RECENTLY ( ) continued Beirao Da Veiga, NM 2007 (residual error estimator) C.-Manzini, CMAME 2008 (post-processing) C.-Manzini-Russo, SINUM accepted (convection-reaction-diffusion) IN PROGRESS Convection-dominated diffusion Higher-order MFD Nodal MFD Mimetic curl operator Mimetic discretization of Stokes Andrea Cangiani (IAC CNR) mimetic finite difference methods 3 / 30
4 Features RECENTLY Generalisation of finite differences to hexahedral meshes. Discrete differential operators defined so as to mimic the properties of the underlying continuum operators (e.g. vector calculus identities, conservation laws, solution symmetries) Applied to wide range of problems. MORE RECENTLY Generalisation of (low order) Mixed Finite Elements/finite volumes to general polyhedral meshes. Constructs a family of methods. Discrete differential operators defined so as to mimic the properties of the underlying continuum operators. May gain from extra freedom given by method construction Currently limited to linear diffusion problems Andrea Cangiani (IAC CNR) mimetic finite difference methods 4 / 30
5 Compatible Spatial Discretizations We define compatible spatial discretizations as those that inherit or mimic fundamental properties of the PDE such as topology, conservation, symmetries, and positivity structures and maximum principles. (IMA Hot Topics Workshop Compatible Spatial Discretizations for Partial Differential Equations, May 11-15, 2004) Advantages: Conserve crucial features of physical, geometrical, and mathematical model: Conservation laws Symmetry Positivity and monotonicity Duality properties of differential operators Provide reliability and accuracy Andrea Cangiani (IAC CNR) mimetic finite difference methods 5 / 30
6 Compatible Spatial Discretizations The following can be cast as compatible spatial discretizations: mixed finite element methods, mimetic finite differences, support operator methods, control volume methods, discrete differential forms, Whitney forms, conservative differencing, discrete Hodge operators, discrete Helmholtz decomposition, finite integration techniques, staggered grid and dual grid methods, etc. Support operator method: Most PDEs are written in terms of invariant differential operators (div, grad, and curl) Define discrete analogues of these invariant operators that satisfy exacly the discrete analogs of the identities satisfied by the continuum operators (e.g. Gauss, Stokes, Hodge orthogonal decomposition). Andrea Cangiani (IAC CNR) mimetic finite difference methods 6 / 30
7 Linear diffusion in mixed form Consider the linear diffusion equation bvp div(kgrad p) = b in Ω IR 2,3 p = 0 on Ω with K strongly elliptic full symmetric tensor. Mixed formulation: F = K grad p div F = b (Constitutive Equation) (Conservation Equation) We also introduce the flux operator Gp = Kgrad p. Andrea Cangiani (IAC CNR) mimetic finite difference methods 7 / 30
8 Linear diffusion in mixed form The operators div and grad satisfy the integral identity p divf dv + F gradp dv = 0 which we rewrite: p divf dv Ω Introducing the inner-products: Ω the identity (1) becomes: Ω Ω (p, q) W = Ω p q dv (F, G) V = Ω K 1 F G dv (divf, p) W (F, Gp) V = 0 K 1 F Gp dv = 0 (1) expressing the fact that divergence and flux are adjoint: G = grad Andrea Cangiani (IAC CNR) mimetic finite difference methods 8 / 30
9 Linear diffusion in mixed form Consider the linear diffusion equation bvp div(kgrad p) = b in Ω IR 2,3 p = 0 on Ω with K strongly elliptic full symmetric tensor. Mixed formulation: F = K grad p div F = b (Constitutive Equation) (Conservation Equation) Mixed variational formulation: find (p, F ) W V s.t. (K 1 F, G) (p, divg) = 0 G V (divf, q) = (b, q) q W W = L 2 (Ω), V = H(div; Ω) Andrea Cangiani (IAC CNR) mimetic finite difference methods 9 / 30
10 Mixed Finite Element RT0-P0 Discretisation Ω E Andrea Cangiani (IAC CNR) mimetic finite difference methods 10 / 30
11 Mixed Finite Element RT0-P0 Discretisation Let Ω h be a triangularization of Ω IR 2. For E Ω h, let P 0 (E) = Polynomials of degree zero on E RT 0 (E) = {a + bx : a (P 0 (E)) 2, b P 0 (E)} from which we define the conforming FE spaces: W h = {q h : q h E P 0 (E), E Ω h } V h = {G h : G h E RT 0 (E) and G h n cont. on Ω h } MFE RT0-P0 discretization: find (p h, F h ) W h V h s.t. (K 1 F h, G h ) (p h, divg h ) = 0 G h V h (divf h, q h ) = (b, q h ) q h W h Andrea Cangiani (IAC CNR) mimetic finite difference methods 11 / 30
12 RT0-P0 nodal variables Finite Element Nodal variables Lifting (E, P 0, N 0 ) N 0 (q h ) = q E L E 0 q E P 0 q E := 1 E E q LE 0 q E q E (E, RT 0, N RT ) N RT (G h ) = (GE e )e = G E L E RT G E RT 0 GE e := 1 e e G ne E L E RT G E e n e E = GE e Global FE space Global nodal spaces Global lifting W h Q h := {q = (q E ) E Ωh } L 0 : Q h W h { } G = (GE ) E Ωh : V h X h := GE e + GE e L RT : X h V h. + = 0 e Ω h Where E and E + are the two elements sharing the edge. Andrea Cangiani (IAC CNR) mimetic finite difference methods 12 / 30
13 Discrete divergence Let L RT G = G h. Then, E div G h E = E Gauss divg h = divl E RT G = divg h E = 1 E E G n E = e GE e e E e GE e =: div hg E e E We may now define the discrete divergence operator satisfying divl RT G = div h G. div h : X h Q h G div h G = (div h G E ) E (p h, div G h ) = E E p E div h G E = [p, div h G] Qh Andrea Cangiani (IAC CNR) mimetic finite difference methods 13 / 30
14 Scalar product in X h (K 1 F h, G h ) = Ω = E K 1 L RT F L RT G K 1 L E RT F E L E RT G E E = E [F E, G E ] E =: [F, G] Xh MFE in terms of the spaces Q h and X h : find (p, F) Q h X h : [F, G] Xh [p, div h G] Qh = 0 G X h [div h F, q] Qh = [b, q] Qh q Q h. with b = b I, I =interpolation operator defined for q V by (q I ) E := 1 q E Ω h E E Andrea Cangiani (IAC CNR) mimetic finite difference methods 14 / 30
15 Mimetic finite difference discretization Ω h partition of Ω into polygonal (polyhedral) elements. E E q e E G E Ω Why polyhedral meshes? They naturally arise in the treatment of complex solution domains and heterogeneous materials (e.g. reservoir models) They facilitate adaptive mesh refinement/de-refinement Andrea Cangiani (IAC CNR) mimetic finite difference methods 15 / 30
16 Discrete MFD pressure and flux spaces Define the discrete pressure space Q h and flux space X h. q Q h q = (q E ) E Ωh ( piecewise constant) ) e E G X h G = ( GE e ( constant normal component) E Ω h e E E + FE e + FE e + = 0 We have the interpolation operators: For q W, (q I ) E = 1 E E q E Ω h, ( I) e For G V, G E = 1 e e G ne E E Ω h e E And the discrete divergence operator: div h : X h Q h G div h G = (div h G E ) E div h G E = 1 E e GE e e E Andrea Cangiani (IAC CNR) mimetic finite difference methods 16 / 30
17 Discretisation of the conservation equation Let b = b I. The mimetic discrete conservation equation reads: i.e. Equivalently, 1 E div h F = b e F e E = (b) E E (2) e E [div h F, q] Qh = [b, q] Qh q Q h where the scalar product for the pressure space Q h is given by [p, q] Qh = E Ω h E p E q E This MFD discrete conservation equation (2) coincides with the MFE RT0-P0 and also with many finite volume discretisation of the conservation equation. Andrea Cangiani (IAC CNR) mimetic finite difference methods 17 / 30
18 Finite volume discretisation of the constitutive equation Directly discretize the constitutive equation defining the flux: F = K grad p (with K = k I) B B K L n AB K L A A n AB F k p K p L K L F k ( pk p L n AB + p ) B p A n KL K L γ AB B A γ KL Andrea Cangiani (IAC CNR) mimetic finite difference methods 18 / 30
19 Mimetic discretisation of the constitutive equation Introduce a scalar product for the flux space X h [F, G] Xh = [F, G] E. ([, ] E to be defined!). E Ω h Moreover, define a discrete flux operator G h : Q h X h as the adjoint of div h : [G, G h q] Xh = [q, div h G] Qh, q Q h G X h. This definition naturally establishes a discrete Green formula with respect to the discrete scalar products. We say that the discrete operators mimic the continuous ones. Andrea Cangiani (IAC CNR) mimetic finite difference methods 19 / 30
20 Discretisation of the constitutive equation (cont.) The mimetic discretisation of the constitutive equation reads: F = G h p or, equivalently, [F, G] Xh = [G h p, G] Xh =[p, div h G] Qh G X h The (family of) Mimetic Finite Difference (MFD) schemes: [F, G] Xh [p, div h G] Qh = 0 G X h [div h F, q] Qh = [b, q] Qh q Q h or, equivalently, [ A B t B 0 ] [ F p ] = [ 0 b ] Andrea Cangiani (IAC CNR) mimetic finite difference methods 20 / 30
21 MFE and MFD As we already saw, if Ω h is made of triangles (tetrahedrons), The MFD conservation equation discretisation is equivalent to RT0-P0 If we define the X h -scalar product to be [F, G] Xh = K 1 L RT F L RT G the constitutive equation disretisation is equivalent....and the two methods coincide! Which properties define L RT? Locally, on each E Ω h, 1 L E RT G e n e E = G e E e E, 2 divl E RT G = div hg E 3 L E RT G I = G if G is constant Ω Andrea Cangiani (IAC CNR) mimetic finite difference methods 21 / 30
22 The flux scalar product: P 0 -compatible liftings IDEA: any reasonable lifting can be used to define a scalar product yielding a reasonable MFD formulation. Set X E = X h E. We define P 0 -compatible lifting a linear map L E : X E L 2 (E) such that for all G X E, 1. L E G e n e E = G E e, e E, 2. divl E G = div h G E, and, for all constant vector C, 3. L E C I = C. The following defines a MFD scalar product for the flux space: [F, G] E = K 1 L E F L E G ds. Ω Andrea Cangiani (IAC CNR) mimetic finite difference methods 22 / 30
23 The flux scalar product: P 0 -compatible liftings Assume that K is constant on E. A crucial property. For any linear function q 1 and G X E [(K q 1 ) I, G] E = K 1 L E ( (K q 1 ) I ) L E G ds E 3 = q 1 L E G ds E Green = q 1 divl E G ds + q 1 (L E G n ext ) dl E E 2,1 = q 1 div h G ds + q 1 dl E e E G e E e (local consistency) the scalar product is independent from the lifting if one of the factors is the interpolant of a constant vector field. the scalar product is exact on the interpolant of constant vector fields: [(e i ) I, (e j ) I ] E = (K 1 ) i,j E Andrea Cangiani (IAC CNR) mimetic finite difference methods 23 / 30
24 A general locally consistent scalar product IDEA: Use local consistency to construct a scalar product without a lifting! Let k E be the number of edges (faces) of E. Set: [F, G] E = k E s,t=1 M s,t E F es E G et E Imposing local consistency and symmetry we get a set of positive semi-definite matrices Imposing stability: exists s, s > 0 independent of h s.t. k E s (G e i E )2 E [G, G] E s i=1 k E i=1 (G e i E )2 E G X E. we restrict to a set of acceptable (symmetric & pos. def.) matrices defining a family of reasonable (e.g. exact on constants) scalar products. Andrea Cangiani (IAC CNR) mimetic finite difference methods 24 / 30
25 A family of consistent and stable MFD scalar products [Brezzi-Lipnikov-Simoncini, Math. Models Methods Appl. Sci. (2005)] Imposing: local consistency ME 0 symmetric & semidef. stability M E = ME 0 +CUC t symmetric & pos. def., where C is a given (computable) full-rank k E (k E d) matrix, U is any symmetric & pos. def. (k E d) (k E d) matrix Proposition: Under a reasonable assumption on the minimum eigenvalue of U, the matrix M E is induced by a lifting. In practice (so far) we use: U = u ( E Trace(K 1 )) II ke d, u IR and the proposition requires that u is sufficiently large. Andrea Cangiani (IAC CNR) mimetic finite difference methods 25 / 30
26 MFD error in function of u pressure. -error pressure L2-error linear pressure L2-error flux. -error Quadrilaterals 20x20 BLS Example Andrea Cangiani (IAC CNR) mimetic finite difference methods 26 / 30
27 Back to Mixed finite elements C. & Manzini, CMAME, 2008 Example. If Ω h is made of triangles, M E = M 0 E + u E [ ei e j ] i,j. with M 0 E symmetric and positive semidefinite. x e x B u E = ke ( ) T i=1 xei x B K 1 ( ) x ei x B 12 E RT 0 P0 Example Rectangles central finite differences Andrea Cangiani (IAC CNR) mimetic finite difference methods 27 / 30
28 Assumptions on the partition The number of edges per face and the number of faces per element is unif. bounded. The length l of each edge is such that l Ch (shape regularity). Every element is uniformly strictly starshaped (if Ω IR 3 ) Every face is uniformly strictly starshaped Andrea Cangiani (IAC CNR) mimetic finite difference methods 28 / 30
29 Error estimates [Brezzi-Lipnikov-Shashkov, SIAM J. Num. Anal. (2005)] If Ω is a convex polyhedron (polygon) with Lipschitz continuous boundary, Ω h is a partition satisfying the previous slide assumptions, the scalar product satisfies local consistency and stability. p H 2 (Ω) Then, p I p Qh + F I F Xh C h, p L 0 p L 2 C h where 2 Q h = [, ] Qh and 2 X h = [, ] Xh. If, moreover, The scalar product is induced by (some) lifting. Then, p I p Qh C h 2 (pressure super convergence) Andrea Cangiani (IAC CNR) mimetic finite difference methods 29 / 30
Local 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 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 informationREMARKS ON DISCRETIZATIONS OF CONVECTION TERMS IN HYBRID MIMETIC MIXED METHODS
NETWORS AND HETEROGENEOUS MEDIA c American Institute of Mathematical Sciences Volume 5, Number 3, September 2010 doi:10.3934/nhm.2010.5.xx pp. 1 xx REMARS ON DISCRETIZATIONS OF CONVECTION TERMS IN HYBRID
More informationMimetic Finite Difference Methods
Mimetic Finite Difference Methods Applications to the Diffusion Equation Asbjørn Nilsen Riseth May 15, 2014 2014-05-15 Mimetic Finite Difference Methods Mimetic Finite Difference Methods Applications to
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 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 Element Methods for general second order elliptic problems
Virtual Element Methods for general second order elliptic problems Alessandro Russo Department of Mathematics and Applications University of Milano-Bicocca, Milan, Italy and IMATI-CNR, Pavia, Italy Workshop
More informationThe Virtual Element Method: an introduction with focus on fluid flows
The Virtual Element Method: an introduction with focus on fluid flows L. Beirão da Veiga Dipartimento di Matematica e Applicazioni Università di Milano-Bicocca VEM people (or group representatives): P.
More 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 informationA Mixed Nonconforming Finite Element for Linear Elasticity
A Mixed Nonconforming Finite Element for Linear Elasticity Zhiqiang Cai, 1 Xiu Ye 2 1 Department of Mathematics, Purdue University, West Lafayette, Indiana 47907-1395 2 Department of Mathematics and Statistics,
More informationarxiv: v1 [math.na] 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 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 informationarxiv: v1 [math-ph] 25 Apr 2013
Higher-order compatible discretization on hexahedrals Jasper Kreeft and Marc Gerritsma arxiv:1304.7018v1 [math-ph] 5 Apr 013 Abstract We derive a compatible discretization method that relies heavily on
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 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 informationSolving the curl-div system using divergence-free or curl-free finite elements
Solving the curl-div system using divergence-free or curl-free finite elements Alberto Valli Dipartimento di Matematica, Università di Trento, Italy or: Why I say to my students that divergence-free finite
More informationFrom the Boundary Element DDM to local Trefftz Finite Element Methods on Polyhedral Meshes
www.oeaw.ac.at From the Boundary Element DDM to local Trefftz Finite Element Methods on Polyhedral Meshes D. Copeland, U. Langer, D. Pusch RICAM-Report 2008-10 www.ricam.oeaw.ac.at From the Boundary Element
More informationFrom the Boundary Element Domain Decomposition Methods to Local Trefftz Finite Element Methods on Polyhedral Meshes
From the Boundary Element Domain Decomposition Methods to Local Trefftz Finite Element Methods on Polyhedral Meshes Dylan Copeland 1, Ulrich Langer 2, and David Pusch 3 1 Institute of Computational Mathematics,
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 informationDivergence-free or curl-free finite elements for solving the curl-div system
Divergence-free or curl-free finite elements for solving the curl-div system Alberto Valli Dipartimento di Matematica, Università di Trento, Italy Joint papers with: Ana Alonso Rodríguez Dipartimento di
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 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 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 informationThe Convergence of Mimetic Discretization
The Convergence of Mimetic Discretization for Rough Grids James M. Hyman Los Alamos National Laboratory T-7, MS-B84 Los Alamos NM 87545 and Stanly Steinberg Department of Mathematics and Statistics University
More informationAn Approach to Higher- Order Mimetic Finite Volume Schemes
An Approach to Higher- Order Mimetic Finite Volume Schemes Blair Perot Context Mimetic Methods Finite Dierence Mimetic SOM Finite Element Edge/Face Staggered Finite Volume Meshless Natural Neighbors DC
More informationLecture 2: Reconstruction and decomposition of vector fields on the sphere with applications
2013 Dolomites Research Week on Approximation : Reconstruction and decomposition of vector fields on the sphere with applications Grady B. Wright Boise State University What's the problem with vector fields
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 informationFinite volume method for nonlinear transmission problems
Finite volume method for nonlinear transmission problems Franck Boyer, Florence Hubert To cite this version: Franck Boyer, Florence Hubert. Finite volume method for nonlinear transmission problems. Proceedings
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 Multigrid Method for Two Dimensional Maxwell Interface Problems
A Multigrid Method for Two Dimensional Maxwell Interface Problems Susanne C. Brenner Department of Mathematics and Center for Computation & Technology Louisiana State University USA JSA 2013 Outline A
More 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 informationA Discrete Vector Calculus in Tensor Grids
COMPUTATIONAL METHODS IN APPLIED MATHEMATICS Vol. (00,, pp. 44 c 00 Institute of Mathematics, National Academy of Sciences A Discrete Vector Calculus in Tensor Grids Nicolas Robidoux Stanly Steinberg Abstract
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 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 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 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 informationWeighted Regularization of Maxwell Equations Computations in Curvilinear Polygons
Weighted Regularization of Maxwell Equations Computations in Curvilinear Polygons Martin Costabel, Monique Dauge, Daniel Martin and Gregory Vial IRMAR, Université de Rennes, Campus de Beaulieu, Rennes,
More informationA non-standard Finite Element Method based on boundary integral operators
A non-standard Finite Element Method based on boundary integral operators Clemens Hofreither Ulrich Langer Clemens Pechstein June 30, 2010 supported by Outline 1 Method description Motivation Variational
More informationInstitut de Recherche MAthématique de Rennes
LMS Durham Symposium: Computational methods for wave propagation in direct scattering. - July, Durham, UK The hp version of the Weighted Regularization Method for Maxwell Equations Martin COSTABEL & Monique
More informationOverlapping Schwarz preconditioners for Fekete spectral elements
Overlapping Schwarz preconditioners for Fekete spectral elements R. Pasquetti 1, L. F. Pavarino 2, F. Rapetti 1, and E. Zampieri 2 1 Laboratoire J.-A. Dieudonné, CNRS & Université de Nice et Sophia-Antipolis,
More informationFinite Elements. Colin Cotter. February 22, Colin Cotter FEM
Finite Elements February 22, 2019 In the previous sections, we introduced the concept of finite element spaces, which contain certain functions defined on a domain. Finite element spaces are examples of
More informationFinite element exterior calculus framework for geophysical fluid dynamics
Finite element exterior calculus framework for geophysical fluid dynamics Colin Cotter Department of Aeronautics Imperial College London Part of ongoing work on UK Gung-Ho Dynamical Core Project funded
More informationarxiv: v1 [math.na] 15 Mar 2015
A Finite Element Framework for Some Mimetic Finite Difference Discretizations arxiv:1503.04423v1 [math.na] 15 Mar 2015 C. Rodrigo a,, F.J. Gaspar a, X. Hu b, L. Zikatanov c,d a Applied Mathematics Department
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 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 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 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 informationFind (u,p;λ), with u 0 and λ R, such that u + p = λu in Ω, (2.1) div u = 0 in Ω, u = 0 on Γ.
A POSTERIORI ESTIMATES FOR THE STOKES EIGENVALUE PROBLEM CARLO LOVADINA, MIKKO LYLY, AND ROLF STENBERG Abstract. We consider the Stokes eigenvalue problem. For the eigenvalues we derive both upper and
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 informationFlux reconstruction and solution post-processing in mimetic finite difference methods
Available online at www.sciencedirect.com Comput. Methods Appl. Mech. ngrg. 197 (008) 933 945 www.elsevier.com/locate/cma Flux reconstruction and solution post-processing in mimetic finite difference methods
More informationINTERGRID OPERATORS FOR THE CELL CENTERED FINITE DIFFERENCE MULTIGRID ALGORITHM ON RECTANGULAR GRIDS. 1. Introduction
Trends in Mathematics Information Center for Mathematical Sciences Volume 9 Number 2 December 2006 Pages 0 INTERGRID OPERATORS FOR THE CELL CENTERED FINITE DIFFERENCE MULTIGRID ALGORITHM ON RECTANGULAR
More informationc 2013 Society for Industrial and Applied Mathematics
SIAM J. NUMER. ANAL. Vol. 51, No. 1, pp. 654 675 c 2013 Society for Industrial and Applied Mathematics HIERARCHICAL A POSTERIORI ERROR ESTIMATORS FOR THE MIMETIC DISCRETIZATION OF ELLIPTIC PROBLEMS PAOLA
More informationU M,j M,j. U 1,j. U i,j. U i,1 (1,N) (M,N) i,j (M,1) (M,1) i,1 (1,1) (1,N) i,n (M,N) i,j. i,1 (1,1) i,n. M,j+1/2 M,j. 1,j. i+1/2,n.
Adjoint Operators for the Natural Discretizations of the Divergence, Gradient, and Curl on Logically Rectangular Grids James M. Hyman and Mikhail Shashkov Los Alamos National Laboratory T-7, MS-B284, Los
More informationComparison of V-cycle Multigrid Method for Cell-centered Finite Difference on Triangular Meshes
Comparison of V-cycle Multigrid Method for Cell-centered Finite Difference on Triangular Meshes Do Y. Kwak, 1 JunS.Lee 1 Department of Mathematics, KAIST, Taejon 305-701, Korea Department of Mathematics,
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 posteriori error estimates for non conforming approximation of eigenvalue problems
A posteriori error estimates for non conforming approximation of eigenvalue problems E. Dari a, R. G. Durán b and C. Padra c, a Centro Atómico Bariloche, Comisión Nacional de Energía Atómica and CONICE,
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 informationA multipoint flux mixed finite element method on hexahedra
A multipoint flux mixed finite element method on hexahedra Ross Ingram Mary F. Wheeler Ivan Yotov Abstract We develop a mixed finite element method for elliptic problems on hexahedral grids that reduces
More informationA priori error analysis of the BEM with graded meshes for the electric eld integral equation on polyhedral surfaces
A priori error analysis of the BEM with graded meshes for the electric eld integral equation on polyhedral surfaces A. Bespalov S. Nicaise Abstract The Galerkin boundary element discretisations of the
More informationTowards Discrete Exterior Calculus and Discrete Mechanics for Numerical Relativity
Towards Discrete Exterior Calculus and Discrete Mechanics for Numerical Relativity Melvin Leok Mathematics, University of Michigan, Ann Arbor. Joint work with Mathieu Desbrun, Anil Hirani, and Jerrold
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 informationHigh order, finite volume method, flux conservation, finite element method
FLUX-CONSERVING FINITE ELEMENT METHODS SHANGYOU ZHANG, ZHIMIN ZHANG, AND QINGSONG ZOU Abstract. We analyze the flux conservation property of the finite element method. It is shown that the finite element
More 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 informationMultiplicable Discrete Operators In Finite Volume Solvers
Multiplicable Discrete Operators In Finite Volume Solvers DRAGAN VIDOVIĆ Institute For The Development Of Water Resources Jaroslav Černi Jaroslava Černog 80, 11226 Pinosava, Belgrade SERBIA draganvid@gmail.com
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 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 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 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: 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 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 informationFinite Element Methods for Maxwell Equations
CHAPTER 8 Finite Element Methods for Maxwell Equations The Maxwell equations comprise four first-order partial differential equations linking the fundamental electromagnetic quantities, the electric field
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 informationTime-Varying Systems; Maxwell s Equations
Time-Varying Systems; Maxwell s Equations 1. Faraday s law in differential form 2. Scalar and vector potentials; the Lorenz condition 3. Ampere s law with displacement current 4. Maxwell s equations 5.
More informationMimetic Methods and why they are better
Mimetic Methods and why they are better Blair Perot Theoretical and Computational Fluid Dynamics Laboratory February 25, 2015 Numerical Methods: Are more than Mathematics Math: Accuracy Stability Convergence
More informationLattice electromagnetic theory from a topological viewpoint a
JOURNAL OF MATHEMATICAL PHYSICS VOLUME 40, NUMBER 1 JANUARY 1999 Lattice electromagnetic theory from a topological viewpoint a F. L. Teixeira b) and W. C. Chew Center for Computational Electromagnetics,
More informationError analysis for a new mixed finite element method in 3D
Trabalho apresentado no CNMAC, Gramado - RS, 2016. Proceeding Series of the Brazilian Society of Computational and Applied Mathematics Error analysis for a new mixed finite element method in 3D Douglas
More informationThe Relationship between Discrete Calculus Methods and other Mimetic Approaches
The Relationship between Discrete Calculus Methods and other Mimetic Approaches Blair Perot Theoretical and Computational Fluid Dynamics Laboratory Woudschoten Conerence Oct. 2011 October 7, 2011 Background
More informationGeometric Multigrid Methods
Geometric Multigrid Methods Susanne C. Brenner Department of Mathematics and Center for Computation & Technology Louisiana State University IMA Tutorial: Fast Solution Techniques November 28, 2010 Ideas
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 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 informationContents as of 12/8/2017. Preface. 1. Overview...1
Contents as of 12/8/2017 Preface 1. Overview...1 1.1 Introduction...1 1.2 Finite element data...1 1.3 Matrix notation...3 1.4 Matrix partitions...8 1.5 Special finite element matrix notations...9 1.6 Finite
More informationA spacetime DPG method for acoustic waves
A spacetime DPG method for acoustic waves Rainer Schneckenleitner JKU Linz January 30, 2018 ( JKU Linz ) Spacetime DPG method January 30, 2018 1 / 41 Overview 1 The transient wave problem 2 The broken
More informationBrief Review of Vector Algebra
APPENDIX Brief Review of Vector Algebra A.0 Introduction Vector algebra is used extensively in computational mechanics. The student must thus understand the concepts associated with this subject. The current
More informationIntroduction to finite element exterior calculus
Introduction to finite element exterior calculus Ragnar Winther CMA, University of Oslo Norway Why finite element exterior calculus? Recall the de Rham complex on the form: R H 1 (Ω) grad H(curl, Ω) curl
More informationAn abstract Hodge Dirac operator and its stable discretization
An abstract Hodge Dirac operator and its stable discretization Paul Leopardi Mathematical Sciences Institute, Australian National University. For presentation at ICCA, Tartu,. Joint work with Ari Stern,
More informationA Posteriori Error Estimation Techniques for Finite Element Methods. Zhiqiang Cai Purdue University
A Posteriori Error Estimation Techniques for Finite Element Methods Zhiqiang Cai Purdue University Department of Mathematics, Purdue University Slide 1, March 16, 2017 Books Ainsworth & Oden, A posteriori
More informationA Multiscale Mortar Multipoint Flux Mixed Finite Element Method
A Multiscale Mortar Multipoint Flux Mixed Finite Element Method Mary F. Wheeler Guangri Xue Ivan Yotov Abstract In this paper, we develop a multiscale mortar multipoint flux mixed finite element method
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 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 posteriori error estimates applied to flow in a channel with corners
Mathematics and Computers in Simulation 61 (2003) 375 383 A posteriori error estimates applied to flow in a channel with corners Pavel Burda a,, Jaroslav Novotný b, Bedřich Sousedík a a Department of Mathematics,
More informationfield using second order edge elements in 3D
The current issue and full text archive of this journal is available at http://www.emerald-library.com using second order edge elements in 3D Z. Ren Laboratoire de GeÂnie Electrique de Paris, UniversiteÂs
More informationarxiv: v3 [math.na] 8 Sep 2015
A Recovery-Based A Posteriori Error Estimator for H(curl) Interface Problems arxiv:504.00898v3 [math.na] 8 Sep 205 Zhiqiang Cai Shuhao Cao Abstract This paper introduces a new recovery-based a posteriori
More informationA posteriori error estimates in FEEC for the de Rham complex
A posteriori error estimates in FEEC for the de Rham complex Alan Demlow Texas A&M University joint work with Anil Hirani University of Illinois Urbana-Champaign Partially supported by NSF DMS-1016094
More 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 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 informationDISCRETE MAXIMUM PRINCIPLES in THE FINITE ELEMENT SIMULATIONS
DISCRETE MAXIMUM PRINCIPLES in THE FINITE ELEMENT SIMULATIONS Sergey Korotov BCAM Basque Center for Applied Mathematics http://www.bcamath.org 1 The presentation is based on my collaboration with several
More 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 informationCompatible Discretization Schemes
Compatible Discretization Schemes Marc Gerritsma 1 1 Collaborators Artur Palha, Guido Oud, Jasper Kreeft & Mick Bouman Marc Gerritsma Burgers Course CFD I, 25-29 January 2010 1 / 81 What if... Why mimetic
More informationCE-570 Advanced Structural Mechanics - Arun Prakash
Ch1-Intro Page 1 CE-570 Advanced Structural Mechanics - Arun Prakash The BIG Picture What is Mechanics? Mechanics is study of how things work: how anything works, how the world works! People ask: "Do you
More informationA Framework for Analyzing and Constructing Hierarchical-Type A Posteriori Error Estimators
A Framework for Analyzing and Constructing Hierarchical-Type A Posteriori Error Estimators Jeff Ovall University of Kentucky Mathematics www.math.uky.edu/ jovall jovall@ms.uky.edu Kentucky Applied and
More information