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.

L'objet de ce travail est de présenter une partie des travaux entrepris au laboratoire CELIA (CEA, CNRS, Université Bordeaux I) dans le domaine de la modélisation numérique des écoulements fortement compressibles. Cette activité au sein de l'équipe Interaction-Fusion par Confinement Inertiel-Astrophysique, a eu pour objectif principal la mise au point et le développement de schémas numériques robustes dédiés à la simulation numérique des plasmas à haute densité d'énergie appliquée à la production d'énergie par fusion. Ces travaux se sont concrétisés par l'écriture du code CHIC (Code d'Hydrodynamique et d'Implosion du CELIA), logiciel permettant de concevoir et de restituer des expériences dans le domaine de la Fusion par Confinement Inertiel (FCI). Le modèle théorique numérique décrivant l'implosion d'une cible laser est un système d'équations aux dérivées partielles au centre duquel on trouve les équations d'Euler écrites dans le formalisme lagrangien, couplées à des équations de diffusion non linéaires modélisant le transport de l'énergie par les électrons et les photons. Dans cet exposé, après un bref rappel du contexte physique, nous décrirons les deux méthodes originales qui constituent l'ossature numérique du code CHIC. Il s'agit de deux schémas numériques d'ordre élevé du type volumes finis dédiés respectivement à la résolution des équations de l'hydrodynamique lagrangienne et à la résolution d'équations de diffusion anisotrope sur des maillages bi-dimensionnels non-structurés. Le premier schéma, dénommé EUCCLHYD (Explicit Unstructured Lagrangian HYDrodynamics), permet de résoudre les équations de la dynamique des gaz sur un maillage mobile qui se déplace à la vitesse du fluide. Il est obtenu à partir d'un formalisme général basé sur le concept de forces de sous-mailles. Dans ce cadre, les flux numériques sont exprimés en fonction des forces de sous-mailles et de la vitesse des noeuds. Leur détermination repose sur les trois principes fondamentaux suivants : compatibilité géométrique entre le mouvement des noeuds et la variation de volume des mailles (loi de conservation géométrique), compatibilité avec le second principe de la thermodynamique et conservation de l'énergie totale et de la quantité de mouvement. L'extension de ce schéma à l'ordre deux est mise en place à l'aide d'une méthode basée sur la résolution d'un problème de Riemann généralisé dans l'approximation acoustique. Le second schéma, appelé CCLAD (Cell-Centered LAgrangian Diffusion), concerne la résolution de l'équation de la chaleur anisotrope non-linéaire. La discrétisation correspondante s'appuye sur une formulation variationnelle locale au niveau des sous-mailles qui permet de construire une approximation multi-points du flux de chaleur. Cette discrétisation d'ordre élevé rend possible la résolution des équations de la diffusion anisotrope avec une précision satisfaisante sur des maillages lagrangiens fortement déformés. La précision et la robustesse de ces méthodes numériques sont démontrées sur des cas-tests représentatifs.