Conditional quasi-optimal error estimate for a finite element discretization of the $p$-Navier–Stokes equations: The case $p>2$

Abstract

Report issue for preceding element

In this paper, we derive quasi-optimal a priori error estimates for the kinematic pressure for a Finite Element (FE) approximation of steady systems of $p$-Navier–Stokes type in the case of shear-thickening, i.e., in the case $p>2$, imposing a new mild Muckenhoupt regularity condition.

Report issue for preceding element

1 Introduction

Report issue for preceding element

In the present paper, we examine a Finite Element (FE) approximation of steady systems of $p$-Navier–Stokes type, i.e.,

Report issue for preceding element

$$-\mathrm{div}\,\mathbf{S}({\bf D}{\bf v})+[\nabla{\bf v}]{\bf v}+\nabla q ={\bf f}\qquad \text{in }\Omega\,,$$ (1.1)

$$\mathrm{div}\,{\bf v} =0\qquad \text{in }\Omega\,,$$

$${\bf v} =\mathbf{0} \text{on }\partial\Omega\,{\color[rgb]{0,0,0}\definecolor[named]{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@gray@stroke{0}\pgfsys@color@gray@fill{0},}$$

for quasi-optimal a priori error estimates for the kinematic pressure in the case of shear-thickening fluids, i.e., $p>2$. The system (1.1) describes the steady motion of a homogeneous, incompressible fluid with shear-dependent viscosity. More precisely, for a given vector field ${\bf f}\colon\Omega\to\mathbb{R}^{d}$ describing external forces, an incompressibility constraint (1.1)₂, and a non-slip boundary condition (1.1)₃, the system (1.1) seeks for a velocity vector field ${\bf v}=(v_{1},\ldots,v_{d})^{\top}\colon\Omega\to\mathbb{R}^{d}$ and a kinematic pressure $q\colon\Omega\to\mathbb{R}$ solving (1.1). Here, $\Omega\subseteq\mathbb{R}^{d}$, $d\in{\{{2,3}\}}$, is a bounded, polygonal (if $d=2$) or polyhedral (if ${d=3}$) Lipschitz domain. The extra stress tensor $\mathbf{S}({\bf D}{\bf v})\colon\Omega\to\smash{\mathbb{R}^{d\times d}_{\textup{sym}}}$ depends on the strain rate tensor ${\bf D}{\bf v}\coloneqq\frac{1}{2}(\nabla{\bf v}+\nabla{\bf v}^{\top})\colon\Omega\to\smash{\mathbb{R}^{d\times d}_{\textup{sym}}}$, i.e., the symmetric part of the velocity gradient $\nabla{\bf v}\coloneqq(\partial_{j}v_{i})_{i,j=1,\ldots,d}\colon\Omega\to\mathbb{R}^{d\times d}$. The convective term $\smash{[\nabla{\bf v}]{\bf v}\colon\Omega\to\mathbb{R}^{d}}$ is defined via $\smash{([\nabla{\bf v}]{\bf v})_{i}\coloneqq\sum_{j=1}^{d}{v_{j}\partial_{j}v_{i}}}$ for all $i=1,\ldots,d$.

Report issue for preceding element

Throughout the paper, we assume that the extra stress tensor $\mathbf{S}$ has $(p,\delta)$-structure (cf. Assumption 2.4). The relevant example falling into this class is $\mathbf{S}\colon\mathbb{R}^{d\times d}\to\mathbb{R}^{d\times d}_{\textup{sym}}$, for every ${\bf A}\in\mathbb{R}^{d\times d}$ defined via

Report issue for preceding element

$$\mathbf{S}({\bf A})\coloneqq\mu_{0}\,(\delta+|{\bf A}|)^{p-2}{\bf A}\,,$$ (1.2)

where $p\in(1,\infty)$, $\delta\geq 0$, and $\mu_{0}>0$.

Report issue for preceding element

The a priori error analysis of the steady $p$-Navier–Stokes problem (1.1) using FE approximations is by now well-understood: recently, in [20], a priori error estimates in the case of shear-thickening, i.e., in the case $p>2$, were derived, which are optimal for the velocity vector field, but sub-optimal for the kinematic pressure. More precisely, this lacuna is mainly due to the following technical hurdle, here, exemplified by the error analysis of the FE approximations of the steady $p$-Navier–Stokes problem (1.1) used in [20]: in it, it turns out that the error of the kinematic pressure measured in the squared $L^{p^{\prime}}(\Omega)$-norm is bounded by the squared $L^{(\varphi_{\smash{|{\bf D}{\bf v}|}})^{*}}(\Omega)$-norm of the extra-stress errors, i.e., it holds that

Report issue for preceding element

$$\begin{aligned} \|q_{h}-q\|_{p^{\prime},\Omega}^{2}&\leq\|{\bf S}({\bf D}{\bf v}_{h})-{\bf S}({\bf D}{\bf v})\|_{(\varphi_{\smash{|{\bf D}{\bf v}|}})^{*},\Omega}^{2}+\textup{(h.o.t.)}\,.\end{aligned}$$ (1.3)

The relation (1.3) is a consequence of the squared discrete inf-sup stability result (cf. [9, Lem. 6.10])

Report issue for preceding element

$$c\,\|z_{h}\|_{p^{\prime},\Omega}^{2}\leq\sup_{{\bf z}_{h}\in{\mathaccent 23{V}}_{h}\;:\;\|\nabla{\bf z}_{h}\|_{p,\Omega}\leq 1}{\big{(}(z_{h},\mathrm{div}\,{\bf z}_{h})_{\Omega}\big{)}^{2}}\,,$$ (1.4)

valid for each $z_{h}\in{\mathaccent 23{Q}}_{h}$, where ${\mathaccent 23{V}}_{h}\subseteq(W^{1,p}_{0}(\Omega))^{d}$ denotes a discrete velocity vector space and and ${\mathaccent 23{Q}}_{h}\subseteq L^{p^{\prime}}_{0}(\Omega)$ a discrete pressure space jointly forming a discretely inf-sup-stable finite element couple. Then, using the estimates (cf. [22, Lem. 4.4(ii)])

Report issue for preceding element

$$\|{\bf S}({\bf D}{\bf v}_{h})-{\bf S}({\bf D}{\bf v})\|_{(\varphi_{\smash{|{\bf D}{\bf v}|}})^{*},\Omega}^{2}\leq c\,\|{\bf F}({\bf D}{\bf v}_{h})-{\bf F}({\bf D}{\bf v})\|_{2,\Omega}^{2}\,,$$

from the relation (1.3), it follows that

Report issue for preceding element

$$\|q_{h}-q\|_{p^{\prime},\Omega}^{2}\leq c\,\|{\bf F}({\bf D}{\bf v}_{h})-{\bf F}({\bf D}{\bf v})\|_{2,\Omega}^{2}+\textup{(h.o.t.)}\,,$$ (1.5)

i.e., the kinematic pressure error measured in the squared $L^{p^{\prime}}(\Omega)$-norm is bounded by the squared velocity vector field error. However, numerical experiments (cf. [22]) indicate that the relation (1.5) is potentially sub-optimal and suggest instead the relation

Report issue for preceding element

$$\rho_{(\varphi_{|{\bf D}{\bf v}|})^{*},\Omega}(q_{h}-q)\leq c\,\|{\bf F}({\bf D}{\bf v}_{h})-{\bf F}({\bf D}{\bf v})\|_{2,\Omega}^{2}+\textup{(h.o.t.)}\,.$$ (1.6)

In [23], in the case of a Local Discontinuous Galerkin (LDG) approximation, such a discrete analogue has been established under the additional assumption that the viscosity of the fluid is a Muckenhoupt weight of class $2$, i.e., if we have that

Report issue for preceding element

$$\mu_{{\bf D}{\bf v}}\coloneqq(\delta+|\overline{{\bf D}{\bf v}}|)^{p-2}\in A_{2}(\mathbb{R}^{d})\,,$$ (1.7)

where $\overline{{\bf D}{\bf v}}\coloneqq{\bf D}{\bf v}$ a.e. in $\Omega$ and $\overline{{\bf D}{\bf v}}\coloneqq{\mathbf{0}}$ a.e.  in $\mathbb{R}^{d}\setminus\Omega$. In [27, 28], it turned out that the Muckenhoupt regularity assumption (1.7) can not be expected, in general, in the three-dimensional case. However, in the two-dimensional case, regularity results (cf. [24]) suggest that (1.7) is satisfied under mild assumptions. Finally, if the Muckenhoupt regularity condition (1.7) is satisfied, from the relation (1.6), we are able to derive an alternative a priori error estimate for the pressure which turns out to be quasi-optimal and, thus, is the first quasi-optimal a priori error estimate for the pressure in the case of shear-thickening, i.e., $p>2$.

Report issue for preceding element

This paper is organized as follows: In Section 2, we introduce the employed notation, define relevant function spaces, basic assumptions on the extra stress tensor $\mathbf{S}$ and its consequences, the weak formulations Problem (Q) and Problem (P) of the system (1.1), and the discrete operators. In Section 3, we introduce the discrete weak formulations Problem (Q_h) and Problem (P_h), recall known a priori error estimates and prove the main result of the paper, i.e., a quasi-optimal (with respect to the Muckenhoupt regularity condition (1.7)) a priori error estimate for the kinematic pressure in the case $p>2$ (cf. Theorem 3.15, Corollary 3.16). In Section 4, we review the theoretical findings via numerical experiments.

Report issue for preceding element

2 Preliminaries

Report issue for preceding element

2.1 Basic Notation

Report issue for preceding element

We employ $c,C>0$ to denote generic constants, that may change from line to line, but do not depend on the crucial quantities. Moreover, we write $u\sim v$ if and only if there exist constants $c,C>0$ such that $c\,u\leq v\leq C\,u$.

Report issue for preceding element

For $k\in\mathbb{N}$ and $p\in[1,\infty]$, we employ the customary Lebesgue spaces $(L^{p}(\Omega),\|\cdot\|_{p,\Omega})$ and Sobolev spaces $(W^{k,p}(\Omega),\|\cdot\|_{k,p,\Omega})$, where $\Omega\subseteq\mathbb{R}^{d}$, $d\in\{2,3\}$, is a bounded, polygonal (if $d=2$) or polyhedral (if $d=3$) Lipschitz domain. The space $\smash{W^{1,p}_{0}(\Omega)}$ is defined as those functions from $W^{1,p}(\Omega)$ whose traces vanish on $\partial\Omega$. The Hölder dual exponent is defined via $p^{\prime}=\tfrac{p}{p-1}\in[1,\infty]$. We denote vector-valued functions by boldface letters and tensor-valued functions by capital boldface letters. The Euclidean inner product between two vectors ${\bf a}=(a_{1},\ldots,a_{d})^{\top},{\bf b}=(b_{1},\ldots,b_{d})^{\top}\in\mathbb{R}^{d}$ is denoted by ${\bf a}\cdot{\bf b}\coloneqq\sum_{i=1}^{d}{a_{i}b_{i}}$, while the Frobenius inner product between two tensors ${\bf A}=(A_{ij})_{i,j\in\{1,\ldots,d\}},$ ${\bf B}=(B_{ij})_{i,j\in\{1,\ldots,d\}}\in\mathbb{R}^{d\times d}$ is denoted by ${\bf A}:{\bf B}\coloneqq\sum_{i,j=1}^{d}{A_{ij}B_{ij}}$. Moreover, for a (Lebesgue) measurable set $M\subseteq\mathbb{R}^{n}$, $n\in\mathbb{N}$, and (Lebesgue) measurable functions, vector or tensor fields $\mathbf{u},\mathbf{w}\in(L^{0}(M))^{\ell}$, $\ell\in\mathbb{N}$, we write

Report issue for preceding element

$$(\mathbf{u},\mathbf{w})_{M}\coloneqq\int_{M}\mathbf{u}\odot\mathbf{w}\,\mathrm{d}x\,,$$

whenever the right-hand side is well-defined, where $\odot\colon\mathbb{R}^{\ell}\times\mathbb{R}^{\ell}\to\mathbb{R}$ either denotes scalar multiplication, the Euclidean inner product or the Frobenius inner product. The integral mean of an integrable function, vector or tensor field $\mathbf{u}\in(L^{0}(M))^{\ell}$, $\ell\in\mathbb{N}$, over a (Lebesgue) measurable set $M\subseteq\mathbb{R}^{n}$, $n\in\mathbb{N}$, with $|M|>0$ is denoted by $\langle{\mathbf{u}}\rangle_{M}\coloneqq\smash{\frac{1}{|M|}\int_{M}\mathbf{u}\,\mathrm{d}x}$.

Report issue for preceding element

2.2 N-functions and Orlicz spaces

Report issue for preceding element

A convex function $\psi\colon\mathbb{R}^{\geq 0}\to\mathbb{R}^{\geq 0}$ is called N-function if it holds that ${\psi(0)=0}$, ${\psi(t)>0}$ for all ${t>0}$, $\lim_{t\rightarrow 0}\psi(t)/t=0$, and $\lim_{t\rightarrow\infty}\psi(t)/t=\infty$. A Carathéodory function $\psi\colon M\times\mathbb{R}^{\geq 0}\to\mathbb{R}^{\geq 0}$, where $M\subseteq\mathbb{R}^{n}$, $n\in\mathbb{N}$, is a (Lebesgue) measurable set, such that $\psi(x,\cdot)$ is an N-function for a.e. $x\in M$, is called generalized N-function. The modular is defined via $\rho_{\psi}(f)\coloneqq\rho_{\psi,\Omega}(f)\coloneqq\int_{\Omega}\psi({\lvert{f}\rvert})\,\textup{d}x$ if $\psi$ is an N-function, and via $\rho_{\psi}(f)\coloneqq\rho_{\psi,\Omega}(f)\coloneqq\int_{\Omega}\psi(x,{\lvert{f(x)}\rvert})\,\textup{d}x$, if $\psi$ is a generalized N-function. We define the (convex) conjugate N-function We define the (convex) conjugate N-function $\psi^{*}\colon M\times\mathbb{R}^{\geq 0}\to\mathbb{R}^{\geq 0}$ via ${\psi^{*}(t,x)\coloneqq\sup_{s\geq 0}(st-\psi(s,x))}$ for all $t\geq 0$ and a.e. $x\in M$, which satisfies $(\partial_{t}(\psi^{*}))(x,t)=(\partial_{t}\psi)^{-1}(x,t)$ for all $t\geq 0$ and a.e. $x\in M$. A (generalized) N-function $\psi$ satisfies the $\Delta_{2}$-condition (in short, $\psi\in\Delta_{2}$), if there exists $K>2$ such that ${\psi(2\,t,x)\leq K\,\psi(t,x)}$ for all $t\geq 0$ and a.e. $x\in M$. The smallest such constant is denoted by $\Delta_{2}(\psi)>0$. We need the following version of the $\varepsilon$-Young inequality: for every ${\varepsilon>0}$, there exists a constant $c_{\varepsilon}>0$, depending only on $\Delta_{2}(\psi),\Delta_{2}(\psi^{*})<\infty$, such that for every $s,t\geq 0$ and a.e. $x\in M$, it holds that

Report issue for preceding element

$$s\,t \leq c_{\varepsilon}\,\psi^{*}(s,x)+\varepsilon\,\psi(t,x)\,.$$ (2.1)

2.3 Basic properties of the extra stress tensor

Report issue for preceding element

Throughout the entire paper, we will always assume that the extra stress tensor $\mathbf{S}$ has $(p,\delta)$-structure. A detailed discussion and full proofs can be found, e.g., in [10, 31]. For a given tensor ${\bf A}\in\mathbb{R}^{d\times d}$, we denote its symmetric part by ${{\bf A}^{\textup{sym}}\coloneqq\frac{1}{2}({\bf A}+{\bf A}^{\top})\in\mathbb{R}^{d\times d}_{\textup{sym}}\coloneqq\{{\bf A}\in\mathbb{R}^{d\times d}\mid{\bf A}={\bf A}^{\top}\}}$.

Report issue for preceding element

For $p\in(1,\infty)$ and $\delta\geq 0$, we define the special N-function $\varphi=\varphi_{p,\delta}\colon\mathbb{R}^{\geq 0}\to\mathbb{R}^{\geq 0}$ via

Report issue for preceding element

$$\varphi(t)\coloneqq\int_{0}^{t}\varphi^{\prime}(s)\,\mathrm{d}s,\quad\text{where}\quad\varphi^{\prime}(t)\coloneqq(\delta+t)^{p-2}t\,,\quad\textup{ for all }t\geq 0\,.$$ (2.2)

An important tool in our analysis plays shifted N-functions $\{\psi_{a}\}_{\smash{a\geq 0}}$ (cf. [11, 31]). For a given N-function $\psi\colon\mathbb{R}^{\geq 0}\to\mathbb{R}^{\geq 0}$, we define the family of shifted N-functions ${\psi_{a}\colon\mathbb{R}^{\geq 0}\to\mathbb{R}^{\geq 0}}$, ${a\geq 0}$, via

Report issue for preceding element

$$\psi_{a}(t)\coloneqq\int_{0}^{t}\psi_{a}^{\prime}(s)\,\mathrm{d}s\,,\quad\text{where }\quad\psi^{\prime}_{a}(t)\coloneqq\psi^{\prime}(a+t)\frac{t}{a+t}\,,\quad\textup{ for all }t\geq 0\,.$$ (2.3)

Assumption 2.4 (Extra stress tensor).

Report issue for preceding element

We assume that the extra stress tensor $\mathbf{S}\colon\mathbb{R}^{d\times d}\to\mathbb{R}^{d\times d}_{\textup{sym}}$ belongs to $C^{0}(\mathbb{R}^{d\times d};\mathbb{R}^{d\times d}_{\textup{sym}})\cap C^{1}(\mathbb{R}^{d\times d}\setminus\{\mathbf{0}\};\mathbb{R}^{d\times d}_{\textup{sym}})$ and satisfies $\mathbf{S}({\bf A})=\mathbf{S}({\bf A}^{\textup{sym}})$ for all ${\bf A}\in\mathbb{R}^{d\times d}$ and $\mathbf{S}(\mathbf{0})=\mathbf{0}$. Moreover, we assume that the tensor $\mathbf{S}=(S_{ij})_{i,j=1,\ldots,d}$ has $(p,\delta)$-structure, i.e., for some $p\in(1,\infty)$, $\delta\in[0,\infty)$, and the N-function $\varphi=\varphi_{p,\delta}$ (cf. (2.2)), there exist constants $C_{0},C_{1}>0$ such that

Report issue for preceding element

$$\sum\limits_{i,j,k,l=1}^{d}\partial_{kl}S_{ij}({\bf A})B_{ij}B_{kl} \geq C_{0}\,\frac{\varphi^{\prime}(|{\bf A}^{\textup{sym}}|)}{|{\bf A}^{\textup{sym}}|}\,|{\bf B}^{\textup{sym}}|^{2}\,,$$ (2.5)

$$\big{|}\partial_{kl}S_{ij}({{\bf A}})\big{|} \leq C_{1}\,\frac{\varphi^{\prime}(|{\bf A}^{\textup{sym}}|)}{|{\bf A}^{\textup{sym}}|}$$ (2.6)

are satisfied for all ${\bf A},{\bf B}\in\mathbb{R}^{d\times d}$ with ${\bf A}^{\textup{sym}}\neq\mathbf{0}$ and all $i,j,k,l=1,\ldots,d$. The constants $C_{0},C_{1}>0$ and $p\in(1,\infty)$ are called the characteristics of $\mathbf{S}$.

Report issue for preceding element

Remark 2.7.

Report issue for preceding element

• [(ii)]

  Report issue for preceding element

• Assume that $\mathbf{S}$ satisfies Assumption 2.4 for some $\delta\in[0,\delta_{0}]$. Then, if not otherwise stated, the constants in the estimates depend only on the characteristics of $\mathbf{S}$ and $\delta_{0}\geq 0$, but are independent of $\delta\geq 0$.

  Report issue for preceding element

• Let $\varphi$ be defined in (2.2) and $\{\varphi_{a}\}_{a\geq 0}$ be the corresponding family of shifted N-functions. Then, the operators $\mathbf{S}_{a}\colon\mathbb{R}^{d\times d}\to\smash{\mathbb{R}_{\textup{sym}}^{d\times d}}$, $a\geq 0$, for every $a\geq 0$ and ${\bf A}\in\mathbb{R}^{d\times d}$ defined via

  Report issue for preceding element

  $$\mathbf{S}_{a}({\bf A})\coloneqq\frac{\varphi_{a}^{\prime}({\lvert{{\bf A}^{\textup{sym}}}\rvert})}{{\lvert{{\bf A}^{\textup{sym}}}\rvert}}\,{\bf A}^{\textup{sym}}\,,$$ (2.8)

  have $(p,\delta+a)$-structure. In this case, the characteristics of $\mathbf{S}_{a}$ depend only on ${p\in(1,\infty)}$ and are independent of $\delta\geq 0$ and $a\geq 0$.

  Report issue for preceding element

Closely related to the extra stress tensor $\mathbf{S}\colon\mathbb{R}^{d\times d}\to\mathbb{R}^{d\times d}_{\textup{sym}}$ with $(p,\delta)$-structure is the non-linear mapping ${\bf F}\colon\mathbb{R}^{d\times d}\to\mathbb{R}^{d\times d}_{\textup{sym}}$, for every ${\bf A}\in\mathbb{R}^{d\times d}$ defined via

Report issue for preceding element

$$\begin{aligned} {\bf F}({\bf A})&\coloneqq(\delta+|{\bf A}^{\textup{sym}}|)^{\frac{p-2}{2}}{\bf A}^{\textup{sym}}\,.\end{aligned}$$ (2.9)

The connections between $\mathbf{S},{\bf F}\colon\mathbb{R}^{d\times d}\to\mathbb{R}^{d\times d}_{\textup{sym}}$ and $\varphi_{a},(\varphi_{a})^{*}\colon\mathbb{R}^{\geq 0}\to\mathbb{R}^{\geq 0}$, ${a\geq 0}$, are best explained by the following proposition (cf. [10, 31, 14]).

Report issue for preceding element

Proposition 2.10.

Report issue for preceding element

Let $\mathbf{S}$ satisfy Assumption 2.4, let $\varphi$ be defined in (2.2), and let ${\bf F}$ be defined in (2.9). Then, uniformly with respect to ${\bf A},{\bf B}\in\mathbb{R}^{d\times d}$, we have that

Report issue for preceding element

$$\big{(}\mathbf{S}({\bf A})-\mathbf{S}({\bf B})\big{)}:({\bf A}-{\bf B}) \sim|{\bf F}({\bf A})-{\bf F}({\bf B})|^{2}$$

$$\sim\varphi_{|{\bf A}^{\textup{sym}}|}(|{\bf A}^{\textup{sym}}-{\bf B}^{\textup{sym}}|)$$ (2.11)

$$\sim(\varphi_{|{\bf A}^{\textup{sym}}|})^{*}({\lvert{\mathbf{S}({\bf A})-\mathbf{S}({\bf B})}\rvert})\,,$$

$${\lvert{\mathbf{S}({\bf A})-\mathbf{S}({\bf B})}\rvert} \sim\smash{\varphi^{\prime}_{{\lvert{{\bf A}}\rvert}}({\lvert{{\bf A}-{\bf B}}\rvert})}\,.$$ (2.12)

The constants depend only on the characteristics of ${\mathbf{S}}$.

Report issue for preceding element

Remark 2.13.

Report issue for preceding element

For the operators $\mathbf{S}_{a}\colon\mathbb{R}^{d\times d}\to\smash{\mathbb{R}_{\textup{sym}}^{d\times d}}$, $a\geq 0$, defined in (2.8), the assertions of Proposition 2.10 hold with $\varphi\colon\mathbb{R}^{\geq 0}\to\mathbb{R}^{\geq 0}$ replaced by $\varphi_{a}\colon\mathbb{R}^{\geq 0}\to\mathbb{R}^{\geq 0}$, $a\geq 0$.

Report issue for preceding element

The following results can be found in [11, 31].

Report issue for preceding element

Lemma 2.14 (Change of Shift).

Report issue for preceding element

Let $\varphi$ be defined in (2.2) and let ${\bf F}$ be defined in (2.9). Then, for each $\varepsilon>0$, there exists $c_{\varepsilon}\geq 1$ (depending only on $\varepsilon>0$ and the characteristics of $\varphi$) such that for every ${\bf A},{\bf B}\in\smash{\mathbb{R}^{d\times d}_{\textup{sym}}}$ and $t\geq 0$, it holds that

Report issue for preceding element

$$\smash{\varphi_{|{\bf B}|}(t)} \leq\smash{c_{\varepsilon}\,\varphi_{|{\bf A}|}(t)+\varepsilon\,|{\bf F}({\bf B})-{\bf F}({\bf A})|^{2}\,,}$$ (2.15)

$$\smash{\varphi_{|{\bf B}|}(t)} \leq\smash{c_{\varepsilon}\,\varphi_{|{\bf A}|}(t)+\varepsilon\,\varphi_{|{\bf A}|}(||{\bf B}|-|{\bf A}||)\,,}$$ (2.16)

$$\smash{(\varphi_{|{\bf B}|})^{*}(t)} \leq\smash{c_{\varepsilon}\,(\varphi_{|{\bf A}|})^{*}(t)+\varepsilon\,|{\bf F}({\bf B})-{\bf F}({\bf A})|^{2}}\,,$$ (2.17)

$$\smash{(\varphi_{|{\bf B}|})^{*}(t)} \leq\smash{c_{\varepsilon}\,(\varphi_{|{\bf A}|})^{*}(t)+\varepsilon\,\varphi_{|{\bf A}|}(||{\bf B}|-|{\bf A}||)}\,.$$ (2.18)

2.4 The $p$-Navier–Stokes system

Report issue for preceding element

Let us briefly recall some well-known facts about the $p$-Navier–Stokes equations (1.1). For $p\in(1,\infty)$, we define the function spaces

Report issue for preceding element

$${\mathaccent 23{V}}\coloneqq(W^{1,p}_{0}(\Omega))^{d}\,,\qquad{\mathaccent 23{Q}}\coloneqq L_{0}^{p^{\prime}}(\Omega)\coloneqq\big{\{}z\in L^{p^{\prime}}(\Omega)\;|\;\langle{z}\rangle_{\Omega}=0\big{\}}\,.$$

With this notation, assuming that $p\geq\frac{3d}{d+2}$, the weak formulation of the $p$-Navier–Stokes equations (1.1) as a non-linear saddle point problem is the following:

Report issue for preceding element

Problem (Q). For given ${\bf f}\in(L^{p^{\prime}}(\Omega))^{d}$, find $({\bf v},q)\in{\mathaccent 23{V}}\times{\mathaccent 23{Q}}$ such that for all $({\bf z},z)^{\top}\in{\mathaccent 23{V}}\times Q$, it holds     that

Report issue for preceding element

$$(\mathbf{S}({\bf D}{\bf v}),{\bf D}{\bf z})_{\Omega}+([\nabla{\bf v}]{\bf v},{\bf z})_{\Omega}-(q,\mathrm{div}\,{\bf z})_{\Omega} =({\bf f},{\bf z})_{\Omega}\,,$$ (2.19)

$$(\mathrm{div}\,{\bf v},z)_{\Omega} =0\,.$$ (2.20)

Alternatively, we can reformulate Problem (Q) “hiding” the kinematic pressure.

Report issue for preceding element

Problem (P). For given ${\bf f}\in(L^{p^{\prime}}(\Omega))^{d}$, find ${\bf v}\in{\mathaccent 23{V}}(0)$ such that for all ${\bf z}\in{\mathaccent 23{V}}(0)$, it holds that

Report issue for preceding element

$$(\mathbf{S}({\bf D}{\bf v}),{\bf D}{\bf z})_{\Omega}+([\nabla{\bf v}]{\bf v},{\bf z})_{\Omega} =({\bf f},{\bf z})_{\Omega}\,,$$ (2.21)

where ${\mathaccent 23{V}}(0)\coloneqq\{{\bf z}\in{\mathaccent 23{V}}\mid\mathrm{div}\,{\bf z}=0\}$.
The names Problem (Q) and Problem (P) are traditional in the literature (cf.  [7, 4]). The well-posedness of Problem (Q) and Problem (P) is usually established in two steps: first, using pseudo-monotone operator theory (cf. [29]), the well-posedness of Problem (P) is shown; then, given the well-posedness of Problem (P), the well-posedness of Problem (Q) follows using DeRham’s lemma. There holds the following regularity property of the pressure if the velocity satisfies a natural regularity assumption.

Report issue for preceding element

Lemma 2.22.

Report issue for preceding element

Let $\mathbf{S}$ satisfy Assumption 2.4 with $p\in(2,\infty)$ and $\delta>0$, and let $({\bf v},q)^{\top}\in{\mathaccent 23{V}}(0)\times{\mathaccent 23{Q}}$ be a weak solution of Problem (Q) with ${\bf F}({\bf D}{\bf v})\in(W^{1,2}(\Omega))^{d\times d}$. Then, the following statements apply:

Report issue for preceding element

• [(ii)]

  Report issue for preceding element

• If ${\bf f}\in(L^{p^{\prime}}(\Omega))^{d}$, then $\nabla q\in(L^{p^{\prime}}(\Omega))^{d}$.

  Report issue for preceding element

• If ${\bf f}\in(L^{2}(\Omega))^{d}$, then $(\delta+|{\bf D}{\bf v}|)^{2-p}|\nabla q|^{2}\in L^{1}(\Omega)$.

  Report issue for preceding element

Proof.

Report issue for preceding element

See [21, Lem. 2.6]. ∎

Report issue for preceding element

2.5 Discussion of Muckenhoupt regularity condition

Report issue for preceding element

In this subsection, we examine the Muckenhoupt regularity condition

Report issue for preceding element

$$\mu_{{\bf D}{\bf v}}\coloneqq(\delta+|\overline{{\bf D}{\bf v}}|)^{p-2}\in A_{2}(\mathbb{R}^{d})\,,$$ (2.23)

for a solution ${\bf v}\in(W^{1,p}_{0}(\Omega))^{d}$ of Problem (P) (or Problem (Q), respectively), where $\overline{{\bf D}{\bf v}}\in(L^{p}(\mathbb{R}^{d}))^{d\times d}$ is defined via

Report issue for preceding element

$$\overline{{\bf D}{\bf v}}\coloneqq\begin{cases}{\bf D}{\bf v}&\text{ a.e.\ in }\Omega\,,\\ {\mathbf{0}}&\text{ a.e.\ in }\mathbb{R}^{d}\setminus\Omega\,.\end{cases}$$

In this connection, recall that for given $p\in[1,\infty)$, a weight $\sigma\colon\mathbb{R}^{d}\to(0,+\infty)$, i.e., $\sigma\in L^{1}_{\textup{loc}}(\mathbb{R}^{d})$ and $0<\sigma(x)<+\infty$ for a.e. $x\in\mathbb{R}^{d}$, is said to satisfy the $A_{p}$-condition, if

Report issue for preceding element

$$[\sigma]_{A_{p}(\mathbb{R}^{d})}\coloneqq\sup_{B\subseteq\mathbb{R}^{d}\,:\,B\text{ is a ball}}{\langle\sigma\rangle_{B}(\langle\sigma^{1-p^{\prime}}\rangle_{B})^{p-1}}<\infty\,.$$

We denote by $A_{p}(\mathbb{R}^{d})$ the class of all weights satisfying the $A_{p}$-condition and use weighted Lebesgue spaces $L^{p}(\Omega;\sigma)$ equipped with the norm $\|\cdot\|_{p,\sigma,\Omega}\coloneqq(\int_{\Omega}|\cdot|^{p}\,\sigma\,\mathrm{d}x)^{1/p}$.

Report issue for preceding element

In two dimensions, the Muckenhoupt regularity condition (2.23) is satisfied under mild assumptions.

Report issue for preceding element

Theorem 2.24.

Report issue for preceding element

Let $\Omega\subseteq\mathbb{R}^{2}$ be a bounded domain with $C^{2}$-boundary, $p>2$, and ${\bf f}\in(L^{s}(\Omega))^{2}$, where ${s>2}$. Then, there exist $q>2$, $\alpha>0$, and a solution $({\bf v},q)^{\top}\in{\mathaccent 23{V}}(0)\times{\mathaccent 23{Q}}$ of Problem (Q) with the following properties:

Report issue for preceding element

• [(iii)]

  Report issue for preceding element

• $({\bf v},q)^{\top}\in(W^{2,q}(\Omega))^{2}\times W^{1,q}(\Omega)$;

  Report issue for preceding element

• $({\bf v},q)^{\top}\in(C^{1,\alpha}(\overline{\Omega}))^{2}\times C^{0,\alpha}(\overline{\Omega})$.

  Report issue for preceding element

Proof.

Report issue for preceding element

See [24, Thm. 6.1]. ∎

Report issue for preceding element

Remark 2.25.

Report issue for preceding element

• [(ii)]

  Report issue for preceding element

• For $|{\bf D}{\bf v}|\in L^{\infty}(\Omega)$, $\delta>0$, and $p>2$, it follows $\delta^{p-2}\leq\mu_{{\bf D}{\bf v}}\leq(\delta+\|{\bf D}{\bf v}\|_{\infty,\Omega})^{p-2}$ a.e. in $\Omega$, so that the Muckenhoupt regularity condition (2.23) is trivially satisfied. As a result, under the assumptions of Theorem 2.24, the regularity assumption (2.23) is satisfied.

  Report issue for preceding element

• We believe that it is possible to prove Theorem 2.24 without the $C^{2}$-boundary assumption for polygonal, convex domains.

  Report issue for preceding element

The following result implies that, in three dimensions, one cannot hope for the regularity assumption to be satisfied, in general.

Report issue for preceding element

Theorem 2.26.

Report issue for preceding element

Let $\Omega\subseteq\mathbb{R}^{d}$, $d\geq 2$, be a bounded domain with $C^{1,1}$-boundary. Then, there exists a vector field ${\bf v}\colon\Omega\to\mathbb{R}^{d}$ with the following properties:

Report issue for preceding element

• [(iii)]

  Report issue for preceding element

• ${\bf v}\in(W^{2,r}(\Omega))^{d}\cap(W^{1,s}_{0}(\Omega))^{d}$ for all $r\in(1,d)$ and $s\in(1,\infty)$;

  Report issue for preceding element

• $\textup{tr}_{\partial\Omega}(\nabla{\bf v})\not\equiv{\mathbf{0}}$;

  Report issue for preceding element

• $\mathrm{div}\,{\bf v}=0$ a.e. in $\Omega$;

  Report issue for preceding element

• $\mu_{{\bf D}{\bf v}}\coloneqq(\delta+|\overline{{\bf D}{\bf v}}|)^{p-2}\notin A_{2}(\mathbb{R}^{d})$ for all $p\in(1,\infty)\setminus\{2\}$ and $\delta\in[0,1]$.

  Report issue for preceding element

Proof.

Report issue for preceding element

See [27, Thm. 2.1] and [28, Thm. 2]. ∎

Report issue for preceding element

Remark 2.27.

Report issue for preceding element

Let ${\bf v}\colon\Omega\to\mathbb{R}^{3}$ be the vector field from Theorem 2.26 in the case $d=3$. Setting ${\bf f}\coloneqq\mathrm{div}\,({\bf v}\otimes{\bf v})-\mathrm{div}\,{\bf S}({\bf D}{\bf v})$, where ${\bf S}({\bf D}{\bf v})\coloneqq\mu_{{\bf D}{\bf v}}{\bf D}{\bf v}$ and $\mu_{{\bf D}{\bf v}}=(\delta+|{\bf D}{\bf v}|)^{p-2}\in A_{2}(\mathbb{R}^{3})$, we see that ${\bf v}$ and $q\equiv 0$ are a weak solution of the $p$-Navier–Stokes system (1.1) for any $p\in(1,\infty)$ and $\delta\in[0,1]$. Due to [27, Rem. 2.3] and [28], we have that ${\bf F}({\bf D}{\bf v})\in(W^{1,2}(\Omega))^{3\times 3}$ and ${\bf f}\in(L^{2}(\Omega))^{3}$ for all $p\in(1,\infty)$.

Report issue for preceding element

3 Finite Element (FE) approximation

Report issue for preceding element

3.1 Triangulations

Report issue for preceding element

Throughout the entire paper, we denote by $\{\mathscr{T}_{h}\}_{h>0}$ a family of triangulations of $\Omega\subseteq\mathbb{R}^{d}$, $d\in\{2,3\}$, consisting of $d$-dimensional simplices (cf. [15]). Here, $h>0$ refers to the maximal mesh-size, i.e., if we set $h_{K}\coloneqq\textup{diam}(K)$ for all $K\in\mathscr{T}_{h}$, then $h\coloneqq\max_{K\in\mathscr{T}_{h}}{h_{K}}$. For every $K\in\mathscr{T}_{h}$, we denote by $\rho_{K}>0$, the supremum of diameters of inscribed balls contained in $K$. We assume that there is a constant $\kappa_{0}>0$, independent of $h>0$, such that $\max_{K\in\mathscr{T}_{h}}{h_{K}}{\rho_{K}^{-1}}\leq\kappa_{0}$. The smallest such constant is called the chunkiness of $\{\mathscr{T}_{h}\}_{h>0}$. For every $K\in\mathscr{T}_{h}$, the element patch is defined via $\omega_{K}\coloneqq\bigcup\{K^{\prime}\in\mathscr{T}_{h}\mid K^{\prime}\cap K\neq\emptyset\}$.

Report issue for preceding element

3.2 Finite element spaces and projectors

Report issue for preceding element

Given $m\in\mathbb{N}_{0}$ and $h>0$, we denote by $\mathbb{P}^{m}(\mathscr{T}_{h})$ the space of (possibly discontinuous) scalar functions that are polynomials of degree at most $m$ on each simplex $T\in\mathscr{T}_{h}$, and set $\mathbb{P}^{m}_{c}(\mathscr{T}_{h})\coloneqq\mathbb{P}^{m}(\mathscr{T}_{h})\cap C^{0}(\overline{\Omega})$. Then, given $k\in\mathbb{N}$ and $\ell\in\mathbb{N}_{0}$, we denote by

Report issue for preceding element