On s'intéresse dans ce travail à la simulation des écoulements diphasiques. Différents modèles, tous hyperboliques, sont considérés suivant les configurations étudiées. Dans un premier temps, plusieurs schémas Volumes Finis sont comparés pour l'approximation du modèle HEM (Homogeneous Equilibrium Model), notamment en présence de faibles densités. Ensuite on démontre l'existence et l'unicité de la solution faible entropique d'une loi de conservation scalaire gouvernant l'évolution de la saturation d'un écoulement diphasique dans un milieu poreux. On propose alors deux schémas Volumes Finis tenant compte du caractère résonnant de cette équation. La troisième partie concerne les écoulements en eaux peu profondes et l'approximation des termes sources raides. Une méthode permettant le maintien d'états au repos ainsi que le recouvrement et l'apparition de zones sèches, est présentée et comparée aux méthodes habituellement utilisées dans l'industrie. Enfin, une classe de modèles hyperboliques non conservatifs se basant sur l'approche bifluide à deux vitesses et deux pressions est proposée. Une étude des solutions discontinues du système convectif permet d'exhiber une classe de fermetures sur la vitesse interfaciale et sur la pression interfaciale, tout en permettant de définir de manière unique les produits non conservatifs. L'approximation se fait à l'aide d'une méthode de splitting d'opérateur. On utilise deux schémas Volumes Finis, le schéma de Rusanov et le schéma de Godunov approché VFRoe-ncv pour l'étape de convection. Plusieurs cas tests sont présentés et commentés : tubes à choc, conditions limites de paroi, robinet d'eau, sédimentation.

Cette thèse concerne le traitement numérique des interfaces pour des probèmes de propagation d'ondes dans des fluides parfaits et des solides élastique isotropes. La présence d'interfaces induit trois types de difficultés. Numériquement, on observe une chute de l'ordre de convergence des schémas, et des problèmes de stabilité. Géométriquement, la représentation en « marches d'escalier » conduit à des diffractions parasites. Physiquement, les schémas ne décrivent pas la nature des contacts. Nous résolvons ces trois problèmes via une méthode d'interface (schéma utilisé aux points de calcul voisins des interfaces, imposant à la solution numérique le respect des différentes conditions vérifiées par la solution exacte). L'étude est divisée en trois parties. La première partie débute par un état de l'art sur les méthodes d'interface classiques, comme l' « Immersed Interface Method » IIM), et sur leurs limitations. Les équations de l'acoustique et de l'élastodynamique sont écrites sous forme de systèmes hyperboliques du premier ordre. Différents schémas numériques, de complexité et qualité croissantes, sont présentés (Lax-Wendroff, volumes finis à limiteur de flux, WENO). La deuxième partie commence par le calcul des conditions de saut vérifiées par la solution exacte et par ses dérivées spatiales successives, pour différentes interfaces : fluide-fluide, fluide-solide et solide-solide, en contact parfait ou imparfait (conditions de masse-ressort). Nous proposons alors une nouvelle méthode d'interface, l' « Explicit Simplified Interface Method ». Cette méthode conserve, en présence d'interfaces, des propriétés des schémas en miieu homogène. L'ESIM est de mise en œuvre aisée, s'adapte aux différents schémas, et conduit à un surcoût informatique négligeable. La troisième partie concerne la validtion numérique de la méthode. Les comparaisons de solutions analytiques et de solutions calculées avec l'ESIM permettent de vérifier les propriétés attendues.

On s'intéresse dans ce travail à la simulation et l'approximation d'écoulements diphasiques avec transition de phase. Il s'agit de modéliser la dynamique d'une bulle de cavitation. Le modèle repose sur les équations d'Euler en coordonnées sphériques, l'interface gaz-liquide étant indiquée par une fonction de couleur. Dans une première partie, aucun transfert de masse n'est supposé. Un schéma Lagrange-projection est d'abord proposé : seule l'interface est déplacée à la vitesse du fluide. La projection repose sur un algorithme de suppression-découpage qui assure que les volumes ne deviennent pas négatifs. Le second chapitre traite du terme source géométrique. On construit un schéma équilibre articulé autour du schéma VFRoe-ncv pour lequel on propose une correction entropique non paramétrique. Un méthode d'ordre élevé de type Galerkin discontinu est ensuite étudiée dans le cadre de la magnétohydrodynamique. L'intégration en temps est réalisée par une méthode de type Adams-Bashforth, qui s'avère bien adaptée aux algorithmes de pas de temps local. La deuxième partie de la thèse traite de la modélisation du changement de phase liquide-vapeur. L'inf-convolution et la transformée de Legendre définissent une structure naturelle pour la construction de loi de pression de mélange. En particulier on montre que la construction de Maxwell est équivalente à la construction de l'enveloppe convexe de l'énergie de van der Waals. Un algorithme de transformée de Legendre rapide, adapté à la prise en compte correcte des bords du domaine de calcul de loi d'état, est développé. La méthode est appliquée à la construction de lois de pression tabulées de mélanges binaires miscibles ou immiscibles de gaz raides.

L'objet de cette thèse est d'étudier l'applicabilité de la méthode d'approximation rationnelle en norme de Hankel à des systèmes dynamiques linéaires de dimension d'état infinie. On illustre par trois exemples concrets les possibilités d'utilisation des techniques d'approximation développées ces dernières années, notamment par Curtain, Glover et Partington. Les exemples choisis représentent des phénomènes d'évolution décrits par des équations aux dérivées partielles, par rapport au temps et aux variables d'espace. Il s'agit: d'un problème de diffusion de chaleur, de type parabolique, pour lequel les techniques d'approximation s'adaptent assez directement ; de deux problèmes hyperboliques décrivant l'évolution d'une poutre en flexion et en torsion, pour lesquels une méthode originale appelée ``relaxation'' a été mise au point: préalable à l'approximation de Hankel, elle permet son application lorsque les pôes associés au système hyperbolique croissent suffisamment rapidement.

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.

La Commande prédictive ou Commande Optimale à Horizon Glissant (COHG) devient de plus en plus populaire dans de nombreuses applications pratiques en raison de ses avantages importants tels que la stabilisation et la prise en compte des contraintes. Elle a été bien étudiée pour des systèmes en dimension finie même dans le cas non linéaire. Cependant, son extension aux systèmes en dimension infinie n'a pas retenu beaucoup d'attention de la part des chercheurs. Ce travail de thèse apporte des contributions à l'application de cette approche aux systèmes de lois de conservation. Nous présentons tout d'abord une preuve de stabilité complète de la COHG pour certaines classes de systèmes en dimension infinie. Ce résultat est ensuite utilisé pour les systèmes hyperboliques 2x2 commandés aux frontières et appliqué à un problème de contrôle de canal d'irrigation. Nous proposons aussi l'extension de cette stratégie au cas de réseaux de systèmes hyperboliques 2x2 en cascade avec une application à un ensemble de canaux d'irrigation connectés. Nous étudions également les avantages de la COHG dans le contexte des systèmes non linéaires et semi-linéaires notamment vis-à-vis des chocs. Toutes les analyses théoriques sont validées par simulation afin d'illustrer l'efficacité de l'approche proposée.

Les systèmes hyperboliques interviennent dans de nombreuses branches des sciences : théorie cinétique, mécanique des fluides non visqueux, magnéto hydrodynamique, dynamique des gaz non visqueux, trafic routier, flux d’une rivière ou d’un glacier, processus de sédimentation, processus d’échanges chimiques, etc. Et souvent, les systèmes qui régissent ces évènements font intervenir des petits paramètres, dont l’étude asymptotique permet d’envisager des simplifications mathématiques et/ou informatiques notoires. L’existence locale et l’existence globale de solutions, uniformément par rapport à ces paramètres, sont des questions fondamentales. Cette thèse regroupe à la fois des résultats généraux sur l’existence locale uniforme de solutions pour des systèmes hyperboliques quasi-linéaires du premier ordre ; et sur l’existence globale uniforme de solutions autour d’un équilibre constant pour ces mêmes systèmes. Le cas du système d’Euler-Maxwell ne satisfaisant pas les conditions requises pour l’existence uniforme globale, nous le traitons à part. / Hyperbolic systems arise in a large field of sciences : kinetic theory, inviscid reactive flow, magnetohydrodynamics, inviscid gas dynamics, traffic flow, river or glacier flow, sedimentation processes, chemical exchange processes, etc. In these kind of systems, small paramaters often appear, and an asymptotic study may lead to mathematical or computational simplifications. One fundamental problem that we may work on is local and global existence of solutions for these systems, uniformly with respect to these parameters. This Ph.D. thesis includes, on one hand, general results on uniform local existence of solutions for first order quasi-linear hyperbolic systems ; and on the other hand, results on uniform global existence of solutions near constant equilibriums for these same systems. In the case of Euler-Maxwell systems, required conditions are not fulfilled for uniform global existence, then we treat it separately.

Cette thèse est essentiellement composée de deux parties. Dans la première partie, on étudie le système d'Euler-Maxwell. En utilisant la méthode d'intégration de l'énergie classique, on montre l'existence et l'unicité de solutions régulières du système avec données initiales petites. Ensuite, on étudie la limite de relaxation en montrant que, le sytème d'Euler-Maxwell converge vers les équations de dérive-diffusion quand le temps de relaxation tend vers zéro. Dans la deuxième partie, on cherche la contrôlabilité et l'observabilité exactes frontières de systèmes hyperboliques quasi-linéaires dans un réseau du type d'arbre. On établit des résultats d'existences de la contrôlabilité et l'observabilité par des méthodes constructives qui sont basées sur la théorie de la solution C1 semi-globale du système hyperbolique quasi-linéaire du premier ordre avec conditions initiales et frontières. Ensuite, on trouve des dualités de la contrôlabilité et l'observabilité.

Cette thèse concerne l'étude de modèles mathématiques et méthodes numériques motivés par la filtration des globules blancs du sang. <br /><br />Dans la première partie, nous définissons des modèles mathématiques qui réprésentent les principaux phénomènes physiques qui entrent en jeu dans le procédé de la filtration.<br /><br />La deuxième partie est dédiée à l'analyse mathématique de systèmes d'équations aux dérivées partielles modélisant le procédé de la filtration. Tout d'abord, nous considérons un système d'équations semi-linéaires de type hyperbolique-parabolique avec une diffusion anisotrope dégénérée. Nous étudions ce problème avec une théorie $L^{1}$; nous considérons en particulier l'existence et l'unicité de solutions faibles ainsi que d'autres propriétés comme le principe du maximum; puis nous établissons la limite quand la constante de réaction devient grande. Nous montrons que le système converge vers une équation non linéaire parabolique-hyperbolique qui généralise le problème de Stefan. Nous étudions également, par des techniques de l'homogénéisation, la filtration au travers de milieux poreux fibrés. Le réseau des fibres étudié est celui utilisé par M. Briane dans le cadre d'une étude sur la conduction thermique des tissus biologiques. Nous dérivons et justifions l'équation de Darcy ainsi que la forme du tenseur de perméabilité pour un tel milieu fibreux. Les résultats théoriques concernant la perméabilité sont illustrés par quelques simulations numériques. Finalement, nous considérons le cas où le diamètre des fibres tend vers zéro. En appliquant des résultats de G. Allaire à notre cas, nous justifions rigoureusement la forme du terme dominant dans les formules de perméabilité efficace utilisées en ingénierie. Ces résultats sont également confirmés par un calcul numérique direct de la perméabilité, dans lequel la petitesse du diamètre des fibres rend nécessaire le recours à des approximations de précision élevée.<br /><br />La définition des méthodes numériques efficaces pour approximer la solution des modèles mathématiques est envisagée dans la troisième partie. Précisément, concernant les équations de Darcy, nous avons utilisé la méthode des éléments finis mixtes hybrides. Pour la résolution de l'équation du transport, nous avons implémenté une méthode numérique utilisant des volumes finis pour la discrétisation du terme convection/réaction associé à une approximation mixte hybride pour la discrétisation du terme dispersif.

L'objectif global de ce travail est d'étudier une modélisation multi-échelle et de développer une méthode adaptée pour la simulation numérique du processus d'érosion à l'échelle du bassin versant. Après avoir passé en revue les différents modèles existants, nous dérivons une solution analytique non triviale pour le système couplé modélisant le transport de sédiments par charriage. Ensuite, nous étudions l'hyperbolicité de ce système avec diverses lois de sédimentation proposées dans la littérature. Concernant le schéma numérique, nous présentons le domaine de validité de la méthode de splitting, pour les équations modélisant l'écoulement et celle décrivant l'évolution du fond. Pour la modélisation du transport en suspension à l'échelle de la parcelle, nous présentons un système d'équations couplant les mécanismes d'infiltration, de ruissellement et le transport de plusieurs classes de sédiments. L'implémentation et des tests de validation d'un schéma d'ordre élevé et de volumes finis bien équilibré sont également présentés. Ensuite, nous discutons sur l'application et la calibration du modèle avec des données expérimentales sur dix parcelles 1m2 au Niger. Dans le but d'aboutir la simulation à l'échelle du bassin versant, nous développons une modélisation multi échelle dans laquelle nous intégrons le taux d'inondation dans les équations d'évolution afin de prendre en compte l'effet à petite échelle de la microtopographie. Au niveau numérique, nous étudions deux schémas bien équilibrés : le schéma de Roe basé sur un chemin conservatif, et le schéma avec reconstruction hydrostatique généralisée. Enfin, nous présentons une première application du modèle avec les données expérimentales du bassin versant de Ganspoel qui nécessite la parallélisation du code.