Penalty method with Crouzeix–Raviart approximation for the Stokes equations under slip boundary condition

• Published in: ESAIM Mathematical Modelling and Numerical Analysis Vol. 53; no. 3; pp. 869 - 891
• Main Authors: Kashiwabara, Takahito, Oikawa, Issei, Zhou, Guanyu
• Format: Journal Article
• Language: English
• Published: Les Ulis EDP Sciences 01.05.2019
• Subjects: 35Q30; 65N30; Approximation; discrete H1/2-norm; domain perturbation; Navier-Stokes equations; Nonconforming FEM; slip boundary condition; Stokes equations
• ISSN: 0764-583X; 1290-3841; 1290-3841
• DOI: 10.1051/m2an/2019008

The Stokes equations subject to non-homogeneous slip boundary conditions are considered in a smooth domain Ω ⊂ ℝN (N = 2,3). We propose a finite element scheme based on the nonconforming P1/P0 approximation (Crouzeix–Raviart approximation) combined with a penalty formulation and with reduced-order numerical integration in order to address the essential boundary condition u · n∂Ω = g on ∂Ω. Because the original domain Ω must be approximated by a polygonal (or polyhedral) domain Ωh before applying the finite element method, we need to take into account the errors owing to the discrepancy Ω ≠ Ωh, that is, the issues of domain perturbation. In particular, the approximation of n∂Ω by n∂Ωh makes it non-trivial whether we have a discrete counterpart of a lifting theorem, i.e., continuous right inverse of the normal trace operator H1 (Ω)N → H1/2(∂Ω); u ↦ u⋅n∂Ω. In this paper we indeed prove such a discrete lifting theorem, taking advantage of the nonconforming approximation, and consequently we establish the error estimates O(hα + ε) and O(h2α + ε) for the velocity in the H1- and L2-norms respectively, where α = 1 if N = 2 and α = 1/2 if N = 3. This improves the previous result [T. Kashiwabara et al., Numer. Math. 134 (2016) 705–740] obtained for the conforming approximation in the sense that there appears no reciprocal of the penalty parameter ϵ in the estimates.