▼ 

Nous présentons de nouveaux schémas numériques pour l'approximation de problèmes de diffusion hétérogène et anisotrope sur des maillages généraux. Sous des hypothèses correspondant aux cas industriels, nous montrons qu'un premier schéma, qui est centré sur les mailles, possède un petit stencil et converge dans le cas de tenseurs discontinus. La preuve de la convergence repose sur des propriétés de consistance des gradients discrets issus du schéma. Dans une seconde partie, nous proposons des méthodes de correction non linéaire du schéma initial pour obtenir le principe du maximum. L'efficacité de ces schémas est étudiée sur des tests numériques ayant fait l'objet de bancs d'essais d'une grande variété de schémas de volumes finis. Les comparaisons avec les schémas volumes finis classiques montrent l'apport de ces schémas en termes de précision. Nous montrons ainsi le bon comportement de ces schémas sur des maillages déformés, et le maintien de la précision des schémas non-linéaires, alors que les oscillations ont été supprimées

We present a new scheme for the discretization of heterogeneous anisotropic diffusion problems on general meshes. With light assumptions, we show that the algorithm can be written as a cell-centered scheme with a small stencil and that it is convergent for discontinuous tensors. The key point of the proof consists in showing both the strong and the weak consistency of the method. Besides, we study non-linear corrections to correct the FECC scheme, in order to satisfy the discrete maximum principle (DMP).The efficiency of the scheme is demonstrated through numerical tests of the 5th & 6th International Symposium on Finite Volumes for Complex Applications - FVCA 5 & 6. Moreover, the comparison with classical finite volume schemes emphasizes the precision of the method. We also show the good behaviour of the algorithm for nonconforming meshes. In addition, we give some numerical tests to check the existence for the non-linear FECC schemes

Advisors/Committee Members: Eymard, Robert (thesis director).

▼ 

Dans cette thèse, nous proposons un modèle pour le calcul des écoulements le long des failles. Ce modèle, baptisé modèle double interface permet de traiter deux difficultés majeures rencontrées lors de la modélisation des failles. Tout d'abord, l'utilisation d'un modèle interface, dans lequel les failles sont représentées par des éléments de dimension inférieure permet de s'affranchir du problème d'échelle spatiale. Ensuite, l'utilisation de deux interfaces pour représenter chaque faille permet de traiter naturellement les maillages non-conformes apparaissant dans ce type de problème. Les questions de failles non-planes et de réseaux de failles sont aussi abordées. Ce modèle est validé numériquement sur différents cas tests académiques et un cas synthétique inspiré du stockage du CO2 a aussi été réalisé. Finalement, une étude théorique a été menée afin de confirmer mathématiquement l'approche retenue.

In this thesis, we are interested in the modelisation of fluid flow along conductive faults. This model, so-called double interface model tackles two majors difficulties encountered when modelising faults. First of all, the use of an interface model, in which the faults are represented by lower dimension elements allows to treat the problem of space scale. Then, the use of two interfaces to modelise each fault allows to handle quite naturally the non-matching grid problem arising from this kind of problem. The question of non-planar fault and fault networks is also addressed. This model is then validated on several academic test cases and a synthetic case inspire by CO2 storage is also performed. Finally, a theoric study is also conducted in order to validate our approach.

Advisors/Committee Members: Gallouët, Thierry (thesis director).

▼ 

Dans cette thèse, nous nous intéressons à la résolution des équations SPN du transport de neutrons au sein des cœurs de réacteurs nucléaires à eau pressurisée. Ces équations forment un problème aux valeurs propres généralisé. Dans notre étude nous commençons par le problème source associé et ensuite nous étudions le problème aux valeurs propres. Un cœur de réacteur est composé de différents milieux: le combustible, le fluide caloporteur, le modérateur... à cause de ces hétérogénéités de la géométrie, le flux solution du problème source peut être peu régulier. Nous proposons l’analyse numérique de l’approximation de la solution par la méthode des éléments finis du problème source dans le cas où la solution est peu régulière. Pour le problème aux valeurs propres, dans le cas mixte, les théories déjà développées ne s’appliquent pas. Nous proposons ici une nouvelle méthode pour étudier la convergence de la méthode des éléments finis mixtes pour les problèmes aux valeurs propres. Pour les solutions peu régulières, la montée en ordre de la méthode des éléments finis n’améliore pas l’approximation du problème, il faut raffiner le maillage aux alentours des singularités de la solution. La géométrie des cœurs de réacteur se prête bien aux maillages cartésiens, mais leur raffinement augmente vite leur nombre de degrés de liberté. Pour palier à cette augmentation, nous proposons ici une méthode de décomposition de domaine qui permet d’utiliser des maillages globalement non-conformes.

In this thesis, we investigate the resolution of the SPN neutron transport equations in pressurized water nuclear reactor. These equations are a generalized eigenvalue problem. In our study, we first considerate the associated source problem and after we concentrate on the eigenvalue problem. A nuclear reactor core is composed of different media: the fuel, the coolant, the neutron moderator... Due to these heterogeneities of the geometry, the solution flux can have a low-regularity. We propose the numerical analysis of its approximation with finite element method for the low regular case. For the eigenvalue problem under its mixed form, we can not rely on the theories already developed. We propose here a new method for studying the convergence of the SPN neutron transport eigenvalue problem approximation with mixed finite element. When the solution has low-regularity, increasing the order of the method does not improve the approximation, the triangulation need to be refined near the singularities of the solution. Nuclear reactor cores are well-suited for Cartesian grids, but the refinement of these sort of triangulations increases rapidly their number of degrees of freedom. To avoid this drawback, we propose domain decomposition method which can handle globally non-conforming triangulations.

Advisors/Committee Members: Ciarlet, Patrick (thesis director), Jamelot, Erell (thesis director).

▼ 

Cette thèse propose l’étude d’une méthode numérique permettant de simuler un plasma. On considère un ensemble de particules, dont le mouvement est régi par l’équation de Vlasov, et qui est sensible aux forces électromagnétiques, qui proviennent des équations de Maxwell. La résolution numérique des équations de Vlasov-Maxwell est réalisée par une méthode Particle In Cell (PIC). La résolution des équations de Maxwell nécessite un maillage suffisamment fin afin de modéliser correctement les problémes multi-échelles que nous souhaitons traiter. Cependant, mailler finement tout le domaine de calcul a un coût. La nouveauté de cette thèse est de proposer un solveur PIC sur des maillages cartésiens localement raffinés, des maillages non conformes, afin de garantir la bonne modélisation du phénomène physique et d’éviter une trop forte pénalisation des temps de calcul.Nous utilisons une méthode Galerkin Discontinue en domaine temporelle (GDDT), qui offre l’avantage d’être d'une grande flexibilité dans le choix du maillage et qui est une méthode d’ordre élevé. Un point fondamental dans l’étude des solveurs PIC concerne le respect de la conservation de la charge. Nous proposons deux approches afin de traiter cet aspect. La première concerne les méthodes utilisant un système de Maxwell augmenté, dont la nouveauté a été de les étendre aux maillages non conformes. La seconde approche repose sur une méthode originale de pré-traitement du calcul du terme source de courant.

This thesis deals with the study of a numerical method to simulate a plasma. We consider a set of particles whose displacement is governed by the Vlasov equation and which creates an electromagnetic field thanks to Maxwell equations. The numerical resolution of the Vlasov-Maxwell system is performed by a Particle In Cell (PIC) method. The resolution of Maxwell equations needs a sufficiently fine mesh to correctly simulate the multi scaled problems that we have to face. Yet, a uniform fine mesh of the whole domain has a prohibitive cost. The novelty of this thesis is a PIC solver on locally refined Cartesian meshes : non conforming meshes, to guarantee the good modeling of the physical phenomena and to avoid too large CPU time. We use the Discontinuous Galerkin in Time Domain (DGTD) method which has the advantage of a great flexibility in the choice of the mesh and which is a high order method. A fundamental point in the study of PIC solvers is the respect of the charge conserving law. We propose two approaches to tackle this point. The first one deals with augmented Maxwell systems, that we have adapted to non conforming meshes. The second one deals with an original method of preprocessing of the calculation of the current source term.

Advisors/Committee Members: Sonnendrücker, Eric (thesis director).

▼ 

Cette thèse s'intéresse à la conception de méthodes de discrétisation non-conforme pour un modèle de poromécanique. Le but de ce travail est de simplifier les couplages liant la géomécanique d'un milieu poreux à l'écoulement polyphasique compositionnel ayant cours en son sein tels qu'ils sont réalisés actuellement dans l'industrie pétrolière, en discrétisant sur un même maillage, typiquement non-conforme car à l'image de la lithologie, la mécanique et l'écoulement. La nouveauté consiste donc à traiter la mécanique par une méthode d'approximation non-conforme sur maillages généraux. Dans cette thèse, nous nous concentrons sur un modèle d'élasticité linéaire. Les difficultés inhérentes à son approximation non-conforme sont son manque de coercivité (se traduisant par la nécessité de satisfaire une inégalité de Korn sur un espace discret discontinu), ainsi que le phénomène de verrouillage numérique lorsque le matériau tend à devenir incompressible. Dans une première partie, nous construisons un espace d'approximation sur maillages généraux, s'apparentant à une extension de l'espace de Crouzeix-Raviart. Nous explicitons ses propriétés d'approximation et de conformité, et montrons que ce dernier est adapté à une discrétisation primale coercive et robuste au locking du modèle d'élasticité sur maillages généraux. La méthode proposée est moins coûteuse que son équivalent éléments finis (en termes de propriétés) P2. Nous nous intéressons dans une deuxième partie à l'approximation non-conforme d'un modèle couplé de poroélasticité. Nous étudions la convergence d'une famille de schémas numériques dont la discrétisation en espace utilise le formalisme des schémas Gradient, auquel appartient la méthode développée pour la mécanique. Nous prouvons la convergence de telles approximations vers la solution de régularité minimale du problème continu, indépendamment des paramètres physiques du système

This manuscript focuses on the conception of nonconforming discretization methods for a poromechanical model. The aim of this work is to ease the coupling between the geomechanics and the multiphase compositional Darcy flow in porous media by discretizing mechanics and flow on the same mesh, typically nonconforming as it represents the lithology. Hence, the novelty hinges on a nonconforming treatment of mechanics on general meshes. In this work, we focus on a linear elasticity model. The nonconforming approximation of such a model is not straightforward owing to its lack of coercivity (meaning that a discrete Korn's inequality must hold on a discontinuous discrete space) and to the numerical locking phenomenon occurring as the material becomes incompressible. In a first part, we design an approximation space on general meshes, which can be viewed as an extension of the so-called Crouzeix-Raviart space. We study its approximation and conformity properties, and prove that this latter is well-adapted to the design of a primal, coercive, and locking-free discretization of the elasticity model on general meshes. The proposed method is less…

Advisors/Committee Members: Eymard, Robert (thesis director).

▼ Les travaux de cette thèse portent sur des méthodes de volumes finis sur maillages quelconque pour la discrétisation de problèmes d'évolution non linéaires modélisant le transport de contaminants en milieu poreux et les écoulements diphasiques.Au Chapitre 1, nous étudions une famille de schémas numériques pour la discrétisation d'une équation parabolique dégénérée de convection-reaction-diffusion modélisant le transport de contaminants dans un milieu poreux qui peut être hétérogène et anisotrope. La discrétisation du terme de diffusion est basée sur une famille de méthodes qui regroupe les schémas de volumes finis hybrides, de différences finies mimétiques et de volumes finis mixtes. Le terme de convection est traité à l'aide d'une famille de méthodes qui s'appuient sur les inconnues hybrides associées aux interfaces du maillage. Cette famille contient à la fois les schémas centré et amont. Les schémas que nous étudions permettent une discrétisation localement conservative des termes d'ordre un et d'ordre deux sur des maillages arbitraires en dimensions d'espace deux et trois. Nous démontrons qu'il existe une solution unique du problème discret qui converge vers la solution du problème continu et nous présentons des résultats numériques en dimensions d'espace deux et trois, en nous appuyant sur des maillages adaptatifs.Au Chapitre 2, nous proposons un schéma de volumes finis hybrides pour la discrétisation d'un problème d'écoulement diphasique incompressible et immiscible en milieu poreux. On suppose que ce problème a la forme d'une équation parabolique dégénérée de convection-diffusion en saturation couplée à une équation uniformément elliptique en pression. On considère un schéma implicite en temps, où les flux diffusifs sont discrétisés par la méthode des volumes finis hybride, ce qui permet de pouvoir traiter le cas d'un tenseur de perméabilité anisotrope et hétérogène sur un maillage très général, et l'on s'appuie sur un schéma de Godunov pour la discrétisation des flux convectifs, qui peuvent être non monotones et discontinus par rapport aux variables spatiales. On démontre l'existence d'une solution discrète, dont une sous-suite converge vers une solution faible du problème continu. On présente finalement des cas test bidimensionnels.Le Chapitre 3 porte sur un problème d'écoulement diphasique, dans lequel la courbe de pression capillaire admet des discontinuité spatiales. Plus précisément on suppose que l'écoulement prend place dans deux régions du sol aux propriétés très différentes, et l'on suppose que la loi de pression capillaire est discontinue en espace à la frontière entre les deux régions, si bien que la saturation de l'huile et la pression globale sont discontinues à travers cette frontière avec des conditions de raccord non linéaires à l'interface. On discrétise le problème à l'aide d'un schéma, qui coïncide avec un schéma de volumes finis standard dans chacune des deux régions, et on démontre la convergence d'une solution approchée vers une solution faible du problème continu. Les test numériques… Advisors/Committee Members: Hilhorst, Danielle (thesis director).