On the Determination of the Yield Surface within the Flow of Yield Stress Fluids using Computational Fluid Dynamics

TL;DR

This paper uses Computational Fluid Dynamics with OpenFOAM to model and analyze the three-dimensional flow behavior of yield stress fluids in open channels, comparing results with experiments and theoretical predictions.

Contribution

It introduces a 3D CFD model for yield stress fluids in open channels, validating it against experimental and theoretical data, and identifies flow zones within the fluid.

Findings

The 3D model accurately reproduces flow patterns.

Flow zones such as rigid and sheared regions are detected.

Numerical results align with experimental and theoretical data.

Abstract

A part of non-Newtonian fluids are yield stress fluids. They require a minimum stress to flow. Below this minimum value, yield stress fluids remain solid. To date, 1D and 2D numerical models have been used predominantly to study free surface flows. However, some phenomena have three-dimensional behaviour such as the appearance of the limit between the liquid regime and the solid regime. Here the aim is to use a Computational Fluid Dynamics (CFD) to reproduce the properties of the free surface flow of yield stress fluids in an open channel. Modelling the behaviour of the yield stress fluid is also expected. The numerical study is driven with the software OpenFOAM. Numerical outcomes are compared with experimental results from model experiment and theorical predictions based on the rheological constitutive law. The 3D model is validated by evaluating its capacity to reproduce reliably flow patterns. The depth, the local velocity and the stress are quantified for different numerical configurations (grid level, rheological parameters). Then numerical results are used to detect the presence of rigid and sheared zones within the flow.