In this dissertation we present an analysis for the discontinuous Galerkin discretization error of multi-dimensional first-order linear symmetric and symmetrizable hyperbolic systems of conservation laws. We explicitly write the leading term of the local DG error, which is spanned by Legendre polynomials of degree p and p+1 when p-th degree polynomial spaces are used for the solution. For special hyperbolic systems, where the coefficient matrices are nonsingular, we show that the leading term of the error is spanned by (p+1)-th degree Radau polynomials. We apply these asymptotic results to observe that projections of the error are pointwise O(h^p+2)-superconvergent in some cases and establish superconvergence results for some integrals of the error. We develop an efficient implicit residual-based a posteriori error estimation scheme by solving local finite element problems to compute estimates of the leading term of the discretization error. For smooth solutions we obtain error estimates that converge to the true error under mesh refinement. We first show these results for linear symmetric systems that satisfy certain assumptions, then for general linear symmetric systems. We further generalize these results to linear symmetrizable systems by considering an equivalent symmetric formulation, which requires us to make small modifications in the error estimation procedure. We also investigate the behavior of the discretization error when the Lax-Friedrichs numerical flux is used, and we construct asymptotically exact a posteriori error estimates. While no superconvergence results can be obtained for this flux, the error estimation results can be recovered in most cases. These error estimates are used to drive h- and p-adaptive algorithms and assess the numerical accuracy of the solution. We present computational results for different fluxes and several linear and nonlinear hyperbolic systems in one, two and three dimensions to validate our theory. Examples include the wave equation, Maxwell's equations, and the acoustic equation. / Ph. D.

Cette thèse porte sur l’étude des systèmes de lois de conservation couplés par une interface mobile. Un modèle dynamique d’un procédé d’extrusion obtenu à partir des bilans de masse, de taux d’humidité et d’énergie est proposé. Ce modèle exprime le transport de la matière et de la chaleur dans une extrudeuse par des systèmes d’équations hyperboliques définis sur deux domaines complémentaires variant dans le temps. L’évolution des domaines est dictée par une Equation aux Dérivées Ordinaires (EDO) issue du bilan de masse total dans une extrudeuse. Par le principe des applications contractantes l’existence et l’unicité de la solution pour cette classe de système sont prouvées. Le problème de stabilisation de l’interface mobile est aussi abordé en utilisation le formalisme des systèmes à retard. La méthode des caractéristiques permet de représenter le système composé des équations issues du bilan de masse par un système à retard sur l’entrée. Au moyen d’un contrôleur prédictif la position de l’interface est stabilisée autour d’un point équilibre. La dernière partie de ce travail est dédiée à l’étude des systèmes Hamiltoniens à ports frontière couplés par une interface mobile. Ces systèmes augmentés de variables couleur qui sont des fonctions caractéristiques du domaine peuvent s’exprimer comme des systèmes Hamiltoniens à ports frontière / This thesis is devoted to the analysis of Partial Differential Equations (PDEs) which are coupled through a moving interface. The motion of the interface obeys to an Ordinary Differential Equation (ODE) which arises from a conservation law. The first part of this thesis concerns the modelling of an extrusion process based on mass, moisture content and energy balances. These balances laws express heat and homogeneous material transport in an extruder by hyperbolic PDEs which are defined in complementary time-varying domains. The evolution of the coupled domains is given by an ODE which is derived from the conservation of mass in an extruder. In the second part of the manuscript, a mathematical analysis has been performed in order to prove the existence and the uniqueness of solution for such class of systems by mean of contraction mapping principle. The third part of the thesis concerns the transformation of an extrusion process mass balance equations into a particular input delay system framework using characteristics method. Then, the stabilization of the moving interface by a predictor-based controller has been proposed. Finally, an extension of the analysis of moving interface problems to a particular class of systems of conservations laws has been developed. Port-Hamiltonian formulation of systems of two conservation laws defined on two complementary time-varying intervals has been studied. It has been shown that the coupled system is a port-Hamiltonian system augmented with two variables being the characteristic functions of the two spatial domains

Le modèle multi-fluide permet de décrire par une approche Eulérienne les sprays polydispersés et apparaît donc comme une méthode indiquée pour les applications de combustion diphasique. Sa pertinence pour la simulation à l’échelle d’applications industrielles est évaluée dans ce travail, par sa mise en oeuvre dans des configurations bi-dimensionnelle et tri-dimensionnelle plus représentatives de ce type de simulations. Cette évaluation couple une étude de faisabilité en terme de coût de calcul avec une analyse de la précision obtenue, par des comparaisons avec les résultats de méthodes de références pour la description des sprays. Afin de définir une telle référence, une hiérarchisation des modèles de spray est proposée dans ce travail, soulignant les niveaux de modélisation associée aux diverses méthodes. Une première configuration d’écoulements tourbillonnaires est utilisée pour caractériser la méthode multi-fluide. L’étude de la structure mathématique du système de lois de conservation permet d’analyser la formation de singularités et de fournir les outils permettant d’évaluer leur impact sur la modélisation. Cette étude permet également de dériver un schéma numérique robuste et efficace pour des configurations bi- et tri-dimensionnelle. La description des dynamiques de gouttes conditionnées par la taille est évaluée dans ces configurations tourbillonnaires au moyen de comparaisons quantitatives, sur des champs instantanés, où le multi-fluide est confronté à une méthode Lagrangienne, ainsi qu’à des résultats expérimentaux. Afin d’évaluer le comportement de la méthode multi-fluide dans des configurations plus représentatives des problématiques industrielles, le solveur MUSES3D est développé, permettant, entre autres, une évaluation fine des méthodes de résolution des sprays. Une implémentation originale de la méthode multi-fluide, conciliant généricité et efficacité pour le calcul parallèle, est réalisée. Le couplage de ce solveur avec le code ASPHODELE, développé au CORIA, permet d’effectuer une évaluation opérationnelle des approches Euler/Lagrange et Euler/Euler pour la description des écoulements diphasiques à inclusions dispersées. Finalement, le comportement de la méthode multi-fluide dans des jets bi-dimensionnels et dans une turbulence homogène isotrope tri-dimensionnelle permet de montrer sa précision pour la description de la dynamique de sprays évaporant dans des configurations plus complexes. La résolution de la polydispersion du spray permet de décrire précisément la fraction massique de combustible en phase vapeur, un élément clé pour les applications de combustion. De plus, l’efficacité du calcul parallèle par décomposition de domaine avec la méthode multi-fluide permet d’envisager son utilisation à l’échelle d’applications industrielles. / The multi-fluid model, providing a Eulerian description of polydisperse sprays, appears as an interesting method for two-phase combustion applications. Its relevance as a numerical tool for industrial device simulations is evaluated in this work. This evaluation assesses the feasibility of multi-fluid simulations in terms of computational cost and analyzes their precision through comparisons with reference methods for spray resolution. In order to define such a reference, the link between the available methods for spray resolution is provided, highlighting their corresponding level of modeling. A first framework of 2-D vortical flows is used to assess the mathematical structure of the multi-fluid model governing system of equations. The link between the mathematical peculiarities and the physical modeling is provided, and a robust numerical scheme efficient for 2-D/3-D configurations is designed. This framework is also used to evaluate the multi-fluid description of evaporating spray sizeconditioned dynamics through quantitative, time-resolved, comparisons with a Lagrangian reference and with experimental data. In order to assess the multi-fluid efficiency in configurations more representative of industrial devices, a numerical solver is designed, providing a framework devoted to spray method evaluation. An original implementation of the multifluid method, combining genericity and efficiency in a parallel framework, is achieved. The coupling with a Eulerian/Lagrangian solver for dispersed two-phase flows, developed at CORIA, is conducted. It allows a precise evaluation of Euler/Lagrange versus Euler/Euler approaches, in terms of precision and computational cost. Finally, the behavior of the multi-fluid model is assessed in 2D-jets and 3-D Homogeneous Isotropic Turbulence. It illustrates the ability of the method to capture evaporating spray dynamics in more complex configurations. The method is shown to describe accurately the fuel vapor mass fraction, a key issue for combustion applications. Furthermore, the method is shown to be efficient in domain decomposition parallel computing framework, a key issue for simulations at the scale of industrial devices.