Ce manuscrit est dédié à l'approximation numérique de plusieurs modèles du transfert radiatif. Dans un premier temps, l'attention est portée sur le modèle cinétique d'ordonnées discrètes. Dans le but de coupler ce modèle avec d'autres phénomènes plus lents, il est nécessaire d'avoir des méthodes numériques performantes et précises sur des temps longs. À partir d'une double approximation polynomiale de la solution en temps et en espace, on développe un schéma de type GRP d'ordre élevé sans restriction sur le pas de temps pour un système hyperbolique linéaire sur des maillages non structurés. Ce schéma est ensuite étendu pour le modèle d'ordonnées discrètes. Dans un second temps, on s'intéresse à des modèles aux moments issus du transfert radiatif. En effet, dans certaines applications, les modèles aux moments de type M1 conservent de nombreuses propriétés de l'ETR et fournissent une approximation suffisante de la solution. Après avoir résolu le problème de Riemann associé au modèle M1 gris, on considère l'approximation numérique du modèle M1 multigroupe. Une attention particulière est portée sur le calcul des moyennes d'opacités et des lois de fermeture. Un algorithme de précalculs est alors mis en place. La dernière application traitée dans ce mémoire porte sur une extension du transfert radiatif pour estimer des doses de radiothérapie. À la différence du M1 gris usuel, les flux dépendent ici de fonctions peu régulières en espace. Grâce à des changements de variables, un schéma HLL rétrograde est développé. De nombreux exemples numériques illustrent l'intérêt des schémas obtenus dans cette étude.

Cette thèse étudie essentiellement des phénomènes touchant aux domaines aéronautique et spatial. Elle traite par simulations numériques des écoulements dans les tuyères de moteur-fusée, les écoulements d'arrière-corps d'avions, et les écoulements supersoniques en présence d'interactions onde de choc/couche limite. La motivation principale de ce travail a été de comprendre les différents facteurs qui gouvernent ces écoulements. Ce travail, soutenu par le groupe SNECMA, est consacré plus particulièrement à l'étude et à la modélisation des phénomènes en proche paroi en utilisant le code industriel N3S-Natur. L'objectif visé est d'améliorer les chaînes de conception et de développement industriels (RANS) et de développer de nouvelles méthodes de calcul (LES/DES) pour une meilleure maîtrise des systèmes énergétiques. Ainsi, les modèles mis au point dans cette étude ont été d'abord validés sur une variété d'écoulement simple (couche limite compressible, jet supersonique,...) avant d'être appliqués aux écoulements complexes (tuyères avec film de refroidissement, arrière-corps 3D, interaction visqueuse et décollement tridimensionnel). Une partie de l‘étude a été consacrée aux phénomènes transitoires d'amorçage rapide des tuyères propulsives. En ce qui concerne les aspects instationnaires, il a été montré, à travers l'étude d'une interaction onde de choc/couche limite, que la DES (Detached Eddy Simulation), tout comme la LES (Large Eddy Simulation), constitue un outil pertinent, permettant une prédiction fine des caractéristiques moyennes et fluctuantes des écoulements supersoniques décollés. Les résultats obtenus dans ce travail confirment la portée et l'intérêt scientifiques des études en aérodynamique supersonique.

Dans ce travail nous nous intéressons à la résolution du problème d'évolution de Navier-Stokes incompressible totalement discrétisé en temps et en espace, en dimension deux par une méthode à deux grilles. Dans un premier temps, nous étendons la méthode à deux grilles, appliquée par V. Girault et J.-L. Lions au problème de Navier-Stokes instationnaire semi-discrétisé au problème de Navier-Stokes totalement discrétisé en temps (par un schéma d'ordre un) et en espace (par une méthode d'éléments finis d'ordre un). Dans la première étape, le problème non-linéaire est discrétisé en espace et en temps sur une grille grossière de pas d'espace H avec un pas de temps Delta t. Puis dans la deuxième étape, le problème, linéarisé autour de la vitesse u_H calculée à l'étape précédente, est discrétisé en espace sur une grille fine de pas d'espace h et le même pas de temps. L'idée de la méthode à deux grilles est que, sous des hypothèses adéquates, la contribution de u_H à l'erreur dans le terme non-linéaire en espace, est mesurée en norme L^2 en espace et en temps et a un ordre plus élevé que si elle était mesurée en norme H^1. Dans un deuxième temps, vu que le but est de gagner en ordre de convergence de l'erreur totale du schéma ainsi qu'en complexité nous étudions un schéma à deux grilles d'ordre deux en temps du problème totalement discrétisé en temps et en espace de Navier-Stokes. Nous présentons les résultats suivants: dans le cas de la résolution du schéma d'ordre un en temps, si h = H^2 = Delta t, alors l'erreur globale de l'algorithme à deux grilles est de l'ordre de h. Dans le cas du schéma d'ordre deux en temps, si h^2 = (Delta t)^2 = H^3, alors l'erreur globale de l'algorithme à deux grilles est de l'ordre de h^2: résultats identiques à ceux de la résolution directe du problème non-linéaire sur une grille fine.

Cette thèse traite des équations de Maxwell en domaine temporel. Le principal objectif est de proposer des méthodes de type éléments finis d'ordre élevé pour les équations de Maxwell et des schémas d'intégration en temps efficaces sur des maillages localement raffinés. Nous considérons des méthodes GDDT (Galerkine Discontinues en Domaine Temporel) s'appuyant sur une interpolation polynomiale d'ordre arbitrairement élevé des composantes du champ électromagnétique. Les méthodes GDDT pour les équations de Maxwell s'appuient le plus souvent sur des schémas d'intégration en temps explicites dont la condition de stabilité peut être très restrictive pour des maillages raffinés. Pour surmonter cette limitation, nous considérons des schémas en temps qui consistent à appliquer un schéma implicite localement, dans les régions raffinées, tout en préservant un schéma explicite sur le reste du maillage. Nous présentons une étude théorique complète et une comparaison de deux méthodes GDDT localement implicites. Des expériences numériques en 2D et 3D illustrent l'utilité des schémas proposés. Le traitement numérique de milieux de propagation complexes est également l'un des objectifs. Nous considérons l'interaction des ondes électromagnétiques avec les tissus biologiques qui est au cœur de nombreuses applications dans le domaine biomédical. La modélisation numérique nécessite alors de résoudre le système de Maxwell avec des modèles appropriés de dispersion. Nous formulons une méthode GDDT localement implicite pour le modèle de Debye et proposons une analyse théorique et numérique complète du schéma.

Les sujets abordés dans cette thèse s'articulent autour de la modélisation numérique du couplage entre l'hydrodynamique et la biologie pour la culture industrielle de microalgues dans des raceways. Ceci est fait au moyen d'un modèle multicouches qui disrétise verticalement les équations de Navier-Stokes hydrostatiques couplé avec un modèle de Droop photosensible pour représenter la croissance des algues, notamment la production de carbone. D'un point de vue numérique, une méthode volumes finis avec schémas cinétiques est appliquée. Elle permet d'obtenir un schéma équilibre qui préserve la positivité de la hauteur d'eau et des quantités biologiques et qui satisfait une inégalité d'énergie. Des simulations sont effectuées en 2D et en 3D, au moyen d'un code C++ développé à cet effet. Du point de vue de l'intérêt pratique de ce travail, ces simulations ont permis de mettre en évidence l'utilité de la roue à aube présente dans les raceways, mais aussi d'exhiber les trajectoires lagrangiennes réalisées par les microalgues, qui permettent de connaitre l'historique lumineux des algues, information d'une grande importance pour les biologistes car elle leur permet d'adapter leurs modèles de croissance phytoplanctoniques à ce contexte très particulier et non naturel. Afin de valider les modèles et les stratégies numériques employées, deux pistes on été explorées. La première consiste à proposer des solutions analytiques pour les équations d'Euler à surface libre, ainsi qu'un modèle biologique spécifique permettant un couplage analytique. La deuxième consiste à faire de l'assimilation de données. Afin de tirer partie de la description cinétique des lois de conservation hyperboliques, une méthode innovante basée sur la construction d'un observateur de Luenberger au niveau cinétique est développée. Elle permet d'obtenir un cadre théorique intéressant pour les lois de conservation scalaires, pour lesquelles on étudie les cas d'observations complètes, partielles en temps, en espace, et bruitées. Pour les systèmes, on se concentre particulièrement sur le système de Saint-Venant, système hyperbolique non linéaire et un observateur basé sur l'observation des hauteurs d'eau uniquement est construit. Des simulations numériques dans les cas scalaires et systèmes, en 1D et 2D sont effectuées et valident l'efficacité de la méthode.

Cette thèse est consacrée à l'étude de certains aspects numériques et mathématiques liés à la simulation d'écoulements incompressibles triphasiques à l'aide d'un modèle à interfaces diffuses de type Cahn-Hilliard/Navier-Stokes. La discrétisation spatiale est effectuée par éléments finis. La présence d'échelles très différentes dans le système suggère l'utilisation d'une méthode de raffinement local adaptatif. La procédure mise en place permet de tenir compte implicitement des non conformités des maillages générés, pour produire in fine des espaces d'approximation conformes. Nous montrons, en outre, qu'il est possible d'exploiter cette méthode pour construire des préconditionneurs multigrilles. Concernant la discrétisation en temps, notre étude a commencé par celle du système de Cahn-Hilliard. Pour remédier aux problèmes de convergence de la méthode de Newton utilisée pour résoudre ce système (non linéaire), nous proposons un schéma semi-implicite permettant de garantir la décroissance de l'énergie. Nous montrons l'existence et la convergence des solutions discrètes. Nous poursuivons ensuite cette étude en donnant une discrétisation en temps inconditionnellement stable du modèle complet Cahn-Hilliard/Navier-Stokes ne couplant pas fortement les deux systèmes. Nous montrons l'existence des solutions discrètes et, dans le cas où les trois fluides ont la même densité, nous montrons leur convergence. Nous étudions, pour terminer cette partie, diverses problématiques liées à l'utilisation de la méthode de projection incrémentale. Enfin, la dernière partie présente plusieurs exemples de simulations numériques, diphasiques et triphasiques, en deux et trois dimensions.

Ce travail s'inscrit dans l'étude mathématique d'écoulements de fluides environnementaux. Nous en abordons deux aspects, à travers deux contextes distincts d'application.En lien avec la simulation des écoulements en milieux poreux, on s'intéresse dans une première partie à la discrétisation d'opérateurs de diffusion anisotropes hétérogènes par des méthodes de volumes finis sur des maillages généraux. Dans le but d'obtenir des solutions approchées qui respectent les bornes physiques des modèles, notre attention se porte sur la conservation du principe du maximum pour les opérateurs elliptiques. Nous présentons des mécanismes généraux permettant de corriger tout schéma volumes finis afin de garantir un principe du maximum discret tout en préservant certaines de ses propriétés principales. On étudie en particulier les propriétés de coercivité et de convergence des schémas corrigés.La deuxième partie est consacrée à la construction de modèles approchés pour la propagation des vagues en eaux peu profondes et sur des topographies irrégulières. A cet effet, nous proposons tout d'abord une adaptation de la démarche d'étude classique à des écoulements bidimensionnels sur des topographies polygonales. Dans un cadre plus général, nous développons ensuite une démarche formelle qui débouche sur des alternatives non locales à quelques modèles classiques (équations de Saint-Venant, équations de Serre, système de Boussinesq). Ces nouveaux modèles contiennent des termes régularisants pour les contributions du fond.

La Simulation aux Grandes Echelles (SGE) est de plus en plus utilisée dans les processus de développement et la conception des réacteurs aéronautiques industriels. L'une des raisons pour ce besoin résulte dans la capacité de la SGE à fournir des informations instantanées d'un écoulement turbulent augmentant la quantité des prédictions de la composition des gaz d'échappement. Ce manuscrit de thèse aborde deux sujets récurrents de la SGE. D'une part, les schémas numériques pour la SGE nécessitent certaines propriétés, notamment une précision élevée avec une diffusivité faible pour ne pas nuire aux modèles de turbulence. Afin de répondre à ce pré requis, une famille de schémas d'intégration temporelle d'ordre élevée est proposée, permettant de modifier la diffusion numérique du schéma. D'autre part, la SGE étant intrinsèquement instationnaire, elle est très consommatrice en temps CPU. De plus, une géométrie complexe prend beaucoup de temps de simulation même avec les super calculateurs d'aujourd'hui. Dans le cas particulier d'intérêt et souvent rencontré dans les applications industrielles, l'approche bas-Mach est constitue une alternative intéressante permettant de réduire le coût et le temps de retour d'une simulation LES. L'impact et la comparaison des formalismes compressible et incompressible sont toutefois rarement quantifiés, ce qui est proposé dans ce travail pour une configuration représentative d'un brûleur swirlé industriel mesuré au CORIA.

Cette recherche a pour objet la compréhension des énoncés de problèmes arithmétiques complexes et leur résolution. Les problèmes complexes choisis combinent des problèmes simples de types Changement et Combinaison. Ce travail s'appuie sur le modèle de la Construction-Intégration de Kintsch. Les résultats montrent qu'il existe une relation entre le niveau d'expertise en compréhension de textes narratifs et la résolution des problèmes arithmétiques complexes. Comprendre un texte narratif ou un énoncé de problème complexe exige de la part des lecteurs la construction d'un réseau propositionnel hiérarchisé et les résultats suggèrent, entre autres, une sensibilité des élèves aux propositions textuelles et aux ellipses contenues dans les textes. La formation des macropropositions est un processus fondamental et les résultats montrent une relation entre le nombre d'objets contenus dans les énoncés de problème et la procédure préférentiellement choisie par les élèves. Ils suggèrent d'une part, la mise en oeuvre du processus de catégorisation au cours du processus de compréhension et d'autre part, l'affaiblissement des liaisons entre les macropropositions élaborées et le schéma de problème Parties-Tout qui leur sont liés. D'un point de vue pédagogique, les résultats montrent que les questions relatives à l'activation d'une part des concepts superordonnés et d'autre part des schémas de problèmes Parties-Tout ne sont pas à privilégier pour aider les élèves. Finalement, les connaissances du lecteur sont essentielles à la compréhension. Cet élément est confirmé ici et la compréhension des problèmes complexes nécessite des connaissances solides relativement aux problèmes arithmétiques simples.

Cette thèse traite des équations de Maxwell en domaine temporel. Le principal objectif est de proposer des méthodes de type éléments finis d'ordre élevé pour les équations de Maxwell et des schémas d'intégration en temps efficaces sur des maillages localement raffinés. Nous considérons des méthodes GDDT (Galerkine Discontinues en Domaine Temporel) s'appuyant sur une interpolation polynomiale d'ordre arbitrairement élevé des composantes du champ électromagnétique. Les méthodes GDDT pour les équations de Maxwell s'appuient le plus souvent sur des schémas d'intégration en temps explicites dont la condition de stabilité peut être très restrictive pour des maillages raffinés. Pour surmonter cette limitation, nous considérons des schémas en temps qui consistent à appliquer un schéma implicite localement, dans les régions raffinées, tout en préservant un schéma explicite sur le reste du maillage. Nous présentons une étude théorique complète et une comparaison de deux méthodes GDDT localement implicites. Des expériences numériques en 2D et 3D illustrent l'utilité des schémas proposés. Le traitement numérique de milieux de propagation complexes est également l'un des objectifs. Nous considérons l'interaction des ondes électromagnétiques avec les tissus biologiques qui est au cœur de nombreuses applications dans le domaine biomédical. La modélisation numérique nécessite alors de résoudre le système de Maxwell avec des modèles appropriés de dispersion. Nous formulons une méthode GDDT localement implicite pour le modèle de Debye et proposons une analyse théorique et numérique complète du schéma.