Skip to content
Licensed Unlicensed Requires Authentication Published by De Gruyter March 9, 2024

Unifying a posteriori error analysis of five piecewise quadratic discretisations for the biharmonic equation

  • Carsten Carstensen EMAIL logo , Benedikt Gräßle and Neela Nataraj

Abstract

An abstract property (H) is the key to a complete a priori error analysis in the (discrete) energy norm for several nonstandard finite element methods in the recent work [Lowest-order equivalent nonstandard finite element methods for biharmonic plates, Carstensen and Nataraj, M2AN, 2022]. This paper investigates the impact of (H) to the a posteriori error analysis and establishes known and novel explicit residual-based a posteriori error estimates. The abstract framework applies to Morley, two versions of discontinuous Galerkin, C0 interior penalty, as well as weakly over-penalized symmetric interior penalty schemes for the biharmonic equation with a general source term in H−2 (Ω).

JEL Classification: 65N30; 65N12; 65N50

Distinguished Visiting Professor, Department of Mathematics, Indian Institute of Technology Bombay, Powai, Mumbai-400076, India.

Funding statement: The research of the first two authors have been supported by the Deutsche Forschungsgemeinschaft in the Priority Program 1748 under the project ‘Foundation and Application of Generalized Mixed FEM Towards Nonlinear Problems in Solid Mechanics’ (CA 151/22-2). This paper has been supported by SPARC project (id 235) ‘The Mathematics and Computation of Plates’ and SERB POWER Fellowship SPF/2020/000019. The second author is supported by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under Germany’s Excellence Strategy – The Berlin Mathematics Research Center MATH+ (EXC-2046/1, project ID: 390685689).

References

[1] G. A. Baker, Finite element methods for elliptic equations using nonconforming elements, Math. Comp. 31 (1977), No. 137, 45–59.10.1090/S0025-5718-1977-0431742-5Search in Google Scholar

[2] L. Beirão da Veiga, J. Niiranen, and R. Stenberg, A posteriori error estimates for the Morley plate bending element, Numer. Math. 106 (2007), No. 2, 165–179.10.1007/s00211-007-0066-1Search in Google Scholar

[3] D. Braess, Finite Elements, Theory, Fast Solvers, and Applications in Elasticity Theory, 3rd ed., Cambridge, 2007.10.1017/CBO9780511618635Search in Google Scholar

[4] S. C. Brenner, T. Gudi, and L.-Y. Sung, An a posteriori error estimator for a quadratic C0-interior penalty method for the biharmonic problem, IMA J. Numer. Anal. 30 (2010), No. 3, 777–798.10.1093/imanum/drn057Search in Google Scholar

[5] S. C. Brenner, T. Gudi, and L.-Y. Sung, A weakly over-penalized symmetric interior penalty method for the biharmonic problem, Electron. Trans. Numer. Anal. 37 (2010), 214–238.10.1007/s10915-009-9278-0Search in Google Scholar

[6] S. C. Brenner and L. R. Scott, The Mathematical Theory of Finite Element Methods, 3rd ed., Springer, 2007.10.1007/978-0-387-75934-0Search in Google Scholar

[7] S. C. Brenner and L.-Y. Sung, C0 interior penalty methods for fourth order elliptic boundary value problems on polygonal domains, J. Sci. Comput. 22/23 (2005), 83–118.10.1007/s10915-004-4135-7Search in Google Scholar

[8] S. C. Brenner, L.-Y. Sung, H. Zhang, and Y. Zhang, A Morley finite element method for the displacement obstacle problem of clamped Kirchhoff plates, J. Comput. Appl. Math. 254 (2013), 31–42.10.1016/j.cam.2013.02.028Search in Google Scholar

[9] S. C. Brenner, Convergence of nonconforming multigrid methods without full elliptic regularity, Math. Comp. 68 (1999), No. 225, 25–53.10.1090/S0025-5718-99-01035-2Search in Google Scholar

[10] C. Carstensen, S. Bartels, and S. Jansche, A posteriori error estimates for nonconforming finite element methods, Numer. Math. 92 (2002), No. 2, 233–256.10.1007/s002110100378Search in Google Scholar

[11] C. Carstensen, M. Eigel, R. H.W. Hoppe, and C. Löbhard, A review of unified a posteriori finite element error control, Numer. Math. Theory Methods Appl. 5 (2012), No. 4, 509–558.10.4208/nmtma.2011.m1032Search in Google Scholar

[12] C. Carstensen and D. Gallistl, Guaranteed lower eigenvalue bounds for the biharmonic equation, Numer. Math. 126 (2014), No. 1, 33–51.10.1007/s00211-013-0559-zSearch in Google Scholar

[13] C. Carstensen, D. Gallistl, and J. Hu, A posteriori error estimates for nonconforming finite element methods for fourth-order problems on rectangles, Numer. Math. 124 (2013), No. 2, 309–335.10.1007/s00211-012-0513-5Search in Google Scholar

[14] C. Carstensen, D. Gallistl, and J. Hu, A discrete Helmholtz decomposition with Morley finite element functions and the optimality of adaptive finite element schemes, Comput. Math. Appl. 68 (2014), No. 12, part B, 2167–2181.10.1016/j.camwa.2014.07.019Search in Google Scholar

[15] C. Carstensen, D. Gallistl, and N. Nataraj, Comparison results of nonstandard P2 finite element methods for the biharmonic problem, ESAIM Math. Model. Numer. Anal. (2015), 977–990.10.1051/m2an/2014062Search in Google Scholar

[16] C. Carstensen, D. Gallistl, and M. Schedensack, Adaptive nonconforming Crouzeix–Raviart FEM for eigenvalue problems, Math. Comp. 84 (2015), 1061–1087.10.1090/S0025-5718-2014-02894-9Search in Google Scholar

[17] C. Carstensen, J. Gedicke, and D. Rim, Explicit error estimates for Courant, Crouzeix–Raviart and Raviart–Thomas finite element methods, J. Comput. Math. 30 (2012), No. 4, 337–353.10.4208/jcm.1108-m3677Search in Google Scholar

[18] C. Carstensen, B. Gräßle, and N. Nataraj, A posteriori error control for fourth-order semilinear problems with quadratic nonlinearity, SIAM J. Numer. Anal. 62 (2024) (to appear).Search in Google Scholar

[19] C. Carstensen and J. Hu, A unifying theory of a posteriori error control for nonconforming finite element methods, Numer. Math. 107 (2007), No. 3, 473–502.10.1007/s00211-007-0068-zSearch in Google Scholar

[20] C. Carstensen, J. Hu, and A. Orlando, Framework For The A Posteriori Error Analysis Of Nonconforming Finite Element, SIAM J. Nu-mer. Anal. 45 (2007), No. 1, 68–82.10.1137/050628854Search in Google Scholar

[21] C. Carstensen, G. Mallik, and N. Nataraj, A priori and a posteriori error control of discontinuous Galerkin finite element methods for the von Kármán equations, IMA J. Numer. Anal. 39 (2019), 167–200.10.1093/imanum/dry003Search in Google Scholar

[22] C. Carstensen and C. Merdon, Computational survey on a posteriori error estimators for the Crouzeix–Raviart nonconforming finite element method for the Stokes problem, Comput. Methods Appl. Math. 14 (2014), No. 1, 35–54.10.1515/cmam-2013-0021Search in Google Scholar

[23] C. Carstensen and N. Nataraj, A priori and a posteriori error analysis of the Crouzeix–Raviart and Morley FEM with original and modified right-hand sides, Comput. Methods Appl. Math. 21 (2021), No. 2, 289–315.10.1515/cmam-2021-0029Search in Google Scholar

[24] C. Carstensen and N. Nataraj, Lowest-order equivalent nonstandard finite element methods for biharmonic plates, ESAIM: M2AN 56 (2022), No. 1, 41–78.10.1051/m2an/2021085Search in Google Scholar

[25] C. Carstensen and S. Puttkammer, How to prove the discrete reliability for nonconforming finite element methods, J. Comput. Math 38 (2020), No. 1, 142–175.10.4208/jcm.1908-m2018-0174Search in Google Scholar

[26] C. Carstensen and S. Puttkammer, Direct guaranteed lower eigenvalue bounds with optimal a priori convergence rates for the bi-Laplacian, SIAM Journal on Numerical Analysis, 61 (2023), No. 2, 812–836.Search in Google Scholar

[27] P. G. Ciarlet, The Finite Element Method for Elliptic Problems, North-Holland, Amsterdam, 1978.10.1115/1.3424474Search in Google Scholar

[28] P. Clément, Approximation by finite element functions using local regularization, Rev. Française Automat. Informat. Recherche Opéra-tionnelle Sér. Rouge Anal. Numér. 9 (1975), No. R-2, 77–84.10.1051/m2an/197509R200771Search in Google Scholar

[29] E. Dari, R. Duran, C. Padra, and V. Vampa, A posteriori error estimators for nonconforming finite element methods, RAIRO Modél. Math. Anal. Numér. 30 (1996), No. 4, 385–400.10.1051/m2an/1996300403851Search in Google Scholar

[30] A. Ern and J. L. Guermond, Finite Elements, I: Approximation and Interpolation, Texts in Applied Mathematicsi, Vol. 72, Springer International Publishing, Cham, 2021.10.1007/978-3-030-56341-7Search in Google Scholar

[31] L. C. Evans, Partial Differential Equations, 2nd ed., Graduate Studies in Mathematics, Vol. 19, American Mathematical Society, Providence, R.I, 2010.Search in Google Scholar

[32] X. Feng and O. A. Karakashian, Fully discrete dynamic mesh discontinuous Galerkin methods for the Cahn–Hilliard equation of phase transition, Math. Comp. 76 (2007), No. 259, 1093–1117.10.1090/S0025-5718-07-01985-0Search in Google Scholar

[33] D. Gallistl, Morley finite element method for the eigenvalues of the biharmonic operator, IMA J. Numer. Anal. 35 (2015), No. 4, 1779– 1811.10.1093/imanum/dru054Search in Google Scholar

[34] E. H. Georgoulis, P. Houston, and J. Virtanen, An a posteriori error indicator for discontinuous Galerkin approximations of fourth-order elliptic problems, IMA J. Numer. Anal. 31 (2011), No. 1, 281–298.10.1093/imanum/drp023Search in Google Scholar

[35] J. Hu, Z. Shi, and J. Xu, Convergence and optimality of the adaptive Morley element method, Numer. Math. 121 (2012), No. 4, 731–752.10.1007/s00211-012-0445-0Search in Google Scholar

[36] J. Hu and Z. C. Shi, A new a posteriori error estimate for the Morley element, Numer. Math. 112 (2009), No. 1, 25–40.10.1007/s00211-008-0205-3Search in Google Scholar

[37] D. Kim, A. K. Pani, and E.-J. Park, Morley finite element methods for the stationary quasi-geostrophic equation, Computer Methods in Applied Mechanics and Engineering 375 (2021), 113639.10.1016/j.cma.2020.113639Search in Google Scholar

[38] L. R. Scott, C Piecewise Polynomials Satisfying Boundary Conditions, Research Report UC/CS TR-2019-18, Dept. Comp. Sci., Univ. Chicago, Report, 2019.Search in Google Scholar

[39] L. R. Scott and S. Zhang, Finite element interpolation of nonsmooth functions satisfying boundary conditions, Math. Comp. 54 (1990), No. 190, 483–493.10.1090/S0025-5718-1990-1011446-7Search in Google Scholar

[40] E. Süli and I. Mozolevski, hp-version interior penalty DGFEMs for the biharmonic equation, Comput. Methods Appl. Mech. Engrg. 196 (2007), No. 13-16, 1851–1863.10.1016/j.cma.2006.06.014Search in Google Scholar

[41] R. Vanselow, New results concerning the DWR method for some nonconforming FEM, Appl. Math. 57 (2012), No. 6, 551–568.10.1007/s10492-012-0033-8Search in Google Scholar

[42] A. Veeser and P. Zanotti, Quasi-optimal nonconforming methods for symmetric elliptic problems, I—Abstract theory, SIAM J. Nu-mer. Anal. 56 (2018), No. 3, 1621–1642.10.1137/17M1116362Search in Google Scholar

[43] A. Veeser and P. Zanotti, Quasi-optimal nonconforming methods for symmetric elliptic problems, III—Discontinuous Galerkin and other interior penalty methods, SIAM J. Numer. Anal. 56 (2018), No. 5, 2871–2894.10.1137/17M1151675Search in Google Scholar

[44] A. Veeser and P. Zanotti, Quasi-optimal nonconforming methods for symmetric elliptic problems, II—Overconsistency and classical nonconforming elements, SIAM J. Numer. Anal. 57 (2019), No. 1, 266–292.10.1137/17M1151651Search in Google Scholar

[45] R. Verfürth, A Review of a Posteriori Error Estimation and Adaptive Mesh-Refinement Techniques, Wiley-Teubner, 1996.Search in Google Scholar

A A posteriori error control of a piecewise polynomial source in H–2(Ω)

This appendix provides an alternative view on the reliable and efficient estimator from Section 7 as lower and upper bounds for the dual norm of a piecewise polynomial source in H−2 (Ω). Suppose the piecewise polynomials Λ0Pk (T), Λ1Pk (T;R2), and Λ2Pk (T; S) define the linear functional Λ ∈ H−2 (Ω) by

(A.1) Λ(v):=ΩΛ0v+Λ1v+Λ2:D2vdxvH02(Ω).

Recall the transfer operators IM, Ih , Jh for the five quadratic discretization schemes of Section 4 listed in Table 1. A reliable and efficient estimator μ2 (T) := μ12(T)+μ22(T)+μ32(T) of the functional Λ is given by

μ12(T):=hT2Λ0divpwΛ1+divpw2Λ22μ22(T):=EE(Ω)hE3Λ1divpwΛ2Λ2τE/sEvEL2(E)2μ32(T):=EE(Ω)hE1ΠE,0Λ2vEEvEL2(E)2 if Ih=idhEΛ2vEEvEL2(E)2 if Ih=IC.

Theorem A.1 (reliability and efficiency). There exist positive constants Crel , Ceff > 0 that exclusively depend on the shape regularity of T and on the polynomial degree k ∈ N0 such that

Crel1|Λ1JhIhIM|μ(T)Ceff|Λ|.

Proof. The discussion in Subsection 7.2 applies to Res := F := Λ and uh := 0 with apx(F, T) = 0. In this particular case, Proposition 7.1 provides the first inequality

Λ1JhIhIMCrelμ(T)

with Crel = C4. Let uH02 (Ω) denote the Riesz representation of a(u, ·) = Λ ∈ H−2 (Ω) with the isometry |||Λ|||* = |||u||| in the Hilbert space H02 ((Ω), a) and ||| · ||| ≡ a(·, ·)1/2. Then the efficiency estimate

μ(T)Ceff |Λ|

follows from Proposition 7.2 with Ceff = C8.

Theorem A.1 allows for a direct application to the linearization of semilinear problems in [18]. It can be further generalized in various directions, e.g., in the spirit of Section 7 that considers the a posteriori error analysis of the linear biharmonic problem for a more general class of functionals in H−2 (Ω) including line and point loads. The reliability requires only piecewise smoothness of Λ0, Λ1, Λ2 so that the traces and derivatives in μ1, μ2, μ3 exist, while the efficiency may require extra oscillation terms (as in (7.5)).

Received: 2022-09-24
Revised: 2023-01-13
Accepted: 2023-01-23
Published Online: 2024-03-09
Published in Print: 2024-03-25

© 2024 Walter de Gruyter GmbH, Berlin/Boston

Downloaded on 28.4.2024 from https://www.degruyter.com/document/doi/10.1515/jnma-2022-0092/html
Scroll to top button