Skip to content
Licensed Unlicensed Requires Authentication Published by De Gruyter October 18, 2021

Two mixed finite element formulations for the weak imposition of the Neumann boundary conditions for the Darcy flow

  • Erik Burman and Riccardo Puppi EMAIL logo

Abstract

We propose two different discrete formulations for the weak imposition of the Neumann boundary conditions of the Darcy flow. The Raviart–Thomas mixed finite element on both triangular and quadrilateral meshes is considered for both methods. One is a consistent discretization depending on a weighting parameter scaling as 𝒪(h−1), while the other is a penalty-type formulation obtained as the discretization of a perturbation of the original problem and relies on a parameter scaling as 𝒪(hk−1), k being the order of the Raviart–Thomas space. We rigorously prove that both methods are stable and result in optimal convergent numerical schemes with respect to appropriate mesh-dependent norms, although the chosen norms do not scale as the usual L2-norm. However, we are still able to recover the optimal a priori L2-error estimates for the velocity field, respectively, for high-order and the lowest-order Raviart–Thomas discretizations, for the first and second numerical schemes. Finally, some numerical examples validating the theory are exhibited.

MSC 2010: 65M60

Funding statement: Erik Burman was partially supported by the EPSRC grants EP/P01576X/1 and EP/T033126/1. Riccardo Puppi was partially supported by ERC AdG project CHANGE No. 694515.

References

[1] R. Adams and J. Fournier, Sobolev Spaces, 2nd ed., Pure and Applied Mathematics, Academic Press, New York, 2003.Search in Google Scholar

[2] D. N. Arnold, D. Bofl, and R. S. Falk, Quadrilateral H(div) finite elements, SIAM J. Numer. Anal. 42 (2004), No. 6, 2429–2451.10.1137/S0036142903431924Search in Google Scholar

[3] I. Babuška, The finite element method with Lagrangian multipliers, Numer. Math. 20 (1973), No. 3, 179–192.10.1007/BF01436561Search in Google Scholar

[4] I. Babuška, The finite element method with penalty, Math. Comput. 27 (1973), No. 122, 221–228.10.1090/S0025-5718-1973-0351118-5Search in Google Scholar

[5] H. J. C. Barbosa and T. J. R. Hughes, The finite element method with Lagrange multipliers on the boundary: circumventing the Babuška–Brezzi condition, Comp. Methods Appl. Mech. Engrg. 85 (1991), No. 1, 109–128.10.1016/0045-7825(91)90125-PSearch in Google Scholar

[6] D. Boffi, F. Brezzi, and M. Fortin, Mixed Finite Element Methods and Applications, Springer Series in Computational Mathematics, 2013.10.1007/978-3-642-36519-5Search in Google Scholar

[7] E. Burman, S. Claus, P. Hansbo, M. G. Larson, and A. Massing, CutFEM: Discretizing geometry and partial differential equations, Int. J. Numer. Meth. Engrg. 104 (2015), No. 7, 472–501.10.1002/nme.4823Search in Google Scholar

[8] E. Burman, M. G. Larson, and L. Oksanen, Primal–dual mixed finite element methods for the elliptic Cauchy problem, SIAM J. Numer. Anal. 56 (2018), No. 6, 3480–3509.10.1137/17M1163335Search in Google Scholar

[9] C. D’Angelo and A. Scotti, A mixed finite element method for Darcy flow in fractured porous media with non-matching grids, ESAIM: Math. Modelling Numer. Anal. 46 (2011), No. 2, 465–489.10.1051/m2an/2011148Search in Google Scholar

[10] A. Ern and J.-L. Guermond, Evaluation of the condition number in linear systems arising in finite element approximations, ESAIM: Math. Modelling Numer. Anal. 40 (2006), No. 1, 29–48.10.1051/m2an:2006006Search in Google Scholar

[11] P. Grisvard, Elliptic Problems in Nonsmooth Domains, Society for Industrial and Applied Mathematics, 2011.10.1137/1.9781611972030Search in Google Scholar

[12] C. Johnson, Numerical Solution of Partial Differential Equations by the Finite Element Method, Cambridge Cambridge University Press, 1987.Search in Google Scholar

[13] J. Könnö, D. Schötzau, and R. Stenberg, Mixed Finite Element Methods for Problmes with Robin Boundary Conditions, SIAM J. Numer. Anal. 49 (2011), No. 1/2, 285–308.10.1137/09077970XSearch in Google Scholar

[14] P. Monk, Finite Element Metods for Maxwell’s Equations, Oxford University Press, 2003.10.1093/acprof:oso/9780198508885.001.0001Search in Google Scholar

[15] J. Nitsche, Über ein Variationsprinzip zur Lösung von Dirichlet-Problemen bei Verwendung von Teilräumen, die keinen Randbedingungen unterworfen sind, Abhandlungen aus dem Mathematischen Seminar der Universität Hamburg 36 (1971), No. 1, 9–15.10.1007/BF02995904Search in Google Scholar

[16] A. Quarteroni and A. Valli, Numerical Approximation of Partial Differential Equations, Springer, Berlin–Heidelberg, 1994.10.1007/978-3-540-85268-1Search in Google Scholar

[17] R. Stenberg, On some techniques for approximating boundary conditions in the finite element method, J. Comp. Appl. Math. 63 (1995), No. 1, 139–148.10.1016/0377-0427(95)00057-7Search in Google Scholar

Received: 2021-04-03
Revised: 2021-06-08
Accepted: 2021-06-18
Published Online: 2021-10-18
Published in Print: 2022-06-27

© 2022 Walter de Gruyter GmbH, Berlin/Boston

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