Pulsatile Flows for Simplified Smart Fluids with Variable Power-Law: Analysis and Numerics

TL;DR

This paper analyzes the flow of variable-exponent non-Newtonian fluids in pipes, proving existence of time-periodic solutions, providing explicit benchmarks, and supporting practical applications with numerical methods.

Contribution

It introduces a unified approach for analyzing time-periodic flows of $p(ar{x})$-fluids, including explicit solutions and a constructive numerical proof, extending prior results for Navier-Stokes and $p$-fluids.

Findings

Existence of time-periodic solutions with specified flow-rate or pressure-drop.

Explicit benchmark solutions for smart fluids like electro-rheological fluids.

Numerical experiments confirming theoretical results.

Abstract

We study the fully-developed, time-periodic motion of a shear-dependent non-Newtonian fluid with variable exponent rheology through an infinite pipe $\Omega:= \mathbb{R}\times \Sigma\subseteq \mathbb{R}^d$, $d\in \{2,3\}$, of arbitrary cross-section $\Sigma\subseteq \mathbb{R}^{d-1}$. The focus is on a generalized $p(\cdot)$-fluid model, where the power-law index is position-dependent (with respect to $\Sigma$), $\textit{i.e.}$, a function $p\colon \Sigma\to (1,+\infty)$. We prove the existence of time-periodic solutions with either assigned time-periodic flow-rate or pressure-drop, generalizing known results for the Navier-Stokes and for $p$-fluid equations. In addition, we identify explicit solutions, relevant as benchmark cases, especially for electro-rheological fluids or, more generally, $\textit{`smart fluids'}$. To support practical applications, we present a fully-constructive existence proof for variational solutions by means of a fully-discrete finite-differences/-elements discretization, consistent with our numerical experiments. Our approach, which unifies the treatment of all values of $p(\overline{x})\in (1,+\infty)$, $\overline{x}\in \Sigma$, without requiring an auxiliary Newtonian term, provides new insights even in the constant exponent case. The theoretical findings are reviewed by means of numerical experiments.