Cette étude concerne la physique des écoulements convectifs résultant d'une instabilité d'évaporation de fluides binaires. Ce problème a de nombreuses applications, l'enrobage par centrifugation, le dépôt de films, les caloducs, etc, pour lesquels le changement de phase et la convection jouent un rôle prépondérant dans la conception et la qualité des procédés. Le système physique étudié est un mélange liquide sous sa propre vapeur, confiné par deux plaques conductrices de chaleur et des bords latéraux isolants. Les plaques sont utilisées pour appliquer un gradient thermique. Aucun gradient de concentration n'est imposé au système. Ces gradients sont induits par les différentes vitesses d'évaporation des composés. Dans ce système, il est important de comprendre comment la dynamique des fluides et les transferts de masse et de chaleur entrent en compétition pour la formation de structures. Le principal objectif de ce travail est d'identifier les conditions pour que le système évolue d'un état conductif vers un état de convection lorsque le gradient vertical de température dépasse une certaine valeur critique.Dans le système, la convection s'installe par trois mécanismes distincts : évaporation, gradients de densité et gradients de tension interfaciale. Trois forces convectives s'opposent aux effets de diffusion qui tendent à garder le système en état conductif. Le seuil d'apparition de la convection dépend de quelques variables, comme les dimensions du contenant, les propriétés thermophysiques des phases liquide et vapeur, la fraction massique, et les caractéristiques de perturbations. L'effet de chacune de ces variables sur le seuil est étudié en présence ou non de gravité.Pour représenter la physique, un modèle mathématique non linéaire complet est développé, basé sur les conservations de quantité de mouvement, d'énergie et de masse dans chaque phase avec les conditions aux limites appropriées. Le fluide binaire est composé de deux alcools légers comme l'éthanol et le sec-butanol. Dans les équations du modèle, la masse volumique ainsi que la tension interfaciale sont fonctions à le fois de la température et de la concentration. Pour la recherche du seuil de transition, les équations sont linéarisées autour d'un état de base connu. Dans notre cas, il s'agit de l'état conductif. Le système d'équations linéaires résultant est résolu par une méthode de collocation spectrale Chebyshev.Nous obtenons quatre résultats principaux. Premièrement, dans un système multi-composants sans gravitation, une instabilité n'apparaît que lorsque le système est chauffé du côté de la phase vapeur contrairement à un système mono-composant. Cela implique que, si on souhaite éviter les instabilités, il vaut mieux un apport de chaleur par la phase liquide en cas de processus d'évaporation en couches minces ou en micro-gravité.Deuxièmement, en présence de gravité, un système multi-composants peut devenir instable quelle que soit la direction du chauffage. Si la convection thermique est négligeable, alors nous montrons que le chauffage par la phase vapeur est la configuration la plus instable. Sinon, les deux modes de chauffage sont à même de produire une instabilité. Ce résultat implique que le gradient thermique appliqué doit être inférieur à une valeur seuil pour éviter les instabilités quelle que soit la direction du chauffage.Troisièmement, lorsque l'instabilité apparaît en absence de gravité, des structures n'apparaitront pas dans le cas de fluide pur mais apparaitront dans le cas d'un fluide multi-composants. De même, des structures apparaitront en présence de gravité en fonction du facteur d'aspect du confinement. Les facteurs d'aspect peuvent être choisis pour éviter des structures multi-cellulaires même en cas d'apparition d'instabilités durant l'évaporation.Enfin, des structures oscillantes ne sont pas prédites de façon générale malgré les effets opposés des convections solutale et thermique dans le problème d'évaporation.

Cette étude concerne la physique des écoulements convectifs résultant d’une instabilité d’évaporation de fluides binaires. Ce problème a de nombreuses applications, l’enrobage par centrifugation, le dépôt de films, les caloducs, etc, pour lesquels le changement de phase et la convection jouent un rôle prépondérant dans la conception et la qualité des procédés. Le système physique étudié est un mélange liquide sous sa propre vapeur, confiné par deux plaques conductrices de chaleur et des bords latéraux isolants. Les plaques sont utilisées pour appliquer un gradient thermique. Aucun gradient de concentration n’est imposé au système. Ces gradients sont induits par les différentes vitesses d’évaporation des composés. Dans ce système, il est important de comprendre comment la dynamique des fluides et les transferts de masse et de chaleur entrent en compétition pour la formation de structures. Le principal objectif de ce travail est d’identifier les conditions pour que le système évolue d’un état conductif vers un état de convection lorsque le gradient vertical de température dépasse une certaine valeur critique.Dans le système, la convection s’installe par trois mécanismes distincts : évaporation, gradients de densité et gradients de tension interfaciale. Trois forces convectives s’opposent aux effets de diffusion qui tendent à garder le système en état conductif. Le seuil d’apparition de la convection dépend de quelques variables, comme les dimensions du contenant, les propriétés thermophysiques des phases liquide et vapeur, la fraction massique, et les caractéristiques de perturbations. L’effet de chacune de ces variables sur le seuil est étudié en présence ou non de gravité.Pour représenter la physique, un modèle mathématique non linéaire complet est développé, basé sur les conservations de quantité de mouvement, d’énergie et de masse dans chaque phase avec les conditions aux limites appropriées. Le fluide binaire est composé de deux alcools légers comme l’éthanol et le sec-butanol. Dans les équations du modèle, la masse volumique ainsi que la tension interfaciale sont fonctions à le fois de la température et de la concentration. Pour la recherche du seuil de transition, les équations sont linéarisées autour d’un état de base connu. Dans notre cas, il s’agit de l’état conductif. Le système d’équations linéaires résultant est résolu par une méthode de collocation spectrale Chebyshev.Nous obtenons quatre résultats principaux. Premièrement, dans un système multi-composants sans gravitation, une instabilité n’apparaît que lorsque le système est chauffé du côté de la phase vapeur contrairement à un système mono-composant. Cela implique que, si on souhaite éviter les instabilités, il vaut mieux un apport de chaleur par la phase liquide en cas de processus d’évaporation en couches minces ou en micro-gravité.Deuxièmement, en présence de gravité, un système multi-composants peut devenir instable quelle que soit la direction du chauffage. Si la convection thermique est négligeable, alors nous montrons que le chauffage par la phase vapeur est la configuration la plus instable. Sinon, les deux modes de chauffage sont à même de produire une instabilité. Ce résultat implique que le gradient thermique appliqué doit être inférieur à une valeur seuil pour éviter les instabilités quelle que soit la direction du chauffage.Troisièmement, lorsque l’instabilité apparaît en absence de gravité, des structures n’apparaitront pas dans le cas de fluide pur mais apparaitront dans le cas d’un fluide multi-composants. De même, des structures apparaitront en présence de gravité en fonction du facteur d’aspect du confinement. Les facteurs d’aspect peuvent être choisis pour éviter des structures multi-cellulaires même en cas d’apparition d’instabilités durant l’évaporation.Enfin, des structures oscillantes ne sont pas prédites de façon générale malgré les effets opposés des convections solutale et thermique dans le problème d’évaporation. / This study focuses on understanding the physics of the convective flow resulting from evaporative instability in binary mixtures. This problem has wide applications in spin coating, film deposition, heat pipes, etc. where phase change and convection play a very important role in the design process and also final quality of the product. The physical system of interest consists of a liquid mixture underlying its own vapor sandwiched between two conducting plates with insulated sidewalls in a closed container. The conducting plates are used to apply a vertical temperature gradient while there is no applied concentration gradient in the system. Concentration gradients are induced by the different evaporation rate of the components. In this system it is important to understand how the fluid dynamics and the heat and mass transfer interact competitively to form patterns. The main goal of this work is to identify the conditions for the system going from the conductive no-flow state to a convection state when the applied vertical temperature gradient exceeds a certain value called the critical value.In the system convection arises due to three distinct phenomena; evaporation, density gradients, and interfacial tension gradients. These convective forces are opposed by the diffusion effects that try to keep the system in the conductive no-flow state. The onset point depends upon several variables such as the dimensions of the container, thermo-physical properties of both liquid and vapor phases, mass fraction, and the characteristic of the disturbance given to the system. The effects of each of these variables on the onset point are investigated both in the presence and in the absence of gravity. To represent the physics a complete non-linear mathematical model is developed including momentum, energy, and mass balances in both phases with appropriate boundary conditions. The binary mixture is assumed to be made up of two low weight alcohols such as ethanol and sec-butanol. In the modeling equations the density and the interfacial tension are taken to be function of both temperature and concentration. To identify the onset point the non-linear equations are linearized around a known base state. In this case the base state is the conductive no-flow state. The resulting set of linear equations is solved using a spectral Chebyshev collocation method. Four major results arise from this work. First, in a multi-component system in the absence of gravity, an instability arises only when the system is heated from the vapor side as opposed to evaporation in a single-component. The implication is that evaporative processes in thin layers or in micro-gravity are best conducted with heat from the liquid side if instabilities are to be avoided.Second, in the presence of gravity, a multi-component system may become unstable no matter the direction of heating. If thermal buoyancy is negligible then it is shown in this study that heating from the vapor side is the unstable arrangement. Otherwise either heating style can produce an instability. This result means that the applied temperature difference must be kept below a threshold in order to avoid flow instabilities no matter the heating direction.Third, whenever instability occurs in the absence of gravity, patterns will not result in the case of a pure component but may result in the case of multi-components. Likewise, patterns will result when gravity is taken into account provided the aspect ratio of the container lies in a suitable range. As a result, aspect ratios can be chosen to avoid multi-cellular patterns even if convective flow instabilities arise during evaporation.Lastly, oscillations are not ordinarily predicted despite opposing effects of solutaland thermal convection in the evaporation problem.

Nous étudions numériquement la stabilité vis-à-vis de perturbations 3D d'un écoulement 2D d'un liquide maintenu par capillarité entre deux barreaux cylindriques coaxiaux isothermes et soumis à un flux thermique latéral. Les solutions numériques sont obtenues par méthode de collocation spectrale. Les écoulements stationnaires sont obtenus par méthode de Newton et une méthode d'Arnoldi est utilisée pour l'étude de la stabilité linéaire. La recherche a été menée sur une large gamme de nombres de Prandtl, allant de 0.001 à 100 et pour un rapport d'aspect égal à 2. Le mécanisme de déstabilisation est analysé en observant le taux de croissance de l'énergie de la perturbation. L'utilisation d'un nouvel outil d'analyse, le système adjoint, permet d'identifier les zones sensibles de l'écoulement à des perturbations impulsionnelles. La localisation des zones sensibles permet d'identifier les structures sensibles des écoulements stationnaires 2D. La structure d'écoulements 3D faiblement non linéaires a aussi été décrite.

Ce travail de thèse est consacré à l’étude des écoulements thermogravitaires et thermocapillaires induits par chauffage laser dans des couches liquides. Le chauffage d’un système à deux liquides donne naissance à deux effets thermiques. Le premier est dû à la variation de la masse volumique avec la température et le second à la variation de la tension interfaciale avec la température. Nous avons dans un premier temps étudié ces deux contributions séparément. En confrontant expériences, théorie et simulations numériques nous avons démontré que la morphologie des écoulements thermogravitaires dépend de l’épaisseur de la couche liquide ainsi que de la largeur du champ de température. Ensuite nous nous sommes intéressés à l’étude théorique et numérique de l’effet thermocapillaire. Cette étude a révélé qu’il est possible d’étudier les écoulements à partir des déformations d’interface induites par ces mêmes écoulements. Nous avons montré qu’il existe deux régimes de déformations en fonction du rapport de hauteurs et du signe de la variation de la tension interfaciale avec la température. Enfin, nous nous sommes intéressés à la compétition entre ces deux mécanismes. L’analyse des déformations d’interface et la comparaison avec un modèle à une dimension a montré qu’en fonction du rapport de hauteurs des couches liquides nous avons une transition d’un régime d’écoulements thermocapillaires vers une régime d’écoulements thermogravitaires. De plus, nous avons montré expérimentalement et numériquement qu’il est possible d’induire une instabilité d’origine thermogravitaire conduisant à la formation d’un pont liquide. / This thesis work is dedicated to thermocapillary and thermogravitary flows induced bylaser-heating in liquid layers. The laser-heating of a two fluids system induces two thermaleffects. The first effect come from the variation of the density with temperature andthe second one is due to the variation of the interfacial tension with the temperature.In a first part, we study separately these two mechanism. With a comparison betweenexperiments, theory and numerical simulations we demonstrated that the morphology ofthe flows depends on the thickness of the layer and the width of the temperature field.Then, we studied numerically and theoretically the thermocapillary effect. That revealedit is possible to understand the flows from the interface deformations. We highligth theseare two deformation regimes depending the thickness ratio and the sign of the interfacialtension gradient. Finally, we studied the competition between the two mechanism and theexperiments revealed a transition between two flow regimes. The comparison of these resultsand a theoritical model showed there is a transition from a thermocapillary regime to athermogravitary regime. In addition, we showed the possibility to induce an instibility bythermogravitary effect which can lead to the formation of a liquid bridge.

Cette étude se concentre sur la compréhension de la physique des instabilités dans différents écoulements de cisaillement, particulièrement la cavité entraînée et la cavité thermocapillaire, où l'écoulement d'un fluide incompressible est assuré soit par le mouvement d’une ou plusieurs parois, soit par des contraintes d’origine thermique.Un code spectral a été validé sur le cas très étudié de la cavité entrainée par une paroi mobile. Il est démontré dans ce cas que l'écoulement transit d'un régime stationnaire à un instationnaire au-delà d'une valeur critique du nombre de Reynolds. Ce travail est le premier à donner une interprétation physique de l'évolution non monotonique du nombre de Reynolds critique en fonction du facteur d'aspect. Lorsque le fluide est entraîné par deux parois mobiles, la cavité entraînée possède un plan de symétrie particulièrement sensible. Des solutions asymétriques peuvent être observés en plus de la solution symétrique au-dessus d'une certaine valeur du nombre de Reynolds. La transition oscillatoire entre la solution symétrique et les solutions asymétriques est expliquée physiquement par les forces en compétition. Dans le cas asymétrique, l'évolution de la topologie permet à l'écoulement de rester stationnaire avec l'augmentation du nombre de Reynolds. Lorsque l'équilibre est perdu une instabilité se manifeste par l'apparition d'un régime oscillatoire dans l'écoulement asymétrique.Dans une cavité thermocapillaire rectangulaire avec une surface libre, Smith et Davis prévoient deux types d'instabilités convectives thermiques: des rouleaux longitudinaux stationnaires et des ondes hydrothermales instationnaires. L'apparition de ses instabilités a été mis en évidence à plusieurs reprises expérimentalement et numériquement. Alors que les applications impliquent souvent plus d'une surface libre, il semble qu'il y ait peu de connaissances sur l'écoulement thermocapillaire entraînée avec deux surfaces libres. Un film liquide libre soumis à des contraintes thermocapillaires possède un plan de symétrie particulier comme dans le cas de la cavité entrainée par deux parois mobiles. Une étude de stabilité linéaire avec deux profils de vitesse pour le film liquide libre est présentée avec différents nombres de Prandtl. Au-delà d'un nombre de Marangoni critique, il est découvert que ces états de base sont sensibles à quatre types d'instabilités convectives thermiques qui peuvent conserver ou briser la symétrie du système. Les mécanismes qui permettent de prédire ces instabilités sont également découverts et interpréter en fonction de la valeur du nombre de Prandtl du fluide. La comparaison avec les travaux de Smith et Davis est faite. Une simulation numérique directe permet de valider les résultats obtenus avec l'étude de stabilité de linéaire. / This study focuses on the understanding of the physics of different instabilities in driven cavities, specifically the lid-driven cavity and the thermocapillarity driven cavity where flow in an incompressible fluid is driven either due to one or many moving walls or due to surface stresses that appear from surface tension gradients caused by thermal gradients. A spectral code is benchmarked on the well-studied case of the lid-cavity driven by one moving wall. In this case, It is shown that the flow transit form a steady regime to unsteady regime beyond a critical value of the Reynolds number. This work is the first to give a physical interpretation of the non-monotonic evolution of the critical Reynolds number versus the size of the cavity. When the fluid is driven by two facing walls moving in the same direction, the cavity possesses a plane of symmetry particularly sensitive. Thus, asymmetrical solutions can be observed in addition to the symmetrical solution above a certain value of the Reynolds number. The oscillatory transition between the symmetric solution and asymmetric solutions is explained physically by the forces in competition. In the asymmetric case, the change of the topology allows the flow to remain steady with increasing the Reynolds number. When the equilibrium is lost, an instability manifests by the appearance of an oscillatory regime in the asymmetric flow. In a rectangular cavity thermocapillary with a free surface, Smith and Davis found two types of thermal convective instabilities: steady longitudinal rolls and unsteady hydrothermal waves. The appearance of its instability has been highlighted repeatedly experimentally and numerically. While applications often involve more than a free surface, it seems that there is little knowledge about the thermocapillary driven flow with two free surfaces. A free liquid film possesses a particular plane of symmetry as in the case of the two-sided lid-driven cavity. A linear stability analysis for the free liquid film with two velocity profiles is presented with various Prandtl numbers. Beyond a critical Marangoni number, it is observed that these basic states are sensitive to four types of thermal convective instabilities, which can keep or break the symmetry of the system. Mechanisms that predict these instabilities are discovered and interpreted according to the value of the Prandtl number of the fluid. Comparison with the work of Smith and Davis is made. A direct numerical simulation is done to validate the results obtained with the linear stability analysis.

Ce travail porte sur l'étude de la thermodiffusion, ou effet Soret, par simulation numérique à l'échelle microscopique. Ce processus de transport croisé couple flux de masse et gradient thermique et est encore largement incompris. Pour cette étude, nous avons appliqué un algorithme de dynamique moléculaire hors équilibre à des mélanges de sphères de Lennard-Jones libres ou confinées. Après avoir testé la validité de nos simulations, nous avons montré que les résultats obtenus permettaient d'estimer la thermodiffusion dans ces fluides modèles à partir de corrélations simples sur les paramètres moléculaires. Cette démarche a été également validée sur des mélanges ternaires. Par ailleurs, les résultats de l'influence du milieu poreux sur la thermodiffusion ont montré la prépondérance des effets d'adsorption sur ceux liés au confinement. Cette influence restant faible dans la majorité des cas, exceptée pour les pores les plus fins et les plus attractifs.

Les ondes hydrothermales constituent un système modèle d'ondes non-linéaires dont nous étudions expérimentalement la transition vers le chaos spatio-temporel. A une dimension d'espace, par comparaison entre une cellule annulaire périodique et une cellule rectangulaire de taille finie, nous montrons que l'instabilité primaire en ondes ne survient dans ce dernier cas que lorsque la transition convectif/absolu est franchie. Nous isolons de même pour l'instabilité secondaire d'Eckhaus les régimes d'instabilité convective et absolue. A deux dimensions d'espace dans une cellule cylindrique, nous étudions la structuration de l'écoulement de base en rouleaux corotatifs, puis les différents modes d'instabilités en fonction de la différence de température appliquée et de la hauteur de fluide~: spirales d'Archimède, cibles, fleurs et rayons apparaissent chacuns par bifurcation de Hopf supercritique. Deux modes d'ondes hydrothermales sont alors distingués. Nous tentons d'expliquer l'existence de ces deux modes différents par une étude théorique et numérique: nous quantifions les effets de la courbure sur les instabilités à deux dimensions.

L’évaporation de gouttes de liquides contenant des particules non volatiles représente un phénomène largement présent dans la vie quotidienne, à l’image des traces laissées par le marc de café après séchage. L’étude de la morphologie des dépôts de particules présente un grand intérêt dans les domaines de la biologie et trouve de nombreuses applications dans l’industrie. De ce fait, elle a fait l’objet de nombreuses recherches durant les dernières décennies. Malgré les nombreuses récentes recherches sur les morphologies des dépôts de particules, les mécanismes les contrôlant restent encore non complétement expliqués. Certains facteurs influençant les morphologies des dépôts sont nombreux (température de substrats…) mais restent encore peu documentés dans la littérature. Cette étude expérimentale s’intéresse à l’influence de la température du substrat sur la morphologie des dépôts de nanoparticules après séchage de gouttes sessiles de liquides. L’augmentation de la température du substrat accélère le processus d’évaporation et entraine des morphologies de dépôts très différentes de celles obtenues sur des substrats à température ambiante. Dans cette étude, la microscopie combinée à la thermographie infrarouge et à l’interférométrie ont permis d’expliquer la dynamique de formation de dépôts. De plus, l’étude a permis d’analyser les effets d’autres paramètres sur la morphologie des dépôts, tel que la composition chimique du liquide composant les gouttes. / Evaporation of liquid droplets containing non-volatile solutes is an omnipresent phenomenon in daily life, e.g., coffee stains on solid surfaces. The study of pattern formation of the particles left after the evaporation of a sessile droplet has attracted the attention of many researchers during the past two decades due to the wide range of biological and industrial applications. Despite the significance of controlling the deposition morphology of droplets, the underlying mechanisms involved in pattern formation are not yet fully understood. There is a varied range of factors that affect the final deposition patterns and some, e.g., substrate temperature, are poorly studied in the literature. This experimental study investigates the effect of a wide range of substrate temperatures on the deposition patterns of nanoparticles from drying sessile droplets. Increasing substrate temperature and accelerating the drying process lead to the formation of the patterns not observed on non-heated substrates. This research elucidates the formation mechanisms of these patterns by optical microscopy, infrared thermography, and white light interferometry techniques. Furthermore, the combined effects of substrate temperature and other factors such as chemical composition of base fluid and particle size on the dried patterns are studied. The underlying mechanisms involved in the formation of the patterns influenced by the combined factors are also discussed and presented.

Plusieurs études tant numériques qu’expérimentales font état de la présence d’instabilités thermiques dans des films liquides chauffés uniformément par le bas pour des conditions aux limites et d’écoulements particuliers. La présence de ces instabilités modifiera les transferts thermiques associés. Le sujet de ce travail de thèse consiste à étudier numériquement les instabilités thermoconvectives d’un écoulement laminaire tridimensionnel de convection mixte dans un canal horizontal à surface libre. Les variations de la tension de surface avec la température (effet Marangoni ou effet thermocapillaire) sont prises en compte. Bien que d’un intérêt certain pour de nombreuses applications industrielles, cette situation a été très peu étudiée d’un point de vue académique dans la configuration considérée ici. Dans cette de configuration plusieurs types de structures thermoconvectives sont susceptibles d’apparaître. Lorsque les forces induites par les courants de convection naturelle, forcée et thermocapillaire sont du même ordre de grandeur, les premiers résultats montrent un développement des instabilités sous forme de rouleaux convectifs longitudinaux stationnaires semblables à ceux rencontrés pour des écoulements de type Poiseuille-Rayleigh-Bénard. A notre connaissance, c’est la première fois que l’écoulement de convection de type Poiseuille-Rayleigh-Bénard associé aux effets Marangoni est étudié. Le nombre et la distribution spatiale des rouleaux convectifs le long du canal dépendent des conditions de l’écoulement. Nous proposons une étude numérique pour ces conditions particulières d’écoulement pouvant conduire à des instabilités avec une évaluation de leur effet sur les transferts de chaleur. Les équations de Navier-Stokes et de l’énergie sont résolues numériquement par la méthode de volumes finis en prenant en compte les effets thermocapillaires. Les résultats présentés concernent l’influence des paramètres contrôlant l’écoulement (nombres de Reynolds, de Rayleigh, de Biot, de Marangoni et le rapport de forme) sur les motifs de l’écoulement et les échanges thermiques. Dans une seconde partie du travail, complémentaire à la première, une analyse de stabilité linéaire de l’écoulement dans un canal ouvert à surface libre d’extension latérale infinie est réalisée en utilisant la méthode spectrale de type collocation Chebyshev pour résoudre un système aux valeurs propres. Les diagrammes de stabilité déterminant les seuils des paramètres conduisant à l’instabilité thermoconvective ont été obtenus et analysés, ainsi que les structures spatiales associées. / Several studies both numerical and experimental have reported the presence of thermal instabilities in liquid films uniformly heated from below for specific boundary conditions and flows. The presence of these instabilities modifies the associated heat transfer. The subject of this PhD thesis is to study numerically the instability of three-dimensional laminar mixed convection within a liquid flowing on a horizontal channel heated uniformly from below. The upper surface is free and assumed to be flat. The variations of the surface tension with the temperature (Marangoni effect or thermocapillary effect) are taken into account. Although of great interest for many industrial applications, this problem has received little attention from an academic point of view. In this configuration, several types of thermoconvective structures may appear. When the strength of the buoyancy, thermocapillary effects and forced convective currents are comparable, the results show the development of instabilities in the form of steady longitudinal convective rolls similar to those encountered in the Poiseuille-Rayleigh-Bénard flow. To our knowledge, this is the first time that the Poiseuille-Rayleigh-Bénard flow associated to the Marangoni effects has been investigated. The number and spatial distribution of the convective rolls along the channel depend on the flow conditions. We propose a numerical study on the flow conditions that could lead to thermal instabilities with an evaluation of their effect on the heat transfer. The coupled Navier-Stokes and energy equations are solved numerically by the finite volume method taking into account the thermocapillary effects. The results presented concern the influence of several control parameters (the Reynolds, Rayleigh, Biot and Marangoni numbers and the aspect ratio of the channel) on the flow patterns and heat transfer characteristics. In the second part of this work, complimentary to the first, a linear stability analysis of a horizontal liquid film flowing in an open channel, with infinite lateral extension and uniform heating from below, is carried out. An eigenvalue problem is obtained in the course of this analysis which is solved numerically using the Chebyshev collocation spectral method. The stability diagrams determining the threshold parameters leading to thermoconvective instabilities were obtained and analyzed as well as the associated spatial patterns.