On the fractal dimension of non-Newtonian Hele-Shaw flow subject to Saffman-Taylor instability

TL;DR

This paper presents a novel numerical method to simulate non-Newtonian Hele-Shaw flow with shear-thinning fluids, revealing that shear-thinning alters pattern morphology by reducing fractal dimension compared to Newtonian fluids.

Contribution

The study introduces a diffusion-limited aggregation-based simulation incorporating shear-rate-dependent viscosity for non-Newtonian fluids in Hele-Shaw flow.

Findings

Shear-thinning reduces the fractal dimension of flow patterns.

Pattern morphology is significantly affected by non-Newtonian fluid properties.

The method effectively models history-dependent shear-thinning effects.

Abstract

We introduce a discrete numerical method based on the diffusion-limited aggregation (DLA) approach to simulate two-fluid Hele-Shaw flow subject to the Saffman-Taylor interfacial instability, in the case where the displaced fluid is non-Newtonian. Focusing on fluids for which the most relevant non-Newtonian aspect of the thin-gap flow is shear-thinning, we introduce a history-dependent aspect into the algorithm, modeling shear-rate-dependent fluid viscosity. The main finding is that the morphology of the emerging patterns, characterized by the fractal dimension, is modified in a nontrivial manner by the shear-thinning nature of the displaced fluid. In particular, we consistently find that shear-thinning leads to the formation of patterns characterized by a smaller fractal dimension, compared to the corresponding Newtonian fluid.