- About this Journal
- Abstracting and Indexing
- Aims and Scope
- Annual Issues
- Article Processing Charges
- Articles in Press
- Author Guidelines
- Bibliographic Information
- Citations to this Journal
- Contact Information
- Editorial Board
- Editorial Workflow
- Free eTOC Alerts
- Publication Ethics
- Reviewers Acknowledgment
- Submit a Manuscript
- Subscription Information
- Table of Contents
Mathematical Problems in Engineering
Volume 2013 (2013), Article ID 829820, 9 pages
Superconvergence of a New Nonconforming Mixed Finite Element Scheme for Elliptic Problem
1Department of Mathematics and Physics, Luoyang Institute of Science and Technology, Luoyang 471003, China
2School of Mathematics and Statistics, Zhengzhou University, Zhengzhou 450001, China
Received 3 January 2013; Revised 30 June 2013; Accepted 5 July 2013
Academic Editor: Trung Nguyen Thoi
Copyright © 2013 Lifang Pei and Dongyang Shi. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
A new nonconforming mixed finite element scheme for the second-order elliptic problem is proposed based on a new mixed variational form. It has the lowest degrees of freedom on rectangular meshes. The superclose property is proven by employing integral identity technique. Then global superconvergence result is derived through interpolation postprocessing operators. At last, some numerical experiments are carried out to verify the theoretical analysis.
Mixed finite element method (MFEM) is an important branch of FEMs and has been used widely in numerical computation of practical problems. A lot of studies on this aspect have been devoted to the second-order elliptic problems [1–5], in which two approximation spaces of MFEM should satisfy the famous B-B condition . However, since the variable and flux belong to and , respectively, it is not easy to construct a stable MFE space pair. In order to circumvent or ameliorate this deficiency, many approaches have been proposed, such as the least squares FEM , stabilization FEM , and -Galerkin FEM . Recently,  presented a new MFEM, in which B-B condition is automatically satisfied when and meet a relation of inclusion; that is, , where and are finite element approximation spaces of and flux , respectively. This advantage makes the construction of stable MFE space pair extremely simple and convenient. A family of triangular and rectangular conforming MFE space pairs with lower degrees of freedom is constructed in , in which the total degrees of freedom of first-order and second-order MFE schemes are about and , respectively; herein denotes the number of nodal points in subdivision. Later,  derived the similar results with  for conforming MFEM and gave a numerical example.
In this paper, motivated by the idea of [10, 11], we first construct a new nonconforming MFEM (NMFEM). The original variable is approximated by the constrained element space  and the flux by piecewise constant vectors space, respectively. Note that the total degrees of freedom of the NMFEM are only about , the lowest on rectangular meshes. We prove that it satisfies B-B condition. Then by the use of integral identity technique, we derive the superclose property for in energy norm and flux in norm. Furthermore, the global superconvergence result with order is obtained through interpolation postprocessing operators. Finally, some numerical results are provided to verify the theoretical analysis. It is observed that, compared with FEM using 4-node quadrilateral (FEM-Q4) and MFEM ( and flux are approximated by piecewise constants and the Raviart-Thomas element, resp.), NMFEM behaves well and has higher rates than MFEM, and it has almost the same rate as FEM-Q4 for in norm, but a higher rate than FEM-Q4 for in energy norm. Moreover, NMFEM is effective and accurate for the diffusion problem studied in .
The remainder of this paper is organized as follows. In Section 2, we introduce NMFEM and derive the superclose property and superconvergence results. In Section 3, we carry out some numerical experiments to show the performance of NMFEM.
We will use standard notations for the Sobolev spaces with norm and seminorm , with norm and seminorm , where is an integer. Besides, let and be the norm and norm, respectively. Throughout the paper, denotes a positive constant independent of the mesh parameter and may be different at each appearance.
2. Superconvergence Analysis for NMFEM
Consider the following elliptic problem: where is a bounded convex polygon domain, .
Obviously, and are continuous bilinear functionals, is a continuous linear functional, and for all . Moreover, for , we have . So there exists a constant such that that is, the B-B condition is satisfied, and therefore (3) has a unique solution .
Let , and then (3) can be written as
Let be a rectangular partition of the domain . For a given element , its four vertices are denoted by , , , , and four edges by , . Let be the reference element with nodes , , , and edges .
Define the affine mapping by Let be the space of polynomials with degrees defined on , and then the constrained element space is defined by : where denotes the jump value of across the boundary , and if .
Let denote the number of interior nodes. It has been proven in  that and , where are defined associated with nodes of as
We choose the following FE spaces and to approximate and , respectively: where is the space of constants on . Obviously, the total degree of freedoms of the nonconforming MFE space pair is , and is a norm over .
Then the MFE approximation of problem (3) is to find such that where .
For , from the affine mapping and the definition of , we can get , and then , that is, . Thus that is, the discrete B-B condition holds, and (10) has a unique solution .
For all , let , and then we can prove the following two important lemmas.
Lemma 1. For all , the following inequality holds:
Proof. In order to prove (12), we will construct a pair satisfying
On the other hand, for a given arbitrary but fixed , by the discrete B-B condition, there exists a such that
So, for , there holds
As a result, setting , we have
which implies that
Note that and we have which follows the desired result (12).
Let and denote the associated interpolation operators of and conforming bilinear element space, respectively. Let be the element interpolation operator ; that is, for all , satisfying , where ,. It has been shown in  that, for all , because is a bubble function for .
Lemma 2. Assume that , , we have where is a interpolation operator satisfying .
Proof. Since (23) has been proven by one of the authors in , we only need to prove (21) and (22).
In fact, because and are constants on each , we have which is (21).
On the other hand, note that Since is a constant vector on and , it follows from integration by parts that Furthermore, let . Note that and . Employing integral identity technique , we have Similarly, Therefore, from (25)–(28), we have The proof is completed.
Now we start to state the following superclose property.
Proof. For , from (2) and (10), we have
Applying (31) yields
Using (12) in Lemma 1 and (32), we can obtain
Hence the desired result follows from the interpolation theorem and Lemma 2.
In order to derive global superconvergence for and flux , we introduce the following postprocessing operators and as , , and , , where are the value of on the nodes and are nodes of on macroelement , while consists of the four small elements in , and and are bilinear and biquadratic piecewise polynomials spaces, respectively. It can be checked that the following properties hold: Then we can have the following superconvergence result.
Theorem 4. Under the assumptions in Theorem 3, there holds
3. Numerical Experiments
In this section, some numerical examples and comparison with other methods are presented to confirm theoretical analysis and good performance of NMFEM.
We consider the problem (1) with and the exact solution is , and then the flux field can be expressed as . We divide the domain into a family of quasiuniform rectangles with number of . The figures of exact solution of problem (1) and finite element approximation with are plotted in Figures 1, 2, and 3, respectively.
In Tables 1 and 2, we present the superclose and superconvergence results of the original variable in energy norm and flux in norm with , respectively. It is clearly that , , , and are converged at order 2 with respect to , which coincide with our theoretical analysis in Theorems 3 and 4. In order to describe the results more intuitively, we plot the errors in the logarithm scales in Figure 4.
Moreover, we compare the results of NMFEM with those of FEM using 4-node quadrilateral (FEM-Q4) and MFEM ( and flux are approximated by piecewise constants and the Raviart-Thomas element, resp.).
The convergence rates of errors of in and energy norm are shown in Figures 5 and 6. The comparison of the flux in norm of NMFEM with MFEM is also given in Figure 7. It is observed that (a) the convergence rates of and flux in norm of NMFEM are better than those of MFEM; (b) NMFEM has almost the same rate as FEM-Q4 for in norm, but a higher rate than FEM-Q4 for in energy norm.
Furthermore, a numerical experiment is carried out to demonstrate the effectiveness and accuracy of NMFEM for the following diffusion problem: where is the source term and , is the boundary data, the permeability is a symmetric tensor-valued function such that (a) is piecewise Lipschitz-continuous on and (b) the set of the eigenvalues of is included in (with ) for all .
As , we consider the problem (37) with , and the analytical solution Let , we can get the superclose and superconvergence results of in energy norm and flux in norm. The errors are listed in Tables 3 and 4 and plotted in the logarithm scales in Figure 8, respectively.
The authors would like to express their sincere thanks to the anonymous referee for his many helpful suggestions, which contribute significantly to the improvement of the paper. The research is supported by the NSF of China (no. 10971203; no. 11271340), Research Fund for the Doctoral Program of Higher Education of China (no. 20094101110006), and Foundation of He’nan Educational Committee (no. 13B110144).
- P.-A. Raviart and J. M. Thomas, “A mixed finite element method for 2nd order elliptic problems,” in Mathematical Aspects of Finite Element Methods, vol. 606 of Lecture Notes in Mathematics, pp. 292–315, Springer, Berlin, Germany, 1977.
- M. Farhloul and M. Fortin, “A nonconforming mixed finite element for second-order elliptic problems,” Numerical Methods for Partial Differential Equations, vol. 13, no. 5, pp. 445–457, 1997.
- B. Cockburn, G. Kanschat, I. Perugia, and D. Schötzau, “Superconvergence of the local discontinuous Galerkin method for elliptic problems on Cartesian grids,” SIAM Journal on Numerical Analysis, vol. 39, no. 1, pp. 264–285, 2001.
- Z.-D. Luo, “Method of nonconforming mixed finite element for second order elliptic problems,” Journal of Computational Mathematics, vol. 18, no. 5, pp. 449–456, 2000.
- D. Y. Shi and C. X. Wang, “A new low-order non-conforming mixed finite-element scheme for second-order elliptic problems,” International Journal of Computer Mathematics, vol. 88, no. 10, pp. 2167–2177, 2011.
- F. Brezzi and M. Fortin, Mixed and Hybrid Finite Element Methods, vol. 15 of Springer Series in Computational Mathematics, Springer, New York, NY, USA, 1991.
- P. B. Bochev and M. D. Gunzburger, “Finite element methods of least-squares type,” SIAM Review, vol. 40, no. 4, pp. 789–837, 1998.
- P. B. Bochev, C. R. Dohrmann, and M. D. Gunzburger, “Stabilization of low-order mixed finite elements for the Stokes equations,” SIAM Journal on Numerical Analysis, vol. 44, no. 1, pp. 82–101, 2006.
- D.-Y. Shi and H.-H. Wang, “Nonconforming -Galerkin mixed FEM for Sobolev equations on anisotropic meshes,” Acta Mathematicae Applicatae Sinica, vol. 25, no. 2, pp. 335–344, 2009.
- S. C. Chen and H. R. Chen, “New mixed element schemes for a second-order elliptic problem,” Mathematica Numerica Sinica, vol. 32, no. 2, pp. 213–218, 2010.
- F. Shi, J. P. Yu, and K. T. Li, “A new mixed finite element scheme for elliptic equations,” Chinese Journal of Engineering Mathematics, vol. 28, no. 2, pp. 231–237, 2011.
- J. Hu, H. Y. Man, and Z. C. Shi, “Constrained nonconforming rotated element for Stokes flow and planar elasticity,” Mathematica Numerica Sinica, vol. 27, no. 3, pp. 311–324, 2005.
- J. Droniou and C. le Potier, “Construction and convergence study of schemes preserving the elliptic local maximum principle,” SIAM Journal on Numerical Analysis, vol. 49, no. 2, pp. 459–490, 2011.
- R. Rannacher and S. Turek, “Simple nonconforming quadrilateral Stokes element,” Numerical Methods for Partial Differential Equations, vol. 8, no. 2, pp. 97–111, 1992.
- D. Y. Shi, Y. C. Peng, and S. C. Chen, “Error estimates for rotated element approximation of the eigenvalue problem on anisotropic meshes,” Applied Mathematics Letters, vol. 22, no. 6, pp. 952–959, 2009.
- Q. Lin and N. N. Yan, The Construction and Analysis of High Accurate Finite Element Methods, Hebei University Press, Baoding, China, 1996.