Surfactant and gravity dependent instability of two-layer channel flows: Linear theory covering all wave lengths

TL;DR

This paper presents a comprehensive linear stability analysis of two-layer Couette flows with surfactant and gravity effects, deriving explicit growth rate formulas and exploring how parameters influence flow stability.

Contribution

It provides a detailed analytical and numerical investigation of the stability landscape, including explicit growth rates and bifurcation points, for flows with surfactant and gravity effects.

Findings

Identified two normal mode branches: robust and surfactant-dependent.

Derived explicit formulas for growth rates in the inertia-less limit.

Mapped stability thresholds in the (Ma, Bo)-plane and analyzed parameter effects.

Abstract

A linear stability analysis of a two-layer plane Couette flow of two immiscible fluid layers with different densities, viscosities and thicknesses, bounded by two infinite parallel plates moving at a constant relative velocity to each other, with an insoluble surfactant along the interface and in the presence of gravity is carried out. The normal modes approach is applied to the equations governing flow disturbances. These equations, together with boundary conditions at the plates and the interface, yield a linear eigenvalue problem. When inertia is neglected velocity amplitudes are linear combinations of hyperbolic functions, and a quadratic dispersion equation for the complex growth rate is obtained where coefficients depend on the aspect ratio, the viscosity ratio, the basic velocity shear, the Marangoni number Ma that measures the effects of surfactant, and the Bond number Bo that measures the influence of gravity. An extensive investigation is carried out that examines the stabilizing or destabilizing influences of these parameters. There are two continuous branches of the normal modes: a robust branch that exists even with no surfactant, and a surfactant branch that vanishes when Ma $\downarrow 0$. Due to the availability of the explicit forms for the growth rates, in many instances the numerical results are corroborated with analytical asymptotics. For the less unstable branch, a mid-wave interval of unstable wavenumbers (Halpern and Frenkel (2003)) sometimes co-exists with a long-wave one. We study the instability landscape, determined by the threshold curve of the long-wave instability and the critical curve of the mid-wave instability in the (Ma, Bo)-plane. The changes of the extremal points of the critical curves with the variation of the other parameters, such as the viscosity ratio, and the extrema bifurcation points are investigated.