Cette thèse étudie les aspects mathématiques et numériques de phénomènes d'ondes dans les plasmas magnétiques. La réflectométrie, une technique de sonde des plasmas de fusion, est modélisée par les équations de Maxwell. Le tenseur de permittivité présente dans ce modèle des valeurs propres ainsi que des termes diagonaux qui s'annulent. La relation de dispersion met en évidence deux phénomènes cruciaux : coupures et résonances, lorsque le nombre d'onde s'annule ou tend vers l'infini. La partie I rassemble les résultats numériques. La grande nouveauté réside dans la définition d'une solution résonante. En effet, à cause des coefficients s'annulant continument en changeant de signe, la solution peut être singulière, i.e. avoir une composante non intégrable. Cependant, grâce au principe d'absorption limite, une solution résonante est explicitement définie comme la limite de solutions intégrables du problème régularisé. L'expression théorique de la singularité est validée par des tests numériques du passage à la limite. La partie II concerne l'approximation numérique. Elle comprend la mise en place d'une nouvelle méthode numérique adaptée aux coefficients réguliers. Celle-ci est basée sur la formulation variationnelle Ultra Faible mais nécessite des fonctions de base spécifiques, construites comme approximations locales du problème adjoint. L'analyse de convergence est effectuée en dimension un, en dimension deux la construction des fonctions de base et leur propriété d'interpolation sont détaillées. La méthode d'ordre élevé obtenue permet de simuler le phénomène de coupure tandis que simuler le phénomène de résonance en dimension deux reste un défi.

Ce travail a permis de développer une technique non invasive d’identification de la distribution du courant dans une pile à combustible à partir du champ magnétique externe. La mesure du champ s’effectue sur un ensemble de points de mesures choisis spécialement pour détecter les variations du champ par rapport à un fonctionnement optimal de la pile. Les deux composantes du champ magnétique sensibles aux variations sont utilisées. La mesure du champ exploitable est la différence entre un mode considéré sain et un mode quelconque de fonctionnement. Autour de ces mesures de champ magnétique, un problème inverse est modélisé en explorant plusieurs approches de paramétrisation de la distribution du courant. Le caractère mal posé du problème s’est traduit par la non-unicité de la solution et sa sensibilité au bruit. L’affranchissement de ces problèmes est atteint par la régularisation du modèle inverse. L’outil développé permis de reconstruire la distribution du courant indépendamment de la taille du défaut dans la limite de sensibilité. La validation est faite sur un simulateur électrique de pile à combustible et sur une pile de type GENEPAC dans un environnement de laboratoire. / A noninvasive technique for identifying the current distribution in a fuel cell from the external magnetic field is developed. The magnetic field measurements are carried out on a set of points chosen in order to detect only the variations of the magnetic field in a case of a faulty fuel cell. Two components of the magnetic field are used because they are sensitive to current heterogeneities. The exploitable measurement of the magnetic field is the difference between a healthy mode and a faulty mode.An ill-posed inverse problem is modeled by investigating several parametrization approaches of the current distribution. The ill-posed property of the problem conducts to a non-uniqueness of the solution and a high sensitivity to the noise. A regularization method is used in order to get a stable solution. The developed tool allows identifying the current distribution independently of the size of the fault within the sensitivity limit. A validation is done on a fuel cell simulator and on a GENEPAC fuel cell in a laboratory conditions.

L’objectif de cette thèse est d’étudier et de donner des outils pour la compréhension des problèmes de perturbations singulières pour des modèles épidémiques et des problèmes de dynamiques de populations. Les modèles considérés sont des équations structurées en âge qui peuvent dans certains cas se réécrire comme des équations à retard. L’étude de ces classes d’exemples s’est faite avec succès et a permis de comprendre et de mettre en évidence toute la complexité et l’étendue de ces problèmes. Comme on peut le remarquer dans la littérature, l’une des clés fondamentales à la compréhension de ces problèmes est l’étude des variétés normalement hyperboliques en dimension infinie que nous avons largement étudiées dans cette thèse. L’approche utilisée est la méthode de Lyapunov-Perron. Ce qui nous a amené à étudier les problèmes de persistance et d’existence de trichotomie (dichotomie) exponentielle qui sont des éléments fondamentaux dans l’utilisation de cette méthode. / In this thesis we aim to give tools to understand singular perturbations in epidemic model sand population dynamic models. We study some singularly perturbed delay differential equation which does not enter into the class frame work of geometric singular perturbation for delay differential equations. An example of singularly perturbed age structured model is also studied. The study of these examples allowed us to understand and highlight some complexities of these problems. One of the main tools in understanding such questions is the normally hyperbolic manifolds theory which is our central focus in this thesis. The approach used here is the Lyapunov-Perron method. Therefore the problems of persistence and existence of exponential trichotomy (dichotomy) are also stressed since there are one of the mainingredients of this method.

Les lipides des microalgues et les glucides de cyanobactéries peuvent être transformés en biodiesel et en bioéthanol, respectivement. L'amélioration de la production de ces molécules doit prendre en compte les entrées périodiques (principalement la lumière) forçant le réseau métabolique de ces organismes photosynthétiques. Il est donc nécessaire de tenir compte de la dynamique du réseau métabolique en réduisant sa dimension pour assurer la maniabilité mathématique. Le but de ce travail est de concevoir une approche originale pour réduire les réseaux métaboliques dynamiques tout en conservant la dynamique de base. Cette méthode est basée sur une séparation en échelles de temps. Pour une classe de modèles de réseaux métaboliques décrits par des ODE, la dynamique des systèmes réduits est calculée à l'aide du théorème de Tikhonov pour les systèmes singulièrement perturbés. Cette approximation quasi-stationnaire coïncide avec la dynamique du réseau d'origine, avec une erreur bornée. L'approche est d'abord développée pour les systèmes de réaction pouvant être linéarisés autour d'un point de travail et forcés par des entrées continues. Ensuite, une généralisation de cette méthode est donnée pour les réseaux à réactions rapides de cinétiques de Michaelis-Menten et tout type de cinétiques lentes, prenant également en compte un nombre fini d'entrées continues externes. La méthode de réduction met en évidence une relation entre la grandeur de la concentration des métabolites et la gamme des vitesses de réaction : les métabolites consommés par les réactions rapides ont une concentration inférieure d'un ordre de grandeur à celle des métabolites consommés à faible vitesse. Cette propriété est satisfaite pour les métabolites à dynamique rapide ne se trouvant pas dans un piège de flux, concept introduit dans ce travail. Le système réduit peut être calibré avec des données expérimentales à l'aide d'une procédure d'identification dédiée basée sur la minimisation. L'approche est illustrée par un réseau métabolique de microalgues autotrophes, comprenant le métabolisme central et représentant la dynamique des glucides et des lipides. Cette approche permet de bien ajuster les données expérimentales de Lacour et al. (2012) avec la microalgue Tisochrysis lutea. Enfin, un schéma visant à optimiser la production de molécules cibles est proposé en utilisant le système réduit. / Lipids from microalgae and carbohydrates from cyanobacteria can be transformed into biodiesel and bioethanol, respectively. Enhancing the production of these molecules must account for the periodic inputs (mainly light) forcing the metabolic network of these photosynthetic organisms. It is therefore necessary to account for the dynamics of the metabolic network, while reducing its dimension to ensure mathematical tractability. The aim of this work is to design an original approach to reduce dynamic metabolic networks while keeping the core dynamics. This method is based on time-scale separation. For a class of metabolic network models described by ODE, the dynamics of the reduced systems are computed using the theorem of Tikhonov for singularly perturbed systems. This Quasi Steady State Approximation accurately coincides with the original network dynamics, with a bounded error. The approach is first developed for reaction systems that can be linearized around a working point and that are forced by external continuous inputs. Then, a generalization of this method is given for networks with fast reactions of Michaelis-Menten kinetics and any type of slow kinetics, also considering a finite number of external continuous inputs. The reduction method highlights a relation between the concentration magnitude of the metabolites and the range of the reaction rates: the metabolites that are consumed by fast reactions have concentration one order of magnitude lower than metabolites consumed at slow rates. This property is satisfied for metabolites with fast dynamics that are not in a flux trap, a concept introduced in this work. The reduced system can be calibrated with experimental data using a dedicated identification procedure based on minimization. The approach is illustrated with an autotrophic microalgae metabolic network, including the core metabolism and representing the carbohydrates and lipids dynamics. The approach efficiently fits the experimental data from Lacour et al. (2012) with the microalgae Tisochrysis lutea. Finally, a scheme to optimize the production of target molecules is proposed using the reduced system.

Ce travail porte sur de nouvelles explosions de cycles (orbites périodiques), l'étude de leur structure par l'analyse qualitative, leur mise en évidence par simulation numérique (Auto, Xpp) et la discussion de leur pertinence dans des modèles mathématiques dans les neurosciences. De telles explosions se produisent dans les systèmes dynamiques lents-rapides. La plupart des neurones sont excitables, dès 1940, Hodgkin identifia trois classes fondamentales d'axones excitables distinguées par leurs réponses à un courant injecté d'amplitude variable. A l'aide de la fonction de Lambert, nous étudions la transition entre les types I et II par des explosions de cycle incomplètes, initiées par une bifurcation de Hopf singulière et qui se terminent dans une bifurcation homocline dans des systèmes une variable rapide/une variable lente. Vient ensuite une étude poussée du système de Hindmarsh-Rose. Il s'agit d'un système deux variables rapides/une variable lente qui produit des oscillations en salves (ou bursting). Nous généralisons la notion d'ensembles candidats-limites-périodiques (clp) aux systèmes tridimensionnels, il s'agit des ensembles invariants du système à la limite singulière. A l'aide de ces derniers, nous obtenons une description très fine de la déformation du cycle limite jusqu'à l'addition d'un nouveau spike au burst. Nous finissons par une étude de la minimalité du modèle de F. Clément et J.-P. Françoise. Ce dernier est un système 4D qui modélise l¿activité des neurones à GnRH. Nous étudions un système une variable rapide/deux variables lentes qui reproduit certaines des caractéristiques du modèle 4D, notamment des Mixed-Modes oscillations. / This thesis is focussed on the analysis of novel explosions of limit cycles (periodic orbits). We provide a study of their structure by qualitative analysis, exhibit evidences of their existence by numerical simulations (Auto, Xpp) and propose a discussion of their relevance in mathematical modeling for neurosciences. Such explosions occur in the slow-fast dynamical systems. Most of neurons are excitable, Hodgkin (1940) identified three fundamental classes of excitable axon distinguished by their responses to a current of variable amplitude injected. Using the Lambert function, we study the transition between types I and II by incomplete explosion of cycle. This explosion, produced by a planar vector field with one fast/one slow variable, is initiated by a singular Hopf bifurcation and ends via a homoclinic bifurcation. The next chapter proposed a study of the Hindmarsh-Rose system. This system, composed of one fast/ two slow variables, is well known to produce square wave bursting oscillation. We generalize the notion of candidate-limit-perodic sets (CLP-sets) to three-dimensional systems. A CLP-set is an invariant set of the system in the singular limit. Using these, we get a very acurate description of the limit cycle deformation under the variation of a parameter until the addition of a new spike to burst. Finally, we propose a study fot the minimality of the model introduced by F. Clement and J.-P. Françoise. The latter is a 4D system that models the activity of GnRH neurons. We study a system composed by one fast /two slow variables that reproduces some of the features of the 4D model, including Mixed-Modes oscillations.

Afin d'exploiter au mieux son parc nucléaire, la R&D d'EDF est en train de développer une nouvelle chaîne de calcul pour simuler le cœur des réacteurs nucléaires avec des outils à l'état de l'art. Ces calculs nécessitent une grande quantité de données physiques, en particulier les sections efficaces. Dans la simulation d'un cœur complet, le nombre de valeurs des sections efficaces est de l'ordre de plusieurs milliards. Ces sections efficaces peuvent être représentées comme des fonctions multivariées dépendant de plusieurs paramètres physiques. La détermination des sections efficaces étant un calcul complexe et long, nous pouvons donc les précalculer en certaines valeurs des paramètres (caluls hors ligne) puis les évaluer en tous points par une interpolation (calculs en ligne). Ce processus demande un modèle de reconstruction des sections efficaces entre les deux étapes. Pour réaliser une simulation plus fidèle du cœur dans la nouvelle chaîne d'EDF, les sections efficaces nécessitent d'être mieux représentées en prenant en compte de nouveaux paramètres. Par ailleurs, la nouvelle chaîne se doit d'être en mesure de calculer le réacteur dans des situations plus larges qu'actuellement. Le modèle d'interpolation multilinéaire pour reconstruire les sections efficaces est celui actuellement utilisé pour répondre à ces objectifs. Néanmoins, avec ce modèle, le nombre de points de discrétisation augmente exponentiellement en fonction du nombre de paramètres ou de manière considérable quand on ajoute des points sur un des axes. Par conséquence, le nombre et le temps des calculs hors ligne ainsi que la taille du stockage des données deviennent problématique. L'objectif de cette thèse est donc de trouver un nouveau modèle pour répondre aux demandes suivantes : (i)-(hors ligne) réduire le nombre de précalculs, (ii)-(hors ligne) réduire le stockage de données pour la reconstruction et (iii)-(en ligne) tout en conservant (ou améliorant) la précision obtenue par l'interpolation multilinéaire. D'un point de vue mathématique, ce problème consiste à approcher des fonctions multivariées à partir de leurs valeurs précalculées. Nous nous sommes basés sur le format de Tucker - une approximation de tenseurs de faible rang afin de proposer un nouveau modèle appelé la décomposition de Tucker . Avec ce modèle, une fonction multivariée est approchée par une combinaison linéaire de produits tensoriels de fonctions d'une variable. Ces fonctions d'une variable sont construites grâce à une technique dite de décomposition en valeurs singulières d'ordre supérieur (une « matricization » combinée à une extension de la décomposition de Karhunen-Loève). L'algorithme dit glouton est utilisé pour constituer les points liés à la résolution des coefficients dans la combinaison de la décomposition de Tucker. Les résultats obtenus montrent que notre modèle satisfait les critères exigés sur la réduction de données ainsi que sur la précision. Avec ce modèle, nous pouvons aussi éliminer a posteriori et à priori les coefficients dans la décomposition de Tucker. Cela nous permet de réduire encore le stockage de données dans les étapes hors ligne sans réduire significativement la précision. / In order to optimize the operation of its nuclear power plants, the EDF's R&D department iscurrently developing a new calculation chain to simulate the nuclear reactors core with state of the art tools. These calculations require a large amount of physical data, especially the cross-sections. In the full core simulation, the number of cross-section values is of the order of several billions. These cross-sections can be represented as multivariate functions depending on several physical parameters. The determination of cross-sections is a long and complex calculation, we can therefore pre-compute them in some values of parameters (online calculations), then evaluate them at all desired points by an interpolation (online calculations). This process requires a model of cross-section reconstruction between the two steps. In order to perform a more faithful core simulation in the new EDF's chain, the cross-sections need to be better represented by taking into account new parameters. Moreover, the new chain must be able to calculate the reactor in more extensive situations than the current one. The multilinear interpolation is currently used to reconstruct cross-sections and to meet these goals. However, with this model, the number of points in its discretization increases exponentially as a function of the number of parameters, or significantly when adding points to one of the axes. Consequently, the number and time of online calculations as well as the storage size for this data become problematic. The goal of this thesis is therefore to find a new model in order to respond to the following requirements: (i)-(online) reduce the number of pre-calculations, (ii)-(online) reduce stored data size for the reconstruction and (iii)-(online) maintain (or improve) the accuracy obtained by multilinear interpolation. From a mathematical point of view, this problem involves approaching multivariate functions from their pre-calculated values. We based our research on the Tucker format - a low-rank tensor approximation in order to propose a new model called the Tucker decomposition . With this model, a multivariate function is approximated by a linear combination of tensor products of one-variate functions. These one-variate functions are constructed by a technique called higher-order singular values decomposition (a « matricization » combined with an extension of the Karhunen-Loeve decomposition). The so-called greedy algorithm is used to constitute the points related to the resolution of the coefficients in the combination of the Tucker decomposition. The results obtained show that our model satisfies the criteria required for the reduction of the data as well as the accuracy. With this model, we can eliminate a posteriori and a priori the coefficients in the Tucker decomposition in order to further reduce the data storage in online steps but without reducing significantly the accuracy.

Le domaine de l’écologie des réseaux est encore limité dans sa capacité à faire des inférences mondiales à grande échelle. Ce défi est principalement dû à la difficulté d’échantillonnage des interactions sur le terrain, entraînant de nombreuses « lacunes » en ce qui concerne la couverture mondiale des données. Cette thèse adopte une approche « centrée sur les méthodes » de l’écologie des réseaux et se concentre sur l’idée de développer des outils pour aider à combler les lacunes en matière de données en présentant la prédiction comme une alternative accessible à l’échantillonnage sur le terrain et introduit deux « outils » différents qui sont prêts à poser des questions à l’échelle mondiale. Le chapitre 1 présente les outils que nous pouvons utiliser pour faire des prédictions de réseaux et est motivé par l’idée selon laquelle avoir la capacité de prédire les interactions entre les espèces grâce à l’utilisation d’outils de modélisation est impératif pour une compréhension plus globale des réseaux écologiques. Ce chapitre comprend une preuve de concept (dans laquelle nous montrons comment un simple modèle de réseau neuronal est capable de faire des prédictions précises sur les interactions entre espèces), une évaluation des défis et des opportunités associés à l’amélioration des prédictions d’interaction et une feuille de route conceptuelle concernant l’utilisation de modèles prédictifs pour les réseaux écologiques. Les chapitres 2 et 3 sont étroitement liés et se concentrent sur l’utilisation de l’intégration de graphiques pour la prédiction de réseau. Essentiellement, l’intégration de graphes nous permet de transformer un graphe (réseau) en un ensemble de vecteurs, qui capturent une propriété écologique du réseau et nous fournissent une abstraction simple mais puissante d’un réseau d’interaction et servent de moyen de maximiser les informations disponibles. dispo- nibles à partir des réseaux d’interactions d’espèces. Parce que l’intégration de graphes nous permet de « décoder » les informations au sein d’un réseau, elle est conçue comme un outil de prédiction de réseau, en particulier lorsqu’elle est utilisée dans un cadre d’apprentissage par transfert. Elle s’appuie sur l’idée que nous pouvons utiliser les connaissances acquises en résolvant un problème connu. et l’utiliser pour résoudre un problème étroitement lié. Ici, nous avons utilisé le métaweb européen (connu) pour prédire un métaweb pour les espèces canadiennes en fonction de leur parenté phylogénétique. Ce qui rend ce travail particulière- ment passionnant est que malgré le faible nombre d’espèces partagées entre ces deux régions, nous sommes capables de récupérer la plupart (91%) des interactions. Le chapitre 4 approfondit la réflexion sur la complexité des réseaux et les différentes ma- nières que nous pourrions choisir de définir la complexité. Plus spécifiquement, nous remet- tons en question les mesures structurelles plus traditionnelles de la complexité en présentant l’entropie SVD comme une mesure alternative de la complexité. Adopter une approche phy- sique pour définir la complexité nous permet de réfléchir aux informations contenues dans un réseau plutôt qu’à leurs propriétés émergentes. Il est intéressant de noter que l’entropie SVD révèle que les réseaux bipartites sont très complexes et ne sont pas nécessairement conformes à l’idée selon laquelle la complexité engendre la stabilité. Enfin, je présente le package Julia SpatialBoundaries.jl. Ce package permet à l’utili- sateur d’implémenter l’algorithme de wombling spatial pour des données disposées de manière uniforme ou aléatoire dans l’espace. Étant donné que l’algorithme de wombling spatial se concentre à la fois sur le gradient et sur la direction du changement pour un paysage donné, il peut être utilisé à la fois pour détecter les limites au sens traditionnel du terme ainsi que pour examiner de manière plus nuancée la direction des changements. Cette approche pourrait être un moyen bénéfique de réfléchir aux questions liées à la détection des limites des réseaux et à leur relation avec les limites environnementales. / The field of network ecology is still limited in its ability to make large-scale, global inferences. This challenge is primarily driven by the difficulty of sampling interactions in the field, leading to many ‘gaps’ with regards to global coverage of data. This thesis takes a ’methods-centric’ approach to network ecology and focuses on the idea of developing tools to help with filling in the the data gaps by presenting prediction as an accessible alternative to sampling in the field and introduces two different ’tools’ that are primed for asking questions at global scales. Chapter 1 maps out tools we can use to make network predictions and is driven by the idea that having the ability to predict interactions between species through the use of modelling tools is imperative for a more global understanding of ecological networks. This chapter includes a proof-of-concept (where we show how a simple neural network model is able to make accurate predictions about species interactions), an assessment of the challenges and opportunities associated with improving interaction predictions, and providing a conceptual roadmap concerned with the use of predictive models for ecological networks. Chapters 2 and 3 are closely intertwined and are focused on the use of graph embedding for network prediction. Essentially graph embedding allows us to transform a graph (net- work) into a set of vectors, which capture an ecological property of the network and provides us with a simple, yet powerful abstraction of an interaction network and serves as a way to maximise the available information available from species interaction networks. Because graph embedding allows us to ’decode’ the information within a network it is primed as a tool for network prediction, specifically when used in a transfer learning framework, this builds on the idea that we can take the knowledge gained from solving a known problem and using it to solve a closely related problem. Here we used the (known) European metaweb to predict a metaweb for Canadian species based on their phylogenetic relatedness. What makes this work particularly exciting is that despite the low number of species shared between these two regions we are able to recover most (91%) of interactions. Chapter 4 delves into thinking about the complexity of networks and the different ways we might choose to define complexity. More specifically we challenge the more traditional structural measures of complexity by presenting SVD entropy as an alternative measure of complexity. Taking a physical approach to defining complexity allows us to think about the information contained within a network as opposed to their emerging properties. Interest- ingly, SVD entropy reveals that bipartite networks are highly complex and do not necessarily conform to the idea that complexity begets stability. Finally, I present the Julia package SpatialBoundaries.jl. This package allows the user to implement the spatial wombling algorithm for data arranged uniformly or randomly across space. Because the spatial wombling algorithm focuses on both the gradient as well as the direction of change for the given landscape it can be used both for detecting boundaries in the traditional sense as well as a more nuanced look at at the direction of changes. This approach could be a beneficial way with which to think about questions which relate to boundary detection for networks and how these relate to environmental boundaries.

Une nouvelle structure de dictionnaire adaptés aux décompositions itératives de type poursuite, appelée un Iteration-Tuned Dictionary (ITD), est présentée. Les ITDs sont structurés en couche, chaque couche se composant d'un ensemble de dictionnaires candidats. Les décompositions itératives basées ITD sont alors réalisées en sélectionnant, à chaque itération i, l'un des dictionnaires de la i-ième couche. Une structure générale des ITDs est proposée, ainsi qu'une variante structurée en arbre appelée Tree-Structured Iteration-Tuned Dictionary (TSITD) et une version contrainte de cette dernière, appelée Iteration-Tuned and Aligned Dictionary (ITAD). Ces structures sont comparées à plusieurs méthodes de l'état de l'art et évaluées dans des applications de débruitage et de compression d'images. Un codec basé sur le schéma ITAD est également présenté et comparé à JPEG2000 dans des évaluations qualitatives et quantitatives. Dans le contexte de l'indexation d'images, un nouveau système de recherche approximative des plus proches voisins est également introduit, qui utilise les représentations parcimonieuses pour réduire la complexité de la recherche. La méthode traite l'instabilité dans la sélection des atomes lorsque l'image est soumise à de faibles transformations affines. Un nouveau système de conditionnement des données est également introduit, permettant de mieux distribuer les données sur la sphère unitaire tout en préservant leurs distances angulaires relatives. Il est montré que cette méthode améliore le compromis complexité/performance de la recherche approximative basée décompositions parcimonieuses.

Les questions de sûreté et d'incertitudes sont au centre des études de faisabilité pour un site de stockage souterrain de déchets nucléaires, en particulier l'évaluation des incertitudes sur les indicateurs de sûreté qui sont dues aux incertitudes sur les propriétés du sous-sol et des contaminants. L'approche globale par les méthodes probabilistes de type Monte Carlo fournit de bons résultats, mais elle demande un grand nombre de simulations. La méthode déterministe étudiée ici est complémentaire. Reposant sur la décomposition en valeurs singulières de la dérivée du modèle, elle ne donne qu'une information locale, mais elle est beaucoup moins coûteuse en temps de calcul. Le modèle d'écoulement suit la loi de Darcy et le transport des radionucléides autour du site de stockage est modélisé par une équation de diffusion-convection linéaire. Différentiation à la main et différentiation automatique sont comparées sur ces modèles en mode direct et en mode adjoint. Une étude comparée des deux approches probabiliste et déterministe pour l'analyse de la sensibilité des flux de contaminants aux exutoires par rapport aux variations des paramètres d'entrée est menée sur des données réalistes fournies par l'ANDRA. Des outils génériques d'analyse de sensibilité et de couplage de code sont développés en langage Caml. Ils permettent à l'utilisateur de ces plates-formes génériques de ne fournir que la partie spécifique de l'application dans le langage de son choix. Une étude sur les écoulements diphasiques eau/air partiellement saturés en hydrogéologie porte sur les limitations des approximations de Richards et de la formulation en pression globale issue du domaine pétrolier.

La quantification par déformation et la correspondance de McKay forment les grands thèmes de l'étude qui porte sur des variétés algébriques singulières, des quotients d'algèbres de polynômes et des algèbres de polynômes invariants sous l'action d'un groupe fini. Nos principaux outils sont les cohomologies de Poisson et de Hochschild et la théorie des représentations. Certains calculs formels sont effectués avec Maple et GAP. Nous calculons les espaces d'homologie et de cohomologie de Hochschild des surfaces de Klein, en développant une généralisation du Théorème de HKR au cas de variétés non lisses et utilisons la division multivariée et les bases de Gröbner. La clôture de l'orbite nilpotente minimale d'une algèbre de Lie simple est une variété algébrique singulière sur laquelle nous construisons des star-produits invariants, grâce à la décomposition BGS de l'homologie et de la cohomologie de Hochschild, et à des résultats sur les invariants des groupes classiques. Nous explicitons les générateurs de l'idéal de Joseph associé à cette orbite et calculons les caractères infinitésimaux. Pour les algèbres de Lie simples B, C, D, nous établissons des résultats généraux sur l'espace d'homologie de Poisson en degré 0 de l'algèbre des invariants, qui vont dans le sens de la conjecture d'Alev et traitons les rangs 2 et 3. Nous calculons des séries de Poincaré à 2 variables pour des sous-groupes finis du groupe spécial linéaire en dimension 3, montrons que ce sont des fractions rationnelles, et associons aux sous-groupes une matrice de Cartan généralisée pour obtenir une correspondance de McKay algébrique en dimension 3. Toute l'étude a donné lieu à 4 articles / Deformation quantization and McKay correspondence form the main themes of the study which deals with singular algebraic varieties, quotients of polynomial algebras, and polynomial algebras invariant under the action of a finite group. Our main tools are Poisson and Hochschild cohomologies and representation theory. Certain calculations are made with Maple and GAP. We calculate Hochschild homology and cohomology spaces of Klein surfaces by developing a generalization of HKR theorem in the case of non-smooth varieties and use the multivariate division and the Groebner bases. The closure of the minimal nilpotent orbit of a simple Lie algebra is a singular algebraic variety : on this one we construct invariant star-products, with the help of the BGS decomposition of Hochschild homology and cohomology, and of results on the invariants of the classical groups. We give the generators of the Joseph ideal associated to this orbit and calculate the infinitesimal characters. For simple Lie algebras of type B, C, D, we establish general results on the Poisson homology space in degree 0 of the invariant algebra, which support Alev's conjecture, then we are interested in the ranks 2 and 3. We compute Poincaré series of 2 variables for the finite subgroups of the special linear group in dimension 3, show that they are rational fractions, and associate to the subgroups a generalized Cartan matrix in order to obtain a McKay correspondence in dimension 3. All the study comes from 4 papers