Analysis of a stabilized penalty-free Nitsche method for the Brinkman, Stokes, and Darcy problems

Blank, Laura; Caiazzo, Alfonso; Chouly, Franz; Lozinski, Alexei; Mura, Joaquin

ESAIM: Mathematical Modelling and Numerical Analysis , Tome 52 (2018) no. 6, pp. 2149-2185.

In this paper we study the Brinkman model as a unified framework to allow the transition between the Darcy and the Stokes problems. We propose an unconditionally stable low-order finite element approach, which is robust with respect to the whole range of physical parameters, and is based on the combination of stabilized equal-order finite elements with a non-symmetric penalty-free Nitsche method for the weak imposition of essential boundary conditions. In particular, we study the properties of the penalty-free Nitsche formulation for the Brinkman setting, extending a recently reported analysis for the case of incompressible elasticity (Boiveau and Burman, IMA J. Numer. Anal. 36 (2016) 770-795). Focusing on the two-dimensional case, we obtain optimal a priori error estimates in a mesh-dependent norm, which, converging to natural norms in the cases of Stokes or Darcy ows, allows to extend the results also to these limits. Moreover, we show that, in order to obtain robust estimates also in the Darcy limit, the formulation shall be equipped with a Grad-Div stabilization and an additional stabilization to control the discontinuities of the normal velocity along the boundary. The conclusions of the analysis are supported by numerical simulations.

[1] R.A. Adams, Sobolev Spaces. In Vol. 65 of Pure and Applied Mathematics. Academic Press [A subsidiary of Harcourt Brace Jovanovich, Publishers], New York-London (1975). | MR | Zbl

[2] I. Babuška, The finite element method with Lagrangian multipliers. Numer. Math. 20 (1972/1973) 179–192. | DOI | MR | Zbl

[3] S. Badia and R. Codina, Unified stabilized finite element formulations for the Stokes and the Darcy problems. SIAM J. Numer. Anal. 47 (2009) 1971–2000. | DOI | MR | Zbl

[4] T. Boiveau, Penalty-free Nitsche method for interface problems in computational mechanics, Ph.D. thesis, University College London (2016).

[5] T. Boiveau and E. Burman, A penalty-free Nitsche method for the weak imposition of boundary conditions in compressible and incompressible elasticity. IMA J. Numer. Anal. 36 (2016) 770–795. | DOI | MR | Zbl

[6] S.C. Brenner and L. Ridgway Scott, The mathematical theory of finite element methods, 3rd edition. In Vol. 15 of Texts in Applied Mathematics. Springer, New York, NY (2008). | MR | Zbl

[7] F. Brezzi, On the existence, uniqueness and approximation of saddle-point problems arising from Lagrangian multipliers. Rev. Française Automat. Informat. Recherche Opérationnelle Sér. 8 (1974) 129–151. | Numdam | MR | Zbl

[8] F. Brezzi and Pitkäranta J., On the stabilization of finite element approximations of the Stokes equations. In: Efficient Solutions of Elliptic Systems (Kiel, 1984). In Vol. 10 of Notes Numer. Fluid Mech. Friedr. Vieweg, Braunschweig (1984) 11–19. | DOI | MR | Zbl

[9] F. Brezzi and M. Fortin, Mixed and hybrid finite element methods. In Vol. 15 of Springer Series in Computational Mathematics. Springer-Verlag, New York, NY (1991). | DOI | MR | Zbl

[10] H.C. Brinkman, A calculation of the viscous force exerted by a flowing fluid on a dense swarm of particles. Flow Turbul. Combust. 1 (1949) 27. | DOI | Zbl

[11] E. Burman, A penalty-free nonsymmetric Nitsche-type method for the weak imposition of boundary conditions. SIAM J. Numer. Anal. 50 (2012) 1959–1981. | DOI | MR | Zbl

[12] E. Burman and P. Hansbo, A unified stabilized method for Stokes’ and Darcy’s equations. J. Comput. Appl. Math. 198 (2007) 35–51. | DOI | MR | Zbl

[13] P.G. Ciarlet, The finite element method for elliptic problems. Reprint of the 1978 original [North-Holland, Amsterdam]. In Vol. 40 of Classics in Applied Mathematics. Society for Industrial and Applied Mathematics (SIAM), Philadelphia, PA (2002). | MR | Zbl

[14] P. Clément, Approximation by finite element functions using local regularization. Rev. Française Automat. Informat. Recherche Opérationnelle Sér. 9 (1975) 77–84. | Numdam | MR | Zbl

[15] C. D’Angelo and P. Zunino, Robust numerical approximation of coupled Stokes’ and Darcy’s flows applied to vascular hemodynamics and biochemical transport. ESAIM: M2AN 45 (2011) 447–476. | DOI | Numdam | MR | Zbl

[16] M. Discacciati, Domain Decomposition Methods for the Coupling of Surface and Groundwater Flows, Ph.D. thesis, EPFL, Lausanne)2004).

[17] J. Douglas Jr and J.P. Wang, An absolutely stabilized finite element method for the Stokes problem. Math. Comp. 52 (1989) 495–508. | DOI | MR | Zbl

[18] A. Ern and J.-L. Guermond, Theory and practice of finite elements. In Vol. 159 of Applied Mathematical Sciences. Springer-Verlag, New York, NY (2004). | DOI | MR | Zbl

[19] L.P. Franca and T.J.R. Hughes, Two classes of mixed finite element methods. Comput. Methods Appl. Mech. Eng. 69 (1988) 89–129. | DOI | MR | Zbl

[20] Freund J. and Stenberg R., On weakly imposed boundary conditions for second order problems. In: Proceedings of the international Conference on Finite Elements in Fluids – New trends and Applications, Venezia (1995).

[21] V. Girault and P.-A. Raviart, Finite element methods for Navier–Stokes equations, Theory and Algorithms. In Vol. 5 of Springer Series in Computational Mathematics. Springer-Verlag, Berlin (1986). | DOI | MR | Zbl

[22] A. Hannukainen, M. Juntunen, R. Stenberg, Computations with finite element methods for the Brinkman problem. Comput. Geosci. 15 (2011) 155–166. | DOI | Zbl

[23] P. Hansbo and M. Juntunen, Weakly imposed Dirichlet boundary conditions for the Brinkman model of porous media flow. Appl. Numer. Math. 59 (2009) 1274–1289. | DOI | MR | Zbl

[24] Q. Hu, F. Chouly, P. Hu, G. Cheng, S.P.A. Bordas, Skew-symmetric Nitsches formulation in isogeometric analysis: Dirichlet and symmetry conditions, patch coupling and frictionless contact. Comput. Methods Appl. Mech. Eng. 341 (2018) 188–220. | DOI | MR | Zbl

[25] T.J.R. Hughes and L.P. Franca, A new finite sed bouelement formulation for computational fluid dynamics. VII. The Stokes problem with various well-posed boundary conditions: symmetric formulations that converge for all velocity/pressure spaces. Comput. Methods Appl. Mech. Eng. 65 (1987) 85–96. | DOI | MR | Zbl

[26] E.W. Jenkins, V. John, A. Linke, L.G. Rebholz, On the parameter choice in grad-div stabilization for the Stokes equations. Adv. Comput. Math. 40 (2014) 491–516. | DOI | MR | Zbl

[27] M. Juntunen and R. Stenberg, Nitsche’s method for general boundary conditions. Math. Comput. 78 (2009) 1353–1374. | DOI | MR | Zbl

[28] M. Juntunen and R. Stenberg, Analysis of finite element methods for the Brinkman problem. Calcolo 47 (2010) 129–147. | DOI | MR | Zbl

[29] K.A. Mardal, X.-C. Tai, R. Winther, A robust finite element method for Darcy-Stokes ow. SIAM J. Numer. Anal. 40 (2002) 1605–1631. | DOI | MR | Zbl

[30] A. Masud and T.J.R. Hughes, A stabilized mixed finite element method for Darcy flow. Comput. Methods Appl. Mech. Eng. 191 (2002) 4341–4370. | DOI | MR | Zbl

[31] J. Nitsche, über ein Variationsprinzip zur Lösung von Dirichlet-Problemen bei Verwendung von Teilräumen, die keinen Randbedingungen unterworfen sind. Collection of articles dedicated to Lothar Collatz on his sixtieth birthday. Abh. Math. Sem. Univ. Hamburg 36 (1971) 9–15. | DOI | MR | Zbl

[32] L. Ridgway Scott and S. Zhang, Finite element interpolation of nonsmooth functions satisfying boundary conditions. Math. Comput. 54 (1990) 483–493. | DOI | MR | Zbl

[33] R. Stenberg, On some techniques for approximating boundary conditions in the finite element method. International Symposium on Mathematical Modelling and Computational Methods Modelling 94 (Prague, 1994). J. Comput. Appl. Math. 63 (1995) 139–148. | DOI | MR | Zbl

[34] V. Thomée, Galerkin finite element methods for parabolic problems, 2nd edition. In Vol. 25 of Springer Series in Computational Mathematics. Springer-Verlag, Berlin (2006). | MR | Zbl

[35] U. Wilbrandt, C. Bartsch, N. Ahmed, N. Alia, F. Anker, L. Blank, A. Caiazzo, S. Ganesan, S. Giere, G. Matthies, R. Meesala, A. Shamim, J. Venkatesan, V. John, ParMooN – A modernized program package based on mapped finite elements. Comput. Math. Appl. 74 (2017) 74–88. | DOI | MR | Zbl