Roll waves in mud

The stability of a viscoplastic fluid film falling down an inclined plane is explored, with the aim of determining the critical Reynolds number for the onset of roll waves. The Herschel-Bulkley constitutive law is adopted and the fluid is assumed two-dimensional and incompressible. The linear stability problem is described for an equilibrium in the form of a uniform sheet flow, when perturbed by introducing an infinitesimal stress perturbation. This flow is stable for very high Reynolds numbers because the rigid plug riding atop the fluid layer cannot be deformed and the free surface remains flat. If the flow is perturbed by allowing arbitrarily small strain rates, on the other hand, the plug is immediately replaced by a weakly yielded 'pseudo-plug' that can deform and reshape the free surface. This situation is modelled by lubrication theory at zero Reynolds number, and it is shown how the fluid exhibits free-surface instabilities at order-one Reynolds numbers. Simpler models based on vertical averages of the fluid equations are evaluated, and one particular model is identified that correctly predicts the onset of instability. That model is used to describe nonlinear roll waves.