UPThèses

L'enveloppe de vol des hélicoptères est limitée par le décrochage dynamique, que le contrôle d'écoulement permet de surmonter. Malgré de nombreuses études, aucun dispositif de contrôle ne s'est révélé utilisable sur un rotor. Conçu par l'ONERA, un Vortex Generator déployable (DVG) s'est révélé efficace pour limiter les effets du décrochage dynamique d'un profil OA209 en oscillation de tangage. Mais l'étude de cet actionneur est nécessaire avant son application. Le travail exposé dans le présent mémoire s'attache à reproduire par simulations numériques l'écoulement du décrochage dynamique contrôlé par DVG afin d'identifier les phénomènes physiques en jeu. Après une première partie bibliographique dédiée au contrôle de décrochage dynamique et aux Vortex Generators, une seconde partie a été consacrée à valider et étudier la simulation numérique du contrôle du décrochage statique. L'écoulement induit a pu être ainsi caractérisé par les interactions tourbillonnaires présentes du fait de l'épaisseur du DVG et qui réduisent l'efficacité du contrôle. Cette analyse a permis d'élaborer une modélisation préliminaire du DVG. Dans une troisième partie, la validation et l'étude de la simulation numérique du contrôle du décrochage dynamique a mis en évidence un effet de contrôle similaire au cas statique, et la décomposition en modes propres orthogonaux de l'écoulement contrôlé a montré une altération du mode lié au tourbillon de décrochage dynamique. Ce travail valide la simulation numérique mise en place, laisse entrevoir des perspectives d'amélioration du dispositif de contrôle et permet la simulation numérique ultérieure du contrôle de décrochage dynamique de voilure tournante.

Nous présentons une méthode numérique particulaire pour résoudre les équations de Navier-Stockes en formulation vitesse-pression pour la modélisation bidimensionnelle des écoulements incompressibles en présence d'obstacles. Cette méthode s'appuie sur la formulation Smoothed particles hydrodynamics pour le calcul du transport des moments. L'absence de maillage structuré permet de traiter des domaines d'écoulement avec des frontières complexes. Par ailleurs, le calcul du champ de pression est réalisé par la résolution d'une équation de poisson assurant l'incompressibilité de l'écoulement et les conditions aux limites sont renforcées par l'utilisation de la méthode des intégrales de frontières. Cette dernière méthode est connue pour être pénalisante en termes de temps CFU. Pour remédier à ce problème, les contributions des termes sources de l'équation de poisson pour la pression et de l'équation de Helmhotz généralisée pour la vitesse sont calculées en superposant une grille cartésienne au domaine de l'écoulement et en utilisant une méthode de différences finies. Les différentes étapes de construction de la méthode que nous proposons ont été validées par l'étude de plusieurs cas académiques parmi lesquels l'écoulement dans une cavité, la rupture de barrage ou encore l'écoulement derrière un cylindre. En plus de son utilisation classique pour la modélisation des écoulements, notre méthode a été utilisée pour reconstruire les champs de vitesse et de pression à partir d'un champ de vitesse mesurée expérimentalement par PIV et appliquée au cas de l'écoulement autour d'un profil NACA en mouvement.

Ce travail propose un modèle mathématique permettant de décrire le comportement thixotrope des suspensions de laponite. Ce modèle de comportement est élaboré afin de contrôler les propriétés physicochimiques de la suspension. Dans un premier temps, nous avons utilisé un modèle simple cinétique découlant directement du modèle d'Herschel-Bulkley. Par la suite nous avons mis en avant le modèle à double cinétique développé au cours de ce travail. Par la suite grâce aux techniques de mesures rhéologiques classiques, les protocoles de détermination des paramètres de ces deux modèles ont été mis en place. Une étude à l'échelle locale du comportement du fluide permettait de mieux comprendre comment s'organise la structure de la suspension en écoulement. Dans cette perspective nous avons, grâce à la technique de vélocimétrie par images de particules (PIV), réalisé des mesures de champs de vitesses locales. Ceci nous a permis d'observer le comportement de la suspension dans ses différents états, en écoulement et au repos. De plus, cette étude nous a permis de mettre en évidence dans le cas de la géométrie employée "disques coaxiaux", la formation d'un cône statique dont la taille dépend de la vitesse imposée au disque tournant. L'étude de la validité et de la bonne adéquation du modèle développé passe par une étude comparative mettant en jeux d'une part le modèle déjà existant d'Herschel-Bulkley et notre modèle. Une étude numérique faisant appel au code de calcul commercial Star-CD a été effectuée. Les deux modèles utilisés ont été introduits dans le code et des simulations ont été effectuées. La comparaison entre les résultats des mesures de PIV et les simulations ont mis en évidence que le modèle mis en place est plus représentatif de la physique de l'écoulement de la suspension de laponite que le modèle d'Herschel-Bulkley.

Cette thèse a pour objectif de reconstruire la conductivité isotrope d'un domaine borné à partir des données sur le bord. Pour cela, nous considérons le problème de conductivité inverse et on utilise les méthodes du Dbar de Nachman et de Brown-Uhlmann. A partir des données aux bords, ces méthodes consistent à calculer une fonction complexe dite la transformée scattering, ensuite calculer la conductivité en résolvant une équation (qui contient la transformée scattering) dite l'équation Dbar.Pour des raisons de stabilité, nous approchons les transformées scattering par plusieurs façons et nous étudions l'erreur de ces approximations. Nous montrons la stabilité des méthodes Dbar via les approximations des transformées scattering. Nous étudions en détails le cas des conductivités radiales et nous obtenons des expressions explicites des approximations des transformées scattering. Nous utilisons la méthode de Vainikko pour la résolution numérique de l'équation Dbar et nous introduisons un schéma de point fixe et nous étudions sa convergence. Les résultats numériques obtenus montrent que la méthode est efficace et justifient nos résultats théoriques.