Le but du travail présenté est la compréhension puis la simulation numérique des mécanismes d'érosion des augets de turbine Pelton par impacts répétés de gouttes d'eau dans le but de prédire la durée de vie des composants. Tout d'abord, les phénomènes de propagation d'ondes dans les milieux fluide et solide sont étudiés. Cela permet de mettre en lumière l'évolution temporelle et la distribution spatiale des pressions de contact, et l'apparition de microjets par éjection supersonique du fluide au contact. Les études expérimentales de l'érosion par gouttes d'eau traduisent un dommage basé sur la fissuration par fatigue. Des simulations numériques en dynamique rapide couplées fluide-structure sont alors effectuées. Le domaine solide est discrétisé par la Méthode des éléments Finis (MEF), et le domaine fluide par la méthode Smoothed Particle Hydrodynamics (SPH), qui est une méthode particulaire (sans maillage) particulièrement adaptée aux grandes distorsions et au suivi des surfaces libres. L'analyse des états de contraintes vient corroborer la nature cyclique de l'endommagement. La simulation d'érosion est alors réalisée à l'aide de critères de fatigue multiaxiaux. Le choix se porte vers un premier critère général de l'American Society of Mechanical Engineers (ASME), utilisant les valeurs principales des différences de contraintes au cours du temps. Le second choix concerne un critère à plan critique : le critère de Dang Van 2. Il traite séparément la contrainte hydrostatique et le cisaillement alterné maximal local. Ces critères permettent de définir les régions érodées du solide au bout d'un nombre d'impact donné, ce qui fait de cette démarche une méthode prédictive. Une étude paramétrique pour différentes tailles de gouttes et vitesses d'impact est ensuite réalisée, puis on évalue l'influence de la présence d'une couche de coating. / The goal of this work is the comprehension and the numerical simulation of erosion caused by repeated droplet impact on Pelton turbine buckets, to predict the lifetime of these components. First, waves propagation phenomenon inside fluid and solid domains are presented, which allows determining the time evolution and spatial distribution of contact pressure, and the birth of lateral microjets by supersonic ejection of the fluid on the contact. Experimental studies of erosion by droplet impact highlight a fatigue cracking-based erosion mechanism. Then, coupled FSI computation are performed. The solid subdomain is discretized by the Finite Element Method (FEM), and the fluid subdomain by the Smoothed Particle Hydrodynamics (SPH), which is a particle method (meshless) effectively recommended for large distortions and free surface tracking. Stress analysis confirms the cyclic nature of the damage mechanism, and erosion simulation is performed using multiaxial fatigue criteria. The first selected criterion is a general one from the American Society of Mechanical Engineers (ASME) using principal values of stress differences over time. The second one is the Dang van 2 criterion, belonging to the family of critical plane criteria. This criterion considers separately the effects due to hydrostatic stress on one hand, and the ones induced by maximum local shear on the other. These two criteria are used to building the equivalent eroded zones of the solid subdomain for a given number of impacts, which allows to qualify this procedure as a predictive predictive. Finally, a parametric study for different droplet sizes and velocites is computed, and the effects of a coating layer are investigated.

Dans cette thèse, nous étudions des modèles mathématiques concernant certains problèmes d'écoulement de fluide à densité variable. Le premier chapitre résume l'ensemble de la thèse et se concentre sur les résultats obtenus, la nouveauté et la comparaison avec la littérature existante. Dans le deuxième chapitre, nous étudions la stabilisation locale des équations non homogènes de Navier-Stokes dans un canal 2d autour du flot de Poiseuille. Nous concevons un contrôle feedback de la vitesse qui agit sur l'entrée du domaine de sorte que la vitesse et la densité du fluide soient stabilisées autour du flot de Poiseuille, à condition que la densité initiale soit donnée par une constante additionnée d'une perturbation dont le support se situe loin du bord latéral du canal. Dans le troisième chapitre, nous étudions un système couplant les équations de Navier-Stokes compressibles à une structure élastique située à la frontière du domaine fluide. Nous prouvons l'existence locale de solutions solides pour ce système couplé. Dans le quatrième chapitre, notre objectif est d'étudier la nulle- contrôlabilité d'un problemè d'interaction fluide-structure linéarisé dans un canal bi dimensional. L'écoulement du fluide est ici modélisé par les équations de Navier-Stokes compressibles. En ce qui concerne la structure, nous considérons une poutre de type Euler-Bernoulli amortie située sur une partie du bord. Dans ce chapitre, nous établissons une inégalité d'observabilité pour le problème considéré d'interaction fluid-structure linéarisé qui constitue le premier pas vers la preuve de la nulle contrôlabilité du système. / In this thesis we study mathematical models concerning some fluid flow problems with variable density. The first chapter is a summary of the entire thesis and focuses on the results obtained, novelty and comparison with the existing literature. In the second chapter we study the local stabilization of the non-homogeneous Navier-Stokes equations in a 2d channel around Poiseuille flow. We design a feedback control of the velocity which acts on the inflow boundary of the domain such that both the fluid velocity and density are stabilized around Poiseuille flow provided the initial density is given by a constant added with a perturbation, such that the perturbation is supported away from the lateral boundary of the channel. In the third chapter we prove the local in time existence of strong solutions for a system coupling the compressible Navier-Stokes equations with an elastic structure located at the boundary of the fluid domain. In the fourth chapter our objective is to study the null controllability of a linearized compressible fluid structure interaction problem in a 2d channel where the structure is elastic and located at the fluid boundary. In this chapter we establish an observability inequality for the linearized fluid structure interaction problem under consideration which is the first step towards the direction of proving the null controllability of the system.

Dans cette thèse on a développé des algorithmes semi-implicites procédures partagées et monolithiques pour l'interaction entre un fluide gouverné par le modèle de Navier Stokes et une structure. Dans le premier chapitre, on présente un algorithme semi-implicite procédures partagées pour l'interaction entre un fluide et une structure gouvernée soit par les équations d'élasticité linéaire ou soit par le modèle de Saint-Venant Kirchhoff non linéaire. Dans le second chapitre, on propose un algorithme semi-implicite procédures partagées pour l'interaction entre un fluide et une structure de modèle linéaire et on montre un résultat de stabilité inconditionnelle en temps de l'algorithme. Un problème d'optimisation est résolu dans les deux algorithmes précédents, afin de satisfaire les conditions de continuité des vitesses et d'égalité des contraintes à l'interface. Durant les itérations de BFGS pour résoudre le problème d'optimisation, le maillage fluide reste fixe et la matrice fluide n'est factorisée qu'une seule fois, ce qui réduit l'effort de calcul. Dans le troisième chapitre, un algorithme semi-implicite monolithique pour l'interaction entre un fluide et une structure de modèle linéaire est proposé. L'algorithme utilise un maillage global pour le domaine fluide structure. La condition de continuité des vitesses à l'interface est automatiquement satisfaite et celle de l'égalité des contraintes n'apparaît pas explicitement dans la formulation faible. A chaque pas de temps on résout un système monolithique d'inconnues vitesse et pression définies sur le domaine global. Le temps CPU est réduit quand l'approche monolithique est utilisée à la place des procédures partagées. / Our aim was to develop some partitioned procedures and monolithic semi-implicit algorithms for solving the interaction between a fluid governed by Navier Stokes equations and a structure. In the first chapter, we propose a partitioned procedures semi-implicit algorithm for solving fluid-structure interaction problems, with a structure governed either by linear elasticity equations or by the non-linear Saint-Venant Kirchhoff model. In the second chapter, we present a partitioned procedures semi-implicit algorithm for solving fluid-structure interaction problem with a linear model for the structure and we prove an unconditional stability result of the algorithm. In the above algorithms, an optimization problem must be solved in order to get the continuity of the velocity as well as the continuity of the stress at the interface. During the iterations of BFGS for solving the optimization problem, the fluid mesh does not move and the fluid matrix is only factorized once, which reduces the computational effort. In the fast chapter, we present a monolithic semi-implicit algorithm for solving fluid-structure interaction problem with linear model for the structure. The algorithm uses one global mesh for the fluid-structure domain. The continuity of velocity at the interface is automatically satisfied and the continuity of stress does not appear explicitly in the global weak form due to the action and reaction principle. At each time step, we have to solve a monolithic system of unknowns velocity and pressure defined on the global fluid-structure domain. When the monolithic approach is used the CPU time is reduced compared to a particular partitioned procedures strategy.

L'os humain est à la fois un matériau (tissu osseux) et une structure (e.g., le fémur). Les traumatismes engendrent des défaillances (fractures) structurelles évidentes de l'os, mais l'intégrité mécanique de celui-ci peut aussi être atteinte d'une manière insidieuse, et non moins dangereuse, par certaines maladies. Très schématiquement, on peut dire que les traumatismes sont la cause de macrofractures (à la structure-os), et les maladies la cause de microfractures (du matériau-os). Celles-ci se développent progressivement en macrofractures, et si des soins ne sont pas prodigués ou efficaces, l'os perd sa fonction de soutien. Cette thèse concerne essentiellement une des maladies de l'os : l'ostéoporose. Pour traiter cette maladie, et/ou prévenir les macrofractures dont elle est à terme la cause, il faut d'abord en faire le diagnostic. Ce problème est compliqué du fait qu'il s'agit d'une caractérisation (essentiellement mécanique) de matériau vivant et qu'il est impératif que le sondage soit de type non-destructif, surtout si l'examen doit se répéter souvent (notamment pour suivre les progrès d'un traitement). Le travail de cette thèse concerne donc le développement de nouvelles méthodes, et/ou l'amélioration de méthodes plus anciennes relatives à l'évaluation non-destructive (END) d'altérations mécaniques du tissu osseux en rapport avec l'ostéoporose. Comme l'END est aussi un terme employé pour la caractérisation de matériaux inertes, le domaine d'application de ce travail dépasse celui des matériaux vivants.

L'exploration pétrolière en eaux profondes crée des nouveaux défis technologiques. Certains de ces problèmes sont liés à des très longues structures servant à relier la plate-forme à la tête de puits. L'augmentation de la profondeur implique des structures plus longues et, par conséquent, plus souple. Deux problèmes liés à ce type de structure sont étudiés dans ce rapport : le positionnement de la structure pour sa connexion à la tête de puits (opération de rentrée), et les vibrations induites par vortex. La première partie de ce rapport est consacrée au contrôle de l'opération de rentrée. Deux stratégies différentes de planification du mouvement sont proposées. La première est basée sur un modèle obtenu par approche modale. L'autre stratégie approxime le comportement de la structure par l'équation d'un câble de Bernoulli avec amortissement. Les solutions approximées de cette EDP sont directement utilisées dans la conception du contrôleur. Deux types de contrôle en boucle fermée sont proposées: un système de suivi de trajectoire en utilisant une fonction de Lyapunov, et un autre qui utilise l'inverse du système. La deuxième partie de ce rapport est consacrée à l'étude des vibrations induites par vortex (VIV) et à leur contrôle. Elle présente une première stratégie de contrôle afin de réduire les VIV. Cette stratégie est basée sur une analyse modale des équations du système. La loi de commande génère un déplacement en opposition de phase par rapport au VIV, l'atténuation résultant d'un effet d'antirésonance. Les résultats de simulation sont confirmés par des expériences réalisées sur une maquette en modèle réduit.

La méthode des mobilités énergétiques a été développée dans le domaine des moyennes fréquences, afin de prévoir le comportement vibratoire des structures par bandes de fréquence. Dans ce travail, on étend l'approche du cas du couplage vibroacoustique entre un fluide et une structure. Les mobilités acoustiques et structurales sont étudiées pour le milieu fluide et la structure respectivement. La surface du couplage est maillée en pavés, afin de réduire le problème en mobilités énergétiques ponctuelles. Une simulation numérique a été réalisée dans le cas d'une plaque couplée à une cavité parallélépipédique, où les influences du nombre de pavés et de l'amortissement des sous-systèmes sont étudiées. L'énergie cinétique de la plaque, la puissance transmise au milieu acoustique et le facteur de rayonnement de la structure sont calculés et comparés avec un calcul de référence (méthode modale) et avec des résultats expérimentaux. Les critères de maillage (taille de pavé par rapport aux longueurs d'ondes de la structure et du milieu fluide) sont déterminés avant et après la fréquence critique, pour assurer une modélisation correcte du couplage dans une bande de fréquence donnée. Afin de tenir compte d'une structure compliquée à modéliser, un calcul mixte (utilisation des mobilités mécaniques mesurées et des mobilités acoustiques calculées) est également présenté.

Ce travail de recherche traite du couplage entre un liquide incompressible, irrotationnel et son contenant : une structure élastique. Cette interaction fluide-structure est traitée dans le cadre des petites déformations autour d'un état d'équilibre.Dans un premier temps, on présente une méthode d'introduction des sources dissipatives visqueuses dans le liquide à partir des équations du système couplé conservatif en s'appuyant sur une approche de type fluide potentiel généralement utilisée pour traiter les problèmes de couplage fluide-structure linéarisés non amortis. Un modèle d'amortissement diagonal est alors choisi pour le liquide et les effets dissipatifs de celui-ci sont pris en compte en calculant les coefficients d'amortissement modaux. Seuls les effets dissipatifs liées à la viscosité du liquide sont alors pris en compte. Le système couplé dissipatif obtenu possède une matrice d'amortissement non symétrique. Une résolution de ce système à amortissement non classique est alors présentée et les expressions des réponses fréquentielle et temporelle linéarisées sont données pour différents types d'excitations.Dans un deuxième temps, le liquide est supposé non visqueux et les forces de tension surfacique sont prises en compte. Cette configuration concerne principalement les satellites où le système couplé est en situation de microgravité. Une formulation du problème conservatif permettant de prendre en compte l'incompressibilité du fluide, la condition de continuité à l'interface fluide structure, les effets de capillarité du fluide ainsi que les effets éventuels de précontraintes statiques est alors établie. On se propose pour cela d'utiliser une méthode énergétique via le Principe de Moindre Action. La démarche est alors décomposée en deux étapes : une étude statique afin de déterminer la position de référence, puis une étude dynamique linéarisée autour de cette position d'équilibre. Cette formulation forme notamment une base pour l'introduction des sources dissipatives liées aux effets de capillarité via la méthode précédemment introduite.

Mémoire numérisé par la Division de la gestion de documents et des archives de l'Université de Montréal

Les travaux présentés dans ce mémoire s'articulent essentiellement autour de deux thèmes : l'étude mathématique et numérique des phénomènes d'interaction fluide-structure et la modélisation de l'appareil respiratoire humain. C'est au cours de ma thèse que j'ai commencé à travailler sur les phénomènes d'interaction fluide-structure. Ces travaux se sont poursuivis par la suite au CEREMADE (Université Paris Dauphine) puis au sein du projet REO à l'INRIA. Je me suis essentiellement intéressée au cas des fluides newtoniens, visqueux, incompressibles satisfaisant aux équations de Navier-Stokes. La structure est, quant à elle, mobile : rigide ou déformable. On se place dans le cas ou les déplacements de la structure ont une amplitude telle que l'on ne peut pas les supposer infinitésimaux, ce qui va faire apparaitre des non linéarités géométriques. J'ai étudie ces problèmes aussi bien d'un point de vue théorique (existence de solutions faibles ou fortes pour des systèmes stationnaires ou instationnaires) que d'un point de vue numérique (étude de la stabilité et de la convergence des schémas, conception d'algorithmes performants). Les applications ou de tels phénomènes apparaissent sont nombreuses : tant en aérodynamique qu'en biomécanique. Nous pourrions citer beaucoup d'exemples : écoulement autour d'un bateau, écoulement sanguin dans les artères, pour l'hydrodynamique (fluide en phase liquide), écoulement autour d'ailes d'avion, étude de l'influence des vents sur le tablier d'un pont, pour l'aéroélasticité (fluide en phase gazeuse). En particulier, ils apparaissent dans l'étude des écoulements biologiques : écoulement sanguin dans les artères, écoulement de l'air dans l'appareil respiratoire. Ces dernières années je me suis particulièrement intéressée à ce type d'applications et plus spécifiquement à la modélisation mathématique et numérique de l'appareil respiratoire, tout d'abord dans le cadre d'une ACI nouvelle interface des mathématiques le poumon vous dis-je, puis au sein du projet REO. Je m'intéresse à la modélisation de l'écoulement de l'air dans les voies aériennes, la modélisation du transport et dépôt de particules ainsi qu'à la modélisation des tissus élastiques pulmonaire. Ces applications nécessitent le développement de modèles mathématiques nouveaux et de méthodes numériques spécifiques et adaptées.

Une méthode innovante pour simuler des interactions fluide-structure complexes tout en gardant un bon compromis temps de calcul/précision est présenté.Pour réduire le temps de simulation des modèles d’ordre réduits sont utilisés au lieu des modèles complets aussi bien pour les modèles structuraux que pour les modèles aérodynamiques. Un des challenges de base était d'utiliser des modèles industrielles hautes fidélités. La technique de condensation dynamique est utilisée pour réduire la taille du modèle éléments finis structures et la décomposition aux valeurs propres est utilisé sur une base de données aérodynamiques construite à partir de simulations CFD.Les non-linéarités structurelles sont réintroduites à posteriori.Une comparaison poussée des méthodes classique d'interpolation comme des méthodes de spline, d’interpolation sur des Manifold de Grassmann avec des méthodes innovantes d'apprentissage statistiques a été amené.Afin de valider complètement la méthodologie développée, une maquette expérimentale visant à imiter le comportement du carénage au sol avant le décollage a été conçue.Ce cas a pu être assimilé à une plaque avec des raideurs de liaisons dans une couche de mélange.La validation de cette méthode est réalisée en comparant les résultats des simulations numériques avec les données enregistrées pendant des essaies en soufflerie. On pourra ainsi comparer aussi bien des champs que des mesures locales. L'ensemble des essais a permis d'améliorer la compréhension de ce phénomène vibratoire qui mène à des problèmes récurrents de fatigue dans cette sous structures.Cette méthode est enfin appliquée à une structure aéronautique: les carénages de volet hypersustentateur / An innovative method for numerical simulating complex problems of fluid structure interaction, such as non-linear transients, characterized by good performances and high precision is presented in this manuscript. To cut down the simulation time, reduced order models are used for both the aerodynamic and structural modules. High fidelity industrial models have been used. A technique of dynamic condensation is employed to reduce the size of the finite element model while the technique of Proper Orthogonal Decomposition is used on a database of aerodynamic pressures built from CFD simulations. Structural non-linearities are reintroduced a posteriori. Different interpolation techniques such as the classic spline interpolation, interpolation on a Grassmann Manifold with more innovative methods of statistical learning have been compared. In order to validate the developed methodology a test campaign has been designed to reproduce a simplified mechanism of interaction inspired by a flap track fairing in take-off configuration. A plate whose stiffness depends on the springs at its attachment to the wind tunnel test section floor is immersed in a mixing layer. In parallel to the test activities a numerical model of the test rig has been developed. The validation of the methodology of fluid structure interaction is done through direct comparison between test data and simulation results. The testing activities have granted a deeper comprehension of the vibratory phenomenon that has led to recurrent fatigue problems on the impacted structures. The methodology is ultimately applied to an industrial problem: the load prediction on flap track fairings excited by engine exhaust.