Abstract
We consider a mixed finite element method for an obstacle problem with the p-Laplace differential operator for , where the obstacle condition is imposed by using a Lagrange multiplier. In the discrete setting the Lagrange multiplier basis forms a biorthogonal system with the standard finite element basis so that the variational inequality can be realized in the point-wise form. We provide a general a posteriori error estimate for adaptivity and prove an a priori error estimate. We present numerical results for the adaptive scheme (mesh-size adaptivity with and without polynomial degree adaptation) for the singular case and the degenerated case . We also present numerical results on the mesh independency and on the polynomial degree scaling of the discrete inf-sup constant when using biorthogonal basis functions for the dual variable defined on the same mesh with the same polynomial degree distribution.
Funding statement: The visit of the third author to the University of Newcastle, Australia was partially supported by the priority research centre of the University of Newcastle for Computer-Assisted Research Mathematics and its Applications.
A Basic Results
For the convenience of the reader, let us recall some, for the p-Laplace fundamental, inequalities. We refer to [1] for the proofs of these results.
Lemma 13 (cf. [23, Lemmas 5.1–5.3]).
For all , , and ξ, there holds
For all a, , , , there holds
where
For all a, , , , and ,
where β is such that .
Lemma 14.
Let a, and . Then there holds
and for ,
with if and if .
Lemma 15 ([15, Proposition 2.1]).
Let .
It holds , and, when , if and only if .
There holds for any , .
Furthermore, for , there holds
(A.1)For , , , there holds
(A.2)
The constant in (A.1) and (A.2) can be stated explicitly and on the right-hand side of (A.1) can be eliminated.
Lemma 16.
For there holds
and for there holds
for all v, .
Corollary 17.
For there exists a constant such that
for all v, .
Lemma 18.
For there exists a constant such that
and for there exists a constant such that
for all , , v, .
Acknowledgements
Ernst P. Stephan expresses his sincere thanks to Bishnu Lamichhane for his hospitality during the visit.
References
[1] L. Banz, B. P. Lamichhane and E. P. Stephan, Higher order FEM for the obstacle problem of the p-Laplacian – A variational inequality approach, preprint, (2017). 10.1016/j.camwa.2018.07.016Search in Google Scholar
[2] L. Banz and A. Schröder, Biorthogonal basis functions in hp-adaptive FEM for elliptic obstacle problems, Comput. Math. Appl. 70 (2015), no. 8, 1721–1742. 10.1016/j.camwa.2015.07.010Search in Google Scholar
[3] L. Banz and E. P. Stephan, A posteriori error estimates of hp-adaptive IPDG-FEM for elliptic obstacle problems, Appl. Numer. Math. 76 (2014), 76–92. 10.1016/j.apnum.2013.10.004Search in Google Scholar
[4] L. Banz and E. P. Stephan, hp-adaptive IPDG/TDG-FEM for parabolic obstacle problems, Comput. Math. Appl. 67 (2014), 712–731. 10.1016/j.camwa.2013.03.003Search in Google Scholar
[5] J. W. Barrett and W. Liu, Quasi-norm error bounds for the finite element approximation of a non-Newtonian flow, Numer. Math. 68 (1994), no. 4, 437–456. 10.1007/s002110050071Search in Google Scholar
[6] S. Bartels and C. Carstensen, Averaging techniques yield reliable a posteriori finite element error control for obstacle problems, Numer. Math. 99 (2004), 225–249. 10.1007/s00211-004-0553-6Search in Google Scholar
[7] D. Braess, A posteriori error estimators for obstacle problems–another look, Numer. Math. 101 (2005), no. 3, 415–421. 10.1007/s00211-005-0634-1Search in Google Scholar
[8] D. Braess, C. Carstensen and R. Hoppe, Convergence analysis of a conforming adaptive finite element method for an obstacle problem, Numer. Math. 107 (2007), 455–471. 10.1007/s00211-007-0098-6Search in Google Scholar
[9] S. C. Brenner and L. R. Scott, The Mathematical Theory of Finite Element Methods, 3 ed., Springer, New York, 2008. 10.1007/978-0-387-75934-0Search in Google Scholar
[10] C. Carstensen and J. Hu, An optimal adaptive finite element method for an obstacle problem, Comput. Methods Appl. Math. 15 (2015), 259–277. 10.1515/cmam-2015-0017Search in Google Scholar
[11] C. Carstensen and R. Klose, A posteriori finite element error control for the p-Laplace problem, SIAM J. Sci. Comput. 25 (2003), no. 3, 792–814. 10.1137/S1064827502416617Search in Google Scholar
[12] C. Carstensen, W. Liu and N. Yan, A posteriori FE error control for p-Laplacian by gradient recovery in quasi-norm, Math. Comp. 75 (2006), no. 256, 1599–1616. 10.1090/S0025-5718-06-01819-9Search in Google Scholar
[13] M. Chipot, Elliptic Equations: An Introductory Course, Birkhäuser, Basel, 2009. 10.1007/978-3-7643-9982-5Search in Google Scholar
[14] L. Diening and C. Kreuzer, Linear convergence of an adaptive finite element method for the p-Laplacian equation, SIAM J. Numer. Anal. 46 (2008), no. 2, 614–638. 10.1137/070681508Search in Google Scholar
[15] C. Ebmeyer and W. Liu, Quasi-norm interpolation error estimates for the piecewise linear finite element approximation of p-Laplacian problems, Numer. Math. 100 (2005), no. 2, 233–258. 10.1007/s00211-005-0594-5Search in Google Scholar
[16] J. Gwinner, hp-FEM convergence for unilateral contact problems with Tresca friction in plane linear elastostatics, J. Comput. Appl. Math. 254 (2013), 175–184. 10.1016/j.cam.2013.03.013Search in Google Scholar
[17] S. Hüeber and B. Wohlmuth, A primal–dual active set strategy for non-linear multibody contact problems, Comput. Methods Appl. Mech. Engrg. 194 (2005), no. 27–29, 3147–3166. 10.1016/j.cma.2004.08.006Search in Google Scholar
[18] G. Jouvet and E. Bueler, Steady, shallow ice sheets as obstacle problems: Well-posedness and finite element approximation, SIAM J. Appl. Math. 72 (2012), no. 4, 1292–1314. 10.1137/110856654Search in Google Scholar
[19] R. Krause, B. Müller and G. Starke, An adaptive least-squares mixed finite element method for the signorini problem, Numer. Methods Partial Differential Equations 33 (2017), 276–289. 10.1002/num.22086Search in Google Scholar
[20] A. Krebs and E. Stephan, A p-version finite element method for nonlinear elliptic variational inequalities in 2D, Numer. Math. 105 (2007), no. 3, 457–480. 10.1007/s00211-006-0035-0Search in Google Scholar
[21] B. Lamichhane and B. Wohlmuth, Biorthogonal bases with local support and approximation properties, Math. Comp. 76 (2007), no. 257, 233–249. 10.1090/S0025-5718-06-01907-7Search in Google Scholar
[22] M. Lewicka and J. J. Manfredi, The obstacle problem for the p-Laplacian via optimal stopping of tug-of-war games, Probab. Theory Related Fields 167 (2017), no. 1–2, 349–378. 10.1007/s00440-015-0684-ySearch in Google Scholar
[23] W. Liu and N. Yan, Quasi-norm local error estimators for p-Laplacian, SIAM J. Numer. Anal. 39 (2001), no. 1, 100–127. 10.1137/S0036142999351613Search in Google Scholar
[24] W. Liu and N. Yan, On quasi-norm interpolation error estimation and a posteriori error estimates for p-Laplacian, SIAM J. Numer. Anal. 40 (2002), no. 5, 1870–1895. 10.1137/S0036142901393589Search in Google Scholar
[25] M. Maischak and E. P. Stephan, Adaptive hp-versions of BEM for Signorini problems, Appl. Numer. Math. 54 (2005), no. 3, 425–449. 10.1016/j.apnum.2004.09.012Search in Google Scholar
[26] D. Malkus, Eigenproblems associated with the discrete LBB condition for incompressible finite elements, Internat. J. Engrg. Sci. 19 (1981), no. 10, 1299–1310. 10.1016/0020-7225(81)90013-6Search in Google Scholar
[27] J. M. Melenk, hp-interpolation of nonsmooth functions and an application to hp-a posteriori error estimation, SIAM J. Numer. Anal. 43 (2005), no. 1, 127–155. 10.1137/S0036142903432930Search in Google Scholar
[28] N. Ovcharova and L. Banz, Coupling regularization and adaptive hp-BEM for the solution of a delamination problem, Numer. Math. 137 (2017), 303–337. 10.1007/s00211-017-0879-5Search in Google Scholar
[29] J. Qin, On the convergence of some low order mixed finite elements for incompressible fluids, Ph.D. thesis, The Pennsylvania State University, 1994. Search in Google Scholar
[30] A. Veeser, Efficient and reliable a posteriori error estimators for elliptic obstacle problems, SIAM J. Numer. Anal. 39 (2001), no. 1, 146–167. 10.1137/S0036142900370812Search in Google Scholar
[31] T. Wick, An error-oriented Newton/inexact augmented Lagrangian approach for fully monolithic phase-field fracture propagation, SIAM J. Sci. Comput. 39 (2017), no. 4, B589–B617. 10.1137/16M1063873Search in Google Scholar
© 2019 Walter de Gruyter GmbH, Berlin/Boston
Abstract
We consider a mixed finite element method for an obstacle problem with the p-Laplace differential operator for , where the obstacle condition is imposed by using a Lagrange multiplier. In the discrete setting the Lagrange multiplier basis forms a biorthogonal system with the standard finite element basis so that the variational inequality can be realized in the point-wise form. We provide a general a posteriori error estimate for adaptivity and prove an a priori error estimate. We present numerical results for the adaptive scheme (mesh-size adaptivity with and without polynomial degree adaptation) for the singular case and the degenerated case . We also present numerical results on the mesh independency and on the polynomial degree scaling of the discrete inf-sup constant when using biorthogonal basis functions for the dual variable defined on the same mesh with the same polynomial degree distribution.
Funding statement: The visit of the third author to the University of Newcastle, Australia was partially supported by the priority research centre of the University of Newcastle for Computer-Assisted Research Mathematics and its Applications.
A Basic Results
For the convenience of the reader, let us recall some, for the p-Laplace fundamental, inequalities. We refer to [1] for the proofs of these results.
Lemma 13 (cf. [23, Lemmas 5.1–5.3]).
For all , , and ξ, there holds
For all a, , , , there holds
where
For all a, , , , and ,
where β is such that .
Lemma 14.
Let a, and . Then there holds
and for ,
with if and if .
Lemma 15 ([15, Proposition 2.1]).
Let .
It holds , and, when , if and only if .
There holds for any , .
Furthermore, for , there holds
(A.1)For , , , there holds
(A.2)
The constant in (A.1) and (A.2) can be stated explicitly and on the right-hand side of (A.1) can be eliminated.
Lemma 16.
For there holds
and for there holds
for all v, .
Corollary 17.
For there exists a constant such that
for all v, .
Lemma 18.
For there exists a constant such that
and for there exists a constant such that
for all , , v, .
Acknowledgements
Ernst P. Stephan expresses his sincere thanks to Bishnu Lamichhane for his hospitality during the visit.
References
[1] L. Banz, B. P. Lamichhane and E. P. Stephan, Higher order FEM for the obstacle problem of the p-Laplacian – A variational inequality approach, preprint, (2017). 10.1016/j.camwa.2018.07.016Search in Google Scholar
[2] L. Banz and A. Schröder, Biorthogonal basis functions in hp-adaptive FEM for elliptic obstacle problems, Comput. Math. Appl. 70 (2015), no. 8, 1721–1742. 10.1016/j.camwa.2015.07.010Search in Google Scholar
[3] L. Banz and E. P. Stephan, A posteriori error estimates of hp-adaptive IPDG-FEM for elliptic obstacle problems, Appl. Numer. Math. 76 (2014), 76–92. 10.1016/j.apnum.2013.10.004Search in Google Scholar
[4] L. Banz and E. P. Stephan, hp-adaptive IPDG/TDG-FEM for parabolic obstacle problems, Comput. Math. Appl. 67 (2014), 712–731. 10.1016/j.camwa.2013.03.003Search in Google Scholar
[5] J. W. Barrett and W. Liu, Quasi-norm error bounds for the finite element approximation of a non-Newtonian flow, Numer. Math. 68 (1994), no. 4, 437–456. 10.1007/s002110050071Search in Google Scholar
[6] S. Bartels and C. Carstensen, Averaging techniques yield reliable a posteriori finite element error control for obstacle problems, Numer. Math. 99 (2004), 225–249. 10.1007/s00211-004-0553-6Search in Google Scholar
[7] D. Braess, A posteriori error estimators for obstacle problems–another look, Numer. Math. 101 (2005), no. 3, 415–421. 10.1007/s00211-005-0634-1Search in Google Scholar
[8] D. Braess, C. Carstensen and R. Hoppe, Convergence analysis of a conforming adaptive finite element method for an obstacle problem, Numer. Math. 107 (2007), 455–471. 10.1007/s00211-007-0098-6Search in Google Scholar
[9] S. C. Brenner and L. R. Scott, The Mathematical Theory of Finite Element Methods, 3 ed., Springer, New York, 2008. 10.1007/978-0-387-75934-0Search in Google Scholar
[10] C. Carstensen and J. Hu, An optimal adaptive finite element method for an obstacle problem, Comput. Methods Appl. Math. 15 (2015), 259–277. 10.1515/cmam-2015-0017Search in Google Scholar
[11] C. Carstensen and R. Klose, A posteriori finite element error control for the p-Laplace problem, SIAM J. Sci. Comput. 25 (2003), no. 3, 792–814. 10.1137/S1064827502416617Search in Google Scholar
[12] C. Carstensen, W. Liu and N. Yan, A posteriori FE error control for p-Laplacian by gradient recovery in quasi-norm, Math. Comp. 75 (2006), no. 256, 1599–1616. 10.1090/S0025-5718-06-01819-9Search in Google Scholar
[13] M. Chipot, Elliptic Equations: An Introductory Course, Birkhäuser, Basel, 2009. 10.1007/978-3-7643-9982-5Search in Google Scholar
[14] L. Diening and C. Kreuzer, Linear convergence of an adaptive finite element method for the p-Laplacian equation, SIAM J. Numer. Anal. 46 (2008), no. 2, 614–638. 10.1137/070681508Search in Google Scholar
[15] C. Ebmeyer and W. Liu, Quasi-norm interpolation error estimates for the piecewise linear finite element approximation of p-Laplacian problems, Numer. Math. 100 (2005), no. 2, 233–258. 10.1007/s00211-005-0594-5Search in Google Scholar
[16] J. Gwinner, hp-FEM convergence for unilateral contact problems with Tresca friction in plane linear elastostatics, J. Comput. Appl. Math. 254 (2013), 175–184. 10.1016/j.cam.2013.03.013Search in Google Scholar
[17] S. Hüeber and B. Wohlmuth, A primal–dual active set strategy for non-linear multibody contact problems, Comput. Methods Appl. Mech. Engrg. 194 (2005), no. 27–29, 3147–3166. 10.1016/j.cma.2004.08.006Search in Google Scholar
[18] G. Jouvet and E. Bueler, Steady, shallow ice sheets as obstacle problems: Well-posedness and finite element approximation, SIAM J. Appl. Math. 72 (2012), no. 4, 1292–1314. 10.1137/110856654Search in Google Scholar
[19] R. Krause, B. Müller and G. Starke, An adaptive least-squares mixed finite element method for the signorini problem, Numer. Methods Partial Differential Equations 33 (2017), 276–289. 10.1002/num.22086Search in Google Scholar
[20] A. Krebs and E. Stephan, A p-version finite element method for nonlinear elliptic variational inequalities in 2D, Numer. Math. 105 (2007), no. 3, 457–480. 10.1007/s00211-006-0035-0Search in Google Scholar
[21] B. Lamichhane and B. Wohlmuth, Biorthogonal bases with local support and approximation properties, Math. Comp. 76 (2007), no. 257, 233–249. 10.1090/S0025-5718-06-01907-7Search in Google Scholar
[22] M. Lewicka and J. J. Manfredi, The obstacle problem for the p-Laplacian via optimal stopping of tug-of-war games, Probab. Theory Related Fields 167 (2017), no. 1–2, 349–378. 10.1007/s00440-015-0684-ySearch in Google Scholar
[23] W. Liu and N. Yan, Quasi-norm local error estimators for p-Laplacian, SIAM J. Numer. Anal. 39 (2001), no. 1, 100–127. 10.1137/S0036142999351613Search in Google Scholar
[24] W. Liu and N. Yan, On quasi-norm interpolation error estimation and a posteriori error estimates for p-Laplacian, SIAM J. Numer. Anal. 40 (2002), no. 5, 1870–1895. 10.1137/S0036142901393589Search in Google Scholar
[25] M. Maischak and E. P. Stephan, Adaptive hp-versions of BEM for Signorini problems, Appl. Numer. Math. 54 (2005), no. 3, 425–449. 10.1016/j.apnum.2004.09.012Search in Google Scholar
[26] D. Malkus, Eigenproblems associated with the discrete LBB condition for incompressible finite elements, Internat. J. Engrg. Sci. 19 (1981), no. 10, 1299–1310. 10.1016/0020-7225(81)90013-6Search in Google Scholar
[27] J. M. Melenk, hp-interpolation of nonsmooth functions and an application to hp-a posteriori error estimation, SIAM J. Numer. Anal. 43 (2005), no. 1, 127–155. 10.1137/S0036142903432930Search in Google Scholar
[28] N. Ovcharova and L. Banz, Coupling regularization and adaptive hp-BEM for the solution of a delamination problem, Numer. Math. 137 (2017), 303–337. 10.1007/s00211-017-0879-5Search in Google Scholar
[29] J. Qin, On the convergence of some low order mixed finite elements for incompressible fluids, Ph.D. thesis, The Pennsylvania State University, 1994. Search in Google Scholar
[30] A. Veeser, Efficient and reliable a posteriori error estimators for elliptic obstacle problems, SIAM J. Numer. Anal. 39 (2001), no. 1, 146–167. 10.1137/S0036142900370812Search in Google Scholar
[31] T. Wick, An error-oriented Newton/inexact augmented Lagrangian approach for fully monolithic phase-field fracture propagation, SIAM J. Sci. Comput. 39 (2017), no. 4, B589–B617. 10.1137/16M1063873Search in Google Scholar
© 2019 Walter de Gruyter GmbH, Berlin/Boston
Articles in the same Issue
- Frontmatter
- Presentation of the Special Issue on Recent Advances in PDE: Theory, Computations and Applications
- Higher Order Mixed FEM for the Obstacle Problem of the p-Laplace Equation Using Biorthogonal Systems
- Edge Patch-Wise Local Projection Stabilized Nonconforming FEM for the Oseen Problem
- Multigrid Methods Based on Hodge Decomposition for a Quad-Curl Problem
- Quasi-Optimality of Adaptive Mixed FEMs for Non-selfadjoint Indefinite Second-Order Linear Elliptic Problems
- Analysis of FEAST Spectral Approximations Using the DPG Discretization
- Stabilizability of Infinite Dimensional Systems by Finite Dimensional Control
- An Explicit-Implicit Splitting Method for a Convection-Diffusion Problem
- Positivity Preserving Gradient Approximation with Linear Finite Elements
- An Elementary Method of Deriving A Posteriori Error Equalities and Estimates for Linear Partial Differential Equations
- Estimation of the Time-Dependent Body Force Needed to Exert on a Membrane to Reach a Desired State at the Final Time
- Mixed Schemes for Fourth-Order DIV Equations
- Improved Analysis and Simulation of a Time-Domain Carpet Cloak Model
- The Boundary Effect in the Accuracy Estimate for the Grid Solution of the Fractional Differential Equation
Articles in the same Issue
- Frontmatter
- Presentation of the Special Issue on Recent Advances in PDE: Theory, Computations and Applications
- Higher Order Mixed FEM for the Obstacle Problem of the p-Laplace Equation Using Biorthogonal Systems
- Edge Patch-Wise Local Projection Stabilized Nonconforming FEM for the Oseen Problem
- Multigrid Methods Based on Hodge Decomposition for a Quad-Curl Problem
- Quasi-Optimality of Adaptive Mixed FEMs for Non-selfadjoint Indefinite Second-Order Linear Elliptic Problems
- Analysis of FEAST Spectral Approximations Using the DPG Discretization
- Stabilizability of Infinite Dimensional Systems by Finite Dimensional Control
- An Explicit-Implicit Splitting Method for a Convection-Diffusion Problem
- Positivity Preserving Gradient Approximation with Linear Finite Elements
- An Elementary Method of Deriving A Posteriori Error Equalities and Estimates for Linear Partial Differential Equations
- Estimation of the Time-Dependent Body Force Needed to Exert on a Membrane to Reach a Desired State at the Final Time
- Mixed Schemes for Fourth-Order DIV Equations
- Improved Analysis and Simulation of a Time-Domain Carpet Cloak Model
- The Boundary Effect in the Accuracy Estimate for the Grid Solution of the Fractional Differential Equation