Cette travail de thèse est lié au besoin de modélisation des écoulements diphasiques en générateurs de vapeur (entrée liquide et sortie vapeur). La démarche proposée consiste à faire le choix d'une modélisation hybride de l'écoulement, en scindant la phase gaz en deux champs, modélisés de manières différentes. Ainsi, les petites bulles sphériques sont modélisées avec une approche dispersée classique avec le modèle eulérien à deux fluides, et les bulles déformées sont simulées à l'aide d'une méthode de localisation d'interface. Le travail effectué porte sur la mise en place, la vérification et la validation du modèle dédié aux larges bulles déformées, ainsi que le couplage entre les deux approches pour le gaz gaz, permettant des premiers calculs de démonstration utilisant l'approche hybride complète.

Les vibrations induites par vortex (VIV) d'un solide immergé interagissant fortement avec son sillage sont des oscillations auto-entretenues. Les VIV peuvent atteindre des amplitudes significatives lorsque la fréquence du détachement tourbillonnaire est proche d'une des fréquences propres de la dynamique du solide, et constituent donc un mécanisme d'interactions fluide/structure intéressant pour extraire de l'énergie d'un écoulement. Les VIV de structures flexibles sont par conséquent analysées dans cette thèse à la fois numériquement, grâce au concept d'oscillateur de sillage, et expérimentalement. Nous démontrons ainsi qu'il est possible d'extraire efficacement de l'énergie d'un écoulement grâce aux VIV d'un câble. Nous prouvons également que la flexibilité de la structure lui permet d'adapter sa dynamique aux sollicitations fluctuantes du sillage pour extraire efficacement de l'énergie de l'écoulement, grâce à (i) l'excitation de ses différents modes de vibrations, dans le cas d'un câble tendu équipé de récupérateurs d'énergie périodiquement distribués, chapitre 2, ou (ii) l'apparition d'ondes progressives de VIV transportant l'énergie lorsque l'on considère un câble suspendu par un unique récupérateur d'énergie, chapitre 3.

Ce travail s'inscrit dans une démarche visant à quantifier l'impact de la diagenèse sur les ratios élémentaires Sr/Ca, Mg/Ca et U/Ca caractéristiques de l'exosquelette (aragonite biogénique) de coraux Porites sp. Ces ratios élémentaires sont très utilisés comme paléo thermomètres en paléo climatologie tropicale pour reconstruire les paléo variables environnementales comme la température de surface de l'océan. On considère dans un premier temps un corail moderne prélevé in vivo et un corail fossile daté au 14C de l'Holocène Moyen (5445 ans BP). Les deux coraux ne présentent aucune trace d'altération diagénétique. On mesure Sr/Ca, Mg/Ca, U/Ca dans l'exosquelette par spectrométrie de masse couplée à un plasma inductif et à la microsonde de Castaing pour Sr/Ca. Les températures de l'océan du temps du vivant des coraux sont alors reconstruites à partir de la géochimie des échantillons et validées. L'approche est ensuite étendue à des coraux pléistocènes (~125000 ans BP) altérés par la diagénèse e.g., apparition de calcite au détriment de l'aragonite initiale sous l'action de l'eau douce percolant dans les récifs. Cette calcitisation affecte les ratios Sr/Ca et Mg/Ca originels et est susceptible d'entacher les reconstructions d'un biais qu'il est important de bien connaître. En s'appuyant sur la spectrométrie Raman et la microanalyse X, on confirme les tendances qualitatives rapportées dans la littérature sur l'impact de la calcite sur les proxies du climat : la présence de calcite diminue (resp. rehausse) le ratio Sr/Ca (resp. Mg/Ca) créant donc des artéfacts chauds dans les températures reconstruites.

L'éolienne Darrieus connaît un intérêt accru ces dernières années parce qu'elle représente une solution alternative potentielle de production d'électricité dans les milieux urbains. En particulier, une éolienne de forte solidité peut être choisie car certaines de ses propriétés peuvent être avantageuses pour son implantation proche de zones habitées. A l'inverse, certaines difficultés aérodynamiques émergent. Ce type d'éolienne fonctionne à de faibles vitesses réduites pour lesquelles le décrochage dynamique a un rôle très significatif. L'objectif de ce travail de thèse consiste à compléter la connaissance du phénomène de décrochage dynamique sur une éolienne à axe vertical afin d'améliorer les modèles numériques de prédiction existants. Cette étude s'appuie sur une analyse combinée de résultats numériques et expérimentaux. Les simulations numériques sont produites avec une méthode des panneaux bidimensionnelle instationnaire. Les effets de la viscosité sont introduits par des corrections utilisant notamment un modèle semi-empirique de décrochage dynamique. Le travail expérimental s'est concentré sur la dynamique tourbillonnaire à proximité immédiate du rotor résultante du décrochage dynamique. Le montage se compose d'une éolienne à pale droite placée dans une soufflerie. Des mesures instationnaires de la répartition de pression pariétale le long de la corde et des mesures de champ de vitesse par vélocimétrie par images de particules ont été accomplies. Les résultats révèlent la manière dont les caractéristiques du décrochage dynamique sont conditionnées par la vitesse réduite. Le retard au décrochage, l'intensité de l'effet du tourbillon de décrochage dynamique et sa convection ont été quantifiés. Enfin, un examen critique de l'applicabilité du modèle de Leishman-Beddoes pour simuler efficacement les effets du décrochage dynamique a été réalisé.

Le développement de techniques de mesure optiques de surfaces libres permettant la mesure des ondes générées par un modèle de navire en bassin des carènes est proposé dans ce mémoire. Trois méthodes, toutes basées sur un système de prises de vues stéréoscopiques, sont présentées. La première repose sur l'emploi de procédures d'inter-corrélation des images des caméras composant le système de stéréovision. La seconde s'appuie sur la théorie de la réfraction de la lumière au travers de l'interface air/eau. La troisième méthode est basée sur l'identification des projections de faisceaux laser sur la surface libre vue par un système stéréoscopique. Ces méthodes sont appliquées à la mesure du champ de vague autour d'un modèle de navire tracté dans le bassin des carène de l'Institut Pprime et dans un bassin des carènes de grande taille (B600 de la DGA Techniques hydrodynamique) et permettent la reconstruction et l'analyse du sillage. Les principales caractéristiques des champs de vagues et leurs dépendances au nombre de Froude sont étudiées en détail. Afin de compléter l'étude de l'écoulement généré par le modèle de navire, des mesures de PIV stéréoscopique autour de la carène, des mesures spécifiques de la vague d'étrave et des simulations numériques sont réalisées. L'ensemble de ces données permettent l'analyse de champs de vitesse tridimensionnelle proche de la coque et l'identification des caractéristiques du champ de vague de l'étrave jusqu'au champ lointain.

La dispersion dans des milieux poreux homogènes (empilements de grains) a été étudiée par des mesures par résonance magnétique nucléaire (RMN) et des simulations de marches aléatoires dans un réseau de pores. La RMN permet de mesurer l'ensemble des déplacements des molécules d'eau durant un temps tΔ, et d'obtenir propagateurs et moments caractéristiques. L'évolution temporelle du second moment σ (σ2 ∝ taΔ) permet de caractériser de manière précise le régime de dispersion des molécules (Gaussien ou anormal). Des mesures pour des écoulements de 15 < Pe < 45 dans un empilement de grains de 30μm ont permis d'observer une dispersion anormale faiblement super-dispersive (a = 1.17) en écoulement saturé et une augmentation progressive du caractère super-dispersif avec la diminution de la saturation en eau (jusqu'à a = 1.5 pour 42 %)lors d'une co-injection stationnaire eau-huile. En écoulement saturé, les propagateurs et courbes de percée sont quasi-gaussiennes, tandis qu'en écoulement insaturé, les propagateurs sont asymétriques et les courbes de percée présentent des trainées aux grands temps. Dans ces conditions, on montre que la dispersion anormale observée est mieux décrite par des lois stables de Lévy que par des lois gaussiennes. Des simulations de marche aléatoire ont été réalisées dans un réseau de pores extrait d'un milieu poreux réel par imagerie microscanner.Elles permettent d'obtenir les mêmes informations que la RMN, les marcheurs se déplaçant par advection et diffusion. Ces simulations montrent l'existence d'une stagnation non observée dans les expériences, montrant que la simplification du réseau poreux est trop importante et empêche de reproduire certains aspects du champ de vitesses détecté par la RMN. Toutefois, l'évolution temporelle du second moment a également un caractère super-dispersif à temps long à 100 % de saturation

Une étude expérimentale des sillages turbulents autour de géométries tridimensionnelles est conduite en augmentant leur complexité depuis les corps axisymétriques jusqu'aux véhicules automobiles. Quelle que soit la géométrie, deux types de mouvement cohérent sont susceptibles d'être observés dans le sillage. Tout d'abord, sur des échelles de temps de l'ordre de 5D/U, D et U étant respectivement la taille caractéristique du corps et la vitesse de l'écoulement, le sillage peut présenter des oscillations périodiques. Elles sont associées aux interactions entre deux couches cisaillées, se faisant face avec une vorticité opposée. Au premier ordre, la fréquence de ces oscillations dépend de la distance séparant les deux couches de mélange ; par conséquent, deux fréquences différentes peuvent être trouvées lorsque le corps a un rapport d'aspect différent de 1. Ces modes globaux oscillants semblent s'atténuer à mesure que le nombre de Reynolds et la complexité géométrique du corps augmentent. Le second type de mouvement cohérent est associé au développement d'instabilités stationnaires. Elles sont reliées aux brisures de symétrie observées en régime laminaire et leur domaine d'apparence est définit à partir de considérations géométriques dans le cas de corps parallélépipédiques avec effet de sol. Ces instabilités génèrent de fortes asymétries dans la topologie instantanée du sillage et notamment peuvent produire des écoulements bistables avec un temps caractéristique supérieur à 100H/U. L'étude de ces phénomènes, associée à des analyses de sensibilité à de petites perturbations, montre que la diminution des asymétries dans le sillage instantané est une stratégie intéressante de réduction de traînée. En particulier, il est montré que les gradients de pression à échelles locale et globale sur les cotés du corps sont sources de structures tourbillonnaires longitudinales qui pénalisent la traînée. Les dépendances quadratiques entre la traînée et les efforts transverses obtenues rappellent les mécanismes de traînée induite bien connus en aéronautique. Ainsi, comme elles sont souvent attribuées à d'importantes asymétries du sillage, les instabilités stationnaires sont identifiées comme des sources importantes de traînée. Par contre, la partie de la traînée associée aux mouvements périodiques du sillage semble négligeable, surtout pour les géométries complexes à grands nombres de Reynolds.

La Simulation aux Grandes Echelles (SGE) est de plus en plus utilisée dans les processus de développement et la conception des réacteurs aéronautiques industriels. L'une des raisons pour ce besoin résulte dans la capacité de la SGE à fournir des informations instantanées d'un écoulement turbulent augmentant la quantité des prédictions de la composition des gaz d'échappement. Ce manuscrit de thèse aborde deux sujets récurrents de la SGE. D'une part, les schémas numériques pour la SGE nécessitent certaines propriétés, notamment une précision élevée avec une diffusivité faible pour ne pas nuire aux modèles de turbulence. Afin de répondre à ce pré requis, une famille de schémas d'intégration temporelle d'ordre élevée est proposée, permettant de modifier la diffusion numérique du schéma. D'autre part, la SGE étant intrinsèquement instationnaire, elle est très consommatrice en temps CPU. De plus, une géométrie complexe prend beaucoup de temps de simulation même avec les super calculateurs d'aujourd'hui. Dans le cas particulier d'intérêt et souvent rencontré dans les applications industrielles, l'approche bas-Mach est constitue une alternative intéressante permettant de réduire le coût et le temps de retour d'une simulation LES. L'impact et la comparaison des formalismes compressible et incompressible sont toutefois rarement quantifiés, ce qui est proposé dans ce travail pour une configuration représentative d'un brûleur swirlé industriel mesuré au CORIA.

Dans le cadre de l'amélioration des méthodes de modélisation du comportement hydraulique des moteurs à ergols liquides, il est important de disposer des méthodes d'identification expérimentale adaptées au besoin de caractérisation des fonctions de transfert des pompes en régime cavitant, notamment pour prédire les risques d'occurrence des phénomènes de couplage fluide-structure comme le Pogo. La matrice de transfert dynamique d'une pompe relie dans le domaine fréquentiel les fluctuations complexes (module + phase) de débit et de pression à l'entrée avec celles de sortie de la pompe. Par ailleurs, les mécanismes physiques qui régissent la réponse d'une pompe soumise à des fluctuations de pression et/ou de débit restent méconnus, ce qui rend aujourd'hui très difficile les prévisions théoriques. Dans le cadre des travaux développés, une nouvelle boucle d'essais a été mise en place permettant la mesure de la pression fluctuante et l'évaluation du débit fluctuant par une approche d'intensimétrie hydroacoustique. Une méthodologie d'essais a été mise en place sur le banc de test pour caractériser les fonctions de transfert de différents organes hydrauliques, notamment en régime cavitant, comme une pompe centrifuge. Ces travaux expérimentaux ont permis également d'alimenter et de valider les modèles de type " système " RLC utilisés pour la simulation du comportement acoustique du circuit. Les méthodologies d'essais, de post-traitement et de modélisation développés dans le cadre de cette thèse seront appliquées par l'industriel lors des essais dédiés sur le banc de test à des composants moteur.

La pollution de l'environnement par les pesticides reste une préoccupation sociale et environnementale importante. Lors de l'application, une partie des pesticides peut contaminer l'environnement (dérive). La diminution des pollutions repose sur la maîtrise de la taille et la vitesse des gouttes en sortie de buse. L'objectif principal de cette thèse est de développer une approche Eulérienne afin d'estimer les caractéristiques initiales des gouttes produites telles que la taille et la vitesse à la sortie de buse. Le modèle Eulérien d'atomisation considère l'écoulement diphasique d'un liquide et d'un gaz comme un écoulement turbulent d'un seul fluide avec une masse volumique variable, variant entre celle du gaz et celle du liquide pulvérisé. Une équation de transport pour la fraction massique liquide moyenne permet de décrire la dispersion du liquide dans la phase gazeuse. La turbulence est modélisée par une approche aux tensions de Reynolds (RSM) en résolvant les équations de transport de chacune des six composantes du tenseur de Reynolds. Par ailleurs, une équation de transport de l'interface volumique liquide/gaz est considérée. La production de l'interface volumique qui exprime la création des gouttes est fonction, d'une part, à grande échelle, du gradient de vitesse moyenne et, d'autre part, à petite échelle, de la turbulence. La destruction de l'interface volumique qui exprime la coalescence des gouttes, quant à elle, prend en compte la tension de surface qui s'oppose à la désintégration de la surface liquide. L'évolution spatiale des rayons et des vitesses des gouttes produites en sortie de buse est décrite en couplant l'équation de la fraction massique liquide moyenne et celle de l'interface volumique à celles de la conservation de la masse, de la quantité de mouvement et de la turbulence. L'étude est faite avec le code de calculs FLUENT V.12 en utilisant les Fonctions Définies par l'Utilisateur (UDF "User Defined Function") pour adapter le code à nos conditions d'injection très particulières puisque le rapport des masses volumiques liquide et gaz est de l'ordre de 800. Les résultats obtenus par le modèle indiquent la formation d'une nappe conique creuse constituée de grosses gouttelettes et la présence d'une zone de recirculation près de l'axe du spray constituée de gouttelettes plus petites, conformément aux expérimentations. La comparaison du Diamètre Moyen de Sauter calculé par le modèle et mesuré expérimentalement par l'Anémométrie Phase Doppler montre un bon accord.