This dissertation considers the genesis and dynamics of interfacial instability in vertical laminar gas-liquid flows, using as a model the two-dimensional channel flow of a thin falling film sheared by counter-current gas. The methodology is linear stability theory by means of Orr-Sommerfeld analysis together with direct numerical simulation of the two-phase flow in the case of nonlinear disturbances. The influence of two main flow parameters on the interfacial dynamics, namely the film thickness and pressure drop applied to drive the gas stream, is investigated. To make contact with existing studies in the literature, the effect of various density and viscosity contrasts as well as surface tension is also examined. Energy budget analyses based on the Orr-Sommerfeld theory reveal various coexisting unstable modes (interfacial, shear, internal) in the case of high density contrasts, which results in mode coalescence and mode competition, but only one dynamically relevant unstable interfacial mode for low and intermediate density contrast. Furthermore, high viscosity contrast and increases in surface tension lead to some amount of mode competition for thin film. A study of absolute and convective instability for low density contrast shows that the system is absolutely unstable for all but two narrow regions of the investigated parameter space. These regions are extended at intermediate density contrast and exhibit only small changes with increased viscosity contrast or surface tension. Direct numerical simulations of the system with low density contrast show that linear theory holds up remarkably well upon the onset of large-amplitude waves as well as the existence of weakly nonlinear waves. For high density contrasts corresponding more closely to an air-water-type system, linear stability theory is also successful at determining the most-dominant features in the interfacial wave dynamics at early-to-intermediate times. Nevertheless, the short waves selected by the linear theory undergo secondary instability and the wave train is no longer regular but rather exhibits chaotic motion. Furthermore, linear stability theory also predicts when the direction of travel of the waves changes - from downwards to upwards. The practical implications of this change in terms of loading and flooding is discussed. The change in direction of the wave propagation is represented graphically for each investigated system in terms of a flow map based on the liquid and gas flow rates and the prediction carries over to the nonlinear regime with only a small deviation. Besides the semi-analytical and numerical analyses, experiments with an practically relevant setup and flow system have been carried out to benchmark and validate the models developed in this work.

Les interfaces entre deux fluides sont le siège de nombreuses instabilités de forme de l’interface si un champ électrique intense est appliqué : génération de gouttelettes, jets micrométriques, etc. Le contrôle de telles instabilités est indispensable pour une fabrication optimale de microsphères ou microfibres : taille, propriétés physico-chimiques, dispersion et structuration spatiale macroscopique d’un agrégat de tels objets. Cette diversité provient de la compétition entre la tension de surface et la gravité avec l’électrodynamique des fluides sous champ électrique induite par les charges électriques, les charges de polarisation, les décharges électriques et/ou le vent ionique. La thèse expérimentale s’articule autour de deux thèmes. Le premier, une compréhension intime des phénomènes spatio-temporels observables lorsqu’un injecteur métallique à la haute tension placé au-dessus d’un bain. Une instabilité originale menant à la formation d’une cloche fluide macroscopique connectant les deux électrodes a été mise en évidence et caractérisée non linéairement. La bifurcation est sous-critique et imparfaite. Le second thème propose une méthode originale de fabrication de microfibres modifiées en une étape par électrofilage au mouillé. Le polymère électrofilé choisi est le PSMA et celui permettant la modification, le PEGDA. Cette étude a été réalisée dans un contexte d’applications de type catalyse. Pour cela les fibres ont été fonctionnalisées à l’aide de la peroxydase (HRP) comme protéine modèle. Les résultats montrent notamment une meilleure stabilité temporelle avec la possibilité de réutilisation du matériau en comparaison à la catalyse utilisant des méthodes standards. / Interfaces between two fluids can lead to various interfacial shape instabilities if an electrical field is applied. Leading, for instance, to micrometric droplets or jets formation. Controlling those instabilities is much-needed for an optimal fabrication of microspheres or microfibers : size, physicochemical properties, dispersion and macroscopic spatial structuring of aggregates of those kind of objects. This diversity is based in the competition between surface tension and gravity forces with gravity during the electrodynamics of fluids under electric field induced by electrical charges, polarization charges, electrical discharges and ionic wind. The experimental thesis deals with two main topics. The first one is a precise understanding of spatiotemporals phenomena occurring in a configuration made of a metallic injector raised to high voltage placed above a liquid bath. We present the formation of an original instability leading to a macroscopic bell-shaped link between both electrodes and its non linear characterization. The bifurcation is subcritical and imperfect. The second topic, based on the experience gained with the first one, is an original method of fabrication of microfibers modified in only one step by wet electrospinning. The chosen electrospun polymer is PSMA and the one used for modification is PEGDA. This study has been realized with a catalyze application context. To do so, fibers has been functionalized with peroxydase (HRP) as the model protein. The results especially show a better temporal stability and possible reuse compared to catalysis with standard methods.

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.