UPThèses

Cette thèse est une contribution expérimentale à l’étude des résonances aéroacoustiques des jets sous-détendus : le Screech. Diverses méthodes expérimentales sont utilisées à ces fins, telles que la mesure de pression acoustique, la strioscopie et Vélocimétrie par Image des Particules, et associées à des techniques classiques de post-traitement comme les décompositions en mode de Fourier et aux valeurs propres. Ces Techniques permettent d’évaluer les effets d’épaisseur de la lèvre de la buse sur l’écoulement, et fournissent des informations sur les différences de comportement d’un même jet montrant des modes oscillatoires différents. Enfin, on entreprend d’étudier la présence de divers mécanismes de fermeture de la boucle de résonance pour divers modes de screech. La présence d’ondes intrinsèques du jet, se propageant vers l’aval pour les modes axisymétrique (A2) et hélicoïdal (C) suggèrent que ces ondes puissent jouer un rôle dans la résonance. La signature de ces ondes n’est en revanche pas attestée pour les modes battants (B). Ces résultats semblent donc indiquer que plusieurs mécanismes de rétroaction différents puissent être à l’oeuvre dans la résonance du jet sous-détendus.

Cette thèse est consacrée à l’étude des bulles d’hélium dans le silicium. Dans un premier temps, des simulations atomistiques de type dynamique moléculaire (DM) ont été réalisées afin d’étudier les propriétés physiques d’une bulle nanométrique : densité et pression d’hélium en fonction du diamètre de la bulle et de la quantité d’hélium initialement introduite dans la bulle. Nos calculs ont montré que la déformation plastique de la matrice ne permettait pas de dépasser une valeur de densité d’hélium nettement inférieure aux résultats expérimentaux. En modélisant un système plus réaliste, un excellent accord avec les mesures expérimentales a été retrouvé. Nous avons aussi montré que le choix du potentiel interatomique était capital afin de modéliser des bulles sous fortes pressions. En outre, nos résultats ont montré que la loi de Laplace-Young n’était pas adaptée pour prédire la pression dans les bulles d’hélium dans du silicium. Enfin, l’étude de l’état de l’hélium dans la bulle indique que l’hélium peut être sous forme solide à 300 K, du fait des pressions élevées. Dans un second temps, une approche de type dynamique d’amas (DA) a été employée afin de comprendre les premiers stades d’évolution des bulles. Nous avons construit un modèle applicable à la DA à l’aide d’un grand nombre de calculs DM. Une étude paramétrique a ensuite été réalisée. Ceci nous a permis d’étudier les mécanismes de croissance des bulles, en prenant en compte l’influence de l’hélium et la présence ou non d’une source continue de défauts.

Cette étude vise à caractériser en soufflerie les mécanismes aérodynamiques à l’origine de la génération puis la transmission des fluctuations de bruit dans un modèle d’habitacle automobile. Le banc d’essai conçu en soufflerie anéchoïque consiste en un écoulement dont la vitesse incidente est modulée par un volet mobile, et qui par interaction avec une marche montante rayonne un bruit aéroacoustique transmis à travers une vitre dans un caisson anéchoïque. L’approche retenue consiste, pendant le temps de maniement du volet, à mesurer et relier le champ de vitesse externe mesuré à l’aide de la technique de vélocimétrie laser par images de particules (TR-PIV échantillonnée à 20 kHz) à la pression pariétale d’une part, puis au champ acoustique interne obtenu par transmission d’autre part. Des outils de corrélation spatio-temporelle sont alors utilisés pour mettre en évidence les zones de l’écoulement les plus corrélées avec les fluctuations d’énergie de la pression pariétale et celles du niveau de bruit intérieur. La fluctuation du chargement aérodynamique de la vitre sous la bulle de recirculation est logiquement liée à l’activité instationnaire de cette dernière, puis plus en aval, au lâcher tourbillonnaire. Quant au bruit transmis dans le modèle d’habitacle, il semble principalement lié aux fluctuations de vitesse dans la couche de cisaillement. Enfin, une procédure spécifique a permis d’évaluer le caractère quasi-stationnaire des variations temporelles des quantités fluctuantes ainsi que la réponse acoustique de la vitre.

Ce travail de recherche présente une étude expérimentale de l’influence d’un écoulement tourbillonnaire sur un lit sédimentaire. Les expériences se déroulent dans un canal hydro-sédimentaire à surface libre, à l’Institut P’. Deux dispositifs permettant de générer des tourbillons isolés, d’axe vertical ou horizontal, ont été conçus lors de cette thèse. Ce travail se concentre sur le premier de ces montages : une étude paramétrique par mesures optiques (PIV, PIV stéréoscopique) est menée, notamment avec deux conditions de rugosité de paroi. Des lois de comportements sont ainsi établies. La caractérisation des écoulements tourbillonnaires révèle une forte tridimensionnalité, de même que la présence de structures turbulentes particulières, liées aux phénomènes de transport sédimentaire. Un second axe d’étude s’attache à la description de la mise en suspension de sédiments par un tourbillon d’axe vertical. Les structures turbulentes liées au transport sont étudiées, ainsi que le comportement du panache de sédiments. La reconstruction du lit sédimentaire sur un temps long, par l’utilisation d’une méthode de stéréo-corrélation, met en avant les régions d’affouillement et d’accrétion induites par le passage d’un seul tourbillon isolé. Enfin, une étude numérique basée sur OpenFOAM est proposée. La méthodologie exposée reproduit bien les comportements du fluide et des sédiments : la validation de cet outil permettrait d’aborder de nouveaux paramètres d’étude par la connaissance précise des quantités de sédiments mobilisés.

Les paliers à feuilles sont des organes de guidage en rotation adaptés pour des rotors légers fonctionnant à des grandes vitesses de rotation. Leur fonctionnement ne nécessite aucun apport d’huile ni de graisse. Ces paliers sont réalisés avec une structure compliante. La présence de frottements dans la structure apporte l’amortissement nécessaire au système rotor-palier. Après une présentation de la technologie, un tour d’horizon des modèles existants et des limitations de chacun d’eux est exposé. Par la suite, un modèle non-linéaire du coussinet compliant est introduit. Dans ce modèle les feuilles sont considérées comme des solides élastiques. Les forces de frottement et les jeux entre les feuilles sont pris en compte par l’utilisation de la condition de non-interférence de Moreau-Signorini. Les efforts normaux et les forces de frottement sont calculés respectivement par la méthode des multiplicateurs de Lagrange augmentés et la méthode des pénalités. Ensuite, cette structure compliante est couplée au film fluide, traité par la résolution de l’équation de Reynolds en conditions de lubrification mixte. L’intérêt est ensuite accordé au palier dans sa globalité : structure déformable, film fluide, rotor. En régime statique, cette étude passe par l’analyse des démarrages du rotor puis de son fonctionnement à des hautes vitesses de rotation et sous de très fortes charges statiques. L’effet des erreurs d’usinage est mis en évidence. En régime dynamique, l’étude se fait par l’analyse non-linéaire d’un rotor de Jeffcott supporté par des paliers à feuilles. Les résultats montrent les limites de stabilité du système rotor-paliers et l’influence du balourd.

Les joints sans contact sont des solutions d’étanchéité optimales pour les systèmes mécaniques fonctionnant à des vitesses relativement élevées. Ils ont une durée de vie importante et ont été proposés pour une utilisation dans les moteurs spatiaux. Une type de joints d'étanchéité sans contact, le joint visqueux, est étudié en détail dans ce travail de thèse au moyen d'une analyse numérique et expérimentale. Ce travail présente un modèle numérique développé pour prédire le comportement du joint visqueux en régime laminaire et turbulent. L’interface liquide-air dans le joint est également étudiée en utilisant une approche CFD basée sur la méthode VOF. La conception et la réalisation d’un dispositif expérimental a permis de confronter les résultats numériques et expérimentaux, l’écart entre les deux approches n’excède pas 10% dans 95% des cas étudiés.

Les matériaux absorbants acoustiques, qui sont d’un intérêt stratégique en aéronautique pour la diminution passive du bruit des réacteurs d’avion, conduisent à une physique complexe où l’écoulement turbulent, des ondes acoustiques, et l’absorbant interagissent. Cette thèse porte sur la simulation de cette interaction dans le problème modèle d’un écoulement de canal turbulent avec des parois impédantes, par le biais de simulations numériques aux grandes échelles implicites, dans un contexte de calcul haute performance. Une étude est d’abord faite des grandes échelles dans un canal turbulent avec des parois rigides, en s’intéressant plus particulièrement à l’effet d’une faible compressibilité (Mach <3) sur les caractéristiques de ces échelles. Un canal turbulent avec une paroi de type impédance est ensuite simulé, avec une condition habituelle de périodicité dans le sens de l’écoulement. On observe que pour des faibles valeurs de la résistance et des fréquences de résonance basses, l’écoulement est instable, ce qui engendre une onde le long de l’absorbant, qui modifie la turbulence et augmente la trainée. Enfin, on se tourne vers une simulation de canal spatial en levant la condition de périodicité dans la direction de l’écoulement, ce qui permet d’introduire une onde acoustique en entrée de domaine. L’atténuation de l’onde dans l’écoulement turbulent est étudiée avec des parois rigides, puis un absorbant acoustique est introduit. Dans cette configuration plus réaliste, il est confirmé que l’écoulement peut devenir instable au bord amont de l’absorbant, ce qui empêche l’atténuation de l’onde acoustique incidente.

Les traitements antireflets (AR) sont très largement utilisés pour améliorer la transmission de systèmes optiques composés de hublots, lentilles, de lames séparatrices,… Dans cette thèse les gammes spectrales visées sont le visible étendu [400-1800nm] et le moyen infrarouge [3,7-4,8µm]. La méthode de nanostructuration par dépôts de films minces utilisant des techniques PVD en incidence oblique (Oblique Angle Deposition) a été choisie car elle permet d’envisager des AR hautes performances sur une large gamme de longueur d’onde, via un procédé industrialisable. L’introduction de porosité via le contrôle des angles de dépôt est utilisée pour nanostructurer l’architecture de chaque couche et de l’empilement ; méthode permettant de modifier et d’optimiser les propriétés optiques des couches constituantes en vue d’un design complet optimal. Une cartographie des indices effectifs accessibles par OAD a été dégagée concernant les trois matériaux déposés (TiO2, SiO2 et Ge). Mais les propriétés optiques de ces couches nanostructurées diffèrent largement de celles des couches denses du fait de la présence d’anisotropie, de gradient d’indice, de diffusion et d’absorption. A partir de caractérisations microstructurales, chimiques et optiques poussées (AFM, MEB, MET, tomographie FIB, tomographie MET, EDX, EELS, spectrophotométrie et ellipsométrie généralisée) un modèle optique analytique plus complexe et couplé à des analyses par éléments finis (FDTD) est présenté. L’ensemble du travail a permis d’élaborer par OAD de simples antireflet bicouches démontrant déjà de hauts niveaux de transmission, supérieurs aux traitements AR existants (interférentiel) ou en développement (Moth-eyes).

Les microfils métalliques polycristallins produits par étirage à froid présentent une résistance mécanique significative en faisant des candidats idéaux pour les renforts de composites. Des études antérieures sur des fils de nickel polycristallin pur ont montré une dépendance importante par rapport à la taille de la limite d'élasticité et de la résistance à la traction, ainsi que de la ductilité. Le but de cette étude est de comprendre cet effet de la taille dans les microfils de nickel pur polycristallin par analyse de diffraction des rayons X in-situ (DRX) et simulations de la plasticité cristalline par éléments finis (CPFE). Des essais de traction monotone et cyclique in-situ sous rayonnement synchrotron ont été réalisés sur des microfils de diamètres allant de 100 à 40 μm. Les fils étirés à 100 micromètres obtenus dans le commerce présentent une architecture cœur-coquille avec une texture de fibre <111> dominante dans le cœur et une texture à double fibre hétérogène <111> et <100> dans la coquille. La réduction de la taille de l'échantillon par polissage électrolytique conduit à des fils ayant une microstructure homogène, tandis que la réduction de la taille de l'échantillon par un étirage à froid supplémentaire conduit à des fils avec une texture plus intense tout en conservant l'architecture cœur-coquille. La limite d'élasticité et la résistance à la traction des fils électropolis augmentent avec la diminution du diamètre, tandis que la ductilité diminue avec la réduction du diamètre. Dans le cas des fils étirés à froid, on observe que la limite d'élasticité et la résistance à la traction, ainsi que la ductilité, augmentent avec la diminution du diamètre. L'analyse DRX indique une plasticité successive des familles de grains sous iso-déformation. Nous avons observé que le gradient de la texture du microfil active des mécanismes de déformation qui ne sont pas observés pour les microfils à texture homogène. Pour comprendre l'influence de différents paramètres microstructuraux, notamment l'influence de la texture cristallographique, une microstructure représentative 3D a été générée et des simulations CPFE ont été réalisées. Le comportement simulé moyen des différentes familles de grains (<111>, <100>) concorde bien avec les résultats expérimentaux. La simulation CPFE indique une hétérogénéité du champ de contrainte à travers la microstructure en présence d'un gradient de texture cristallographique. Nous montrons que la micro-texture (texture simple ou double texture) et leur dispersion spatiale (homogène ou architecturée) peuvent être utilisées comme stratégie de conception pour obtenir une microstructure optimale en fonction de l’ensemble désiré de propriétés mécaniques.

Une étude de l'impact des mascarets sur le transport des sédiments à l'aide de la simulation numérique a été réalisée dans ce travail. En utilisant le logiciel OpenFOAM CFD, nous avons généré 17 simulations numériques de mascaret avec divers nombres de Froude Fr, allant de 0,99 à 1,66. Deux types de mascarets, ondulant et déferlant, ont été couverts dans ces 17 simulations numériques. Pour les particules sédimentaires non cohésives, nous avons utilisé les équations de Maxey et Riley pour déterminer la trajectoire des particules sédimentaires non cohésives sous l’influence d’un mascaret ondulant. En utilisant le schéma Runge-Kutta du quatrième ordre, une méthode tracker résout les équations de Maxey et Riley qui nécessitent l’information des champs de vitesse au temps t. Pour les particules sédimentaires cohésives, nous avons calculé la distribution des particules sédimentaires cohésives en utilisant un modèle de transport de flocs, présenté par Winterwerp (2001). Dans ce modèle, la concentration volumique solide des sédiments et le diamètre des flocs D sont estimés. Les équations de transport de et D sont résolues en utilisant la méthode des moments présentée par Beaudoin et al. (2002 et 2004). La méthode des moments permet de réduire le temps CPU rendant possible une étude paramétrique. De ce travail, nous avons trouvé une classification du mascaret en fonction du nombre de Froude Fr. Cette classification est également basée sur l’étude menée par Furgerot (2014). Pour un nombre de Froude 1,04 < Fr < 1,43, le mascaret est ondulant. Pour un nombre de Froude 1,43 < Fr < 1,57, le mascaret est partiellement déferlant, similaire à la transition de mascaret définie par Furgerot (2014). Pour un nombre de Froude Fr > 1,57, le mascaret est totalement déferlant. Une analyse de la distribution de la pression a été effectuée par Baddour et Song (1990). Nous avons trouvé que les pressions totales et hydrostatiques d’un mascaret ondulant ont de grandes valeurs sous la crête et le creux. Dans le cas d’un mascaret ondulant, les pressions totales ne sont pas égales à les pressions hydrostatiques. Cela provoque la présence de pressions dynamiques. Dans le cas de mascaret déferlant, les pressions totales deviennent égales aux pressions hydrostatiques. La turbulence réduit les pressions dynamiques. L'impact des mascarets sur le transport de particules sédimentaires non cohésives et cohésives a été étudié dans ce travail. Pour les particules sédimentaires non cohésives, nous avons observé que la trajectoire utilisant l’écoulement généré par OpenFOAM est similaire à la trajectoire de type e proposée par Chen et al. (2010). Des modifications du modèle de Chen ont été faites en incluant les effets de la gravité, l’élévation et l’atténuation pour reproduire des trajectoires de particules non cohésives sous un mascaret ondulant. Nous avons obtenu des relations linéaires entre les paramètres du modèle de Chen modifié (β1 , β2 et β3) et le nombre de Froude Fr. C’est parce que le niveau de la turbulence du mascaret ondulant est faible. L’écoulement induit par mascaret ondulant n’est pas complexe. Ce phénomène physique est quasi linéaire. Le paramètre β1, lié à la célérité avant du mascaret ondulant, diminue lorsque le nombre de Froude Fr augmente. Les paramètres β2 et β3, liés respectivement à l’élévation et à l’atténuation du mascaret déferlant, augmentent lorsque le nombre de Froude Fr augmente. Enfin, pour les particules sédimentaires cohésives, nous avons calculé la distribution de la taille des flocs D sous deux types de mascaret, ondulant et déferlant. Nous avons utilisé le diamètre initial de la particule sédimentaire cohésive d = 4 μm. La taille du floc initial D est égale à 10 μm avec la concentration du floc c = 0,5 kg/m3 . Et nous avons limité la taille maximale du floc à 2000 μm. Nous avons observé que la valeur maximal de la taille des flocs Dmax augmente de façon exponentielle par rapport au nombre de Froude Fr.