ENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED FINITE ELEMENT METHODS
|
|
- Sherilyn Cannon
- 5 years ago
- Views:
Transcription
1 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 for second order elliptic equations It is shown that a reliable and efficient error estimator can be constructed using a postprocessed solution of the method The analysis is performed in two different ways; under a saturation assumption and using a Helmholtz decomposition for vector fields 1 Introduction We consider the mixed finite element approximation of second order elliptic equations with the Poisson problem as a model: u = f in Ω R n, u = 0 The problem is written as the system on Ω σ u = 0, div σ + f = 0, which is approximated with the Mixed method Find σ h, u h S h V h Hdiv:Ω L 2 Ω such that σ h, τ + div τ, u h = 0 τ S h, div σ h, v + f, v = 0 v V h In the method the polynomial used for approximating the flux σ is of higher degree than that used for the displacement u, which is counterintuitive in view of 13 As a consequence, the mixed method has to be carefully designed in order to satisfy the Babuška-Brezzi conditions, cf eg [8] There are two ways of posing these conditions, both yielding the same a priori estimates The more common one is to use the Hdiv : Ω norm for the flux and the L 2 Ω norm for the displacement The other one is to use so called mesh dependent norms [3] which are close to the energy norm of the continuous problem The a posteriori error analysis of mixed methods has been performed in [1], [10] and [5] In [10] the estimate is for the Hdiv : Ω norm This is in a way unsatisfactory since the div part of the norm is trivially computable and also may dominate the error, see Remark 34 below In [5] an estimate for the L 2 -norm 1991 Mathematics Subject Classification 65N30 Key words and phrases mixed finite element methods, a-posteriori error estimates, postprocessing This work has been supported by the European Project HPRN-CT New Materials, Adaptive Systems and their Nonlinearities Modelling, Control and Numerical Simulation 1
2 2 C LOVADINA AND R STENBERG of the flux is derived but it is, however, not optimal The reason for this is that the estimator includes the element residual in the constitutive relation 13 As the polynomial degree of approximation for the displacement is lower than that for the flux, it is clear that this residual is large The purpose of this paper is to point out a simple remedy to this Since the work of Arnold and Brezzi [2] it is known that the mixed finite element solution can be locally postprocessed in order to obtain an improved displacement Later other postprocessing has been proposed [6, 9, 7, 17, 16] On each element the postprocessed displacement is of one degree higher than the flux, which is in accordance with 13 Hence, it is natural to use it in the a posteriori estimate In this paper, we will focus on the postprocessing introduced in [17, 16] In Section 2 we develop an a-priori error analysis by recognizing that the postprocessed output can be viewed as the direct solution of a suitable modified method In Section 3 we introduce our estimator based on the postprocessed solution, and we prove its efficiency and reliability Throughout the paper we will use standard notations for Sobolev norms and seminorms Moreover, we will denote with C and C i i = 1, 2, generic constants independent of the mesh parameter h, which may take different values in different occurrences 2 A-priori estimates and postprocessing In this section we will consider the mixed methods, their postprocessing and error analysis We will also give the stability and error analysis by treating the method and the postprocessing as one method This will be useful for the a posteriori analysis We will use standard notation used in connection with mixed FE methods By C h we denote the finite element regular partitioning and by Γ h the collection of edges or faces of C h The subspaces σ h, u h S h V h Hdiv : Ω L 2 Ω are piecewise polynomial spaces defined on C h In this paper we will consider the following families of elements The results are, however, easily applicable for other families as well RTN elements the triangular elements of Raviart-Thomas [15] and their tetrahedral counterparts of Nedelec [14]; BDM elements the triangular elements of Brezzi-Douglas-Marini [9] and their tetrahedral counterparts of Brezzi-Douglas-Duran-Fortin [7] Accordingly, given an integer k 1, we define: 21 S RT N h = { τ Hdiv:Ω τ K [P k 1 K] n x P k 1 K K C h } 22 S BDM h = { τ Hdiv:Ω τ K [P k K] n K C h } 23 V RT N h = V BDM h = { v L 2 Ω v K P k 1 K K C h }, where P k 1 K denotes the homogeneous polynomials of degree k 1 For quadrilateral and hexahedral meshes there exist a wide choice of different alternatives, cf [8] By defining the following bilinear form 24 Bϕ, w; τ, v = ϕ, τ + div τ, w + div ϕ, v
3 ENERGY NORM A POSTERIORI ESTIMATES FOR MIXED METHODS 3 the mixed method can compactly be defined as: Find σ h, u h S h V h such that 25 Bσ h, u h ; τ, v + f, v = 0 τ, v S h V h For the displacement and the flux we will use the following norms: 26 v 2 1,h = v 2 0,K + h 1 E [[v]] 2 0,E, and 27 τ 2 0,h = τ h E τ n 2 0,E, where n is the unit normal to E Γ h and [[v ] is the jump in v along interior edges/faces and v on edges/faces on Ω By an element by element partial integration we have 28 div τ, v τ 0,h v 1,h τ, v S h V h In the FE subspace the norm for the flux is equivalent to the L 2 -norm: 29 C τ 0,h τ 0 τ 0,h τ S h Hence, it also holds 210 div τ, v C τ 0 v 1,h τ, v S h V h With this choice of norms the Babuška-Brezzi stability condition is the following Lemma 21 There is a positive constant C such that 211 sup τ S h div τ, v τ 0 C v 1,h v V h Proof We first point out that since Vh RT N = Vh BDM the result for BDM is a consequence of that for RTN Therefore, we focus on the RTN family, first recalling that the local degrees of freedom for the flux variable are the following: and S RT N h τ n, z E z P k 1 E, E K, τ, z K z [P k 2 K] n S BDM h Above and in the rest of the paper, we use the notation, K and, E for the L 2 inner product on the element K and on the edge/face E, respectively Hence, given v V h we can define τ S h by τ n, z E = h 1 E [[v]], z E z P k 1 E, E Γ h, τ, z K = v, z K z [P k 2 K] n, K C h Noting that v K [P k 2 K] n, [[v]] E P k 1 E, from we obtain It follows that cf also 26 τ n, [[v]] E = h 1 E [[v ] 2 0,E, τ, v K = v 2 0,K
4 4 C LOVADINA AND R STENBERG 218 div τ, v = Using scaling arguments imply τ, v K + τ n, [[v]] E = v 2 0,K + h 1 E [v]] 2 0,E = v 2 1,h 219 τ 0,h C v 1,h The assertion now follows from 218 and 219 From this stability estimate, the following full stability result holds Lemma 22 There is a positive constant C such that Bϕ, w; τ, v 220 sup C ϕ 0 + w 1,h τ,v S h V h τ 0 + v 1,h ϕ, w S h V h In our analysis we will exploit the interpolation operator R h : Hdiv : Ω [L s Ω] n S h, with s > 2, such that 221 div τ R h τ, v = 0 v V h, which can be constructed by using the degrees of freedom for S h, cf [15, 14, 9, 7] In addition, we will use the equilibrium property 222 div S h V h When denoting by P h : L 2 Ω V h the L 2 -projection, this implies that 223 div τ, u P h u = 0 τ S h The projection and interpolation operators satisfy the following commuting property: 224 div R h = P h div Theorem 23 There is a positive constant C such that 225 σ σ h 0 + P h u u h 1,h C σ R h σ 0 Proof By Lemma 22 there is a pair τ, v S h V h, with τ 0 + v 1,h C, such that 226 σ h R h σ 0 + u h P h u 1,h Bσ h R h σ, u h P h u; τ, v Next, 221, 223 and 224 give 227 Bσ h R h σ, u h P h u; τ, v = σ h R h σ, τ + div τ, u h P h u + div σ h R h σ, v = σ R h σ, τ σ R h σ 0 τ 0 C σ R h σ 0 The assertion then follows from the triangle inequality This gives assuming full regularity: σ σ h 0 + P h u u h 1,h Ch k+1 σ k+1 σ σ h 0 + P h u u h 1,h Ch k σ k for BDM, for RTN
5 ENERGY NORM A POSTERIORI ESTIMATES FOR MIXED METHODS 5 We note that these estimates contain a superconvergence result for P h u u h 1,h This, together with the fact that σ h is a good approximation of u, implies that an improved approximation for the displacement can be constructed by local postprocessing Below we will consider the method introduced in [17, 16] The postprocessed displacement is sought in a FE space V h V h For our choices, the spaces are V BDM h = { v L 2 Ω v K P k+1 K K C h }, Vh RT N = { v L 2 Ω v K P k K K C h } Postprocessing method Find u h V h 232 P h u h = u h and such that 233 u h, v K = σ h, v K v I P h V h K The error analysis of this postprocessing is done in [17, 16] Here we proceed in a slightly different way by considering the method and the postprocessing as one method To this end we define the bilinear form 234 B h ϕ, w ; τ, v =ϕ, τ + div τ, w + div ϕ, v + w ϕ, I P h v K Then we have the following equivalence to the original problem Lemma 24 Let σ h, u h S h V h be the solution to the problem 235 B h σ h, u h; τ, v + P h f, v = 0 τ, v S h V h, and set u h = P h u h V h Then σ h, u h S h V h coincides with the solution of Conversely, let σ h, u h S h V h be the solution of 15 16, and let u h V h be the postprocessed displacement defined by Then σ h, u h S h Vh is the solution to 235 Proof Testing by τ, 0 S h V h in 235 gives 236 σ h, τ + div τ, u h = 0 τ S h The equilibrium property 222 implies 237 div τ, u h = div τ, u h Hence, 15 is satisfied Next, for a generic v Vh set v = P hv V h and observe that V h = P h Vh Testing in 235 with 0, v, and using the fact that P h f, v = f, v, we obtain 238 div σ h, v + f, v = 0 v V h, ie the equation 16 Conversely, let σ h, u h S h V h be the solution of 15 16, and let u h V h be defined by Splitting a generic v V h as v = P h v + I P h v we have
6 6 C LOVADINA AND R STENBERG 239 B h σ h, u h; τ, v = B h σ h, u h; τ, P h v + B h σ h, u h; 0, I P h v = σ h, τ + div τ, u h + div σ h, P h v + + div σ h, I P h v + u h σ h, I P h P h v K u h σ h, I P h I P h v K = σ h, τ + div τ, u h P h f, P h v = P h f, v τ, v S h V h Therefore, σ h, u h S h V h solves 235 Next, we prove the stability In the proof we will use the following norm equivalence Lemma 25 There are positive constants C 1 and C 2, such that 240 w 1,h P h w 1,h + I P h w 1,h C 2 w 1,h and 1/2 241 C 1 w 1,h P h w 1,h + I P h w 2 0,K C2 w 1,h, for every w V h Proof We first prove 240 The estimate w 1,h P h w 1,h + I P h w 1,h follows immediately from the triangle inequality To continue, we notice that 242 P h w 1,h + I P h w 1,h 2 P h w 1,h + w 1,h We now fix an interior edge/face E, and we consider the elements K 1 and K 2 such that E = K 1 K 2 A scaling argument shows that 243 h 1 E [[P hw ]] 2 0,E + 2 P h w 2 0,K i i=1 C h 1 E [[w ]] 2 0,E + If E K is an edge/face lying in Ω, a similar argument gives 2 w 2 0,K i 244 h 1 E P hw 2 0,E + P h w 2 0,K C h 1 E w 2 0,E + w 2 0,K The estimate 245 P h w 1,h + I P h w 1,h C 2 w 1,h easily follows from cf also 26 Hence, 240 is proved To prove 241 we first notice that 245 implies 1/2 246 P h w 1,h + I P h w 2 0,K C2 w 1,h Next, scaling arguments lead to i=1
7 ENERGY NORM A POSTERIORI ESTIMATES FOR MIXED METHODS h 1 E [[w ] 2 0,E w 2 0,K i i=1 C h 1 E [[P hw ]] 2 0,E 2 Ph w 2 0,K i + I P h w 2 0,K i, i=1 for an interior edge/face E, and to 248 h 1 E w 2 0,E + w 2 0,K C h 1 E P hw 2 0,E + P h w 2 0,K + I P h w 2 0,K, for a boundary edge/face E The estimate 1/2 C 1 w 1,h P h w 1,h + I P h w 2 0,K is a consequence of The proof is complete Lemma 26 There is a positive constant constant C such that 249 sup τ,v S h V h B h ϕ, w ; τ, v τ 0 + v 1,h C ϕ 0 + w 1,h ϕ, w S h V h Proof Let ϕ, w S h V h be arbitrary By choosing v = v V h and using the equilibrium condition 222 we then get 250 B h ϕ, w ; τ, v = ϕ, τ + div τ, w + div ϕ, v = ϕ, τ + div τ, P h w + div ϕ, v = Bϕ, P h w ; τ, v, Hence, the stability of Lemma 22 implies that we can choose τ, v such that 251 B h ϕ, w ; τ, v ϕ P h w 2 1,h and 252 τ 0 + v 1,h C 1 ϕ 0 + P h w 1,h
8 8 C LOVADINA AND R STENBERG Next, 210 and Schwarz inequality give 253 B h ϕ, w ; 0, I P h w = div ϕ, I P h w + w ϕ, I P h w K C 2 ϕ 0 I P h w 1,h + w, I P h w K = C 2 ϕ 0 I P h w 1,h + + I P h w 2 0,K P h w, I P h w K C 2 ϕ 0 + P h w 1,h I Ph w 1,h + I P h w 2 0,K We now notice that I P h w is L 2 orthogonal to the piecewise constant functions; therefore, a scaling argument shows that 1/2 254 I P h w 1,h C 3 I P h w 2 0,K For α > 0, we obtain from 253 and B h ϕ, w ; 0, I P h w 1 C2 ϕ 0 + P h w 2 α 1,h 2α 2 I P hw 2 1,h + I P h w 2 0,K C2 ϕ 0 + P h w 2 1,h 1 2α Choosing α > 0 sufficiently small, we get 256 B h ϕ, w ; 0, I P h w C αc2 3 I P h w 2 2 0,K I P h w 2 0,K ϕ 2 0 P h w 2 1,h Combining 251 and 256, with δ > 0 to be chosen, we have 257 B h ϕ, w ; τ, v + δi P h w 1 δc 4 ϕ P h w 2 1,h + δc4 I P h w 2 0,K Next, by 241 we have 258 P h w 2 1,h + δ I P h w 2 0,K C 5 w 2 1,h
9 ENERGY NORM A POSTERIORI ESTIMATES FOR MIXED METHODS 9 From 252 and 240 we have 259 τ 0 + v + δi P h w 1,h τ 0 + v 1,h + δ I P h w 1,h C 1 ϕ 0 + P h w 1,h + δ I Ph w 1,h C 6 ϕ 0 + w 1,h Choosing δ = 1/2C 4, estimate 249 is proved by combining Theorem 27 The following a priori error estimate holds σ σ h 0 + u u h 1,h C σ R h σ 0 + inf v V h u v 1,h Proof From Lemma 26 it follows that there is ϕ, w S h Vh, with ϕ 0 + w 1,h C, such that 260 σh R h σ 0 + u h v 1,h Bh σ h R h σ, u h v ; ϕ, w Next, from the definition of B h and the equations it follows that 261 B h σ, u; ϕ, w + f, w = 0 Hence it holds 262 B h σ h R h σ, u h v ; ϕ, w = B h σ R h σ, u v ; ϕ, w + f P h f, w Writing out the right hand side we have 263 B h σ R h σ, u v ; ϕ, w + f P h f, w = σ R h σ, ϕ + div ϕ, u v + div σ R h σ, w + u v σ R h σ, I P h w K + f P h f, w The commuting property 224 gives 264 div σ R h σ, w = f P h f, w Hence, the third and the last term on the right hand side of 263 cancel The other terms are directly estimated and using σ R h σ, ϕ σ R h σ 0 ϕ 0 C σ R h σ 0, div ϕ, u v C ϕ 0 u v 1,h C u v 1,h u v σ R h σ, I P h w K C u v 1,h + σ R h σ 0 w 1,h C u v 1,h + σ R h σ 0 The assertion then follows by collecting the above estimate and using the triangle inequality For our choices of spaces we obtain the estimates with the assumption of a sufficiently smooth solution
10 10 C LOVADINA AND R STENBERG Corollary 28 There are positive constants C such that 268 σ σ h 0 + u u h 1,h Ch k+1 u k σ σ h 0 + u u h 1,h Ch k u k+1 for BDM, for RTN 3 A-posteriori estimates We define the following local error indicators on the elements 31 η 1,K = u h σ h 0,K, η 2,K = h K f P h f 0,K, and on the edges 32 η E = h 1/2 E [[u h]] 0,E Using these quantities, the global estimator is 33 η = η 2 1,K + η2,k 2 + η 2 E 1/2 The efficiency of the estimator is given by the following lower bounds, which directly follow from 13 using the triangle inequality, and from 32 noting that [[u]] = 0 on each edge E Theorem 31 It holds η 1,K u u h 0,K + σ σ h 0,K, 34 η E = h 1/2 E [[u u h]] 0,E As far as the estimator reliability is concerned, below we will use two different techniques 31 Reliability via a saturation assumption The first technique to prove the upper bound is based on the following saturation assumption We let C h/2 be the mesh obtained from C h by refined each element into 2 n n = 2, 3 elements For clarity all variables in the spaces defined on C h will be equipped with the subscript h whereas h/2 will be used for those defined on C h/2 Accordingly, we let σ h/2, u h/2 S h/2 Vh/2 be the solution to 35 B h/2 σ h/2, u h/2 ; τ h/2, vh/2 + P h/2f, vh/2 = 0 τ h/2, vh/2 S h/2 Vh/2 As already done in [5], we make the following assumption for the solutions of 235 and 35 Saturation assumption There exists a positive constant β < 1 such that 36 σ σ h/2 0 + u u h/2 1,h/2 β σ σ h 0 + u u h 1,h Since it holds 37 u u h 1,h u u h 1,h/2 we also have 38 σ σ h/2 0 + u u h/2 1,h/2 β σ σ h 0 + u u h 1,h/2 Using the triangle inequality we then get 39 σ σ h 0 + u u h 1,h/2 1 σh/2 σ h 0 + u h/2 1 β u h 1,h/2
11 ENERGY NORM A POSTERIORI ESTIMATES FOR MIXED METHODS 11 By again using 37 we obtain 310 σ σ h 0 + u u h 1,h 1 σh/2 σ h 0 + u h/2 1 β u h 1,h/2 We now prove the following result Theorem 32 Suppose that the saturation assumption 36 holds Then there exists a positive constant C such that 311 σ σ h 0 + u u h 1,h Cη Proof By 310 it is sufficient to prove the following bound 312 σ h/2 σ h 0 + u h/2 u h 1,h/2 Cη By Lemma 26 applied to the finer mesh C h/2, there is τ h/2, v h/2 S h/2 V h/2, with τ h/2 0 + v h/2 1,h/2 C, such that σh σ h/2 0 + u h u h/2 1,h/2 313 Using the fact that B h/2 σ h σ h/2, u h u h/2 ; τ h/2, v h/2 314 σ h/2, τ h/2 + div τ h/2, u h/2 = 0 we have 315 B h/2 σ h σ h/2, u h u h/2 ; τ h/2, v h/2 = σ h σ h/2, τ h/2 + div τ h/2, u h u h/2 + div σ h σ h/2, v h/2 + /2 u h u h/2 σ h σ h/2, I P h/2 v h/2 K = σ h, τ h/2 + div τ h/2, u h + div σ h σ h/2, v h/2 + /2 u h σ h, I P h/2 v h/2 K, We now notice that it holds cf C τ h/2 0,h τ h/2 0 τ h/2 0,h τ h/2 S h/2 Therefore, using 316 and 31 33, we obtain 317 σ h, τ h/2 + div τ h/2, u h = σ h u h, τ h/2 K + τ h/2 n, [[u h ] E σ h u h 0,K τ h/2 0,K + τ h/2 n 0,E [[u h ] 0,E η τ h/2 0,h ηc τ h/2 0 Cη Similarly for the last term in 315 we get using u h σ h, I P h/2 vh/2 K Cη I P h/2 vh/2 1,h/2 /2 Cη v h/2 1,h/2 Cη
12 12 C LOVADINA AND R STENBERG When estimating the term div σ h σ h/2, vh/2 in 315 we recall that div σ h = P h f and div σ h/2 = P h/2 f, and that P h, P h/2 are L 2 -projection operators Therefore, we have div σ h σ h/2, v h/2 = P h/2f P h f, v h/2 319 = P h/2f f, v h/2 + f P hf, v h/2 = P h/2 f f, v h/2 P h/2v h/2 + f P hf, v h/2 P hv h/2 Next, we use the following interpolation estimates, which are easily proved by standard scaling arguments cf [5, Lemma 31]: where v h/2 P hv h/2 0,K Ch K v h/2 1,h/2,K, K C h, v h/2 2 1,h/2,K = K i v h/2 2 0,K i + E i h 1 E i [[v h/2 ] 2 0,E i Here K i K are the elements of C h/2 and E i are the edges of Γ h/2 lying in the interior of K This gives f P h f, vh/2 P hvh/2 C h 2 K f P h f 2 1/2 v 320 0,K h/2 1,h/2 C h 2 K f P h f 2 1/2 0,K Cη We also have 321 P h/2 f f, v h/2 P h/2v h/2 f P h/2 f 0,K vh/2 P h/2vh/2 0,K /2 C h K f P h/2 f 0,K vh/2 0,K /2 C h 2 K f P h/2 f 2 1/2 v 0,K h/2 1,h/2 /2 C h 2 K f P h/2 f 2 1/2 0,K /2 C h 2 K f P h/2 f 2 1/2 0,K Since, by the properties of L 2 -projection operators, it holds from 321 we obtain f P h/2 f 0,K f P h f 0,K K C h, 322 P h/2 f f, vh/2 P h/2vh/2 C h 2 K f P h f 2 1/2 0,K Cη By collecting the estimates and 322, from 315 we get 323 B h/2 σ h σ h/2, u h u h/2 ; τ h/2, vh/2 Cη The assertion now follows from 313
13 ENERGY NORM A POSTERIORI ESTIMATES FOR MIXED METHODS 13 We have presented the above proof since this is rather general and can be used for other problems as well In [13] we use it for a plate bending method 32 Reliability via a Helmholtz decomposition Now, let us give another proof of the estimator reliability, not relying on the saturation assumption Theorem 33 Suppose that Ω R 2 is a simply connected domain Then there exists a positive constant C such that 324 σ σ h 0 + u u h 1,h Cη Proof We use the techniques of [11] and [10] We first notice that σ σ h, ϕ 325 σ σ h 0 = sup ϕ L 2 Ω ϕ 0 For a generic ϕ L 2 Ω, we consider the L 2 -orthogonal Helmholtz decomposition see, eg [12]: 326 ϕ = ψ + curl q, ψ H 1 0 Ω, q H 1 Ω/R, with 327 ϕ 0 = ψ curl q 2 0 1/2 Therefore, from we see that it holds σ σ h, ψ σ σ h, curl q 328 σ σ h 0 sup + sup ψ H0 1Ω ψ 1 q H 1 Ω/R q 1 Given ψ H 1 0 Ω, from 14 and 16 it follows that 329 div σ σh, P h ψ = 0 Hence, we have 330 σ σ h, ψ = div σ σ h, ψ As a consequence, we get cf 31 = div σ σ h, ψ P h ψ 1/2 ψ 1 C h 2 K div σ σ h 2 0,K 1/2 ψ 1 C h 2 K f P h f 2 0,K 331 σ σ h, ψ 1/2 sup C h 2 ψ H0 1Ω ψ K f P h f 2 0,K = C 1 η 2 2,K 1/2
14 14 C LOVADINA AND R STENBERG To continue, let I h q be the Clément interpolant of q in the space of continuous piecewise linear functions see [4], for instance satisfying 1/2 332 q I h q 1 + h 1 E q I hq 0,E 2 C q 1 Noting that curl I h q S h, and div curl I h q = 0, from 13 and 15 we get 333 σ σ h, curl I h q = 0 Therefore, using 332, one has 334 σ σ h, curl q = σ σ h, curlq I h q = u σ h, curlq I h q = σ h, curlq I h q = σh u h, curlq I h q + K u h, curlq I h q K 1/2 q 1 C σ h u h 2 0,K + u h, curlq I h q K Furthermore, an integration by parts and standard arguments and 332 give u h, curlq I h q K = u h t, q I h q K = [[ u h t]], q I h q E 335 1/2 1/2 h E [ u h t]] 2 0,E h 1 E q I hq 2 0,E 1/2 q 1 C h 1 E [[u h ] 0,E 2 From 334 and 335 we obtain see 31 and σ σ h, curl q C q 1 sup q H 1 Ω/R Using 331 and 336 we deduce σ h u h 2 0,K + 1/2 = C η1,k 2 + ηe σ σ h 0 C η1,k 2 + η2,k 2 + η 2 E We now estimate the term u u h 1,h We first recall that h 1 E [[u h]] 2 0,E 1/2 338 u u h 1,h = u u h 2 0,K + h 1 E [u u h ] 2 0,E 1/2 1/2
15 ENERGY NORM A POSTERIORI ESTIMATES FOR MIXED METHODS 15 and we notice that cf 32 1/2 339 h 1 E [u 1/2 u h ] 2 0,E = h 1 E [[u h]] 2 0,E = η 2 E 1/2 We have u u h 2 0,K = u u h, u u h K = σ u h, u u h K 340 by which we obtain = σ σ h, u u h + σ K h u h, u u h K σ σ h 0,K + σ h u h 0,K u u h 0,K, 341 u u h 0,K σ σ h 0,K + σ h u h 0,K Hence we infer 1/2 342 u u h 2 0,K σ σh 0 + σ h u h 2 0,K Using 337 and recalling 31, from 342 we get 1/2 343 u u h 2 0,K C η1,k 2 + η2,k 2 + Therefore, joining 339 and 343 we obtain 344 u u h 1,h C η1,k 2 + η2,k 2 + From 337 and 344 we finally deduce see 33 η 2 E 345 σ σ h 0 + u u h 1,h C η1,k 2 + η2,k 2 + We end the paper by the following 1/2 η 2 E η 2 E 1/2 1/2 1/2 = Cη Remark 34 On the estimate in the Hdiv:Ω-norm In the paper we have repeatedly used the fact that by the equilibrium property 222 we have div σ σ h = P h f f and hence div σ σ h 0 = f P h f 0 is a quantity that is directly computable from the data to the problem For the BDM spaces it furthermore holds that for a general loading and a smooth solution it holds f P h f 0 = Oh k, whereas σ σ h 0 = Oh k+1, and hence this trivial component in the Hdiv:Ω norm can dominate the whole estimate References [1] A Alonso, Error estimators for a mixed method, Numer Math, , pp [2] D N Arnold and F Brezzi, Mixed and nonconforming finite element methods: implementation, postprocessing and error estimates, RAIRO Modél Math Anal Numér, , pp 7 32 [3] I Babuška, J Osborn, and J Pitkäranta, Analysis of mixed methods using mesh dependent norms, Math Comp, , pp
16 16 C LOVADINA AND R STENBERG [4] D Braess, Finite Elements: Theory, Fast Solvers, and Applications in Solid Mechanics, Cambridge University Press, 2001 [5] D Braess and R Verfürth, A posteriori error estimators for the Raviart-Thomas element, SIAM J Numer Anal, , pp [6] J H Bramble and J Xu, A local post-processing technique for improving the accuracy in mixed finite-element approximations, SIAM J Numer Anal, , pp [7] F Brezzi, J Douglas, Jr, R Durán, and M Fortin, Mixed finite elements for second order elliptic problems in three variables, Numer Math, , pp [8] F Brezzi and M Fortin, Mixed and Hybrid Finite Element Methods, Springer-Verlag, 1991 [9] F Brezzi, J Douglas Jr, and L Marini, Two families of mixed finite elements for second order elliptic problems, Numer Math, , pp [10] C Carstensen, A posteriori error estimate for the mixed finite element method, Math Comp, , pp [11] E Dari, R Durán, C Padra, and V Vampa, A posteriori error estimators for nonconforming finite element methods, RAIRO Modél Math Anal Numér, , pp [12] V Girault and P Raviart, Finite Element Methods for Navier-Stokes Equations Theory and Algorithms, Springer-Verlag, 1986 [13] C Lovadina and R Stenberg, A posteriori error analysis of the linked interpolation technique for plate bending problems, SIAM J Numer Anal, Accepted for publication [14] J-C Nédélec, A new family of mixed finite elements in R 3, Numer Math, , pp [15] P Raviart and J Thomas, A mixed finite element method for second order elliptic problems, in Mathematical Aspects of the Finite Element Method Lecture Notes in Math 606, Springer- Verlag, 1977, pp [16] R Stenberg, Some new families of finite elements for the Stokes equations, Numer Math, , pp [17] R Stenberg, Postprocessing schemes for some mixed finite elements, RAIRO Modél Math Anal Numér, , pp Dipartimento di Matematica, Università di Pavia and IMATI-CNR, VIa Ferrata 1, Pavia 27100, Italy address: carlolovadina@unipvit Institute of Mathematics, Helsinki University of Technology, PO Box 1100, TKK, Finland address: rolfstenberg@tkkfi
ENERGY 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 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 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 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 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 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 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 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 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 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 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 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 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 informationA-priori and a-posteriori error estimates for a family of Reissner-Mindlin plate elements
A-priori and a-posteriori error estimates for a family of Reissner-Mindlin plate elements Joint work with Lourenco Beirão da Veiga (Milano), Claudia Chinosi (Alessandria) and Carlo Lovadina (Pavia) Previous
More informationNONSTANDARD NONCONFORMING APPROXIMATION OF THE STOKES PROBLEM, I: PERIODIC BOUNDARY CONDITIONS
NONSTANDARD NONCONFORMING APPROXIMATION OF THE STOKES PROBLEM, I: PERIODIC BOUNDARY CONDITIONS J.-L. GUERMOND 1, Abstract. This paper analyzes a nonstandard form of the Stokes problem where the mass conservation
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 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 informationLecture Note III: Least-Squares Method
Lecture Note III: Least-Squares Method Zhiqiang Cai October 4, 004 In this chapter, we shall present least-squares methods for second-order scalar partial differential equations, elastic equations of solids,
More informationarxiv: v1 [math.na] 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 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 informationError estimates for the Raviart-Thomas interpolation under the maximum angle condition
Error estimates for the Raviart-Thomas interpolation under the maximum angle condition Ricardo G. Durán and Ariel L. Lombardi Abstract. The classical error analysis for the Raviart-Thomas interpolation
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 MULTIGRID METHOD FOR THE PSEUDOSTRESS FORMULATION OF STOKES PROBLEMS
SIAM J. SCI. COMPUT. Vol. 29, No. 5, pp. 2078 2095 c 2007 Society for Industrial and Applied Mathematics A MULTIGRID METHOD FOR THE PSEUDOSTRESS FORMULATION OF STOKES PROBLEMS ZHIQIANG CAI AND YANQIU WANG
More informationANALYSIS OF H(DIV)-CONFORMING FINITE ELEMENTS FOR THE BRINKMAN PROBLEM
Helsinki University of Technology Institute of Mathematics Research Reports Espoo 2010 A582 ANALYSIS OF H(DIV)-CONFORMING FINITE ELEMENTS FOR THE BRINMAN PROBLEM Juho önnö Rolf Stenberg ABTENILLINEN OREAOULU
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 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 informationEnergy norm a-posteriori error estimation for divergence-free discontinuous Galerkin approximations of the Navier-Stokes equations
INTRNATIONAL JOURNAL FOR NUMRICAL MTHODS IN FLUIDS Int. J. Numer. Meth. Fluids 19007; 1:1 [Version: 00/09/18 v1.01] nergy norm a-posteriori error estimation for divergence-free discontinuous Galerkin approximations
More informationAn 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 informationA u + b u + cu = f in Ω, (1.1)
A WEIGHTED H(div) LEAST-SQUARES METHOD FOR SECOND-ORDER ELLIPTIC PROBLEMS Z. CAI AND C. R. WESTPHAL Abstract. This paper presents analysis of a weighted-norm least squares finite element method for elliptic
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 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 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 informationAdaptive approximation of eigenproblems: multiple eigenvalues and clusters
Adaptive approximation of eigenproblems: multiple eigenvalues and clusters Francesca Gardini Dipartimento di Matematica F. Casorati, Università di Pavia http://www-dimat.unipv.it/gardini Banff, July 1-6,
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 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 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 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 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 informationDownloaded 07/25/13 to Redistribution subject to SIAM license or copyright; see
SIAM J. NUMR. ANAL. Vol. 39, No. 6, pp. 2034 2044 c 2002 Society for Industrial and Applied Mathematics RSIDUAL-BASD A POSRIORI RROR SIMA FOR A NONCONFORMING RISSNR MINDLIN PLA FINI LMN CARSN CARSNSN Abstract.
More informationA Posteriori Error Estimates for Weak Galerkin Finite Element Methods for Second Order Elliptic Problems
DOI 10.1007/s10915-013-9771-3 A Posteriori Error Estimates for Weak Galerkin Finite Element Methods for Second Order Elliptic Problems Long Chen Junping Wang Xiu Ye Received: 29 January 2013 / Revised:
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 informationc 2008 Society for Industrial and Applied Mathematics
SIAM J. NUMER. ANAL. Vol. 46, No. 3, pp. 640 65 c 2008 Society for Industrial and Applied Mathematics A WEIGHTED H(div) LEAST-SQUARES METHOD FOR SECOND-ORDER ELLIPTIC PROBLEMS Z. CAI AND C. R. WESTPHAL
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 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 informationDISCRETE EXTENSION OPERATORS FOR MIXED FINITE ELEMENT SPACES ON LOCALLY REFINED MESHES
DISCRETE EXTENSION OPERATORS FOR MIXED FINITE ELEMENT SPACES ON LOCALLY REFINED MESHES MAR AINSWORTH, JOHNNY GUZMÁN, AND FRANCISCO JAVIER SAYAS Abstract. The existence of uniformly bounded discrete extension
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 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 informationUNIFIED A POSTERIORI ERROR ESTIMATOR FOR FINITE ELEMENT METHODS FOR THE STOKES EQUATIONS
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 10, Number 3, Pages 551 570 c 013 Institute for Scientific Computing and Information UNIFIED A POSTERIORI ERROR ESTIMATOR FOR FINITE ELEMENT
More 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 informationANALYSIS OF A LINEAR LINEAR FINITE ELEMENT FOR THE REISSNER MINDLIN PLATE MODEL
Mathematical Models and Methods in Applied Sciences c World Scientific Publishing Company ANALYSIS OF A LINEAR LINEAR FINITE ELEMENT FOR THE REISSNER MINDLIN PLATE MODEL DOUGLAS N. ARNOLD Department of
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 informationOn an Approximation Result for Piecewise Polynomial Functions. O. Karakashian
BULLETIN OF THE GREE MATHEMATICAL SOCIETY Volume 57, 010 (1 7) On an Approximation Result for Piecewise Polynomial Functions O. arakashian Abstract We provide a new approach for proving approximation results
More informationICES REPORT Analysis of the DPG Method for the Poisson Equation
ICES REPORT 10-37 September 2010 Analysis of the DPG Method for the Poisson Equation by L. Demkowicz and J. Gopalakrishnan The Institute for Computational Engineering and Sciences The University of Texas
More informationA SECOND ELASTICITY ELEMENT USING THE MATRIX BUBBLE WITH TIGHTENED STRESS SYMMETRY
A SECOND ELASTICITY ELEMENT USING THE MATRIX BUBBLE WITH TIGHTENED STRESS SYMMETRY J. GOPALAKRISHNAN AND J. GUZMÁN Abstract. We presented a family of finite elements that use a polynomial space augmented
More informationc 2007 Society for Industrial and Applied Mathematics
SIAM J. NUMR. ANAL. Vol. 45, No. 1, pp. 68 82 c 2007 Society for Industrial and Applied Mathematics FRAMWORK FOR TH A POSTRIORI RROR ANALYSIS OF NONCONFORMING FINIT LMNTS CARSTN CARSTNSN, JUN HU, AND ANTONIO
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 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 informationFinite element approximation on quadrilateral meshes
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2001; 17:805 812 (DOI: 10.1002/cnm.450) Finite element approximation on quadrilateral meshes Douglas N. Arnold 1;, Daniele
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 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 informationDerivation and Analysis of Piecewise Constant Conservative Approximation for Anisotropic Diffusion Problems
Derivation Analysis of Piecewise Constant Conservative Approximation for Anisotropic Diffusion Problems A. Agouzal, Naïma Debit o cite this version: A. Agouzal, Naïma Debit. Derivation Analysis of Piecewise
More informationDiscontinuous Galerkin Methods: Theory, Computation and Applications
Discontinuous Galerkin Methods: Theory, Computation and Applications Paola. Antonietti MOX, Dipartimento di Matematica Politecnico di Milano.MO. X MODELLISTICA E CALCOLO SCIENTIICO. MODELING AND SCIENTIIC
More 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 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 informationFIRST-ORDER SYSTEM LEAST SQUARES (FOSLS) FOR PLANAR LINEAR ELASTICITY: PURE TRACTION PROBLEM
SIAM J. NUMER. ANAL. c 1998 Society for Industrial Applied Mathematics Vol. 35, No. 1, pp. 320 335, February 1998 016 FIRST-ORDER SYSTEM LEAST SQUARES (FOSLS) FOR PLANAR LINEAR ELASTICITY: PURE TRACTION
More 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 informationANALYSIS OF MIXED FINITE ELEMENTS FOR LAMINATED COMPOSITE PLATES
ANALYSIS OF MIXED FINITE ELEMENTS FOR LAMINATED COMPOSITE PLATES F. Auricchio Dipartimento di Meccanica Strutturale Università di Pavia, Italy C. Lovadina Dipartimento di Ingegneria Meccanica e Strutturale
More informationTwo Nonconforming Quadrilateral Elements for the Reissner-Mindlin Plate
Two Nonconforming Quadrilateral Elements for the Reissner-Mindlin Plate Pingbing Ming and Zhong-ci Shi Institute of Computational Mathematics & Scientific/Engineering Computing, AMSS, Chinese Academy of
More informationSchur Complements on Hilbert Spaces and Saddle Point Systems
Schur Complements on Hilbert Spaces and Saddle Point Systems Constantin Bacuta Mathematical Sciences, University of Delaware, 5 Ewing Hall 976 Abstract For any continuous bilinear form defined on a pair
More informationA UNIFORMLY ACCURATE FINITE ELEMENT METHOD FOR THE REISSNER MINDLIN PLATE*
SIAM J. NUMER. ANAL. c 989 Society for Industrial and Applied Mathematics Vol. 26, No. 6, pp. 5, December 989 000 A UNIFORMLY ACCURATE FINITE ELEMENT METHOD FOR THE REISSNER MINDLIN PLATE* DOUGLAS N. ARNOLD
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 informationMixed Finite Element Methods. Douglas N. Arnold, University of Minnesota The 41st Woudschoten Conference 5 October 2016
Mixed Finite Element Methods Douglas N. Arnold, University of Minnesota The 41st Woudschoten Conference 5 October 2016 Linear elasticity Given the load f : Ω R n, find the displacement u : Ω R n and the
More informationA PROJECTION-BASED ERROR ANALYSIS OF HDG METHODS BERNARDO COCKBURN, JAYADEEP GOPALAKRISHNAN, AND FRANCISCO-JAVIER SAYAS
A PROJECTION-BASED ERROR ANALYSIS OF HDG METHODS BERNARDO COCBURN, JAYADEEP GOPALARISHNAN, AND FRANCISCO-JAVIER SAYAS Abstract. We introduce a new technique for the error analysis of hybridizable discontinuous
More 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 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 Second Elasticity Element Using the Matrix Bubble
Portland State University PDXScholar Mathematics and Statistics Faculty Publications and Presentations Fariborz Maseeh Department of Mathematics and Statistics 2011 A Second Elasticity Element Using the
More informationETNA Kent State University
Electronic Transactions on Numerical Analysis. Volume 11, pp. 1-24, 2000. Copyright 2000,. ISSN 1068-9613. ETNA NEUMANN NEUMANN METHODS FOR VECTOR FIELD PROBLEMS ANDREA TOSELLI Abstract. In this paper,
More informationWEAK GALERKIN FINITE ELEMENT METHOD FOR SECOND ORDER PARABOLIC EQUATIONS
INERNAIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 13, Number 4, Pages 525 544 c 216 Institute for Scientific Computing and Information WEAK GALERKIN FINIE ELEMEN MEHOD FOR SECOND ORDER PARABOLIC
More informationLeast-squares Finite Element Approximations for the Reissner Mindlin Plate
NUMERICAL LINEAR ALGEBRA WITH APPLICATIONS Numer. Linear Algebra Appl., 6, 479 496 999 Least-squares Finite Element Approximations for the Reissner Mindlin Plate Zhiqiang Cai, Xiu Ye 2 and Huilong Zhang
More informationNumerische Mathematik
Numer. Math. (2003) 94: 195 202 Digital Object Identifier (DOI) 10.1007/s002110100308 Numerische Mathematik Some observations on Babuška and Brezzi theories Jinchao Xu, Ludmil Zikatanov Department of Mathematics,
More informationBUBBLE STABILIZED DISCONTINUOUS GALERKIN METHOD FOR STOKES PROBLEM
BUBBLE STABILIZED DISCONTINUOUS GALERKIN METHOD FOR STOKES PROBLEM ERIK BURMAN AND BENJAMIN STAMM Abstract. We propose a low order discontinuous Galerkin method for incompressible flows. Stability of the
More informationPolynomial-degree-robust a posteriori estimates in a unified setting for conforming, nonconforming, discontinuous Galerkin, and mixed discretizations
Polynomial-degree-robust a posteriori estimates in a unified setting for conforming, nonconforming, discontinuous Galerkin, and mixed discretizations Alexandre Ern, Martin Vohralík To cite this version:
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 informationA NEW CLASS OF MIXED FINITE ELEMENT METHODS FOR REISSNER MINDLIN PLATES
A NEW CLASS OF IXED FINITE ELEENT ETHODS FOR REISSNER INDLIN PLATES C. LOVADINA Abstract. A new class of finite elements for Reissner indlin plate problem is presented. The family is based on a modified
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 informationA posteriori error estimation for elliptic problems
A posteriori error estimation for elliptic problems Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in
More 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 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 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 information1. Introduction. We consider the model problem that seeks an unknown function u = u(x) satisfying
A SIMPLE FINITE ELEMENT METHOD FOR LINEAR HYPERBOLIC PROBLEMS LIN MU AND XIU YE Abstract. In this paper, we introduce a simple finite element method for solving first order hyperbolic equations with easy
More informationInterior superconvergence in mortar and non-mortar mixed finite element methods on non-matching grids
Interior superconvergence in mortar and non-mortar mixed finite element methods on non-matching grids Gergina Pencheva a,, Ivan Yotov b a Center for Subsurface Modeling, Institute for Computational Engineering
More informationINNOVATIVE FINITE ELEMENT METHODS FOR PLATES* DOUGLAS N. ARNOLD
INNOVATIVE FINITE ELEMENT METHODS FOR PLATES* DOUGLAS N. ARNOLD Abstract. Finite element methods for the Reissner Mindlin plate theory are discussed. Methods in which both the tranverse displacement and
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 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 informationNumerische Mathematik
Numer. Math. (2012) 122:61 99 DOI 10.1007/s00211-012-0456-x Numerische Mathematik C 0 elements for generalized indefinite Maxwell equations Huoyuan Duan Ping Lin Roger C. E. Tan Received: 31 July 2010
More informationA LAGRANGE MULTIPLIER METHOD FOR ELLIPTIC INTERFACE PROBLEMS USING NON-MATCHING MESHES
A LAGRANGE MULTIPLIER METHOD FOR ELLIPTIC INTERFACE PROBLEMS USING NON-MATCHING MESHES P. HANSBO Department of Applied Mechanics, Chalmers University of Technology, S-4 96 Göteborg, Sweden E-mail: hansbo@solid.chalmers.se
More 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 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 informationHybridized DG methods
Hybridized DG methods University of Florida (Banff International Research Station, November 2007.) Collaborators: Bernardo Cockburn University of Minnesota Raytcho Lazarov Texas A&M University Thanks:
More 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 information