La recherche expérimentale, présentée dans cette thèse, a pour objectif l'étude de l'influence d'un écoulement amont turbulent ou non sur les tourbillons coniques (primaire et secondaire) engendrés à partir d'une configuration de montant de baie. Cette configuration d'écoulement est rencontrée sur les véhicules terrestres et l'interaction des tourbillons avec les vitres latérales du véhicule est le principal responsable du bruit aérodynamique transmis dans le véhicule. Une maquette du type dièdre a été conçue et testée en soufflerie afin de générer et d'étudier ces tourbillons coniques. Les propriétés spatiales et temporelles des champs de vitesse et de pression pariétale ont été étudiées en associant simultanément la mesure de la vitesse par PIV ou Stéréo-PIV à la mesure de pression pariétale fluctuante par capteurs déportés. Nous avons montré qu'un niveau modéré de turbulence amont (intensité : 4%) se traduit par une modification très importante de la dynamique du tourbillon et de la pression pariétale fluctuante associée (Cp'). Notre conjecture - basée sur l'étude des moments statistiques et des corrélations spatio-temporelles - est qu'en l'absence de turbulence externe, la contribution principale au Cp' est liée à l'empreinte du coeur tourbillonnaire et du décollement secondaire induit. Au contraire, en présence de turbulence amont, la forte augmentation du Cp' semble associée à une réponse globale de la structure tourbillonnaire et à la modification associée de l'aérodynamique autour du corps. L'inclusion de dispositifs géométriques simples de contrôle passif sur la maquette a fait objet d'essais supplémentaires. Nous avons montré que ces dispositifs permettent de modifier très significativement le développement de la structure tourbillonnaire et son empreinte en terme de pression pariétale fluctuante.

Par l'exploitation de nos capacités visuelles, la visualisation scienti que entend proposer aux ingénieurs et aux chercheurs un outil visant à les assister dans l'acquisition de connaissance à partir de phénomènes complexes. Ce travail de thèse s'inscrit dans cette lignée, à ceci près qu'il s'intéresse plutôt à l'utilisation des technologies de la Réalité Virtuelle dans le but d'amener un utilisateur expert au coeur du processus d'exploration et d'analyse des données, ce que nous appelons Exploration de données scienti fiques. Dans le but d'arriver à des processus d'exploration efficaces, notre recherche s'est portée sur la mise en place de techniques d'interactions intuitives, susceptibles d'exploiter au mieux les capacités sensorimotrices de l'être humain. Pour atteindre cet objectif, deux conditions nous paraissent être essentielles. D'une part, il faut que les informations transmises via différents canaux sensorimoteurs aient une certaine cohérence à être délivrées ensemble : il est primordial que l'exploitation d'un canal pour véhiculer une information ne soit pas en concurrence avec ce qui est fait sur un autre canal. D'autre part, il est souhaitable que le potentiel de chaque canal soit utilisé au meilleur de sa capacité. Dans ce contexte, ce travail a débuté par une analyse de l'utilisation de l'haptique dans l'Exploration de données scientifi ques. Pour ce type d'usage, il a été identifié quatre tâches fondamentales (Sélectionner, Localiser, Relier et Arranger) pour lesquelles l'haptique semble présenter un réel avantage par rapport aux autres canaux sensoriels. Pour chacune de ces tâches, nous avons montré, au travers d'une large étude bibliographique, comment l'interaction haptique pouvait être exploitée afi n d'off rir des méthodes d'exploration efficaces. Sur la base de cette analyse organisée autour de ces quatre catégories, nous avons ensuite mis en évidence les problématiques liées aux tâches identifiées. Ainsi, nous avons souligné, d'une part que l'haptique pouvait faciliter la sélection de données scienti fiques dans des contextes où celles-ci sont massives, et d'autre part nous avons montré le besoin de mettre en place de nouvelles méthodes de suivi de structures d'intérêts (iso-surfaces, lignes de courant etc.). Notre problématique ayant ainsi été posée, nous avons d'une part étudié l'utilisation de retour multimodaux non visuels pour la recherche et la sélection de cibles dans un environnent virtuel 3d. Les situations impliquant une ou plusieurs cibles furent analysées, et plusieurs paradigmes d'interaction ont été proposés. Dans cet ordre d'idées, nous sommes arrivés à défi nir et valider un principe directeur pour l'usage de retours haptico-sonores pour la recherche et la sélection d'une cible donn ée située dans une scène 3d pouvant en contenir plusieurs autres. Nous avons en eff et montré que, pour une telle tâche, il était préférable d'exploiter la spatialisation sonore a n de localiser la cible désirée dans l'espace, tandis que le retour haptique permettait une sélection précise de la cible. D'autre part, nous nous sommes attaqués aux problèmes liés au rendu haptique d'ensembles de données pouvant présenter de fortes variations. A cet eff et, suite à un rappel de l'apport de l'haptique pour le rendu de surfaces dans le domaine médical, nous avons analysé certains besoins pouvant être comblés par l'ajout de ce canal sensorimoteur dans l'analyse d'iso-surfaces issues de simulation de Mécanique des Fluides Numérique (MFN). Par la suite nous avons proposé et évalué, par l'intermédiaire d'expériences de perception et de mesures de performance, de nouvelles méthodes de rendu haptique d'iso-surfaces dont l'une des originalités est de pouvoir se passer d'une représentation polygonale intermédiaire. En fin, nous avons appliquécette nouvelle approche d'exploration de données scientifi ques à l'analyse des résultats d'une simulation d'un écoulement dans une cavité ouverte. Ainsi, nous avons proposé deux méthodes d'analyse multi-sensorielle, dédiées à l'exploration d'un ensemble de données issu d'une simulation de MFN, en exploitant les approches génériques développées précédemment. La première méthode concerne une analyse interactive de la géométrie de l'écoulement, alors que la seconde se rapporte à une analyse multi-sensorielle de la topologie de l'écoulement. Les premières évaluations menées ont indiqué que les méthodes proposées tendaient à favoriser une meilleure compréhension du phénomène analysé et qu'elles pouvaient diminuer la charge cognitive habituellement requise par une telle tâche. A titre de conclusion, soulignons que cette thèse de doctorat ouvre la voie à un certain nombre de perspectives de recherches. A court terme, il s'agit de compléter les travaux relatifs à l'analyse de l'aspect dynamique de la simulation d'un écoulement dans une cavité ouverte. Ces travaux consisteront à proposer et valider différentes interactions multimodales visant à examiner les échanges/recirculations pouvant exister entre/dans les parties de la dite cavité. Mais au-delà des travaux relatifs à la MFN, l'expérience acquise à travers ces travaux pluridisciplinaires (informatique, mécanique des fluides, ergonomie cognitive) me permettra à moyen terme d'élargir mes travaux de recherche à l'exploration de données médicales, météorologiques, ou géologiques.

Dans l'étude de faisabilité d'un stockage géologique profond des déchets radioactifs, les questions essentielles posées aux ingénieurs et chercheurs en géomécanique sont la stabilité des ouvrages pendant la période de réversibilité et la variation des propriétés de confinement de la roche hôte induite par le creusement des cavités et par le stockage lui-même. Le présent travail porte sur l'étude des propriétés de confinement vis-à-vis des gaz des argilites du Callovo-Oxfordien du site de Meuse/Haute-Marne, étudié par l'ANDRA depuis 1995, au voisinage d'un puits d'accès. Un protocole de mesure de la perméabilité au gaz a été élaboré en prenant en compte la très faible perméabilité des argilites et les effets de désaturation. Les trajets de chargement mécanique et de pression de gaz ont été déterminés par rapport aux chargements in situ. Deux méthodes d'interprétation, une fondée sur une solution analytique simplifiée et l'autre sur une solution numérique, ont été utilisées pour déterminer la valeur de perméabilité au gaz. Des mesures de la perméabilité au gaz sur des éprouvettes prélevées à différentes distances, de 0,1 à 12,5 m, de la paroi du puits d'accès au Laboratoire Souterrain de Meuse/Haute- Marne ont été effectuées pour caractériser l'endommagement mécanique induit suite au fonçage du puits. Les résultats obtenus sous une contrainte isotrope de 11 MPa sont compris entre 10-21 et 10-22 m² et ils ne présentent pas de variations significatives entre la zone potentiellement endommagée (éprouvette proximité à la paroi) et la zone présumée intacte (éprouvette à 12,5 m de la paroi). Ce résultat est cohérent avec les caractérisations de l'endommagement menées in situ. Sous chargement déviatorique, la variation des mesures de perméabilité est inférieure d'un ordre de grandeur par rapport à celle mesurée sous contrainte isotrope ; elle est donc peu significative compte tenu de la précision du dispositif de mesure. Des éprouvettes jumelles d'argilites supposés intactes ont subi différents paliers hydriques à l'aide des solutions salines sursaturées (HR variant de 25 % à 98 %) pour constituer un cycle complet désaturation - resaturation. Pour chacune des valeurs du degré de saturation, des mesures de perméabilité ont été effectuées. La cinétique de désaturation, la courbe de rétention, les déformations hydriques ainsi que la fissuration hydrique des éprouvettes ont été mesurées et analysées. La valeur du coefficient de diffusivité hydrique déterminée d'après l'évolution de la masse des éprouvettes est estimée à 5.10-10 m².s-1. La perméabilité au gaz mesurée des éprouvettes, sous un confinement de 5 MPa, augmente de 10-21 m² à 10-18 m² lorsque l'humidité relative à laquelle les éprouvettes ont été conditionnées varie de 98 % à 25 %. La perméabilité au gaz (k) et le degré de saturation (Sr) de la roche présentent une relation quasi-linéaire dans un repère log(k) - Sr. L'évolution de la perméabilité au gaz en fonction d'un déviateur de contraintes a aussi été étudiée pour différents degrés de saturation. Les résultats confirment que l'effet du déviateur de contraintes sur la perméabilité n'est pas significatif, même si le déviateur de contraintes dépasse le seuil d'endommagement et s'il est proche du seuil de rupture. Ces résultats ont été comparés à ceux obtenus par d'autres laboratoires utilisant des techniques différentes (Zhang 2007). La déformabilité et la vitesse des ondes des éprouvettes d'argilites ont été suivies pendant les essais. La microstructure des argilites a été étudiée à l'aide de la porosimétrie au mercure. Les facteurs influençant la mesure de perméabilité tels que : l'effet Klinkenberg, la variation du degré de saturation pendant l'essai, la dissolution du gaz et le déplacement d'eau par le gaz, et surtout l'étanchéité du système, ont été abordés dans cette étude.

La connaissance du comportement thermique des fours et des pièces est un problème difficile et cependant essentiel dans les thématiques de recherche industrielles actuelles. L'industrie cherche à se doter de moyens numériques de plus en plus efficaces tout en réduisant sans cesse le temps de calcul afin de modéliser des pièces et des assemblages de plus en plus réalistes. Les cinétiques de chauffe ainsi que la distribution de la température dans l'enceinte et dans les pièces traitées doivent également être maîtrisées, afin d'améliorer la qualité des produits chauffées en terme de résistance et dureté. Le travail présenté dans cette thèse porte donc sur le développement de méthodes de résolution numérique pour la simulation du transfert thermique. Ces méthodes permettent le calcul couplé de la température des pièces avec l'environnement du four afin d'optimiser la géométrie de ces pièces, leur position dans le four mais également les conditions initiales de chauffe. Le calcul de dynamique des fluides (CFD) peut être considéré comme un puissant outil technique de prédiction de l'écoulement des fluides et du transfert thermique dans des cas industriels réalistes. La première partie de la thèse porte sur l'étude de schémas numériques avancés relatifs aux méthodes éléments finis stabilisés. Ces méthodes de calculs efficaces ont étaient utilisées pour simuler des écoulements instationnaires et des transferts thermiques conjugués. Par la suite, deux modèles de turbulence (modèle k-epsilon et modèle Large Eddy Simulation) sont introduits et utilisés pour prendre en compte les écoulements complexes et turbulents dans l'enceinte du four. Le transfert radiatif est assuré par la résolution du modèle P1 tout en calculant un terme source volumique qui sera intégré dans l'équation de la chaleur. La deuxième partie de cette thèse se consacre à la mise en place d'une méthodologie de discrétisation robuste qui permet aux utilisateurs de générer de façon entièrement automatique un seul maillage. Ce maillage contient à la fois des domaines axés sur la résolution d'un problème fluide (air, eau, ...) mais également solides spécifiques aux structures, et ce quel que soit le niveau de détail et donc de complexité du cas étudié. Cette approche, connue sous le nom d' « immersion de volume », garantit un maillage anisotrope précis aux interfaces fluide-solide afin de capturer plus précisément les gradients thermiques et la forte discontinuité des propriétés physiques. Cette méthode offre donc une grande flexibilité dans la mise en données du problème mixte fluide-structure et aussi dans la prise en compte de plusieurs géométries (four, pièces, supports) et elle est également bien adaptée aux solveurs thermomécaniques développés.

A bas nombre de Reynolds, le concept de voilure battante est récemment apparu comme une alternative aux concepts conventionnels de voilure fixe et voilure tournante. Dans le cadre d'une application pratique (micro-drones ou MAVs), l'évaluation de l'adaptabilité d'un tel mode de sustentation requiert la compréhension fine des principaux mécanismes aérodynamiques mis en jeu et de leur impact sur les efforts résultants. Ces derniers se caractérisent par une instationnarité forte et des comportements complexes. Le travail réalisé s'articule autour de trois points. Dans un premier temps, l'établissement d'une corrélation précise entre comportement spatio-temporel des structures cohérentes et effort résultant passe par la définition de méthodes d'investigation adaptées à l'étude d'écoulements bi et tridimensionnels résolus en temps. En ce sens, l'analyse repose d'une part sur une approche numérique DNS par maillage mobile, et d'autre part sur les approches expérimentales TRPIV- 2D2C/3D3C conduisant à la mesure non-intrusive des efforts aérodynamiques par application du bilan de quantité de mouvement. Les difficultés associées à ces techniques novatrices sont ici mises en exergue. En second lieu, nous abordons deux types de configurations bidimensionnelles caractérisées par des cinématiques de battement respectivement symétriques et asymétriques. Nous nous limitons à l'étude du vol stationnaire ou sur place pour lequel les effets instationnaires sont prépondérants vis- à-vis des effets quasi-stationnaires. D'abord, l'analyse détaillée d'un cas de référence permet d'affiner les résultats reportés dans la littérature. Ensuite, la confrontation systématique de cas asymétriques rév`ele l'occurrence de phénom`enes aérodynamiques particuliers contribuant favorablement à l'efficacité de l'aile. L'approximation bidimensionnelle occulte toutefois l'influence du tourbillon marginal ou du gradient transverse de vitesse incidente établi le long de l'envergure d'une aile en révolution autour de son emplanture. L'une des spécificités des configurations tridimensionnelles provient effectivement du caract`ere stabilisateur de certains mécanismes aérodynamiques. Ainsi, nous présentons les effets tridimensionnels et leur dépendance vis-à-vis de l'allongement dans le cas de mouvements de battement en translation rectiligne et en rotation autour d'un axe vertical.

Un écoulement de coextrusion peut être considéré comme un écoulement de Poiseuille multicouche de fluides viscoélastiques. Nous avons étudié expérimentalement l'influence des conditions opératoires sur la stabilité de l'écoulement polyéthylène/polystyrène. Nous avons ensuite développé une analyse de stabilité linéaire aux ondes longues pour une loi de comportement de White-Metzner. Cette analyse a donné des résultats cohérents avec l'expérience dans certaines situations et des résultats mitigés voire contradictoires dans d'autres. Nous avons alors montré que l'instabilité était de nature convective. Nous avons utilisé à cet effet deux montages de coextrusion originaux: une filière d'essais industriels de très grande longueur démontable et une filière de laboratoire transparente qui nous permet d'observer en temps réel le développement de l'instabilité le long de l'écoulement. Nous avons alors développé deux outils numériques pour rendre compte de ces instabilités convectives: une modélisation directe, utilisant pour les deux polymères une loi de comportement de Maxwell et une étude de stabilité linéaire spatiale pour une loi de comportement de White-Metzner. Cette dernière étude nous a permis de rendre compte de l'essentiel des phénomènes observés. Nous avons alors montré qu'une telle étude permet de déterminer, pour une couple de produits donnés, l'incidence des conditions opératoires sur l'apparition d'instabilités interfaciales.

La compréhension des écoulements turbulents en aval d'un cylindre circulaire horizontal, soumis à un écoulement d'air transverse et en présence d'importants effets de flottabilité, constitue l'objectif principal de ce présent mémoire. Les approches expérimentale et numérique (DNS) sont utilisées afin de compléter les connaissances des mécanismes d'instabilité tourbillonnaire en régime de convection mixte. En dehors de la perte de symétrie du sillage induit par les effets de pesanteur, l'écoulement étudié ici est caractérisé par le développement de deux instabilités principales. Au niveau de la couche de cisaillement inférieure on assiste au développement d'instabilités intermittentes de Kelvin-Helmholtz. Nous montrons qu'en présence d'effets de gravité importants la transition vers la turbulence est accélérée dans cette région de l'écoulement. A l'arrière du cylindre un écoulement secondaire apparaît sous l'effet des forces de pesanteur. Cet écoulement est le lieu d'importants mouvements tridimensionnels tant sur le plan moyen que fluctuant et donne naissance à des dipôles de vorticité caractéristiques. Ces structures se forment dans la partie supérieure du sillage et évoluent librement au niveau de l'écoulement principal ou s'apparient avec les structures tourbillonnaires du sillage turbulent. Finalement, en présence d'effets de gravité non négligeables, on assiste à un enrichissement du spectre des structures turbulentes et à une augmentation des composantes tridimensionnelles du champ dynamique. Dans ces conditions la température agit comme un scalaire actif capable d'accélérer la transition vers la turbulence.

Un lit de particules, par exemple le lit d'une rivière ou les dépots d'hy drates dans un pipeline, peut évoluer de différentes manières quand il est soumis à un écoulement de fluide. Nous étudions expérimentalement l'interaction fluide/particules dans la géométrie d'un tube. En faisant varier les paramètres expérimentaux, nous avons déterminé l'existence de cinq régimes différents : particules immobiles, mise en mouvement du lit plat, formation de dunes laminaires, formation de dunes à vortex et apparition de dunes sinueuses. Ce travail s'articule principalement en quatre parties. Après l'étude qualitative des différents régimes existant, nous avons montré que le lit de particules se met en mouvement quand les forces hydrodynamiques deviennent supérieures à une fraction du poids apparent des particules, c'est-à-dire pour un nombre de Shields critique constant sur une large gamme de Reynolds particulaire. Un modèle continu à deux phases, dans lequel nous avons utilisé une rhéologie de type frottement solide pour modéliser la contrainte solide, nous permet de déterminer le flux de particules et de prédire l'évolution expérimentale du lit plat. En réalisant une étude de stabilité linéaire dans laquelle le flux de particules obtenu est couplé avec l'écoulement de fluide sur un fond fixe, nous avons déterminé le seuil d'apparition des dunes laminaires. Ce seuil correspond à un nombre de Reynolds constant et présente un bon accord avec les résultats expérimentaux dans la limite du modèle. Les expériences et les modèles réalisés nous permettent de construire un diagramme de phase représentant les zones d'existence des différents régimes observés dans le tube.

Le but de cette thèse est de développer une méthode de volumes finis qui s'applique à une classe de maillages beaucoup plus grande que celle des méthodes classiques, limitées par des conditions d'orthogonalité très restrictives. On construit des opérateurs différentiels discrets agissant sur les trois maillages décalés nécessaires à la construction de la méthode. Ces opérateurs vérifient des propriétés discrètes analogues à celles des opérateurs continus. La méthode est tout d'abord appliquée au problème divergence-rotationnel qui peut etre considéré comme une brique du problème de Stokes. Ensuite, le problème de Stokes est discrétisé avec diverses conditions aux limites. Par ailleurs, il est bien connu que lorsque le domaine est polygonal et non-convexe, l'ordre de convergence des méthodes numériques se dégrade. Par conséquent, nous avons étudié sous quelles conditions un raffinement local approprié permet de restaurer l'ordre de convergence optimal. Enfin, nous avons discrétisé le problème non-linéaire de Navier-Stokes, en utilisant la formulation rotationnelle du terme de convection, associée à la pression de Bernoulli. Par un algorithme itératif, nous sommes amenés à résoudre un problème de point-selle à chaque itération, pour lequel nous testons quelques préconditionneurs issus des éléments finis, que l'on adapte (quand c'est possible) à la méthode. Chaque problème est illustré par des cas tests numériques sur des maillages "arbitraires", tels que des maillages fortement non-conformes.

Nous présentons dans cette étude le développement de la convection à partir de différentes perturbations de l'état conductif d'une couche fluide confinée dans une cavité cylindrique, chauffée par le bas et avec une surface supérieure libre. La discrétisation spatiale du domaine repose sur la méthode des éléments spectraux et les itérations temporelles sont assurées par une méthode splitting.Au déclenchement de la convection, les structures convectives correspondent à des modes de Fourier, et les seuils critiques dépendent du rapport de forme de la cavité, et des nombres de Biotet de Marangoni qui caractérisent la surface libre. Les transitions d'écoulements au-delà du seuil primaire sont caractérisées quantitativement en fonction du nombre de Rayleigh pour différentes valeurs du nombre de Biot et Ma = 0. Les résultats présentés sont obtenus en résolvant l'ensemble des équations non-linéaires de conservation à travers une méthode de continuation. Lorsque la convection se déclenche sous la forme d'un mode axisymétrique m = 0, l'évolution non-linéaire montre la coexistence de différentes structures convectives, des structures axisymétriques avec écoulement montant ou descendant au centre de la cavité et des structures correspondant à des combinaisons de modes qui apparaissent sur des branches secondaires sous-critiques.L'action d'un champ magnétique constant est ensuite étudiée pour des fluides conducteurs dans une même configuration comprenant une surface supérieure libre. Nous montrons l'effet stabilisateur du champ magnétique sur les seuils primaires ainsi que son action sélective sur les différents modes de convection. Nous analysons l'évolution des structures convectives au delà de ces seuils et montrons comment le champ magnétique modifie les transitions entre ces structures.En soumettant le bain fondu à un champ magnétique tournant, le mouvement de rotation du fluide se superpose aux mouvements de convection thermique et on observe une diminution des fluctuations de température et un retard du déclenchement de l'instabilité de Rayleigh-Bénard(lorsque les deux parois haut/bas du bain sont rigides). La rotation influe sur ce déclenchement qui de stationnaire devient oscillatoire, à l'exception du mode m = 0 de Fourier, pour qui la transition reste stationnaire jusqu'à une certaine valeur critique du nombre de Taylor magnétique.La dynamique de l'écoulement axisymétrique de part et d'autre de cette valeur critique sera étudiée en détail.