Error analysis for a fully-discrete finite element approximation of the unsteady $p(\cdot,\cdot)$-Stokes equations
Luigi C. Berselli, Alex Kaltenbach, Seungchan Ko
TL;DR
This work develops and analyzes a fully discrete finite element method for the unsteady $p(\cdot,\cdot)$-Stokes equations, using backward Euler in time and conforming, discretely inf-sup stable spaces in space. It leverages variable-exponent and generalized-function frameworks, including $L^{p(\cdot)}$ and $W^{1,p(\cdot)}$, Bochner–Nikolskii and Calderón spaces, to establish a priori error estimates under fractional regularity assumptions on the velocity and pressure; crucial toolchains include the non-linear operator pair $(\mathbf{S},\mathbf{F},\mathbf{F}^*)$ and shifted $N$-functions. The paper proves stability and interpolation results for FE projections, derives fractional interpolation error bounds, and culminates in a main theorem giving rate-decay for velocity and nonlinear stresses with respect to $\tau$ and $h$, plus auxiliary results under a pressure regularity framework. Numerical experiments corroborate the theoretical rates, particularly demonstrating quasi-optimal convergence for $p^{-}\ge 2$, and highlight the practical viability of the approach for non-Newtonian, variable-exponent fluids described by the unsteady $p(\cdot,\cdot)$-Stokes model. The results provide rigorous guidance for accurately simulating smart fluids with spatially and temporally varying rheology in computational contexts.
Abstract
In this paper, we examine a fully-discrete finite element approximation of the unsteady $p(\cdot,\cdot)$-Stokes equations ($i.e.$, $p(\cdot,\cdot)$ is time- and space-dependent), employing a backward Euler step in time and conforming, discretely inf-sup stable finite elements in space. More precisely, we derive error decay rates for the vector-valued velocity field imposing fractional regularity assumptions on the velocity and the kinematic pressure. In addition, we carry out numerical experiments that confirm the optimality of the derived error decay rates in the case $p(\cdot,\cdot)\ge 2$.
