Les matériaux granulaires possèdent des propriétés rhéologiques et mécaniques peu communes. Ils peuvent en effet s'écouler comme des liquides mais, sous certaines conditions, se bloquer et résister à des contraintes extérieures sans se déformer. Dans ce mémoire, nous étudions cette dualité solide/liquide à travers deux expériences dont la caractéristique majeure est de permettre la comparaison des comportements individuels des grains à ceux de l'empilement. La première partie du mémoire est consacrée aux écoulements à la surface d'un empilement incliné. Nous dérivons tout d'abord une description continue à partir des équations de conservations intégrées sur l'épaisseur en mouvement en étudiant les écoulements stationnaires dans un tambour tournant. Nous nous intéressons ensuite à la rhéologie du matériau et analysons les corrélations spatiales du champ de vitesse et les fluctuations de compacité. Nous montrons ainsi l'existence d'agrégats de grains trop denses pour pouvoir se déformer. La taille de ces agrégats est distribuée en loi de puissance, sans échelle caractéristique. Nous en discutons les implications sur les approches théoriques proposées pour décrire les écoulements denses. La deuxième partie porte sur les réarrangements internes dans un empilement apparemment statique soumis à une petite sollicitation. Une perturbation thermique de quelques degrés engendre des fluctuations électriques géantes dans un empilement de grains conducteurs. L'analyse statistique de ces fluctuations nous permet de relier cette sensibilité électrique à des fluctuations tribologiques locales et non pas à des réorganisations collectives des contacts comme cela a pu être suggéré par le passé. Nous visualisons dans un deuxième temps les micro-déplacements des grains engendrés par l'ajout d'une petite surcharge en surface. Cela nous permet de discuter la validité des différents modèles décrivant la répartition des contraintes dans un empilement statique.

Cette thèse est dediée à l'étude des mécanismes d'instabilité non-modaux dans les écoulements cisaillés, principalement des couches de mélange. On se concentre sur les perturbations linéarisées qui ont la plus grande croissance d'énergie à un temps donné, les 'perturbations optimales', différentes du mode propre le plus instable pour les systémes non-normaux. Une description originale de la non-normalité et ses conséquences est donnée dans le chapitre 2. Le chapitre 3 traite du mécanisme de 'lift-up' dans le cas des perturbations longitudinales non-visqueuses sur un écoulement de base parallèle et arbitraire. On trouve une nouvelle équation 1D qui détermine l'ensemble des perturbations orthogonales, dont l'optimale, et permet de trouver des nouveaux résultats exacts et asymptotiques. Dans le chapitre 4 on s'intéresse aux instabilités secondaires d'une couche de mélange en utilisant comme état de base l'écoulement 2D instationnaire et non-linéaire pour calculer les perturbations optimales 3D. Selon le nombre d'onde et les temps d'optimisation (initial et final), on retrouve comme réponses optimales les perturbations de types 'elliptique' et 'hyperbolique'. Dans les chapitres 5 et 6 on considère une stratification en densité dans la direction orthogonale au plan de l'état de base, et on utilise une décomposition de 'Craya-Herring' pour analyser les perturbations en termes de leur contenu en ondes internes et vorticité verticale. Les équations d'évolution des perturbations linéarisées autour d'un état de base 2D général sont obtenues et analysées au niveau de l'énergie onde/tourbillon. Ces résultats sont appliqués dans le chapitre 6 pour analyser la génération et l'émission d'ondes dans les perturbations optimales sur une couche de cisaillement horizontale.

La tomographie de résistivité électrique (ERT) est une méthode courante géophysique de terrain, souvent utilisée pour détecter et l'évolution suivre les panaches de polluants en zone saturée. L'ERT est cependant une méthode intégratrice dont la fiabilité des modèles est confrontée aux problèmes de non-unicité des solutions du problème inverse. Ces contraintes limitent l'interprétation des modèles à un aspect qualitatif de la distribution des contrastes de résistivité modélisés en 2D ou 3D, résultant du choix des paramètres d'inversion et de l'association de paramètres du milieu non identifiables à l'échelle du volume poreux. Cette thèse propose de tester la faisabilité de la méthode pour quantifier les paramètres de transport de polluants et de solutés miscibles au contact des eaux souterraines, ainsi que la sensibilité des paramètres d'inversion les plus influents sur la modélisation. Les tests expérimentaux sont réalisés en laboratoire sur des empillements 2D de billes de verre sphériques (de l'ordre de la 100ène de μm) dans un réservoir en plexiglas transparent. Deux réseaux verticaux de 21 électrodes sont disposés sur les bords latéraux du réservoir pour effectuer le suivi ERT du traceur salé (NaCl dissout) à partir de 210 points de mesures en dipôle-dipôle transverse acquis toutes les 5 minutes afin d'optimiser la résolution temporelle. Le dispositif est également disposé face à un panneau lumineux permettant de réalisé un suivi vidéo simultané du colorant. L'analyse vidéo révèle une propagation plus rapide du colorant sur les bords latéraux qui reste somme toute négligeable pour les débits à la pompe les plus faibles. En revanche les mesures ERT sont fortement perturbées par les effets résistant de la cellule plexiglas qui se répercutent sur les modèles. La normalisation des mesures de résistivité apparente à partir d'une série acquise à l'état initial permet de les atténuer fortement. La modélisation est particulièrement sensible au choix du maillage, aux normes appliqués (L1 ou L2) sur les données et les paramètres, et au facteur d'amortissement. Des valeurs trop élevées du facteur d'acceptance tendent à lisser les contrastes au niveau du front de dispersion et augmentent l'impact des effets des bords horizontaux sur D et . A l'inverse, une modélisation contrainte par de faibles valeurs du facteur d'acceptance donnent des résultats plus proches de l'analyse vidéo, mais produit des effets de bosses à l'avant et à l'arrière du front. La vitesse interstitielle u est indépendante du choix des paramètres d'inversion pour l'ERT. Pour les deux méthodes u est toujours inférieure au débit imposé par la pompe, dont le décalage est exprimé par le facteur retard Rf . Les effets de retard résultent de l'adsorption du Na+ sur les surfaces des billes de verre chargées négativement qui retarde le front de dispersion du suivi ERT. Pour le suivi vidéo, la taille importante de la molécule du colorant favorise son piégeage dans les zones où la perméabilité est plus faible, en plus d'une éventuelle affinité avec la surface solide. Les contrastes de conductivité et la stabilité de l'interface créent par la différence de densité entre les fluides testés ici n'ont pas d'influences significatives sur la dispersion qui est dominée par le débit imposé à la pompe. Les estimations du coefficient de dispersion D en fonction du nombre de Péclet sont cohérentes avec la courbe théorique de Bachmat (1968). Cependant la dispersivité augmente pour les vitesses d'écoulement les plus élevées. Les premières expérimentations de terrain réalisées en 2D sur des sables de Fontainebleau présentent l'avantage de s'affranchir des effets de bords inhérents au laboratoire. En revanche la recalibration des données normalisée par la loi de Archie est plus complexe puisqu'il est nécessaire de tenir compte de l'état de saturation de la résistivité des fluides initialement présents. De plus l'erreur importante sur les modèles ne permet pas de déduire une estimation fiable des paramatères de transport u, (ou D), et Rf .

Pour comprendre les mécanismes de déstabilisation, d'arrêt et de coexistence des phases statiques (solide) et coulantes (liquide), nous réalisons expériences et simulations numériques d'empilements granulaires modèles. Trois sujets ont été abordés au cours de cette thèse.<br />L'étude expérimentale de la transition d arrêt d'un empilement après un écoulement de surface met en évidence l'existence de relaxations de durée bien supérieure au temps de relaxation d'un grain sous l'action de son poids. Celle-ci est constituée de phases de relaxation et de réactivations liées à des déplacements corrélés des grains, prises en compte dans un modèle statistique.<br />L'étude numérique d'un empilement incliné en deçà de l'angle d'avalanche met en évidence l'influence du domaine métastable –au-delà de l'angle de repos– sur ses propriétés hystérétiques au cours de cycles quasi-statiques. Le réseau des contacts faibles est très affecté par le passage dans le domaine métastable. Les corrélations entre micro-structure, contrainte et déformation sont discutées.<br />Enfin, l'expérience d'un écoulement non confiné sur plan incliné permet d'étudier la loi d'écoulement dans le cas de la coexistence solide-liquide. Les caractéristiques s´electionnées par l'écoulement (épaisseur, largeur, vitesse) évoluent lentement avec le temps. La prise en compte de l'existence d'une couche basale statique permet de retrouver la corrélation entre épaisseur coulante et vitesse, comme pour un écoulement confiné, et de prédire la morphologie des dépôts (présence ou non de levées).

Ce travail porte sur l'étude par simulation numérique de jets ronds coaxiaux à fort rapport de vitesse avec une attention particulière portée sur trois aspects : la transition de ces jets vers la turbulence, le mélange dans cette configuration d'écoulement et l'influence du nombre de Reynolds sur ces deux précédents points.<br /><br />Dans un premier temps, nous réalisons des simulations numériques directes de jets coaxiaux à un nombre de Reynolds modéré. Ces jets générent une région de recirculation lorsque le rapport des vitesses entre les jets externe et interne dépasse une valeur critique. Les anneaux de Kelvin-Helmholtz intérieurs et extérieurs ont un développement couplé. Ils sont convectés avec une même fréquence de passage contrôlée par la couche cisaillée extérieure. Ensuite, comme pour les jets ronds, des tourbillons longitudinaux contra-rotatifs apparaissent initiant la tridimensionnalisation de l'écoulement. La région de recirculation influence les anneaux de Kelvin-Helmholtz internes en les ralentissant et en les étirant longitudinalement. Deux modèles théoriques prédisant des grandeurs globales du jet (la longueur du cône potentiel interne et la valeur du rapport de vitesse critique au-delà duquel la région de recirculation apparaît) montrent l'influence de l'épaisseur de quantité de mouvement intérieure initiale.<br /><br />Les propriétés de mélange ont ensuite été étudiées en résolvant l'équation de transport d'un traceur simultanément aux équations de Navier-Stokes. Les structures cohérentes de l'écoulement contrôlent le processus de mélange. Les tourbillons longitudinaux augmentent le mélange par un phénomène d'éjection du traceur en périphérie du jet. Cependant, la configuration initiale du jet montre que des poches de traceur non mélangé persistent à la fin de la transition. Les modifications des conditions d'entrée du jet diminuant ces poches sont celles qui permettent une génération plus précoce ou plus intense de structures longitudinales. C'est le cas de la région de recirculation qui étire longitudinalement les structures. De la même façon, un forçage azimutal de la couche cisaillée externe (qui domine la dynamique) améliore nettement le mélange et semble être plus performant qu'un forçage axisymétrique en ce qui concerne le mélange en champ proche.<br /><br />Pour finir, nous avons réalisé des simulations des grandes échelles de jets coaxiaux à hauts nombres de Reynolds. L'auto-similitude des jets coaxiaux en turbulence pleinement développée a permis une validation sur les données expérimentales. Les quantités globales des jets coaxiaux sont fortement dépendantes du nombre de Reynolds jusqu'à ce que celui-ci atteigne une valeur de l'ordre de 10000. Au-delà de cette valeur, ces quantités sont quasi-indépendantes du nombre de Reynolds en raison du phénomène de ``mixing transition'' qui implique une déstabilisation tridimensionnelle immédiate du jet. Cette déstabilisation précoce des couches cisaillées conduit la région de recirculation à un comportement instationnaire. Enfin, cela permet une nette amélioration du mélange avec un phénomène d'éjection du traceur proche de l'entrée du jet.

De nombreuses hypothèses sont classiquement utilisées pour décrire le comportement du fluide dans un contact lubrifié : film continu, viscosité constante dans l'épaisseur, film mince, fluide newtonien... Or, certaines s'avèrent erronées dès lors que l'on s'intéresse aux contacts Elasto- HydroDynamiques fortement glissants ou à la répartition du lubrifiant en sortie de contact. Une approche numérique originale, basée sur un retour aux équations de la mécanique des fluides générale et prenant en compte le couplage fluide/solide et les effets thermiques sont proposés ici dans le but d'apporter des éléments physiques supplémentaires aux modélisations usuelles. Dans un premier temps, l'influence des effets thermiques sur l'évolution du frottement dans les contacts Thermo-EHD est mise en évidence. La présence d'un minimum de frottement pour le cas du glissement pur est expliquée par l'analyse des transferts thermiques entre le lubrifiant et les solides. L'origine des modifications locales d'épaisseur de film observées et l'existence même d'une épaisseur de film lubrifiant pour les cas de vitesse d'entraînement nulle sont alors reliées à la présence d'un fort gradient de viscosité dans l'épaisseur de film. Une comparaison qualitative avec des éléments expérimentaux de la littérature est réalisée, validant les tendances obtenues. Dans un second temps, l'écoulement à surface libre du lubrifiant en périphérie du contact est étudié expérimentalement puis numériquement par une méthode à interface diffuse. Le rôle des effets capillaires est analysé et les résultats numériques confrontés à des résultats issus de la littérature. Un bon accord est obtenu tant qualitativement que quantitativement. Validé par l'étude numérique diphasique (air/lubrifiant) réalisée, un modèle analytique simplifié est alors développé, prédisant une loi de répartition du lubrifiant en sortie de contact. La zone de sortie des contacts EHD est ensuite traitée par un modèle de cavitation vaporeuse et la prise en compte nécessaire de l'air environnant est discutée. Enfin, une première modélisation tridimensionnelle de l'écoulement à surface libre du lubrifiant autour d'un contact ponctuel est réalisée mettant en avant l'influence des effets capillaires et la faisabilité d'une telle approche.

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.

Nous étudions dans cette thèse, composée de deux parties, la modélisation des écoulements en eaux peu profondes par la méthode de Boltzmann sur réseau et l'analyse des propriétés de commandabilité et d'observabilité des modèles obtenus. Dans la première partie, nous nous consacrons d'abord à la modélisation par la méthode de Boltzmann sur réseau des équations de Saint-Venant. En utilisant une linéarisation autour d'un profil d'équilibre, une représentation sous forme d'état des modèles de Boltzmann sur réseau est définie. Cette représentation incorpore les termes de force, et permet une définition complète des entrées (commandes) et des sorties (mesures). Nous représentons ensuite les phénomènes de sédimentation dans les écoulements en eaux peu profondes avec la méthode de Boltzmann sur réseau. Ce modèle défini en une dimension est validé numériquement en le comparant avec un modèle de volumes finis qui résout les équations de Saint-Venant-Exner. Le modèle LB défini est moins gourmand en temps de calcul et plus facile à manipuler que les modèles traditionnels. Dans la deuxième partie, nous traitons de l'analyse des propriétés de commandabilit é et d'observabilité des modèles LB obtenus. La première analyse est faite sur les critères algébriques de Kalmann et permet d'établir la non conservation des propriétés de commandabilité et d'observabilité lorsque l'ordre de réduction du système est augmenté. Une analyse plus approfondie basée sur la détermination les grammiens de commandabilit é et d'observabilité montre également que le constat reste valide pour les méthodes de discrétisation classique. La résolution des grammiens est faite avec des méthodes particulièrement adaptées aux structures creuses et de grande dimension que sont les matrices de la dynamique, de commande et/ou d'observation des modèles LB. Enfin, nous établissons que pour une commande aux frontières classique des canaux d'irrigation en débit et hauteur, la famille de systèmes des modèles LB d'ordre réduit n'est pas uniformément commandable alors qu'avec l'utilisation des variables de scattering comme variables de commande, cette famille devient uniformément commandable.

La démolition de bâtiments par foudroyage est une technique devenue classique. Elle génère des nuisances qui doivent être prises en compte, sous peine de rendre caduc ce type d'activité. Les principales nuisances générées par le foudroyage sont : - les projections générées par l'explosion et la dislocation du bâtiment ; - l'onde de surpression générée par l'explosif ; - la vibration du sol induite par l'effondrement du bâtiment ; - la création et la propagation d'un nuage de poussière. Les deux premiers points sont maitrisés. En ce qui concerne les vibrations, aucune méthode n'existe pour les limiter. Seuls, des fossés creusés permettent de bénéficier d'un effet d'ombre vis-à-vis de certains bâtiments. Les nuisances générées par les poussières et les manières de les limiter commencent à être une question fréquemment posée. Dans le cadre de la destruction d'installation classée pour la protection de l'environnement, cette question devient critique, car ce nuage peut déplacer des particules nocives. Actuellement, aucune information sur les mécanismes de production ou sur les mécanismes de propagation des poussières n'a été élaborée. Cette thèse propose un protocole de mesure pour l'étude des caractéristiques du nuage de particule. Elle présente les premières informations sur la production et la propagation des particules. Et pour finir, elle présente les premières simulations numériques du phénomène.

Une méthode de frontières immergées est développée pour la simulation d'écoulements compressibles et validée au travers de cas-tests spécifiques (réflexion d'ondes acoustiques et quantification de la conservation de la masse dans des canaux inclinés). La simulation aux grandes échelles (LES) d'une cavité transsonique est ensuite présentée. Le bouclage aéro-acoustique, très sensible aux conditions aux limites, est reproduit avec précision par la LES dans le cas où les parois sont immergées dans un maillage structurée. La comparaison des stratégies de modélisation de sous-maille pour cet écoulement transsonique et l'adaptation des filtres en présence de frontières immergées sont également discutées. Le rôle, souvent sous-estimé, du schéma de viscosité artificiel, est quantifié.Dans la dernière partie du manuscrit, des études sont réalisées pour aider au dimensionnement d'un nouveau concept de chambre de combustion où la flamme est stabilisée par la recirculation de gaz brûlés dans une cavité (chambre TVC pour Trapped Vortex Combustor). La modélisation de la combustion turbulente est basée sur une chimie tabulée, couplée à une fonction densité de probabilité présumée (PCM-FPI). L'étude de la dynamique de la flamme est réalisée pour diverses conditions de fonctionnement (débit de l'écoulement principal et présence ou non d'un swirl). Les spécificités de mise en œuvre de la simulation d'un écoulement de ce type sont discutées et un soin particulier est apporté au traitement de la condition de sortie, qui constitue un point sensible de la chaîne de modélisation. Les phénomènes d'instabilités et de retour de la flamme sont mis en évidence ainsi que les modifications à apporter au dispositif afin de minimiser ces effets. L'existence d'un cycle limite acoustique est souligné et une formule permettant d'anticiper le niveau des fluctuations de pression est proposée et validée. Une correction au modèle PCM-FPI est présentée afin de préserver la vitesse de flamme et d'assurer une reproduction plus précise de la dynamique de flamme.