Le cadre général de ce travail est l'étude des problématiques de couplage de modèles dans le contexte applicatif de l'hydraulique et de l'océanographie. Ce sont des domaines dans lesquels les besoins sont nombreux, utilisant une vaste gamme de modèles tels que les équations de Navier-Stokes, de Saint-Venant ou des modèles intermédiaires tels que les équations primitives ou de Navier-Stokes hydrostatiques. Pour aborder ces questions de couplage, on a choisi d'utiliser les méthodes de Schwarz, pour leur caractère non-intrusif vis-à-vis des codes de calcul. On a travaillé tout d'abord sur les équations de Navier-Stokes, en envisageant deux discrétisations : une méthode de projection où la vitesse et la pression sont calculées dans deux étapes différentes, et un schéma couplé où toutes les inconnues sont calculées lors de la même étape. À chaque pas de temps, on est ramené à une équation de Stokes, pour laquelle différents opérateurs d'interface ont été considérés (Dirichlet, absorbants exacts et approchés...). La convergence de la méthode de Schwarz a été étudiée et testée numériquement dans la configuration classique de la cavité entraînée. Ces conditions ont ensuite été étendues, sans être testées numériquement, à des équations proches des équations de Navier-Stokes, à savoir dans un premier temps aux équations d'Oseen (qui sont une version linéarisée des équations de Navier-Stokes), et dans un deuxième temps en vue d'un couplage entre les équations de Navier-Stokes et les équations de Navier-Stokes hydrostatiques.

Les schémas de calcul déterministes permettent une modélisation à moindre coût du comportement de la population de neutrons en réacteur, mais sont traditionnellement construits sur des approximations (décomposition réseau/cœur, homogénéisation spatiale et énergétique…). La thèse revient sur une partie de ces sources d’erreur, de façon à rapprocher la méthode déterministe d’un schéma de référence. L’objectif est de profiter des architectures informatiques modernes (HPC) pour résoudre le problème neutronique à l’échelle du cœur 3D, tout en préservant l’opérateur de transport et une partie des hétérogénéités de la géométrie. Ce travail est réalisé au sein du solveur cœur Sn Minaret de la plateforme de calcul Apollo3® pour des réacteurs à neutrons rapides.Une méthode de décomposition de domaine en espace, est retenue. L'idée consiste à décomposer un problème de grande dimension en sous-problèmes "indépendants" de taille réduite. La convergence vers la solution globale est assurée par échange de flux angulaires entre sous-domaines au cours d'un processus itératif. En favorisant un recours massif au parallélisme, les méthodes de décomposition de domaine contribuent à lever les contraintes en mémoire et temps de calcul. La mise en place d'un parallélisme hybride, couplant les technologies MPI et OpenMP, est en particulier propice au passage sur supercalculateur. Une méthode d'accélération de type Coarse Mesh Rebalance  est ajoutée pour pallier à la pénalité de convergence constatée sur la méthode de décomposition de domaine. Le potentiel du nouveau schéma est finalement mis en évidence sur un coeur CFV 3D, construit en préservant l'hétérogénéité des assemblages absorbants. / Deterministic calculation schemes are devised to numerically solve the neutron transport equation in nuclear reactors. Dealing with core-sized problems is very challenging for computers, so much that the dedicated core codes have no choice but to allow simplifying assumptions (assembly- then core-scale steps…). The PhD work aims to correct some of these ‘standard’ approximations, in order to get closer of reference calculations: thanks to important increases in calculation capacities (HPC), nowadays one can solve 3D core-sized problems, using both high mesh refinement and the transport operator. Developments were performed inside the Sn core solver Minaret, from the new CEA neutronics platform Apollo3® for fast neutrons reactors of the CFV-kind.This work focuses on a Domain Decomposition Method in space. The fundamental idea involves splitting a core-sized problem into smaller and 'independent' subproblems. Angular flux is exchanged between adjacent subdomains. In doing so, all combined subproblems converge to the global solution at the outcome of an iterative process. Domain decomposition is well-suited to massive parallelism, allowing much more ambitious computations in terms of both memory requirements and calculation time. An hybrid MPI/OpenMP parallelism is chosen to match the supercomputers architecture. A Coarse Mesh Rebalance accelration technique is added to balance the convergence penalty observed using Domain Decomposition. The potential of the new calculation scheme is demonstrated on a 3D core of the CFV-kind, using an heterogeneous description of the absorbent rods.

Les calculs de réactivité constituent une brique fondamentale dans la simulation des coeurs des réacteurs nucléaires. Ceux-ci conduisent à la résolution de problèmes aux valeurs propres généralisées résolus par l'algorithme de la puissance inverse. A chaque itération, on est amené à résoudre un système linéaire de manière approchée via un algorithme d'itérations imbriquées. Il est difficile de traiter les modélisations très fines avec le solveur développé à EDF, au sein de la plate-forme Cocagne, en raison de la consommation mémoire et du temps de calcul. Au cours de cette thèse, on étudie une méthode de décomposition de domaine de type Schur dual. Plusieurs placements de l'algorithme de décomposition de domaine au sein du système d'itérations imbriquées sont envisageables. Deux d'entre eux ont été implémentés et les résultats analysés. Le deuxième placement, utilisant les spécificités des éléments finis de Raviart-Thomas et de l'algorithme des directions alternées, conduit à des résultats très encourageants. Ces résultats permettent d'envisager l'industrialisation de la méthodologie associée. / The reactivity computations are an essential component for the simulation of the core of a nuclear plant. These computations lead to generalized eigenvalue problems solved by the inverse power iteration algorithm. At each iteration, an algebraic linear system is solved through an inner/outer process. With the solver Cocagne developed at EDF, it is difficult to take into account very fine discretisation, due to the memory requirement and the computation time. In this thesis, a domain decomposition method based on the Schur dual technique is studied. Several placement in the inner/outer process are possible. Two of them are implemented and the results analyzed.The second one, which uses the specificities of the Raviart Thomas finite element and of the alternating directions algorithm, leads to very promising results. From these results the industrialization of the method can be considered.

La première partie de cette thèse (articles I et II) est consacrée à l'étude du comportement asymptotique des solutions d'un système dynamique du second ordre avec dissipation évanescente. Le système dynamique est étudié dans sa version continue et dans sa version discrète via un algorithme.La deuxième partie de cette thèse (articles III à VI) est consacrée à l'étude de plusieurs algorithmes de type proximal. Nous montrons que ces algorithmes convergent vers des solutions de certains problèmes de minimisation. Dans chaque cas, une application est donnée dans le cadre de la décomposition de domaine pour les EDP. / The first part of this thesis is devoted to the study of the asymptotic behavior of solutions of a second order dynamic system with vanishing dissipation. The dynamic system is studied in its continuous version and in its discrete version via an algorithm.The second part is about the study of several proximal-type algorithms. We show that these algorithms converge to solutions of some minimization problems. In each case, an application is given in the area of domain decomposition for PDE's.

Le contrôle de la qualité des calculs de structure suscite un intérêt croissant dans les processus de conception et de certification. Il repose sur l'utilisation d'estimateurs d'erreur, dont la mise en pratique entraîne un sur-coût numérique souvent prohibitif sur des calculs de grande taille. Le présent travail propose une nouvelle procédure permettant l'obtention d'une estimation garantie de l'erreur de discrétisation dans le cadre de problèmes linéaires élastiques résolus au moyen d'approches par décomposition de domaine. La méthode repose sur l'extension du concept d'erreur en relation de comportement au cadre des décompositions de domaine sans recouvrement, en s'appuyant sur la construction de champs admissibles aux interfaces. Son développement dans le cadre des approches FETI et BDD permet d'accéder à une mesure pertinente de l'erreur de discrétisation bien avant convergence du solveur lié à la décomposition de domaine. Une extension de la procédure d'estimation aux problèmes hétérogènes est également proposée. Le comportement de la méthode est illustré et discuté sur plusieurs exemples numériques en dimension 2. / The control of the quality of mechanical computations arouses a growing interest in both design and certification processes. It relies on error estimators the use of which leads to often prohibitive additional numerical costs on large computations. The present work puts forward a new procedure enabling to obtain a guaranteed estimation of discretization error in the setting of linear elastic problems solved by domain decomposition approaches. The method relies on the extension of the constitutive relation error concept to the framework of non-overlapping domain decomposition through the recovery of admissible interface fields. Its development within the framework of the FETI and BDD approaches allows to obtain a relevant estimation of discretization error well before the convergence of the solver linked to the domain decomposition. An extension of the estimation procedure to heterogeneous problems is also proposed. The behaviour of the method is illustrated and assessed on several numerical examples in 2 dimension.

Au cours des dernières décennies, les plasmas créés par une décharge micro-ondes ont de plusen plus attiré l’attention de la communauté scientifique aérospatiale sur le sujet du contrôled’écoulements. En effet, il a été démontré expérimentalement que le dépôt d’énergie obtenu parle plasma peut modifier les propriétés aérodynamiques de l’écoulement autour d’un objet telleque la trainée de frottement. Or, la conception et l’optimisation de ces actionneurs plasma entant que technique de contrôle d’écoulements nécessitent une compréhension approfondie de laphysique sous-jacente que les seules expériences sont incapables de fournir.Dans ce contexte, nous nous intéressons à la modélisation numérique de l’interaction desondes électromagnétiques avec un plasma et le gaz afin de mieux comprendre la nature desdécharges micro-ondes et leur applicabilité. La modélisation de ces phénomènes présente desdifficultés importantes en raison du couplage multi-physique et donc de la multitude des échellesspatiales et temporelles qui apparaissent. Ce travail de thèse traite des questions de physiqueet de mathématiques appliquées soulevées par la modélisation numérique de ces plasmas.La première partie du travail se focalise sur les questions de validité du modèle physique duclaquage micro-onde fondé sur l’approximation de champ effectif local. En raison des gradientsde densité du plasma très élevés, la validité du concept de champ effectif local peut être misen doute. Pour cela, un modèle fluide du second ordre est développé en incluant une equationd’énergie électronique non-locale. Cette modélisation permet de décrire de façon plus précisele dépôt d’énergie par plasma induisant la formation d’ondes de choc dans le gaz. Une analysedimensionnelle du système d’équations fluide permet de caractériser la non-localité en espace dubilan d’énergie électronique en fonction du champ électrique réduit et de la fréquence de l’onderéduite. Une discussion est également menée sur d’autres approximations des coefficients detransport électronique. Dans une deuxième partie, la construction et l’analyse d’une méthode multi-échelles derésolution numérique du problème de propagation des ondes électromagnétiques dans le plasmasont réalisées. Il s’agit du couplage entre les équations de Maxwell dans le domaine temporel avecune équation de quantité de mouvement pour les électrons. L’approche s’appuie sur la méthodede décomposition de domaine de type Schwartz, basée sur une formulation variationnelle duschéma de Yee et utilisant deux niveaux de grilles Cartésiennes emboitées. Une grille locale,appelée patch, est utilisée pour calculer de manière itérative la solution dans la région du plasmaoù une meilleure précision est requise. La méthode proposée permet le raffinement local etdynamique du maillage spatial tout en conservant l’énergie du système. Une analyse théorique dela convergence de l’algorithme pour les résolutions temporelles explicite et implicite est égalementréalisée. Dans la dernière partie, des simulations numériques sur le claquage micro-ondes et la formation de structures filamentaires de plasma sont conduites. Les effets de différents types d’approximations sur le modèle physique du plasma sont analysés. Puis, ces expériences numériques démontre la précision et l’efficacité, en terme de temps de calcul, de la méthode multi-échelleproposée. Enfin, on étudie les effets de chauffage du gaz sur la formation et l’entretien de structures filamentaires dans l’air à pression atmosphérique. Pour cela, le modèle micro-onde-plasma développé est couplé avec les équations de Navier-Stokes instationnaires pour les écoulements compressibles. Les simulations montrent des caractéristiques intéressantes de la dynamique deces structures plasma pendant le processus de chauffage du gaz, qui sont en accord étroit avec les données expérimentales. / In recent decades, microwave discharge plasmas have attracted increasing attention of aerospace scientific community to the subject of aerodynamic flow control because of their capability of sub- stantially modifying the properties of the flow around bodies by effective energy deposition. The design and optimization of these plasma actuators as flow control technique require a compre- hensive understanding of the complex physics involved that the sole experiments are incapable to provide.In this context, we have interest in the numerical modeling of the mutual interaction of elec- tromagnetic waves with plasma and gas in order to better understand the nature of microwave discharges and their applicability. A challenging problem arises when modeling such phenomena because of the coupling of different physics and therefore the multiplicity of spatial and tempo- ral scales involved. A solution is provided by this thesis work which addresses both physics and applied mathematics questions related to microwave plasma modeling.The first part of this doctorate deals with validity matters of the physical model of microwave breakdown based on the local effective field concept. Because of large plasma density gradients, the local effective field approximation is questionable and thus a second-order plasma fluid model is developed, where the latter approximation is replaced by the local mean energy approximation. This modeling approach enables to take into account the non-locality in space of the electron energy balance that provides a more accurate description of the energy deposition by microwave plasma leading to the shock waves formation into the gas. A dimensionless analysis of the plasma fluid system is performed in order to theoretically characterize the non-locality of the introduced electron energy equation as function of the reduced electric field and wave frequency. It also discusses other approximations related to the choice and method of calculation of electron transport coefficients.Concerning the mathematical aspects, the thesis work focuses on the design and the analysis of a multiscale method for numerically solving the problem of electromagnetic wave propagation in microwave plasma. The system of interest consists of time-dependent Maxwell’s equations coupled with a momentum transfer equation for electrons. The developed approach consists of a Schwartz type domain decomposition method based on a variational formulation of the standard Yee’s scheme and using two levels of nested Cartesian grids. A local patch of finite elements is used to calculate in an iterative manner the solution in the plasma region where a better precision is required. The proposed technique enables a conservative local and dynamic refinement of the spatial mesh. The convergence behavior of the iterative resolution algorithm both in an explicit and implicit time-stepping formulation is then analyzed.In the last part of the doctorate, a series of numerical simulations of microwave breakdown and the filamentary plasma array formation in air are performed. They allow to study in detail the consequences of the different types of physical approximations adopted in the plasma fluid model. Then, these numerical experiments demonstrate the accuracy and the computational efficiency of the proposed patch correction method for the problem of interest. Lastly, a numerically investigation of the effects of gas heating on the formation and sustaining of the filamentary plasma array in atmospheric-pressure air is carried out. For doing this, the developed microwave-plasma model is coupled with unsteady Navier-Stokes equations for compressible flows. The simulations provide interesting features of the plasma array dynamics during the process of gas heating, in close agreement with experimental data.

Lors d'un hypothétique accident grave de réacteur à eau sous pression, un mélange de matériaux fondus, appelé corium, issu de la fusion du cœur peut se relocaliser dans le puits de cuve constitué par un radier en béton. Les codes d'évaluation réacteur pour simuler la phénoménologie de l'interaction corium-béton sont basés sur une description à grande échelle des échanges qui soulève de nombreuses questions, tant sur la prise en compte des phénomènes multi-échelles mis en jeu que sur la structure adoptée de la couche limite au voisinage du front d'ablation. Dans ce contexte, l'objectif principal de ce travail consiste à aborder le problème de la structure de la couche limite par simulation numérique directe. Ce travail s'inscrit dans le cadre plus général d'une description et d'une modélisation multi-échelle des échanges, c'est-à-dire de l'échelle locale associée au voisinage du front d'ablation jusqu'à l'échelle du code d'évaluation réacteur. Une telle description multi-échelle des échanges soulève le problème de la description locale de l'écoulement multiphasique multiconstituant mais aussi le problème du changement d'échelle et en particulier le passage de l'échelle locale à l'échelle de description supérieure dite macroscopique associée aux mouvements convectifs dans le bain de corium. Parmi les difficultés associées au changement d'échelle, nous nous intéressons à la problématique de la construction de conditions aux limites effectives ou lois de parois pour les modèles macroscopiques. Devant la complexité du problème multiphasique multiconstituant posé au voisinage du front, cette contribution a été abordée sur un problème modèle. Des conditions aux limites dites effectives ont été construites dans le cadre d'une méthode de décomposition de domaine puis testées pour un problème d'écoulement laminaire de convection naturelle sur parois rugueuses. Mˆeme si le problème traité reste encore éloigné des applications visées, cette contribution offre de nombreuses perspectives et constitue une première étape d'une modélisation multiéchelle des échanges pour la problématique de l'interaction corium-béton. Dans le cas plus complexe des écoulements multiphasiques multiconstituants et devant les difficultés expérimentales associées, le développement de lois de parois pour les outils existants aux échelles de description supérieures nécessite, au préalable, de disposer d'un outil de simulation numérique directe de l'écoulement au voisinage du front d'ablation. L'outil développé dans ce travail correspond à un modèle de Cahn-Hilliard/Navier-Stokes pour un mélange diphasique (liquide-gaz) compositionnel (corium-béton fondu) s'appuyant sur une description du système selon trois paramètres d'ordre associés respectivement aux fractions volumiques du gaz et aux deux espèces miscibles de la phase liquide ainsi que sur une décomposition de l'énergie libre selon une contribution diphasique et compositionnelle. Les équations de transport sont dérivées dans le cadre de la thermodynamique des processus irréversibles et résolues sur la base d'une application éléments finis de la plate-forme PELICANS. Plusieurs expériences numériques illustrent la validité et les potentialités d'application de cet outil sur des problèmes diphasiques et/ou compositionnels. Enfin, à partir de l'outil développé, nous abordons par simulation numérique directe une étude de la structure de la couche limite au voisinage du front d'ablation pour des bétons siliceux et silico-calcaire. / In the late phases of some scenario of hypothetical severe accident in Pressurized Water Reactors, a molten mixture of core and vessel structures, called corium, comes to interact with the concrete basemat. The safety numerical tools are lumped parameter codes. They are based on a large averaged description of heat and mass transfers which raises some uncertainties about the multi-scale description of the exchanges but also about the adopted boundary layer structure in the vicinity of the ablation front. In this context, the aim of this work is to tackle the problem of the boundary layer structure by means of direct numerical simulation. This work joins within the more general framework of a multi-scale description and a multi-scale modeling, namely from the local scale associated with the vicinity of the ablation front to the scale associated with the lumped parameter codes. Such a multi-scale description raises not only the problem of the local description of the multiphase multicomponent flow but also the problem of the upscaling between the local- and the macro-scale which is associated with the convective structures within the pool of corium. Here, we are particularly interested in the building of effective boundary conditions or wall laws for macro-scale models. The difficulty of the multiphase multicomponent problem at the local scale leads us to consider a relatively simplified problem. Effective boundary conditions are built in the frame of a domain decomposition method and numerical experiments are performed for a natural convection problem in a stamp shaped cavity to assess the validity of the proposed wall laws. Even if the treated problem is still far from the target applications, this contribution can be viewed as a first step of a multi-scale modeling of the exchanges for the molten core concrete issue. In the more complicated case of multiphase multicomponent flows, it is necessary to have a direct numerical simulation tool of the flow at the local scale to build wall laws for macro-scale models. Here, the developed tool corresponds to a Cahn-Hilliard/Navier-Stokes model for a two-phase compositional system. It relies on a description of the system by three volume fractions and on a free energy composed by a two-phase part and a compositional part. The governing equations are derived in the frame of the thermodynamic of irreversible processes. They are solved on the basis of a finite element application of the object-oriented software component library PELICANS. Several numerical experiments illustrate the validity and the potentialities of application of this tool on two-phase compositional problems. Finally, using the developed tool, we tackle by means of direct numerical simulation the problem boundary layer structure in the vicinity of the ablation front for limestone-sand and siliceous concretes.

Cette thèse présente une méthode de résolution parallèle de systèmes linéaires creux qui combine efficacement les techniques de résolutions directes et itératives en utilisant une approche de type complément de Schur. Nous construisons une décomposition de domaine. L'intérieur des sous-domaines est éliminé de manière directe pour se ramener à un problème sur l'interface. Ce problème est résolu grâce à une méthode itérative préconditionnée par une factorisation incomplète. Un réordonnancement de l'interface permet la construction d'un préconditionneur global du complément de Schur. Des algorithmes minimisant le pic mémoire de la construction du préconditionneur sont proposés. Nous exploitons un schéma d'équilibrage de charge utilisant une répartition de multiples sous-domaines sur les processeurs. Les méthodes sont implémentées dans le solveur HIPS et des résultats expérimentaux parallèles sont présentés sur de grands cas tests industriels. / This thesis presents a parallel resolution method for sparse linear systems which combines effectively techniques of direct and iterative solvers using a Schur complement approach. A domain decomposition is built ; the interiors of the subdomains are eliminated by a direct method in order to use an iterative method only on the interface unknowns. The system on the interface (Schur complement) is solved thanks to an iterative method preconditioned by a global incomplete factorization. A special ordering on the Schur complement allows to build a scalable preconditioner. Algorithms minimizing the memory peak that appears during the construction of the preconditioner are presented. The memory is balanced thanks to a multiple domains per processors parallelization scheme. The methods are implemented in the HIPS solver and parallel experimental results are presented on large industrial test cases.

Ce travail de thèse est consacré à la mise en œuvre d’une méthode de décomposition de domaine appliquée à l’équation du transport. L’objectif de ce travail est l’accès à des solutions déterministes haute-fidélité permettant de correctement traiter les hétérogénéités des réacteurs nucléaires, pour des problèmes dont la taille varie d’un motif d’assemblage en 3 dimensions jusqu’à celle d’un grand cœur complet en 3D. L’algorithme novateur développé au cours de la thèse vise à optimiser l’utilisation du parallélisme et celle de la mémoire. La démarche adoptée a aussi pour but la diminution de l’influence de l’implémentation parallèle sur les performances. Ces objectifs répondent aux besoins du projet APOLLO3, développé au CEA et soutenu par EDF et AREVA, qui se doit d’être un code portable (pas d’optimisation sur une architecture particulière) permettant de réaliser des modélisations haute-fidélité (best estimate) avec des ressources allant des machines de bureau aux calculateurs disponibles dans les laboratoires d’études. L’algorithme que nous proposons est un algorithme de Jacobi Parallèle par Bloc Multigroupe. Chaque sous domaine est un problème multigroupe à sources fixes ayant des sources volumiques (fission) et surfaciques (données par les flux d’interface entre les sous domaines). Le problème multigroupe est résolu dans chaque sous domaine et une seule communication des flux d’interface est requise par itération de puissance. Le rayon spectral de l’algorithme de résolution est rendu comparable à celui de l’algorithme de résolution classique grâce à une méthode d’accélération non linéaire par la diffusion bien connue nommée Coarse Mesh Finite Difference. De cette manière une scalabilité idéale est atteignable lors de la parallélisation. L’organisation de la mémoire, tirant parti du parallélisme à mémoire partagée, permet d’optimiser les ressources en évitant les copies de données redondantes entre les sous domaines. Les architectures de calcul à mémoire distribuée sont rendues accessibles par un parallélisme hybride qui combine le parallélisme à mémoire partagée et à mémoire distribuée. Pour des problèmes de grande taille, ces architectures permettent d’accéder à un plus grand nombre de processeurs et à la quantité de mémoire nécessaire aux modélisations haute-fidélité. Ainsi, nous avons réalisé plusieurs exercices de modélisation afin de démontrer le potentiel de la réalisation : calcul de cœur et de motifs d’assemblages en 2D et 3D prenant en compte les contraintes de discrétisation spatiales et énergétiques attendues. / This thesis is devoted to the implementation of a domain decomposition method applied to the neutron transport equation. The objective of this work is to access high-fidelity deterministic solutions to properly handle heterogeneities located in nuclear reactor cores, for problems’ size ranging from colorsets of assemblies to large reactor cores configurations in 2D and 3D. The innovative algorithm developed during the thesis intends to optimize the use of parallelism and memory. The approach also aims to minimize the influence of the parallel implementation on the performances. These goals match the needs of APOLLO3 project, developed at CEA and supported by EDF and AREVA, which must be a portable code (no optimization on a specific architecture) in order to achieve best estimate modeling with resources ranging from personal computer to compute cluster available for engineers analyses. The proposed algorithm is a Parallel Multigroup-Block Jacobi one. Each subdomain is considered as a multi-group fixed-source problem with volume-sources (fission) and surface-sources (interface flux between the subdomains). The multi-group problem is solved in each subdomain and a single communication of the interface flux is required at each power iteration. The spectral radius of the resolution algorithm is made similar to the one of a classical resolution algorithm with a nonlinear diffusion acceleration method: the well-known Coarse Mesh Finite Difference. In this way an ideal scalability is achievable when the calculation is parallelized. The memory organization, taking advantage of shared memory parallelism, optimizes the resources by avoiding redundant copies of the data shared between the subdomains. Distributed memory architectures are made available by a hybrid parallel method that combines both paradigms of shared memory parallelism and distributed memory parallelism. For large problems, these architectures provide a greater number of processors and the amount of memory required for high-fidelity modeling. Thus, we have completed several modeling exercises to demonstrate the potential of the method: 2D full core calculation of a large pressurized water reactor and 3D colorsets of assemblies taking into account the constraints of space and energy discretization expected for high-fidelity modeling.

La simulation numérique des systèmes multicorps en présence d'interactions complexes, dont le contact frottant, pose de nombreux défis, tant en terme de modélisation que de temps de calcul. Dans ce manuscrit de thèse, nous étudions deux familles de décomposition de domaine adaptées au formalisme de la dynamique non régulière des contacts (NSCD). Cette méthode d'intégration implicite en temps de l'évolution d'une collection de corps en interaction a pour caractéristique de prendre en compte le caractère discret et non régulier d'un tel milieu. Les techniques de décomposition de domaine classiques ne peuvent de ce fait être directement transposées. Deux méthodes de décomposition de domaine, proches des formalismes des méthodes de Schwarz et de complément de Schur sont présentées. Ces méthodes se révèlent être de puissants outils pour la parallélisation en mémoire distribuée des simulations granulaires 2D et 3D sur un centre de calcul haute performance. Le comportement de structure des milieux granulaires denses est de plus exploité afin de propager rapidement l'information sur l'ensemble des sous domaines via un schéma semi-implicite d'intégration en temps.