Cette thèse développe des méthodes d'ondelettes originales en vue de simuler des écoulements incompressibles.<br />Nous commencerons par présenter une certaine manière de concevoir le phénomène de la turbulence dans les fluides, puis nous ferons une introduction à la théorie des ondelettes.<br />Dans le but de construire des ondelettes 2D et 3D adaptées aux écoulements fluides, nous reprenons en les enrichissant les travaux de P-G Lemarié-Rieusset et K. Urban sur les ondelettes à divergence nulle. Nous mettons en évidence l'existence d'algorithmes rapides associés.<br />Par la suite, nous démontrons qu'il est possible d'utiliser ces ondelettes à divergence nulle pour définir la décomposition de Helmholtz d'un champ de vecteurs 2D ou 3D quelconque. Cette décomposition est définie par un algorithme itératif dont nous prouvons la convergence pour des ondelettes particulières. L'optimisation de la convergence fait ensuite l'objet d'une étude poussée.<br />Tous ces ingrédients permettent de définir une nouvelle méthode de résolution des équations de Navier-Stokes incompressible, dont nous prouvons la faisabilité sur un cas test.<br />On applique également la décomposition en ondelettes à divergence nulle à l'analyse de champs d'écoulements turbulents 2D et 3D, ainsi qu'à la compression dans une méthode d'Extraction de Structures Cohérentes.

Dans cette thèse, nous développons et mettons en oeuvre une méthode numérique de discrétisation spatiale entièrement tridimensionnelle afin d'étudier la transition des écoulements visqueux et incompressibles dans un canal infiniment long, d'un état stable et bidimensionnel à un état tridimensionnel. Le principe de stabilité linéaire permet de déterminer l'apparition d'une telle bifurcation et la théorie des systèmes dynamiques montre que l'écoulement au voisinage de la solution stable bidimensionnelle considérée tend alors à suivre une direction privilégiée. Dans certains cas il en résulte un écoulement tridimensionnel et périodique qu'il n'est possible de décrire qu'en adoptant une approche entièrement tridimensionnelle. Nous avons adopté une approche combinant la robustesse de la méthode des éléments finis à la précision des méthodes de Fourier. La théorie de la méthode de discrétisation est expliquée, un code est validé en utilisant plusieurs bancs d'essai et la description qualitative du comportement local de l'écoulement après bifurcation est présentée.

Les équations fondamentales décrivant la dynamique de l'océan sont en théorie les équations de Navier-Stokes sur une sphère en rotation, auxquelles il faut a jouter une équation d'état pour la densité, et des équations de transport-diffusion pour les traceurs. Toutefois, un certain nombre de considérations physiques et de limitations pratiques ont nécessité le développement de modèles plus simples. En effet, un certain nombre d'hypothèses simplificatrices sont pleinement justifiées du point de vue de la physique des mouvements océaniques, dont les principales sont les approximations de couche mince et de Boussinesq. D'autre part, étant donné les dimensions des bassins océaniques (plusieurs centaines à plusieurs milliers de kilomètres), les coûts de calculs sont un facteur pratique extrêmement limitant. On est, à l'heure actuelle, capable de simuler l'océan mondial avec une résolution de l'ordre de dix kilomètres, en utilisant des modèles dits aux équations primitives, dont le coût de calcul est bien inférieur à celui des équations de Navier-Stokes. On est donc bien loin d'une modélisation complète des phénomènes décrits par ces équations, qui nécessiterait en théorie de considérer des échelles de l'ordre du millimètre. Les équations primitives sont issues des équations complètes de la mécanique des fluides en effectuant l'approximation hydrostatique, justifiée par la faible profondeur des domaines considérés au regard de leur dimension horizontale. Dans cette thèse, nous considérerons les équations de Navier-Stokes (NS) qui sont le coeur du modèle complet évoqué ci-dessus, sans prendre en compte les équations de la densité et des traceurs (salinité, température, etc.). Nous utiliserons l'approximation hydrostatique dans le chapitre 10, et le modèle sera naturellement appelé Navier-Stokes hydrostatique (NSH). Il correspond aux équations primitives dans lesquelles on ne prendrait pas en compte la densité et les traceurs. C'est dans ce cadre que se situe le travail présenté dans cette thèse, avec l'objectif à moyen terme de pouvoir coupler de façon rigoureuse et efficace les équations de Navier-Stokes avec les équations primitives. Dans une première partie, on présentera quelques rappels sur les équations de Navier-Stokes, leur discrétisation, ainsi que le cas-test de la cavité entraînée qui sera utilisé dans tout ce document. Dans une deuxième partie, on met en œuvre les méthodes de Schwarz sur les équations de Stokes et Navier-Stokes, en dérivant notamment des conditions absorbantes exactes et approchées pour ces systèmes. Enfin, dans une troisième partie, on proposera des pistes vers le couplage Navier-Stokes/Navier-Stokes hydrostatique décrit ci-dessus. / Fundamental equations describing the ocean dynamic are theoretically Navier-Stokes equations over a rotating sphere, whom need to add a state equation for the fluid density, and advection-diffusion equations for tracers. However, some physical considerations and practical limitations required to developped more simple models. Indeed, some simplifying hypotheses are well justified from a ocean dynamic point of view, whose principal ones are thin layer and Boussinesq approximations. On the other hand, considering the dimensions of oceans (from serveral hundreds to serveral thousands kilometers), computations costs are a very practical limitating factor. We are, by now, able to simulate the global ocean with about ten kilometers large grid mesh. This is very far from a complete modelisation of all phenomenes decribed by the Navier-Stokes equations, which require to consider scales of milimeters order. Primitives equations derive from complete equations describing fluid mecanics, by doing the hydrostatic approximations, which is justified by the low deepness of considered domains with regard to their horizontal dimension. In this thesis, we are considering Navier-Stokes equations (NS) which are the heart of the complete modele mentionned previously, without holding in account density and tracers equations. We will use the hydrostatic approximations, and the resulting equations will be named as hydrostatic Navier-Stokes equations (NSH).The mid term objective is to couple carefully Navier-Stokes equations with primitive equation. In a first part, we will remind few results for Navier-Stokes equations, their discretization, and the lid-driven cavity which wil be used as a test-case. In a second part, we will use Schwarz method with Stokes and Navier-Stokes equations, deriving in particular exact and approched absorbing interface conditions for these systems. Finally, in a third part, we shall propose first results towards coupling Navier-Stokes and hydrostatic Navier-Stokes equations.

La simulation de fluide fait l’objet de recherches actives en infographie. Largement utilisée dans le domaine des jeux vidéos ou de l’animation, elle permet de simuler le comportement des liquides, des gaz et autres phénomènes pouvant être apparentés à un fluide. Pour cela, la simulation de fluide dispose d’outils de calcul numériques adaptés, permettant de produire des animations visuellement réalistes pour un temps de calcul raisonnable. Ce mémoire décrit les deux principales approches utilisées en simulation de fluide : l’approche eulérienne et l’approche lagrangienne, ainsi que certains outils numériques associés, que sont les différences finies et les fonctions de lissage. Chaque approche et chaque outil numérique possède ses avantages et ses inconvénients. Les noyaux constants par morceaux constituent un nouvel outil de calcul numérique et ouvrent de nouvelles possibilités à la simulation de fluide. Ils seront étudiés en détails puis intégrés dans une simulation de fluide eulérienne. L’atout notable qu’apportent les noyaux constants par morceaux est la possibilité d’augmenter la précision des calculs là où cela est jugé nécessaire dans la simulation. En augmentant la précision des calculs aux endroits clés, où sont susceptibles d’apparaitre des effets visuellement attrayants comme les tourbillons ou les remous, nous améliorons la qualité des animations.

.

Nous proposons un modèle mathématique intégré du poumon dont l'approche globale repose sur une modélisation multibloc. En effet, on décompose en trois niveaux l'arbre bronchique qui s'étend sur vingt quatre générations de bronches allant de la trachée aux alvéoles. Au premier niveau (les six premières générations), a lieu un écoulement de Navier-Stokes, qui est simulé directement. Au deuxième niveau (de la génération sept à la génération dix sept), les flux à travers les bronches sont régis par la loi de Poiseuille. La linéarité de cette loi nous permet de condenser cette partie de l'arbre et de proposer des conditions aux bords dissipatives adaptées à la similation de la ventilation et permettant d'éviter le maillage de cette partie géométriquement complexe. Le dernier niveau du modèle, prend en compte la partie distale de l'arbre qui est la zone alvéolaire. Elle est composée des acini, qui agissent comme un ensemble de petites pompes et dont l'effet macroscopique est le moteur même de la respiration. A ce niveau, on propose les déplacements d'un piston comme modèle simplifié des mouvements du diaphragme pulmonaire. Dans un premier temps, on se place dans le cadre particulier des équations de Stokes et on s'intéresse au couplage des deux premiers compartiments, dont la validité est illustrée par des tests numériques. On explique également le calcul de la résistance globale équivalente qui intervient dans le calcul de la condition aux limites qui remplace la zone condensée. L'étude est ensuite généralisée au cas des équations de Navier-Stokes. La difficulté réside dans le contrôle du flux d'énergie cinétique, on introduit alors une classe de conditions aux limites, qu'on désigne par dissipatives essentielles, pour lesquelles la trace du champ de vitesse sur les sections d'entrée et de sorties vit dans un espace de dimension fini, et pour lesquelles on prouve des résultats d'existence de solutions faibles locales en temps pour données quelconques et globales en temps pour données petites. Pour le cas de conditions dites dissipatives naturelles, c'est à dire sans contrainte sur la trace du champ de vitesse, on a existence de solutions faibles locales en temps pour données petites et globales en temps pour données plus petites, mais seulement en dimension deux. Cependant, on prouve pour ces conditions aux limites, que pour une classe de solutions plus régulières on a l'existence d'une unique solution locale en temps ainsi que l'existence d'une solution globale en temps pour données petites. Pour le couplage global, incluant le piston, on prouve l'existence de solutions faibles locales en temps pour des données quelconques en ce qui concerne les conditions aux limites dissipatives essentielles, tandis que pour les conditions dissipatives naturelles, on obtient l'existence de solutions locales en temps pour données petites et toujours seulement en dimension deux. Finalement, on propose une discrétisation en temps du problème global et on établit un bilan énergétique à l'ordre 1 pour le problème régulier en espace et discrétisé en temps. Nous présentons ainsi plusieurs simulations numériques bi-dimensionnelles correspondants aussi bien à un poumon sain que pathologique et notamment asthmatique.

La thèse est composée de deux parties indépendantes.<br />Dans la première partie, on étudie le comportement asymptotique de problèmes elliptiques et paraboliques à données $L^1+W^{-1,p'}$ (respectivement $L^1+L^p(0,T;W^{-1,p'})$ dans le cas parabolique), dans des domaines devenant infiniment grands. En utilisant le cadre des solutions renormalisées et les résultats de régularité des solutions pour de telles données, on prouve, sous certaines hypothèses structurelles sur les variables d'espace, des résultats de convergence dans les espaces de régularité des solutions.<br />Dans la seconde partie, dans le cas de la dimension $2$, on étudie des systèmes de type Boussinesq. Ces systèmes dérivent de modèles de mécanique des fluides et consistent en un couplage des équations de Navier-Stokes incompressibles et de l'équation de la chaleur. On s'intéresse essentiellement aux questions d'unicité de la solution, particulièrement délicate à prouver du fait du couplage très non linéaire entre les équations. On travaille dans le cadre des solutions faibles pour les équations de Navier-Stokes et dans le cadre des solutions renormalisées pour des problèmes paraboliques pour l'équation de la chaleur. On établit tout d'abord des résultats de régularité pour ces équations puis on montre plusieurs résultats d'existence et d'unicité de la solution du système pour de petites données.

Ce travail de thèse concerne les méthodes numériques à base d'ondelettes pour la simulation de la turbulence incompressible. L'objectif principal est la prise en compte de conditions aux limites physiques dans la résolution des équations de Navier-Stokes. Contrairement aux travaux précédents où la <i>vorticité</i> était décomposée sur base d'ondelettes <i>classiques</i>, le point de vue qui est adopté ici vise à calculer le champ de <i>vitesse</i> de l'écoulement sous la forme d'une série d'ondelettes à divergence nulle. On est alors dans le cadre des équations de Navier-Stokes incompressibles en formulation <i>vitesse-pression</i>, pour lesquelles les conditions aux limites sur la <i>vitesse</i> s'écrivent explicitement, ce qui diffère de la formulation <i>vitesse-tourbillon</i>. Le principe de la méthode développée dans cette thèse consiste à injecter directement les conditions aux limites sur la base d'ondelettes. Ce travail prolonge la thèse de E. Deriaz réalisée dans le cas périodique. La première partie de ce travail a donc été la définition et la mise en œuvre de nouvelles bases d'ondelettes à divergence nulle ou à rotationnel nul sur $[0,1]^n$, permettant la prise en compte de conditions aux limites, à partir des travaux originaux de P. G. Lemarié-Rieusset, K. Urban, E. Deriaz et V. Perrier. Dans une deuxième partie, des méthodes numériques efficaces utilisant ces nouvelles ondelettes sont proposées pour résoudre différents problèmes classiques : équation de la chaleur, problème de Stokes et calcul de la décomposition de Helmholtz-Hodge en non périodique. L'existence d'algorithmes rapides associés rend les méthodes compétitives. La dernière partie est consacrée à la définition de deux nouveaux schémas de résolution des équations de Navier-Stokes incompressibles par ondelettes, qui utilisent les ingrédients précédents. Des expériences numériques menées pour la simulation d'écoulement en cavité entraînée en dimension deux ou le problème de la reconnection de tubes de vortex en dimension trois montrent le fort potentiel des algorithmes développés.

L'objectif des travaux présentés dans la thèse concerne la résolution par itérations de systèmes algébriques à grande échelle. Ces systèmes sont issus de la discrétisation par éléments finis de problèmes aux limites. Dans la majorité des cas en 3D, la phase de résolution s'avère l'étape la plus exigeante en terme de ressources informatiques. Ainsi, il est impératif de développer des méthodes itératives efficaces et robustes pour un large éventail de problèmes aux limites. Dans cette thèse, nous nous plaçons dans le cadre des méthodes itératives de Krylov à préconditionneur variable, c'est-à-dire autorisant une flexibilité au niveau du choix du préconditionnement en cours d'itération. Nous visons principalement des problèmes de type convection-diffusion, d'élasticité et de Navier-Stokes discrétisés par des éléments finis quadratiques. Afin de réduire les coûts inhérents aux éléments quadratiques, nous proposons une méthode de résolution multi-niveaux basée sur la hiérarchie naturelle entre les éléments finis linéaires et quadratiques d'où le nom de méthode hiérarchique. Elle possède plusieurs points en commun avec les méthodes multi-grilles mais a l'avantage de s'appliquer aux géométries complexes et aux maillages non-structures. L'utilisation de cette méthode comme préconditionneur à une méthode de Krylov à préconditionneur variable permet d'obtenir une méthode très efficace. L'autre partie de la thèse, concerne la résolution globale et itérative des systèmes de type point selle. Ces systèmes proviennent de la discrétisation des équations linéarisées du problème de Navier-Stokes. La résolution efficace de ces systèmes joue un rôle majeur dans le traitement numérique des équations de Navier-Stokes. Pour cela, nous avons mis en place un préconditionneur adroite de format triangulaire par bloc. Pour rendre ce préconditionneur efficace, nous avons fait appel à trois ingrédients : l'ajout du terme rV(div(u)) aux équations continues de Navier-Stokes, une résolution efficace en vitesse par la méthode hiérarchique et une bonne approximation du complément de Schur. Nos tests numériques montrent l'efficacité des méthodes présentées dans ce travail.

Cette thèse porte sur la simulation numérique des écoulements bi-fluides par l'approche à interface diffuse. La description mathématique d'un écoulement bi-fluide par l'approche à interface diffuse consiste en une équation de Navier-Stokes modifiée couplée à un modèle de capture de l'interface mobile entre les deux fluides. Dans cette thèse, pour la capture de l'interface mobile nous avons porté notre attention sur le modèle de Allen-Cahn. En premier lieu nous nous sommes intéressés de prime abord à la résolution numérique du système semi-discrétisé en temps et totalement implicite Navier-Stokes/Allen-Cahn (NSAC). Pour ce faire nous avons développé un algorithme itératif basé sur la méthode du point fixe. Nous avons montré qu'à chaque itération, le système d'équations (contenu dans l'algorithme) est bien défini ; de plus, la solution de l'équation de Allen-Cahn satisfait le principe du maximum. Par la suite nous avons montré que l'algorithme de point fixe converge et que sa limite est la solution du système NSAC semi-discret. Si cette première méthode itérative nous a donné une méthode de résolution, elle n'est pas satisfaisante quant à la performance. En second lieu nous proposons un nouveau schéma de discrétisation en temps à pas fractionnaire inconditionnellement stable. Utilisant une approche de type point fixe (une projection couplée) nous avons montré la convergence à chaque pas de temps et que la limite correspond à la solution du système semi-discret (donnant ainsi l'existence et l'unicité de la solution). Nous concluons enfin avec des applications numériques aux fins d'illustrer la pertinence et les potentielles limitations du modèle d'une part, puis les performances de notre méthode de résolution du système Navier-Stokes/Allen-Cahn d'autre part.