Dans les procédés industriels mettant en œuvre des systèmes solide liquide, le contrôle du processus d'agglomération des particules est fondamental compte tenu de son influence déterminante sur l'évolution de la distribution granulométrique des produits finaux. Ce travail a pour objectif l'étude du processus d'agglomération de suspensions d'alumine, agitées, peu concentrées en s'intéressant tout particulièrement aux aspects physico-chimiques et hydrodynamiques. La principale technique expérimentale utilisée est la turbidimétrie spectrale, méthode optique qui permet de suivre in-situ, l'évolution dans le temps de l'ensemble de la distribution granulométrique d'une suspension. Nous avons ainsi pu identifier expérimentalement trois phénomènes responsables de cette évolution: l'agglomération, la sédimentation et la rupture des agglomérats formés et étudier l'influence des paramètres hydrodynamiques et physico-chimiques. Nous proposons alors une modélisation par l'intermédiaire du calcul d'un noyau d'agglomération original en régime turbulent qui tient compte des interactions hydrodynamiques et des forces d'interaction interparticulaires. Les résultats expérimentaux et ceux du modèle sont confrontés au moyen de simulations numériques de l'évolution dans le temps de la turbidité d'une suspension.

Afin de déterminer la conduction thermique de sols en climat continental humide, l'étude des paramètres minéralogiques, physiques, physico-chimiques et géophysiques est essentielle. Le diffractomètre à rayons X, le microscope électronique à balayage et le terramètre sont utilisés pour la caractérisation et le calcul de la conduction thermique, soit la capacité volumétrique, la conductivité et la diffusivité thermique. La teneur en eau gravimétrique est la caractéristique qui interfère le plus fortement sur le bilan thermique. La minéralogie, quant à elle, joue deux rôles importants: (1) le quartz, minéral dominant dans la plupart des échantillons, possède une forte conductivité thermique et (2) la minéralogie des particules fines inférieures à 63 microns influencent la perméabilité des sols. L'interaction qui existe entre la température du sol, le contenu en eau et la minéralogie est clairement définie. Un lien direct unit les quantités de vermiculite et d'interstratifiés chlorite-smectite, potentiellement absorbantes, et les quantités d'eau retrouvées dans certains échantillons. Ainsi, la définition du bilan thermohydrodynamique représentatif des sols doit inclure une étude minéralogique. ______________________________________________________________________________ MOTS-CLÉS DE L’AUTEUR : Conduction thermique, Sols agricoles, Vermiculite, Diffraction aux rayons X, Microscopie électronique à balayage, Résistivité électrique.

La microscopie à force atomique est utilisée dans de nombreuses applications. Ce travail a ici pour but d'étudier, à l'échelle du nanomètre, les propriétés mécaniques de monocouches et de multicouches de phospholipides avec un microscope à force atomique en mode dynamique. Le comportement non linéaire du système pointe-levier, proche des couches de phospholipides, est étudié avec un modèle de granulation. Le MFA est aussi utilisé pour étudier l'influence d'un peptide membranaire, la syntaxine, et la structure des couches de phospholipides autoassemblées sur du mica. Dans un second temps, nous étudions l'oscillation d'un microlevier de MFA clampé en milieu liquide. Les résultats sont comparés aux prédictions théoriques des simulations numériques de la résolution des équations de Navier-Stokes à trois dimensions, pour mesurer l'influence du mouvement du fluide sur le comportement d'un levier oscillant de MFA. Nous montrons aussi comment la force hydrodynamique peut être réduite d'au moins un ordre de grandeur en usinant, à l'aide d'un faisceau d'ions focalisés (FIB), la largeur du levier. Enfin, nous présentons une description analytique qui détermine le mouvement d'un levier excité acoustiquement en milieu liquide. La troisième partie est une étude du comportement dynamique d'un nanoménisque avec une nanopointe oscillante. La nanoaiguille, découpée avec un<br />faisceau d'ions focalisés, est approchée d'une interface air-liquide et seule son extrémité oscille dans le liquide (Eau-Glycérol-PDMS). Il est montré qu'il est possible d'imager une interface de liquide avec une résolution à l'échelle du nanomètre.

Reproduction de : Thèse de doctorat : Dynamique des fluides : Toulouse, INPT : 2007. / Titre provenant de l'écran-titre. Bibliogr. 58 réf.

Reproduction de : Thèse de doctorat : Génie des procédés et de l'environnement : Toulouse, INPT : 2008. / Titre provenant de l'écran-titre. Bibliogr. 98 réf.

Cette thèse a pour objectif de détecter les structures turbulentes aux grandes échelles présentes dans une couche de mélange de type Rayleigh-Taylor incompressible à faible nombre d'Atwood. Diverses grandeurs statistiques conditionnées par la présence de ces structures ont été obtenues, et il est désormais possible de les comparer avec les résultats des modèles statistiques turbulents dits bi-structure, tel le modèle 2SFK développé au CEA. Afin de réaliser les simulations numériques directes du mélange turbulent, un code numérique tridimensionnel incompressible à densité variable a été développé. Ce code a été parallélisé dans les trois directions. Plusieurs méthodes de détection de structure ont été conçues et testées. Bien que toutes ces méthodes présentent différents intérêts, seule la plus efficace vis-à-vis de nos critères de détection a été gardée pour faire des simulations à forte résolution (plus d'un milliard de mailles, 1024^3). Un filtrage temporel de la vitesse verticale est utilisé dans cette méthode de détection afin de : 1) corriger les distorsions dues aux points d'arrêt et zones de recirculation dans l'écoulement, 2) minimiser l'effet de la turbulence aux petites échelles et mieux mettre en évidence les grandes échelles, 3) introduire un effet mémoire permettant de prolonger la bimodalité du champ de détection depuis les zones laminaires extérieures jusqu'au centre de la zone de mélange turbulent. Plusieurs simulations numériques directes 1024^3 ont été effectuées. Les résultats viennent conforter ceux obtenus avec le modèle bi-structure 2SFK et justifient une étude plus poussée des grandeurs statistiques en vue de sa validation.

Les supernovæ de type Ia (SNe Ia) sont devenues un outil important pour retracer l'expansion de notre Univers, leur étude est donc importante pour la cosmologie. Le modèle le plus populaire est celui de l'explosion d'une naine blanche (NB) accrétante dont la contraction relance la combustion sous la forme d'une déflagration subsonique, qui transiterait ensuite en une détonation supersonique. Ce scénario de détonation retardée repose sur un mécanisme physique de Transition Déflagration-Détonation (TDD) encore très mal compris, que nous étudions dans cette thèse.Les modèles actuels de détonation retardée reproduisent les observations en se fondant sur le mécanisme des gradients de Zel'dovich. Cependant, les échelles d'ignition n'étant pas résolues, ces simulations n'expliquent pas à elles seules la TDD, phénomène mal compris, même sur Terre, lorsqu'il s'agit de milieux non-confinés. D'autre part, ce mécanisme requiert une turbulence trop intense et impose des conditions probablement trop restrictives.C'est dans ce contexte que nous avons proposé un nouveau mécanisme de TDD: le chauffage acoustique de l'enveloppe du progéniteur. Un modèle simplifié, en géométrie plane, permet de mettre en évidence l'amplification d'ondes acoustiques (générés par une flamme turbulente) dans un gradient de densité similaire à ceux d'une NB. Selon leur fréquence et leur amplitude, leur amplification peut aller jusqu'à la formation d'un choc suffisamment fort pour initier une détonation. Ensuite, ce mécanisme est analysé en géométrie sphérique dans le cadre plus réaliste d'une NB en expansion. Une étude paramétrique montre la validité de notre mécanisme sur une gamme raisonnable de fréquences et d'amplitudes acoustiques.Finalement, quelques simulations MHD 2D et 3D, où l'on recherche une source de perturbations acoustiques, sont présentées pour démontrer le caractère réaliste de notre nouveau mécanisme de TDD.

L'ablation laser dans l'air à la pression atmosphérique est souvent appliquée dans l'industrie, l'analyse chimique, la chirurgie, ... Pour un développement plus approfondi des technologies laser basées sur l'effet d'ablation, il est nécessaire de mieux comprendre les phénomènes à l'origine de l'interaction laser-matière. Lors d'un impact du faisceau laser sur la surface d'un matériau, une plume plasma se forme au dessus de la cible. Ce plasma contient des électrons, des atomes et des ions du matériau évaporé, ainsi que du gaz ambiant, s'il est présent. Lors de l'impulsion laser, cette plume absorbe une grande partie d'énergie du faisceau laser, réduisant ainsi la quantité du rayonnement laser qui arrive à la surface de la cible. Les dimensions, ainsi que les paramètres de cette plume plasma évoluent très rapidement avec le temps. L'étude de la dynamique et des paramètres de cette plume est très importante, parce qu'ils influent sur tous les processus physiques ayant lieu à la surface du matériau traité. Nous avons étudié l'expansion de la plume plasma formée lors de l'ablation des échantillons métalliques (Al, Ti, Fe) par faisceau laser Nd :YAG nanoseconde (durée d'impulsion : 5.1 ns, longueur d'onde : 1064 nm, irradiance de l'ordre de grandeur de GW/cm2) dans l'air à la pression atmosphérique en utilisant la technique d'imagerie rapide. Cette technique nous a permis d'observer l'évolution spatio-temporelle de la plume au début de son expansion. Les résultats obtenus indiquent que la plume d'ablation laser a une structure : deux régions ont été distinguées - le cœur et la périphérie de la plume. La dynamique de ces deux régions de la plume a été étudiée et les vitesses de leur expansion ont été déterminées. En plus, nous avons examiné l'influence de l'irradiance laser, ainsi que l'influence de la composition de la cible sur la dynamique de la plume. D'autre part, nous avons développé des modèles numériques sous COMSOL Multiphysics pour simuler le processus d'ablation laser. Un modèle thermique a été utilisé pour simuler l'interaction laser - cible. Les résultats de ce modèle ont été employés comme conditions à la limite du modèle hydrodynamique, utilisé pour simuler l'expansion de la plume du plasma dans l'air ambiant. Deux approches ont été proposées : approche microscopique et macroscopique. Les résultats de simulation basée sur l'approche macroscopique donnent l'évolution temporelle et distribution spatiale (1D) des paramètres de la plume : masse volumique, vitesse, pression.

Cette thèse s'intéresse à la turbulence dans l'hélium (4He) superfluide, à des températures comprises entre 1.15 K et la température de transition superfluide, Tlambda = 2.17 K, ce qui correspond à une fraction de superfluide rho_s/rho comprise entre 97.6% et 0%. Il s'agit d'un travail essentiellement expérimental dont le but est de comparer la turbulence classique et la turbulence quantique, à l'aide de mesures locales de fluctuations de vitesse et de vorticité. Ces mesures sont complétées par l'analyse de champs de vitesse issus de simulations numériques.Nous avons développé une instrumentation spécifique, adaptée aux écoulements cryogéniques: des tubes de Pitot miniatures, dont les dimensions effectives ont pu être rendues sub-millimétriques, et un capteur original, basé sur la déflection d'une micro-poutrelle (300 microns x 100 microns x 1 microns) mesurée à l'aide d'un micro-résonateur supraconducteur dans la gamme de fréquence du GHz. Un premier prototype de ce capteur, micro-fabriqué à partir d'une galette de silicium, a été réalisé en salle blanche puis validé dans une conduite cryogénique. La résolution spatiale obtenue est du même ordre que celle des meilleurs anémomètres en He II, et il devrait être possible de l'améliorer d'une décade.Les tubes de Pitot ont été placés dans les souffleries superfluides TSF et TOUPIE. La première, fruit d'une collaboration nationale, a fourni un écoulement stationnaire de grille (Rlambda= 250, 1.65 K < T < 2.6 K, rho_s/rho < 80%) et un sillage proche. La seconde soufflerie, refroidie pour la première fois dans le cadre de cette thèse, a fourni un écoulement de sillage lointain (Rlambda = 1100, 1.55 K < T < 4.2 K, rho_s/rho < 86%).Les mesures ont mis en évidence dans ces écoulements des similarités fortes avec les écoulements classiques, aux échelles inertielles: spectre de vitesse en k^{-5/3}, constante de Kolmogorov et taux de turbulence identiques aux écoulements classiques, loi des 4/5, exposants anormaux pour les fonctions de structures des incréments de vitesse (intermittence). À plus petite échelle, les simulations numériques (1.15 K < T < 2.1565 K) mettent en évidence un comportement exotique : l'énergie s'accumule et tend vers l'équipartition, ce qui se traduit par un spectre de vitesse simulé en k^2. Ce phénomène s'accentue à basse température. Enfin, des mesures locales de fluctuation de vorticité ont été réalisées à l'aide de pinces à second son sur une gamme de température comprise entre 1.69 K (rho_s/rho=77%) et 2.01 K (rho_s/rho = 42%). Nous avons observé un raidissement de la pente des spectres de vorticité lorsque la température diminue. Ce résultat peut être interprété comme une conséquence du phénomène d'équipartition mis en évidence dans les simulations numériques à petite échelle.

Le travail présenté ici avait pour but le développement d'un schéma de type Galerkin discontinu (GD) d'ordre élevé pour la résolution des équations de la dynamique des gaz écrites dans un formalisme Lagrangien total, sur des maillages bi-dimensionnels totalement déstructurés. A cette fin, une méthode progressive a été utilisée afin d'étudier étape par étape les difficultés numériques inhérentes à la discrétisation Galerkin discontinue ainsi qu'aux équations de la dynamique des gaz Lagrangienne. Par conséquent, nous avons développé dans un premier temps des schémas de type Galerkin discontinu jusqu'à l'ordre trois pour la résolution des lois de conservation scalaires mono-dimensionnelles et bi-dimensionnelles sur des maillages déstructurés. La particularité principale de la discrétisation GD présentée est l'utilisation des bases polynomiales de Taylor. Ces dernières permettent, dans le cadre de maillages bi-dimensionnels déstructurés, une prise en compte globale et unifiée des différentes géométries. Une procédure de limitation hiérarchique, basée aux noeuds et préservant les extrema réguliers a été mise en place, ainsi qu'une forme générale des flux numériques assurant une stabilité globale L2 de la solution. Ensuite, nous avons tâché d'appliquer la discrétisation Galerkin discontinue développée aux systèmes mono-dimensionnels de lois de conservation comme celui de l'acoustique, de Saint-Venant et de la dynamique des gaz Lagrangienne. Nous avons noté au cours de cette étude que l'application directe de la limitation mise en place dans le cadre des lois de conservation scalaires, aux variables physiques des systèmes mono-dimensionnels étudiés provoquait l'apparition d'oscillations parasites. En conséquence, une procédure de limitation basée sur les variables caractéristiques a été développée. Dans le cas de la dynamique des gaz, les flux numériques ont été construits afin que le système satisfasse une inégalité entropique globale. Fort de l'expérience acquise, nous avons appliqué la discrétisation GD mise en place aux équations bi-dimensionnelles de la dynamique des gaz, écrites dans un formalisme Lagrangien total. Dans ce cadre, le domaine de référence est fixe. Cependant, il est nécessaire de suivre l'évolution temporelle de la matrice jacobienne associée à la transformation Lagrange-Euler de l'écoulement, à savoir le tenseur gradient de déformation. Dans le travail présent, la transformation résultant de l'écoulement est discrétisée de manière continue à l'aide d'une base Éléments Finis. Cela permet une approximation du tenseur gradient de déformation vérifiant l'identité essentielle de Piola. La discrétisation des lois de conservation physiques sur le volume spécifique, le moment et l'énergie totale repose sur une méthode Galerkin discontinu. Le schéma est construit de sorte à satisfaire de manière exacte la loi de conservation géométrique (GCL). Dans le cas du schéma d'ordre trois, le champ de vitesse étant quadratique, la géométrie doit pouvoir se courber. Pour ce faire, des courbes de Bézier sont utilisées pour la paramétrisation des bords des cellules du maillage. Nous illustrons la robustesse et la précision des schémas mis en place à l'aide d'un grand nombre de cas tests pertinents, ainsi que par une étude de taux de convergence.