Stability of gravity-driven free-surface flow past a deformable solid: The role of depth-dependent modulus

Phys Rev E. 2020 Apr;101(4-1):043107. doi: 10.1103/PhysRevE.101.043107.

Abstract

The linear stability of a Newtonian liquid layer flowing down an inclined plane lined with a deformable linear elastic solid characterized by a continuously varying modulus is analyzed in this study. A low-wave-number asymptotic analysis is performed to obtain an analytical expression for the complex wavespeed which shows striking similarity with the earlier results of Sahu and Shankar [Sahu and Shankar, Phys. Rev. E 94, 013111 (2016)10.1103/PhysRevE.94.013111] for gravity-driven flow of Newtonian fluid past a solid bilayer having constant shear modulus (in each layer) lined on a rigid inclined plane. This shows that a deformable solid layer having a continuously varying shear modulus can be treated as a generalization of a system having multiple solid layers of constant shear modulus. Also, in the low-wave-number limit, we show that the stability of the free surface is governed by the value of effective shear modulus G_{eff}, and not by the detailed spatial variation of the modulus. Here the effective shear modulus (H/G_{eff}=∫_{1}^{1+H}1/[E_{0}+E[over ¯](z)]dz, where [E_{0}+E[over ¯](z)] represents the modulus gradient function) characterizes the overall modulus of the elastic solid, which is obtained by treating the continuous variation to be the limit of the arrangement of solid layers of infinitesimal thickness (each having a constant shear modulus) in a series. At finite wave numbers, we show that the free-surface and the liquid-solid interface become unstable as we increase the value of Γ, where Γ indicates the ratio of viscous stresses in the fluid to elastic stresses in the solid. When the system is analysed for different types of spatial modulus variations, we find results similar to those of Gkanis and Kumar [Gkanis and Kumar, Phys. Rev. E 73, 026307 (2006)10.1103/PhysRevE.73.026307], i.e., when we have two different configurations of the shear modulus function that have the same spatially averaged modulus, but have different values at the interface, the system is more stable for the configuration having higher shear modulus at the liquid-solid interface. In a similar manner, when we examined systems having the same shear modulus at the liquid-solid interface and same average modulus, the more stable case is the one which has a higher value of shear modulus just below the interface. Thus the use of deformable solids with a depth-dependent modulus potentially offers more control in the passive manipulation of the instabilities.