A comparative study of dynamic models for gravity-driven particle-laden flows

TL;DR

This study compares two dynamic models for viscous particle-laden flows down inclined surfaces, analyzing their similarities, differences, and applicability through theoretical derivations and numerical simulations.

Contribution

It provides a direct comparison of diffusive flux and suspension balance models using lubrication theory and simulations, highlighting their similarities and differences.

Findings

Fairly good agreement between models in flux profiles

Differences mainly at the free surface

Solutions differ by about 10% in certain conditions

Abstract

The dynamics of viscous thin-film particle-laden flows down inclined surfaces are commonly modeled with one of two approaches: a diffusive flux model or a suspension balance model. The diffusive flux model assumes that the particles migrate via a diffusive flux induced by gradients in both the particle concentration and the effective suspension viscosity. The suspension balance model introduces non-Newtonian bulk stress with shear-induced normal stresses, the gradients of which cause particle migration. Both models have appeared in the literature of particle-laden flow with virtually no comparison between the two models. For particle-laden viscous flow on an incline, in a thin-film geometry, one can use lubrication theory to derive a compact dynamic model in the form of a $2\times 2$ system of conservation laws. We can then directly compare the two theories side by side by looking at similarities and differences in the flux functions for the conservation laws, and in exact and numerical simulations of the equations. We compare the flux profiles over a range of parameters, showing fairly good agreement between the models, with the biggest difference involving the behavior at the free surface. We also consider less dense suspensions at lower inclination angles where the dynamics involve two shock waves that can be clearly measured in experiments. In this context the solutions differ by no more than about 10%, suggesting that either model could be used for this configuration.