La complexification de la modélisation multi-physique conduit d’une part à devoir simuler des systèmes d’équations différentielles ordinaires et d’équations différentielles algébriques de plus en plus grands en nombre d’inconnues et sur des temps de simulation longs. D’autre part l’évolution des architectures de calcul parallèle nécessite d’autres voies de parallélisation que la décomposition de système en sous-systèmes. Dans ce travail, nous proposons de concevoir des méthodes de décomposition de domaine pour la résolution d’EDO en temps. Nous reformulons le problème à valeur initiale en un problème aux valeurs frontières sur l’intervalle de temps symétrisé, sous l’hypothèse de réversibilité du flot. Nous développons deux méthodes, la première apparentée à une méthode de complément de Schur, la seconde basée sur une méthode de type Schwarz dont nous montrons la convergence pouvant être accélérée par la méthode d’Aitken dans le cadre linéaire. Afin d’accélérer la convergence de cette dernière dans le cadre non-linéaire, nous introduisons les techniques d’extrapolation et d’accélération de la convergence des suites non-linéaires. Nous montrons les avantages et les limites de ces techniques. Les résultats obtenus nous conduisent à développer l’accélération de la méthode de type Schwarz par une méthode de Newton. Enfin nous nous intéressons à l’étude de conditions de raccord non-linéaires adaptées à la décomposition de domaine de problèmes non-linéaires. Nous nous servons du formalisme hamiltonien à ports, issu du domaine de l’automatique, pour déduire les conditions de raccord dans le cadre l’équation de Saint-Venant et de l’équation de la chaleur non-linéaire. Après une étude analytique de la convergence de la DDM associée à ces conditions de transmission, nous proposons et étudions une formulation de Lagrangien augmenté sous l’hypothèse de séparabilité de la contrainte. / Complexification of multi-physics modeling leads to have to simulate systems of ordinary differential equations and algebraic differential equations with increasingly large numbers of unknowns and over large times of simulation. In addition the evolution of parallel computing architectures requires other ways of parallelization than the decomposition of system in subsystems. In this work, we propose to design domain decomposition methods in time for the resolution of EDO. We reformulate the initial value problem in a boundary values problem on the symmetrized time interval, under the assumption of reversibility of the flow. We develop two methods, the first connected with a Schur complement method, the second based on a Schwarz type method for which we show convergence, being able to be accelerated by the Aitken method within the linear framework. In order to accelerate the convergence of the latter within the non-linear framework, we introduce the techniques of extrapolation and of acceleration of the convergence of non-linear sequences. We show the advantages and the limits of these techniques. The obtained results lead us to develop the acceleration of the method of the type Schwarz by a Newton method. Finally we investigate non-linear matching conditions adapted to the domain decomposition of nonlinear problems. We make use of the port-Hamiltonian formalism, resulting from the control field, to deduce the matching conditions in the framework of the shallow-water equation and the non-linear heat equation. After an analytical study of the convergence of the DDM associated with these conditions of transmission, we propose and study a formulation of augmented Lagrangian under the assumption of separability of the constraint.

Cette thèse développe des estimations d’erreur a posteriori et critères d’arrêt pour les méthodes de décomposition de domaine avec des conditions de transmission de Robin optimisées entre les interfaces. Différents problèmes sont considérés: l’équation de Darcy stationnaire puis l’équation de la chaleur, discrétisées par les éléments finis mixtes avec un schéma de Galerkin discontinu de plus bas degré en temps pour le second cas. Pour l’équation de la chaleur, une méthode de décomposition de domaine globale en temps, avec mêmes ou différents pas de temps entre les différents sous domaines, est utilisée. Ce travail est finalement étendu à un modèle diphasique en utilisant une méthode de volumes finis centrés par maille en espace. Pour chaque modèle, un problème d’interface est résolu itérativement, où chaque itération nécessite la résolution d’un problème local dans chaque sous-domaine, et les informations sont ensuite transmises aux sous-domaines voisins. Pour les modèles instationnaires, les problèmes locaux dans les sous-domaines sont instationnaires et les données sont transmises par l’interface espace-temps. L’objectif de ce travail est, pour chaque modèle, de borner l’erreur entre la solution exacte et la solution approchée à chaque itération de l’algorithme de décomposition de domaine. Différentes composantes d’erreur en jeu de la méthode sont identifiées, dont celle de l’algorithme de décomposition de domaine, de façon à définir un critère d’arrêt efficace pour cette méthode. En particulier, pour l’équation de Darcy stationnaire, on bornera l’erreur par un estimateur de décomposition de domaine ainsi qu’un estimateur de discrétisation en espace. On ajoutera à la borne de l’erreur un estimateur de discrétisation en temps pour l’équation de la chaleur et pour le modèle diphasique. L’estimation a posteriori répose sur des techniques de reconstructions de pressions et de flux conformes respectivement dans les espaces H1 et H(div) et sur la résolution de problèmes locaux de Neumann dans des bandes autour des interfaces de chaque sous-domaine pour les flux. Ainsi, des critères pour arrêter les itérations de l’algorithme itératif de décomposition de domaine sont développés. Des simulations numériques pour des problèmes académiques ainsi qu’un problème plus réaliste basé sur des données industrielles sont présentées pour illustrer l’efficacité de ces techniques. En particulier, différents pas de temps entre les sous-domaines sont considérés pour cet exemple. / This work contributes to the developpement of a posteriori error estimates and stopping criteria for domain decomposition methods with optimized Robin transmission conditions on the interface between subdomains. We study several problems. First, we tackle the steady diffusion equation using the mixed finite element subdomain discretization. Then the heat equation using the mixed finite element method in space and the discontinuous Galerkin scheme of lowest order in time is investigated. For the heat equation, a global-in-time domain decomposition method is used for both conforming and nonconforming time grids allowing for different time steps in different subdomains. This work is then extended to a two-phase flow model using a finite volume scheme in space. For each model, the multidomain formulation can be rewritten as an interface problem which is solved iteratively. Here at each iteration, local subdomain problems are solved, and information is then transferred to the neighboring subdomains. For unsteady problems, the subdomain problems are time-dependent and information is transferred via a space-time interface. The aim of this work is to bound the error between the exact solution and the approximate solution at each iteration of the domain decomposition algorithm. Different error components, such as the domain decomposition error, are identified in order to define efficient stopping criteria for the domain decomposition algorithm. More precisely, for the steady diffusion problem, the error of the domain decomposition method and that of the discretization in space are estimated separately. In addition, the time error for the unsteady problems is identified. Our a posteriori estimates are based on the reconstruction techniques for pressures and fluxes respectively in the spaces H1 and H(div). For the fluxes, local Neumann problems in bands arround the interfaces extracted from the subdomains are solved. Consequently, an effective criterion to stop the domain decomposition iterations is developed. Numerical experiments, both academic and more realistic with industrial data, are shown to illustrate the efficiency of these techniques. In particular, different time steps in different subdomains for the industrial example are used.