arxiv: v1 [math.na] 19 Dec 2017
|
|
- Cleopatra McCoy
- 5 years ago
- Views:
Transcription
1 Arnold-Winther Mixed Finite Elements for Stokes Eigenvalue Problems Joscha Gedicke Arbaz Khan arxiv: v [math.na] 9 Dec 207 Abstract This paper is devoted to study the Arnold-Winther mixed finite element method for two dimensional Stokes eigenvalue problems using the stress-velocity formulation. A priori error estimates for the eigenvalue and eigenfunction errors are presented. To improve the approximation for both eigenvalues and eigenfunctions, we propose a local post-processing. With the help of the local post-processing, we derive a reliable a posteriori error estimator which is shown to be empirically efficient. We confirm numerically the proven higher order convergence of the post-processed eigenvalues for convex domains with smooth eigenfunctions. On adaptively refined meshes we obtain numerically optimal higher orders of convergence of the post-processed eigenvalues even on nonconvex domains. Keywords a priori analysis, a posteriori analysis, Arnold-Winther finite element, mixed finite element, Stokes eigenvalue problem AMS subject classification 65N5, 65N25, 65N30 Introduction Over the last decade, the numerical analysis of the finite element method for eigenvalue problems has been of increasing interest because of various practical applications. Specifically, the numerical analysis of the Stokes eigenvalue problem is a broad research area. Huang et al. [6] discuss the numerical analysis of several stabilized finite element methods for the Stokes eigenvalue problem. In [2], Meddahi et al. proposed a finite element analysis of a pseudo-stress formulation for the Stokes eigenvalue problem. In [25], Türk et al. introduced a stabilized finite element method for two-field and three-field Stokes eigenvalue problems. From [22], one can also study a variety of mixed or hybrid finite element methods for eigenvalue problems. In the literature, most of the results based on the a posteriori error analysis for the finite element method (see [, 26] and the references therein) consider the source problem. In comparison there are only a few results for the a posteriori error analysis of the Stokes eigenvalue problem available. In [20], Lovadina et al. presented the numerical analysis for a residual-based a posteriori error estimator for the finite element discretization of the Stokes eigenvalue problem. In [9], Liu et al. proposed the finite element approximation of the Stokes eigenvalue problem based on a projection method, and derive some superconvergence results and the related recovery type a posteriori error estimators. A posteriori error estimators for stabilized low-order mixed The first author has been funded by the Austrian Science Fund (FWF) through the project P 2997-N32. The second author has been funded by the Mathematics Center Heidelberg (Match) at University of Heidelberg, Germany. Faculty of Mathematics, University of Vienna, 090 Vienna, Austria (joscha.gedicke@univie.ac.at). Corresponding author: School of Mathematics, University of Manchester, M3 9PL Manchester, UK (arbaz.khan@manchester.ac.uk).
2 finite elements for the Stokes eigenvalue problem are presented by Armentano et al. [2]. A new adaptive mixed finite element method based on residual type a posteriori error estimators for the Stokes eigenvalue problem is proposed by Han et al. [4]. In [5], Huang presented two stabilized finite element methods for the Stokes eigenvalue problem based on the lowest equal-order finite element pair and also discussed a posteriori lower and upper bounds of Stokes eigenvalues. Arnold and Winther introduced the strongly symmetric Arnold-Winther mixed finite element (MFEM) for linear elastic problems in [3] and proved its stability for any material parameters. Hence, the proposed Arnold-Winther MFEM is also stable for the Stokes problem as a limit case of linear elasticity. In [], Carstensen et al. presented the Arnold-Winther mixed finite element formulation for the Stokes source problem. In this paper, we are presenting the Arnold-Winther mixed finite element formulation for the two dimensional Stokes eigenvalue problem using the stress-velocity formulation. In principal, the stress-velocity formulation is originated from a physical model where incompressible Newtonian flows are modeled by the conservation of momentum and the constitute law. The stress-velocity formulation for the Stokes eigenvalue problem reads as follows: find a symmetric stress tensor σ, a nonzero eigenfunction u, and an eigenvalue λ such that div σ = λu in Ω, Aσ ɛ(u) = 0 in Ω, u = 0 on Ω, for a bounded Lipschitz domain Ω R 2. Here σ, Aσ and ɛ(u) are the stress tensor, the deviatoric stress tensor and the deformation rate tensor, respectively. For simplicity, we restrict ourselves in this paper to the analysis for simple eigenvalues λ. This simplifies the a priori error analysis in that the eigenspace to λ is one-dimensional. Note that the presented a posteriori error estimator can directly be extended to multiple eigenvalues, in the way that the sum over the squared error estimators, for all the discrete eigenvectors that approximate eigenvectors to the same multiple eigenvalue, provides an upper bound of the eigenvalue error. The remaining parts of this paper are organized as follows: Section 2 presents the necessary notation, the formulation of the problem, the discretization of the domain and the mixed finite element formulation. Section 3 is devoted to study the a priori error analysis of the eigenvalue problem. Section 4 presents the local post-processing and its higher order convergence analysis. The a posteriori error analysis is introduced in Section 5. Finally, we verify the reliability and efficiency of the a posteriori error estimator and the higher order convergence of post-processed eigenvalues in three numerical experiments in Section 6. 2 Preliminaries Let H s (ω) be the standard Sobolev space with the associated norm s,ω for s 0. In case of ω = Ω, we use s instead of s,ω. Let H s (ω) := (H s (ω)) be the dual space of H s (ω). Now we extend the definitions for vector and matrix-valued function. Let H s (ω) = H s (ω; R 2 ) and H s (ω, R 2 2 ) be the Sobolev spaces over the set of 2-dimensional vector and 2 2 matrix-valued function, respectively. Define v = (v, v 2 ) t R 2, τ = (τ ij ) 2 2 and σ = (σ ij ) 2 2 R 2 2, then v := τ v := ( v x v 2 x v y v 2 y ), div(τ) := ( τ x τ 2 x ( ) τ v + τ 2 v 2, τ : σ := τ 2 v + τ 22 v 2 i,j + τ2 y + τ22 y τ ij σ ij, ), tr τ := τ + τ 22, δ := 2 2 unit matrix. Let the divergence conforming stress space H(div, Ω, R 2 2 ) be defined as H(div, Ω, R 2 2 ) := {τ L 2 (Ω; R 2 2 ) div τ L 2 (Ω)}, 2
3 equipped with the norm τ 2 H(div,Ω,R 2 2 ) := (τ, τ ) + (div τ, div τ ). In this paper, we are using the notation (, ) ω for the L 2 (ω; R 2 2 ) inner product τ : τ dx as ω well as the L 2 (ω) inner product ω τ τ dx. For w = Ω, the notation (, ) is used instead of (, ) Ω. The symbols and are used throughout the paper to denote inequalities which are valid up to positive constants that are independent of the local mesh size h but may depend on the size of the eigenvalue λ and the coefficient ν. Let Ω R 2 be a bounded and connected Lipschitz domain. Consider the Stokes eigenvalue problem ν u + p = λu in Ω, div u = 0 in Ω, u = 0 on Ω, with the compatibility condition on the pressure p dx = 0. (2) Define σ = (σ i,j ) 2 2 as a stress tensor and the deformation rate tensor as Ω ɛ(u) := 2 ( u + ( u)t ). From () we can derive the stress-velocity-pressure formulation for the Stokes eigenvalue problem, which is the set of original physical equations for incompressible Newtonian flow, div σ = λu in Ω, σ + pδ 2νɛ(u) = 0 in Ω, div u = 0 in Ω, u = 0 on Ω. Next, we define the deviatoric operator A : S S, where S is the space of symmetric tensors S = {τ R 2 2 τ = τ t }. Then we can define the deviator Aτ of τ by Aτ := 2 ( τ 2 (tr τ )δ) for all τ S. Here, Ker(A) = {qδ S q R} and Aτ is a trace-free tensor. In addition, A satisfies for all τ, σ S the following properties, (Aτ, σ) = (τ, Aσ), (Aτ, 2νAσ) = (Aσ, τ ) = ((σ, τ ) ) 2ν 2 (tr σ, tr τ ), Aτ τ. 2 Using the deviatoric tensor A, we arrive at the stress-velocity formulation of the Stokes eigenvalue problem div σ = λu in Ω, Aσ ɛ(u) = 0 in Ω, u = 0 on Ω. () (3) 3
4 Using the compatibility condition (2), we have tr σ dx = 0. Defining V := L 2 (Ω) and Φ := H(div, Ω; S)/R Ω { } τ H(div, Ω, S) tr τ dx = 0, Ω we can now derive the following weak form for problem (3): find σ Φ, u V with u 0 =, and λ R + such that (Aσ, τ ) + (div τ, u) = 0 for all τ Φ, (div σ, v) = λ(u, v) for all v V. (4) Let {T h } denote a family of regular triangulations of Ω into triangles K of diameter h K. For each T h, we define E h as the set of all edges of T h and h E as the length of the edge E E h. Furthermore, let [w] denote the jump of w, [w] E := (w K+ ) E (w K ) E if E = K + K. Let Th and ɛ Th denote the piecewise gradient and the piecewise symmetric gradient, i.e. ( Th ) K = ( K ) and (ɛ Th ) K = ɛ( K ), for all K T h. Finally, we define the following finite element spaces associated with the triangulation T h, AW k (K) := {τ P k+2 (K; S) div τ P k (K; R 2 )}, Φ h := {τ Φ τ K AW k (K)}, V h := {v L 2 (Ω) v K P k (K; R 2 )}. Here, AW k (K) denotes the Arnold-Winther MFEM of index k of [3] and P k (K) denotes the set of all polynomials of total degree up to k on the domain K. The space Φ h consists of all symmetric polynomial matrix fields of degree at most k + together with the divergence free matrix fields of degree k + 2. Note that Φ h Φ. When τ h Φ h we have that τ h has continuous normal components and τ h satisfies Ω tr τ hdx = 0. The Arnold-Winther MFEM of the Stokes eigenvalue problem becomes the following: find σ h Φ h, u h V h with u h 0 =, and λ h R + such that (Aσ h, τ h ) + (div τ h, u h ) = 0 for all τ h Φ h, (div σ h, v h ) = λ h (u h, v h ) for all v h V h. (5) 3 A priori error analysis Our main aim is to show that the solutions of the Arnold-Winther MFEM of the Stokes eigenvalue problem converge to the solution of the corresponding spectral problem which comes to apply the classical spectral approximation theory presented in [4] to the associated source problem. Consider the stress-velocity formulation for the associated source problem div σ f = f in Ω, Aσ f ɛ(u f ) = 0 in Ω, u f = 0 on Ω. (6) 4
5 Using the compatibility condition (2), we have tr σ f dx = 0. Ω For the above problem (6) there exist well-known regularity results for convex domains with sufficiently smooth boundary Ω. If f L 2 (Ω) then the solution of problem (6) satisfies u H 2 (Ω) H ( Ω), p H (Ω)/R, σ H (Ω; S), and u 2 + p + σ f 0. (7) Using the well-posedness of problem (6), the operators S : V Φ and T : V V are well defined for any f V such that Sf = σ f, and T f = u f are the stress and velocity solutions, respectively. We now define the following weak form for problem (3): for given f V, (Sf, T f) (Φ, V ) is the solution of (A(Sf), τ ) + (div τ, T f) = 0 for all τ Φ, (div(sf), v) = (f, v) for all v V. The above problem (8) has unique solution from the well known inf-sup condition of the mixed formulation and [, Lemma 2.]. For the discrete solution operators S h : V Φ h and T h : V V h, the Arnold-Winther MFEM of the Stokes source problem becomes the following: find S h f Φ h and T h f V h such that (A(S h f), τ h ) + (div τ h, T h f) = 0 for all τ h Φ h, (div(s h f), v h ) = (f, v h ) for all v h V h. From [], the discrete source problem is well-posed and has a unique solution. From [3] we have the following a priori estimates for (σ f, u f ) (Φ H k+2 (Ω; S)) H k+2 (Ω) and f H k (Ω) Sf S h f 0 h m σ f m, m k + 2, (0) div(sf S h f) 0 h m div σ f m, 0 m k +, () T f T h f 0 h m u f m+, m k +. (2) Hence, we can state the following convergence results by (0) and (2) T T h L(V,V ) 0 if h 0, (3) S S h L(V,Φ) 0 if h 0. (4) The above results (3) and (4) are equivalent to the convergence of eigenvalues and eigenfunctions. Thus, using the abstract theory from [6, 22] and the a priori results (0) and (2), we have u u h 0 h m, m k +, (5) σ σ h 0 h m, m k + 2. (6) Note that p = tr σ/2, hence the approximation of the pressure is defined by p h = tr σ h /2 and satisfies the following estimate p p h 0 = 2 tr σ tr σ h 0 σ σ h 0 h m, m k + 2. The next lemma establishes a connection between the errors in the eigenvalues and in the eigenfunctions. (8) (9) 5
6 Lemma 3.. Let (σ, u, λ) and (σ h, u h, λ h ) be solutions of the continuous eigenvalue problem (4) and the discrete eigenvalue problem (5), respectively. Then, we have the identity λ λ h = 2ν A(σ σ h ) 2 0 λ h u u h 2 0. (7) Proof. From (4) and (5), using (u, u) = and (u h, u h ) =, we have Moreover, it follows Then we obtain (2νAσ, Aσ) = λ(u, u) = λ, (2νAσ h, Aσ h ) = λ h (u h, u h ) = λ h, (2νAσ, Aσ h ) = λ h (u h, u). 2ν A(σ σ h ) 2 0 = 2ν(A(σ σ h ), A(σ σ h )) Using u u h 2 0 = 2 2(u h, u) yields = 2ν(Aσ, Aσ) + 2ν(Aσ h, Aσ h ) 4ν(Aσ, Aσ h ) = λ + λ h 2λ h (u h, u). λ λ h = 2ν A(σ σ h ) 2 0 2λ h + 2λ h (u h, u). λ λ h = 2ν A(σ σ h ) 2 0 λ h u u h 2 0. An immediate consequence of Lemma 3. are the following two lemmas. Lemma 3.2. For sufficiently smooth u H k+2 (Ω), σ H k+2 (Ω; S), the following a priori error estimate holds Proof. The assertion follows from (5) (7). λ λ h h 2m, 0 m k +. (8) Lemma 3.3. For sufficiently smooth u H k+2 (Ω), σ H k+2 (Ω; S), the following a priori error estimate holds div(σ σ h ) 0 h m, 0 m k +. (9) Proof. From (4) and (5), we have div(σ) = λu and div(σ h ) = λ h u h. Therefore we get div(σ σ h ) 0 = λ h u h λu 0 (λ h λ) u h 0 + λ u h u 0. The assertion follows from (5) (7), and u h 0 =. Let P h denote the L 2 projection onto V h with the well known approximation property P h v v 0 h m v m, for all v H m (Ω). (20) In the following we relate the discrete eigenfunctions (σ h, u h ) to the discrete approximations (σ λ,h, u λ,h ) of the associated source problem (9) with right hand side f = λu V. Then σ λ,h = λs h u and u λ,h = λt h u and we have the following lemma. 6
7 Lemma 3.4 ([, Theorem 3.4]). For sufficiently smooth boundary Ω, σ H k+2 (Ω; S), and div σ H k+ (Ω), it holds that P h u u λ,h 0 h k+3 ( σ k+2 + div σ k+ ). (2) Next, we prove a similar estimate for the eigenfunction u h which is used in Section 4 to derive an error estimate for the post-processed eigenfunction. Theorem 3.5. For sufficiently smooth boundary Ω, σ H k+2 (Ω; S), div σ H k+ (Ω), it holds that P h u u h 0 h k+3 (h k + σ k+2 + div σ k+ ). (22) The proof of Theorem 3.5 does not follow directly from Lemma 3.4, because the eigenvalue problem does not satisfy an orthogonality condition, i.e. (div(σ σ h ), v h ) 0, therefore we need the following lemma. Lemma 3.6. For sufficiently smooth boundary Ω, σ H k+2 (Ω; S), and div σ H k+ (Ω), the difference between the discrete eigenfunction u h and the discrete solution of the associated source problem u λ,h can be estimated by Proof. First we define some relations u λ,h u h h k+3 (h k + σ k+2 + div σ k+ ). λt u = u, λ h T h u h = u h and λt h u = u λ,h. Since λ is a simple eigenvalue, the eigenspace of the eigenvalue λ is spanned by u. The operator T : V V is self adjoint. Therefore the orthogonal complement of u is a T -invariant subspace denoted by U,V. Moreover the eigenvalue λ does not belong to the spectrum of T U,V which is defined as T U,V : U,V U,V. Thus, we can define the following invertible map (I λt ) : U,V U,V. Here (I λt ) is bounded. Define δ h := u λ,h u h (u λ,h u h, u)u. Using u 0 =, it follows Hence we have (δ h, u) = (u λ,h u h (u λ,h u h, u)u, u) = 0. δ h U,V and δ h 0 (I λt )δ h 0. Since (I λt )u = 0, the following estimate holds Moreover, we obtain the following inequality δ h 0 (I λt )(u λ,h u h ) 0. (23) (I λt )(u λ,h u h ) = (λ h T h λt )(u λ,h u h ) + u λ,h u h λ h T h (u λ,h u) λ h T h (u u h ) = (λ h T h λt )(u λ,h u h ) + (λ λ h )T h u λ h T h (u λ,h u). (24) 7
8 Inserting the result of (24) into (23) and using the triangle inequality, implies where δ h 0 (λ h T h λt h )(u λ,h u h ) 0 + (λt h λt )(u λ,h u h ) 0 + (λ λ h )T h u 0 + λ h T h (u λ,h u) 0, I + I 2 + I 3 + I 4, I = (λ h T h λt h )(u λ,h u h ) 0, I 2 = (λt h λt )(u λ,h u h ) 0, I 3 = (λ λ h )T h u 0, I 4 = λ h T h (u λ,h u) 0. Using the boundedness of T h, the estimates of I and I 3 read I λ λ h u λ,h u h 0 and I 3 λ λ h. (25) Using f = (u λ,h u h ) in (2) and inserting the regularity result (7), we obtain I 2 = (λt h λt )(u λ,h u h ) 0 λh (u λ,h u h ) 0. (26) Adding and subtracting T h (P h u) in I 4 leads to I 4 λ h T h (u λ,h P h u) 0 + λ h T h (P h u u) 0. (27) Applying Lemma 3.4 in first term of the estimate (27), implies T h (u λ,h P h u) 0 h k+3 ( σ k+2 + div σ k+ ). Note that the second term in the estimate (27) is equal to zero due to definition of T h, since the right hand side of problem (9) vanishes for f = P h u u. Hence Using the definition of δ h, leads to I 4 λ h h k+3 ( σ k+2 + div σ k+ ). u λ,h u h 0 δ h 0 + (u λ,h u h, u)u 0. (28) The first term of the estimate (28) is already estimated. To estimate the second term, observe that with u 0 =, we get Moreover, we have Since u 0 = and u h 0 =, we get (u λ,h u h, u)u 0 = (u λ,h u h, u). (u λ,h u h, u)u 0 (u λ,h u, u) + (u u h, u). (u u h, u) = (u h, u) = 2 u u h 2 0 h 2k+2. (29) Next, we estimate (u λ,h u, u). Choosing the test functions τ = σ λ,h σ and v = u λ,h u in (4), gives (Aσ, σ λ,h σ) + (div(σ λ,h σ), u) = 0, (div σ, u λ,h u) = λ(u, u λ,h u). 8
9 From that, we obtain Moreover, we have λ(u λ,h u, u) = (div σ, u λ,h u) + (Aσ, σ λ,h σ) + (div(σ λ,h σ), u). λ(u λ,h u, u) = (div(σ σ λ,h ), u λ,h u) + (div σ λ,h, u λ,h u) + (A(σ σ λ,h ), σ λ,h σ) + (Aσ λ,h, σ λ,h σ) + (div(σ λ,h σ), u u λ,h ) + (div(σ λ,h σ), u λ,h ). (30) Choosing τ h = σ λ,h, v h = u λ,h, f = λu in (9) and τ = σ λ,h, v = u λ,h in (4), and subtracting (4) from (9), we obtain (A(σ λ,h σ), σ λ,h ) + (div σ λ,h, u λ,h u) = 0, (div(σ λ,h σ), u λ,h ) = 0. (3) Inserting the result from (3) in (30), implies λ(u λ,h u, u) = (A(σ σ λ,h ), σ λ,h σ) + 2(div(σ σ λ,h ), u λ,h u) σ σ λ,h div(σ σ λ,h ) u u λ,h 2 0 h 2k+4 σ 2 k+2 + h 2k+2 div σ 2 k+ + h 2k+2 u 2 k+2 h 2k+2. (32) Summing up the estimates (25) (27) and (29), (32), we obtain u λ,h u h 0 λ λ h u λ,h u h 0 + h u λ,h u h 0 + λ λ h + h k+3 ( σ k+2 + div σ k+ ) + h 2k+2. Choosing h small enough and using λ λ h h 2k+2, the following estimate holds u λ,h u h 0 h k+3 (h k + σ k+2 + div σ k+ ). Proof of Theorem 3.5. Combining Lemma 3.4 and 3.6, implies the desired result. 4 Local post-processing In this section we present an improved approximation of the eigenfunction by a local postprocessing similar to the one in [] for the source problem. Define for m k + 2 W h = {v L 2 (Ω) v K P m (K; R 2 ) for all K T h }. Choose u h W h on each K T h with P K = P h K as a solution of the system P K u h = u h, (ɛ(u h), ɛ(v)) K = (Aσ h, ɛ(v)) K for all v (δ P K )W h K. (33) Here, u h is the Riesz representation of the linear functional (Aσ h, ɛ( )) K in the Hilbert space (δ P K )W h K {v m P m (K; R 2 ) (v m, w k ) K = 0 for all {w k P k (K; R 2 )}, 9
10 which is equipped with the scalar product (ɛ( ), ɛ( )) K. We define the local post-processing on each triangle with Lagrange multiplier µ k P k (K; R 2 ) as follows (ɛ(u h), ɛ(v m )) K + (µ k, v m ) K = (Aσ h, ɛ(v m )) K for all v m P m (K; R 2 ), (u h, w k ) K = (u h, w k ) K for all w k P k (K; R 2 ). (34) From Korn s inequality, we get positive definiteness of (ɛ( ), ɛ( )) K on (δ P K )W h K, and since P k (K; R 2 ) P m (K; R 2 ), we obtain (v m, µ k ) K sup C µ k 0,K for all µ k P k (K; R 2 ). 0 v m P m(k;r 2 ) v m,k Hence, there exist unique solutions on each triangle, c.f. [8]. Combining the identity Aσ = ɛ(u) and (33), gives the following error identity (ɛ(u u h), ɛ(v)) K = (A(σ σ h ), ɛ(v)) K for all v (δ P K )W h K. Theorem 4.. With sufficiently smooth boundary Ω, if u H m+ (Ω), σ H k+2 (Ω; S), and div σ H k+ (Ω) solve problem (), then the following estimates hold for the post-processed eigenfunction u h W h u u h 0 h k+3 (h k + σ k+2 + div σ k+ ) + h m+ u m+, Th (u u h) 0 h k+2 (h k + σ k+2 + div σ k+ ) + h m u m+. Proof. Let û denote the L 2 projection of u onto Wh. From the triangle inequality we get u u h 0 u û 0 + P h (û u h) 0 + (δ P h )(û u h) 0. (35) Applying (20) in first part of the right hand side of (35) leads to Using P K u h = u h on each K T h shows u û 0 h m+ u m+ for all u H m+ (Ω). (36) P K (û u h) 0,K = P K û u h 0,K. Here V h Wh, so we can apply the result of Theorem 3.5 to deduce P h (û u h) 0 = P h u u h 0 h k+3 (h k + σ k+2 + div σ k+ ). (37) The estimates for the last term in (35) are identical to the estimates for the post-processing for the source problem, hence we omit the details and quote the final estimate from the proof of [, Theorem 4.] (δ P h )(û u h) 0 h k+3 ( σ k+2 + div σ k+ ) + h m+ u m+. (38) Combining (35) with (36) (38) proves the first estimate. To prove the second estimate, we use the triangle inequality u u h,k u û,k + P h (û u h),k + (δ P h )(û u h),k. Applying an discrete inverse inequality leads to u u h,k u û,k + h P h (û u h) 0,K + h (δ P h )(û u h) 0,K. 0
11 Combining (37) and (38), with the estimate u û,k h m u m+,k obtained from interpolation type arguments and inverse estimates as in the proof of [, Theorem 4.], proves the second estimate. Definition 4.2. For (u h, u h ) 0, we define the post-processed eigenvalue as the value of the Rayleigh quotient of the post-processed eigenfunction λ h := (div σ h, u h ) (u h, u h ). (39) Theorem 4.3. Let (σ, u, λ) (Φ H k+2 (Ω; S)) H k+2 (Ω) R + be the solution of (4), with u 0 =. For sufficiently small h, the following a priori estimate for the post-processed eigenvalue λ h holds λ λ h h 2k+4. Proof. From (4), (5), (34), (39), and (u, u) =, we have Moreover, it follows Therefore we have (2νAσ, Aσ) = λ(u, u) = λ, and (2νAσ h, Aσ h ) = λ h(u h, u h). 2ν A(σ σ h ) 2 0 = 2ν(A(σ σ h ), A(σ σ h )) = 2ν(Aσ, Aσ) + 2ν(Aσ h, Aσ h ) 4ν(Aσ, Aσ h ) = λ + λ h(u h, u h) 2λ h(u, u h) + 2(div σ h, u) + 2λ h(u, u h). λ λ h = 2ν A(σ σ h ) 2 0 λ h λ h(u h, u h) + 2λ h(u, u h) 2(div σ h + λ hu h, u). Using (u, u) + (u h, u h ) 2(u, u h ) = u u h 2 0, and (div σ h + λ h u h, u h ) = 0 we get which leads to λ λ h =2ν A(σ σ h ) 2 0 λ h u u h 2 0 2(div σ h + λ hu h, u u h), (40) λ λ h =2ν A(σ σ h ) 2 0 λ h u u h 2 0 2(div σ h div σ, u u h) 2(λ hu h λu, u u h) =2ν A(σ σ h ) 2 0 λ h u u h 2 0 2(div σ h div σ, u u h) 2(λ h(u h u), u u h) 2((λ h λ)u, u u h) =2ν A(σ σ h ) λ h u u h (div(σ σ h ), u u h) + 2(λ λ h)(u, u u h). From Cauchy-Schwarz inequality it follows This results in λ λ h 2ν A(σ σ h ) λ h u u h div(σ σ h ) 0 u u h λ λ h u u h 0. λ λ h σ σ h u u h div(σ σ h ) 0 u u h 0 + λ λ h 2. Using the estimates (6) and (9), and Theorem 4., we obtain for h small enough λ λ h h 2k+4.
12 5 A posteriori error analysis First, we define ũ h H = H 0 (Ω). Here ũ h is closely related to the discontinuous approximation u h. The primal mixed formulation of the Stokes source problem with right hand side λu reads as follows: find σ L and u H such that (σ, ɛ(v)) = λ(u, v) for all v H, (2νAσ, Aτ ) (τ, ɛ(u)) = 0 for all τ L. (4) Following [0, Section.3], we have the following lemma Lemma 5.. The operator A : X X, defined for (σ, u) X := L H by is linear, bounded, and bijective. From the above lemma, we obtain (A(σ, u))(τ, v) := (2νAσ, Aτ ) (σ, ɛ(v)) (τ, ɛ(u)), A(σ σ h ) 0 + ɛ(u ũ h ) 0 Res L L + Res H H, (42) for any approximation (σ h, ũ h ) L H of the primal source problem with right hand side λu, where Res H (v) := λ(u, v) + (σ h, ɛ(v)) for all v H, Res L (τ ) := (2νAσ h, Aτ ) (τ, ɛ(ũ h )) for all τ L. Lemma 5.2. Let (σ, u, λ) L H R + be a solution of (4). Then the approximation (σ h, ũ h ) L H satisfies with the higher order term A(σ σ h ) 0 + ɛ Th (u ũ h ) 0 ( ) /2 Aσ h ɛ(ũ h ) K,0 + h 2 K λu + div σ h Θ, K T h Proof. Gauss theorem implies for any v H Θ := λ h u u h 0 + λ λ h. Res H (v) = (λu + div σ h, v). Let v h denote the Scott-Zhang interpolation [23] of v, then it holds that Res H (v) = (λu + div σ h, v v h ) + (λu + div σ h, v h ). For the second term on the right hand side, (5) and (33) show (λu + div σ h, v h ) = (λu λ h u h, v h ) = (λu λ h u h, v h ) = (λ λ h )(u, v h ) + λ h (u u h, v h ). 2
13 Applying Cauchy-Schwarz inequality, it follows Res H (v) ( K T h h 2 K λu + div σ h 2 K,0 + (λ h u u h + λ λ h ) v h 0. ) /2 ( K T h h 2 K (v v h) 2 K,0 Poincare s inequality, and stability and approximation properties of the Scott-Zhang interpolation show ( ) /2 Res H (v) h 2 K λu + div σ h 2 K,0 + Θ v. K T h Finally, Res L L is estimated as Res L (τ ) = (Aσ h ɛ(ũ h )) : τ dx Aσ h ɛ(ũ h ) 0 τ 0. Ω ) /2 Now, we present an a posteriori error estimator that involves the discontinuous post-processed approximation u h H (T h ) := {v L 2 (Ω) v T H (T ) for all T T h }. In the following let ũ h H be the conforming approximation to u, that is obtained from the discontinuous post-processed function u h by taking the arithmetic mean value ũ h (z) := {K T h : z K} K T h :z K u h(z) K for each vertex and edge degree of freedom in z R 2. A discrete scaling argument and [7, Theorem 2.2] show K T h h 2 K u h ũ h 2 K,0 + K T h ɛ Th (u h) ɛ(ũ h ) 2 K,0 E E h h E [u h] 2 E,0. (43) Theorem 5.3. Let (σ, u, λ) L H R + be a solution of (4). The post-processed eigenfunction u h H (T h ) satisfies the following reliability estimate A(σ σ h ) 0 + ɛ Th (u u h) 0 η + Υ, (44) for the a posteriori error estimator η 2 := Aσ h ɛ Th (u h) h 2 K λ hu h + div σ h 2 K,0 + h E [u h] 2 E,0, K T h E E h and the higher order term Υ := ( T T h h 2 K λu λ hu h 2 K,0 3 ) /2 + Θ.
14 Proof. Adding and subtracting ɛ(ũ h ) in the second term of (44) and using the triangle inequality, we have A(σ σ h ) 0 + ɛ Th (u u h) 0 A(σ σ h ) 0 + ɛ Th (u ũ h ) 0 + ɛ Th (u h) ɛ(ũ h ) 0. Applying Theorem 5.2, implies A(σ σ h ) 0 + ɛ Th (u u h) 0 Aσ h ɛ(ũ h ) 0 + ɛ Th (u h) ɛ(ũ h ) 0 ( ) /2 + h 2 K λu + div σ h 2 K,0 + Θ. K T h Another triangle inequality yields A(σ σ h ) 0 + ɛ Th (u u h) 0 Aσ h ɛ Th (u h) 0 + ɛ Th (u h) ɛ(ũ h ) 0 ( ) /2 + h 2 K λ hu h + div σ h 2 K,0 + Υ. K T h Using (43), the desired estimate holds. Theorem 5.4. The a posteriori error estimator η 2 provides an upper bound of the post-processed eigenvalue error for sufficiently small mesh size h, up to the higher order term Υ 2, λ λ h η 2 + Υ 2. Proof. Adding and subtracting ũ h in the last term of (40) yields λ λ h = 2ν A(σ σ h ) 2 0 λ h u u h 2 0 2(div σ h + λ hu h, u ũ h ) Applying Gauss theorem, implies 2(div σ h + λ hu h, ũ h u h) = 2ν A(σ σ h ) 2 0 λ h u u h 2 0 2(divσ h divσ, u ũ h ) 2(λ hu h λu, u ũ h ) 2(div σ h + λ hu h, ũ h u h). λ λ h = 2ν A(σ σ h ) λ h u u h 2 0 2(σ σ h, (u ũ h )) + 2λ h(u u h, u ũ h ) + 2(λ λ h)(u, u ũ h ) 2(divσ h + λ hu h, ũ h u h). Using Cauchy-Schwarz inequality, we have λ λ h A(σ σ h ) u u h σ σ h (u ũ h ) u ũ h (λ λ h) 2 + h 2 K λ hu h + div σ h 2 0,K + h 2 K ũ h u h 0,K. K T h K T h From the trace estimate [9, Proposition 3., IV.3], we get σ σ h 0 A(σ σ h ) 0 + div(σ σ h ). 4
15 Note that div(σ σ h ) = sup (λu + div σ h, v) = Res H H. v H v = Therefore, (42), Korn s and triangle inequalities lead to λ λ h A(σ σ h ) ɛ Th (u u h) h 2 K λ hu h + div σ h 2 0,K K T h + u h ũ h ɛ Th (u h ũ h ) h 2 K u h ũ h 2 0,K K T h + u u h (λ λ h) 2. For h small enough such that λ λ h < /2, we conclude with (43) and (44) λ λ h η 2 + Υ 2. 6 Numerical experiments In this section, we present numerical results for the square domain, the L-shaped domain and the slit domain. We verify the proven (asymptotic) reliability of the a posteriori error estimator of Section 5 and show empirically its efficiency. We present numerical results that show sixth order convergence of the post-processed eigenvalues of Section 4 on adaptively refined meshes even for the (nonconvex) L-shaped and slit domains. In all experiments, we take to the lowest order Arnold-Winther finite element (k = ), the parameter ν = and the polynomial order P 3 (T h, R 2 ) for the post-processing. Since the exact eigenvalues for all three domains are unknown, we compare the computed eigenvalues to some reference values with high accuracy. Note that the eigenvalues of the Stokes eigenvalue problem are related to the eigenvalues of the buckling eigenvalue problem of clamped plates via the stream function formulation. Hence, we can use known [5] or computed reference eigenvalues for the plate eigenvalue problem. We consider the standard adaptive finite element loop Solve Estimate Mark Refine, that creates a sequence of adaptively refined (nested) regular meshes (T l ) level index l. For the algebraic eigenvalue solver we use the Matlab implementation of ARPACK [8]. To estimate the error, we compute the a posteriori error estimator of Section 5 η 2 l = K T l Aσ l ɛ Tl (u l ) 2 0,K + E E l h E [u l ] 2 0,E + K T l h 2 K λ l u l + div σ l 2 0,K. In the above formula, u l W l is the solution of the local post-processing which is given in Section 4. Since the conforming function ũ l of Section 5 is closely related to u l we also compare the error estimator ηl 2 to the heuristical estimator µ 2 l = Aσ l ɛ(ũ l ) 2 0,K + h 2 K λ l ũ l + div σ l 2 0,K, K T l K T l 5
16 η 2 l (uniform) λ λl (uniform) λ λ l (uniform) η 2 l (adaptive) λ λl (adaptive) λ λ l (adaptive) order 2 order µ 2 l (uniform) λ λl (uniform) λ λl (uniform) µ 2 l (adaptive) λ λl (adaptive) λ λl (adaptive) order 2 order Nl Nl (a) (b) Figure : Convergence history of (a) λ λ l, λ λ l, η2 l, (b) λ λ l, λ λ l and µ 2 l on uniformly and adaptively refined meshes for the square domain. where we replaced u l by ũ l, and λ l is computed from (39) with u l replaced by ũ l. We numerically demonstrate reliability and efficiency of µ 2 l for the eigenvalue error λ λ l. We mark triangles of the triangulation T l in a minimal set of marked triangles M l according to the bulk marking strategy [3], such that θηl 2 η2 l (M l) for the bulk parameter θ = /2, and refine the mesh with the red-green-blue refinement strategy [26]. Let N l denote the degrees of freedom N l := dim(φ l ) + dim(v l ). Note that for uniform meshes, we have the relationship O(N r l 6. Square domain ) O(h 2r l ), r > 0. In the first example, we consider the square domain Ω = (0, ) 2. The reference value for the first eigenvalue λ = is taken from [5, 2]. Figures (a) and (b) are devoted to the convergence history of the eigenvalue errors and a posteriori error estimators. Due to the smoothness of the eigenfunction, the error of the post-processed eigenvalue λ l and the corresponding a posteriori error estimator ηl 2 achieve optimal third order of convergence for both uniform and adaptive meshes. Moreover, the error for the eigenvalue λ l and the estimator µ 2 l also achieve optimal convergence of O(N 3 ) for both uniform and adaptive meshes. For both uniform and adaptive meshes, the convergence rate of the eigenvalue error of λ l is equal to O(N 2 ), which confirms the theoretical result (8). The streamline plot of the discrete eigenfunction u l and a plot of the discrete pressure p l = trσ l /2 are displayed in 2(a) and 2(b), respectively. 6.2 L-shaped domain The second example is for the L-shaped domain with Ω = (, ) 2 \ [0, ] 2. Here, the domain in nonconvex and has a re-entrant corner at the origin, which causes a singularity in the first eigenfunction. To compute the eigenvalue error of the first eigenvalue, we take λ = as a reference value, cf. to the values in [7, 24] scaled by /4 due to the differently sized domain. The convergence results for uniform meshes in Figure 3(a) and 3(b) show reduced orders of convergence O(N ) for all eigenvalue errors and both error estimators. Furthermore, the convergence results based on adaptive refinement recover optimal higher order convergence 6
17 (a) (b) Figure 2: (a) Streamline plot of the discrete eigenfunction u l. p l = trσ l /2. (b) Plot of discrete pressure η 2 l (uniform) λ λl (uniform) λ λ l (uniform) η 2 l (adaptive) λ λl (adaptive) λ λ l (adaptive) order order 2 order µ 2 l (uniform) λ λl (uniform) λ λl (uniform) µ 2 l (adaptive) λ λl (adaptive) λ λl (adaptive) order order 2 order Nl Nl (a) (b) Figure 3: Convergence history of (a) λ λ l, λ λ l, η2 l (b) λ λ l, λ λ l and µ 2 l on uniformly and adaptively refined meshes for the L-shaped domain. 7
18 (a) Figure 4: (a) Adaptive refined meshes for ηl 2 with 392 nodes. µ 2 l (b) with 386 nodes, (b) Adaptively refined meshes for (a) (b) Figure 5: (a)streamline plot of the discrete eigenfunction u l. p l = trσ l /2. (b) Plot of discrete pressure 8
19 η 2 l (uniform) λ λl (uniform) λ λ l (uniform) η 2 l (adaptive) λ λl (adaptive) λ λ l (adaptive) order /2 order 2 order µ 2 l (uniform) λ λl (uniform) λ λl (uniform) µ 2 l (adaptive) λ λl (adaptive) λ λl (adaptive) order /2 order 2 order Nl Nl (a) (b) Figure 6: Convergence history of (a) λ λ l, λ λ l, η2 l, (b) λ λ l, λ λ l and µ 2 l on uniformly and adaptively refined meshes for the slit domain. O(N 3 ) of the post-processed eigenvalues λ l and λ l. In both cases the corresponding a posteriori error estimators ηl 2 and µ2 l are reliable and efficient and close to the true error. Moreover, the eigenvalue errors for adaptive refinement are several orders of magnitude below the eigenvalue errors for uniform refinement, which illustrates the importance of adaptive mesh refinement. Figures 4(a) and 4(b) show two adaptively refined meshes for the a posteriori error estimators ηl 2 and µ2 l, respectively. Both meshes show strong refinement toward the origin. Figures 5(a) and 5(b) show the discrete velocity and pressure as a streamline plot computed on an adaptive mesh. 6.3 Slit domain Finally, we consider the slit domain Ω = (, ) 2 \ {0 x, y = 0}, with maximal re-entrant corner of angle 2π. We take λ = as a reference value for the first eigenvalue. Again the first eigenfunction is singular. For uniform meshes in Figure 6(a) and 6(b), the convergence results show suboptimal convergence O(N /2 ) for the eigenvalue errors and error estimators. On the contrary, the convergence results, which are based on adaptive refinement, achieve optimal convergence O(N 3 ) for the post-processed eigenvalues λ l and λ l, and for the a posteriori error estimators ηl 2 and µ2 l. In Figures 6(a) and 6(b), the graphs for the estimators η2 l and µ2 l are parallel to the eigenvalue errors of the eigenvalues λ l and λ l, thus it confirms that both error estimator are numerically reliable and efficient. Note that the efficiency index corresponding to ηl 2 in Figure 6(a) and the efficiency index corresponding µ2 l in Figure 6(b) are close to one in case of adaptive mesh refinement. Figures 7(a) and 7(b) display adaptively refined meshes for the proposed error estimators ηl 2 and µ2 l, which both show strong refinement toward the re-entrant corner. The computed velocity and discrete pressure are shown in Figures 8(a) and 8(b) as a streamline plot on an adaptive mesh. 9
20 (a) (b) Figure 7: (a) Adaptive refined meshes for η 2 l with 375 nodes, (b) Adaptive refined meshes for µ2 l with 336 nodes (a) (b) Figure 8: (a)streamline plot of the discrete eigenfunction u l. p l = trσ l /2. (b) Plot of discrete pressure 20
21 References [] M. Ainsworth and J. T. Oden. A posteriori error estimation in finite element analysis. Pure and Applied Mathematics (New York). John Wiley & Sons, New York, [2] M. G. Armentano and V. Moreno. A posteriori error estimates of stabilized low-order mixed finite elements for the Stokes eigenvalue problem. J. Comput. Appl. Math., 269:32 49, 204. [3] D. N. Arnold and R. Winther. Mixed finite elements for elasticity. Numer. Math., 92(3):40 49, [4] I. Babuška and J. Osborn. Eigenvalue problems. In Handbook of numerical analysis. Volume II: Finite element methods (Part ), pages Amsterdam etc.: North-Holland, 99. [5] P. E. Bjørstad and B. P. Tjøstheim. High precision solutions of two fourth order eigenvalue problems. Computing, 63(2):97 07, 999. [6] D. Boffi. Finite element approximation of eigenvalue problems. Acta Numer., 9: 20, 200. [7] S. C. Brenner, P. Monk, and J. Sun. C 0 interior penalty Galerkin method for biharmonic eigenvalue problems. In Spectral and high order methods for partial differential equations ICOSAHOM 204, volume 06 of Lect. Notes Comput. Sci. Eng., pages 3 5. Springer, Cham, 205. [8] F. Brezzi. On the existence, uniqueness and approximation of saddle-point problems arising from Lagrangian multipliers. Rev. Française Automat. Informat. Recherche Opérationnelle Sér. Rouge, 8(R-2):29 5, 974. [9] F. Brezzi and M. Fortin. Mixed and hybrid finite element methods, volume 5 of Springer Series in Computational Mathematics. Springer-Verlag, New York, 99. [0] C. Carstensen. A unifying theory of a posteriori finite element error control. Numer. Math., 00(4):67 637, [] C. Carstensen, J. Gedicke, and E.-J. Park. Numerical experiments for the Arnold-Winther mixed finite elements for the Stokes problem. SIAM J. Sci. Comput., 34(4):A2267 A2287, 202. [2] W. Chen and Q. Lin. Approximation of an eigenvalue problem associated with the Stokes problem by the stream function-vorticity-pressure method. Appl. Math., 5():73 88, [3] W. Dörfler. A convergent adaptive algorithm for Poisson s equation. SIAM J. Numer. Anal., 33(3):06 24, 996. [4] J. Han, Z. Zhang, and Y. Yang. A new adaptive mixed finite element method based on residual type a posterior error estimates for the Stokes eigenvalue problem. Numer. Methods Partial Differential Equations, 3():3 53, 205. [5] P. Huang. Lower and upper bounds of Stokes eigenvalue problem based on stabilized finite element methods. Calcolo, 52():09 2,
22 [6] P. Huang, Y. He, and X. Feng. Numerical investigations on several stabilized finite element methods for the Stokes eigenvalue problem. Math. Probl. Eng., pages Art. ID , 4, 20. [7] O. A. Karakashian and F. Pascal. A posteriori error estimates for a discontinuous Galerkin approximation of second-order elliptic problems. SIAM J. Numer. Anal., 4(6): , [8] R. Lehoucq, D. Sorensen, and C. Yang. ARPACK Users Guide: Solution of Large-Scale Eigenvalue Problems with Implicitly Restarted Arnoldi Methods. Society for Industrial and Applied Mathematics, Philadelphia, PA. USA, 998. [9] H. Liu, W. Gong, S. Wang, and N. Yan. Superconvergence and a posteriori error estimates for the Stokes eigenvalue problems. BIT, 53(3): , 203. [20] C. Lovadina, M. Lyly, and R. Stenberg. A posteriori estimates for the Stokes eigenvalue problem. Numer. Methods Partial Differential Equations, 25(): , [2] S. Meddahi, D. Mora, and R. Rodríguez. A finite element analysis of a pseudostress formulation for the Stokes eigenvalue problem. IMA J. Numer. Anal., 35(2): , 205. [22] B. Mercier, J. Osborn, J. Rappaz, and P.-A. Raviart. Eigenvalue approximation by mixed and hybrid methods. Math. Comp., 36(54): , 98. [23] L. R. Scott and S. Zhang. Finite element interpolation of nonsmooth functions satisfying boundary conditions. Math. Comp., 54(90): , 990. [24] J. Sun and A. Zhou. Finite element methods for eigenvalue problems. Monographs and Research Notes in Mathematics. CRC Press, Boca Raton, FL, 207. [25] O. Türk, D. Boffi, and R. Codina. A stabilized finite element method for the two-field and three-field Stokes eigenvalue problems. Comput. Methods Appl. Mech. Engrg., 30: , 206. [26] R. Verfürth. A review of a posteriori error estimation and adaptive mesh-refinement techniques. Chichester: John Wiley & Sons; Stuttgart: B. G. Teubner,
NUMERICAL EXPERIMENTS FOR THE ARNOLD WINTHER MIXED FINITE ELEMENTS FOR THE STOKES PROBLEM
NUMERICAL EXPERIMENTS FOR THE ARNOLD WINTHER MIXED FINITE ELEMENTS FOR THE STOKES PROBLEM CARSTEN CARSTENSEN, JOSCHA GEDICKE, AND EUN-JAE PARK Abstract. The stress-velocity formulation of the stationary
More 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 informationarxiv: v1 [math.na] 15 Nov 2017
Noname manuscript No. (will be inserted by the editor) An HDG method with orthogonal projections in facet integrals Issei Oikawa arxiv:1711.05396v1 [math.na] 15 Nov 2017 Received: date / Accepted: date
More informationarxiv: v1 [math.na] 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 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 informationGeometric Multigrid Methods
Geometric Multigrid Methods Susanne C. Brenner Department of Mathematics and Center for Computation & Technology Louisiana State University IMA Tutorial: Fast Solution Techniques November 28, 2010 Ideas
More 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 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 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 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 informationYongdeok Kim and Seki Kim
J. Korean Math. Soc. 39 (00), No. 3, pp. 363 376 STABLE LOW ORDER NONCONFORMING QUADRILATERAL FINITE ELEMENTS FOR THE STOKES PROBLEM Yongdeok Kim and Seki Kim Abstract. Stability result is obtained for
More informationA 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 informationC 0 IPG Method for Biharmonic Eigenvalue Problems
C IPG Method for Biharmonic Eigenvalue Problems Susanne C. Brenner, Peter Monk 2, and Jiguang Sun 3 Department of Mathematics and Center for Computation & Technology, Louisiana State University, Baton
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 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 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 informationTo link to this article:
This article was downloaded by: [Wuhan University] On: 11 March 2015, At: 01:08 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 1072954 Registered office: Mortimer
More informationUnified A Posteriori Error Control for all Nonstandard Finite Elements 1
Unified A Posteriori Error Control for all Nonstandard Finite Elements 1 Martin Eigel C. Carstensen, C. Löbhard, R.H.W. Hoppe Humboldt-Universität zu Berlin 19.05.2010 1 we know of Guidelines for Applicants
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 informationELLIPTIC RECONSTRUCTION AND A POSTERIORI ERROR ESTIMATES FOR PARABOLIC PROBLEMS
ELLIPTIC RECONSTRUCTION AND A POSTERIORI ERROR ESTIMATES FOR PARABOLIC PROBLEMS CHARALAMBOS MAKRIDAKIS AND RICARDO H. NOCHETTO Abstract. It is known that the energy technique for a posteriori error analysis
More informationarxiv: v1 [math.na] 2 Sep 2015
The lower bound property of the Morley element eigenvalues arxiv:1509.00566v1 [math.na] Sep 015 Yidu Yang, Hao Li, Hai Bi School of Mathematics and Computer Science, Guizhou Normal University, Guiyang,
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 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 informationAxioms of Adaptivity (AoA) in Lecture 3 (sufficient for optimal convergence rates)
Axioms of Adaptivity (AoA) in Lecture 3 (sufficient for optimal convergence rates) Carsten Carstensen Humboldt-Universität zu Berlin 2018 International Graduate Summer School on Frontiers of Applied and
More informationEnhancing eigenvalue approximation by gradient recovery on adaptive meshes
IMA Journal of Numerical Analysis Advance Access published October 29, 2008 IMA Journal of Numerical Analysis Page 1 of 15 doi:10.1093/imanum/drn050 Enhancing eigenvalue approximation by gradient recovery
More 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 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 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 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 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 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 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 informationA Posteriori Existence in Adaptive Computations
Report no. 06/11 A Posteriori Existence in Adaptive Computations Christoph Ortner This short note demonstrates that it is not necessary to assume the existence of exact solutions in an a posteriori error
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 informationA gradient recovery method based on an oblique projection and boundary modification
ANZIAM J. 58 (CTAC2016) pp.c34 C45, 2017 C34 A gradient recovery method based on an oblique projection and boundary modification M. Ilyas 1 B. P. Lamichhane 2 M. H. Meylan 3 (Received 24 January 2017;
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 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 informationIMPROVED LEAST-SQUARES ERROR ESTIMATES FOR SCALAR HYPERBOLIC PROBLEMS, 1
Computational Methods in Applied Mathematics Vol. 1, No. 1(2001) 1 8 c Institute of Mathematics IMPROVED LEAST-SQUARES ERROR ESTIMATES FOR SCALAR HYPERBOLIC PROBLEMS, 1 P.B. BOCHEV E-mail: bochev@uta.edu
More informationA NEW FAMILY OF STABLE MIXED FINITE ELEMENTS FOR THE 3D STOKES EQUATIONS
MATHEMATICS OF COMPUTATION Volume 74, Number 250, Pages 543 554 S 0025-5718(04)01711-9 Article electronically published on August 31, 2004 A NEW FAMILY OF STABLE MIXED FINITE ELEMENTS FOR THE 3D STOKES
More 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 informationFinite element approximation of eigenvalue problems
Acta NumeTica (2010), pp. 1-120 doi: 10.1017/80962492910000012 Cambridge University Press, 2010 Finite element approximation of eigenvalue problems Daniele Boffi Dipartimento di Matematica, Universita
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 informationAxioms of Adaptivity (AoA) in Lecture 2 (sufficient for optimal convergence rates)
Axioms of Adaptivity (AoA) in Lecture 2 (sufficient for optimal convergence rates) Carsten Carstensen Humboldt-Universität zu Berlin 2018 International Graduate Summer School on Frontiers of Applied and
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 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 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 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:25. Fictitious domain finite element methods using cut elements: II. A stabilized Nitsche method ERIK BURMAN PETER HANSBO
PREPRINT 2010:25 Fictitious domain finite element methods using cut elements: II. A stabilized Nitsche method ERIK BURMAN PETER HANSBO Department of Mathematical Sciences Division of Mathematics CHALMERS
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 informationFULLY COMPUTABLE ERROR BOUNDS FOR EIGENVALUE PROBLEM
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 15, Number 1-2, Pages 260 276 c 2018 Institute for Scientific Computing and Information FULLY COMPUTABLE ERROR BOUNDS FOR EIGENVALUE PROBLEM
More 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 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 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 informationChapter 2 Finite Element Spaces for Linear Saddle Point Problems
Chapter 2 Finite Element Spaces for Linear Saddle Point Problems Remark 2.1. Motivation. This chapter deals with the first difficulty inherent to the incompressible Navier Stokes equations, see Remark
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 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 informationHybridized Discontinuous Galerkin Methods
Hybridized Discontinuous Galerkin Methods Theory and Christian Waluga Joint work with Herbert Egger (Uni Graz) 1st DUNE User Meeting, Stuttgart Christian Waluga (AICES) HDG Methods October 6-8, 2010 1
More 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 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 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 DECOMPOSITION RESULT FOR BIHARMONIC PROBLEMS AND THE HELLAN-HERRMANN-JOHNSON METHOD
Electronic ransactions on Numerical Analysis. Volume 45, pp. 257 282, 2016. Copyright c 2016,. ISSN 1068 9613. ENA A DECOMPOSIION RESUL FOR BIHARMONIC PROBLEMS AND HE HELLAN-HERRMANN-JOHNSON MEHOD WOLFGANG
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 informationSUPERCONVERGENCE PROPERTIES FOR OPTIMAL CONTROL PROBLEMS DISCRETIZED BY PIECEWISE LINEAR AND DISCONTINUOUS FUNCTIONS
SUPERCONVERGENCE PROPERTIES FOR OPTIMAL CONTROL PROBLEMS DISCRETIZED BY PIECEWISE LINEAR AND DISCONTINUOUS FUNCTIONS A. RÖSCH AND R. SIMON Abstract. An optimal control problem for an elliptic equation
More informationIntroduction. J.M. Burgers Center Graduate Course CFD I January Least-Squares Spectral Element Methods
Introduction In this workshop we will introduce you to the least-squares spectral element method. As you can see from the lecture notes, this method is a combination of the weak formulation derived from
More 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 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 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 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 informationA NONCONFORMING PENALTY METHOD FOR A TWO DIMENSIONAL CURL-CURL PROBLEM
A NONCONFORMING PENALTY METHOD FOR A TWO DIMENSIONAL CURL-CURL PROBLEM SUSANNE C. BRENNER, FENGYAN LI, AND LI-YENG SUNG Abstract. A nonconforming penalty method for a two-dimensional curl-curl problem
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 informationFinite Elements. Colin Cotter. February 22, Colin Cotter FEM
Finite Elements February 22, 2019 In the previous sections, we introduced the concept of finite element spaces, which contain certain functions defined on a domain. Finite element spaces are examples of
More 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 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 informationA POSTERIORI ERROR ESTIMATORS FOR CONVECTION-DIFFUSION EIGENVALUE PROBLEMS
A POSERIORI ERROR ESIMAORS FOR CONVECION-DIFFUSION EIGENVALUE PROBLEMS JOSCHA GEDICKE AND CARSEN CARSENSEN Abstract. A posteriori error estimators for convection-diffusion eigenvalue model problems are
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 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 informationA POSTERIORI ERROR ESTIMATES OF TRIANGULAR MIXED FINITE ELEMENT METHODS FOR QUADRATIC CONVECTION DIFFUSION OPTIMAL CONTROL PROBLEMS
A POSTERIORI ERROR ESTIMATES OF TRIANGULAR MIXED FINITE ELEMENT METHODS FOR QUADRATIC CONVECTION DIFFUSION OPTIMAL CONTROL PROBLEMS Z. LU Communicated by Gabriela Marinoschi In this paper, we discuss a
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 informationPSEUDO-COMPRESSIBILITY METHODS FOR THE UNSTEADY INCOMPRESSIBLE NAVIER-STOKES EQUATIONS
PSEUDO-COMPRESSIBILITY METHODS FOR THE UNSTEADY INCOMPRESSIBLE NAVIER-STOKES EQUATIONS Jie Shen Department of Mathematics, Penn State University University Par, PA 1680, USA Abstract. We present in this
More informationNew Discretizations of Turbulent Flow Problems
New Discretizations of Turbulent Flow Problems Carolina Cardoso Manica and Songul Kaya Merdan Abstract A suitable discretization for the Zeroth Order Model in Large Eddy Simulation of turbulent flows is
More 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 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 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 informationMultigrid Methods for Maxwell s Equations
Multigrid Methods for Maxwell s Equations Jintao Cui Institute for Mathematics and Its Applications University of Minnesota Outline Nonconforming Finite Element Methods for a Two Dimensional Curl-Curl
More information