Ce travail de thèse concerne l'étude des écoulements de fluides gazeux compressibles laminaires subsonique-supersonique dans une tuyère de type convergent-divergent. Les écoulements étudiés sont à nombres de Reynolds modérés et s'affranchissent de l'hypothèse de condition adiabatique de paroi couramment utilisée afin de mieux prendre en compte les phénomènes de transfert de chaleur par convection et rayonnement avec le milieu extérieur. Cette étude des phénomènes de transferts conjugués a permis de déterminer le comportement dynamique simultané du fluide et de la paroi de la tuyère. Enfin, compte tenu des niveaux élevés de températures mis en jeu dans ces systèmes, une analyse concernant le comportement thermomécanique de l'ensemble de la structure de paroi avec des matériaux monocouches et multicouches a été réalisé. De nombreuses configurations géométriques, propriétés physiques et conditions aux limites sur le fluide et la paroi ont été analysées. Les résultats présentés montrent, la structure des écoulements à travers les iso-contours de vitesses, des nombres de Mach, des pressions dans le fluide, des températures dans le fluide et dans la paroi ainsi que les déformations et les contraintes de la paroi qui résultent des couplages thermomécaniques. Une analyse des performances de la tuyère, en termes de force de poussée et de coefficient de débit spécifique, est largement discutée dans ce travail.

Cette thèse est consacrée è la modélisation et a la simulation de l'ensemble PISTON /SEGMENT/ CHEMISE, d'un moteur thermique. Dans un moteur, le piston muni de segments est en mouvement relatif dans la chemise. Ce système est lubrifié par projection d'huile, les segments et la surface de la chemise "contrôlent" la quantité de lubrifiant restant sur la paroi. Pour modéliser ce phénomène, le modèle conservatif parabolique-hyperbolique d'Elrod-Adams (P − θ), prenant en compte la cavitation (présence de gaz dans le fluide, aspect multiphasique de l'écoulement), est couplé a un modèle de contact, celui de Greenwood-Tripp (basé sur une approche statistique de la surface). En se restreignant a l'écoulement et au contact, nous négligeons beaucoup de phénomènes physiques pouvant intervenir comme le grippage ou les effets thermiques. Dans une première partie, les phénomènes mécaniques sont étudiés, de la réalisation de la surface des chemises è la modélisation de la physique de la lubrification du contact segment/chemise. Cette étude expose les simplifications effectuées et le choix des conditions aux limites pour la prise en compte de la cavitation en tant que problème a frontière libre, notamment au niveau du débit d'entrée pour simuler un fonctionnement "normal". Dans une deuxième partie, la mise en oeuvre et la modification de l'algorithme "classique" pour des surfaces mesurées sont décrites en détail. Nous présentons les performances simulées de différentes surfaces, et nous les comparons è des essais effectués par Renault. Dans une troisième partie, des modifications sont introduites au modèle (P − θ) utilisé précédemment afin de rendre l'écoulement cohérent avec le modèle de Navier-Stokes. Bien que ce nouveau problème puisse être posé en dimension 2, l'étude mathématique ne concerne que la dimension 1. On se ramène a un système dynamique dont les inconnues sont les frontières libres. L'unicité de la solution impose l'introduction d'un paramètre artificiel. Le comportement local et/ou global en temps est déterminé par les paramètres géométriques et les conditions aux limites. L'introduction de deux paramètres supplémentaires permet de démontrer l'existence et l'unicité de la solution. Dans une dernière partie, nous comparons le modèle de cavitation proposé avec le modèle P − θ.

Les ondes gravito-inertielles jouent un rôle majeur dans les échanges d'énergie globaux sur la planète. Si la génération des ondes est bien connue dans l'atmosphère et l'océan, le devenir de ces ondes au cours de leur propagation n'est pas complètement défini aujourd'hui. Ces ondes peuvent interagir de façon non-linéaire avec elles-mêmes et créer des structures de plus petite échelle qui vont se dissiper plus facilement. Ainsi, le phénomène d'instabilité paramétrique sous-harmonique (PSI), a été étudié de façon expérimentale. Nous avons effectué la première mise en évidence expérimentale de l'interaction de trois ondes planes inertielles bi-dimensionnelles, sous la forme d'une triade résonnante. Cette étude améliore en outre la compréhension de la turbulence en rotation. Les ondes internes peuvent aussi créer, ou interagir avec des écoulements lents de grande échellequi peuvent modifier la biodiversité au fond des océans. Nous avons mis en évidence une situation expérimentale à l'origine d'un tel écoulement moyen induit par les ondes et, à l'aide d'un modèle théorique simplifié, nous avons expliqué la formation de ces écoulements. Enfin, on étudie également des tourbillons en fluide stratifié pour permettre de futures études sur l'interaction d'ondes gravito-inertielles avec des tourbillons.

Lors de sa détente dans une turbine, la vapeur subit une chute d'enthalpie qui entraîne sa condensation spontanée sous forme d'un nuage de gouttelettes submicroniques. Ces gouttes vont se déposer sur les aubes aval et y former un film d'eau. Ce dernier est arraché sous l'effet de l'écoulement de vapeur environnant ce qui crée des gouttes de quelques dizaines de microns qui peuvent se redéposer sur les aubes aval. Ces phénomènes sont à l'origine de pertes, généralement regroupées sous le terme générique de "pertes par humidité", estimées grâce à la loi de Baumann. Le but de cette thèse est de mettre en place dans un code méridien des modèles permettant la prise en compte des phénomènes de condensation, de déposition et d'écoulement des films liquides afin de pouvoir estimer les pertes liées à chacun d'entre eux. Dans cette optique un modèle diphasique homogène permettant d'avoir accès à la fraction massique de liquide et au nombre de gouttes est implanté dans le code méridien. Ce modèle est validé sur un cas test expérimental de détente en tuyère et alimente le modèle de déposition. Les contributions de la diffusion, de la turbophorèse, de la thermophorèse, de la gravité et de l'inertie des gouttes à la déposition sont étudiées. Il apparaît nécessaire de prendre en compte la diffusion, la turbophorèse et la déposition inertielle sur les bords d'attaque des aubes. Un modèle permettant d'avoir accès à l'épaisseur et à la vitesse d'un film liquide soumis au cisaillement d'un écoulement environnant, à la force de frottement sur la paroi et aux effets de la rotation est mis en place. Ce modèle est validé par rapport à des résultats expérimentaux d'écoulement de film liquide sur une plaque plane dans des conditions proches de celles rencontrées en turbine à vapeur basse pression. Finalement, un calcul réalisé sur une géométrie réelle de turbine basse pression de 8 étages permet de démontrer l'applicabilité de la méthodologie mise en place sur un cas industriel. Les contributions des phénomènes étudiés aux pertes par humidité sont explicitées.

L'écoulement dans les turbomachines est tri-dimensionnel et instationnaire. Actuellement, les concepteurs de moteurs cherchent à réduire l'encombrement et le poids des machines. En conséquence, les interactions entre les roues, appelées interactions rotor-stator, sont renforcées. Parmi elles, l'effet potentiel remontant n'est désormais plus négligeable malgré sa rapide atténuation spatiale. Dans cette étude, cet effet potentiel remontant a été analysé sur une configuration spécialement conçue : une grille linéaire d'aubes de turbine, suivie de barreaux défilants en aval à une distance de 20% de corde axiale, simulant des aubes de rotor en aval. La grande échelle du banc d'essais facilite l'étude du comportement de la couche limite des aubes de la grille. Des mesures de pression et d'anémométrie laser à deux composantes, synchronisées avec le défilement des barreaux aval sont réalisées. Le nombre de Reynolds, basé sur la corde, est 1.6 × 105. Une grille de turbulence placée en amont de la grille d'Aube afin de pouvoir augmenter le taux de turbulence amont a été utilisée. Des résultats de mesures en absence de cette grille (faible taux de turbulence amont) sont également présentés et analysés. Une modélisation numérique, basée sur un calcul laminaire avec un préconditionnement basse vitesse pour la même configuration, a montré la déformation des lignes de courant de l'écoulement dans le canal inter-aubes, en fonction de la position du barreau aval. La distribution de pression autour de l'aube est également périodiquement modifiée. Les résultats stationnaires expérimentaux, en absence de tout barreau aval, ont révélé un décollement de la couche limite à l'extrados de l'aube à bas taux de turbulence amont(Tu−am = 1.2%) qui est supprimé à haut taux de turbulence amont (Tu−am = 4.2%) ; la couche limite commence alors sa transition par un mode by-pass. Aucun effet instationnaire dans la couche limite n'a été observé à l'intrados, quel que soit le taux de turbulence amont. L'étude instationnaire, avec le défilement des barreaux en aval, a permis de mettre en évidence un décollement périodique de la couche limite à l'extrados à bas taux de turbulence amont (Tu−am = 1.8%). Dans ce cas, la couche limite suit deux modes de transition au cours d'une période : une transition par décollement et une transition bypass. Au contraire, dans le cas à fort taux de turbulence amont (Tu−am = 4.0%), aucun décollement de la couche limite n'a été décelé. La couche limite est sujette à l'effet instationnaire à l'extrados. Elle est devenue turbulente au bord de fuite à tout instant par un mode by-pass. Cette étude a montré que l'effet potentiel issu d'un roue en aval est du même ordre de grandeur que les effets de sillage et doit être pris en compte dans l'analyse des phénomènes. Par des méthodes d'indexation de roues, le décollement de la couche limite pourrait être supprimé.

Cette thèse porte sur la caractérisation numérique et la modélisation stochastique de l'accélération du fluide pour l'écoulement en canal à grand nombre de Reynolds. La motivation concerne l'observation et l'analyse des effets de l'intermittence liés aux interactions à longue portée à travers le canal. Dans la première partie, l'accélération est étudiée par simulation numérique directe pour trois différents nombres de Reynolds (180, 590 et 1000). La lognormalité de la norme de l'accélération est observée quelle que soit la distance à la paroi. Un profil universel de la norme de l'accélération est également recherché par analyse dimensionnelle. La seconde partie présente une modélisation stochastique de l'accélération basée sur la décomposition norme/orientation. Le modèle stochastique pour la norme s'appuie sur un processus de fragmentation afin de représenter les interactions à longue portée à travers le canal. Pour l'orientation, l'évolution vers l'isotropie lorsque la distance à la paroi augmente (observée par la DNS) est reproduite grâce à un modèle de marche aléatoire sur une sphère. Ces modèles ont été appliqués à l'approche LES-SSAM (Stochastic Subgrid Acceleration Model) introduite par Sabel'nikov, Chtab et Gorokhovski. Nos calculs montrent que les estimations de la vitesse moyenne, du spectre d'énergie, des contraintes de l'écoulement et de la non-gaussianité des statistiques de l'accélération peuvent être améliorées de façon significative par rapport à la LES classique. L'intérêt de l'approche LES-SSAM, donnant un accès vers la structure intermittente de sous-maille, est illustré dans la dernière partie, par l'étude du transport de particules inertielles ponctuelles par l'écoulement de canal. Cette étude commence par l'analyse par DNS de l'influence des structures de paroi sur la dynamique des particules

L'hétérogénéité de la zone non saturée joue un rôle important dans le transfert d'eau et de soluté car elle accentue à la fois le développement des zones de stockage temporelles et les écoulements préférentiels. Par conséquent, la validation des modèles prédictifs nécessite le développement des outils expérimentaux spécifiques afin d'observer et de quantifier les mécanismes de transport impliqués dans un système non saturé hétérogène. Cette thèse vise à étudier l'effet combiné de la vitesse d'infiltration, de la barrière capillaire et l'angle de la pente d'interface entre deux matériaux sur les processus de l'écoulement de l'eau et du transport de soluté dans un modèle physique, le lysimètre de laboratoire 1x1x1.6 m3, nommé LUGH (Lysimeter for Urban Groundwater Hydrology) et un modèle numérique 3D de ce lysimètre. Le lysimètre LUGH est rempli par un sable fin et un mélange bimodal (50 % sable fin et 50 % gravier) en deux configurations: un profil uniforme de matériau bimodal ou un profil avec deux couches avec une pente de 14o. Ces agencements figurent l'hétérogénéité structurale et texturale observée sur un des sites expérimentaux de l'OTHU (Observatoire de Terrain en Hydrologie Urbaine) : le bassin d'infiltration d'eaux pluviales Django Reinhardt géré par la ville de Lyon (France). Le lysimètre est alimenté en eau et avec un traceur inerte (bromure de potassium, KBr) sur une partie de la surface par un système d'arrosage automatique. Les effluents ont été recueillis dans quinze sorties différentes en bas du lysimètre. La forte hétérogénéité des flux des sorties et des courbes de percée souligne la mise en place des écoulements préférentiels résultant à la fois de l'effet de barrière capillaire et de l'effet de fond du lysimètre. A partir des résultats expérimentaux, la modélisation numérique à l'aide de logiciel COMSOL MultiphysicsTM a permis de mieux comprendre les mécanismes responsables de ces transferts hétérogènes. Lorsque le modèle numérique validé, un test de sensibilité a été conduit pour étudier les effets de la vitesse d'infiltration et de la pente de l'interface sur l'écoulement. Les résultats montrent que la diminution de la vitesse d'infiltration ou l'augmentation de la pente de l'interface favorisent le développement des écoulements préférentiels. Notre étude a donné également des renseignements pertinents sur le couplage entre les processus hydrodynamiques et le transfert des solutés dans les sols non saturés hétérogènes en soulignant le rôle de la géométrie des interfaces ainsi que celui des conditions aux limites comme des facteurs clés pour la quantification des écoulements préférentiels.

Ce travail est une contribution expérimentale à l'étude rhéologique et en écoulement de fluides complexes (à seuil et thixotropes) transitant dans un circuit hydraulique. Il s'agit notamment de suspensions de bentonite ainsi que des complexes bentonite/polymère. L'étude porte sur l'évolution des pertes de charge et des champs de vitesse et se situe en régime laminaire, transitoire et turbulent. L'étude a été réalisée en utilisant un vélocimètre ultrasonore Doppler pulsé développé au laboratoire. Dans la première partie expérimentale de la thèse, des mesures rhéologiques et en écoulement ont été effectuées sur des suspensions de bentonite pures (sans additifs) à différentes concentrations. A partir des essais sur boucle hydraulique, une étude détaillée est présentée sur l'évolution des coefficients de frottement et des profils de vitesse pour les différents régimes d'écoulement. Dans une seconde partie, une suspension de bentonite pure et des mélanges bentonite/CMC à différentes concentrations massiques ont été étudiées en termes de comportement rhéologique et hydrodynamique en écoulement en conduite. En ce qui concerne les mesures effectuées en boucle hydraulique, il a été montré que le polymère a des propriétés viscosifiantes en régime laminaire. Par contre, en régime turbulent, le polymère agit comme un réducteur de frottement. La dernière partie de la thèse a été consacrée à l'étude de l'influence de la température sur le comportement rhéologique des solutions de polymère et des mélanges argile/polymère. Les mesures rhéologiques à différents paliers de températures ainsi que les balayages en température ont mis en évidence le caractère thermodépendant des dispersions.

La vectorisation de la poussée d'une tuyère propulsive supersonique axisymétrique est étudiée par le biais d'une injection fluidique transversale dans sa partie divergente. Cette étude menée dans le cadre du programme PERSEUS du CNES a été motivée par la recherche d'une solution alternative au pilotage conventionnel de la poussée par actionneurs mécaniques. Le travail de la thèse, tout en s'appuyant sur des approches expérimentale et numérique, comprend essentiellement une large étude paramétrique concernant principalement la position de l'injection, la forme de la tuyère, la nature et le débit du fluide injecté. L'analyse des résultats montre que pour certaines configurations optimales, des angles de déviation pertinents peuvent être obtenus pour des taux d'injections modérés. L'analyse numérique étendue aux écoulements chauds multi-espèces, plus proches des applications réelles, a montré que la vectorisation fluidique reste très performante lors de l'injection de produits de combustion dans le divergent.

Cette thèse est une contribution à l'étude de l'écoulement autour d'un drone UCAV générique. L'écoulement tourbillonnaire se développant à l'extrados d'aile volante à flèche modérée peut générer des phénomènes aérodynamiques non linéaires et ainsi dégrader fortement les qualités de vol, les performances aérodynamiques et limiter le domaine de vol. Plusieurs études expérimentales en soufflerie ont été réalisées sur une maquette de drone de flèche 53° afin de décrire les phénomènes de pitch down/up et d'instabilité latérale. Une analyse détaillée des mesures PIV a permis de mettre en évidence une topologie tourbillonnaire complexe à l'origine du comportement non linéaire du drone. L'utilisation d'un critère topologique a permis la détection de sous structures tourbillonnaires et instationnaires émanant du bord d'attaque et convectées dans la couche de mélange. La compréhension de la phénoménologie de l'écoulement tourbillonnaire a été complétée par différentes simulations numériques. Ces simulations numériques ont permis de prédire globalement le comportement de l'aéronef. La deuxième partie de cette thèse concerne le contrôle de l'écoulement aérodynamique à l'extrados du drone afin d'améliorer les performances aérodynamiques. Une stratégie de contrôle par des dispositifs passifs a été évaluée expérimentalement et numériquement. La mise en place de dispositifs de contrôle aux bords d'attaque a permis une modification significative de la topologie de l'écoulement tourbillonnaire à l'extrados du drone. Ces dispositifs ont ainsi contribué à une réduction de l'influence des phénomènes aérodynamiques non linéaires sur le comportement du drone.