International Journal of Heat and Mass Transfer, Vol.90, 689-701, 2015
Falling liquid films on a slippery substrate with Marangoni effects
This paper investigates the dynamics of liquid films of mean thickness' h(0) flowing down a uniformly heated slippery substrate. The dimensionless slippery length beta is small and is assumed to have the same order as the film parameter epsilon = h(0)/l where l is a typical wave length. A weighted-residual model and a Benney-type model are derived to discuss the influences of slippery length beta, Reynolds number Re, Marangoni number M and Blot number Bi on the dynamics of the film. It is found that, when Reynolds number is small (Re < 1) and the wall is non-slip, the Benney-type model, weighted-residual model are in good agreements with the linearized Navier Stokes equations (LNS). However, when the wall is slip or the Reynolds number is moderate, the Benney-type model does not agree with LNS. Our results show that, for a small slippery effect at the substrate, the weighted-residual model compares well with LNS at moderate Re. It demonstrates that the weighted-residual model is more reasonable than the Benney-type model. A critical Reynolds number obtained by the weighted-residual model for this system is Re-c = 35(1-beta)(4)/(6(7-67 beta+159 beta(2)) (Wek(2) + cot alpha - 3MBi/2(1+Bi-2)) (alpha is the inclined angle and We is the Weber number). Linear stability analysis shows that, the instability is enhanced as the slippery length,beta and the Reynolds number Re increase. In the long-wave range, the linear wave speed c =1 + 2 beta. When the wall is heated, the instability is reinforced. The instability is impeded when the substrate is cooled. Nonlinear analysis by the traveling wave solutions shows that, the wave speed and wave height of traveling waves are promoted by the slippery and fluid inertial effects. Nonlinear study shows that the traveling wave speed is promoted by a heated substrate and is reduced by a cooled substrate. It shows that the speed of fast waves is larger than the linear wave speed and the higher the wave is the faster it travels. Direct numerical simulation of the weighted residual model supports the linear stability analysis and the nonlinear analysis of traveling wave solutions. (C) 2015 Elsevier Ltd. All rights reserved.