A posteriori error estimates applied to flow in a channel with corners
|
|
- Vernon Marvin Porter
- 5 years ago
- Views:
Transcription
1 Mathematics and Computers in Simulation 61 (2003) 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, Faculty of Mechanical Engineering, Czech University of Technology, Karlovo nám. 13, CZ Praha 2, Czech Republic b Institute of Thermomechanics AS CR, Dolejškova 5, CZ Prague 8, Czech Republic Abstract The paper consists of three parts. In the first part, we investigate a posteriori error estimates for the Stokes and Navier Stokes equations on two-dimensional polygonal domains. Special attention is paid to the sources of the constants in the estimates, as these play a crucial role in practical applications to adaptive refinements, as we also show. In the second part, we deal with the problem of determining accurately the constants that appear in the estimates. We present a technique for calculating the constant with high accuracy. In the third part, we apply the a posteriori error estimates with the constants found numerically to the technique of adaptive mesh refinement we solve an incompressible flow problem in a domain with corners that cause singularities in the solution IMACS. Published by Elsevier Science B.V. All rights reserved. MSC: 65M50; 76D05 Keywords: A posteriori error estimate; Stokes problem; Taylor Hood element; Singularities; Adaptive mesh refinement 1. Introduction At present various a posteriori error estimates for the Stokes problem are available. Let us mention, e.g. Ainsworth and Oden [1], Carstensen and Verfürth [7], Johnson et al. [10], Verfürth [13], and other references in [5]. The aim of the first part of the paper is to focus on the aspect of the constant that appears in the estimate because it plays a significant role in the adaptive mesh refinement, cf. also [12]. That is why we derive our own a posteriori estimate and trace carefully the role of different constants and their sources. In [5], we derived an a posteriori error estimate for the Stokes problem in a two-dimensional polygonal domain. In [6], we used similar technique to derive a posteriori error estimates also for three-dimensional domains. Corresponding author. addresses: burda@fsik.cvut.cz (P. Burda), novotny@bivoj.it.cas.cz (J. Novotný), sousedik@student.fsid.cvut.cz (B. Sousedík) /02/$ see front matter 2002 IMACS. Published by Elsevier Science B.V. All rights reserved. PII: S (02)
2 376 P. Burda et al. / Mathematics and Computers in Simulation 61 (2003) The outline of the paper is as follows. In Section 2, we recall the Stokes problem and its finite element approximation with Taylor Hood elements. The theorem on a posteriori error estimate for the Stokes problem is in Section 3. In Sections 4 6, we describe implementation and testing of a posteriori error estimates of the discretization error for Navier Stokes equations. These a posteriori estimates are computed using approximate numerical solution on an initial finite element mesh. On each element they give us information about the discretization error reached on it. This enables us to consider the quality of the mesh, and also to refine the elements of the mesh, where the discretization error is too high and then compute new solution on that new mesh. This way we can continue, until the prescribed accuracy is reached. Numerical results are demonstrated in Section 7 on a computation of fluid flow in domain with corner singularity. 2. The Stokes problem and its finite element approximation Let us consider the Stokes problem on a bounded Lipschitzian domain Ω R d, with either d = 2or d = 3: given f L 2 (Ω), find {v,p} H 1 (Ω) d L 2 0 (Ω) such that, in the weak sense, ν v + p = f in Ω, div v = 0 in Ω, v = 0 on Ω, (1) where v is the velocity vector, p the pressure, and ν>0 is the viscosity. L 2 0 (Ω) is the space of L2 functions having mean value zero. Let us denote (, ) 0 the scalar product in L 2, and let V = H 1 0 (Ω)d L 2 0 (Ω). Problem (1) is known to have a unique solution, cf. [9], and consists in: find {v,p} V such that ν( v, v ) 0 (p, div v ) 0 + (p, div v) 0 = (f, v ) 0, {v,p } V. (2) For the finite element approximation, we suppose Ω to be a polygon in R 2 or a polyhedron in R 3, for simplicity. Let {T h } h 0 be a regular (cf. [9]) family of triangulations of Ω. InR 3 we decompose Ω into tetrahedra. Let P m be polynomials of degrees m>0. We let X h,m h be the finite element spaces of Taylor Hood elements (cf. [3]), i.e. X h ={v H0 1(Ω)d C(Ω),v/K P 2 (K) d,k T h }, M h ={p L 2 0 (Ω), p/k P 1 (K), K T h }, (3) where d = 2ord = 3 as in (1), cf. Fig. 1 for the case d = 2. Velocity values are given in corner nodes and midside nodes of the quadrilateral or the triangle, and pressure values only in corner nodes, in order to satisfy Babuška Brezzi condition [3]. Then velocities and pressure are approximated as continuous functions of spatial variables. The finite element approximation of the Stokes problem consists in finding {v h,p h } X h M h such that, {v h,ph } Xh M h, ν( v h, v h ) 0 (p h, div v h ) 0 + (p h, div vh ) 0 = (f, v h ) 0. (4)
3 P. Burda et al. / Mathematics and Computers in Simulation 61 (2003) Fig. 1. Referential elements. 3. A posteriori error estimate for the Stokes problem Applying the ideas of Eriksson et al. [8], and those of Babuška and Rheinboldt [2] we proved in [5] the estimate for the Taylor Hood elements on the Stokes problem. We define the residual components on the elements K T h, by the relations R 1 (v h,p h ) = f + ν v h p h, R 2 (v h,p h ) = div v h. (5) The error components are defined on Ω by e v = v v h, e p = p p h, where {v,p} is the exact solution defined in Eq. (2), {v h,p h } is the approximate solution, by Eq. (4). The V norm of {e v,e p } is {e v,e p } 2 V = (e v, e v ) 1 + (e p,e p ) 0. Using the Poincaré Friedrichs inequality, Galerkin orthogonality, the Schwarz inequality, the interpolation properties of X h, M h, and the estimate of the solution of the dual problem, we get the following theorem (proof in [5,6]). Theorem 1. Let Ω be a polygon in R 2 or a polyhedron in R 3, with Lipschitz continuous boundary. Let T h be a regular family of triangulations of Ω. Let {v h,p h } be the Hood Taylor approximation of the solution {v,p} of the Stokes problem. Then the error {e v,e p } satisfies the following a posteriori estimate: e v 1 + e p 0 ( 2C P C I C R h K R 1 (v h,p h ) 0,K + R 2 (v h,p h ) 0,K + h 1/2 1 K 2 K T h where C P,C I,C R are positive constants. l K [ν vh n ] l 0,l ), (6)
4 378 P. Burda et al. / Mathematics and Computers in Simulation 61 (2003) Remarks The constants C P,C I, and C R in Theorem 1 come, in turn, from the Poincaré inequality, the interpolation properties of X h,m h, and the regularity of the dual problem, respectively. Our result in Theorem 1 is in agreement with that of Verfürth [13], though the technique of the proof is different, and we do not require any regularity. Let us note that an algorithm that generates local refinements of tetrahedral meshes has been proposed by Křížek and Strouboulis [11]. 4. Steady Navier Stokes equations in 2D Let Ω R 2 be a bounded domain with boundary Ω, f = (f 1,f 2 ) the volumetric loads, and ν the kinematic viscosity coefficient. The steady Navier Stokes problem for the incompressible fluid consists in finding the velocity v = (v 1,v 2 ), and pressure p defined in Ω and satisfying (v )v ν v + p = f, div v = 0 (7) in Ω together with boundary conditions on disjoint parts of the boundary Γ in, Γ wall and Γ out (meaning, in turn, the inlet, the wall, and the outlet parts), v = g on Γ in Γ wall, ν v n pn = 0 on Γ out ( do nothing boundary condition). (8) For the finite element discretization we use again Taylor Hood elements P 2/P A posteriori estimates for the Navier Stokes equations in 2D Suppose, that exact solution of the problem is denoted by (v 1,v 2,p) (v x,v y,p)and the approximate finite element solution by (v h 1,vh 2,p h). The exact solution differs from the approximate solution in the error (e v1,e v2,e p ) = (v 1 v h 2,v 2 v h 2,p p h). For the solution (v 1,v 2,p), we denote (v 1,v 2,p) 2 V = (v 1,v 2 ) 2 1,Ω + p 2 0,Ω ( ( ) 2 = v1 2 + v2 2 + v1 + x Ω ( ) 2 v1 + y ( ) 2 v2 + x The estimate in Theorem 1 can be generalized to the Navier Stokes equations: ( ) ) 2 v2 dω + p 2 dω. y Ω (e v1,e v2 ) 2 1,Ω + e p 2 0,Ω E2 (v h 1,vh 2,ph ), (9) where (cf. [13]) [ E 2 (v1 h,vh 2,ph,Ω)= C h 2 (r1 2 + r2 2 ) dω + ] r3 2 dω, (10) l K Ω l l K Ω l
5 P. Burda et al. / Mathematics and Computers in Simulation 61 (2003) where h denotes the diameter of the finite element and r i stand for the residuals ( r 1 = f x v1 h v1 h x + v v h ) ( 1 2 v1 h 2 + ν y x + 2 v h ) 1 ph 2 y 2 x, ( r 2 = f y v1 h v2 h x + v h ) ( vh 2 2 v2 h 2 + ν y x + 2 v h ) 2 ph 2 y 2 y, r 3 = vh 1 x + vh 2 y. Let us note that according to our practical experience we use only the element residuals in (10). In this respect cf. also [7]. Qualitatively, the value of the constant C is not simple to determine, the sources are seen in Theorem 1.It is important, that this constant does not depend on the mesh size and so can be determined experimentally for general situation, cf. [12]. By computing the estimates (9), we obtain absolute numbers that will depend on given quantities in different problems. We are mainly interested in the error related to the computed solution, i.e. relative error. This is given by the ratio of absolute norm of the solution error, related to unit area of the element Ω l, (1/ Ω l )E 2 (v1 h,vh 2,ph,Ω l ), and the solution norm on the whole domain Ω, related to unit area (1/ Ω l ) (v1 h,vh 2,ph ) 2 V,Ω, i.e. R 2 (v1 h,vh 2,ph,Ω l ) = Ω E2 (v1 h,vh 2,ph,Ω l ) Ω l (v1 h,vh 2,ph ) 2. (11) V,Ω 6. Example and determination of the constant C in the estimates 6.1. Example Steady Stokes equations on one element, Ω = [0; 2] [0; 2]. Let the exact solution and corresponding boundary conditions be given as follows: v 1 = 1 + x 3, v 2 = 3x 2 y, p = 3νx 2 3νy 2, x [0; 2],y [0; 2], Fig. 2. Geometry of the channel.
6 380 P. Burda et al. / Mathematics and Computers in Simulation 61 (2003) and the approximate solution as v h 1 = 3x2 2x + 1, v h 2 = 3x2 y, p h = 0.04xy + bx + cy, where b = ,c = Fig. 3. Initial finite element mesh. Fig. 4. Relative errors on elements of initial mesh. Fig. 5. (a) Finite element mesh after the first refinement; (b) finite element mesh after second refinement; (c) finite element mesh after third refinement.
7 P. Burda et al. / Mathematics and Computers in Simulation 61 (2003) For determining the constant C in estimates (9), (10) we compute error estimate of problem, where exact solution is known and also the norm of the error between exact and approximate solution. The required constant is then given by ratio of these numbers. Then, we have: } E 2 (v1 h,vh 2,ph ) = C {h 2 [(ay+b 6ν) 2 + (ax + c + 6νy) 2 ]dx dy + (3x 2 6x + 2) 2 dx dy Ω Ω = C Fig. 6. (a) Streamlines near the corners after third refinement; (b) velocity v x after third refinement; (c) velocity v y after third refinement; (d) pressure p after third refinement.
8 382 P. Burda et al. / Mathematics and Computers in Simulation 61 (2003) Our computer code gives using three-point Gaussian integration E 2 computed (vh 1,vh 2,ph ) = C For determining the constant C we now have to compute the exact error norm between approximate and exact solution, which is available in this special case. For the error between approximate and exact solution we can write: v 1 v1 h = x3 3x 2 + 2x, v 2 v2 h = 0, p p h = 0.03x y 2 axy bx cy. Substituting, we get (v 1 v h 1,v 2 v h 2,p ph ) analytical = And so we get C = Numerical results and application of estimates to the construction of adaptive meshes Consider two-dimensional flow of viscous, incompressible fluid described by Navier Stokes equations in domain with corner singularity, cf. Fig. 2. Due to symmetry, we solve the problem only on half of the channel, cf. Fig. 3. On the inflow we consider parabolic velocity profile, at the outflow do nothing boundary condition. On the upper wall, no-slip condition and on the lower wall, condition of symmetry (i.e. only y component of velocity equals zero). We consider the following parameters: ν = m 2 /s, v in = 1 m/s. The initial mesh is in Figs. 3 and 4. Elements, where the relative error exceeds 3% are refined, and new solution together with new error estimates is computed. In Fig. 5a c, we show three steps of adaptive refinement of the mesh. In Fig. 6a d, we give the results reached on the third level of the mesh refinement. The singularity in v y component of velocity and in pressure is well seen in Fig. 6c and d, respectively. 8. Conclusions In the estimate (6) only the constant C P can be evaluated directly. The others are not known. But in the application to the adaptive mesh refinement we need the value of the constant. We developed a technique for calculating the constant with high accuracy. Of course, this approach needs some improvement, and it will be a subject of future research. Let us note that another way of local mesh refinement near the singularity has been suggested in [4]. This is based on the asymptotic behavior of the solution near the singularity, not on a posteriori error estimates.
9 Acknowledgements P. Burda et al. / Mathematics and Computers in Simulation 61 (2003) This research has been supported partly by the State Research Project No. J04/98/ and partly by the State Research Project No. J04/98/ References [1] M. Ainsworth, J.T. Oden, A posteriori error estimators for the Stokes and Oseen problems, SIAM J. Numer. Anal. 34 (1997) [2] I. Babuška, W.C. Rheinboldt, A posteriori error estimates for the finite element method, Int. J. Numer. Meth. Eng. 12 (1978) [3] F. Brezzi, M. Fortin, Mixed and Hybrid Finite Element Methods, Springer, Berlin, [4] P. Burda, On the F.E.M. for the Navier Stokes equations in domains with corner singularities, in: M. Křížek, P. Neittaanmäki, R. Stenberg (Eds.), Finite Element Methods, Supeconvergence, Post-Processing and A Posteriori Estimates, Marcel Dekker, New York, 1998, pp [5] P. Burda, An a posteriori error estimate for the Stokes problem in a polygonal domain using Hood Taylor elements, in: P. Neittaanmäki, T. Tiihonen, P. Tarvainen (Eds.), Proceedings of the 3rd European Conference on Numerical Mathematics and Advanced Applications, ENUMATH 99, Jyväskylä, Finland, July 1999, World Scientific, Singapore, 2000, pp [6] P. Burda, A posteriori error estimates for the Stokes flow in 2D and 3D domains, in: P. Neittaanmäki, M. Křížek (Eds.), Finite Element Methods, 3D Problems, GAKUTO Int. Ser., Math. Sci. Appl. 15 (2001) [7] C. Carstensen, R. Verfürth, Edge residuals dominate a posteriori error estimates for low order finite element methods, SIAM J. Numer. Anal. 36 (5) (1999) [8] K. Eriksson, D. Estep, P. Hansbo, C. Johnson, Introduction to adaptive methods for differential equations, Acta Numerica, CUP, 1995, pp [9] V. Girault, P.G. Raviart, Finite Element Method for Navier Stokes Equations, Springer, Berlin, [10] C. Johnson, R. Rannacher, M. Boman, Numerics and hydrodynamic stability: towards error control in computational fluid dynamics, SIAM J. Numer. Anal. 32 (1995) [11] M. Křížek, T. Stroboulis, How to generate local refinements of unstructured tetrahedral meshes satisfying a regularity ball condition, Numer. Meth. Partial Differential Equations 13 (1997) [12] B. Sousedík, P. Burda, A. Szölös, J. Novotný, Application of a posteriori error estimates to fluid flow in a domain with singularity, in: K. Kozel, J. Příhoda, M. Feistauer (Eds.), Proceedings of the 4th Seminar on Euler and Navier Stokes Equations, IT AS CR Praha, 2001, pp [13] R. Verfürth, A Review of A Posteriori Error Estimation and Adaptive Mesh-Refinement Techniques, Wiley and Teubner, Chichester, 1996.
Finite Element Modeling of Incompressible Fluid Flows
PANM 2008 Programs and Algorithms of Numerical Mathematics Dolní Maxov, June 1 6, 2008 Finite Element Modeling of Incompressible Fluid Flows Pavel Burda, Jaroslav Novotný, Jakub Šístek, Alexandr Damašek
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 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 informationChapter 6 A posteriori error estimates for finite element approximations 6.1 Introduction
Chapter 6 A posteriori error estimates for finite element approximations 6.1 Introduction The a posteriori error estimation of finite element approximations of elliptic boundary value problems has reached
More informationNUMERICAL SIMULATION OF INTERACTION BETWEEN INCOMPRESSIBLE FLOW AND AN ELASTIC WALL
Proceedings of ALGORITMY 212 pp. 29 218 NUMERICAL SIMULATION OF INTERACTION BETWEEN INCOMPRESSIBLE FLOW AND AN ELASTIC WALL MARTIN HADRAVA, MILOSLAV FEISTAUER, AND PETR SVÁČEK Abstract. The present paper
More informationSTEADY AND UNSTEADY 2D NUMERICAL SOLUTION OF GENERALIZED NEWTONIAN FLUIDS FLOW. Radka Keslerová, Karel Kozel
Conference Applications of Mathematics 1 in honor of the th birthday of Michal Křížek. Institute of Mathematics AS CR, Prague 1 STEADY AND UNSTEADY D NUMERICAL SOLUTION OF GENERALIZED NEWTONIAN FLUIDS
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 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 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 informationOn angle conditions in the finite element method. Institute of Mathematics, Academy of Sciences Prague, Czech Republic
On angle conditions in the finite element method Michal Křížek Institute of Mathematics, Academy of Sciences Prague, Czech Republic Joint work with Jan Brandts (University of Amsterdam), Antti Hannukainen
More informationUne méthode de pénalisation par face pour l approximation des équations de Navier-Stokes à nombre de Reynolds élevé
Une méthode de pénalisation par face pour l approximation des équations de Navier-Stokes à nombre de Reynolds élevé CMCS/IACS Ecole Polytechnique Federale de Lausanne Erik.Burman@epfl.ch Méthodes Numériques
More informationA Two-Grid Stabilization Method for Solving the Steady-State Navier-Stokes Equations
A Two-Grid Stabilization Method for Solving the Steady-State Navier-Stokes Equations Songul Kaya and Béatrice Rivière Abstract We formulate a subgrid eddy viscosity method for solving the steady-state
More informationA posteriori error estimates for Maxwell Equations
www.oeaw.ac.at A posteriori error estimates for Maxwell Equations J. Schöberl RICAM-Report 2005-10 www.ricam.oeaw.ac.at A POSTERIORI ERROR ESTIMATES FOR MAXWELL EQUATIONS JOACHIM SCHÖBERL Abstract. Maxwell
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 informationNodal O(h 4 )-superconvergence of piecewise trilinear FE approximations
Preprint, Institute of Mathematics, AS CR, Prague. 2007-12-12 INSTITTE of MATHEMATICS Academy of Sciences Czech Republic Nodal O(h 4 )-superconvergence of piecewise trilinear FE approximations Antti Hannukainen
More informationA mixed finite element approximation of the Stokes equations with the boundary condition of type (D+N)
wwwijmercom Vol2, Issue1, Jan-Feb 2012 pp-464-472 ISSN: 2249-6645 A mixed finite element approximation of the Stokes equations with the boundary condition of type (D+N) Jaouad El-Mekkaoui 1, Abdeslam Elakkad
More informationAn interpolation operator for H 1 functions on general quadrilateral and hexahedral meshes with hanging nodes
An interpolation operator for H 1 functions on general quadrilateral and hexahedral meshes with hanging nodes Vincent Heuveline Friedhelm Schieweck Abstract We propose a Scott-Zhang type interpolation
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 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 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 informationOn the equivalence of regularity criteria for triangular and tetrahedral finite element partitions
Computers and Mathematics with Applications 55 (2008) 2227 2233 www.elsevier.com/locate/camwa On the equivalence of regularity criteria for triangular and tetrahedral finite element partitions Jan Brandts
More informationA GAGLIARDO NIRENBERG INEQUALITY, WITH APPLICATION TO DUALITY-BASED A POSTERIORI ESTIMATION IN THE L 1 NORM
27 Kragujevac J. Math. 30 (2007 27 43. A GAGLIARDO NIRENBERG INEQUALITY, WITH APPLICATION TO DUALITY-BASED A POSTERIORI ESTIMATION IN THE L 1 NORM Endre Süli Oxford University Computing Laboratory, Wolfson
More informationELLIPTIC RECONSTRUCTION AND A POSTERIORI ERROR ESTIMATES FOR PARABOLIC PROBLEMS
ELLIPTIC RECONSTRUCTION AND A POSTERIORI ERROR ESTIMATES FOR PARABOLIC PROBLEMS CHARALAMBOS MAKRIDAKIS AND RICARDO H. NOCHETTO Abstract. It is known that the energy technique for a posteriori error analysis
More informationarxiv: v1 [math.na] 19 Dec 2017
Arnold-Winther Mixed Finite Elements for Stokes Eigenvalue Problems Joscha Gedicke Arbaz Khan arxiv:72.0686v [math.na] 9 Dec 207 Abstract This paper is devoted to study the Arnold-Winther mixed finite
More informationPREPRINT 2010:23. A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO
PREPRINT 2010:23 A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO Department of Mathematical Sciences Division of Mathematics CHALMERS UNIVERSITY OF TECHNOLOGY UNIVERSITY OF GOTHENBURG
More informationA UNIFYING THEORY OF A POSTERIORI FINITE ELEMENT ERROR CONTROL
A UNIFYING THEORY OF A POSTERIORI FINITE ELEMENT ERROR CONTROL C. CARSTENSEN Abstract. Residual-based a posteriori error estimates are derived within a unified setting for lowest-order conforming, nonconforming,
More informationA NEW FAMILY OF STABLE MIXED FINITE ELEMENTS FOR THE 3D STOKES EQUATIONS
MATHEMATICS OF COMPUTATION Volume 74, Number 250, Pages 543 554 S 0025-5718(04)01711-9 Article electronically published on August 31, 2004 A NEW FAMILY OF STABLE MIXED FINITE ELEMENTS FOR THE 3D STOKES
More informationOn discontinuity capturing methods for convection diffusion equations
On discontinuity capturing methods for convection diffusion equations Volker John 1 and Petr Knobloch 2 1 Universität des Saarlandes, Fachbereich 6.1 Mathematik, Postfach 15 11 50, 66041 Saarbrücken, Germany,
More informationComplementarity based a posteriori error estimates and their properties
Complementarity based a posteriori error estimates and their properties Tomáš Vejchodský November 3, 2009 Institute of Mathematics, Czech Academy of Sciences Žitná 25, CZ-115 67 Praha 1, Czech Republic
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 informationNUMERICAL SOLUTION OF CONVECTION DIFFUSION EQUATIONS USING UPWINDING TECHNIQUES SATISFYING THE DISCRETE MAXIMUM PRINCIPLE
Proceedings of the Czech Japanese Seminar in Applied Mathematics 2005 Kuju Training Center, Oita, Japan, September 15-18, 2005 pp. 69 76 NUMERICAL SOLUTION OF CONVECTION DIFFUSION EQUATIONS USING UPWINDING
More informationDiscrete Maximum Principle for a 1D Problem with Piecewise-Constant Coefficients Solved by hp-fem
The University of Texas at El Paso Department of Mathematical Sciences Research Reports Series El Paso, Texas Research Report No. 2006-10 Discrete Maximum Principle for a 1D Problem with Piecewise-Constant
More informationON USING ARTIFICIAL COMPRESSIBILITY METHOD FOR SOLVING TURBULENT FLOWS
Conference Applications of Mathematics 212 in honor of the 6th birthday of Michal Křížek. Institute of Mathematics AS CR, Prague 212 ON USING ARTIFICIAL COMPRESSIBILITY METHOD FOR SOLVING TURBULENT FLOWS
More informationA connection between grad-div stabilized FE solutions and pointwise divergence-free FE solutions on general meshes
A connection between grad-div stabilized FE solutions and pointwise divergence-free FE solutions on general meshes Sarah E. Malick Sponsor: Leo G. Rebholz Abstract We prove, for Stokes, Oseen, and Boussinesq
More informationPSEUDO-COMPRESSIBILITY METHODS FOR THE UNSTEADY INCOMPRESSIBLE NAVIER-STOKES EQUATIONS
PSEUDO-COMPRESSIBILITY METHODS FOR THE UNSTEADY INCOMPRESSIBLE NAVIER-STOKES EQUATIONS Jie Shen Department of Mathematics, Penn State University University Par, PA 1680, USA Abstract. We present in this
More informationA STOKES INTERFACE PROBLEM: STABILITY, FINITE ELEMENT ANALYSIS AND A ROBUST SOLVER
European Congress on Computational Methods in Applied Sciences and Engineering ECCOMAS 2004 P. Neittaanmäki, T. Rossi, K. Majava, and O. Pironneau (eds.) O. Nevanlinna and R. Rannacher (assoc. eds.) Jyväskylä,
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 informationYongdeok Kim and Seki Kim
J. Korean Math. Soc. 39 (00), No. 3, pp. 363 376 STABLE LOW ORDER NONCONFORMING QUADRILATERAL FINITE ELEMENTS FOR THE STOKES PROBLEM Yongdeok Kim and Seki Kim Abstract. Stability result is obtained for
More informationA 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 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 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 informationOn Pressure Stabilization Method and Projection Method for Unsteady Navier-Stokes Equations 1
On Pressure Stabilization Method and Projection Method for Unsteady Navier-Stokes Equations 1 Jie Shen Department of Mathematics, Penn State University University Park, PA 1682 Abstract. We present some
More informationA NOTE ON THE LADYŽENSKAJA-BABUŠKA-BREZZI CONDITION
A NOTE ON THE LADYŽENSKAJA-BABUŠKA-BREZZI CONDITION JOHNNY GUZMÁN, ABNER J. SALGADO, AND FRANCISCO-JAVIER SAYAS Abstract. The analysis of finite-element-like Galerkin discretization techniques for the
More informationPolynomial Preserving Recovery for Quadratic Elements on Anisotropic Meshes
Polynomial Preserving Recovery for Quadratic Elements on Anisotropic Meshes Can Huang, 1 Zhimin Zhang 1, 1 Department of Mathematics, Wayne State University, Detroit, Michigan 480 College of Mathematics
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 informationExact a posteriori error analysis of the least squares nite element method 1
Applied Mathematics and Computation 116 (2000) 297±305 www.elsevier.com/locate/amc Exact a posteriori error analysis of the least squares nite element method 1 Jinn-Liang Liu * Department of Applied Mathematics,
More informationSimulating Drag Crisis for a Sphere Using Skin Friction Boundary Conditions
Simulating Drag Crisis for a Sphere Using Skin Friction Boundary Conditions Johan Hoffman May 14, 2006 Abstract In this paper we use a General Galerkin (G2) method to simulate drag crisis for a sphere,
More informationANALYSIS OF THE FEM AND DGM FOR AN ELLIPTIC PROBLEM WITH A NONLINEAR NEWTON BOUNDARY CONDITION
Proceedings of EQUADIFF 2017 pp. 127 136 ANALYSIS OF THE FEM AND DGM FOR AN ELLIPTIC PROBLEM WITH A NONLINEAR NEWTON BOUNDARY CONDITION MILOSLAV FEISTAUER, ONDŘEJ BARTOŠ, FILIP ROSKOVEC, AND ANNA-MARGARETE
More informationChapter 2. General concepts. 2.1 The Navier-Stokes equations
Chapter 2 General concepts 2.1 The Navier-Stokes equations The Navier-Stokes equations model the fluid mechanics. This set of differential equations describes the motion of a fluid. In the present work
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 informationAveraging technique for FE ± a posteriori error control in elasticity. Part II: k-independent estimates
Comput. Methods Appl. Mech. Engrg. 190 2001) 4663±4675 www.elsevier.com/locate/cma Averaging technique for FE ± a posteriori error control in elasticity. Part II: k-independent estimates Carsten Carstensen,
More informationA Two-grid Method for Coupled Free Flow with Porous Media Flow
A Two-grid Method for Coupled Free Flow with Porous Media Flow Prince Chidyagwai a and Béatrice Rivière a, a Department of Computational and Applied Mathematics, Rice University, 600 Main Street, Houston,
More informationInteraction of Incompressible Fluid and Moving Bodies
WDS'06 Proceedings of Contributed Papers, Part I, 53 58, 2006. ISBN 80-86732-84-3 MATFYZPRESS Interaction of Incompressible Fluid and Moving Bodies M. Růžička Charles University, Faculty of Mathematics
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 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 informationOn a Suitable Weak Solution of the Navier Stokes Equation with the Generalized Impermeability Boundary Conditions
Proceedings of the 3rd IASME/WSEAS Int. Conf. on FLUID DYNAMICS & AERODYNAMICS, Corfu, Greece, August -, 5 pp36-41 On a Suitable Weak Solution of the Navier Stokes Equation with the Generalized Impermeability
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 informationPREPRINT 2008:17. A Stokes model with cavitation for the numerical simulation of hydrodynamic lubrication BERTIL NILSSON PETER HANSBO
PREPRINT 28:17 A Stokes model with cavitation for the numerical simulation of hydrodynamic lubrication BERTIL NILSSON PETER HANSBO Department of Mathematical Sciences Division of Mathematics CHALMERS UNIVERSITY
More informationNon-Newtonian Fluids and Finite Elements
Non-Newtonian Fluids and Finite Elements Janice Giudice Oxford University Computing Laboratory Keble College Talk Outline Motivating Industrial Process Multiple Extrusion of Pastes Governing Equations
More informationPseudo-divergence-free element free Galerkin method for incompressible fluid flow
International Workshop on MeshFree Methods 23 Pseudo-divergence-free element free Galerkin method for incompressible fluid flow Y. Vidal () and A. Huerta (2) Abstract: Incompressible modelling in finite
More informationMATHEMATICAL MODELING AND NUMERICAL SOLUTION OF 3D ATMOSPHERIC BOUNDARY LAYER
, Vol, Pt, Special Issue Proceedings of International Conference RDAMM 585 MATHEMATICAL MODELING AND NUMERICAL SOLUTION OF D ATMOSPHERIC BOUNDARY LAYER L. Beneš, K. Kozel Department of Technical Mathematics,
More informationA NOTE ON CONSTANT-FREE A POSTERIORI ERROR ESTIMATES
A NOTE ON CONSTANT-FREE A POSTERIORI ERROR ESTIMATES R. VERFÜRTH Abstract. In this note we look at constant-free a posteriori error estimates from a different perspective. We show that they can be interpreted
More informationNumerical simulation of steady and unsteady flow for generalized Newtonian fluids
Journal of Physics: Conference Series PAPER OPEN ACCESS Numerical simulation of steady and unsteady flow for generalized Newtonian fluids To cite this article: Radka Keslerová et al 2016 J. Phys.: Conf.
More informationHIGHER-ORDER LINEARLY IMPLICIT ONE-STEP METHODS FOR THREE-DIMENSIONAL INCOMPRESSIBLE NAVIER-STOKES EQUATIONS
STUDIA UNIV. BABEŞ BOLYAI, MATHEMATICA, Volume LIII, Number 1, March 2008 HIGHER-ORDER LINEARLY IMPLICIT ONE-STEP METHODS FOR THREE-DIMENSIONAL INCOMPRESSIBLE NAVIER-STOKES EQUATIONS IOAN TELEAGA AND JENS
More informationShape Optimization in Problems Governed by Generalised Navier Stokes Equations: Existence Analysis
Shape Optimization in Problems Governed by Generalised Navier Stokes Equations: Existence Analysis Jaroslav Haslinger 1,3,JosefMálek 2,4, Jan Stebel 1,5 1 Department of Numerical Mathematics, Faculty of
More informationInstitut für Mathematik
U n i v e r s i t ä t A u g s b u r g Institut für Mathematik Dietrich Braess, Carsten Carstensen, and Ronald H.W. Hoppe Convergence Analysis of a Conforming Adaptive Finite Element Method for an Obstacle
More informationON SINGULAR PERTURBATION OF THE STOKES PROBLEM
NUMERICAL ANALYSIS AND MATHEMATICAL MODELLING BANACH CENTER PUBLICATIONS, VOLUME 9 INSTITUTE OF MATHEMATICS POLISH ACADEMY OF SCIENCES WARSZAWA 994 ON SINGULAR PERTURBATION OF THE STOKES PROBLEM G. M.
More informationSome remarks on grad-div stabilization of incompressible flow simulations
Some remarks on grad-div stabilization of incompressible flow simulations Gert Lube Institute for Numerical and Applied Mathematics Georg-August-University Göttingen M. Stynes Workshop Numerical Analysis
More informationAdaptive Finite Element Methods Lecture 1: A Posteriori Error Estimation
Adaptive Finite Element Methods Lecture 1: A Posteriori Error Estimation Department of Mathematics and Institute for Physical Science and Technology University of Maryland, USA www.math.umd.edu/ rhn 7th
More informationThe Role of the Pressure in Finite Element Methods for Incompressible Flow Problems
Weierstrass Institute for Applied Analysis and Stochastics The Role of the Pressure in Finite Element Methods for Incompressible Flow Problems Volker John Mohrenstrasse 39 10117 Berlin Germany Tel. +49
More informationA POSTERIORI ERROR ESTIMATES FOR MAXWELL EQUATIONS
A POSTERIORI ERROR ESTIMATES FOR MAXWELL EQUATIONS JOACHIM SCHÖBERL Abstract. Maxwell equations are posed as variational boundary value problems in the function space H(curl) and are discretized by Nédélec
More informationA Posteriori Estimates for Cost Functionals of Optimal Control Problems
A Posteriori Estimates for Cost Functionals of Optimal Control Problems Alexandra Gaevskaya, Ronald H.W. Hoppe,2 and Sergey Repin 3 Institute of Mathematics, Universität Augsburg, D-8659 Augsburg, Germany
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 informationPOLYNOMIAL PRESERVING GRADIENT RECOVERY AND A POSTERIORI ESTIMATE FOR BILINEAR ELEMENT ON IRREGULAR QUADRILATERALS
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume, Number, Pages 24 c 2004 Institute for Scientific Computing and Information POLYNOMIAL PRESERVING GRADIENT RECOVERY AND A POSTERIORI ESTIMATE
More informationNew Discretizations of Turbulent Flow Problems
New Discretizations of Turbulent Flow Problems Carolina Cardoso Manica and Songul Kaya Merdan Abstract A suitable discretization for the Zeroth Order Model in Large Eddy Simulation of turbulent flows is
More informationEnhancing eigenvalue approximation by gradient recovery on adaptive meshes
IMA Journal of Numerical Analysis Advance Access published October 29, 2008 IMA Journal of Numerical Analysis Page 1 of 15 doi:10.1093/imanum/drn050 Enhancing eigenvalue approximation by gradient recovery
More informationAnalysis of the Oseen-Viscoelastic Fluid Flow Problem
Analysis of the Oseen-Viscoelastic Fluid Flow Problem Vincent J. Ervin Hyesuk K. Lee Louis N. Ntasin Department of Mathematical Sciences Clemson University Clemson, South Carolina 9634-0975 December 30,
More informationA high order adaptive finite element method for solving nonlinear hyperbolic conservation laws
A high order adaptive finite element method for solving nonlinear hyperbolic conservation laws Zhengfu Xu, Jinchao Xu and Chi-Wang Shu 0th April 010 Abstract In this note, we apply the h-adaptive streamline
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 informationRESIDUAL BASED ERROR ESTIMATES FOR THE SPACE-TIME DISCONTINUOUS GALERKIN METHOD APPLIED TO NONLINEAR HYPERBOLIC EQUATIONS
Proceedings of ALGORITMY 2016 pp. 113 124 RESIDUAL BASED ERROR ESTIMATES FOR THE SPACE-TIME DISCONTINUOUS GALERKIN METHOD APPLIED TO NONLINEAR HYPERBOLIC EQUATIONS VÍT DOLEJŠÍ AND FILIP ROSKOVEC Abstract.
More informationConstrained H 1 -interpolation on quadrilateral and hexahedral meshes with hanging nodes
Constrained H 1 -interpolation on quadrilateral and hexahedral meshes with hanging nodes Vincent Heuveline Friedhelm Schieweck Abstract We propose a Scott-Zhang type finite element interpolation operator
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 informationarxiv: v1 [math.na] 28 Feb 2008
BDDC by a frontal solver and the stress computation in a hip joint replacement arxiv:0802.4295v1 [math.na] 28 Feb 2008 Jakub Šístek a, Jaroslav Novotný b, Jan Mandel c, Marta Čertíková a, Pavel Burda a
More informationTraces and Duality Lemma
Traces and Duality Lemma Recall the duality lemma with H / ( ) := γ 0 (H ()) defined as the trace space of H () endowed with minimal extension norm; i.e., for w H / ( ) L ( ), w H / ( ) = min{ ŵ H () ŵ
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 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 informationGeneralized Newtonian Fluid Flow through a Porous Medium
Generalized Newtonian Fluid Flow through a Porous Medium V.J. Ervin Hyesuk Lee A.J. Salgado April 24, 204 Abstract We present a model for generalized Newtonian fluid flow through a porous medium. In the
More informationFLOW OVER A GIVEN PROFILE IN A CHANNEL WITH DYNAMICAL EFFECTS
Proceedings of Czech Japanese Seminar in Applied Mathematics 2004 August 4-7, 2004, Czech Technical University in Prague http://geraldine.fjfi.cvut.cz pp. 63 72 FLOW OVER A GIVEN PROFILE IN A CHANNEL WITH
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 informationEXISTENCE AND REGULARITY OF SOLUTIONS FOR STOKES SYSTEMS WITH NON-SMOOTH BOUNDARY DATA IN A POLYHEDRON
Electronic Journal of Differential Equations, Vol. 2017 (2017), No. 147, pp. 1 10. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu EXISTENCE AND REGULARITY OF SOLUTIONS FOR
More informationAxioms of Adaptivity (AoA) in Lecture 3 (sufficient for optimal convergence rates)
Axioms of Adaptivity (AoA) in Lecture 3 (sufficient for optimal convergence rates) Carsten Carstensen Humboldt-Universität zu Berlin 2018 International Graduate Summer School on Frontiers of Applied and
More informationAdaptive Finite Element Methods Lecture Notes Winter Term 2017/18. R. Verfürth. Fakultät für Mathematik, Ruhr-Universität Bochum
Adaptive Finite Element Methods Lecture Notes Winter Term 2017/18 R. Verfürth Fakultät für Mathematik, Ruhr-Universität Bochum Contents Chapter I. Introduction 7 I.1. Motivation 7 I.2. Sobolev and finite
More informationFINITE ELEMENT APPROXIMATION OF ELLIPTIC DIRICHLET OPTIMAL CONTROL PROBLEMS
Numerical Functional Analysis and Optimization, 28(7 8):957 973, 2007 Copyright Taylor & Francis Group, LLC ISSN: 0163-0563 print/1532-2467 online DOI: 10.1080/01630560701493305 FINITE ELEMENT APPROXIMATION
More informationON LIQUID CRYSTAL FLOWS WITH FREE-SLIP BOUNDARY CONDITIONS. Chun Liu and Jie Shen
DISCRETE AND CONTINUOUS Website: http://aimsciences.org DYNAMICAL SYSTEMS Volume 7, Number2, April2001 pp. 307 318 ON LIQUID CRYSTAL FLOWS WITH FREE-SLIP BOUNDARY CONDITIONS Chun Liu and Jie Shen Department
More informationNUMERICAL EXPERIMENTS FOR THE ARNOLD WINTHER MIXED FINITE ELEMENTS FOR THE STOKES PROBLEM
NUMERICAL EXPERIMENTS FOR THE ARNOLD WINTHER MIXED FINITE ELEMENTS FOR THE STOKES PROBLEM CARSTEN CARSTENSEN, JOSCHA GEDICKE, AND EUN-JAE PARK Abstract. The stress-velocity formulation of the stationary
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 informationAn A Posteriori Error Estimate for Discontinuous Galerkin Methods
An A Posteriori Error Estimate for Discontinuous Galerkin Methods Mats G Larson mgl@math.chalmers.se Chalmers Finite Element Center Mats G Larson Chalmers Finite Element Center p.1 Outline We present an
More 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 informationStabilization of finite element method for solving incompressible viscous flows
Czech Technical University in Prague Faculty of Mechanical Engineering Department of Technical Mathematics Diploma thesis Stabilization of finite element method for solving incompressible viscous flows
More informationFULLY COMPUTABLE ERROR BOUNDS FOR EIGENVALUE PROBLEM
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 15, Number 1-2, Pages 260 276 c 2018 Institute for Scientific Computing and Information FULLY COMPUTABLE ERROR BOUNDS FOR EIGENVALUE PROBLEM
More information