Un des ingrédients principaux pour la simulation numérique des écoulements des fluides (hydrocarbures, gaz naturel) en milieux poreux est la discrétisation des termes elliptiques anisotropes et hétérogènes. Dans l’industrie pétrolière, la nécessité d’améliorer la précision des simulations dans les régions proches des puits a suscité l’utilisation de maillages non structurés généraux et des tenseurs de perméabilité pleins. Notre défi a été de trouver des discrétisations consistantes et robustes des termes elliptiques anisotropes, hétérogènes sur maillages généraux. Notre recherche s’est focalisée sur des méthodes volumes finis qui soient consistantes, stables sur maillages polyédriques généraux, robustes par rapport à l’anisotropie et l’hétérogénéité du tenseur de perméabilité, qui mènent à des systèmes linéaires bien conditionnés pour lesquels des stratégies de préconditionnement optimales peuvent être conçues et qui donnent un stencil compact pour réduire les communications dans les implémentations parallèles. Pour répondre à cette recherche, nous avons proposé plusieurs schémas tels que MPFA O généralisé, G scheme, CG method, VFSYM, DIOPTRE. Nous avons prouvé également que toutes ces méthodes convergent sous des hypothèses adéquates à la fois sur le tenseur de perméabilité et le maillage / One of the key ingredients for the numerical simulation of Darcy flow in heterogeneous porous media is the discretization of anisotropic heterogeneous elliptic terms. In the oil industry, the need to improve accuracy in near wellbore regions has prompted the use of general unstructured meshes and full permeability tensors. Our effort has therefore been devoted to find consistent and robust finite volume discretizations of anisotropic, heterogeneous elliptic terms on general meshes. Our research was focused on finite volume methods which are consistent and coercive on general polyhedral meshes as well as robust with respect to the anisotropy and heterogeneity of the permeability tensor ; yield well-conditioned linear systems for which optimal preconditioning strategies can be devised ; have a narrow stencil to reduce the communications in parallel implementations. To answer to this search, we have proposed several scemes such that generalized MPFA O, G scheme, CG method, VFSYM, DIOPTRE. We proved also the convergence of all these methods under suitable assumptions on both the permeability tensor and the mesh

En permettant le partage à grande échelle de données et d'algorithmes et en fournissant une quantité importante de puissance de calcul et de stockage, les grilles de calcul sont des plates-formes intéressantes pour les applications d'analyse d'images médicales. Dans cette thèse, nous étudions un problème d'analyse d'images médicales qui s'avère être une application dimensionnante pour les grilles, conduisant au développement de nouvelles méthodes et outils pour la description, l'implémentation et l'optimisation de flots de traitements. Le problème applicatif étudié est l'évaluation de la précision d'algorithmes de recalage d'images médicales en l'absence de vérité terrain. Nous faisons passer à l'échelle une méthode statistique d'évaluation de ces algorithmes et nous montrons des résultats de précision sur une base de données liée au suivi de la radiothérapie du cerveau. Ces résultats permettent notamment de détecter des défauts très légers au sein des données. Nous étendons ce schéma pour quantifier l'impact de la compression des images sur la qualité du recalage. Cette application étant représentative de problèmes typiques survenant sur les grilles, nous nous attachons à son déploiement et à son exécution sur ce type d'infrastructures. Pour faciliter une parallélisation transparente, nous adoptons un modèle générique de flots de traitements, dont nous proposons une nouvelle taxonomie. Pour répondre aux limitations de performance des moteurs d'exécution de flots existants, nous présentons MOTEUR, qui permet d'exploiter les différents types de parallélisme inhérents à ces applications. La définition d'un nouvel opérateur de composition de données facilite la description des applications d'analyse d'images médicales sur les grilles. Par une comparaison entre la grille de production EGEE et des grappes dédiées de Grid'5000, nous mettons en évidence l'importance de la variabilité de la latence sur une grille de production. En conséquence, nous proposons un modèle probabiliste du temps d'exécution d'un flot de traitement sur une grille. Ce modèle est centré sur l'utilisateur : il considère la grille toute entière comme une boîte noire introduisant une latence aléatoire sur le temps d'exécution d'une tâche. A partir de ce modèle, nous proposons trois stratégies d'optimisation visant à réduire l'impact de la latence et de sa variabilité : dans un flot de traitement, grouper les tâches séquentiellement liées permet de réduire la latence moyenne rencontrée, optimiser la valeur du délai d'expiration des tâches prémunit contre les valeurs extrêmes de la latence et optimiser la granularité des tâches permet de réduire le risque de rencontrer de fortes latences. Des accélérations significatives sont ainsi obtenues.

Cette thèse est organisé en deux parties. Dans la première partie nous nous adressons à un problème vieux dans la théorie de jauge - le calcul du spectre et des fonctions d'onde. La stratégie que nous proposons est de construire une base d'états stochastiques de liens de Bargmann, construite à partir d'une distribution physique de densité de probabilité. Par la suite, nous calculons les amplitudes de transition entre ces états par une approche analytique, en utilisant des intégrales de chemin standards ainsi que la théorie des groupes. Également, nous calculons numériquement matrices symétrique et hermitienne des amplitudes de transition, via une méthode Monte Carlo avec échantillonnage pondéré. De chaque matrice, nous trouvons les valeurs propres et les vecteurs propres. En appliquant cette méthode â la théorie de jauge U(l) en deux dimensions spatiales, nous essayons d'extraire et de présenter le spectre et les fonctions d'onde de cette théorie pour des grilles de petite taille. En outre, nous essayons de faire quelques ajustement dynamique des fenêtres de spectres d'énergie et les fonctions d'onde. Ces fenêtres sont outiles de vérifier visuellement la validité de l'hamiltonien Monte Carlo, et de calculer observables physiques. Dans la deuxième partie nous étudions le comportement chaotique de deux systèmes de billard classiques, par la théorie des matrices aléatoires. Nous considérons un gaz périodique de Lorentz à deux dimensions dans des régimes de horizon fini et horizon infini. Nous construisons quelques matrices de longueurs de trajectoires de un particule mobile dans ce système, et réalisons des études des spectres de ces matrices par l'analyse numérique. Par le calcul numérique des distributions d'espacement de niveaux et rigidité spectral, nous constatons la statistique des espacements de niveaux suggère un comportement universel. Nous étudions également un tel comportement pour un système optique chaotique. En tant que quasi-système de potentiel, ses fluctuations dans l'espacement de ses niveaux suivent aussi un comportement GOE, ce qui est une signature d'universalité. Dans cette partie nous étudions également les propriétés de diffusion du gaz de Lorentz, par la longueur des trajectoires. En calculant la variance de ce quantité, nous montrons que dans le cas d'horizons finis, la variance de longueurs est linéaire par rapport au nombre de collisions de la particule dans le billard. Cette linéarité permet de définir un coefficient de diffusion pour le gaz de Lorentz, et dans un schéma général, elle est compatible avec les résultats obtenus par d'autres méthodes.

Des patchs en matériaux composites sont classiquement utilisés dans l'aéronautique pour réparer ou renforcer des structures métalliques. Ces deux types de matériaux présentent cependant des coefficients de dilatation très différents. De ce fait, leur assemblage, soumis à de grandes variations de température en cours de fonctionnement, est le siège de contraintes pouvant provoquer la rupture ou le décollement du patch. Afin d'évaluer précisement les distributions de contraintes dans l'assemblage collé, différents modèles de calcul ont été proposés : d'abord un modèle analytique bidimensionnel mettant en évidence des effets liés à la géométrie plane de la structure, puis un modèle thermo-viscoélastique intégrant le comportement non-linéaire de la colle en fonction du temps et de la température. Ces modèles ont été ensuite partiellement validés expérimentalement sur des éprouvettes patchées en utilisant des méthodes de mesure de champs.