Abstract:We propose a finite element discretization for the steady, generalized Navier-Stokes equations for fluids with shear-dependent viscosity, completed with inhomogeneous Dirichlet boundary conditions and an inhomogeneous divergence constraint. We establish (weak) convergence of discrete solutions as well as a priori error estimates for the velocity vector field and the scalar kinematic pressure. Numerical experiments complement the theoretical findings.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is the finite - element discretization problem of the steady - state generalized Navier - Stokes equations for non - Newtonian fluids with shear - dependent viscosity. Specifically, the paper focuses on the finite - element discretization method under inhomogeneous Dirichlet boundary conditions and inhomogeneous divergence constraints. The main objectives include:
1. **Establish the existence and stability of the discrete solution**:
- Prove the existence of the discrete solution.
- Establish a priori energy bounds for the discrete solution.
2. **Convergence analysis**:
- Prove the (weak) convergence of the discrete solution to the solution of the continuous problem in an appropriate weak form.
3. **Error estimation**:
- Derive a priori error estimates for the velocity vector field and the dynamic pressure.
The paper addresses these problems by introducing specific assumptions and mathematical tools, such as the Fortin interpolation operator, the theory of locally coercive operators, etc. In particular, for shear - thickening and shear - thinning fluids, certain small - quantity assumptions need to be made for the data to ensure the existence and convergence of the discrete solution.
### Formulas and symbol explanations
- **Basic equations**:
\[
-\text{div}\, S(Dv)+\text{div}(v\otimes v)+\nabla q = f\quad\text{in }\Omega
\]
where \(\Omega\subseteq\mathbb{R}^d\) is a bounded polyhedral Lipschitz domain, \(f:\Omega\rightarrow\mathbb{R}^d\) represents the external body force, and \(v:\Omega\rightarrow\mathbb{R}^d\) and \(q:\Omega\rightarrow\mathbb{R}\) represent the velocity vector field and the dynamic pressure, respectively.
- **Boundary conditions**:
\[
v = g_2\quad\text{on }\partial\Omega
\]
\[
\text{div}\, v = g_1\quad\text{in }\Omega
\]
- **Extra stress tensor**:
\[
S:\mathbb{R}^{d\times d}\rightarrow\mathbb{R}^{d\times d}_{\text{sym}}
\]
satisfying \((p,\delta)\)-structure, that is, there exist constants \(C_0,C_1 > 0\) such that:
\[
(S(A)-S(B)):(A - B)\geq C_0\phi_{|A^{\text{sym}}|}(|A^{\text{sym}}-B^{\text{sym}}|)
\]
\[
|S(A)-S(B)|\leq C_1\phi'_{|A^{\text{sym}}|}(|A^{\text{sym}}-B^{\text{sym}}|)
\]
where \(\phi=\phi_{p,\delta}\) is defined as:
\[
\phi(t)=\int_0^t\phi'(s)\,ds,\quad\phi'(t)=(\delta + t)^{p - 2}t
\]
- **Discrete weak form**:
\[
(S(Dv_h),Dz_h)_\Omega + b(v_h,v_h,z_h)-(q_h,\text{div}\,z_h)_\Omega=(f,z_h)_\Omega\quad\forall z_h\in V_h
\]
\[
(\text{div}\,v_h,z_h)_\Omega=(g_1,z_h)_\Omega\quad\forall z_h\in Y_h
\]
\[
v_h = g_{2h}