arxiv: v1 [math.na] 1 May 2013

Transcription

1 arxiv: v [mathna] May 03 Approximation Properties of a Gradient Recovery Operator Using a Biorthogonal System Bishnu P Lamichhane and Adam McNeilly May, 03 Abstract A gradient recovery operator based on projecting the discrete gradient onto the standard finite element space is considered We use an oblique projection, where the test and trial spaces are different, and the bases of these two spaces form a biorthogonal system Biorthogonality allows efficient computation of the recovery operator We analyse the approximation properties of the gradient recovery operator Key words Gradient reconstruction, oblique projection, biorthogonal system AMS subject classification 65N30, 65N5, 65N50 Introduction The gradient reconstruction is a popular technique to develop a reliable a posteriori error estimators for approximating the solution of partial differential equations using adaptive finite element methods [, 3,, 9, 3, ] Recently we have presented a gradient reconstruction operator based on an oblique projection [5] The oblique projection operator is constructed by using a biorthogonal system In fact, for the linear finite element in simplicial meshes, this approach reproduces the so-called gradient reconstruction scheme by the weighted averaging [4, 9, ] We proved that the approximation property of the recovered gradient for any finite element space is similar to the one obtained by using the orthogonal projection with respect to L -norm [9, 3] In this article, we aim at analysing the approximation property of the recovered gradient in one School of Mathematical and Physical Sciences, University of Newcastle, Callaghan, NSW 308, bishnulamichhane@newcastleeduau, adammcneilly@uoneduau

2 dimension using the oblique projection This construction is quite useful in extending the weighted average gradient recovery of linear finite elements [9] to quadrilaterals and hexahedras Let = (α, β) with α, β R and α < β Let = {α = < x < < x n = β} be a partition of the interval We define the interior of the grid, denoted int ( ), as int ( ) = {x i : i n } We also define the set of intervals in the partition as {I i } n i=0, where I i = [x i, x i+ ) Two sets A n and B n of indices are also defined as A n = {i N : 0 i n} and B n = A n \ {0, n}, respectively A piecewise linear interpolant of a continuous function u is written as I h u V h with I h u(x) = n u(x i )φ i (x), i=0 where φ i is the standard hat function associated with the point x i, 0 i n We define a discrete space, V h = span {φ 0,, φ n } H () The linear interpolant of u H () is the continuous function defined by I h u = n i=0 u(x i)φ i However, if we compute the derivative of this interpolant I h u, the resulting function will not be continuous To make the derivative continuous we project the derivative of the interpolant, u h x = n i=0 u(x i) φi x, onto the discrete space V h There are two different types of projection One is an orthogonal projection and the other is an oblique projection The orthogonal projection operator, P h, that projects u h x onto V u h is to find a g h = P h h x V h that satisfies: u h g h φ j dx = x φ j dx () Since g h V h, we can represent it as an (n + )-dimensional vector: g = g 0 g n with g h = n g i φ i Now the requirement given in equation () is equivalent to a linear system: M g = f, where M is a mass matrix, and f = f 0 f n with f j = i=0 u h x φ j dx

3 Here the mass matrix M is tridiagonal We can reduce computation time greatly if we have a diagonal mass matrix This can be done if we use a suitable oblique projection instead of an orthogonal projection We consider the projection Q h : L () V h, du which is defined as the problem of finding g h = Q h dx V h such that g h λ h dx = dx λ h dx, λ h M h, where M h is another piecewise polynomial space, not orthogonal to V h, with dim (M h ) = dim (V h ), see [8] In fact, the projection operator Q h is well-defined due to the following stability condition There is a constant β > 0 independent of the mesh-size h such that [5, 6] φ h L () β sup µ h M h \{0} µ hφ h dx, φ h V h µ h L () In order to achieve that the mass matrix M is diagonal we need to define a new set of basis functions for M h, {λ 0,, λ n }, that are biorthogonal to the standard hat basis function we used previously This biorthogonality relation is defined as: λ i φ j dx = c j δ ij, c j 0, i, j n, () where δ ij is the Kronecker delta function: {, if i = j δ ij = 0, otherwise, and c j is a positive scaling factor The basis functions for M h are simply given by { (x x)+(x ) λ 0 (x) = x, x x 0, otherwise { (x xn )+(x x n) x λ n (x) = n x n, x n x x n, 0, otherwise and for i n (x x i )+(x x i) x i x i, x i x x i (x x λ i (x) = i+)+(x x i) x i x i+, x i x x i+ 0, otherwise 3

4 Figure : The hat basis function (left) and biorthogonal basis function (right) with stepsize h = By using an oblique projection Q h the mass matrix will be diagonal We let the diagonal mass matrix be D, so that our system is D g = f The values g i are our estimates of the gradient of u at the point x i So, we estimate the gradient by finding g = D f, where g i = dx λ idx φ i λ i dx We want to calculate the error in this approximation, and find out when g i approximates u (x i ) exactly for each x i As in [0, ] we want to see if g i approximates u (x i ) exactly when u is a quadratic polynomial Superconvergence Theorem Let u C 0 () Then we have g i = u(x i+) u(x i ) x i+ x i, i B n, and Proof: g 0 = u(x ) u( ) x, g n = u(x n) u(x n ) x n x n We note that (x x i )+(x x i) x i x i, x i x x i (x x λ i (x) = i+)+(x x i) x i x i+ x i x x i+, for all i B n 0, otherwise Now, we calculate g i for i B n : g i = dx λ i dx φ i λ i dx, 4

5 where φ i λ i dx = + xi x i xi+ x i ( ) ( x xi (x xi ) + (x x i ) x i x i x i x i ( ) ( x xi+ (x xi+ ) + (x x i ) x i x i+ x i x i+ ) dx ) dx = (x i x i+ ), and dx λ i dx = n j=0 = u(x i ) dφ j u(x j ) dx λ i dx xi x i dφ i dx λ i dx + u(x i ) ( xi x i xi+ dφ i+ + u(x i+ ) x i dx λ i dx (since φ j and λ i overlap only when j {i, i, i + } ) Therefore, g i = u(x i+) u(x i ) x i+ x i ) dφ xi+ i dx λ dφ i i dx + x i dx λ i dx Now we look at the end-points We note that g 0 = Computing as before we get dx λ 0 dx φ 0 λ 0 dx, and g n = dx λ n dx φ n λ n dx g 0 = u(x ) u( ) x, g n = u(x n) u(x n ) x n x n We have the following super-convergence in L -norm This is proved as in [7, 9] Theorem Let h i = x i+ x i for 0 i n, h = ma i n h i, and 0 = (x, x n ) If the point distribution satisfies h i+ h i = O(h ) for 0 i n Then we have the estimate du dx Q h dx h u W 3, ( 0), u W 3 () L ( 0) For the tensor product meshes in two or three dimensions satisfying the above mesh condition this theorem has an easy extension 5

6 Application to quadratic functions Corollary 3 Let u P (R) Then g i reproduces u ( x i ) exactly for all x i, where: +x, i = 0 x x i = i +x i+ x, i n, and g i = 0 dx λ i dx x n +x n x, i = n 0 φ i λ i dx Figure : A non-uniform grid with 8 nodes (vertical lines) The points x i (dots) are also shown Proof: We use the result of the previous theorem to get ( ax i + bx i ax i+ bx i+) g i = (x i x i+ ) = a ( x i ) x i+ + b (xi x i+ ) x i x i+ = a (x i x i+ ) (x i + x i+ ) + b (x i x i+ ) x i x i+ = a (x i + x i+ ) + b On the other hand, ( ) u xi + x i+ ( x i ) = a + b = a (x i + x i+ ) + b = g i 6

7 So, g i reproduces u ( x i ) exactly for i B n Now for i = 0 and i = n, we have ( ) ax g 0 = + bx ax 0 b (x = a ( ) x x 0 + b (x ) ) x = a (x ) (x + ) + b (x ) = a (x + ) + b, x x ( 0 ) ax and g n = n + bx n ax n bx n (x = a ( ) x n x n + b (xn x n ) n x n ) x n x n = a (x n x n ) (x n + x n ) + b (x n x n ) = a (x n + x n ) + b x n x n Since u ( ) = a ( + x ) + b, and u ( x n ) = a (x n + x n ) + b, we have g 0 and g n reproduce u ( ) and u ( x n ), respectively, exactly Remark 4 (Uniform Grid) Let = {α = < x < < x n = β} be a uniform grid on the interval so that x i x i = h, i A n \ {0}, where h is some constant, called the stepsize We note that if our grid is uniform, then x i = x i x i int ( ) So, our gradient recovery operator will reproduce the exact gradient of any quadratic function on the interior of a uniform grid We cannot recover the gradients at the endpoints exactly, however, since and x n x n Corollary 5 Let u P (R) with u(x) = ax +bx+c, and let the grid be uniform with stepsize h and = 0 Then g i u (x i ) = ah for i = 0, n (ie for the endpoints of the grid) Proof: We will start with the case where i = 0 (ie the left endpoint) We know from Theorem 3 that g 0 = a (x ) + b Since our grid is uniform with stepsize h, this simplifies to g 0 = ah + b u ( ) = a + b = b, since = 0 Therefore g 0 u ( ) = (ah + b) b = ah The case for i = n (ie the right endpoint) is proven similarly For a non-uniform grid, we cannot simplify our approximations using the stepsize h, since the spacing between each adjacent node is not always equal We did not make any assumption about the uniformity of the grid in Theorem 3 Thus g i u (x i ), i A n, is not zero for a non-uniform grid This is estimated in the following corollary 7

8 Corollary 6 Let u (x) = ax + bx + c Then, g i u (x i ) = a (x i + x i+ x i ) x i int ( ), g 0 u ( ) = a (x ), and g n u (x n ) = a (x n x n ) Remark 7 For i B n let h i = x i+ x i Then we have g i u (x i ) = a (x i + x i+ x i ) = a(h i+ h i ) We still get superapproximation of the gradient recovery when h i+ h i = O(h ) when i B n Application to cubic functions Corollary 8 Let u P 3 (R) with u(x) = ax 3 + bx + cx + d Then, g i u ( x i ) = a 4 (x i x i+ ) for all i B n, and g 0 u ( ) = a 4 ( x ) g n u ( x n ) = a 4 (x n x n ) where x i is defined as in Theorem 3 Similarly, for all i B n we have g i u (x i ) = a ( x i + x i x i+ + x i+ 3x ) i + b (xi + x i+ x i ), and g 0 u ( ) = a ( x + x x 0) + b (x ), g n u (x n ) = a ( x n + x n x n x n) + b (xn x n ) Proof: The proof of this theorem is similar to Theorem 3 3 Conclusion We have presented an analysis of approximation property of the reconstructed gradient using an oblique projection The reconstruction of the gradient is numerically efficient due to the use of a biorthogonal system It is useful to investigate the extension to higher order finite elements 8

9 References [] M Ainsworth and J T Oden A Posteriori Error Estimation in Finite Element Analysis Wiley Interscience, New York, 000 [] J Chen and D Wang Three-dimensional finite element superconvergent gradient recovery on par6 patterns Numerical Marthematics: Theory, Methods and Applications, 3:78 94, 00 [3] L Chen Superconvergence of tetrahedral linear finite elements International Journal of Numerical Analysis and Modeling,, 3:73 8, 006 [4] J Goodsell Pointwise superconvergence of the gradient for the linear tetrahedral element Numerical Methods for Partial Differential Equations, 0:65 666, 994 [5] BP Lamichhane A gradient recovery operator based on an oblique projection Electronic Transactions on Numerical Analysis, 37:66 7, 00 [6] BP Lamichhane A stabilized mixed finite element method for the biharmonic equation based on biorthogonal systems Journal of Computational and Applied Mathematics, 35: , 0 [7] B Li and Z Zhang Analysis of a class of superconvergence patch recovery techniques for linear and bilinear finite elements Numerical Methods for Partial Differential Equations, 5:5 67, 999 [8] BI Wohlmuth Discretization Methods and Iterative Solvers Based on Domain Decomposition, volume 7 of LNCS Springer Heidelberg, 00 [9] J Xu and Z Zhang Analysis of recovery type a posteriori error estimators for mildly structured grids Mathematics of Computation, 73:39 5, 004 [0] Z Zhang Ultraconvergence of the patch recovery technique Mathematics of Computation, 65:43 437, 996 [] Z Zhang Ultraconvergence of the patch recovery technique ii Mathematics of Computation, 69:4 58, 000 [] OC Zienkiewicz and JZ Zhu The superconvergent patch recovery and a posteriori error estimates part : The recovery technique International Journal for Numerical Methods in Engineering, 33:33 364, 99 [3] OC Zienkiewicz and JZ Zhu The superconvergent patch recovery and a posteriori error estimates part : Error estimates and adaptivity International Journal for Numerical Methods in Engineering, 33:365 38, 99 9