La capacité de simuler la dynamique de la couche limite turbulente représente aujourd'hui un enjeu important pour la prévision de l'aérodynamique instationnaire et de l'aéroacoustique des aéronefs et des véhicules terrestres. Aussi, les travaux présentés dans ce manuscrit proposent une méthode originale de simulation de la dynamique des écoulements turbulents pariétaux, à haut nombre de Reynolds, sur des géométries curvilignes. L'approche ZDES, dans son mode " Wall Modelled Large Eddy Simulation ", est ici retenue. Dans un premier temps, une méthode de réactivation turbulente, par combinaison de l'injection de turbulence synthétique et de l'application de termes de forçage, est développée afin de permettre une résolution WMLES locale des régions d'intérêt au sein d'une simulation majoritairement RANS. Puis l'étude des interactions entre la physique et la méthode de résolution numérique (ZDES), sur le cas d'une couche limite turbulente sans gradient de pression en développement spatial jusqu'à Re13 000, a conduit à une généralisation des conditions d'emploi du mode WMLES de la ZDES. Des post-traitements instationnaires avancés ont permis de démontrer la capacité de la méthode à simuler la dynamique particulière de la zone externe à ces grands nombres de Reynolds. Enfin, les différentes méthodes développées au cours de cette étude ont été appliquées à la simulation d'une manche à air coudée de section rectangulaire. Cette nouvelle méthode a permis de simuler avec succès la dynamique multi-échelles de cet écoulement et des distorsions dynamiques associées, pour un coût environ 50 fois inférieur à celui d'une approche LES classique

Cette thèse est une contribution à la résolution numérique d'une loi de conservation hyperbolique résultante d'un couplage entre les équations de Saint-Venant, associée à la modélisation des écoulements en eaux peu profondes, et l'équation de transport-diffusion d'un polluant non actif. Le modèle mathématique utilisé est bi-dimensionnel, intégrant des termes de friction, de diffusion, des tensions de surface et un terme tenant compte la variation de la bathymétrie. Nous présentons un modèle numérique basé sur un schéma volumes finis bidimensionnel d'ordres deux, conservatif et consistant, sur un maillage non structuré adaptatif. Ce modèle préserve la positivité de la hauteur d'eau et l'état stationnaire associé au lac au repos, il permet de capturer avec précision les ondes de chocs. Dans le temps une extension à l'ordre deux est garantie en utilisant un schéma de Runge-Kutta ce qui permettra de prendre en compte les différentes vitesses de propagation de l'information présentes dans les différents problèmes traités. Nous appliquons le modèle numérique développé sur plusieurs problèmes. Entre autre, la simulation d'une propagation d'une onde de crue, écoulement autour d'une singularité géométrique, écoulement sur des fonds variables et présentant des fronts raides. Et en fin, L'étude numérique s'achève par une application du modèle pour la simulation du transport de polluant dans une géométrie réelle avec une bathymétrie fortement variable telle que présente la baie de Tanger.

Nous etudions dans cette these differentes lois de conservation hyperboliques associees a la modelisation des ecoulements en eaux peu profondes.<br />Nous nous consacrons d'abord a l'analyse numerique du systeme de Saint-Venant avec termes sources. Nous presentons un schema volumes finis bidimensionnel d'ordre 2, conservatif et consistant, qui s'appuie sur une interpretation cinetique du systeme et une methode de reconstruction hydrostatique des variables aux interfaces. Ce schema preserve la positivite de la hauteur d'eau et l'etat stationnaire associe au lac au repos.<br />Nous etendons ensuite l'interpretation cinetique au couplage du systeme avec une equation de transport. Nous construisons un schema volumes finis a deux pas de temps, qui permet de prendre en compte les differentes vitesses de propagation de l'information presentes dans le probleme. Cette approche preserve les proprietes de stabilite du systeme et reduit sensiblement la diffusion numerique et les temps de calcul.<br />Nous proposons egalement un nouveau modele de Saint-Venant multicouche, qui permet de retrouver des profils de vitesse non constants, tout en preservant le caractere invariant et bidimensionnel du domaine de definition. Nous presentons sa derivation a partir des equations de Navier-Stokes et une etude de stabilite - energie, hyperbolicite. Nous etudions egalement ses relations avec d'autres modeles fluides et sa mise en oeuvre numerique, la encore basee sur l'utilisation des schemas cinetiques.<br />Enfin nous etablissons un theoreme d'unicite pour les lois de conservation scalaires avec flux discontinus. La preuve est basee sur l'utilisation d'une nouvelle famille d'entropies, qui constituent une adaptation naturelle des entropies de Kruzkov classiques au cas discontinu. Cette methode permet de lever certaines hypotheses classiques sur le flux - convexite, existence de bornes BV, nombre fini de discontinuites - et ne necessite pas l'introduction d'une condition d'interface.

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 est consacrée à l'étude des systèmes non linéaires hyperboliques de lois de conservation avec des termes sources autorisées à devenir raides. Les applications viennent pour la plupart de problèmes physiques, par exemple l'equation de Burgers ou de Buckley-Leverett et le système d'Euler muni de lois de pression pour un gaz parfait ou une version simplifiée de mélange à deux phases. La première partie traite des aspects théoriques de ce sujet. Nous nous concentrons sur le cas unidimensionnel et considérons l'équation scalaire et d'un cas particulier d'un système d'écoulement diphasique réactif. La question principale à gérer lorsqu'on établit la convergence forte d'une suite d'approximations au moyen de la théorie de compacité par compensation dans de tels cas est de trouver des régions invariantes dans l'espace des phases. Pour le problème scalaire, il suffit que le terme source conserve une borne uniforme sur la norme L ∞ de la séquence. Mais c'est moins évident pour un système car les zones invariantes sont plutôt compliquées et peuvent être perturbées en raison de la forme de certains termes de forçage. La deuxième partie s'intéresse aux critères de stabilité pour les schémas numériques implicites impliquant des quadratures en temps sophistiquées. Nous présentons d'abord l'analyse d'une semi-discrétisation (méthode des lignes) afin d'en tirer des conditions suffisantes pour le système différentiel d'être bien posé. Cela amène des restrictions assez naturelles pour assurer la stabilité d'une méthode à une seule étape implicite: le pas de temps doit être petit près des points de compression et de répulsion du terme source. Une analyse de variation totale est ensuite réalisée pour dériver les régimes du second ordre sans oscillations autour de chocs pour le plus grand pas de temps possible. Des expériences numériques sont présentés sur des problèmes scalaires non linéaires. La troisième partie est une suite directe afin de corriger certaines des difficultés rencontrées auparavant. L'objectif principal est de construire une nouvelle discrétisation qui peut gérer un terme source arbitrairement raide sans aucune influence (sauf la condition habituelle CFL) sur le pas de temps. Suivant les idées de J. Greenberg et A.-Y. LeRoux, nous calculons d'abord un schéma de Godunov scalaire afin d'obtenir un résultat de convergence vers la solution entropique à l'aide des estimations BV. Comme la caractéristique principale est de gérer les sources au moyen des produits non-conservatifs, nous décidons d'utiliser le formalisme proposé par G. Dalmaso, PG LeFloch et F. Murat pour étendre ces idées à des systèmes non-homogènes. Suivant I. Toumi, nous introduisons des matrices de type Roe le long des chemins régularisants afin d'obtenir des approximations qui montrent des propriétés intéressantes. En fait, c'est ce que l'on appelle les régimes bien équilibrés, à savoir les régimes qui préservent les intégrales premières théoriques du mouvement. Des tests numériques sont basés sur le système d'Euler avec termes sources géométriques. Le quatrième (et dernière) partie concerne la mise en place de cette formulation non-conservative dans un schéma de type flux-splitting. L'objectif est de construire une méthode très robuste, bien équilibrée et facile à mettre en oeuvre. Les relations de saut généralisées provenant des sources sont détaillées et le système est testé dans un très large éventail de problèmes: les flux de tuyères, diphasique avec chimie et amortissement, les systèmes de relaxation ... Un chapitre est également consacré à des cas de résonance: en suivant Majda, nous étudions le problème stationnaire pour valider nos résultats. Ce travail conclut avec un écoulement diphasique en deux dimensions qui évolue de manière assez compliquée, les relations non conservatives étant traitées par un processus itératif convergent.

Les machines électriques sont modélisées pour prédire leurs performances et optimiser leur rendement. Cette modélisation peut être faite par des simulations avec la méthode des éléments finis. En particulier, les machines asynchrones nécessitent des simulations 3D pour prendre en compte les courants de Foucault et les têtes de bobines. Dans le logiciel Flux®, des formulations 3D basées sur le potentiel scalaire magnétique sont utilisées avec succès depuis de nombreuses années. Néanmoins, des coupures mathématiques artificielles sont nécessaires, lorsque le domaine n'est pas simplement connexe.Afin de se libérer de ces contraintes de connexité, des formulations en potentiel vecteur magnétique ont été étudiées et développées. En 3D, leur mise en œuvre nécessite l'utilisation d'éléments finis d'arêtes afin de respecter la nature des champs. Avec les éléments d'arêtes, les formulations sont généralement résolues avec une condition de jauge pour les solveurs directs comme pour les solveurs itératifs. De nouvelles formulations en potentiel vecteur magnétique auto-jaugées ont été développées permettant la prise en compte des bobines maillées et des bobines non maillées. La prise en compte du mouvement est relativement simple à mettre en œuvre pour les formulations en potentiel scalaire magnétique avec l'interpolation nodale.Avec les éléments d'arête, l'interpolation est plus délicate. C'est pourquoi la méthode des éléments avec joints a été développée pour prendre en compte le mouvement dans un cas général. / Electric machines are modeled in order to predict their performance and to optimize their output. This modeling can be done by simulation with the finite element method. In particular, induction machines require 3D simulation to take into account eddy currents and coils overhangs. In the Flux® software, 3D formulations based on magnetic scalar potential has been used with success for many years. Nevertheless, artificial mathematical cuts are necessary, when the domain is not simply connected.In order to avoid connection constraints, magnetic vector potential formulations have been studied and developed. In 3D, their implementation requires the use of edge elements to respect the nature of fields. With edge elements, formulations are generally solved with a gauge condition for direct solvers as well as for iterative solvers. New auto-gauged magnetic vector potential formulations have been developed to take into account meshed coils and non-meshed coils. Consideration of movement is relatively simple to implement for magnetic scalar potential formulations with nodal interpolation. With edge elements, the interpolation is more delicate. For this reason, the mortar method has been developed to take into account movement in a general case.

Cette étude porte sur le calcul numérique du champ acoustique rayonné par des écoulements subsoniques turbulents présentant des inhomogénéités de température. Des méthodes hybrides sont développées grâce à un développement de Janzen-Rayleigh des équations de Navier-Stokes. L'écoulement est résolu par un calcul quasi incompressible puis les perturbations acoustiques sont propagées selon deux méthodes : les équations d'Euler linéarisées (EEL) et l'approximation à faible nombre de Mach perturbée(PLMNA). Les méthodes sont validées sur des cas simples puis appliquées à une couche de mélange isotherme et anisotherme en développement spatial. / This study focuses on the numerical calculation of the acoustic field radiated by subsonic turbulent flows with temperature inhomogeneities. Hybrid methods are developed through a Rayleigh-Janzen expansion of the Navier-Stokes equations. The flow is solved in a quasi-incompressible way then the acoustic disturbances are propagated by two methods : the linearized Euler's equations (EEL) and the perturbed low Mach number approximation (PLMNA). The methods are validated on simple cases and then applied to an isothermal and non isothermal spatially evolving mixing layer.

Cette étude porte sur le calcul numérique du champ acoustique rayonné par des écoulements subsoniques turbulents présentant des inhomogénéités de température. Des méthodes hybrides sont développées grâce à un développement de Janzen-Rayleigh des équations de Navier-Stokes. L'écoulement est résolu par un calcul quasi incompressible puis les perturbations acoustiques sont propagées selon deux méthodes : les équations d'Euler linéarisées (EEL) et l'approximation à faible nombre de Mach perturbée (PLMNA). Les méthodes sont validées sur des cas simples puis appliquées à une couche de mélange isotherme et anisotherme en développement spatial.