This thesis presents the author's contribution <p>to the design and implementation of COOLFluiD,<p>an object oriented software platform for <p>the high performance simulation of multi-physics phenomena on unstructured grids. In this context, the final goal has been to provide a reliable tool for handling high speed aerothermodynamic <p>applications. To this end, we introduce a number of design techniques that have been developed in order to provide the framework with flexibility<p>and reusability, allowing developers to easily integrate new functionalities such as arbitrary mesh-based data structures, numerical algorithms (space discretizations, time stepping schemes, linear system solvers, ),and physical models. <p>Furthermore, we describe the parallel algorithms <p>that we have implemented in order to efficiently <p>read/write generic computational meshes involving <p>millions of degrees of freedom and partition them <p>in a scalable way: benchmarks on HPC clusters with <p>up to 512 processors show their effective suitability for large scale computing. <p>Several systems of partial differential equations, <p>characterizing flows in conditions of thermal and <p>chemical equilibrium (with fixed and variable elemental fractions)and, particularly, nonequilibrium (multi-temperature models) <p>have been integrated in the framework. <p>In order to simulate such flows, we have developed <p>two state-of-the-art flow solvers: <p>1- a parallel implicit 2D/3D steady and unsteady cell-centered Finite Volume (FV) solver for arbitrary systems of PDE's on hybrid unstructured meshes; <p>2- a parallel implicit 2D/3D steady vertex-centered Residual Distribution (RD) solver for arbitrary systems of PDE's on meshes with simplex elements (triangles and tetrahedra). <p>The FV~code has been extended to handle all <p>the available physical models, in regimes ranging from incompressible to hypersonic. <p>As far as the RD code is concerned, the strictly conservative variant of the RD method, denominated CRD, has been applied for the first time in literature to solve high speed viscous flows in thermochemical nonequilibrium, yielding some preliminary outstanding results on a challenging double cone flow simulation. <p>All the developments have been validated on real-life testcases of current interest in the aerospace community. A quantitative comparison with experimental measurements and/or literature has been performed whenever possible. / Doctorat en Sciences de l'ingénieur / info:eu-repo/semantics/nonPublished

L'étude des phénomènes fondamentaux de<p>transport et de thermodynamique apparaissant lors de l'évaporation<p>en milieu poreux permet l'investigation d'applications pratiques<p>variées. Dans ce travail, nous développons des modèles fondamentaux<p>d'évaporation en milieu poreux que nous appliquons ensuite au<p>séchage en lit fluidisé de deux matériaux granulaires poreux :le<p>PVC et la levure.<p><p>Les modèles mis au point sont réalisés suivant une approche<p>multiéchelle. Nous nous intéressons tout d'abord aux phénomènes se<p>déroulant à l'échelle d'un pore. Les modèles développés à cette<p>échelle sont ensuite exploités dans le cadre d'une étude à l'échelle<p>d'un grain poreux. Le couplage des modèles de grain avec un modèle à<p>l'échelle du réacteur permet alors l'étude des applications<p>industrielles.<p><p>A l'échelle du pore, nous étudions les phénomènes de transport dans<p>un capillaire initialement rempli de liquide qui s'évapore vers<p>l'atmosphère ambiante. L'objectif est de prédire le taux<p>d'évaporation dans cette configuration. Nous nous intéressons<p>successivement à la modélisation du transport de matière par<p>convection-diffusion en phase gazeuse et la modélisation de l'impact<p>de films liquides adsorbés à la paroi du capillaire sur le transport<p>de matière. Ces deux modèles sont confrontés à des expériences<p>d'évaporation en capillaires cylindriques visualisées à l'aide de<p>deux dispositifs optiques. Le premier offre un suivi d'ensemble au<p>cours du temps du capillaire, alors que le second, un interféromètre<p>de Mach-Zehnder, permet une visualisation locale de la région<p>entourant le ménisque.<p><p>Le modèle portant sur le transport de matière par<p>convection-diffusion mène à la définition d'un critère non<p>dimensionnel permettant d'évaluer si la convection dans la phase<p>gazeuse dans le capillaire doit être prise en compte pour évaluer le<p>taux d'évaporation. Le modèle de film permet de prédire l'impact de<p>celui-ci sur l'évaporation en présence d'un gaz inerte lorsque les<p>mouvements convectifs en phase gazeuse sont négligeables. La<p>confrontation de ce modèle avec les profils d'épaisseur des films<p>obtenus à l'aide de interféromètre de Mach-Zehnder ne permet pas de<p>valider le modèle, et ce, suite à une trop grande incertitude sur<p>l'évaluation des interactions entre la paroi et le liquide.<p><p>A l'échelle d'un grain, nous développons un modèle discret par<p>réseau de pores et deux modèles continus pour tenter de prédire le<p>taux d'évaporation et la distribution des phases dans le milieu<p>poreux. Le modèle par réseau de pores prend en compte les transports<p>de matière par diffusion en phase gazeuse, par convection dans les<p>pores remplis de liquide et par convection dans les films liquides.<p>Les effets visqueux en phase liquide sont également modelisés. Pour<p>la prise en compte de ces derniers, nous montrons l'importance de<p>l'usage d'un algorithme approprié. Nous évaluons ensuite au travers<p>de trois nombres sans dimensions l'impact du transport par film et<p>des effets visqueux sur l'évaporation et la distribution des phases.<p>Cette analyse dimensionnelle est ensuite appliquée à l'étude de<p>réseaux de pores pour lesquels la section des liaisons les<p>constituant est idéalisée par des polygones réguliers. Pour les<p>modèles continus après une vérification simplifiée de<p>l'applicabilité de cette démarche, nous développons deux modèles<p>simples. Dans les deux modèles, l'étape de séchage à vitesse<p>constante est supposée limitée par le transport de matière externe<p>au grain. Le premier modèle, dit à front pénétrant, suppose que<p>l'étape de séchage à vitesse décroissante correspond à l'existence<p>d'un front d'évaporation qui s'enfonce dans la matrice poreuse. Le<p>second modèle, dit à surface d'échange variable, attribue cette même<p>étape du séchage à une diminution progressive de la surface<p>d'évaporation en surface du grain.<p><p>A l'échelle du réacteur, nous présentons deux modèles visant à<p>simuler deux types d'essais différents :le séchage en lit fluidisé<p>et la thermogravimétrie analytique. Ces deux modèles sont couplés<p>aux différents modèles à l'échelle de grain pour étudier le séchage<p>de PVC et de levure tant en lit fluidisé que lors des essais de<p>thermogravimétrie analytique.<p><p>Dans le cas du PVC, le modèle par réseau de pores ne peut pas être<p>appliqué de par la nécessité d'une trop grande puissance de calcul.<p>Au niveau des modèles continus, nous montrons que l'application du<p>modèle à surface d'échange variable permet de reproduire les courbes<p>de séchage expérimentales des essais en lit fluidisé.<p><p>Dans le cas de la levure, nous appliquons le modèle par réseau de<p>pores et le modèle à front pénétrant. L'utilisation du modèle par<p>réseau de pores nécessite une connaissance plus détaillée de la<p>structure poreuse des grains. Le traitement d'une analyse par<p>microtomographie nous permet d'obtenir un réseau de pores<p>expérimental. Celui-ci est utilisé pour montrer que la méthode de<p>caractérisation de la porosité par intrusion de mercure ne semble<p>pas adaptée à un matériau deformable comme la levure. Le même réseau<p>est utilisé pour simuler le séchage de grains de levure à l'aide du<p>modèle par réseau de pores. Les simulations sont réalisées sur des<p>réseaux équivalents à des coupes dans le solide. Le modèle par<p>réseau de pores et le modèle à front pénétrant permettent tous deux<p>de modéliser correctement le séchage de levure en lit fluidisé<p>moyennant l'ajustement de leurs paramètres ajustables,<p>respectivement la conductibilité des films liquide et la tortuosité.<p>Pour l'essai de thermogravimétrie, ils ne parviennent à approcher<p>que la première étape de celui-ci. Les avantages, défauts et<p>complémentarités de ces deux modèles sont discutés.<p><p>Nous évaluons ensuite l'impact du rétrécissement de la levure et des<p>types d'eau sur le séchage de ce matériau. Le rétrécissement est,<p>pour ce faire, mesuré à l'aide d'un stéréomicroscope. Ces premières<p>mesures, exploratoires, mènent à la définition d'un modèle empirique<p>de retrait du solide au cours de son séchage. En le prenant en<p>compte dans les modèles déjà appliqués à la levure, nous montrons<p>que le retrait a un impact significatif sur l'étape de séchage à<p>vitesse décroissante. Cet impact<p> peut cependant être masqué intégralement par la réévaluation<p>des paramètres ajustables des différents modèles. Finalement,<p>l'étude des types d'eau au travers d'un modèle simple appliqué à<p>l'essai de thermogravimétrie analytique montre que les types d'eau<p>ne doivent pas être pris en compte pour modéliser le séchage de<p>levure.<p><p>A l'issue de ce travail, nous disposons donc de modèles fondamentaux<p>d'évaporation en milieu poreux. Ceux-ci peuvent être appliqués à des<p>cas d'intérêt industriel, comme nous le réalisons pour le PVC et la<p>levure. Ils peuvent servir à améliorer la<p>compréhension de ces procédés. Ils représentent donc des<p>outils de choix pour la conception, le dimensionnement et<p>l'optimisation du séchage. / Doctorat en Sciences de l'ingénieur / info:eu-repo/semantics/nonPublished

Several types of instability can affect the interface between two fluids. For instance, a Rayleigh-Taylor instability (or density fingering) is encountered when a heavier fluid is placed upon a lighter one in the gravity field and double diffusive instabilities can be triggered by differential diffusivity of the different species present in the fluid. <p>In this context our work aims to understand theoretically in which way a chemical reaction can induce and influence such instabilities in a fluid initially at rest.<p>To understand the dynamics resulting from the coupling between chemical reactions and hydrodynamical instabilities we use chemical fronts as model systems. These fronts result from the coupling between autocatalytical chemical reactions and diffusion and they allow to create a self-organized interface between the products and the reactants. As during a chemical reaction the density may vary due to solutal and thermal effects, the products and the reactants can have different densities which may trigger convection movements leading to the destabilization of the fronts.<p><p>We have in particular studied the influence of the exothermicity of the reaction on the fingering of chemical fronts, focusing first on the influence of heat losses through the walls of the set-up.<p>These leaks have a marked influence on the dynamics because they affect the temperature profiles and hence the density profiles too. We have also classified the various types of instabilities that may appear dues to solutal and thermal effects. We have found a new type of hydrodynamic instability of statically stable fronts induced by the chemical reaction. <p><p>We have furthermore analyzed an isothermal model with two chemical species. If they diffuse at different rates the front can be subject to diffusive instabilities as well. We have shown that the coupling between such a diffusive instability and fingering can trigger complex dynamics. We have eventually studied the influence of an external electric field on the diffusive instabilities and on fingering underlying the possibility to destabilize otherwise stable fronts./<p>Différents types d'instabilités hydrodynamiques peuvent affecter les interfaces entre deux fluides comme par exemple, une instabilité de Rayleigh-Taylor (ou digitation de densité) quand un fluide plus dense se trouve placé au-dessus d'un fluide moins dense dans le champ de gravité ou des instabilités de double diffusion induites par des différences entre les diffusivités d'un soluté et de la chaleur contenus dans les fluides. Dans ce contexte, notre thèse s'attache à comprendre de manière théorique comment une réaction chimique peut influencer ces instabilités voire les générer dans un fluide initialement au repos. Pour étudier les dynamiques résultant du couplage entre réactions chimiques et instabilités hydrodynamiques, nous utilisons des systèmes modèles: les fronts chimiques de conversion résultant de la compétition entre réactions chimiques autocatalytiques et diffusion créant une interface auto-organisée entre les réactifs et les produits. Comme au cours d'une réaction chimique la densité peut varier par des effets solutaux et thermiques, les produits et les réactifs de densités différentes peuvent générer des mouvements de convection qui conduisent à la déstabilisation des fronts. <p><p>Nous avons en particulier étudié l'influence de l'exothermicité de la réaction sur les instabilités de digitation de fronts chimiques, en nous focalisant dans un premier temps sur l'influence des pertes de chaleur par les parois du réacteur.<p>Ces fuites ont un effet marqué sur les instabilitités car elles affectent les profils de température et donc les profils de densité dans le système. Nous avons également classifié les différentes instabilités qui peuvent apparaître via des changements de densité dûs à des effets thermiques et solutaux et mis en évidence un nouveau type de déstabilisation hydrodynamique de fronts statiquement stables induit par une réaction chimique. <p>Nous avons ensuite analysé un modèle isotherme impliquant deux espèces chimiques. Si ces dernières diffusent a des vitesses différentes le front peut être sujet à une instabilité diffusive. Nous avons montré qu'un couplage entre une telle instabilité diffusive et de la digitation peut être à l'origine de dynamiques complexes. Nous avons ensuite considéré l'influence d'un champ électrique sur les instabilité diffusives et de digitation en soulignant la possibilié de déstabiliser via ce champ des fronts initialement stables. / Doctorat en Sciences / info:eu-repo/semantics/nonPublished

A numerical code for solving the time-dependent incompressible 3D Navier-Stokes equations with finite volumes on overlapping staggered grids in cylindrical and rectangular geometry is developed. In the code, written in FORTRAN, the momentum equation for the velocity is solved by projection method and Poisson equation for the pressure is solved by ADI implicit method in two directions combined with discrete fast Fourier transform in the third direction. A special technique for overcoming the singularity on the cylinder's axis is developed. This code, taking into account dependence upon temperature of the viscosity, density and surface tension of the liquid, is used to study the fluid motion in a cylinder with free cylindrical surface (under normal and zero-gravity conditions); and in a rectangular closed cell with a source of thermocapillary convection (bubble inside attached to one of the cell's faces). They are significant problems in crystal growth and in general experiments in fluid dynamics respectively. Nevertheless, the main study is dedicated to the liquid bridge problem.<p><p>The development of thermocapillary convection inside a cylindrical liquid bridge is investigated by using a direct numerical simulation of the 3D, time-dependent problem for a wide range of Prandtl numbers, Pr = 0.01 - 108. For Pr > 0.08 (e.g. silicon oils), above the critical value of temperature difference between the supporting disks, two counter propagating hydrothermal waves bifurcate from the 2D steady state. The existence of standing and traveling waves is discussed. The dependence of viscosity upon temperature is taken into account. For Pr = 4, 0-g conditions, and for Pr = 18.8, 1-g case with unit aspect ratio an investigation of the onset of chaos was numerically carried out. <p><p>For a Pr = 108 liquid bridge under terrestrial conditions ,the appearance and the development of thermoconvective oscillatory flows were investigated for different ambient conditions around the free surface.<p><p>Transition from 2D thermoconvective steady flow to a 3D flow is considered for low-Prandtl fluids (Pr = 0.01) in a liquid bridge with a non-cylindrical free surface. For Pr < 0.08 (e.g. liquid metals), in supercritical region of parameters 3D but non-oscillatory convective flow is observed. The computer program developed for this simulation transforms the original non-rectangular physical domain into a rectangular computational domain.<p><p>A study of how presence of a bubble in experimental rectangular cell influences the convective flow when carrying out microgravity experiments. As a model, a real experiment called TRAMP is numerically simulated. The obtained results were very different from what was expected. First, because of residual gravity taking place on board any spacecraft; second, due to presence of a bubble having appeared on the experimental cell's wall. Real data obtained from experimental observations were taken for the calculations.<p> / Doctorat en sciences appliquées / info:eu-repo/semantics/nonPublished

Le présent manuscrit contribue à l’étude des écoulements liquides dans des milieux poreux artificiels, plus spécifiquement dans les cas où la matrice poreuse présente des gradients de perméabilité importants, par exemple dans un milieu stratifié ou fracturé. Nous étudions l’influence de tels milieux poreux hétérogènes sur différents types d’écoulements. Ce travail est principalement expérimental, mettant en oeuvre une technique optique non-intrusive appelée effet Christiansen. Cette méthode permet de déterminer quantitativement des distributions soit de température, soit de concentration au sein d’un milieu poreux. <p><p>Trois problèmes physiques sont étudiés: tout d’abord le problème de Horton-Rodgers-Lapwood qui est l’équivalent du très connu problème de Rayleigh-Bénard mais pour un milieu poreux, ensuite les phénomènes de dispersion hydrodynamique que l’on rencontre dans des écoulements multiphasiques. Cette dispersion hydrodynamique est essentiellement envisagée comme un processus macroscopique de diffusion, renforcé par rapport à la diffusion moléculaire que l’on rencontre en milieu fluide libre. Enfin, le troisième problème englobe les écoulements capillaires en milieux poreux en environnement de pesanteur réduite. Dans le cas d’écoulements immiscibles multiphasiques, il faut prendre en considération l’effet de la tension superficielle aux interfaces. Comme les effets capillaires sont partiellement masqués par les effets de pesanteur durant des expériences au sol, une étude précise des effets de mouillage dans ces écoulements en milieu poreux nécessite de les découpler au maximum des autres effets physiques. Un programme de recherche en microgravité a été réalisé, et un nouveau modèle mathématique qui prend en compte l’influence des forces capillaires a été élaboré dans le cadre d’une collaboration entre le Service de Chimie-Physique et le Prof. N.N. Smirnov du Département de Mécanique et de Mathématique de l’Université d’Etat de Moscou.<p><p><p>La structure de ce travail part du Chapitre 1, qui présente essentiellement les milieux poreux et leurs spécificités. Ce dernier introduit le formalisme et les concepts nécessaires au traitement des trois problèmes de recherche envisagés. Le Chapitre 2 présente ensuite une étude bibliographique du problème de Horton-Rodgers-Lapwood et des phénomènes de dispersion hydrodynamique en milieux poreux. Le Chapitre 3 est consacré à l’effet Christiansen. Le Chapitre 4 présente les dispositifs de laboratoire mis au point, ainsi qu’une compilation des résultats expérimentaux obtenus. Les problèmes d’écoulements capillaires sont exposés au Chapitre 5, étant donné que la technique expérimentale est différente de celle basée sur l’effet Christiansen. Ce Chapitre compare le nouveau modèle mathématique aux résultats des expériences menées en microgravité durant de nombreuses campagnes de vols paraboliques. Le Chapitre 6 referme ce travail par ses conclusions et perspectives. / Doctorat en sciences appliquées / info:eu-repo/semantics/nonPublished

The present research is an experimental investigation of the atomization of a superheated pressurized liquid jet that is exposed to the ambient pressure due to a sudden depressurization. This phenomena is called flashing and occurs in several industrial environments.<p><p>Liquid flashing phenomena holds an interest in many areas of science and engineering. Typical examples one can mention: a) the accidental release of flammable and toxic pressure-liquefied gases in chemical and nuclear industry; the failure of a vessel or pipe in the form of a small hole results in the formation of a two-phase jet containing a mixture of liquid droplets and vapor, b) atomisation improvement in the fuel injector technology, c) flashing mechanism occurrence in expansion devices of refrigerator cycles etc. The interest in flashing events is especially true in the safety field where any unexpected event is undesirable. In case of an accident, flammable or toxic gas clouds are anticipated in close regions of the release because of the sudden phase change .Due to the non-equilibrium nature of the flow in these near field regions, conducting accurate data measurements for droplet size and velocity is a challenging task resulting in scarce data in the very close area.<p><p>This research has been carried out at the von Karman Institute (VKI) within the 5th framework of European Commission to fulfill the goal of understanding of source processes in flashing liquids in accidental releases. The program is carried out under name of FLIE (Flashing Liquids in Industrial Environments)(Contract no: EVG1-CT-2000-00025). The specific issues that are presented in this thesis study are the following:a) a comprehensive state of art of the jet break up patterns, spray characteristics and studies related to flashing phenomena; b)flashing jet breakup patterns and accurate characterization of the atomized jet such as droplet diameter size, velocity and temperature evolution through carefully designed laboratory-scale experiments; c) the influence of the initial storage conditions on the final atomized jet; d) a physical model on the droplet transformation and rapid evaporation in aerosol jets.<p><p>In order to characterize the atomization of the superheated liquid jet, laser-based optical techniques like Particle Image Velocimetry (PIV), Phase Doppler Anemometry (PDA) are used to obtain information for particle diameter and velocity evolution at various axial and radial distances. Moreover, a high-speed video photography presents the possibility to understand the break-up pattern changes of the simulating liquid namely R-134A jet in function of driving pressure, superheat and discharge nozzle characteristics. Global temperature measurements with an intrusive technique such as thermocouples, non-intrusive measurements with Infrared Thermography are performed. Cases for different initial pressures, temperatures, orifice diameters and length-to-diameter ratios are studied. The break-up patterns, the evolution of the mean droplet size, velocity, RMS, turbulence<p>intensity and temperature along the radial and axial directions are presented in function of initial parameters. Highly populated drop size and velocity count distributions are provided. Among the initial storage conditions, superheat effect is found to be very important in providing small droplets. A 1-D analytical rapid evaporation model is developed in order to explain the strong temperature decrease during the measurements. A sensitivity analysis of this model is provided.<p> / Doctorat en sciences appliquées / info:eu-repo/semantics/nonPublished

This thesis deals with the dynamics of a thin liquid film falling down a heated plate. The heating yields surface tension gradients that induce thermocapillary stresses on the free surface, thus affecting the stability and the evolution of the film. Accounting for the coherence of the flow due to viscosity, two main approaches that reduce the dimensionality of the original problem are usually considered depending on the flow rate (as measured by the Reynolds number): the `long wave' asymptotic expansion for small Reynolds numbers and the `integral boundary layer' approximation for moderate Reynolds numbers. The former suffers from singularities and the latter from incorrectness of the instability threshold for the occurrence of hydrodynamic waves. Thus, the aim of this thesis is twofold: in a first part, we define quantitatively the validity of the `long wave' evolution equation (Benney equation) for the film thickness h including the thermocapillary effect; and in a second part, we improve the `integral boundary layer' approach by combining a gradient expansion to a weighted residual method. <p>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.<p>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.<p>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.<p>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.<p><p>------------------------------------------------<p><p>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.<p>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. <p>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<p> / Doctorat en sciences appliquées / info:eu-repo/semantics/nonPublished

<p align="justify">This work deals with the development of numerical methods for compressible flow simulation with application to the interaction of fluid flows and structural bodies.</p><p><p><p align="justify">First, we develop numerical methods based on multidimensional upwind residual distribution (RD) schemes. Theoretical results for the stability and accuracy of the methods are given. Then, the RD schemes for unsteady problems are extended for computations on moving meshes. As a second approach, cell centered and vertex centered finite volume (FV) schemes are considered. The RD schemes are compared to FV schemes by means of the 1D modified equation and by the comparison of the numerical results for scalar problems and system of Euler equations. We present a number of two and three dimensional steady and unsteady test cases, illustrating properties of the numerical methods. The results are compared with the theoretical solution and experimental data.</p><p><p><p align="justify">In the second part, a numerical method for fluid-structure interaction problems is developed. The problem is divided into three distinct sub-problems: Computational Fluid Dynamics, Computational Solid Mechanics and the problem of fluid mesh movement. The problem of Computational Solid Mechanics is formulated as a system of partial differential equations for an anisotropic elastic continuum and solved by the finite element method. The mesh movement is determined using the pseudo-elastic continuum approach and solved again by the finite element method. The coupling of the problems is achieved by a simple sub-iterative approach. Capabilities of the methods are demonstrated on computations of 2D supersonic panel flutter and 3D transonic flutter of the AGARD 445.6 wing. In the first case, the results are compared with the theoretical solution and the numerical computations given in the references. In the second case the comparison with experimental data is presented.</p> / Doctorat en Sciences de l'ingénieur / info:eu-repo/semantics/nonPublished