Evolution and stability of falling liquid films with thermocapillary effects

In the first part, we further investigate the Benney equation in its validity domain in the case of periodically inhomogeneous heating in the streamwise direction. It induces steady-state deformations of the free surface with increased transfer rate in regions where the film is thinner, and also in average. The inhomogeneities of the heating also modify the nature of travelling wave solutions at moderate temperature gradients and allows for suppressing wave motion at larger ones.

Moreover, large temperature gradients (for instance positive ones) in the streamwise direction produce large local film thickening that may in turn become unstable with respect to transverse disturbances such that the flow may organize in rivulet-like structures. The mechanism of such instability is elucidated via an energy analysis. The main features of the rivulet pattern are described experimentally and recovered by direct numerical simulations.

In the second part, various models are obtained, which are valid for larger Reynolds numbers than the Benney equation and account for second-order viscous and inertial effects. We then elaborate a strategy to select the optimal model in terms of linear stability properties and existence of nonlinear solutions (solitary waves), for the widest possible range of parameters. This model -- called reduced model -- is a system of three coupled evolution equations for the local film thickness h, the local flow rate q and the surface temperature Ts. Solutions of this model indicate that the interaction of the hydrodynamic and thermocapillary modes is non-trivial, especially in the region of large-amplitude solitary waves.

Finally, the three-dimensional evolution of the solutions of the reduced model in the presence of periodic forcing and noise compares favourably with available experimental data in isothermal conditions and with direct numerical simulations in non-isothermal conditions.

------------------------------------------------

Cette thèse analyse la dynamique d'un film mince s'écoulant le long d'une paroi chauffée. Le chauffage crée des gradients de tension superficielle qui induisent des tensions thermocapillaires à la surface libre, altérant ainsi la stabilité et l'évolution du film. Grâce à la cohérence de l'écoulement assurée par la viscosité, deux approches permettant de réduire la dimensionnalité du problème original sont habituellement considérées suivant le débit (mesuré par le nombre de Reynolds): l'approximation asymptotique dite `longues ondes' pour les faibles nombres de Reynolds et l'approximation `intégrale couche limite' pour les nombres de Reynolds modérés. Cependant, la première approximation souffre de singularités et la dernière de prédictions imprécises du seuil de stabilité des ondes hydrodynamiques à la surface du film. Le but de cette thèse est donc double: dans une première partie, il s'agit de déterminer, de manière quantitative, la validité de l'équation d'évolution `longues ondes' (ou équation de Benney) pour l'épaisseur du film h, en y incluant l'effet thermocapillaire; et dans une seconde partie, il s'agit d'améliorer l'approche `intégrale couche limite' en combinant un développement en gradients avec une méthode aux résidus pondérés.

Dans la première partie, nous étudions l'équation de Benney, dans son domaine de validité, dans le cas d'un chauffage inhomogène et périodique dans la direction de l'écoulement. Cela induit des déformations permanentes de la surface libre avec un accroissement du transfert de chaleur dans les régions où le film est plus mince, mais aussi en moyenne. Un chauffage inhomogène modifie également la nature des solutions d'ondes progressives pour des gradients de températures modérés et conduit même à leur suppression pour des gradients de températures plus importants. De plus, ceux-ci, lorsqu'ils sont par exemple positifs le long de l'écoulement, produisent des épaississements localisés du film qui peuvent à leur tour devenir instables par rapport à des perturbations suivant la direction transverse à l'écoulement. Ce dernier s'organise alors sous forme d'une structure en rivulets. Le mécanisme de cette instabilité est élucidé via une analyse énergétique des perturbations. Les principales caractéristiques des structures en rivulets sont décrites expérimentalement et retrouvées par l'intermédiaire de simulations numériques.

Dans la seconde partie, nous dérivons une famille de modèles valables pour des nombres de Reynolds plus grands que l'équation de Benney, qui prennent en compte les effets visqueux et inertiels du second ordre. Nous élaborons ensuite une stratégie pour sélectionner le modèle optimal en fonction de ses propriétés de stabilité linéaire et de l'existence de solutions non-linéaires (ondes solitaires), et ce pour la gamme de paramètres la plus large possible. Ce modèle -- appelé modèle réduit -- est un système de trois équations d'évolution couplées pour l'épaisseur locale de film h, le débit local q et la température de surface Ts. Les solutions de ce modèle indiquent que l'interaction des modes hydrodynamiques et thermocapillaires n'est pas triviale, spécialement dans le domaine des ondes solitaires de grande amplitude. Finalement, l'évolution tri-dimensionnelle des solutions du modèle réduit en présence d'un forçage périodique ou d'un bruit se compare favorablement aux données expérimentales disponibles en conditions isothermes, ainsi qu'aux simulations numériques directes en conditions non-isothermes