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 est consacrée à la modélisation de la ventilation mécanique chez l'humain et à l'analyse numérique des systèmes en découlant. Des simulations directes d'écoulement d'air dans l'ensemble des voies aériennes étant impossibles (maillages indisponibles et géométrie trop complexe), il est nécessaire de considérer un domaine d'intérêt réduit, qui implique de travailler dans une géométrie tronquée, comportant des frontières artificielles ou encore de considérer des modèles réduits simples mais représentatifs. Si on cherche à effectuer des simulations numériques 3D où l'écoulement du fluide est décrit par les équations de Navier-Stokes, différentes problématiques sont soulevées :- Si on considère que la ventilation est la conséquence de différences de pression, les conditions aux limites associées sont des conditions de type Neumann. Cela aboutit à des questions théoriques en terme d'existence et d'unicité de solution et à des questions numériques en terme de choix de schémas et de méthodes adaptées.- Lorsque l'on travaille dans un domaine tronqué, il peut être nécessaire de prendre en compte les phénomènes non décrits grâce à des modèles réduits appropriés. Ici nous considérons des modèles 0D. Ces couplages 3D/0D sont à l'origine d'instabilités numériques qu'on étudie mathématiquement et numériquement dans ce manuscrit. Par ailleurs, lorsqu'on s'intéresse à des régimes de respiration forcée, les modèles usuels linéaires sont invalidés par les expériences. Afin d'observer les différences entre les résultats expérimentaux et numériques, il est nécessaire de prendre en compte plusieurs types de non linéarités, comme la déformation du domaine ou les phénomènes de type Bernoulli. Une approche par modèles réduits est adoptée dans ce travail.Pour finir, on a cherché à valider les modèles obtenus en comparant des résultats numériques et des résultats expérimentaux dans le cadre d'un travail interdisciplinaire.Parvenir à modéliser et simuler ces écoulements permet de mieux comprendre les phénomènes et paramètres qui entrent en jeu lors de pathologies (asthme, emphysème...). Un des objectifs à moyen terme est d'étudier l'influence du mélange hélium-oxygène sur le dépôt d'aérosol, toujours dans le cadre du travail interdisciplinaire. A plus long terme, l'application de ces modèles à des situations pathologiques pourrait permettre de construire des outils d'aide à la décision dans le domaine médical (compréhension de la pathologie, optimisation de thérapie...). / In this thesis, we study the modelling of the human mecanical ventilation and the numerical analysis of linked systems. Direct simulations of air flow in the whole airways are impossible (complex geometry, unavailable meshes). Then a reduced area of interest can be considered, working with reduced geometries and artificial boundaries. One can also use reduced models, simple but realistic. If one try to make 3D numerical simulations where the fluid flow is described by the Navier-Stokes equations, various issues are raised:- If we consider that ventilation is the result of pressure drops, the associated boundary conditions are Neumann conditions. It leads to theoretical questions in terms of existence and uniqueness of solution and numerical issues in terms of scheme choice and appropriate numerical methods.- When working in a truncated domain, it may be necessary to take into account non-described phenomena with appropriate models. Here we consider 0D models. These 3D/0D couplings imply numerical instabilities that we mathematically and numerically study in this thesis.Furthermore, when we focus on forced breathing, linear usual models are invalidated by experiments. In order to observe the differences between the experimental and numerical results, it is necessary to take into account several types of non-linearities, such as deformation of the domain or the Bernoulli phenomenon. A reduced model approach is adopted in this work. Finally, we sought to validate the obtained models by comparing numerical and experimental results in the context of interdisciplinary work.Achieving model and simulate these flows allow to better understand phenomena and parameters that come into play in diseases (asthma, emphysema ...). A medium-term objective is to study the influence of helium-oxygen mixture in the aerosol deposition. In the longer term, the application of these models to pathological situations could afford to build decision support tools in the medical field (understanding of pathology, therapy optimization ...).

There is no abstract

Les travaux ont traité de l'ajustement des efforts de calcul vibro-acoustique dans le cadre d'une étude de qualité sonore sur un système vibro-acoustique dès le stade de la conception. Cette recherche d'ajustement a été appliquée au système vibro-acoustique constitué d'une plaque vibrante couplée à une cavité amortie. Pour un calcul fréquentiel, un premier pas dans ce sens demandait à ajuster deux paramètres de simulation importants : la fréquence maximale de calcul, puis le pas fréquentiel de calcul. En termes de contraintes perceptives, il a été cherché pour ces paramètres une valeur ajustée permettant de conserver les tendances qualitatives originellement établies à partir d'une évaluation auditive portant sur les sons réels rayonnés par la plaque dans la cavité. Les sons réels ont été collectés via l'expérimentation sur un banc d'essai constitué d'une plaque en acier couplée à une cavité parallélépipédique ; ces sons correspondaient au son rayonné par la plaque vibrante à l'intérieur de la cavité, enregistré dans diverses configurations structurales du système. Ces configurations impliquaient différentes combinaisons des modalités prises par les 3 paramètres structuraux variables retenus : l'épaisseur de la plaque, les conditions de serrage de la plaque et les propriétés en absorption de la cavité. Les sons enregistrés ont ensuite été soumis par paires à un jury d'auditeurs auxquels il était demandé d'évaluer leur dissimilarité et de donner un jugement de préférence. A partir de l'analyse des informations recueillies, des tendances qualitatives de référence ont été établies, en termes d'espace perceptif de dissimilarité, d'espace perceptif de préférence et de classement de préférence. Dès lors, en travaillant à partir des sons enregistrés, le processus d'ajustement a consisté, i) à déterminer, à partir d'un filtrage passe-bas appliqué aux sons, la fréquence de coupure ajustée qui permettait de préserver les tendances qualitatives de référence, ii) à rechercher le pas fréquentiel ajusté, à imposer aux spectres des sons avant synthèse, qui permettait de maintenir les tendances qualitatives précédemment établies.

Les travaux de la thèse portent sur la résolution des équations de Stokes, d'abord avec des conditions au bord portant sur la composante normale du champ de vitesse et la composante tangentielle du tourbillon, ensuite avec des conditions au bord portant sur la pression et la composante tangentielle du champ de vitesse. Dans chaque cas nous démontrons l'existence, l'unicité et la régularité de la solution. Nous traitons aussi le cas de solutions très faibles, par dualité. Le cadre fonctionnel que nous avons choisi est celui des espaces de Banach du type H(div) et H(rot) ou l'intersection des deux, basés sur l'espace Lp , avec 1 < p < ∞. En particulier, on se place dans des domaines non simplement connexes, avec des frontières non connexes. Nous nous intéressons en premier lieu à l'obtention d'inégalités de Sobolev pour des champs de vecteurs u ∈ Lp (Ω). Dans un second temps, nous établissons des résultats d'existence pour les potentiels vecteurs avec diverses conditions aux limites. Ceci nous permet d'abord d'effectuer des décompositions de type Helmholtz et ensuite de démontrer des conditions Inf − Sup lorsque la forme bilinéaire est un produit de rotationnels. Ces conditions aux limites font que l'équation de la pression est indépendante des autres variables. C'est la raison pour laquelle nous sommes naturellement conduit à étudier les problèmes elliptiques qui se traduisent par les systèmes de Stokes sans la pression. La résolution de ces problèmes se fait au moyen des Conditions Inf − Sup qui jouent un rôle clef pour établir l'existence et l'unicité de solutions. Nous donnons une applications aux systèmes de Navier-Stokes, où on obtient l'existence d'une solution en effectuant un point fi xe autour du problème d'Oseen. Enfi n, deux méthodes numériques sont proposées pour approcher le problème de Stokes. Nous analysons d'abord une méthode de Nitsche et puis une méthode de Galerkin discontinu. Quelques résultats numériques de convergence sont décrits qui sont parfaitement cohérents avec l'analyse.

Nous proposons une approche spectrale permettant d'aborder des problèmes d'interpolation à noyaux dont la résolution numérique n'est pas directement envisageable. Un tel cas de figure se produit en particulier lorsque le nombre de données est infini. Nous considérons dans un premier temps le cadre de l'interpolation optimale dans les sous-espaces hilbertiens. Pour un problème donné, un opérateur intégral est défini à partir du noyau sous-jacent et d'une paramétrisation de l'ensemble des données basée sur un espace mesuré. La décomposition spectrale de l'opérateur est utilisée afin d'obtenir une formule de représentation pour l'interpolateur optimal et son approximation est alors rendu possible par troncature du spectre. Le choix de la mesure induit une fonction d'importance sur l'ensemble des données qui se traduit, en cas d'approximation, par une plus ou moins grande précision dans le rendu des données. Nous montrons à titre d'exemple comment cette approche peut être utilisée afin de rendre compte de contraintes de type "conditions aux limites" dans les modèles d'interpolation à noyaux. Le problème du conditionnement des processus gaussiens est également étudié dans ce contexte. Nous abordons enfin dans la dernière partie de notre manuscrit la notion de noyaux conditionnellement positifs. Nous proposons la définition générale de noyaux symétriques conditionnellement positifs relatifs à une espace de référence donné et développons la théorie des sous-espaces semi-hilbertiens leur étant associés. Nous étudions finalement la théorie de l'interpolation optimale dans cette classe d'espaces. / We propose a spectral approach for the resolution of kernel-based interpolation problems of which numerical solution can not be directly computed. Such a situation occurs in particular when the number of data is infinite. We first consider optimal interpolation in Hilbert subspaces. For a given problem, an integral operator is defined from the underlying kernel and a parameterization of the data set based on a measurable space. The spectral decomposition of the operator is used in order to obtain a representation formula for the optimal interpolator and spectral truncation allows its approximation. The choice of the measure on the parameters space introduces a hierarchy onto the data set which allows a tunable precision of the approximation. As an example, we show how this methodology can be used in order to enforce boundary conditions in kernel-based interpolation models. The Gaussian processes conditioning problem is also studied in this context. The last part of this thesis is devoted to the notion of conditionally positive kernels. We propose a general definition of symmetric conditionally positive kernels relative to a given space and exposed the associated theory of semi-Hilbert subspaces. We finally study the optimal interpolation problem in such spaces.

Les écoulements en eau peu profonde se rencontrent dans de nombreuses situations d’intérêts : écoulements de rivières et dans les lacs, mais aussi dans les mers et océans (courants de marée, tsunami, etc.). Ils sont modélisés par un système d’équations aux dérivées partielles, où les inconnues sont la vitesse de l’écoulement et la hauteur d’eau. On peut supposer que la composante verticale de la vitesse est petite devant les composantes horizontales et que ces dernières sont indépendantes de la profondeur. Le modèle est alors donné par les équations de shallow water (SWEs). Cette thèse se concentre sur la conception d’une nouvelle technique d’interaction de plusieurs grilles imbriquées pour modèle en eau peu profonde en utilisant des méthodes numériques. La première partie de cette thèse comprend, La dérivation complète de ces équations à partir des équations de Navier- Stokes est expliquée. Etudier le développement et l’évaluation des méthodes numériques en utilisant des méthodes de différences finies et plusieurs exemples numériques sont appliqués utilisant la condition initiale du niveau gaussien pour 2DSWEs. Dans la deuxième partie de la thèse, nous sommes intéressés à proposer une nouvelle technique d’interaction de plusieurs grilles imbriquées pour résoudre les modèles océaniques en utilisant quatre choix des opérateurs de restriction avec des résultats de haute précision. Notre travail s’est concentré sur la résolution numérique de SWE par grilles imbriquées. A chaque niveau de résolution, nous avons utilisé une méthode classique de différences finies sur une grille C d’Arakawa, avec un schéma de leapfrog complété par un filtre d’Asselin. Afin de pouvoir affiner les calculs dans les régions perturbées et de les alléger dans les zones calmes, nous avons considéré plusieurs niveaux de résolution en utilisant des grilles imbriquées. Ceci permet d’augmenter considérablement le rapport performance de la méthode, à condition de régler efficacement les interactions (spatiales et temporelles) entre les grilles. Dans la troisième partie de cette thèse, plusieurs exemples numéériques sont testés pour 2DSWE avec imbriqués 3:1 et 5:1. Finalement, la quatrième partie de ce travail, certaines applications de grilles imbriquées pour le modèle tsunami sont présentées. / Most flows in the rivers, seas, and ocean are shallow water flow in which the horizontal length andvelocity scales are much larger than the vertical ones. The mathematical formulation of these flows, so called shallow water equations (SWEs). These equations are a system of hyperbolic partial differentialequations and they are effective for many physical phenomena in the oceans, coastal regions, riversand canals. This thesis focuses on the design of a new two-way interaction technique for multiple nested grids 2DSWEs using the numerical methods. The first part of this thesis includes, proposing several ways to develop the derivation of shallow water model. The complete derivation of this system from Navier-Stokes equations is explained. Studying the development and evaluation of numerical methods by suggesting new spatial and temporal discretization techniques in a standard C-grid using an explicit finite difference method in space and leapfrog with Robert-Asselin filter in time which are effective for modeling in oceanic and atmospheric flows. Several numerical examples for this model using Gaussian level initial condition are implemented in order to validate the efficiency of the proposed method. In the second part of our work, we are interested to propose a new two-way interaction technique for multiple nested grids to solve ocean models using four choices of higher restriction operators (update schemes) for the free surface elevation and velocities with high accuracy results. Our work focused on the numerical resolution of SWEs by nested grids. At each level of resolution, we used explicit finite differences methods on Arakawa C-grid. In order to be able to refine the calculations in troubled regions and move them into quiet areas, we have considered several levels of resolution using nested grids. This makes it possible to considerably increase the performance ratio of the method, provided that the interactions (spatial and temporal) between the grids are effectively controlled. In the third part of this thesis, several numerical examples are tested to show and verify twoway interaction technique for multiple nested grids of shallow water models can works efficiently over different periods of time with nesting 3:1 and 5:1 at multiple levels. Some examples for multiple nested grids of the tsunami model with nesting 5:1 using moving boundary conditions are tested in the fourth part of this work.

Ce texte de synthèse a pour but de présenter l'´évolution de mes recherches postèrieures à ma thèse. Ce travail s'articule autour de plusieurs axes de recherche dans le cadre des équations aux dérivées partielles non linéaires et en particulier des lois de conservation.<br />Il s'inscrit dans l'étude des problèmes hyperboliques, des problème mixtes et des équations cinétiques. Les domaines d'application sont la mécanique des fluides ou du solide, la propagation de composants chimiques, l'électromagnétisme, l'optique.<br />Mon activité concerne d'abord la modélisation de phénomènes physiques ou chimiques sous forme d'équations aux dérivées partielles non linéaires telles que les équations de Bloch, Korteweg, Navier-Stokes, Saint-Venant, puis vient l'étude mathématique de ces équations à travers les<br />problèmes d'existence, d'unicité, de régularité avec éventuellement la mise au point de méthodes numériques de résolution.<br /> Ce document est divisé en une introduction générale et trois chapitres qui concernent respectivement les systèmes hyperboliques avec conditions aux limites et la chromatographie, les problèmes d'analyse asymptotique et enfin les méthodes cinétiques.<br />Dans chaque partie, un historique et une présentation des différents résultats mathématiques sont faits et quelques problèmes ouverts sont donnés.

L'objet de cette thèse est de présenter une étude du comportement thermomécanique des disques de frein automobiles pour la prédiction de leur tenue en fatigue. La stratégie de calcul numérique est repose sur le code de calcul Ansys v.11. Ce dernier qui est basé sur la méthode des éléments finis et qui possède des algorithmes de gestion du contact avec frottement est utilisé pour simuler dans l'application du freinage le comportement du mécanisme malgré son complexité. Dans un premier temps, est présentée une analyse des phénomènes thermiques opérant dans un disque de frein en service (flux de chaleur généré par frottement, gradients thermiques élevés, élévation de température). Cette modélisation est effectuée en tenant compte l'influence d'un certains nombre de paramètre tel que le type de freinage, le mode de refroidissement, les matériaux de conceptions. Ensuite, une étude purement mécanique du contact sec entre le disque et plaquettes est développée avec une bonne prédiction devient un enjeu majeur pour les industriels tout en modélisant le chargement et les conditions aux limites autours du disque. Nous avons utilisé le même code de calcul pour visualiser les déplacements, les déformations globales dans le disque, les contraintes de cisaillement, les contraintes de Von Mises et les outils de contact des plaquettes tout en effectuant une étude paramétrique telle que ( le module de Young des plaquettes, le coefficient de frottement , le type de chargement , la vitesse de rotation du disque ,..) pour voir sa sensibilité sur les résultats de calcul. Ainsi, les analyses faites sur le comportement thermique et mécanique que de ces prototypes montrent que ces types de solutions technologiques représentent de réelles pistes d'amélioration qui répond au besoin de l'ingénieur en charge de la conception des disques de frein.