Abstract and Applied Analysis

Volume 2014 (2014), Article ID 234375, 11 pages

http://dx.doi.org/10.1155/2014/234375

Research Article

Convergence Analysis of Incomplete Biquadratic Rectangular Element for Fourth-Order Singular Perturbation Problem on Anisotropic Meshes

Pingli Xie¹ and Meng Hu²

¹School of Sciences, Henan University of Technology, Zhengzhou 450001, China

²School of Mathematics and Statistics, Anyang Normal University, Anyang 455000, China

Received 2 June 2013; Accepted 29 December 2013; Published 13 February 2014

Academic Editor: Youyu Wang

Copyright © 2014 Pingli Xie and Meng Hu. 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.

Abstract.

The convergence analysis of a Morley type rectangular element for the fourth-order elliptic singular perturbation problem is considered. A counterexample is provided to show that the element is not uniformly convergent with respect to the perturbation parameter. A modified finite element approximation scheme is used to get convergent results; the corresponding error estimate is presented under anisotropic meshes. Numerical experiments are also carried out to demonstrate the theoretical analysis.

1. Introduction

The elliptic perturbation problems, which are derived from the stationary formation of parabolic perturbation problems, such as the Cahn-Hilliard type equation, are very important in both theoretical research and applications. The finite element methods are always chosen to be the appropriate way to solve the numerical solutions (cf. [1–5]). Here, we consider the following two-dimensional linear stationary Cahn-Hilliard type equation as our model problem:

where is the standard Laplace operator, is a bounded polygonal domain in , is the boundary of , and is a real parameter such that . Let denote the normal derivative of along the boundary . Particularly, the differential equations (1) formally degenerate to Poisson equations (a plate model degenerates towards an elastic membrane problem) when tends to zero.

Semper considers its conforming finite element methods in [5]. The regularity of the solution is analyzed, quasioptimal global error estimates are presented when , and local analysis is also done by using techniques of Nitsche and Schatz [6] and Schatz and Wahlbin [7]. The author points out that the method behaves poorly when the perturbation parameter is much smaller than the mesh size by some numerical experiments.

On the other hand, it is well known that when fourth-order problems are discretized by a finite element method, the standard variational formulation will require the piecewise smooth functions in space. However, it is very difficult to construct such functions, and even if we can do that, the element will be rather complicated. Hence a common approach to solve this problem is to use nonconforming finite elements which violate the -continuity requirement. In this case, two convergence criteria are generally employed: the Patch-Test [8] is used widely in engineers, but it is neither necessary nor sufficient; the Generalized Patch-Test [9] is proved to be the sufficient and necessary condition, while in practice it is often hard to be verified. To overcome the difficulty, the F-E-M criteria were proposed in [10] to make the test tractable.

Many successful nonconforming plate elements have been constructed (e.g., see [2, 3, 9–17]), but not all of them are convergent uniformly for (1) with respect to perturbation parameter . The very simple nonconforming Morley element (see [18]), which is convergent (cf. [19, 20]) and even has some superconvergent properties under uniform meshes for plate problems, see [21], however, is proved to be not uniformly convergent for (1) in [3] when , that is; it may diverge for second order problem like Poisson equation (see also [22]). It is considered that the main reason for this degeneracy is the fact that the finite element space is not a subspace of . Indeed, it is not of type. A counterexample is given in [3]. For more discussions on this element, we refer to [3, 19, 22]. As an alternative, a new modified element is proposed in [3], which is robust with respect to the parameter .

In [12], the convergence analysis of a nonconforming incomplete biquadratic rectangular plate element with the shape function space and the degrees of freedom and , respectively, is studied, where is the function value at the vertex of element , is the unit outer normal derivative value at the middle point of the edge of , and is the unit outer normal vector to . This element, similar to the famous triangular Morley element [3, 18, 19, 23], is also a non element, and its convergence order was given based on the Generalized Patch-Test. In [24], a modified element is provided by replacing the degrees of freedom and the shape function of [12] with and , respectively. Recently, [2] applied the modified Morley element of [24] to the fourth-order elliptic singular perturbation problem and proved the convergence uniformly in the perturbation parameter. However, all the studies above are based on the traditional regular triangulations.

In this paper we will present another improved element by using the same degrees of freedom of [2] and the same shape function space of [12]. Obviously, the above element is convergent for fourth-order plate bending problems according to FEM test in [10]. However, to our knowledge, there is no literature considering the convergence of this element for fourth-order singular perturbation problems. Here, we will show that this element is not uniformly convergent for fourth-order singular perturbation problems with respect to the perturbation parameter with a counterexample presented. Moreover, the convergence results are presented even under anisotropic meshes when the modified approximation formulation in [3] is employed.

The paper is organized as follows. The next section lists some preliminaries and the construction of the element. In Section 3, a counterexample is presented. In Section 4, the convergence results under the quasiuniform assumption and anisotropic meshes are provided. Numerical experiments are carried out in last section to confirm the theoretical analysis.

2. Premilinaries

Denote the inner product on by , the usual Sobolev space, norm, and seminorm by , , and , respectively. The space is the closure in of . Equivalently, we have

Let be the gradient of and let be the tensor of the second order partial derivatives . Define

Then the weak form of (1) reads: find , such that

By Green’s formula, it is easy to get

However, it does not hold on nonconforming finite element spaces.

Without loss of generality, we assume that the edges of are parallel to the axis. For mesh size , a rectangular triangulation of is then formed by lines also parallel to axis. Let be a rectangle with the central point , and the lengths of edges parallel to axis and axis, respectively, , , , , and the four vertices, (). Let be a reference element in () plane with central point (0,0), four vertices , , , and , and four edges  , (). Then there exists a reversible mapping :

On we define the finite element as:

where denotes the set of quadratic polynomials on element . Then it is easy to check that can be uniquely determined by the degrees of freedom . The degrees of freedom are plotted in Figure 1.

Figure 1: The degrees-of-freedom of incomplete biquadratic element.

For every , we define the interpolation operator as , and satisfies

such that

It can be checked that

Let be the associated finite element space defined by

where is the jump value of on and if .

Then the corresponding finite element approximation of (4) is as follows: find , such that

where for all ,

Next, we will present the modified discretization form of problem (4) in [2].

Let be the interpolation operator of the Lagrange bilinear rectangular element corresponding to the triangulation . The modified finite element method of (4) reads as: find , such that

Note that the problem has a unique solution when , but when , the problem degenerates to

in this case, is uniquely determined, though the solution is not unique.

We introduce the same mesh dependent norm and semi–norm on space :

The energy norm is defined by

3. A Counterexample

In this section, we construct a counterexample to show that the element presented in Section 2 is not convergent uniformly with respect to the perturbation parameter . That means if the element is applied to a nearly second order problem with the form of (1) when , the convergence rate of the method will deteriorate. In fact, like the Morley element, when it is applied to a second order equation like Poisson’s equation, the method will diverge.

As in [3], we consider the slightly modified reduced problem

to simplify some calculations. Assume , where and are disjoint subsets of . This problem is a second order problem with mixed boundary conditions, which can be regarded as the formal limit of the fourth-order problems

Let be a triangulation of , and let be the finite element space of incomplete biquadratic plate element corresponding to the boundary conditions of (18) here; see Figure 2. Then the approximation problem to (18) reads as: find , such that

where , denotes the arc length along . Moreover, the exact solution of (18) satisfies

where

Let be the corresponding energy norm of problem (18), that is, . Then, employing Cauchy-Schwarz inequality (we refer to [3]), we have

In the following, we will choose a suitable exact solution to prove the divergence of the method by virtue of (23).

Figure 2: The triangulation of .

The domain is taken as the unit square. To simplify the analytic process, the uniform triangulation with the mesh size is employed. Assume to be the intersection of with the coordinate axis, while to be the part of on and . Hence, the functions in finite element space are zeros at the vertices on the coordinate axis.

We assume that the exact solution of (18) is given by . Thus, on and on . Obviously, is harmonic, which means . Therefore:

Note that , thus and here is the finite element interpolation on . Similar to the discussion in [3], it is easy to derive that finite element space can be naturally decomposed into two spaces:

where and correspond to the vertex values and the edge values, respectively. In fact, the two spaces can be expressed as

where and represent the sets of the edges and vertices corresponding to the triangulation , respectively. Let be the interpolant of onto , then, according to its definition, we can get

We begin to prove that the limit is strictly positive, so we can show from (23) that the method is divergent.

In [3], the authors proved for Morley element, the decomposition of space is orthogonal, while it does not hold any more for this incomplete biquadratic plate element. However, in the mesh fashion chosen before, for any element , the length of each edge is . Let the center point of be , then by the definition of the space , direct calculation implies the expression of on element :

where

We denote the element on plane by the reference element . Apparently, the mapping from to is affine, hence we have

moreover, , , and then

thus

On the other hand, denoting the edges of on and by and , respectively, we have

We first consider the first term, let those components of lying on be , , and let the corresponding element be , the center point of be , where , . Then

thus

Therefore, by a limit process, we can obtain

Similarly, for , we have

this immediately leads to

Finally, since , and

we can verify that

From this expression we obtain that

This together with (32), (36), and (38) implies

The divergence of the method is therefore a consequence of the basic lower bound (23).

4. Convergence Analysis in a Modified Discretization Form

In the last section, we provide a counterexample to show that the incomplete biquadratic plate element may diverge for a second order problem like Poisson equation, and hence, in the standard finite element approximation, the convergence can not be insured for problem (1). In [2], a new modified approximation form is presented for Morley element and another Morley type rectangular element. In this section, we will show that the incomplete biquadratic plate element is convergent for problem (1) uniformly with respect to the parameter under anisotropic meshes.

To begin with, we introduce the following error estimate regarding the operator on anisotropic meshes, as to its proof, we refer to [25].

Lemma 1. For the bilinear interpolation operator , for all , there holds

where is a constant independent of triangulation.

Denote the quadratic part of the interpolant function by and the corresponding part of the function on the element by , then we can get the following lemma.

Lemma 2. For any , for all , without the regular or quasiuniform assumption, we have the following estimate:

Proof. Apparently, can be considered as an interpolant of on . For simplicity, we denote it by . We first show that on the reference element , for multi-index ,

For the convenience of notations, we denote degrees-of-freedom by , , . Correspondingly, , , . Hence, for all , direct computation provides the expression of :

where

thus, when , we have

Obviously, is a basis of . Moreover, let , we can get

By Hölder’s inequality, () is a bounded linear functional on . Therefore, by the basic theorem in [26], (45) holds. Then,

The conclusion when can be derived similarly, that is,

Equations (50) and (51) immediately imply the desired result.

The following error estimates can be found in [27].

Lemma 3. For all , without the regular or quasiuniform assumption, the following estimates hold:

The following theorem shows that for any fixed the new incomplete biquadratic element method converges linearly with respect to .

Theorem 4. Suppose and are the solutions of (1) and (14), respectively, , then without the quasiuniform assumption and regular condition, there exists a constant independent of and , such that

Proof. The second Strang lemma implies

where is the consistency error given by

Furthermore, from Lemmas 1 and 2,

Hence, it suffices to estimate the consistency error .

Since and are continuous, it follows from (1) that

By using the approximation formulation, the consistency error can therefore be expressed as

On the other hand,

and by Green’s formula, can be rewritten as

It follows from Hölder’s inequality and Lemma 1 that

Together with Lemma 3, we can immediately get the desired estimate.

Remark 5. By Lemma 1 and the estimates above, we can derive the same convergence result as in [2]:

We should mention that the result here does not need the quasiuniform assumptions. Moreover, similar discussions can also lead to the following estimate when the meshes satisfy the quasiuniform assumption

Remark 6. In the last section, a counterexample is presented to show the possible divergence of incomplete biquadratic plate element when applied to second order problem, but the theorem above implicates that, when the approximation formulation (14) is employed, the uniform convergence result can be ensured even without the regular condition or quasiuniform assumption.

5. Numerical Experiments

In this section, numerical experiments are carried out to confirm our theoretical analysis of the incomplete biquadratic element. We calculate several numerical examples for problem (1) in different approximation schemes. We consider problem (1) with and , where . The domain is divided into the following two fashions.

 Mesh 1: square mesh. The mesh obtained in this way for is illustrated in Figure 3(a). Mesh 2: each edge of is divided into segments with points , , . The mesh obtained in this way for is illustrated in Figure 3(b).

(a)

(b)

Figure 3: Mesh fashions: mesh 1 (a) and mesh 2 (b).

We first compute the relative errors in the energy norm under mesh 1 when we use the standard finite element approximation, for different and . For a comparison, we also consider the case when , that is, the Poisson equation with Dirichlet boundary conditions, and the biharmonic problem

for this case, the corresponding relative error is presented by . From Table 1 we can see that if the standard approximation scheme is applied, the method is divergent when .

Table 1: The errors employing approximation form (12) under mesh 1.

Poisson

Biharmonic

We first present the relative errors in the energy norm under mesh 1 when we use the standard finite element approximation.

To see the numerical effect more clearly, we plot the errors and under different meshes in a logarithm scale in Figure 4. Obviously the slope of the curve represents the convergence rate. We immediately get from the figures that when , the errors are no longer descending, which means that the method is divergent. We should point out that, for biharmonic equation, the method is very efficient.

(a)

(b)

Figure 4: The errors and under mesh 1 in the formulation (12).

In Tables 2 and 3, we list the relative errors and the errors for the approximation scheme (14) under different meshes. In this case, for biharmonic problem, is the solution of problem

Table 2: The errors employing approximation form (14) under mesh 1.

Poisson

Biharmonic

Table 3: The errors employing approximation form (14) under mesh 2.

Poisson

Biharmonic

When the formulation (14) is employed, the result shows that the method is uniformly convergent with respect to the parameter . Moreover, we can get that

At the same time, we plot the logarithm figures of the errors under different meshes in Figure 5. The results consist with our analysis.

(a)

(b)

Figure 5: The errors under mesh 1 (a) and mesh 2 (b).

From Figures 5, 6, and 7, we can also see that, with the numerical results derived under mesh 1 and mesh 2 differ slightly, the method under mesh 2 is still very efficient when the approximation scheme (14) is employed. This matches our theoretical analysis.

(a)

(b)

Figure 6: The error under mesh 1 (a) and mesh 2 (b).

(a)

(b)

Figure 7: The error under mesh 1 (a) and mesh 2 (b).

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The work is supported by the National Natural Sciences Foundation of China (no. 11126272 and 11201122), High-Level Personal Foundation of Henan University of Technology (Grant no. 2009BS066), and The basic and frontier project of Henan province (no. 132300410232).

References

S. C. Chen, Y. C. Zhao, and D. Y. Shi, “Non nonconforming elements for elliptic fourth order singular perturbation problem,” Journal of Computational Mathematics, vol. 23, no. 2, pp. 185–198, 2005.

M. Wang, J. C. Xu, and Y. C. Hu, “Modified Morley element method for a fourth order elliptic singular perturbation problem,” Journal of Computational Mathematics, vol. 24, no. 2, pp. 113–120, 2006.

T. K. Nilssen, X. C. Tai, and R. Winther, “A robust nonconforming -element,” Mathematics of Computation, vol. 70, no. 234, pp. 489–505, 2001.

S. Zhang and M. Wang, “A posteriori estimator of nonconforming finite element method for fourth order elliptic perturbation problems,” Journal of Computational Mathematics, vol. 26, no. 4, pp. 554–577, 2008.

B. Semper, “Conforming finite element approximations for a fourth-order singular perturbation problem,” SIAM Journal on Numerical Analysis, vol. 29, no. 4, pp. 1043–1058, 1992.

J. A. Nitsche and A. H. Schatz, “Interior estimates for Ritz-Galerkin methods,” Mathematics of Computation, vol. 28, pp. 937–958, 1974.

A. H. Schatz and L. B. Wahlbin, “On the finite element method for singularly perturbed reaction-diffusion problems in two and one dimensions,” Mathematics of Computation, vol. 40, no. 161, pp. 47–89, 1983.

G. Strang, “Variational crimes in the finite element method,” in The Mathematical Foundations of the Finite Element Method with Applications to Partial Differential Equations, A. R. Aziz, Ed., pp. 689–710, Academic Press, New York, NY, USA, 1972.

F. Stummel, “The generalized patch test,” SIAM Journal on Numerical Analysis, vol. 16, no. 3, pp. 449–471, 1979.

Z. C. Shi, “The F-E-M test for convergence of nonconforming finite elements,” Mathematics of Computation, vol. 49, no. 180, pp. 391–405, 1987.

S. C. Chen, D. Y. Shi, and H. Ichiro, “Generalized estimate formulation of nonconforming finite elements,” Mathematica Numerica Sinica, vol. 22, no. 3, pp. 295–300, 2000.

Z. C. Shi, “On the convergence of the incomplete biquadratic nonconforming plate element,” Mathematica Numerica Sinica, vol. 8, no. 1, pp. 53–62, 1986.

Z. C. Shi, “The generalized patch test for Zienkiewicz's triangles,” Journal of Computational Mathematics, vol. 2, no. 3, pp. 279–286, 1984.

Z. C. Shi and S. C. Chen, “An analysis of a nine-parameter plate element of Specht,” Mathematica Numerica Sinica, vol. 11, no. 3, pp. 312–318, 1989.

D. Y. Shi, S. P. Mao, and S. C. Chen, “On the anisotropic accuracy analysis of ACM's nonconforming finite element,” Journal of Computational Mathematics, vol. 23, no. 6, pp. 635–646, 2005.

S. P. Mao and S. C. Chen, “Convergence analysis of Morley element on anisotropic meshes,” Journal of Computational Mathematics, vol. 24, no. 2, pp. 169–180, 2006.

Z. C. Shi and Q. Y. Chen, “A new high accuracy rectangular element,” Scince in China A, vol. 30, no. 6, pp. 504–515, 2000.

L. S. D. Morley, “The triangular equilibrium element in the solution of plate bending problems,” Aero Quartet, vol. 19, pp. 149–169, 1968.

P. Lascaux and P. Lesaint, “Some nonconforming finite elements for the plate bending problem,” Analyse Numérique, vol. 9, no. 1, pp. 9–53, 1975.

Z. C. Shi, “Error estimates for the Morley element,” Mathematica Numerica Sinica, vol. 12, no. 2, pp. 113–118, 1990.

S. P. Mao and Z. C. Shi, “High accuracy analysis of two nonconforming plate elements,” Numerische Mathematik, vol. 111, no. 3, pp. 407–443, 2009.

M. Wang, “On the necessity and sufficiency of the patch test for convergence of nonconforming finite elements,” SIAM Journal on Numerical Analysis, vol. 39, no. 2, pp. 363–384, 2001.

P. G. Ciarlet, The Finite Element Method for Elliptic Problems, North-Holland, Amsterdam, The Netherlands, 1978.

H. Q. Zhang and M. Wang, The Mathematical Theory of Finite Elements, Science Press, Beijing, China, 1991.

T. Apel and M. Dobrowolski, “Anisotropic interpolation with applications to the finite element method,” Computing, vol. 47, no. 3-4, pp. 277–293, 1992.

S. C. Chen, D. Y. Shi, and Y. C. Zhao, “Anisotropic interpolation and quasi-Wilson element for narrow quadrilateral meshes,” IMA Journal of Numerical Analysis, vol. 24, no. 1, pp. 77–95, 2004.

P. L. Xie, Nonconforming finite element methods for fourth order elliptic problems [Ph.D. thesis], Fudan University, Shanghai, China, 2009.