Dans cette thèse, nous nous intéressons à un modèle simplifié de corrosion issu du modèle ''Diffusion Poisson Coupled Model'' (DPCM). Nous analysons de manière approfondie le schéma numérique qui a été implémenté dans le code CALIPSO utilisé par l'ANDRA. Il est de type Euler implicite en temps et volumes finis en espace, avec des flux de Scharfetter-Gummel. Nous étudions notamment la convergence de ce schéma ainsi que son comportement asymptotique en différentes limites de paramètres. Enfin, nous explorons différentes possibilités pour augmenter l'ordre en temps. / In this PhD thesis, we are interested in a simplified corrosion model derived from the Diffusion Poisson Coupled Model (DPCM). We analyze the numerical scheme implemented in the CALIPSO code used by the French nuclear waste management agency ANDRA. It is a backward Euler scheme in time and a finite volume scheme in space, with Schafetter-Gummel approximation of the convection-diffusion fluxes. We study the convergence of this scheme and its asymptotic behavior for different limits of parameters. Finally, we compare several higher order schemes in time.

Thèse diffusée initialement dans le cadre d'un projet pilote des Presses de l'Université de Montréal/Centre d'édition numérique UdeM (1997-2008) avec l'autorisation de l'auteur.

Thèse diffusée initialement dans le cadre d'un projet pilote des Presses de l'Université de Montréal/Centre d'édition numérique UdeM (1997-2008) avec l'autorisation de l'auteur.

Les écoulements diphasiques dans les centrales de type réacteur à eau pressurisée appartiennent à des régimes très variés allant du faible nombre de Mach jusqu'aux ondes de chocs. Calculer des solutions approchées précises de ces écoulements peut s'avérer délicat dans certains régimes. On s'intéresse dans cette thèse à la conception et à l'étude de méthodes numériques robustes et stables à grand pas de temps, capables de calculer des solutions approchées précises quel que soit le régime d'écoulement, y compris sur maillage grossier. Une stratégie pour construire de tels schémas consiste à : utiliser un schéma semi implicite basé sur un splitting d’opérateurs pour séparer la résolution approchée des phénomènes rapides de celles des phénomènes lents ; corriger les flux numériques afin d’améliorer la précision du schéma dans certains régimes. Deux approches sont utilisées pour analyser la capacité du schéma numérique à gérer plusieurs régimes d'écoulement. L’approche des schémas asymptotic preserving est utilisée pour traiter le système de la dynamique des gaz avec termes sources raides. On utilise ensuite la notion de schéma tout-régime pour le système de la dynamique des gaz et les systèmes diphasiques homogénéisés HRM et HEM à bas nombre de Mach. Des propriétés garantissant la stabilité et la robustesse des schémas ont été obtenues, et en particulier des inégalités d'entropie discrètes. L'implémentation de ces méthodes a permis de mener des expériences numériques en 1D et 2D sur maillage non structuré qui confirment le gain en précision et en temps de calcul des schémas asymptotic preserving et tout-régime ainsi construits par rapport à des schémas numériques classiques. / Two-phase flows in Pressurized Water Reactors belong to a wide range of Mach number flows. Computing accurate approximate solutions of those flows may be challenging from a numerical point of view as classical finite volume methods are too diffusive in the low Mach regime. In this thesis, we are interested in designing and studying some robust numerical schemes that are stable for large time steps and accurate even on coarse meshes for a wide range of flow regimes. An important feature is the strategy to construct those schemes. We use a mixed implicit-explicit strategy based on an operator splitting to solve fast and slow phenomena separately. Then, we introduce a modification of a Suliciu type relaxation scheme to improve the accuracy of the numerical scheme in some regime of interest. Two approaches have been used to assess the ability of our numerical schemes to deal with a wide range of flow regimes. The first approach, based on the asymptotic preserving property, has been used for the gas dynamics equations with stiff source terms. The second approach, based on the all-regime property, has been used for the gas dynamics equations and the homogeneous two-phase flows models HRM and HEM in the low Mach regime. We obtained some robustness and stability properties for our numerical schemes. In particular, some discrete entropy inequalities are shown. Numerical evidences, in 1D and in 2D on unstructured meshes, assess the gain in term of accuracy and CPU time of those asymptotic preserving and all-regime numerical schemes in comparison with classical finite volume methods.

Dans cette thèse, nous développons des méthodes analytiques et numériques pour capturer les dynamiques asymptotiques de problèmes issus de la physique des plasmas et de la modélisation des mouvements collectifs dans les populations animales. Dans une première partie, nous présentons une méthode numérique Particle-In-Cell (PIC) pour le système Vlasov-Poisson préservant l'asymptotique quasi-neutre. Dans une seconde partie, nous étudions la limite macroscopique d'un modèle de Vicsek décrivant des interactions d'alignement entre deux populations, une population à l'arrêt et une population en mouvement. Nous sélectionnons ensuite un schéma numérique pour capturer les solutions du modèle macroscopique de Vicsek correspondant à la dynamique particulaire sous-jacente. La troisième partie est dédiée à l'étude des transitions compressible-incompressible apparaissant sous l'effet d'une contrainte de congestion dans un modèle macroscopique de déplacement collectif. Des schémas numériques préservant l'asymptotique de congestion sont ensuite mis au point pour le système d'Euler avec une contrainte de densité maximale.

Cette étude porte sur la modélisation et la simulation numérique d'un piano, en domaine temporel, par modèles phy- siques. Nous souhaitons rendre compte du comportement vibratoire et acoustique du piano, en prenant en compte les éléments principaux qui contribuent à la production du son. La table d'harmonie est modélisée par une équation bidimensionnelle de plaque épaisse, le système de Reissner Mindlin, pour un matériau orthotrope et hétérogène, dont l'amortissement dépend de la fréquence. Grâce aux équations de la vibroacoustique, la table rayonne dans l'air, dans lequel on souhaite calculer le champ acoustique complet autour de la ceinture du piano, que l'on suppose rigide. La table d'harmonie est d'autre part sollicitée par les cordes, à travers le chevalet où elles présentent un léger angle par rapport au plan horizontal. Chaque corde est modélisée par un système d'équations monodimensionnelles amorties dans lequel on prend en compte non seulement les ondes transversales excitées par le marteau, mais aussi la raideur à travers les ondes de cisaillement, ainsi que le couplage avec les ondes longi- tudinales provenant de la prise en compte des non linéarités géométriques. Le marteau est lancé avec une vitesse initiale vers un chœur de cordes, contre lequel il s'écrase avant d'être repoussé par les cordes. La force d'interaction dépend de façon non linéaire de l'écrasement du marteau.Le modèle complet de piano, que l'on souhaite résoudre numériquement, consiste donc en un système couplé d'équations aux dérivées partielles, dont chacune revêt des difficultés de nature différente : la corde est régie par un système d'équations non linéaires, la table d'harmonie est soumise à un amortissement dépendant de la fréquence, la propagation acoustique requiert un très grand nombre d'inconnues; auxquelles s'ajoute la difficulté inhérente aux couplages. D'une part, la stabilité numérique du schéma discret peut être compromise par la présence d'équations non linéaires et de nombreux couplages. Une méthode efficace pour garantir cette stabilité a priori est de construire un schéma qui conserve, ou dissipe, un équivalent discret de l'énergie physique d'un pas de temps au suivant. Une contribution majeure de ce travail a été de développer des schémas préservant une énergie discrète pour une classe de systèmes non linéaires dans laquelle s'inscrit le modèle de corde. D'autre part, afin d'augmenter l'efficacité de la méthode et de réduire le coût des calculs numériques, il est souhaitable de mettre à jour de façon découplée les inconnues liées aux différentes parties du problème, sur lesquelles la discrétisation en temps est faite de façon différente, afin de s'adapter aux spécificités de chacune. L'introduction de multiplicateurs de Lagrange nous permet de réaliser ce découplage artificiel grâce à des compléments de Schur adaptés. L'utilisation du code de calcul en situation réaliste montre le potentiel d'une telle modélisation d'un piano complet en domaine temporel. Au delà de très bien reproduire les mesures, il est possible d'étudier l'influence de certains phénomènes physiques (corde raide, non linéaire), de la géométrie ou encore des matériaux utilisés sur le comportement vibratoire général du piano, et sur le son en particulier. L'enrichissement spectral, ainsi que l'apparition des " partiels fantômes " et du précurseur non linéaire sont clairement mis en évidence pour les grandes amplitudes de jeu, soulignant l'intérêt de notre approche dans la compréhension du fonctionnement de l'instrument.

Cette étude porte sur la modélisation et la simulation numérique d'un piano, en domaine temporel, par modèles phy- siques. Nous souhaitons rendre compte du comportement vibratoire et acoustique du piano, en prenant en compte les éléments principaux qui contribuent à la production du son. La table d'harmonie est modélisée par une équation bidimensionnelle de plaque épaisse, le système de Reissner Mindlin, pour un matériau orthotrope et hétérogène, dont l'amortissement dépend de la fréquence. Grâce aux équations de la vibroacoustique, la table rayonne dans l'air, dans lequel on souhaite calculer le champ acoustique complet autour de la ceinture du piano, que l'on suppose rigide. La table d'harmonie est d'autre part sollicitée par les cordes, à travers le chevalet où elles présentent un léger angle par rapport au plan horizontal. Chaque corde est modélisée par un système d'équations monodimensionnelles amorties dans lequel on prend en compte non seulement les ondes transversales excitées par le marteau, mais aussi la raideur à travers les ondes de cisaillement, ainsi que le couplage avec les ondes longi- tudinales provenant de la prise en compte des non linéarités géométriques. Le marteau est lancé avec une vitesse initiale vers un chœur de cordes, contre lequel il s'écrase avant d'être repoussé par les cordes. La force d'interaction dépend de façon non linéaire de l'écrasement du marteau.Le modèle complet de piano, que l'on souhaite résoudre numériquement, consiste donc en un système couplé d'équations aux dérivées partielles, dont chacune revêt des difficultés de nature différente : la corde est régie par un système d'équations non linéaires, la table d'harmonie est soumise à un amortissement dépendant de la fréquence, la propagation acoustique requiert un très grand nombre d'inconnues; auxquelles s'ajoute la difficulté inhérente aux couplages. D'une part, la stabilité numérique du schéma discret peut être compromise par la présence d'équations non linéaires et de nombreux couplages. Une méthode efficace pour garantir cette stabilité a priori est de construire un schéma qui conserve, ou dissipe, un équivalent discret de l'énergie physique d'un pas de temps au suivant. Une contribution majeure de ce travail a été de développer des schémas préservant une énergie discrète pour une classe de systèmes non linéaires dans laquelle s'inscrit le modèle de corde. D'autre part, afin d'augmenter l'efficacité de la méthode et de réduire le coût des calculs numériques, il est souhaitable de mettre à jour de façon découplée les inconnues liées aux différentes parties du problème, sur lesquelles la discrétisation en temps est faite de façon différente, afin de s'adapter aux spécificités de chacune. L'introduction de multiplicateurs de Lagrange nous permet de réaliser ce découplage artificiel grâce à des compléments de Schur adaptés. L'utilisation du code de calcul en situation réaliste montre le potentiel d'une telle modélisation d'un piano complet en domaine temporel. Au delà de très bien reproduire les mesures, il est possible d'étudier l'influence de certains phénomènes physiques (corde raide, non linéaire), de la géométrie ou encore des matériaux utilisés sur le comportement vibratoire général du piano, et sur le son en particulier. L'enrichissement spectral, ainsi que l'apparition des " partiels fantômes " et du précurseur non linéaire sont clairement mis en évidence pour les grandes amplitudes de jeu, soulignant l'intérêt de notre approche dans la compréhension du fonctionnement de l'instrument.

Cette thèse porte sur la modélisation des ondes élastodynamiques dans deux situations particulières qui pénalisent les méthodes numériques utilisées pour simuler ces phénomènes. Dans la première partie, on se place dans le cas où les ondes de pression (ondes P) se propagent à une vitesse beaucoup plus grande que celle des ondes de cisaillement (ondes S). Les modèles numériques utilisés habituellement pour traiter cette configuration sont pénalisés par la plus petite vitesse qui dicte le choix du pas du schéma. Nous proposons ici un schéma qui découple numériquement, dans le volume, les ondes P et les ondes S, pour deux types de conditions de bord en utilisant la décomposition du déplacement en potentiels de Lamé, en deux dimensions. Les conditions aux limites de Dirichlet homogènes, qui sont des conditions essentielles pour la formulation classique en déplacement, deviennent des conditions naturelles, mais non standard, pour la formulation en potentiels qui se présente comme un système de deux équations d'ondes couplées par les conditions aux limites. Cette formulation préserve une énergie équivalente à l'énergie élastodynamique. Nous construisons un schéma éléments finis en espace et utilisons un thêta-schéma en temps sur les termes de bord afin de ne pas pénaliser la CFL et mener à une condition sur le pas de temps indépendante des termes de couplage au bord. Ce schéma préserve une énergie discrète. Le cas des conditions de surface libre mène à des instabilités. Nous les avons traitées comme des perturbations des conditions de Dirichlet, ce qui permet d'obtenir de bons résultats dans le domaine fréquentiel mais donne naissance à de sévères instabilités après discrétisation en temps. La seconde partie de la thèse est consacrée à la construction, l'analyse et la validation de conditions de transmission effectives (CTE) à travers une couche mince de matériau homogène et isotrope d'épaisseur constante h. Ici, la finesse de la couche affecte les schémas explicites usuels car le maillage de la couche avec des éléments suffisamment petits entraîne une diminution analogue du pas de temps critique via la condition CFL, tandis que l'on espère avec les CTE obtenir un pas de temps indépendant de l'épaisseur de la couche. Une analyse complète du cas de la bande mince rectiligne est donnée en deux et trois dimensions. Les conditions obtenues sont stables via la conservation d'une énergie et l'ordre de l'erreur d'approximation par rapport à l'épaisseur de la couche pour les conditions d'ordre 2 est de O(h^3). Des résultats numériques sont présentés pour les configurations bi et tridimensionnelles, ils valident les résultats de stabilité, d'estimation d'erreur et de conditions de stabilité de schémas en temps proposés, qui sont des modifications du schéma explicite utilisé en l'absence de couche mince. Enfin, le traitement d'une couche curviligne est effectué dans le cas bidimensionnel. Sa stabilité est à nouveau vérifiée par conservation d'énergie et des résultats numériques sont également présentés.

Nous étudions plusieurs aspects d'équations aux dérivées partielles multi-échelles. Pour trois exemples, la présence de multiples échelles, spatiales ou temporelles, motive un travail de modélisation mathématique ou constitue un enjeu de discrétisation. La première partie est consacrée à la construction et l'étude d'un système multicouche de type Saint-Venant pour décrire un fluide à surface libre (océan). Son obtention s'appuie sur l'analyse des échelles spatiales, précisément l'hypothèse " eau peu profonde ". Nous justifions nos équations à partir du modèle primitif et montrons un résultat d'existence locale de solution. Puis nous proposons un schéma volumes finis et des simulations numériques. Nous étudions ensuite un problème hyperbolique de relaxation, inspiré de la théorie cinétique des gaz. Nous construisons un schéma numérique via une stratégie préservant l'asymptotique : nous montrons sa convergence pour toute valeur du paramètre de relaxation, ainsi que sa consistance avec le problème à l'équilibre local. Des estimations d'erreurs sont établies et des simulations numériques sont présentées. Enfin, nous étudions un problème d'écoulement sanguin dans une artère avec stent, modélisé par un système de Stokes dans un domaine contenant une petite rugosité périodique (géométrie double échelle). Pour éviter une discrétisation coûteuse du domaine rugueux (l'artère stentée), nous formulons un ansatz de développement de la solution type Chapman-Enskog, et obtenons une loi de paroi implicite sur le bord du domaine lisse (artère seule). Nous montrons des estimations d'erreurs et des simulations numériques

La migration de données entre des sources aux schémas hétérogènes est un domaine en pleine croissance avec l'augmentation de la quantité de données en accès libre, et le regroupement des données à des fins d'apprentissage automatisé et de fouilles. Cependant, la description du processus de transformation des données d'une instance source vers une instance définie sur un schéma différent est un processus complexe même pour un utilisateur expert dans ce domaine. Cette thèse aborde le problème de la définition de mapping par un utilisateur non expert dans le domaine de la migration de données, ainsi que la vérification du respect par ce mapping des contraintes d'accès ayant été définies sur les données sources. Pour cela, dans un premier temps nous proposons un système dans lequel l'utilisateur fournit un ensemble de petits exemples de ses données, et est amené à répondre à des questions booléennes simples afin de générer un mapping correspondant à ses besoins. Dans un second temps, nous proposons un système permettant de réécrire le mapping produit de manière à assurer qu'il respecte un ensemble de vues de contrôle d'accès définis sur le schéma source du mapping. Plus précisément, le premier grand axe de cette thèse est la formalisation du problème de la définition interactive de mappings, ainsi que la description d'un cadre formel pour la résolution de celui-ci. Cette approche formelle pour la résolution du problème de définition interactive de mappings est accompagnée de preuves de bonnes propriétés. A la suite de cela, basés sur le cadre formel défini précédemment, nous proposons des algorithmes permettant de résoudre efficacement ce problème en pratique. Ces algorithmes visent à réduire le nombre de questions auxquelles l'utilisateur doit répondre afin d'obtenir un mapping correspondant à ces besoins. Pour cela, les mappings possibles sont ordonnés dans des structures de treillis imbriqués, afin de permettre un élagage efficace de l'espace des mappings à explorer. Nous proposons également une extension de cette approche à l'utilisation de contraintes d'intégrité afin d'améliorer l’efficacité de l'élagage. Le second axe majeur vise à proposer un processus de réécriture de mapping qui, étant donné un ensemble de vues de contrôle d'accès de référence, permet d'assurer que le mapping réécrit ne laisse l'accès à aucune information n'étant pas accessible via les vues de contrôle d'accès. Pour cela, nous définissons un protocole de contrôle d'accès permettant de visualiser les informations accessibles ou non à travers un ensemble de vues de contrôle d'accès. Ensuite, nous décrivons un ensemble d'algorithmes permettant la réécriture d'un mapping en un mapping sûr vis-à-vis d'un ensemble de vues de contrôle d'accès. Comme précédemment, cette approche est complétée de preuves de bonnes propriétés. Afin de réduire le nombre d'interactions nécessaires avec l'utilisateur lors de la réécriture d'un mapping, une approche permettant l'apprentissage des préférences de l'utilisateur est proposée, cela afin de permettre le choix entre un processus interactif ou automatique. L'ensemble des algorithmes décrit dans cette thèse ont fait l'objet d'un prototypage et les expériences réalisées sur ceux-ci sont présentées dans cette thèse / Data exchange between sources over heterogeneous schemas is an ever-growing field of study with the increased availability of data, oftentimes available in open access, and the pooling of such data for data mining or learning purposes. However, the description of the data exchange process from a source to a target instance defined over a different schema is a cumbersome task, even for users acquainted with data exchange. In this thesis, we address the problem of allowing a non-expert user to spec- ify a source-to-target mapping, and the problem of ensuring that the specified mapping does not leak information forbidden by the security policies defined over the source. To do so, we first provide an interactive process in which users provide small examples of their data, and answer simple boolean questions in order to specify their intended mapping. Then, we provide another process to rewrite this mapping in order to ensure its safety with respect to the source policy views. As such, the first main contribution of this thesis is to provide a formal definition of the problem of interactive mapping specification, as well as a formal resolution process for which desirable properties are proved. Then, based on this formal resolution process, practical algorithms are provided. The approach behind these algorithms aims at reducing the number of boolean questions users have to answers by making use of quasi-lattice structures to order the set of possible mappings to explore, allowing an efficient pruning of the space of explored mappings. In order to improve this pruning, an extension of this approach to the use of integrity constraints is also provided. The second main contribution is a repairing process allowing to ensure that a mapping is “safe” with respect to a set of policy views defined on its source schema, i.e., that it does not leak sensitive information. A privacy-preservation protocol is provided to visualize the information leaks of a mapping, as well as a process to rewrite an input mapping into a safe one with respect to a set of policy views. As in the first contribution, this process comes with proofs of desirable properties. In order to reduce the number of interactions needed with the user, the interactive part of the repairing process is also enriched with the possibility of learning which rewriting is preferred by users, in order to obtain a completely automatic process. Last but not least, we present extensive experiments over the open source prototypes built from two contributions of this thesis