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 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.