Le présent travail est consacré à l’approximation numérique des solutions du modèle aux moments M1 pour le transfert radiatif. Il s’agit, ici, de développer des solveurs numériques performants et précis capables de prédire avec précision et robustesse des écoulements où le transfert radiatif joue un rôle essentiel. Dans ce sens, plusieurs méthodes numériques ont été envisagées pour la dérivation des schémas numériques de type solveur de Godunov. Une attention particulière a été portée sur les solveurs préservant les ondes de contact stationnaires. En particulier, un schéma de relaxation et un solveur HLLC sont présentés dans ce travail. Pour chacun de ces solveurs, la robustesse de la méthode a été établie (positivité de l’énergie radiative et limitation du flux radiatif). La validation et l’intérêt des méthodes abordées sont exhibés à travers de nombreuses expériences numériques mono et multidimensionelles. / The present work is dedicated to the numerical approximation of the M1 moments model solutions for radiative transfer. The objective is to develop efficient and accurate numerical solvers, able to provide with precise and robust computations of flows where radiative transfer effects are important. With this aim, several numerical methods have been considered in order to derive numerical schemes based on Godunov type solvers. A particular attention has been paid to solvers preserving the stationary contact waves. Namely, a relaxation scheme and a HLLC solver are presented in this thesis. The robustness of each of these solvers has been established (radiative energy positivity and radiative flux limitation). Several numerical experiments in one and two space dimensions validate the developed methods and outline their interest.

Dans ce mémoire, on s'intéresse à la simulation des écoulements liquide-vapeur en transition de phase. Pour décrire ces écoulements, une approche bifluide moyennée à deux pressions indépendantes est retenue. Cette description du mélange liquide-vapeur s'appuie sur le modèle à sept équations de Baer et Nunziato. On étudie les aptitudes de cette modélisation à simuler les transitions de phase apparaissant en ingénierie nucléaire.<br /><br />Dans un premier temps, on élabore un cadre thermodynamique théorique pour décrire les écoulements liquide-vapeur. Dans ce cadre, on réalise la fermeture du modèle de Baer et Nunziato. De nouvelles modélisations sont proposées pour les termes d'interaction entre les phases. Ces nouvelles modélisations dotent le modèle bifluide à deux pressions d'une inégalité d'entropie. On étudie ensuite les propriétés mathématiques de ce modèle. Sa partie convective hyperbolique se présente sous une forme non-conservative. On étudie tout d'abord la définition de ses solutions faibles. Divers régimes d'écoulement sont alors mis à jour pour le mélange diphasique. Ces différents régimes d'écoulement présentent des analogies avec le comportement fluvial et torrentiel des écoulements en rivière. Les stabilités linéaire et non-linéaire de l'équilibre liquide-vapeur sont ensuite établies. Pour affiner notre description des interactions diphasiques, on étudie pour finir l'implémentation d'un modèle de turbulence, ainsi que l'implémentation d'une procédure de reconstruction pour la densité d'aire interfaciale.<br /><br />On s'intéresse ensuite à la simulation de ce modèle. Suivant une approche à pas fractionnaires, une méthode numérique est élaborée dans un formalisme Volumes Finis. Pour réaliser l'approximation de la partie convective, diverses adaptations non-conservatives de solveurs de Riemann standard sont tout d'abord proposées. A l'inverse du cadre non-conservatif classique, l'ensemble de ces schémas converge vers une unique solution. Un nouveau schéma de relaxation est ensuite proposé pour approcher la dynamique des transferts interfaciaux. L'ensemble de la méthode numérique se caractérise alors par la préservation des équilibres liquide-vapeur. Dans un premier temps, cette méthode numérique est employée à la comparaison des différentes modélisations bifluides à une et deux pressions. On l'applique ensuite à la simulation des écoulements liquide-vapeur dans les circuits hydrauliques des réacteurs à eau sous pression en configuration accidentelle.

Ce travail traite de la modélisation des plasmas de fusion qui est ici abordée à l'aide d'un modèle Euler bi-températures et du modèle de la magnétohydrodynamique (MHD) idéale et résistive. Ces modèles sont tout d'abord établis à partir des équations de la MHD bi-fluide et nous montrons qu'ils correspondent à des régimes asymptotiques différents pour des plasmas faiblement ou fortement magnétisés. Nous décrivons ensuite les méthodes de volumes finis pour des maillages structurés et non-structurés qui ont été utilisées pour approcher les solutions de ces modèles. Pour les trois modèles mathématiques étudiés dans cette thèse, les méthodes numériques reposent sur des schémas de relaxation. Afin d'appliquer ces méthodes aux problèmes de fusion par confinement magnétique, nous décrivons comment modifier les méthodes de volumes finis pour les appliquer à des problèmes formulés en coordonnées cylindriques tout en gardant une formulation conservative forte des équations. Enfin nous étudions une stratégie pour maintenir la contrainte de divergence nulle du champ magnétique qui apparait dans les modèles MHD. Une série de cas tests numériques pour les trois modèles est présentée pour différentes géométries afin de valider les méthodes numériques proposées. / This work deals with the modeling of fusion plasma which is discussed by using a bi-temperature Euler model and the ideal and resistive magnetohydrodynamic (MHD) ones. First, these models are established from the bi-fluid MHD equations and we show that they correspond to different asymptotic regimes for lowly or highly magnetized plasma. Next, we describe the finite volume methods for structured and non-structured meshes which have been used to approximate the solution of these models. For the three mathematical models studied in this thesis, the numerical methods are based on relaxation schemes. In order to apply those methods to magnetic confinement fusion problems, we explain how to modify the finite volume methods to apply it to problem given in cylindrical coordinates while keeping a strong conservative formulation. Finally, a strategy dealing with the divergence-free constraint of the magnetic field is studied. A set of numerical tests for the three models is presented for different geometries to validate the proposed numerical methods.