Nous nous intéressons à l'analyse mathématique et aux simulations numériques d'un modèle de fluides non newtoniens. Nous couplons le modèle d'Hébraud et Lequeux décrivant l'écoulement de Couette, plan de suspensions concentrés à l'échelle mésoscopique avec l'équation de conservation de la quantité de mouvement. Outre le couplage multi-échelle, la difficulté principale de ce modèle vient du fait que l'équation mésoscopique est une équation parabolique de type Fokker-Planck non linéaire qui peut dégénérer en une équation hyperbolique. Le but de cette thèse est d'étudier des deux points de vue théorique et numérique le modèle obtenu. Ainsi, dans la première partie, nous démontrons les théorèmes d'existence et d'unicité de solutions dans des espaces fonctionnels appropriés, et dans la seconde partie, nous développons des schémas numériques pour approcher les solutions du problème. Deux méthodes (déterministe et stochastique) ont été implémentées. Des tests de réduction de variance relatifs à la méthode stochastique ont été réalisés. Finalement, le problème de calibrage des paramètres par la méthode de l'adjoint a été abordé.

La production industrielle d'aluminium met en jeu plusieurs aspects physiques, à la fois chimiques, thermiques et magnétohydrodynamiques (MHD). L'une de ses particularités est la coexistence dans une cuve de deux fluides non miscibles, ce qui conduit à la présence d'une interface libre. Ce procédé consomme près de 2% de l'électricité mondiale, la moitié étant perdue par effet Joule. L'enjeu est de réduire ce coût sans déstabiliser le procédé: il s'agit typiquement d'un problème de contrôle optimal, que nous traitons en considérant une modélisation MHD de la cuve. Deux approches sont utilisées pour effectuer cette optimisation, à savoir considérer une contrainte d'état non linéaire basée sur un couplage entre les équations de Maxwell et de Navier-Stokes multifluides, et une contrainte d'état linéaire résultant d'une approximation shallow water de la précédente. Après une courte introduction à la modélisation du procédé et aux concepts du contrôle optimal basé sur le principe de Pontryagin, nous décrivons dans un premier temps le contrôle de l'évolution de l'interface modélisée par l'approximation shallow water. S'ensuivent un travail de parallèlisation du logiciel de simulation du procédé dans le cadre non linéaire et la recherche numérique d'actionneurs acceptables pour son contrôle. Enfin, un algorithme d'optimisation de la forme de l'interface est proposé sous une contrainte d'état non linéaire simplifié, à savoir les équations de Navier-Stokes bifluides en dimension deux.