We consider the dynamics of a thin liquid film falling down a uniformly heated wall. The heating sets up surface tension gradients that induce thermocapillary stresses on the free surface, thus affecting the evolution of the film. We model this thermocapillary flow by using a gradient expansion combined with a Galerkin projection with polynomial test functions for both velocity and temperature fields. We obtain equations for the evolution of the velocity and temperature amplitudes at first- and second-order in the expansion parameter. These equations are fully compatible close to criticality with the Benney long-wave expansion. Models of reduced dimensionality for the evolution of the local film thickness, flow rate and interfacial temperature only, are proposed.

We analyse the regularized reduced model derived in Part 1 (Ruyer-Quil et al. 2005). Our investigation is two-fold: (i) we demonstrate that the linear stability properties of the model are in good agreement with the Orr–Sommerfeld analysis of the linearized Navier–Stokes/energy equations; (ii) we show the existence of nonlinear solutions, namely single-hump solitary pulses, for the widest possible range of parameters. We also scrutinize the influence of Reynolds, Prandtl and Marangoni numbers on the shape, speed, flow patterns and temperature distributions for the solitary waves obtained from the regularized model. The hydrodynamic and Marangoni instabilities are seen to reinforce each other in a non-trivial manner. The transport of heat by the flow has a stabilizing effect for small-amplitude waves but promotes the instability for large-amplitude waves when a recirculating zone is present. Nevertheless, in this last case, by increasing the shear in the bulk and thus the viscous dissipation, increasing the Prandtl number decreases the amplitude and speed of the waves.

We present new insights and results for the problem of a film falling down a heated wall: (i) treatment of a mixed heat flux boundary condition on the substrate; (ii) development of a long-wave theory for large Péclet numbers; (iii) refined treatment of the energy equation based on a high-order Galerkin projection in terms of polynomial test functions which satisfy all boundary conditions; (iv) time-dependent computations for the free-surface height and interfacial temperature; (v) numerical solution of the full energy equation; (vi) demonstration of the existence of a thermal boundary layer at the front stagnation point of a solitary pulse; (vii) development of models that prevent negative temperatures and are in good agreement with the numerical solution of the full energy equation.

We investigate the dynamics of a thin laminar liquid film flowing under gravity down the lower wall of an inclined channel when turbulent gas flows above the film. The solution of the full system of equations describing the gas–liquid flow faces serious technical difficulties. However, a number of assumptions allow isolating the gas problem and solving it independently by treating the interface as a solid wall. This permits finding the perturbations to pressure and tangential stresses at the interface imposed by the turbulent gas in closed form. We then analyse the liquid film flow under the influence of these perturbations and derive a hierarchy of model equations describing the dynamics of the interface, i.e. boundary-layer equations, a long-wave model and a weakly nonlinear model, which turns out to be the Kuramoto–Sivashinsky equation with an additional term due to the presence of the turbulent gas. This additional term is dispersive and destabilising (for the counter-current case; stabilizing in the co-current case). We also combine the long-wave approximation with a weighted-residual technique to obtain an integral-boundary-layer approximation that is valid for moderately large values of the Reynolds number. This model is then used for a systematic investigation of the flooding phenomenon observed in various experiments: as the gas flow rate is increased, the initially downward-falling film starts to travel upwards while just before the wave reversal the amplitude of the waves grows rapidly. We confirm the existence of large-amplitude stationary waves by computing periodic travelling waves for the integral-boundary-layer approximation and we corroborate our travelling-wave results by time-dependent computations.

We consider the static and dynamic behaviour of two-dimensional droplets on inclined heterogeneous substrates. We utilize an evolution equation for the droplet thickness based on the long-wave approximation of the Stokes equations in the presence of slip. Through a singular perturbation procedure, evolution equations for the location of the two moving fronts are obtained under the assumption of quasi-static dynamics. The deduced equations, which are verified by direct comparisons with numerical solutions to the governing equation, are scrutinized in a variety of dynamic and equilibrium settings. For example, we demonstrate the possibility for stick–slip dynamics, substrate-induced hysteresis, the uphill motion of the droplet for sufficiently strong chemical gradients and the existence of a critical inclination angle beyond which the droplet can no longer be supported at equilibrium. Where possible, analytical expressions are obtained for various quantities of interest, which are also verified by appropriate numerical experiments.