Nous considérons les fibrés à connexion non-singulière ou méromorphe, de rang 2 et sans trace sur les surfaces de Riemann compactes de genre quelconque. <br />En déformant la courbe, la position des pôles et la connexion, nous construisons la déformation isomonodromique universelle d'un tel fibré à connexion. Notre construction spécifique au cas du rang 2 et sans trace est plus élémentaire que la construction en rang quelconque due à B. Malgrange et I. Krichever au sens où elle ne nécessite pas d'analyse de Stokes des singularités irrégulières. De plus, elle englobe le cas des singularités résonantes de manière naturelle.<br />Nous montrons que le fibré vectoriel sous-jacent à la déformation isomonodromique universelle est génériquement 'maximalement' stable, pourvu que le fibré à connexion initial soit irréductible. À cette fin, nous démontrons une version analytique du résultat de semicontinuité de M. Maruyama, puis nous nous ramenons à un problème de transversalité de feuilletages. À l'aide d'exemples explicites, nous montrons que la condition d'irréductibilité est nécessaire et que l'ensemble analytique des paramètres non génériques au sens ci-dessus peut être non algébrique.

Ce mémoire est consacré à la résolution du problème de Plateau à bord polygonal dans l'espace euclidien et dans l'espace de Minkowski de dimension trois. Il s'appuie sur la méthode de résolution proposée par René Garnier dans le cas euclidien dans un article publié en 1928 et qui a été oublié depuis, voire ignoré à l'époque. Plus géométrique et constructive que la méthode variationnelle, l'approche de Garnier est cependant parfois très compliquée, voire obscure et incomplète. On retranscrit sa démonstration dans un formalisme moderne, tout en proposant de nouvelles preuves plus simples, et en en complétant certaines lacunes. Ce travail repose principalement sur l'utilisation plus systématique des systèmes fuchsiens et la mise en évidence du lien entre la réalité de ces systèmes et leur monodromie. Ceci nous permet d'étendre le résultat de Garnier dans l'espace de Minkowski. La méthode de Garnier repose sur le fait que, par la représentation de Weierstrass spinorielle des surfaces minimales, on peut associer une équation fuchsienne réelle du second ordre définie sur la sphère de Riemann à tout disque minimal à bord polygonal. La monodromie de cette équation est déterminée par les directions orientées des côtés du bord. Pour résoudre le problème de Plateau, on est donc amené à résoudre un problème de Riemann-Hilbert. On procède ensuite en deux étapes : on construit d'abord, par déformations isomonodromiques, la famille de tous les disques minimaux dont le bord est un polygone de directions orientées données. Puis on montre, en étudiant les longueurs des côtés des bords polygonaux, qu'on obtient ainsi tout polygone comme bord d'un disque minimal.

Ces travaux portent sur la modélisation et la simulation numérique de la dispersion de matériaux nocifs (pulvérisations liquides ou gazeuses) en milieu urbain ou naturel (attentat ou explosion accidentelle survenant en zone peuplée, fuites de produits toxiques gazeux ou liquides, éclatement de réservoir..). Afin de prédire ces risques un outils de simulation tridimensionnel a été développé. Celui-ci est basé sur un modèle de milieu hétérogène afin de traiter des phénomènes dont la durée et les distances associées peuvent être très grandes. La topographie des milieux étudiées est prise en compte grâce à des données numériques d'´elévation ainsi que les conditions météo permettant l'utilisation de profils de température et de vent complexes. Les transferts de chaleur et de masse sont considérés, notamment au niveau des obstacles. Un schéma numérique d'ordre élevé en temps et en espace est utilisé pour calculer les concentrations massiques de polluants. Par ailleurs, un modèle d'écoulement gaz-particule a été développé et implémenté dans le code de calcul. L'instabilité d'une couche de fluide soumise à un important gradient de pression est également étudiée, ceci afin de mieux comprendre et de caractériser les conditions initiales à utiliser pour ce type d'écoulement, impliquant des couches de particules. / This work deals with the modeling and the numerical simulation of the dispersion of toxic cloud of dropplets or gas in uneven geometry such as urban environment, industrial plants and hilly environment. Examples of phenomena under study are the dispersion of chemical products from damaged vessels, gas diffusion in an urban environment under explosion conditions, shock wave propagation in urban environment etc. A 3D simulation code has been developed in this aim. To simplify the consideration of complex geometries, a heterogeneous discrete formulation has been developed. When dealing with large scale domains, such as hilly natural environment, the topography is reconstructed with the help of numerical elevation data. Meteorological conditions are also considered, concerning temperature and wind velocity profiles. Heat and mass transfers on subscale objects, such as buildings are studied. A high order numerical scheme in space and time is used to compute mass concentration of pollutant. A two-phase model for dilute gas-particles flow has been developed and implemented in the 3D simulation code. The instability of a fluid layer appearing under high pressure gradient is also studied. This analysis allows us a better understanding of initial conditions for similar problems involving particles layer.

Ce travail traite des méthodes de calcul numérique pour les simulations hydrodynamiques appliquées principalement sur des produits développés par ANDRITZ HYDRO. Il s’agit ici de mettre en place des schémas d’ordre élevé pour des simulations CFD en utilisant le code de calcul ASPHODEL développé et utilisé par ANDRITZ HYDRO. Les principales motivations sont l’augmentation de la fiabilité des résultats de calculs numériques avec un coût de calcul raisonnable. Cette fiabilité s’exprime à travers l’augmentation de la précision et de la robustesse des schémas numériques. Le code de calcul ASPHODEL est basé sur la méthode sans maillage SPH-ALE. Mélange entre les volumes finis et la méthode SPH (Smoothed Particle Hydrodynamics), la méthode SPH-ALE emploie un ensemble de points appelés particules servant à la discrétisation du domaine fluide. Elle permet en particulier de par son caractère sans maillage, un suivi des surfaces libres sans effort de calcul supplémentaire. Cet aspect est véritablement attrayant pour bon nombre d’applications industrielles notamment la simulation des écoulements à surface libre se produisant dans une turbine Pelton, mais également le remplissage d’une turbine Francis. Cependant, le bémol à cette méthode est son manque de précision spatiale. En effet les points de calcul étant mobiles, les opérateurs spatiaux doivent être en mesure de conserver leur précision et leur robustesse au cours du temps. La qualité des résultats en est du coup impactée, en particulier le champ de pression souvent excessivement bruité. La montée en ordre et l’amélioration de la consistance des opérateurs pour un vaste panel de configurations géométriques sont donc les enjeux de ce travail. En utilisant des outils inspirés par les volumes finis non-structurés, il est possible d’améliorer les opérateurs spatiaux. En effet, la montée en ordre ou p-raffinement peut notamment se faire avec des reconstructions d’ordres élevés pour évaluer les états aux interfaces des problèmes de Riemann. La sommation des flux numériques résolus par un solveur de Riemann est ensuite retravaillée pour obtenir un schéma numérique d’ordre global cohérent. Le même soucis de cohérence avec les schémas en temps doit d’ailleurs être pensé. Le gain de précision apporté par les schémas numériques d’ordre élevé est comparé avec un raffinement spatial, c’est à dire une augmentation du nombre des particules de taille plus petite, aussi appelé h-raffinement. La méthode SPH-ALE améliorée est ensuite testée sur des cas représentatifs des applications visées. En conclusion, les développements effectués dans cette étude ont été guidés par l’application en turbine Pelton principalement mais il va de soi qu’ils sont applicables à des écoulements sans surface libre dans les turbines Francis par exemple. Ce travail montre les possibilités d’une méthode sans maillage pour des cas d’écoulements complexes autour de géométrie tournantes. / This work deals with numerical methods for hydrodynamic testing applied mainly on products developed by ANDRITZ HYDRO. This is to put in place high order schemes for CFD simulations using the ASPHODEL calculation code developed and used by ANDRITZ HYDRO. The main reasons are the increased reliability of the results of numerical calculations with a reasonable computational cost. This reliability is expressed through increasing the accuracy and robustness of numerical schemes. The ASPHODEL computer code is based on the meshfree method SPH-ALE. Mix between finite volume method and SPH (Smoothed Particle Hydrodynamics), the SPH-ALE method uses a set of points called particles serving as the fluid domain discretization. It allows track free surfaces without additional computational effort. This is truly attractive for many industrial applications including the simulation of free surface flows occurring in a Pelton turbine, but also filling a Francis turbine. However, the downside of this method is its lack of spatial accuracy. Indeed calculation points are mobile, space operators must be able to keep their accuracy and robustness over time. The quality of results is impacted especially the pressure field is often excessively noisy. The rise in order and improving the consistency of the operators for a wide range of geometric configurations are the challenges of this work. Using tools inspired by the unstructured finite volume, it is possible to improve the spatial operators. Indeed, the increasing order or p-refinement particular can be done with reconstructions of high order to assess the conditions at the interfaces of Riemann problems. The summation of discret fluxes solved by Riemann solver is then reworked to obtain a coherent global order scheme. The same concern for consistency with temporal schemes should also be considered. The precision gain provided by numerical schemes of higher orders is compared with a spatial refinement ie an increase in the number of smaller particles ; also called h -refinement . Improved SPH -ALE method is then tested on representative cases of intended applications. In conclusion, the developments made in this study were guided in accordance mainly with the Pelton turbine but it goes without saying that they are applicable to non- free surface flows in Francis turbines for example. This work shows the possibilities of a free mesh method for cases of complex flow around rotating geometry.

Cette thèse est consacrée à l'analyse, à l'application et à l'extension bidimensionnelle, d'un nouveau schéma aux volumes finis (SRNH) proposé récemment pour une classe de système non homogène. L'analyse de stabilité du schéma, d'abord dans le cas scalaire ensuite dans le cas de systèmes, mène à une nouvelle formulation où intervient le signe de la matrice Jacobienne du système de lois de bilan considéré. Pour le système de Saint Venant avec terme de pente, on montre formellement que le schéma SRNHS vérifie la C-propriété exacte introduite pour les schémas équilibres par Bermùdez et Vázquez. Les résultats numériques 1D et 2D, en particulier du cas de rupture de barage sur un fond en forme de marche, montrent le degrés d'efficacité du schéma. Pour le système diphasiques des zones de non hyperbolicité peuvent exister, avec apparition de valeurs propres complexes dans la Jacobienne du système. On montre que pour les configurations faiblement non hyperboliques, on peut calculer le signe de la Jacobienne par l'algorithme de Newton-Schultz. Pour les configurations plus raides, où la méthode précédente ne fonctionne plus, on a recours à la méthode de perturbation par densité. Dans les deux cas évoqués, les tests numériques montrent que l'on approche la solution exacte du problème de Ransom avec une grande précision, et que l'on conserve la stabilité des calculs même avec un maillage de finesse relativement élevée.

Cette thèse étudie une classe des problèmes de Riemann-Hilbert Fuchsiens (à coefficients méromorphes dont tous les pôles sont d'ordre un). Les variétés de Frobenius apparaissent comme une formulation géométrique des structures d'équations de Witten-DijkgraafVerlande-Verlande (WDVV). Nous considérons ces variétés sur les espaces de Hurwitz vus, quant à eux, comme variétés réelles motivés par le fait qu'une variété de Frobenius semisimple peut être construite à partir d'une solution fondamentale du problème de Riemann-Hilbert associé. Une solution au problème Fuchsien de Riemann-Hilbert matriciel (problème de monodromie inverse) correspondant aux structures "réelles doubles" de Frobenius de Dubrovin sur les espaces de Hurwitz, a été construite. La solution est donnée en termes de certaines différentielles méromorphes integrées sur une base appropriée d'homologie relative de la surface de Riemann. La relation avec la solution du problème Fuchsien de Riemann-Hilbert pour les structures de Frobenius Hurwitz de Dubrovin est établie. Une solution du problème de Riemann-Hilbert correspondant aux déformations des "réelles doubles" est aussi donnée.

Dans cette thèse, nous développons des méthodes analytiques et numériques pour capturer les dynamiques asymptotiques de problèmes issus de la physique des plasmas et de la modélisation des mouvements collectifs dans les populations animales. Dans une première partie, nous présentons une méthode numérique Particle-In-Cell (PIC) pour le système Vlasov-Poisson préservant l'asymptotique quasi-neutre. Dans une seconde partie, nous étudions la limite macroscopique d'un modèle de Vicsek décrivant des interactions d'alignement entre deux populations, une population à l'arrêt et une population en mouvement. Nous sélectionnons ensuite un schéma numérique pour capturer les solutions du modèle macroscopique de Vicsek correspondant à la dynamique particulaire sous-jacente. La troisième partie est dédiée à l'étude des transitions compressible-incompressible apparaissant sous l'effet d'une contrainte de congestion dans un modèle macroscopique de déplacement collectif. Des schémas numériques préservant l'asymptotique de congestion sont ensuite mis au point pour le système d'Euler avec une contrainte de densité maximale.

Dans la première partie de cette thèse nous traitons l’écoulement diphasique compositionnel, partiellement miscible et compressible en milieux poreux. Déplacement d'une phase par un autre est analysé. Nous examinons les mélanges non idéals, la pression est variable, et les concentrations de phase, la densité et la viscosité sont les fonctions de la pression. Le processus est décrit par le problème de Riemann qui admet des solutions discontinues. Nous avons développé une méthode numérique-analytique de solution pour déterminer les paramètres à tous les chocs avant résoudre les équations de flux. Cette méthode est basée sur la séparation de thermodynamique et hydrodynamique, proposée dans [Oladyshkin, Panfilov 2006] et qui était inapplicable à problème de Riemann, en raison de manque des conditions d’Hugoniot. Dans cette thèse, nous avons construit les conditions supplémentaires d'Hugoniot. Dans la deuxième partie, nous examinons l'écoulement diphasique lors que les zones monophasique apparaissent, dans cette zone, le fluide est sur/sous-saturés et les équations diphasique dégénèrent.Nous avons proposé de décrire les zones diphasique et sur/sous-saturés avec un système uniforme des équations diphasique classique en étendant le concept de saturation d'être négatif et supérieur à un. Physiquement, cela signifie que les états monophasiques sont considérés comme des états diphasiques consistant une phase imaginaire avec la saturation négative. Une telle extension de la saturation exige développement des conditions de consistance qui sont fait dans cette thèse.La dernière partie est consacrée ensuite à étendre le modèle HT-split pour le cas d’écoulement triphasique compositionnel. Nous avons obtenu le modèle asymptotique, dans lequel la thermodynamique et l'hydrodynamique sont séparées / In the first part of thesis we deal with two-phase multicomponent, partially miscible, compressible flow in porous media. Displacement of one phase by another is analyzed. We examine non ideal solutions, pressure is variable, and phase compositions, densities and viscosities are variable functions of pressure.The process is described by Riemann problem which admits discontinuous solutions.We developed a numerical-analytical method of solution to explicitly determine all shock parameters before solving the flow equations. This method is based on splitting thermodynamics and hydrodynamics, suggested in [Oladyshkin, Panfilov 2006]. Earlier this method was inapplicable to Riemann problem, due to the lack of Hugoniot conditions. In this thesis we have constructed additional Hugoniot conditions.In the second part we examine two-phase flow when the single-phase zones appear, in this zone the fluid is over/under-saturated and two-phase flow equations degenerate and they cannot be used. We proposed to describe two-phase and over/under-saturated single-phase zones by uniform system of classic two-phase equations while extending the concept of phase saturation to be negative and higher than one. Physically it means that the oversaturated single-phase states are considered as pseudo two-phase states consisting an imaginary phase with negative saturation. Such an extension of saturation requires developing some consistence conditions which have developed in this thesis.The last part then is devoted to extend the HT-split model to the case of three-phase compositional flow. We have obtained the general asymptotic model, in which the thermodynamics and hydrodynamics are split

Ce travail de thèse est consacré à deux thèmes de recherche en Calcul Scientifique liés par l'approximation numérique de problèmes en mécanique des fluides. Le premier thème concerne l'approximation numérique des équations de Stokes, modélisant les écoulements de fluides incompressibles à vitesse faible. Ce thème est présent dans plusieurs travaux en Calcul Scientifique. La discrétisation en temps est réalisée à l'aide de la méthode de projection. La discrétisation en espace utilise la méthode des éléments finis mixtes hybrides qui permet d'imposer de façon exacte la contrainte d'incompressibilité. Cette approche est originale: la méthode des éléments mixtes hybrides est couplée avec une méthode d'éléments finis standards. L'ordre de convergence des deux méthodes est préservé. Le second thème concerne la mise au point de méthodes numériques de type volumes finis pour la résolution de l'équation Level Set. Ces équations interviennent de manière essentielle dans la résolution des problèmes de propagation d'interfaces. Dans cette partie, nous avons développé une nouvelle méthode d'ordre 2 de type MUSCL pour résoudre le système hyperbolique résultant de l'équation Level Set. Nous illustrons ces propriétés par des applications numériques. En particulier nous avons regardé le cas du problème des deux demi-plans pour lequel notre schéma donne une approximation pour le gradient de la fonction Level Set. Par ailleurs, l'ordre de précision attendu est obtenu avec les normes $L_1$ et $L_{\infty}$ pour des fonctions régulières. Pour finir, il est à noter que notre méthode peut être facilement étendue aux problèmes d'Hamilton-Jacobi du premier et du second ordre.

Cette thèse est consacrée à la simulation numérique de problèmes d'écoulements de fluides raides régis par des systèmes de lois de bilan non homogènes, dans des configurations monodimensionnelles et bidimensionnelles. La méthode numérique utilisée est une extension d'un schéma à deux pas (SRNH), comportant un paramètre \alpha^n_(j+\frac(1)(2)) ajustable, proposé par le professeur F.Benkhaldoun dans un cadre monodimensionnel. Ainsi, en un premier temps on a introduit une variante SRNHR, obtenue en remplaçant la vitesse numérique (\frac(\Delta x)(\Delta t)) par la vitesse de Rusanov locale, en vue de l'extension du schéma au cas bidimensionnel. Par la suite, une analyse de stabilité du schéma, révèle que celui-ci peut être d'ordre 1 ou 2 selon la valeur du paramètre \alpha^n_(j+\frac(1)(2)). Une stratégie de variation de ce paramètre, basée sur la théorie des limiteurs a alors été adoptée. Le schéma peut ainsi être rendu d'ordre 1 dans les zones à forte variation de l'écoulement, et d'ordre 2, là où l'écoulement est régulier. Ensuite on a établi les conditions pour que ce schéma respecte la C-propriété exacte introduite par Bermùdez et Vazquez. Une étude d'implémentation des conditions aux limites, adaptée à ce schéma, a également été menée en se basant sur les invariants de Riemann. Dans la deuxième partie de la thèse, on a appliqué ce schéma à des systèmes monophasiques homogènes et non homogènes. Par exemple on a réalisé la simulation du problème de rupture de barrage sur une marche, pour des configurations 1D et 2D, en menant en particulier une étude de convergence numérique via la détermination des courbes d'erreurs. Enfin, on a utilisé le schéma pour la simulation numérique de systèmes diphasiques (Ransom 1D et 2D).