Low-Reynolds-number instabilities in three-layer flow down an inclined wall

Abstract

The finite wavelength instability of viscosity-stratified three-layer flow down an inclined wall is examined for small but finite Reynolds numbers. It has previously been demonstrated using linear theory that three-layer zero-Reynolds-number instabilities can have growth rates that are orders of magnitude larger than those that arise in two-layer structures. Although the layer configurations yielding large growth instabilities have been well characterized, the physical origin of the three-layer inertialess instability remains unclear. Using analytic, numerical and experimental techniques, we investigate the origin and evolution of these instabilities. Results from an energy equation derived from linear theory reveal that interfacial shear and Reynolds stresses contribute to the energy growth of the instability at finite Reynolds numbers, and that this remains true in the limit of zero Reynolds number. This is thus a rare example that demonstrates how the Reynolds stress can play an important role in flow instability, even when the Reynolds number is vanishingly small. Numerical solutions of the Navier–Stokes equations are used to simulate the nonlinear evolution of the interfacial deformation, and for small amplitudes the predicted wave shapes are in excellent agreement with those obtained from linear theory. Further comparisons between simulated interfacial deformations and linear theory reveal that the linear evolution equations are surprisingly accurate even when the interfaces are highly deformed and nonlinear effects are important. Experimental results obtained using aqueous gelatin systems exhibit large wave growth and are in agreement with both the theoretical predictions of small-amplitude behaviour and the nonlinear simulations of the large-amplitude behaviour. Quantitative agreement is confounded owing to water diffusion driven by differences in gelatin concentration between the layers in experiments. However, the qualitative agreement is sufficient to confirm that the correct mechanism for the experimental instability has been determined.