Réduire la traînée d'un avion permet de diminuer les coûts du transport aérien, enjeux des plus stratégiques. Dans ce but, le maintien d'une couche limite laminaire étendue permet une forte réduction de la traînée de frottement. Les progrès réguliers dans ce domaine ont débouché sur l'optimisation des profils d'aile, et sur la sensibilisation de la couche limite par aspiration au travers de parois poreuses. Ces techniques ont été testées en vol sur la voilure d'un Falcon 900, avec Dassault Aviation, faisant apparaître des écarts entre prévisions et mesures de transition. Ces écarts ont été attribués à des défauts de géométrie (ondulations) de la voilure, et ont motivé un important travail d'analyse. Il a été décidé de développer une nouvelle génération de critères de tolérance de forme, basés non plus sur des corrélations empiriques globales, mais sur les analyses de stabilité de la couche limite. La théorie de la stabilité permet, en effet, de caractériser la croissance des instabilités, la mesure de cette croissance se faisant au moyen du facteur N. Une modélisation du saut de facteur N en fonction des paramètres caractérisant la déformation géométrique permet de prévoir, à partir de l'étude de la géométrie de base, quels seront les effets sur la transition. La démarche employée consiste en une série d'études expérimentales et numériques, de complexité croissante, permettant d'affiner progressivement les lois de variations physiques et leur modélisation. L'ensemble de ces travaux a permis l'élaboration d'un critère tenant compte de l'ensemble des paramètres étudiés et permettant une bonne estimation de la position de transition pour les différentes configurations. Enfin des déformations tridimensionnelles ont été examinées ; cette étude à mis en évidence, pour certaines configurations, la prépondérance des effets tridimensionnels.

​ Nombreuses stratégies de contrôle ont été récemment proposées par la communauté scientifique afin depouvoir réduire la traînée dans les écoulements pariétaux. Entre autres, les Surfaces Superhydrophobes (SHS) ontmontré leurs capacités de pouvoir réduire considérablement le frottement pariétal d'un écoulement liquide grâce à laprésence de microbulles de gaz piégées dans les nano-rugosités de la surface. Dans des conditions géométrique etthermodynamique données pour lesquelles la transition de mouillage est évitée (condition pour laquelle normalementla taille des rugosités qui caractérise la SHS est de plusieurs ordres de grandeur plus petite que l'échellecaractéristique de l'écoulement principal), on peut atteindre ce qu'on appelle 'l'effet Lotus', pour lequel l'écoulementglisse à la paroi, avec une vitesse différente de zéro.. Dans ce cadre, nous nous sommes proposés d'étudier, à l'aidede simulations numériques l'influence des SHS sur la transition laminaire-turbulent dans un écoulement de canal.Pour cela, nous avons réalisé une série de simulations numériques directes (DNS), allant de l'état laminaire au casturbulent pleinement développé, en traitant la plupart de scénarios de transition connu en littérature. Des analyses destabilité locale et globale ont aussi été réalisées afin de déterminer l'influence de ces surfaces sur la première phasedu processus de transition. Bien que la procédure de déclenchement de la transition contrôlée (type K, H, C,...) soitbien décrite dans la littérature, cela n'est pas le cas pour les transitions naturelles. À cette fin, une nouvelle méthode aété développée pour déclencher puis étudier la transition naturelle dans des écoulements de type canal. Cette méthodeest basée sur des mécanismes de réceptivité de l'écoulement (resolvent global) permettant de construire un forçagevolumique spécifique. Plusieurs approches pour modéliser les SHS ont été utilisées, de complexités croissantes, touten tenant en compte des caractéristiques physiques de ces surfaces. Dans un premier temps, une condition deglissement homogène a été utilisée et son influence analysée. Chaque rugosité a été ensuite discrétisée spatialement,d'abord avec une alternance de condition limite sur une surface plate, ensuite en tenant compte de la dynamique del'interface gaz-liquide par une méthode Lagrangienne-Eulerienne Arbitraire (ALE). Nous avons montré que les SHSpermettent d'efficacement retarder les transitions contrôlées mais qu'en revanche elles ont peu d'influence sur lestransitions naturelles (développant des stries de vitesse). En effet, ce comportement dérive de l'équilibre entre deuxeffets contradictoires. D'un côté, le glissement pariétal nuit au développement des structures cohérentes de typehairpin ​ , en altérant le processus de ​ vortex stretching-tilting ​ . D'autre part, le mouvement de l'interface gaz-liquideinteragit avec les structures cohérentes de l'écoulement, en produisant des vitesses normales à la paroi favorisantdavantage le processus de ​ sweep-ejection et entraînant le développement de structures en forme d'arche. Nous avonsmontré que les interfaces gaz-liquide statiques retardent la transition de façon analogue à une condition aux limiteshomogène (si l'hétérogénéité pariétale est petite). En revanche la prise en compte de leur dynamique limite le retardde la transition, montrant l'importance du modèle de SHS dans les écoulements transitionnels.

La thèse propose une étude numérique et une modélisation de la transition vers la turbulence dans l'écoulement de Couette plan. On se concentre particulièrement sur le régime de bandes oblique laminaire turbulentes, au coeur de la transition. On suit deux types d'approches selon le type d'échelles considérées. D'une part on étudie le détail à petite échelle de l'écoulement. A l'aide de DNS, une instabilité de cisaillement créant de la vorticité transverse est mise en évidence, et la vitesse d'advection des perturbations est mesurée. La modélisation de l'écoulement turbulent et cette instabilité vient confirmer ces mécanismes et met en évidence la transition absolue/convective à l'entrée des bandes. une description du maintien de la turbulence en termes de cycle auto-entretenu est proposée. D'autre part, on étudie l'écoulement 'a l'échelle de la bande. Des outils statistiques permettent d'étudier les bandes en termes de formation de motif. L'effet du bruit intrinsèque issu de la turbulence est mis en évidence, en particulier dans les changements d'orientation. La modélisation du motif en termes d'équation de Landau-Langevin et d'équation de Fokker-Planck correspondante fait le parallèle avec les motifs bruités et les phénomènes critiques. Partant de ce modèle on peut comprendre les retournements en termes de temps de premier passage moyen. Finalement, on regarde deux autres types de problèmes en couplant les approches microscopiques et macroscopiques. L'étude des relaminarisations tant au seuil de disparition de la turbulence qu'au seuil d'apparition du motif permet de mettre en évidence les mécanismes de disparition de la turbulence. L'étude de spots turbulents propose un aperçu de leur zoologie et permet d'approcher les mécanismes petite échelle à plus grand nombre de Reynold

Lorsque des ondes d'amplitudes suffisamment importantes se propagent dans un milieu, elles peuvent interagir et donner naissance à des ondes de longueurs d'onde différentes. L'énergie est transportée par interactions non linéaire et une cascade d'énergie peut alors avoir lieu sur une large gamme d'échelles spatiales, d'une échelle d'injection à une échelle de dissipation. Cette phénoménologie, présentant des analogies formelles avec la turbulence hydrodynamique peut se rencontrer dans de nombreux systèmes physiques tels que les ondes de gravité et de capillarité à la surface de la mer, les ondes internes dans l'océan ou l'atmosphère, les ondes élastiques, les ondes d'Alfvén dans le vent solaire... Dans cette thèse, la dynamique et la statistique d'un ensemble d'ondes de surface en interaction est étudiée dans différents systèmes : les ondes capillaires et de gravité à la surface d'un liquide ainsi que les ondes hydro-élastiques à la surface d'une feuille élastique flottante. Les résultats expérimentaux et numériques sont comparés aux prédictions de la théorie de turbulence d'ondes, aussi appelée turbulence faible. Dans un premier temps, la turbulence d'ondes gravito-capillaires à la surface d'un liquide est étudiée, à la fois en régime stationnaire et instationnaire. Le déclin auto-similaire en temps de la turbulence d'ondes capillaires en déclin est observé expérimentalement, en accord avec le scénario théorique. De plus nous mettons en évidence l'existence de dissipation au sein de la cascade. L'influence de la dissipation à toutes les échelles sur les cascades de turbulence d'ondes est alors étudiée en régime stationnaire : pour une faible dissipation, la dépendance du régime de turbulence d'ondes avec l'échelle est trouvée en accord avec la théorie, cependant nous mettons en évidence un flux d'énergie non constant ; pour une forte dissipation, le régime de turbulence d'ondes est en fort désaccord avec la théorie. Par ailleurs, les premières simulations numériques de turbulence d'ondes capillaires à partir des équations diphasiques de Navier-Stokes (code Gerris) sont présentées. Un spectre d'énergie des vagues en loi de puissance est observé en bon accord avec la théorie. Les cascades inverse et directe de la turbulence d'ondes de gravité sont ensuite étudiées. Nous présentons la première observation expérimentale de la cascade inverse des petites vers les grandes échelles. Les effets de taille finie et l'influence des conditions aux limites sur la cascade directe d'ondes de gravité sont discutés à l'aide d'expérience dans des bassins de grandes tailles. Dans une seconde partie, je m'intéresse aux ondes hydro-élastiques à la surface d'une feuille élastique flottante sur un liquide. Deux régimes d'ondes élastiques sont observés : des ondes de tension et de flexion (respectivement à basse et haute fréquence) en bon accord avec la théorie linéaire. La valeur de la tension de la membrane est contrôlée par la pression hydrostatique imposée sur le système. Lorsque l'amplitude de forçage est augmentée, un décalage significatif de la relation de dispersion non linéaire des ondes est observé, due à une surtension créée par les oscillations lentes du mode fondamental de la feuille. Le processus d'interactions à trois ondes est mis en évidence expérimentalement et un régime de turbulence d'ondes est observé. Un spectre d'énergie des vagues en loi de puissance est toujours observé, avec un exposant en désaccord avec les prédictions théoriques. Différentes hypothèses sont discutées afin d'expliquer ce désaccord, comme l'existence de dissipation à différentes échelles au sein de la zone inertielle, ou l'influence des fortes non linéarités du mode fondamental

Le travail présenté dans ce mémoire, divisé en deux parties principales, traite de la transition laminaire-turbulent de la couche limite en écoulement incompressible sur des ailes en flèche. La première partie concerne les problèmes de contamination et de transition naturelle de la ligne de partage alors que la seconde est relative à la transition par instabilités (longitudinales ou transversales) sur une aile en flèche. Les conclusions dégagées au cours de ce travail s'appuient sur des résultats numériques et expérimentaux obtenus dans la soufflerie subsonique F2 du Fauga-Mauzac. L'étude de la ligne de partage a permis de quantifier les performances d'un Gaster bump et d'une aspiration pariétale en tant que dispositifs anti-contamination. Des mesures par films chauds montrent que la couche limite sur la ligne de partage reste laminaire jusqu'à des nombres de Reynolds très élevés (proche des valeurs sur le bord d'attaque d'un avion de transport). Une étude numérique de la transition naturelle de la théorie de la stabilité linéaire est aussi entreprise. Les résultats obtenus sont en bon accord avec les observations expérimentales. Dans la seconde partie, on s'intéresse aux moyens de prévision de la transition sur une aile en flèche munie d'un dispositif d'aspiration. Les positions de transition et les spectres de fluctuations (observées expérimentalement) fournissent une base de données pour l'exploitation de la méthode du eN (théories liéaires Orr-Sommerfeld et PSE). Deux stratégies d'intégrations sont retenues pour le calcul du facteur N : la méthode enveloppe et la méthode à β fixé. La confrontation entre résultats expérimentaux et numériques montre que la stratégie à β fixé semble la plus appropriée. Enfin, une étude à caractère plus fondamental est entreprise à l'aide du modèle PSE non linéaire.