UNIFIED A POSTERIORI ERROR ESTIMATOR FOR FINITE ELEMENT METHODS FOR THE STOKES EQUATIONS
|
|
- Erik Evan Crawford
- 5 years ago
- Views:
Transcription
1 INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 10, Number 3, Pages c 013 Institute for Scientific Computing and Information UNIFIED A POSTERIORI ERROR ESTIMATOR FOR FINITE ELEMENT METHODS FOR THE STOKES EQUATIONS JUNPING WANG, YANQIU WANG, AND XIU YE Abstract. This paper is concerned with residual type a posteriori error estimators for finite element methods for the Stokes equations. In particular, the authors established a unified approach for deriving and analyzing a posteriori error estimators for velocity-pressure based finite element formulations for the Stokes equations. A general a posteriori error estimator was presented with a unified mathematical analysis for the general finite element formulation that covers conforming, non-conforming, and discontinuous Galerkin methods as examples. The key behind the mathematical analysis is the use of a lifting operator from discontinuous finite element spaces to continuous ones for which all the terms involving jumps at interior edges disappear. Key words. A posteriori error estimate, finite element methods, Stokes equations 1. Introduction A posteriori error estimator refers to a computable formula that offers a measure for judging the reliability and efficiency of a particular numerical scheme employed for approximating the solution of partial differential equations or alike. With a mathematically justified a posteriori error estimator, one would be able to generate a mesh that is tailored at reducing computational errors at places of great need. This process is commonly known as adaptive mesh refinement which has become a useful and important tool in today s scientific and engineering computing. The goal of this paper is to offer a systematic framework for developing and analyzing a posteriori error estimators for finite element methods for model Stokes equations. This paper is concerned with residual type a posteriori error estimators. In other words, the computable formula for judging the efficiency and reliability of numerical schemes shall be given by functions of residuals. Along this avenue, several fine results have been developed for finite element methods for the Stokes equations. For conforming finite element methods, some a posteriori error estimators have been derived for mini-elements by Verfurth [1] and Bank-Welfert [4]. Ainsworth- Oden [3] and Nobile [18] have considered more general conforming finite elements in their study. For nonconforming finite elements, a posteriori error estimation for the Crouzeix-Raviart element [8] has been developed by several researchers such as Verfurth [], Dari-Durán-Padra [9] and Doerfler-Ainsworth [10]. Carstensen, Gudi, and Jensen [5] proposed and analyzed an a posteriori error estimator for discontinuous Galerkin methods by using a stress-velocity-pressure formulation for the Stokes equations. Kay and Silvester[16] established a posteriori error estimation Received by the editors August 1, Mathematics Subject Classification. 65N15, 65N30, 76D07, 35B45, 35J50. The research of Junping Wang was supported by the NSF IR/D program, while working at the Foundation. However, any opinion, finding, and conclusions or recommendations expressed in this material are those of the author and do not necessarily reflect the views of the National Science Foundation. The research of Xiu Ye was supported in part by National Science Foundation Grant DMS
2 55 J. WANG, Y. WANG, AND X. YE for the stabilized finite element formulation. The recovery based a posteriori error estimate for the Stokes equations is investigated in [1]. In both [9] and [10], the analysis for their a posteriori error estimators was based on a Helmholtz decomposition for decomposing the Crouzeix-Raviart element into two parts: an exactly divergence-free part and the second as its orthogonal complement. While the Helmholtz decomposition offers an applaudable approach for analyzing the efficiency and reliability of a posteriori error estimators for the Stokes equations, the method has difficulty in being extended to finite element approximations arising from discontinuous Galerkin methods. The main difficulty comes from the fact that the approximate velocity field from the discontinuous finite element methods is not divergence-free in the classical sense. Therefore, other analytical techniques have been developed for discontinuous finite elements; but most of them requires special and unnecessary properties about the finite element mesh. For example, Houston, Schötzau and Wihler [14] have developed an a posteriori error analysis for the discontinuous Q k Q k 1 element on partitions consisting of parallelograms only. In this paper, we establish a unified approach for deriving and analyzing a posteriori error estimators of residual type for velocity-pressure based formulations of the Stokes equations. In particular, we shall develop a general finite element formulation that covers conforming, non-conforming, and discontinuous Galerkin methods as examples. Then, a general a posteriori error estimator shall be presented with a unified mathematical analysis. The key behind the analysis is the use of a lifting operator from discontinuous finite element spaces to continuous ones for which all the terms involving jumps at interior edges disappear. A similar lifting operator was employed by Karakashian and Pascal [15] for analyzing a posteriori error estimates for a discontinuous Galerkin approximation to second order elliptic equations. The paper is organized as follows. In Section, a model Stokes problem and some notations are introduced. In Section 3, we shall first present a general finite element formulation for the Stokes equations, and then illustrate how most existing conforming, nonconforming, and discontinuous Galerkin methods be represented by the general framework. In Section 4, we establish an analytical tool for analyzing the general a posteriori error estimator of residual type. Finally in Section 5, we present some numerical results to confirm the theory developed in previous sections.. Preliminaries and notations Let Ω be an open bounded domain in R d,d =,3. Denote by Ω the boundary of Ω. The model problem seeks a velocity function u and a pressure function p satisfying (1) () (3) u+ p = f in Ω, u = 0 in Ω, u = 0 on Ω, where,, and denote the Laplacian, gradient, and divergence operators, respectively, and f is the external volumetric force acting on the fluid. For simplicity, the algorithm and its analysis will be presented for the model Stokes problem (1)-(3) only in two-dimensional spaces (i.e.; d = ) with polygonal domains. An extension to the Stokes problem in three dimensions can be made formally for general polyhedral domains.
3 UNIFIED A POSTERIORI ERROR ESTIMATE FOR STOKES EQUATIONS 553 For any given polygon D Ω, we use the standard definition of Sobolev spaces H s (D) with s 0 (e.g., see [1, 6] for details). The associated inner product, norm, and seminorms in H s (D) are denoted by (, ) s,d, s,d, and r,d,0 r s, respectively. When s = 0, H 0 (D) coincides with the space of square integrable functions L (D). In this case, the subscript s is suppressed from the notation of norm, semi-norm, and inner products. Furthermore, the subscript D is also suppressed when D = Ω. Denote by L 0 (D) the subspace of L (D) consisting of functions with mean value zero. The above definition/notation can easily be extended to vector-valued and matrixvalued functions. The norm, semi-norms, and inner-product for such functions shall follow the same naming convention. In addition, all these definitions can be transferred from a polygonal domain D to an edge e, a domain with lower dimension. Similar notation system will be employed. For example, s,e and e would denote the norm in H s (e) and L (e) etc. Throughout the paper, we follow the convention that a bold Latin letter denotes a vector. Let u = [u i ] 1 i, v = [v i ] 1 i be two vectors, and σ = [σ ij ] 1 i,j, τ = [τ ij ] 1 i,j be two matrices, define ) v = u v = v σ = ( v1 x v x It is not hard to see that v 1 y v y, v = v 1 x + v y, ( ) u1 v 1 u 1 v, σ : τ = σ u v 1 u v ij τ ij, ( ) σ11 v 1 +σ 1 v, v σ u = σ 1 v 1 +σ v σ : (u v) = v σ u. i,j=1 σ ij u i v j. i,j=1 Let T h be a geometrically conformal triangulation of the domain Ω; i.e., the intersection of any two triangles in T h is either empty, a common vertex, or a common edge. Denote by h T the diameter of triangle T T h, and h the maximum of all h T. We assume that T h is shape regular in the sense that for each T T h, the ratio between h T and the diameter of the inscribed circle is bounded from above. The shape regularity of T h ensures a validity of the inverse inequality for finite element functions. In addition, shape regularity allows one to apply the routine scaling arguments in finite element analysis. Let us introduce two finite dimensional spaces V s and P t as follows: V s = {v V : v T [P s (T)], for all T T h }, P t = {q L 0 (Ω) : q T P t (T), for all T T h }, wherev = [L (Ω)] or [H 1 0(Ω)], s and tarenon-negativeintegers, and P k (T)is the set of polynomials of degree no more than k on T. A finite element method usually seeks discrete velocity and pressure approximations in some subspaces V h V s and Q h Q t. Certain continuity condition may be imposed on V h and Q h, depending on the type of finite elements. For a posteriori error estimates only, V h and Q h are not required to satisfy the discrete inf-sup condition. However, we do need to make the following assumption (4) {v [H 1 0 (Ω)] : v T [P 1 (T)], for all T T h } V h,
4 554 J. WANG, Y. WANG, AND X. YE in the analysis of the a posteriori error estimates. For a finite element partition T h, with triangles only, (4) is satisfied as long as k, order of the polynomials, is no less than 1. Therefore, the assumption (4) is reasonable and examples that satisfy it will be given later. Denote by E h the set of all edges in T h, and denote by E 0 h := E h\ Ω the collection of all interior edges. Next, we define the average and jump on edges for scalar-valued function q, vector-valued function w, and matrix-valued function τ, respectively. For any interior edge e E 0 h, let T 1, T be two triangles sharing e and n 1 and n be the unit outward normal vectors on e, associated with T 1 and T, respectively. Define the average { } and jump [ ] on e by {q} = 1 (q T 1 +q T ), [q] = q T1 n 1 +q T n, {w} = 1 (w T 1 +w T ), [w] = w T1 n 1 +w T n, {τ} = 1 (τ T 1 +τ T ), [τ] = n 1 τ T1 +n τ T. We also define a matrix-valued jump [ ] for w on e by [w] = w T1 n 1 +w T n. If e is a boundary edge, the above definitions need to be adjusted accordingly so that both the average and the jump are equal to the one-sided values on e. That is, {q} = q e, {w} = w e, {τ} = τ e, [q] = q e n, [w] = w e n, [τ] = n τ e, [w] = w e n, where n is the unit outward normal of Ω. Let q, v and τ be scalar-, vector-, and matrix-valued functions that are regular enough to make all involving terms well-defined, then the following identities are standard []: qv nds = [q] {v}ds+ (5) {q}[v]ds, T T T h e Eh 0 e e E e h n τ vds = [τ] {v}ds+ (6) {τ} : [v]ds. T T T h e Eh 0 e e E e h 3. Finite element formulation We first derive a general weak formulation for the Stokes equations that covers a set of existing numerical methods including discontinuous Galerkin, conforming finite elements, and some nonconforming finite element schemes. To this end, let v and w be vector-valued functions, and we introduce the following notation (v, w) Th := v w dx, K T h K (v, w) Eh := e E h e v w ds. Testing the momentum equation (1) by v V h yields (7) ( u, v)+( p, v) = (f, v).
5 UNIFIED A POSTERIORI ERROR ESTIMATE FOR STOKES EQUATIONS 555 Next, using integration by parts, Equation (6) and the fact that [ u] = 0 on all e Eh 0, we have ( u, v) = ( u, v) Th n u vds K T h K = ( u, v) Th ({ u}, [v]) Eh. Analogously, since [p] = 0 on all e Eh 0, then by using Equation (5) we arrive at ( p, v) = ( v, p) Th + v np ds Thus, (7) can be rewritten as K T h K = ( v, p) Th +({p}, [v]) Eh. (8) ( u, v) Th ({ u}, [v]) Eh ( v, p) Th +({p}, [v]) Eh = (f, v). The mass conservation equation () can be tested by using any q Q h, yielding (9) ( u, q) Th = 0. In order to derive a general weak formulation, we introduce the following bilinear forms: a(u,v) := ( u, v) Th ( (10) γ ({ u}, [v]) Eh +δ({ v}, [u]) Eh α(h 1 e [u], [v]) Eh ), (11) b(v,p) := ( v, p) Th +γ([v],{p}) Eh, where γ = 0 or 1, δ = 1, 1 or 0, and α 0 are parameters with various values. The general weak formulation for the Stokes equations is then given as follows: (1) (13) Algorithm G: Find (u h ;p h ) V h Q h such that a(u h,v)+b(v,p h ) = (f,v), b(u h, q) = 0 for all (v;q) V h Q h. When γ = 1, it is not hard to see that system (1)-(13) is consistent with the system (8)-(9), as the exact solution of the Stokes problem satisfies [u] = 0 and [u] = 0 on all e E h. When γ = 0, these two systems are in general not consistent. It should be pointed out that by choosing different values of γ and the spaces V h and Q h, Algorithm G represents different types of finite element formulations. For illustrative purpose, we shall list three most important examples; all of which satisfy the assumption (4) on V h. Example 1: Discontinuous Galerkin. Set γ = 1 and V h {v [L (Ω)] : v T [P s (T)], for all T T h }, Q h {q L 0(Ω) : q T P t (T), for all T T h }. When δ = 1, the corresponding formulation is symmetric. One example of such a discontinuous Galerkin formulation was discussed in [5], where V h, Q h are chosen to be H(div) conforming elements and the well-known inf-sup condition was inherited to be valid.
6 556 J. WANG, Y. WANG, AND X. YE Example : Conforming finite element method. Set γ = 1 and V h {v [H 1 0 (Ω)] : v T [P s (T)], for all T T h }, Q h {q L 0(Ω) : q T P t (T), for all T T h }. By the continuity requirement on V h, we easily see that the two bilinear forms are simplified to a(v,w) = ( v, w), b(v,q) = ( v, q). The conforming finite element scheme is one of the most well-studied formulations for the Stokes problems. Later in the section for numerical experiments, we shall consider one such element, namely, the Taylor-Hood element [13]. Example 3: Nonconforming finite element method. Set γ = 0 and chose a finite element space V h such that V h [H 1 0 (Ω)]. This gives some nonconforming finite element methods, depending on the selection of V h and Q h. In this paper, we only consider nonconforming finite elements which satisfies the following condition: Assumption (H): ({ u pi}, [v]) Eh = 0 for all u, v V h, p Q h, where I is the identity matrix. It is not hard to see that the Crouzeix- Raviart type nonconforming elements [8, 11, 7] satisfy this assumption, because by its definition, { u pi} is a polynomial one degree lower than [v] and [v] vanishes at all Gaussian points on each edge e E h. The goal of this manuscript is to provide a tool that can be employed to analyze a posteriori error estimators for the unified formulation (1)-(13). To this end, we first derive an orthogonality property of the error between the exact solution and its finite element approximation: e = u u h, ǫ = p p h. By subtracting (1) from (8) and using ({p h }, [v]) Eh = ({p h I}, [v]) Eh, we have ( ) ( e, v) Th ({ e}, [v]) Eh +γ δ({ v}, [u h ]) Eh α(h 1 e [u h], [v]) Eh ( v,ǫ) Th +({ǫ}, [v]) Eh (1 γ)({ u h p h I}, [v]) Eh = 0. Note that for all of the above-mentioned discontinuous Galerkin, conforming and nonconforming cases, we must have (14) (1 γ)({ u h p h I}, [v]) Eh = 0. In fact, (14) follows from γ = 1 for the discontinuous Galerkin and the conforming cases, and from Assumption (H) for the nonconforming case. Combining this and the fact that [u] = 0 on all edges, the above equation can be rewritten as (15) ( ) ( e, v) Th ({ e}, [v]) Eh γ δ({ v}, [e]) Eh α(h 1 e [e], [v]) Eh ( v,ǫ) Th +({ǫ}, [v]) Eh = 0. Moreover, if v V h [H 1 0 (Ω)], then (15) becomes (16) ( e, v) Th +γδ({ v}, [u h ]) Eh ( v,ǫ) Th = 0.
7 UNIFIED A POSTERIORI ERROR ESTIMATE FOR STOKES EQUATIONS 557 Next, we introduce a norm for the space [H0 1(Ω)] +V h as follows (17) v = v 1,T + h 1 e [v] e. e E h For a well-crafted numerical scheme, we usually expect to have an a priori error estimate like the following (18) e + ǫ Ch k+1 ( u k+1 + p k ), wherec isapositiveconstantandk isdeterminedbythe orderofthecorresponding finite elements and the regularity of the exact solution (u; p). For some elements, the a priori error estimate (18) may have variations on the right-hand side, but this does not affect our analysis to be presented. A priori error estimation like (18) is well known for conforming finite element methods, for which v = v 1. Such an estimate is also known for some discontinuous Galerkin formulations; e.g., see [4, 5] for an approach with H(div)-elements. For nonconforming finite elements, such as the Crouzeix-Raviart type elements, similar error estimates hold true due to the fact that the jump term, [u h ], vanishes at all Gaussian points on the edge e. In fact, for any interior edge e shared by two elements T 1 and T, let ē be the average value of e on edge e. Since [e] = [u h ] vanishes at all Gaussian points on the edge e, then ē has the same value when the trace of e was taken from either T 1 or T. Hence, he 1 [e] e = h 1 e [e ē] e Consequently, we have e E h h 1 h 1 e e ē e T i T 1,T C e 1,T i. T 1,T e [e] e C e 1,T. Therefore, e is in the same order as e 1,T for the Crouzeix-Raviart type elements. The analysis here shows that the norm as defined in (17) is a reasonable one to use in the a priori and later the a posteriori error estimates for the unified formulation (1)-(13). 4. A posteriori error estimation The goal of this section is to derive an a posteriori error estimation for the Algorithm G by using a unified framework. For simplicity of notation, we shall use to denote less than or equal to up to a constant independent of the mesh size, variables, or other parameters appearing in the inequality. To this end, define { [ uh p J 1 ( u h p h I) = h I] on e Eh 0, 0 on boundary edges, and { [uh ] on e E J (u h ) = h 0, u h n on boundary edges. Our residual-based global error estimator is given by (19) η = η T,
8 558 J. WANG, Y. WANG, AND X. YE where ηt =h T f h + u h p h T + u h T + 1 ( he J 1 ( u h p h I) +h 1 e J (u h ) ) ds, e E h e with h e being the length of edge e, and f h the L projection of the load function f onto V h. It is also convenient to introduce an oscillation quantity for the load function f on T h as follows ( ) 1/ osc(f) = h T f f h T. Our ultimate goal is to establish the following result. Theorem 1. Let (u;p) be the solution of (1)-(), and (u h ;p h ) be its finite element approximation arising from (1)-(13). Then, one has ( ) 1/ (0) (u u h ) T + p p h η +osc(f) and (1) η ( (u u h ) T ) 1/ + p p h +osc(f). For convenience, the relation (0) is referred to as a reliability estimate and (1) as an efficiency estimate A lifting operator and some technical estimates. Define H 1 (T h ) = H 1 (T) and S k = P k (T). For any triangle T T h, denote by T (T) the set of all triangles in T h having a nonempty intersect with T, including T itself. Denote by E(T) the set of all edges in E h having a nonempty intersection with T, including all three edges of T. Similarly, for any point x in Ω, denote by E(x) the set of all edges that pass through x. Note that E(x) is nonempty only when x lies on E h. Let e be an edge of triangle T. For any v S k, let us define a lifting operator () L k : v L k (v) S k H 1 0 (Ω) that lifts a discontinuous piecewise polynomial to a continuous piecewise polynomial with vanishing boundary trace as follows. Let G k (T) be the set of all Lagrangian nodal points for P k (T). At all internal Lagrangian nodal points x j G k (T), we set L k (v)(x j ) = v(x j ). At boundary Lagrangian points x j G k (T) T, we let L k (v)(x j ) be either a traceof v from anyside oraprescribed weighted averageof all possible traces. At global Lagrangian points x j T Ω, we set L k (v)(x j ) = 0. Note that for any function g H 1 (T) the following estimate holds: (3) g e h 1 T g T +h T g T. Lemma 1. For any v S k, k 1, the lifting operator L k as defined in () satisfies the following estimate: (4) v L k (v) T +h T (v L k(v)) T h e [v] e T T h. e E(T)
9 UNIFIED A POSTERIORI ERROR ESTIMATE FOR STOKES EQUATIONS 559 Proof. The proof follows from a routine scaling argument and the fact that all norms on finite dimensional spaces are equivalent. To this end, we observe that v L k (v) T +h T (v L k (v)) T h T (v L k (v))(x j ), x j G k (T) where G k (T) is the set of all Lagrangian nodal points for P k (T). For simplicity, let L k (v) be defined by taking a random one-sided trace on the boundary Lagrangian points. It follows from the definition of L k that v L k (v) vanishes at all internal Lagrangian points in T. At Lagrangian points on T, there are two possibilities: (1) x j is in the interior of an edge e; () x j is a vertex of T. In the first case, we see that (v L k (v))(x j ) is either 0 or [v] e (x j ), where [v] e denotes the jump on e. In the second case, we can use the triangle inequality to traverse through all edges e E(x j ) and to obtain (v L k (v))(x j ) [v] e (x j ). e E(x j) The above analysis, together with the shape regularity of T h, implies that (v L k (v))(x j ) [v] e (x j ), x j G k (T) e E(T) x j G k (e) where G k (e) was used to denote the corresponding Lagrangian points on edge e. Then, using the routine scaling argument on edges, inequality (4) follows immediately. Thefollowingisanotherresultthatturnsouttobe veryuseful intheforthcoming analysis. Lemma. For any v H 1 (T h ), there exists a v I S k H0(Ω), 1 k 1, satisfying (5) v v I T +h T (v v I ) T h T v T + h e [v] e T T h. T T (T) e E(T) Proof. First of all, there exists a piecewise constant v 0 S k (e.g., the cell average of v) such that (6) v v 0 T +h T (v v 0 ) T h T v T T T h. Furthermore, by the approximation property, inequality (3) and the shape regularity of T h, we have for each e Eh 0 h e [v 0 ] e = h e v 0 T1 n 1 +v 0 T n ds e (7) h e (v 0 v) T1 n 1 ds+h e (v 0 v) T n ds+h e [v] e e h T 1 v T 1 +h T v T +h e [v] e, where T 1 and T are the two triangles sharing e as a common edge. On boundary edges, similar result can obviously be obtained without any difficulty. Since v 0 S k with k = 1, according to Lemma 1, one may lift v 0 to a continuous piecewise linear function L 1 (v 0 ) S 1 H0 1 (Ω) which satisfies (4). By taking e
10 560 J. WANG, Y. WANG, AND X. YE v I = L 1 (v 0 ), we obtain from the usual triangle inequality, (6), and (7) that v v I T +h T (v v I) T v v 0 T +h T (v v 0) T + v 0 L 1 (v 0 ) T +h T (v 0 L 1 (v 0 )) T h T v T + h e [v 0 ] e e E(T) h T v T + h e [v] e, T T(T) e E(T) which completes the proof. 4.. Reliability estimate. We first establish an estimate for the pressure error in terms of the a priori error estimator and the velocity error. The result can be stated as follows. Lemma 3. Let (u;p) be the solution of (1)-() and (u h ;p h ) be its finite element approximation arising from (1)-(13). Then, we have (8) p p h η +osc(f)+ ( (u u h ) T ) 1/. Proof. Let v [H 1 0(Ω)] and v I V h [H 1 0(Ω)] be an interpolation of v such that both components satisfy (5). Observe that such an interpolation v I is possible if V h satisfiesthe assumption(4). Notethat [v] = 0oneveryedgesincev [H 1 0(Ω)]. Let e = u u h and ǫ = p p h be the error for velocity and pressure approximations, respectively. Using integration by parts, (6), (16), (3), (5), and the fact that both v and v I are continuous across each interior, we arrive at ( v, ǫ) = ( (v v I ), ǫ)+( v I, ǫ) = ( (v v I ), ǫ)+( e, v I ) Th +γδ({ v I }, [e]) Eh = ( (v v I ), ǫi e) Th +( e, v) Th +γδ({ v I }, [e]) Eh = (f + u h p h, v v I ) Th +({v v I }, [ u h p h I]) E 0 h γδ({ v I }, [u h ]) Eh +( e, v) Th ( v 1 ( h T f + u h p h T) 1/ +( h e [ u h p h I] e) 1/ e Eh 0 +( e [u h ] e) 1/ +( ) e T) 1/, e E h h 1 which, together with the following inf-sup condition ( v, p p h ) (9) p p h sup v [H0 1 v (Ω)] 1 yields the required estimate (8). The next result is concerned with an estimate for the velocity approximation, which can be stated as follow.
11 UNIFIED A POSTERIORI ERROR ESTIMATE FOR STOKES EQUATIONS 561 Theorem. Let (u;p) and (u h ;p h ) be the solutions of (1)-() and (1)-(13). Then one has the following global reliability estimate: (30) ( (u u h ) T ) 1/ η +osc(f). Substituting (30) into (8) yields the following estimate for the pressure approximation: (31) p p h η +osc(f). Thus, it follows from the definition of that e + ǫ η +osc(f). Proof. Let e I V h [H 1 0 (Ω)] be an interpolation of e satisfying (5). Again, such a choice of e I is possible because V h was assumed to satisfy (4). Observe that [e] = [u h ]. Thus, it follows from (5) and the mesh regularity that (3) e e I T+h T (e e I ) T h T T T (T) e T + e E(T) e J (u h ) e. Using(16),integrationbyparts,(6), andthefactthate I, u, u, parecontinuous across all interior edges, we obtain h 1 ( e, e) Th = ( e, (e e I )) Th +( e, e I ) Th = ( e, (e e I )) Th +( e I, ǫ) γδ({ e I }, [e]) Eh = ( (e e I ), e ǫi) Th +( e, ǫ) Th γδ({ e I }, [e]) Eh = (f + u h p h, e e I ) Th +( e, ǫ) Th γδ({ e I }, [e]) Eh ({e e I }, [ u h p h I]) E 0 h ({ e ǫi}, [u h ]) Eh = (f + u h p h, e e I ) Th +( e, ǫ) Th γδ({ e I }, [e]) Eh (33) ({e e I }, [ u h p h I]) E 0 h ({ e ǫi}, [u h χ]) Eh = I 1 +I +I 3 +I 4 +I 5 ; whereχ = L k (u h ) V h [H 1 0(Ω)] isacontinuousinterpolationofu h suchthat (4) is satisfied, and I j represents the corresponding term from the previous equation. The first term I 1 can be estimated by using the inequality (3) as follows I 1 (34) ( )1 ( h T e e I T h T f + u h p h T ( )1 ( h T f + u h p h T e T + ( ) 1/ (η +osc(f)) η + e T. )1 e E h h 1 e J (u h ) e )1
12 56 J. WANG, Y. WANG, AND X. YE The term I can be handled by using (8) and the fact that e = u h as follows ( ) 1/ I u h T p p h ( ) 1/ (35) η η +osc(f)+ e T. Using [e] = [u h ] and the estimate (3), the term I 3 can be bounded as follows ( )1 ( )1 I 3 γδ h e e I e h 1 e [u h ] e e E h e E h ( )1 ( )1 e I T h 1 e [u h] e e E h ( )1 e T + ( )1 h 1 e [u h] e h 1 e [u h] e e E h e E h ( )1 (36) e T +η η. The same argument can be applied to provide an estimate for the term I 4 : ( )1 (37) I 4 e T +η η. As to I 5, we first use (15) to obtain I 5 =({ e ǫi}, [u h χ]) Eh =({ e}, [u h χ]) Eh ({ǫ}, [u h χ]) Eh =( e, (u h χ)) Th γδ({ (u h χ)}, [e]) Eh +γα(h 1 e [u h χ], [e]) Eh ( (u h χ), ǫ) Th =( e, (u h χ)) Th γδ({ (u h χ)}, [e]) Eh γα(h 1 e [u h], [u h ]) Eh ( (u h χ), ǫ) Th =J 1 +J +J 3 +J 4. Recall that χ = L k (u h ) is a lift of u h so that the estimate (4) holds true. It follows from (4) that ( )1 J 1 η e T Next, we use [e] = [u h ] to obtain J γδ e E h (u h L k (u h )) e [u h ] e, which, with the help of (3) and (4), can be bounded by J η.
13 UNIFIED A POSTERIORI ERROR ESTIMATE FOR STOKES EQUATIONS 563 It is obvious that J 3 η, and it follows from (8) and (4) that ) J 4 η (η +osc(f)+( e, e) 1. Th Collecting all the estimates for J s yields ) (38) I 5 η (η +osc(f)+( e, e) 1. Th Now we substitute all the estimates for I s into (33) to obtain ( ) 1/ ( e, e) Th (η +osc(f)) η +osc(f)+ e T This completes the proof of (30). 3 (η +osc(f)) + 1 ( e, e) T h Efficiency. We first define two bubble functions, which are widely used in a posteriori error estimations [3]. For each triangle T T h, denote by φ T the following bubble function { 7λ 1 λ λ 3 in T, φ T = 0 in Ω\T, where λ i, i = 1,,3 are barycentric coordinates on T. It is clear that φ T H 1 0 (Ω) and satisfies the following properties [3]: (39) (40) For any polynomial q with degree at most m, there exist positive constants c m and C m, depending only on m, such that c m q T q φ T dx q T, (qφ T ) T C m h 1 T q T. T For each e Eh 0, we can analogously define an edge bubble function φ e. Let T 1 and T be two triangles sharing the edge e. To this end, denote by ω e = T 1 T the union of the elements T 1 and T. Assume that in T i, i = 1,, the barycentric coordinates associated with the two ends of e are λ Ti 1 and λ Ti, respectively. The edge bubble function can be defined as follows 4λ T1 1 λt1 in T 1, φ e = 4λ T 1 λt in T, 0 in Ω\ω e. Then φ e H0 1 (Ω) and satisfies the following properties [3]: For any polynomial q with degree at most m, there exist positive constants d m, D m and E m, depending only on m, such that (41) d m q e q φ e ds q e, (4) (43) (qφ e ) ωe D m h 1/ e q e, qφ e ωe E m h 1/ e q e. Then we have the following efficiency bound. e
14 564 J. WANG, Y. WANG, AND X. YE Theorem 3. Let (u;p) and (u h ;p h ) be the solutions of (1)-() and (1)-(13), respectively. Then for all T T h and e Eh 0, we have (44) h T f h + u h p h T e T + ǫ T +h T f f h T, ( ) 1/ (45) h 1/ e J 1 ( u h p h I) e h e f f h ωe + e T + ǫ ωe, T ω e (46) u h T e T. By summing the above estimates over all element T T h and edges e E h, we obtain the following efficiency estimate: η e + ǫ +osc(f). Proof. Let w T = (f h + u h p h )φ T. Since f = u+ p and w T vanishes on T, we clearly have (f f h, w T ) T +(f h + u h p h, w T ) T = ( e, w T ) T ( w T, ǫ ) T. Then, by inequalities (39)-(40), f h + u h p h T (f h + u h p h, w T ) T = ( e, w T ) T ( w T, ǫ ) T (f f h, w T ) T ) ( e T + ǫ T w T T + f f h T w T T ( h 1 T e T +h 1 T ǫ T + f f h T ) f h + u h p h T. This completes the proof of (44). Similarly, for any e Eh 0, let w e = ([ u h p h I])φ e. Using integration by parts and the fact that w e = 0 on ω e, we have (47) ( u h, w e ) ωe = ( u h,w e ) T + [ u h ] w e ds, T ω e e and (48) ( w e, p h ) ωe = ( p h,w e ) T + T ω e e [p h I] w e ds. Testing (1) by using w e over ω e and then using integration by parts give (49) (f, w e ) ωe = ( u, w e ) ωe ( w e, p) ωe. Using the properties of φ e and equations (47)-(49), we have [ u h p h I] e [u h p h I] w e ds e = ) ((f f h, w e ) T +(f h + u h p h, w e ) T ( e, w e ) T +(ǫ, w e ) T T ω e [ u h p h I] e ( h 1/ e f f h ωe +h 1/ e f h + u h p h ωe Combining the above with (44) gives (45). +he 1/ ( ) e T )1/ +he 1/ ǫ ωe T ω e
15 UNIFIED A POSTERIORI ERROR ESTIMATE FOR STOKES EQUATIONS 565 Finally, the estimate (46) holds true as u = 0 and clearly u h T = u h u T = e T e T. This completes the proof of the theorem. 5. Numerical results Some numerical results for such a posteriori error estimators have been reported in [10] for the Crouzeix-Raviart element, and in [4] for an H(div) based discontinuous Galerkin formulation. In this section, we would like to present some computational results for the conforming finite element method in order to verify the theory established in previous Sections. To this end, we consider the Taylor-Hood element [13] when applied to the Stokes problem. We use the BiConjugate Gradient (BICG) iterative solver for solving the resulting linear algebraic equations, with a relative residual = 10 8 being set as a stopping criteria Test problems. Three problems are considered in the numerical test; all are defined on the domain of unit square Ω = (0,1) (0,1). Two of them have exact solutions given by: Test Problem 1: A first Stokes problem with exact solution ( ) x u = y(x 1) (y 1)(y 1) xy (x 1)(x 1)(y 1), π(y x) p = sin, Test Problem : A second Stokes problem with exact solution in the polar coordinate system ( 3 ( )) r cos θ 3θ cos u = 3 ( r 3sin θ sin ), 3θ p = 6r 1/ cos θ. Observe that the solution has a corner singularity of order 0.5 at the origin (0,0). Test Problem 3: D lid driven cavity. This is a Stokes problem which describes the flow of fluid in a rectangular container driven by a uniform motion of the top lid [0]. Note that the boundary condition features a discontinuity at two top corners. The exact solution (u; p), if it exists, should not be sufficiently smooth so that (u; p) [H 1 ] L. However, the discrete problem is still well-posed and should provide an approximation to the actual solution. It must be pointed out that the finite element formulation was only presented for the homogeneous Dirichlet boundary condition in the previous Sections. But the formulation can be easily extended to problems with non-homogeneous Dirichlet boundary conditions such as u = g on Ω. In this case, the term J (u h ) in the a posteriori error estimator needs to be modified as follows { [uh ] on interior edge e E J (u h ) = h 0, (u h g) n on boundary edges. Also, in the computational implementation, we had replaced h T f h+ u h p h T by T f h + u h p h T, where T stands for the area of triangle T.
16 566 J. WANG, Y. WANG, AND X. YE 5.. Results with uniform meshes. In this experiment, we solve test problems 1and on uniform meshesand compute the asymptoticorderofthe errorestimator η. The coarsest mesh is generated by dividing Ω into 4 4 sub-rectangles, and then dividing each sub-rectangle into four triangles by connecting its two diagonal lines. We then use the routine uniform refinement procedure, which divides each triangle into four sub-triangles by connecting the center of its three edges, to generate several levels of fine meshes. For test problem 1, we know theoretically that η = O(h ). For test problem, the order of η is expected to be O(h 0.5 ). Numerical results for these two test problems are reported in tables 1 and. For test problem, the error of pressure is not calculated since the pressure goes to infinity at the point (0,0). Table 1. Error profiles for test problem 1 on uniform triangular meshes. h dofs η (u u h ) u u h p p h e e e e e e e e e-03.71e e e e e-05.13e e-05 Asym. Order O(h k ), k = Table. Error profiles for test problem on uniform triangular meshes. h ndofs η (u u h ) u u h e e e e e e e+00.33e-01.17e e e e-04 Asym. Order O(h k ), k = Results with adaptive refinement. We perform adaptive refinements for test problems and 3, which have corner singularities. Two different refinement strategies are employed in this study. The first one is based on a comparison of each error η T with the maximum value of all the error estimators. The strategy can be described as follows: Local Refinement by Maximum Strategy : (1) Given a current triangular mesh, error indicators η T on each triangle, and a threshold θ (0,1) (e.g., θ = 0.5). One computes the maximum error η max = maxη T. () For each triangle T, if η T θη max, mark this triangle for refinement. (3) The actual refinement is done by the newest node bisection method [17, 19]. It has been proved that this method will not cause mesh degeneration. The only requirement is that the newest nodes for the coarsest mesh must be assigned carefully such that every triangle is compatibly divisible. This can be easily checked.
17 UNIFIED A POSTERIORI ERROR ESTIMATE FOR STOKES EQUATIONS 567 Table 3. Error profiles for test problem using adaptive mesh refinements. Strategy Refinement times dofs η (u u h ) u u h Maximum Local The second refinement strategy is based on a comparison of η T with those for its neighbors. To explain the main idea, let ρ > 0 be a prescribed distance parameter and set T ρ,t = { T : 0 < T T ρ}, where T T stands for the distance of the centers of T and T. With a given threshold θ > 1, we mark a triangle T for refinement if η T θη N(T), where η N(T) is the average of the local error indicator on all the neighboring triangles T T ρ,t. The following is such a refinement strategy that was numerically investigated in this study. Local Refinement by Local Strategy : (1) Given a current triangular mesh, error estimators η T on each triangle, and a threshold θ > 1.0 (e.g., θ = 1.5). One computes an error indicator η N(T) as the average of the local error indicator on neighboring triangles that share a vertex or an edge with T, not including T itself. () For each triangle T, if η T θη N(T), mark this triangle for refinement. (3) The actual refinement is again done by the newest node bisection method [17, 19]. The residual estimator η and errors between the true solution and its finite element approximations for test problem are reported in Table 3. By comparing tables and 3, one clearly sees the power of adaptive refinements in numerical methods for PDEs. For example, the error in L and H 1 norms are given by and with only 1117 degree of freedoms when the local adaptive refinement strategy was used, while the uniform partition requires more than degree of freedoms in order to reach a comparable accuracy. The refined meshes after 0, 8, and 16 refinements are illustrated in Figures 1 and, which show that our residual estimator really captures the corresponding corner singularity correctly, under both refinement strategies. Furthermore, in Figures 3 and 4, we examine the relation of η, (u u h ), and u u h with the degree of freedoms N during the process of the adaptive refinement. The plots start from the coarsest mesh and ends after 16 refinements. For the driven cavity problem, since the exact solution may not be sufficiently as smooth as in [H 1 ] L, the residual error η is not expected to decrease when the mesh is refined. Our numerical experiments show that the error indicator η T is able to locate both corner singularities for this problem. The meshes after 0, 8, and 16 refinements are plotted in Figures 5 and 6. Readers are invited to make their own conclusions for the numerical results illustrated in this Section.
18 568 J. WANG, Y. WANG, AND X. YE Figure 1. Test problem, maximum strategy adaptive refinement. Meshes after 0, 8, and 16 refinements. Figure. Test problem, local strategy adaptive refinement. Meshes after 0, 8, and 16 refinements O(1/N) η H1 semi norm L norm N: degrees of freedom Figure 3. Test problem, maximum strategy adaptive refinement. Plot of η, (u u h ) and u u h, versus the degrees of freedom N. References [1] R. Adams and J. Fournier, Sobolev Spaces, Academic press, 003. [] D. Arnold, F. Brezzi, B. Cockburn and D. Marini, Unified analysis of discontinuous Galerkin methods for elliptic problems, SIAM J. Numer. Anal., 39 (00),
19 UNIFIED A POSTERIORI ERROR ESTIMATE FOR STOKES EQUATIONS O(1/N) η H1 semi norm L norm 10 3 N: degrees of freedom 10 3 Figure 4. Test problem, local strategy adaptive refinement. Plot of η, (u u h ) and u u h, versus the degrees of freedom N. Figure 5. Test problem 3, maximum strategy adaptive refinement. Meshes after 0, 8, and 16 refinements. Figure 6. Test problem 3, local strategy adaptive refinement. Meshes after 0, 8, and 16 refinements. [3] M. Ainsworth and J.T. Oden, A Posteriori Error Estimators for Stokes and Oseen s E- quations. SIAM J. Numer. Anal., 17 (1997), [4] R. Bank and B. Welfert A Posteriori Error Estimates for the Stokes Problem, SIAM J. Numer. Anal. 8, (1991),
20 570 J. WANG, Y. WANG, AND X. YE [5] C. Carstensen, T. Gudi, M. Jensen A Unifying Theory of A Posteriori Control for Discontinuous Galerkin FEM, Numerische Mathematik, 11, (009), [6] P.G. Ciarlet, The finite element method for elliptic problems, North-Holland, Amsterdam, [7] M. Crouzeix and R.S. Falk, Nonconforming Finite Elements for the Stokes Problem, Math. Comp., 5, (1989), [8] M. Crouzeix and P.-A. Raviart, Conforming and nonconforming finite element methods for solving the stationary Stokes equations, I. RAIRO Anal. Numér. 7, (1973), [9] E. Dari, R. G. Durán and C. Padra Error estimators for nonconforming finite element approximations of the Stokes problem, Math. Comp. 64, (1995) [10] W. Doerfler and M. Ainsworth, Reliable a posteriori error control for non-conforming finite element approximation of Stokes flow, Math. Comp. 74 (005), [11] M. Fortin and M. Soulie, A non-conforming piecewise quadratic finite element on triangles, Int. J. Num. Meth. Eng. 19, (1983), [1] A. Hannukainen, R. Stenberg and M. Vohralik. A unified framework for a posteriori error estimation for the Stokes problem, Numerische Mathematik, submitted. [13] P. Hood and C. Taylor, A numerical solution of the Navier-Stokes equations using the finite element technique, Comp. and Fluids, 1, (1973), [14] P. Houston, D. Schötzau and T. P. Wihler Energy Norm A Posteriori Error Estimation for Mixed Discontinuous Galerkin Approximations of the Stokes Problem, J. Sci. Comput., -3, (005), [15] O. A. Karakashian and F. Pascal, A posteriori error estimates for a discontinuous Galerkin approximation of second order elliptic problems, SIAM J. Numer. Anal. 41 (003), [16] David Kay and David Silvester A posteriori error estimation for stabilised mixed approximations of the Stokes equations, SIAM J. Scientific Computing 1, (1999), [17] W. F. Mitchell Unified Multilevel Adaptive Finite Element Methods for Elliptic Problems, Ph.D. dissertation, University of Illinois at Urbana-Champaign, [18] F. Nobile, A posteriori error estimates for the finite element approximation of the Stokes problem, TICAM Report 03-13, April 003. [19] E. G. Sewell, Automatic generation of triangulations for piecewise polynomial approximation, Ph.D. dissertation, Purdue University, 197. [0] P.N. Shankar and M.D. Deshpande, Fluid Mechanics in the driven cavity, Annu. Rev. Fluid Mech., 3 (000), [1] R. Verfurth, A posteriori error estimators for the Stokes equations, Numer. Math., 3, (1989), [] R. Verfurth, A posteriori error estimators for the Stokes equations II. Non-conforming discretizations. Numer. Math. 60, (1991), [3] R. Verfurth, A review of a posteriori error estimation and adaptive mesh-refinement techniques. Teubner Skripten zur Numerik. B.G. Willey-Teubner, Stuttgart, [4] J. Wang, Y. Wang, and X. Ye, A posteriori error estimation for an interior panelty type method employing H(div) elements for the Stokes equations, SIAM J. Scientific Computing, accepted. [5] J. Wang and X. Ye, New finite element methods in computational fluid dynamics by H(div) elements, SIAM J. Numerical Analysis, 45 (007), Division of Mathematical Sciences, National Science Foundation, Arlington, VA 30, USA jwang@nsf.gov Department of Mathematics, Oklahoma State University, Stillwater, OK 74075, USA yqwang@math.okstate.edu Department of Mathematics, University of Arkansas at Little Rock, Little Rock, AR 704, USA xxye@ualr.edu
arxiv: 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 informationMaximum-norm a posteriori estimates for discontinuous Galerkin methods
Maximum-norm a posteriori estimates for discontinuous Galerkin methods Emmanuil Georgoulis Department of Mathematics, University of Leicester, UK Based on joint work with Alan Demlow (Kentucky, USA) DG
More informationBasic Concepts of Adaptive Finite Element Methods for Elliptic Boundary Value Problems
Basic Concepts of Adaptive Finite lement Methods for lliptic Boundary Value Problems Ronald H.W. Hoppe 1,2 1 Department of Mathematics, University of Houston 2 Institute of Mathematics, University of Augsburg
More informationA posteriori error estimator based on gradient recovery by averaging for discontinuous Galerkin methods
A posteriori error estimator based on gradient recovery by averaging for discontinuous Galerkin methods Emmanuel Creusé, Serge Nicaise October 7, 29 Abstract We consider some (anisotropic and piecewise
More informationAn a posteriori error estimator for the weak Galerkin least-squares finite-element method
An a posteriori error estimator for the weak Galerkin least-squares finite-element method James H. Adler a, Xiaozhe Hu a, Lin Mu b, Xiu Ye c a Department of Mathematics, Tufts University, Medford, MA 02155
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 informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University The Residual and Error of Finite Element Solutions Mixed BVP of Poisson Equation
More 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 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 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 informationINTRODUCTION TO FINITE ELEMENT METHODS
INTRODUCTION TO FINITE ELEMENT METHODS LONG CHEN Finite element methods are based on the variational formulation of partial differential equations which only need to compute the gradient of a function.
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 informationOverview. A Posteriori Error Estimates for the Biharmonic Equation. Variational Formulation and Discretization. The Biharmonic Equation
Overview A Posteriori rror stimates for the Biharmonic quation R Verfürth Fakultät für Mathematik Ruhr-Universität Bochum wwwruhr-uni-bochumde/num1 Milan / February 11th, 013 The Biharmonic quation Summary
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 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 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 informationBasic Principles of Weak Galerkin Finite Element Methods for PDEs
Basic Principles of Weak Galerkin Finite Element Methods for PDEs Junping Wang Computational Mathematics Division of Mathematical Sciences National Science Foundation Arlington, VA 22230 Polytopal Element
More 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 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 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 informationScientific Computing WS 2018/2019. Lecture 15. Jürgen Fuhrmann Lecture 15 Slide 1
Scientific Computing WS 2018/2019 Lecture 15 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 15 Slide 1 Lecture 15 Slide 2 Problems with strong formulation Writing the PDE with divergence and gradient
More informationLECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES. Sergey Korotov,
LECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES Sergey Korotov, Institute of Mathematics Helsinki University of Technology, Finland Academy of Finland 1 Main Problem in Mathematical
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 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 informationScientific Computing WS 2017/2018. Lecture 18. Jürgen Fuhrmann Lecture 18 Slide 1
Scientific Computing WS 2017/2018 Lecture 18 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 18 Slide 1 Lecture 18 Slide 2 Weak formulation of homogeneous Dirichlet problem Search u H0 1 (Ω) (here,
More informationA note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations
A note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations S. Hussain, F. Schieweck, S. Turek Abstract In this note, we extend our recent work for
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Nonconformity and the Consistency Error First Strang Lemma Abstract Error Estimate
More 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 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 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 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 informationSimple Examples on Rectangular Domains
84 Chapter 5 Simple Examples on Rectangular Domains In this chapter we consider simple elliptic boundary value problems in rectangular domains in R 2 or R 3 ; our prototype example is the Poisson equation
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 local-structure-preserving local discontinuous Galerkin method for the Laplace equation
A local-structure-preserving local discontinuous Galerkin method for the Laplace equation Fengyan Li 1 and Chi-Wang Shu 2 Abstract In this paper, we present a local-structure-preserving local discontinuous
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 informationNon-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions
Non-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions Wayne McGee and Padmanabhan Seshaiyer Texas Tech University, Mathematics and Statistics (padhu@math.ttu.edu) Summary. In
More 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 informationA Finite Element Method Using Singular Functions for Poisson Equations: Mixed Boundary Conditions
A Finite Element Method Using Singular Functions for Poisson Equations: Mixed Boundary Conditions Zhiqiang Cai Seokchan Kim Sangdong Kim Sooryun Kong Abstract In [7], we proposed a new finite element method
More informationMultigrid Methods for Elliptic Obstacle Problems on 2D Bisection Grids
Multigrid Methods for Elliptic Obstacle Problems on 2D Bisection Grids Long Chen 1, Ricardo H. Nochetto 2, and Chen-Song Zhang 3 1 Department of Mathematics, University of California at Irvine. chenlong@math.uci.edu
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 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 informationarxiv: v1 [math.na] 5 Jun 2018
PRIMAL-DUAL WEAK GALERKIN FINIE ELEMEN MEHODS FOR ELLIPIC CAUCHY PROBLEMS CHUNMEI WANG AND JUNPING WANG arxiv:1806.01583v1 [math.na] 5 Jun 2018 Abstract. he authors propose and analyze a well-posed numerical
More 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 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 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: 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 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 informationThe Mortar Boundary Element Method
The Mortar Boundary Element Method A Thesis submitted for the degree of Doctor of Philosophy by Martin Healey School of Information Systems, Computing and Mathematics Brunel University March 2010 Abstract
More informationA P4 BUBBLE ENRICHED P3 DIVERGENCE-FREE FINITE ELEMENT ON TRIANGULAR GRIDS
A P4 BUBBLE ENRICHED P3 DIVERGENCE-FREE FINITE ELEMENT ON TRIANGULAR GRIDS SHANGYOU ZHANG DEDICATED TO PROFESSOR PETER MONK ON THE OCCASION OF HIS 6TH BIRTHDAY Abstract. On triangular grids, the continuous
More informationConvergence and optimality of an adaptive FEM for controlling L 2 errors
Convergence and optimality of an adaptive FEM for controlling L 2 errors Alan Demlow (University of Kentucky) joint work with Rob Stevenson (University of Amsterdam) Partially supported by NSF DMS-0713770.
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 informationUltraconvergence of ZZ Patch Recovery at Mesh Symmetry Points
Ultraconvergence of ZZ Patch Recovery at Mesh Symmetry Points Zhimin Zhang and Runchang Lin Department of Mathematics, Wayne State University Abstract. The ultraconvergence property of the Zienkiewicz-Zhu
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 informationHigh order, finite volume method, flux conservation, finite element method
FLUX-CONSERVING FINITE ELEMENT METHODS SHANGYOU ZHANG, ZHIMIN ZHANG, AND QINGSONG ZOU Abstract. We analyze the flux conservation property of the finite element method. It is shown that the finite element
More 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 informationarxiv: v1 [math.na] 14 Apr 2011
A WEAK GALERKIN FINIE ELEMEN MEHOD FOR SECOND-ORDER ELLIPIC PROBLEMS JUNPING WANG AND XIU YE arxiv:1104.2897v1 [math.na] 14 Apr 2011 Abstract. In this paper, authors shall introduce a finite element method
More informationA COMPUTATIONAL STUDY OF THE WEAK GALERKIN METHOD FOR SECOND-ORDER ELLIPTIC EQUATIONS
A COMPUTATIONAL STUDY OF THE WEA GALERIN METHOD FOR SECOND-ORDER ELLIPTIC EQUATIONS LIN MU, JUNPING WANG, YANQIU WANG, AND XIU YE Abstract. The weak Galerkin finite element method is a novel numerical
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 informationMaximum norm estimates for energy-corrected finite element method
Maximum norm estimates for energy-corrected finite element method Piotr Swierczynski 1 and Barbara Wohlmuth 1 Technical University of Munich, Institute for Numerical Mathematics, piotr.swierczynski@ma.tum.de,
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 informationEfficient Augmented Lagrangian-type Preconditioning for the Oseen Problem using Grad-Div Stabilization
Efficient Augmented Lagrangian-type Preconditioning for the Oseen Problem using Grad-Div Stabilization Timo Heister, Texas A&M University 2013-02-28 SIAM CSE 2 Setting Stationary, incompressible flow problems
More information1. Introduction. We consider the model problem: seeking an unknown function u satisfying
A DISCONTINUOUS LEAST-SQUARES FINITE ELEMENT METHOD FOR SECOND ORDER ELLIPTIC EQUATIONS XIU YE AND SHANGYOU ZHANG Abstract In tis paper, a discontinuous least-squares (DLS) finite element metod is introduced
More informationDiscontinuous Galerkin Method for interface problem of coupling different order elliptic equations
Discontinuous Galerkin Method for interface problem of coupling different order elliptic equations Igor Mozolevski, Endre Süli Federal University of Santa Catarina, Brazil Oxford University Computing Laboratory,
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 informationHamburger Beiträge zur Angewandten Mathematik
Hamburger Beiträge zur Angewandten Mathematik Numerical analysis of a control and state constrained elliptic control problem with piecewise constant control approximations Klaus Deckelnick and Michael
More informationA POSTERIORI ERROR ESTIMATION FOR NON-CONFORMING QUADRILATERAL FINITE ELEMENTS
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 2, Number, Pages 8 c 2005 Institute for Scientific Computing and Information A POSTERIORI ERROR ESTIMATION FOR NON-CONFORMING QUADRILATERAL
More informationNew class of finite element methods: weak Galerkin methods
New class of finite element methods: weak Galerkin methods Xiu Ye University of Arkansas at Little Rock Second order elliptic equation Consider second order elliptic problem: a u = f, in Ω (1) u = 0, on
More informationGraded Mesh Refinement and Error Estimates of Higher Order for DGFE-solutions of Elliptic Boundary Value Problems in Polygons
Graded Mesh Refinement and Error Estimates of Higher Order for DGFE-solutions of Elliptic Boundary Value Problems in Polygons Anna-Margarete Sändig, Miloslav Feistauer University Stuttgart, IANS Journées
More informationWeak Galerkin Finite Element Methods and Applications
Weak Galerkin Finite Element Methods and Applications Lin Mu mul1@ornl.gov Computational and Applied Mathematics Computationa Science and Mathematics Division Oak Ridge National Laboratory Georgia Institute
More informationA Balancing Algorithm for Mortar Methods
A Balancing Algorithm for Mortar Methods Dan Stefanica Baruch College, City University of New York, NY 11, USA Dan Stefanica@baruch.cuny.edu Summary. The balancing methods are hybrid nonoverlapping Schwarz
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 informationNumerical Experiments Using MATLAB: Superconvergence of Nonconforming Finite Element Approximation for Second-Order Elliptic Problems
Applied Matematics, 06, 7, 74-8 ttp://wwwscirporg/journal/am ISSN Online: 5-7393 ISSN Print: 5-7385 Numerical Experiments Using MATLAB: Superconvergence of Nonconforming Finite Element Approximation for
More informationR T (u H )v + (2.1) J S (u H )v v V, T (2.2) (2.3) H S J S (u H ) 2 L 2 (S). S T
2 R.H. NOCHETTO 2. Lecture 2. Adaptivity I: Design and Convergence of AFEM tarting with a conforming mesh T H, the adaptive procedure AFEM consists of loops of the form OLVE ETIMATE MARK REFINE to produce
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 information20. A Dual-Primal FETI Method for solving Stokes/Navier-Stokes Equations
Fourteenth International Conference on Domain Decomposition Methods Editors: Ismael Herrera, David E. Keyes, Olof B. Widlund, Robert Yates c 23 DDM.org 2. A Dual-Primal FEI Method for solving Stokes/Navier-Stokes
More informationAn Iterative Substructuring Method for Mortar Nonconforming Discretization of a Fourth-Order Elliptic Problem in two dimensions
An Iterative Substructuring Method for Mortar Nonconforming Discretization of a Fourth-Order Elliptic Problem in two dimensions Leszek Marcinkowski Department of Mathematics, Warsaw University, Banacha
More informationDomain Decomposition Algorithms for an Indefinite Hypersingular Integral Equation in Three Dimensions
Domain Decomposition Algorithms for an Indefinite Hypersingular Integral Equation in Three Dimensions Ernst P. Stephan 1, Matthias Maischak 2, and Thanh Tran 3 1 Institut für Angewandte Mathematik, Leibniz
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 informationDISCRETE MAXIMUM PRINCIPLES in THE FINITE ELEMENT SIMULATIONS
DISCRETE MAXIMUM PRINCIPLES in THE FINITE ELEMENT SIMULATIONS Sergey Korotov BCAM Basque Center for Applied Mathematics http://www.bcamath.org 1 The presentation is based on my collaboration with several
More information