Anisotropic Measures of Third Order Derivatives and the Quadratic Interpolation Error on Triangular Elements

The main purpose of this paper is to present a closer look at how the $H^1$- and $L^2$-errors for quadratic interpolation on a triangle are determined by the triangle geometry and the anisotropic behavior of the third order derivatives of interpolated functions. We characterize quantitatively the anisotropic behavior of a third order derivative tensor by its orientation and anisotropic ratio. Both exact error formulas and numerical experiments are presented for model problems of interpolating a cubic function $u$ at the vertices and the midpoints of three sides of a triangle. Based on the study on model problems, we conclude that when an element is aligned with the orientation of $\nabla^3 u$, the aspect ratio leading to nearly the smallest $H^1$- and $L^2$-norms of the interpolation error is approximately equal to the anisotropic ratio of $\nabla^3 u$. With this alignment and aspect ratio taken, the $H^1$-seminorm of the error is proportional to the reciprocal of the anisotropic ratio of $\nabla^3 u$, the $L^2$-norm of the error is proportional to the $-\frac 32$th power of the anisotropic ratio, and both of them are insensitive to the internal angles of the element.

[1]  Weiming Cao,et al.  On the Error of Linear Interpolation and the Orientation, Aspect Ratio, and Internal Angles of a Triangle , 2005, SIAM J. Numer. Anal..

[2]  Robert D. Russell,et al.  A Study of Monitor Functions for Two-Dimensional Adaptive Mesh Generation , 1999, SIAM J. Sci. Comput..

[3]  S. Rippa Long and thin triangles can be good for linear interpolation , 1992 .

[4]  Mark A. Taylor,et al.  An Algorithm for Computing Fekete Points in the Triangle , 2000, SIAM J. Numer. Anal..

[5]  J. Shewchuk What Is a Good Linear Finite Element? Interpolation, Conditioning, Anisotropy, and Quality Measures , 2002 .

[6]  Simona Perotto,et al.  New anisotropic a priori error estimates , 2001, Numerische Mathematik.

[7]  Weizhang Huang,et al.  Variational mesh adaptation II: error estimates and monitor functions , 2003 .

[8]  Weizhang Huang,et al.  Measuring Mesh Qualities and Application to Variational Mesh Adaptation , 2005, SIAM J. Sci. Comput..

[9]  Long Chen,et al.  Optimal anisotropic meshes for minimizing interpolation errors in Lp-norm , 2007, Math. Comput..

[10]  E. F. D'Azevedo,et al.  On optimal triangular meshes for minimizing the gradient error , 1991 .

[11]  T. Apel Anisotropic Finite Elements: Local Estimates and Applications , 1999 .

[12]  Philippe G. Ciarlet,et al.  6. Finite Element Methods for the Plate Problem , 2002 .

[13]  E. F. D’Azevedo,et al.  Optimal Triangular Mesh Generation by Coordinate Transformation , 1991, SIAM J. Sci. Comput..

[14]  R. B. Simpson Anisotropic mesh transformations and optimal error control , 1994 .

[15]  P. George,et al.  Delaunay mesh generation governed by metric specifications. Part I algorithms , 1997 .

[16]  I. Babuska,et al.  ON THE ANGLE CONDITION IN THE FINITE ELEMENT METHOD , 1976 .

[17]  LongChen,et al.  OPTIMAL DELAUNAY TRIANGULATIONS , 2004 .