Étudier les ondes internes est crucial pour comprendre le mélange dans l'océan. Dans cette thèse, un attracteur d'ondes 2D est tout d'abord simulé de manière directe, appuyé par une bonne comparaison avec une expérience pré-existante. Nous dérivons un modèle simple de la largeur de l'attracteur et mettons en évidence des effets non-linéaires. Nous réalisons dans une deuxième partie une étude expérimentale de la réflexion d'ondes planes sur une paroi inclinée. Les résonances prédites entre différents harmoniques n'apparaissent pas mais en revanche, un fort écoulement moyen horizontal apparaît, courbant les caractéristiques des ondes par effet Doppler. 70 à 80% du flux d'énergie incident sont dissipés ou convertis en écoulement moyen, ce dernier semblant alimenté par la dissipation des ondes. La génération d'ondes solitaires consécutive à la réflexion d'ondes sur une pycnocline est ensuite étudiée numériquement dans la troisième partie. Dans un premier temps, une étude académique, 2D est réalisée à l'aide de simulations directes. Nous montrons que des ondes solitaires de différents modes et piégées dans la pycnocline peuvent être générées. Deux critères pour comprendre la sélection d'un mode donné, l'un portant sur les différentes vitesses de phase, l'autre sur des arguments géométriques, sont définis. Ces critères sont dans un second temps comparés aux conditions du Golfe de Gascogne en été. Nous montrons qu'un rayon d'ondes internes seul ne peut générer des ondes solitaires correspondant aux observations, ce qui est corrigé en tenant compte de l'écoulement présent dans la pycnocline et indépendant du rayon d'ondes internes.

La simulation de cas réalistes d'écoulements ou de transferts thermiques implique souvent l'utilisation d'obstacles ou d'interfaces de forme complexe. De part leur manque de flexibilité, les maillages structurés ne sont pas initialement adaptés au traitement d'interfaces irrégulières, ces dernières coïncidant rarement avec les lignes du maillage. Afin de permettre à l'approche structurée de traiter des interfaces complexes avec précision, des méthodes dites de domaines fictifs sont nécessaires. La première contribution de cette thèse est une nouvelle méthode de travail sur maillage curviligne structuré qui permet de réutiliser de nombreuses méthodes fonctionnant initialement sur des maillages cartésiens sur maillages curvilignes. Nous avons ensuite mis au point deux nouvelles méthodes de domaines fictifs : la méthode de pénalisation de sous-maille (PSM) pour la gestion des frontières immergées pour les équations elliptiques et de Navier-Stokes et la méthode d'interface immergée algébrique (IIA) pour les problèmes d'interfaces immergées pour les équations elliptiques. L'un des intérêts de ces deux méthodes à l'ordre deux en espace est leur simplicité. Ces différents développements ont finalement été appliqués à des cas de couplage fluide-structure académiques et réalistes (sédimentation d'un cylindre, hydroplanage d'un pneu, écoulements dans une tête de forage et convection naturelle dans la grotte de Lascaux).

Le principal objectif de cette thèse est de caractériser l'écoulement dans un anévrisme abdominal aortique (AAA) sous différentes conditions physiologiques et à différents stades de son développement. Cette étude est consacrée aux AAA ax- isymétriques, modélisés comme une dilatation de profil gaussien et de section circu- laire. Ainsi, les résultats s'appliquent surtout aux étapes précoces du développement d'un AAA. Le modèle d'AAA est caractérisé par une hauteur maximale H et une largeur W , l'unité de mesure étant le diamètre d'entrée de l'artère. Pour com- mencer, la dynamique est étudiée pour les écoulements stationnaires. La stabilité globale de ces écoulements de base est analysée en calculant les valeurs propres et les fonctions propres pour des perturbations de faible amplitude. Pour comprendre les mécanismes d'instabilité, le transfer d'énergie entre l'écoulement de base et les perturbations est calculé. L'écoulement pour des AAA peu profonds (ou de grande longueur) se déstabilise par un mécanisme de 'lift-up' et les perturbations ampli- fiées sont stationnaires. Des anévrismes plus localisés (ou plus profonds) deviennent instables pour des nombres de Reynolds plus élevés, sans doute par instabilité el- liptique ; dans cette situation, les perturbations sont des modes oscillants. Dans le cas des écoulements pulsés, deux types de profil de débit physiologique ont été considérés dans cette étude, correspondant à une situation de repos ou d'exercice physique. Ces écoulements restent collés aux parois pendant la phase de systole et un écoulement décollé est généralement observé pendant la décélération après le maximum de systole. Dans cette phase, un vortex se forme à l'extrémité aval. Ce vortex s'agrandit au cours du temps et impacte l'extrémité aval de l'AAA, ce qui conduit à de forts gradients de contrainte pariétale, qui ne sont pas observés dans les cas sains. Il a été observé que les conditions d'écoulement varient significative- ment avec les nombre de Womersley (Wo) et de Reynolds (Re); l'écoulement reste attaché aux parois plus longtemps pour des nombres de Womersley croissants. Le principal effet d'une augmentation de Re est un renforcement du vortex primaire qui se forme après le maximum de systole. Les décollements de l'écoulement, l'impact de vortex au bord aval de l'AAA ou encore de faibles contraintes pariétales oscil- lantes (des caractéristiques importantes dans les cas d'anévrismes pathologiques) sont observés même pour des anévrismes de faible profondeur. Pour des anévrismes plus développés, des vortex multiples sont observés tout au long du cycle dans la cavité de l'AAA. Une analyse de stabilité de ces écoulements de base pulsés a mon- tré que le maximum des perturbations se développe vers l'extérmité aval des AAA. Cependant, les perturbations ne sont pas complètement confinées dans la cavité de l'AAA et se développent aussi au-delà en aval. On en déduit qu'une fois qu'un AAA s'est développé, les perturbations affectent aussi les artères saines en aval de l'AAA. Enfin, en considérant deux profils équivalents d'AAA, de formes sinusoï- dale et gaussienne, la sensibilité des résultats aux détails de la géométrie a pu être établie.

Ce travail de thèse, réalisé dans le cadre d'une convention CIFRE entre TURBOMECA et le CERFACS et en partenariat avec l'IVK, se place dans un contexte d'amélioration des performances des turbines axiales équipant les turboréacteurs d'hélicoptère. Un des points critiques du dimensionnement de tels moteurs est la maitrise de la durée de vie des pales de la turbine haute pression qui font face à de très hautes températures provenant de la chambre de combustion. Les prédictions numériques de l'environnement aérothermique des pales (écoulements dans la veine et système de refroidissement) sont réalisées aujourd'hui dans le milieu industriel à l'aide de la modélisation Reynolds Averaged Navier-Stokes (RANS). Grâce à des capacités de calculs grandissantes, l'approche Simulation aux Grandes Echelles (SGE) offre désormais un nouveau potentiel de prédictions d'écoulements. Les travaux de cette thèse s'intéressent ainsi à la capacité de la SGE à prédire l'écoulement du circuit de refroidissement interne d'une pale de turbine. Pour simplifier l'analyse de ce problème ou plusieurs phénomènes physiques sont en jeu, une progression en trois parties est proposée. La première s'intéresse à l'étude aérothermique de géométries simplifiées de canaux de refroidissement (coude à 180° et canal avec promoteurs de turbulence) en configuration statique. Aux régimes d'écoulement considérés, une approche résolue en paroi avec maillage non-structuré hybride est proposée et validée en vue d'une application industrielle facilitée. La seconde partie étend l'analyse de l'écoulement à un cas de canal avec promoteurs de turbulence en rotation utilisant une méthode de résolution numérique dans un repère absolu. Les investigations des résultats de la SGE fournissent des prédictions moyennes et instationnaires en bon accord avec les expériences disponibles et les travaux précédents aussi bien pour la dynamique de l'écoulement que les transferts de chaleur. Enfin, une troisième partie présente une application de la méthode sur un cas de pale réelle avec couplage thermique entre le circuit de refroidissement et le solide de la pale. Cette dernière partie classée confidentielle n'est pas présente dans le manuscrit disponible publiquement. Les résultats de l'approche résolue en paroi et de la rotation dans le repère absolu comparés aux résultats RANS disponibles pour le cas applicatif montrent d'importante différences locales et ainsi le potentiel de la méthode proposée.

Ce travail est une contribution pour mieux comprendre le transport de solutés dans un milieu poreux. Ce phénomène se rencontre dans de nombreux domaines: transport de contaminants dans les eaux souterraines, séquestration du CO2, stockage souterrain des déchets nucléaires, simulations de réservoirs pétroliers. On obtient la dispersion effective de Taylor en tenant compte de la convection, de la diffusion, de la géométrie du milieu poreux et des réactions chimiques. Le but de la théorie d'homogénéisation est, à partir d'équations microscopiques, de dériver un modèle effectif à l'échelle macroscopique. Ici, on applique la méthode de convergence à deux échelles avec dérive pour arriver au comportement effectif. Dans un premier temps, on considère les réactions de type adsorption à la surface des pores. À l'échelle microscopique, le phénomène de transport est modélisé par des équations couplées de type advection-diffusion, une pour la concentration dans le fluide et l'autre pour la concentration à la surface de milieu poreux. Le couplage est fait par les isothermes d'adsorption. Le système microscopique avec des coefficients fortement oscillants est étudié dans un régime de forte convection i.e., dans un régime de grand nombre de Péclet. La présence de forte convection dans le modèle microscopique se traduit par l'apparition d'une large dérive dans les profils de concentrations. On considère à la fois l'isotherme linéaire et l'isotherme non linéaire et les résultats ainsi obtenus sont comparés. Dans la deuxième partie, on généralise nos résultats concernant le transport réactif d'un seul soluté à ceux de plusieurs solutés dans un cadre linéaire. Dans ce cas, les paramètres effectifs sont obtenus en utilisant le principe de Factorisation et la convergence à deux échelles avec dérive. On étudie numériquement le comportement des paramètres effectifs par rapport au nombre de Péclet et au nombre de Damköhler. On utilise Freefem++ pour effectuer des tests numériques en dimension deux.

Cette thèse contribue à la simulation numérique des écoulements d'un fluide compressible modélisé par les équations de Euler et Navier-Stokes: étude d'un schéma d'ordre élévé basé sur une matrice de masse, modélisation des écoulement turbulents compressibles à très haut Reynolds, développement des estimateurs et correcteurs d'erreurs a posteriori et a priori, et adaptation de maillage anisotrope pour les fonctionnelles d'observation. Concernant la prédiction des écoulements turbulents, on s'est intéressé aux modèles hybrides de type RANS/LES comportant les nouveautés suivantes: traitement des tourbillons de grande échelle utilisant la formulation VMS (Variational Multi-Scale) et du RANS employé sur la paroi sur une distance imposée via une zone de protection conçue afin d'éviter le phénomène assez commun apelé "grid induced model depletion". Le niveau de viscosité du modèle VMS-LES est de plus controllé par un procédé de double filtre dynamique. La seconde partie concerne l'adaptation de maillage anisotrope pour mieux observer une fontionnelle d'observation. Les estimations a priori sont réalises pour le modèle des équations d'Euler et Navier-Stokes en instationnaire en 2D et 3D. A partir de ces estimations on sait définir les maillages optimaux au cours du calcul instationnaire, en fonction de l'état et de l'état adjoint. Le système d'optimalité est discrétisé et résolu à l'aide d'une méthode de point fixe instationnaire global, comportant une stratégie de stockage/recalcul pour le couplage état/ état adjoint. Des applications à la propagation d'ondes de choc et d'ondes acoustiques sont présentées.

L'objectif de cette étude est d'analyser l'écoulement traversant le récepteur solaire haute température de la centrale PEGASE. Cet écoulement est turbulent et soumis à un très fort gradient de température. Nous utilisons la simulation des grandes échelles thermiques appliquée aux équations bas-Mach pour étudier l'écoulement turbulent dans une géométrie de canal plan bi-périodique avec températures imposées aux parois. Pour cette étude, nous avons réalisé des simulations pour deux nombres de Reynolds turbulents différents (Re_tau m = 180 et Re _taum = 395) et pour quatre rapports de température (T2=T1 = 1 ; 1; 04 ; 2 et 5). Après avoir validé notre modèle, nous avons réalisé une étude sur la modélisation sous-maille thermique qui conclut sur la nécessité d'utiliser un modèle sous-maille thermique dynamique pour un écoulement ayant une forte intensité turbulente soumis à un très fort gradient de température. Nous analysons l'influence du gradient de température sur les différents profils de vitesse et de température ainsi que sur les spectres d'énergie. Augmenter le gradient de température crée une dissymétrie de tous les profils. Une relaminarisation du côté chaud est visible. On remarque une nouvelle répartition des fluctuations à travers le canal ainsi que de l'énergie en fonction de la taille des échelles. Une création de fluctuations de vitesse et de température a également été mise évidence du côté froid et du côté chaud du canal. Ces effets sont dus, à la fois à une variation de la viscosité cinématique (effet direct de la température) et à une interaction entre les champs turbulents dynamique et thermique.

Cette thèse a pour objet l'étude des couplages électrocinétiques en milieux non-saturés : d'une part les phénomènes d'électrofiltration liés à l'écoulement de l'eau, et d'autre part les couplages sismo-électriques, c'est à dire le champ électrique associé à la propagation d'une onde sismique sur le terrain. Ces couplages existent en raison du mouvement de particules chargées présentes dans l'eau. Le premier axe développé est une approche expérimentale. Des différences de potentiels d'électrofiltration ont été mesurées lors d'expériences de drainage réalisées dans une colonne de sable. Ces mesures combinées aux mesures des conditions hydrodynamiques ont mené aux premiers enregistrements continus de coefficients d'électrofiltration C en fonction de la saturation Sw. Les valeurs de coefficients normalisées ne varient pas de fa¸con monotone avec la saturation Sw. En effet, C augmente lorsque Sw diminue jusqu'à atteindre un maximum pour Sw égale à 65-85%, puis diminue avec la saturation. En outre, les valeurs expérimentales semblent pouvoir dépasser d'un facteur 100 la valeur mesurée à saturation. Ce comportement a été observé de façon similaire lors de trois expériences de drainage. La modélisation de ces signaux constitue le deuxième axe de travail. Une nouvelle expression, intégrant une contribution dynamique liée aux variations temporelles de pression, est proposée. Les résultats illustrent la nécessité d'intégrer cette contribution dynamique pour modéliser les signaux d'une des expériences, particulièrement lorsque C augmente pour Sw =[0. 65-0. 85 ; 1]. La dernière expérience tend à montrer que les conditions d'écoulement, et plus précisément le flux hydrique, pourraient jouer un rôle important dans le comportement de C. Le dernier axe développé est une approche de terrain. Les amplitudes de signaux sismo-électriques ont été étudiées en non-saturé, par la mesure d'un champ électrique Ex et d'un déplacement sismique ux, et comparées au calcul de fonctions de transfert Ex/üx. Une image de saturation du sol a été déduite de l'interprétation conjointe de données radar et de résistivité électrique. Les amplitudes sismo-électriques liées aux arrivées des ondes P et des ondes de surface, traitées indépendamment, ont pu être comparées à la teneur en eau, puis à des mesures de résistances électriques. Les résultats suggèrent que ces amplitudes sont plus faibles en milieu non-saturé. De plus, le traitement des ondes de volume et des ondes de surface a montré que ces dernières doivent être considérées indépendamment.

Dans ce mémoire, nous proposons une méthode permettant de modéliser, de façon globale, tout type de mouvement par des combinaisons linéaires de polynômes orthogonaux. Pour cela, nous projetons chaque champ de déplacement représentant le mouvement étudié sur une base orthogonale composée des polynômes de Legendre. Nous obtenons alors une expression polynomiale du mouvement considéré. Une modélisation d'ordre faible permet d'interpréter physiquement le mouvement dominant de séquences d'images. Cette caractéristique trouve un intérêt particulier dans des applications telles que la vidéo-surveillance, l'indexation ou l'étude du comportement. Nous montrons également que notre méthode permet de modéliser des mouvements complexes tels que des écoulements fluides. La représentation du mouvement sous forme polynomiale, les faibles temps de calcul ainsi que les taux de compression élevés sont alors des atouts importants pour traiter ce type de données.

Les écoulements dans les circuits d'admission et d'échappement des moteurs à combustion interne jouent un rôle majeur pour la détermination des différents rendements (volumétrique, de balayage et mécanique), la puissance indiquée, les performances, les émissions et en établissant le champ d'écoulement dans le cylindre. Lorsque les gaz s'écoulent transitoirement à travers ces systèmes, l'ensemble des forces de frottement, de pression et d'inertie de gaz sont présentes. L'importance relative de ces forces dépend de la vitesse du gaz et les dimensions et les formes de tels systèmes. Traditionnellement, ces écoulements sont étudiés au moyen des équations de dynamique des gaz unidimensionnelles (1D) où, le phénomène tridimensionnel de l'écoulement et la déformation des ondes de pression, la turbulence et la viscosité sont ignorés ou négligés. L'approche thermodynamique (0D) a été aussi utilisée où l'effet important d'inertie de fluide lié aux dimensions des composantes de ces circuits d'admission et d'échappement a été ignoré. Dans cette étude, la méthode de vidange remplissage, a été complètement révisée et une nouvelle méthode (Méthode Inertielle Capacitive MIC), basée sur la formulation thermodynamique de la méthode de vidange remplissage et sur l'équation fondamentale de la conservation de la quantité de mouvement, est développée. L'objectif est ensuite de tenir compte des effets d'inertie de fluide sur le comportement de fluide ou de l'écoulement sans utilisation du code unidimensionnel 1D (à cause du long temps de calcul). Pour cet objectif, des analyses de calcul CFD ont été élaborées afin de calculer les paramètres de calage de " la méthode inertielle capacitive " correspondant au nouveau modèle. Dans cette étude, il est apparu que les effets d'inertie ignorés de la formulation de la méthode de vidange remplissage sont les causes pour que cette dernière soit valable uniquement pour les circuits compacts d'admission et d'échappement de MCI. Pour valider le nouveau modèle, une investigation expérimentale a été réalisée sur un moteur monocylindrique à quatre temps. Le débit volumétrique du moteur est ensuite calculé avec le nouveau modèle. Le résultat est comparé à celui issu de l'expérimental et un bon accord a été obtenu.