La méthode des éléments frontières pour l’équation des ondes (BEM) est utilisée en acoustique eten électromagnétisme pour simuler la propagation d’une onde avec une discrétisation en temps(TD). Elle permet d’obtenir un résultat pour plusieurs fréquences à partir d’une seule résolution.Dans cette thèse, nous nous intéressons à l’implémentation efficace d’un simulateur TD-BEM sousdifférents angles. Nous décrivons le contexte de notre étude et la formulation utilisée qui s’exprimesous la forme d’un système linéaire composé de plusieurs matrices d’interactions/convolutions.Ce système est naturellement calculé en utilisant l’opérateur matrice/vecteur creux (SpMV). Nousavons travaillé sur la limite du SpMV en étudiant la permutation des matrices et le comportementde notre implémentation aidé par la vectorisation sur CPU et avec une approche par bloc surGPU. Nous montrons que cet opérateur n’est pas approprié pour notre problème et nous proposonsde changer l’ordre de calcul afin d’obtenir une matrice avec une structure particulière.Cette nouvelle structure est appelée une matrice tranche et se calcule à l’aide d’un opérateur spécifique.Nous décrivons des implémentations optimisées sur architectures modernes du calculhaute-performance. Le simulateur résultant est parallélisé avec une approche hybride (mémoirespartagées/distribuées) sur des noeuds hétérogènes, et se base sur une nouvelle heuristique pouréquilibrer le travail entre les processeurs. Cette approche matricielle a une complexité quadratiquesi bien que nous avons étudié son accélération par la méthode des multipoles rapides (FMM). Nousavons tout d’abord travaillé sur la parallélisation de l’algorithme de la FMM en utilisant différentsparadigmes et nous montrons comment les moteurs d’exécution sont adaptés pour relâcher le potentielde la FMM. Enfin, nous présentons des résultats préliminaires d’un simulateur TD-BEMaccéléré par FMM . / The time-domain BEM for the wave equation in acoustics and electromagnetism is used to simulatethe propagation of a wave with a discretization in time. It allows to obtain several frequencydomainresults with one solve. In this thesis, we investigate the implementation of an efficientTD-BEM solver using different approaches. We describe the context of our study and the TD-BEMformulation expressed as a sparse linear system composed of multiple interaction/convolutionmatrices. This system is naturally computed using the sparse matrix-vector product (SpMV). Wework on the limits of the SpMV kernel by looking at the matrix reordering and the behavior of ourSpMV kernels using vectorization (SIMD) on CPUs and an advanced blocking-layout on NvidiaGPUs. We show that this operator is not appropriate for our problem, and we then propose toreorder the original computation to get a special matrix structure. This new structure is called aslice matrix and is computed with a custom matrix/vector product operator. We present an optimizedimplementation of this operator on CPUs and Nvidia GPUs for which we describe advancedblocking schemes. The resulting solver is parallelized with a hybrid strategy above heterogeneousnodes and relies on a new heuristic to balance the work among the processing units. Due tothe quadratic complexity of this matrix approach, we study the use of the fast multipole method(FMM) for our time-domain BEM solver. We investigate the parallelization of the general FMMalgorithm using several paradigms in both shared and distributed memory, and we explain howmodern runtime systems are well-suited to express the FMM computation. Finally, we investigatethe implementation and the parametrization of an FMM kernel specific to our TD-BEM, and weprovide preliminary results.

Dans le contexte de cette thèse, nous nous focalisons sur des algorithmes pour l’algèbre linéaire numérique, plus précisément sur la résolution de grands systèmes linéaires creux. Nous mettons au point des méthodes de parallélisation pour le solveur linéaire hybride MaPHyS. Premièrement nous considerons l'aproche MPI+threads. Dans MaPHyS, le premier niveau de parallélisme consiste au traitement indépendant des sous-domaines. Le second niveau est exploité grâce à l’utilisation de noyaux multithreadés denses et creux au sein des sous-domaines. Une telle implémentation correspond bien à la structure hiérarchique des supercalculateurs modernes et permet un compromis entre les performances numériques et parallèles du solveur. Nous démontrons la flexibilité de notre implémentation parallèle sur un ensemble de cas tests. Deuxièmement nous considérons un approche plus innovante, où les algorithmes sont décrits comme des ensembles de tâches avec des inter-dépendances, i.e., un graphe de tâches orienté sans cycle (DAG). Nous illustrons d’abord comment une première parallélisation à base de tâches peut être obtenue en composant des librairies à base de tâches au sein des processus MPI illustrer par un prototype d’implémentation préliminaire de notre solveur hybride. Nous montrons ensuite comment une approche à base de tâches abstrayant entièrement le matériel peut exploiter avec succès une large gamme d’architectures matérielles. À cet effet, nous avons implanté une version à base de tâches de l’algorithme du Gradient Conjugué et nous montrons que l’approche proposée permet d’atteindre une très haute performance sur des architectures multi-GPU, multicoeur ainsi qu’hétérogène. / In the context of this thesis, our focus is on numerical linear algebra, more precisely on solution of large sparse systems of linear equations. We focus on designing efficient parallel implementations of MaPHyS, an hybrid linear solver based on domain decomposition techniques. First we investigate the MPI+threads approach. In MaPHyS, the first level of parallelism arises from the independent treatment of the various subdomains. The second level is exploited thanks to the use of multi-threaded dense and sparse linear algebra kernels involved at the subdomain level. Such an hybrid implementation of an hybrid linear solver suitably matches the hierarchical structure of modern supercomputers and enables a trade-off between the numerical and parallel performances of the solver. We demonstrate the flexibility of our parallel implementation on a set of test examples. Secondly, we follow a more disruptive approach where the algorithms are described as sets of tasks with data inter-dependencies that leads to a directed acyclic graph (DAG) representation. The tasks are handled by a runtime system. We illustrate how a first task-based parallel implementation can be obtained by composing task-based parallel libraries within MPI processes throught a preliminary prototype implementation of our hybrid solver. We then show how a task-based approach fully abstracting the hardware architecture can successfully exploit a wide range of modern hardware architectures. We implemented a full task-based Conjugate Gradient algorithm and showed that the proposed approach leads to very high performance on multi-GPU, multicore and heterogeneous architectures.

Cette thèse porte sur l'étude numérique des écoulements fluide-particules rencontrés dans l'industrie. Ces travaux se situent dans le cadre de la compréhension des phénomènes qui se déroulent dans des tambours tournants et réacteurs à lit fixe en présence de particules de forme non convexe. En effet, la forme des particules influence de manière importante la dynamique de ces milieux. A cet effet, nous nous sommes servis de la plateforme numérique parallèle Grans3D pour la dynamique des milieux granulaires et PeliGRIFF pour les écoulements multiphasiques. Dans la première partie de cette thèse, nous avons développé une nouvelle stratégie numérique qui permet de prendre en compte des particules de forme arbitrairement non convexe dans le solveur Grains3D. Elle consiste à décomposer une forme non convexe en plusieurs formes convexes quelconques. Nous avons nommé cette méthode « glued-convex ». Le modèle a été validé avec succès sur des résultats théoriques et expérimentaux de tambours tournants en présence de particules en forme de croix. Nous avons aussi utilisé le modèle pour simuler le chargement de réacteurs à lits fixes puis des lois de corrélation sur les taux de vide ont été déduites de nos résultats numériques. Dans ces travaux, nous avons aussi testé les performances parallèles de nos outils sur des simulations numériques à grande échelle de divers systèmes de particules convexes. La deuxième partie de cette thèse a été consacrée à l'extension du solveur PeliGRIFF à pouvoir prendre en compte la présence de particules multilobées (non convexes) dans des écoulements monophasiques. Une approche du type Simulation Numérique Directe, basée sur les Multiplicateurs de Lagrange Distribués / Domaine Fictif (DLM/FD), a alors été adoptée pour résoudre l'écoulement autour des particules. Une série d'études de convergence spatiale a été faite basée sur diverses configurations et divers régimes. Enfin, ces outils ont été utilisés pour simuler des écoulements au travers de lits fixes de particules de forme multi-lobée dans le but d'étudier l'influence de la forme des particules sur l'hydrodynamique dans ces lits. Les résultats ont montré une consistance avec les résultats expérimentaux disponibles dans la littérature. / Non convex granular media are involved in many industrial processes as, e.g., particle calcination/drying in rotating drums or solid catalyst particles in chemical reactors. In the case of optimizing the shape of catalysts, the experimental discrimination of new shapes based on packing density and pressure drop proved to be difficult due to the limited control of size distribution and loading procedure. There is therefore a strong interest in developing numerical tools to predict the dynamics of granular media made of particles of arbitrary shape and to simulate the flow of a fluid (either liquid or gas) around these particles. Non-convex particles are even more challenging than convex particles due to the potential multiplicity of contact points between two solid bodies. In this work, we implement new numerical strategies in our home made high-fidelity parallel numerical tools: Grains3D for granular dynamics of solid particles and PeliGRIFF for reactive fluid/solid flows. The first part of this work consists in extending the modelling capabilities of Grains3D from convex to non-convex particles based on the decomposition of a non-convex shape into a set of convex particles. We validate our numerical model with existing analytical solutions and experimental data on a rotating drum filled with 2D cross particle shapes. We also use Grains3D to study the loading of semi-periodic small size reactors with trilobic and quadralobic particles. The second part of this work consists in extending the modelling capabilities of PeliGRIFF to handle poly-lobed (and hence non-convex) particles. Our Particle Resolved Simulation (PRS) method is based on a Distributed Lagrange Multiplier / Fictitious Domain (DLM/FD) formulation combined with a Finite Volume / Staggered Grid (FV/SG) discretization scheme. Due to the lack of analytical solutions and experimental data, we assess the accuracy of our PRS method by examining the space convergence of the computed solution in assorted flow configurations such as the flow through a periodic array of poly-lobed particles and the flow in a small size packed bed reactor. Our simulation results are overall consistent with previous experimental work.

L'émergence de l'informatique dans les nuages met en avant de nombreux défis qui pourraient limiter l'adoption du paradigme Cloud. Tandis que la taille des données traitées par les applications Cloud augmente exponentiellement, un défi majeur porte sur la conception de solutions efficaces pour la gestion de données. Cette thèse a pour but de concevoir des mécanismes d'auto-adaptation pour des systèmes de gestion de données, afin qu'ils puissent répondre aux exigences des services de stockage Cloud en termes de passage à l'échelle, disponibilité et sécurité des données. De plus, nous nous proposons de concevoir un service de données qui soit à la fois compatible avec les interfaces Cloud standard dans et capable d'offrir un stockage de données à haut débit. Pour relever ces défis, nous avons proposé des mécanismes génériques pour l'auto-connaissance, l'auto-protection et l'auto-configuration des systèmes de gestion de données. Ensuite, nous les avons validés en les intégrant dans le logiciel BlobSeer, un système de stockage qui optimise les accès hautement concurrents aux données. Finalement, nous avons conçu et implémenté un système de fichiers s'appuyant sur BlobSeer, afin d'optimiser ce dernier pour servir efficacement comme support de stockage pour les services Cloud. Puis, nous l'avons intégré dans un environnement Cloud réel, la plate-forme Nimbus. Les avantages et les désavantages de l'utilisation du stockage dans le Cloud pour des applications réelles sont soulignés lors des évaluations effectuées sur Grid'5000. Elles incluent des applications à accès intensif aux données, comme MapReduce, et des applications fortement couplées, comme les simulations atmosphériques.

Dans cette thèse, nous nous intéressons à la modélisation et la simulation de dispositifs nanoélectroniques innovants. Premièrement, nous dérivons formellement un modèle avec masse effective pour décrire le transport quantique des électrons dans des nanostructures très fortement confinées. Des simulations numériques illustrent l'intérêt du modèle obtenu pour un dispositif simplifié mais déjà significatif. La deuxième partie est consacrée à l'étude du transport non ballistique dans ces mêmes structures confinées. Nous analysons rigoureusement un modèle de drift-diffusion et puis nous décrivons et implémentons une approche de couplage spatial classique-quantique. Enfin, nous modélisons et simulons un nanodispositif de spintronique. Plus précisement, nous étudions le renversement d'aimantation dans un matériau ferromagnétique multi-couches sous l'effet d'un courant de spin.

Les modèles de transport réactif sont un outil basique pour la modélisation de l'intéraction entre les réactions chimiques et l'écoulement du fluide dans un milieu poreux. Nous présentons un modèle de transport réactif multi-espèces totalement réduit incluant des réactions cinétiques et en équilibre. Une formulation structurée ainsi que différentes approches numériques sont proposées. Les méthodes de décomposition de domaine offrent la possibilité de diviser des problèmes de grande taille dans des problèmes plus petits dont la solution se fait en parallèle. Partant d'un point de vue géométrique, nous présentons la classe des méthodes de Schwarz ayant prouvé une haute performance dans de nombreuses applications. Des questions quant à la réalisation d'une décomposition de domaine et des conditions de transmission au niveau discret sont traitées dans le contexte des volumes finis. Nous proposons et validons numériquement un schéma de volumes finis hybrides pour l'opérateur d'advection-diffusion étant particulièrement adapté à l'utilisation dans le contexte d'une décomposition de domaine. Nous étudions théoriquement et numériquement des méthodes de Schwarz relaxation d'ondes en détail pour un système de deux espèces couplées de type transport réactif avec des termes de couplage linéaire et non-linéaire. Des résultats qualifiant le problème comme bien posé ainsi que la convergence des méthodes de décomposition de domaine sont développés et la sensibilité du comportement de convergence de l'algorithme de Schwarz par rapport au terme de couplage est étudiée. Finalement, nous appliquons une méthode de Schwarz relaxation d'ondes au modèle de transport réactif multi-espèces présenté.

Les enceintes de confinement des centrales nucléaires françaises sont réalisées en béton précontraint de 1,2 m d'épaisseur. L'étanchéité de la structure dépend principalement de la précontrainte résiduelle du béton. Or les dispositifs de surveillance des déformations différées font apparaître des différences sensibles, suivant les sites, que les modèles de calcul réglementaires ne prennent pas en compte d'une façon satisfaisante. Pour améliorer la gestion du parc de centrales, au sens durée de vie, EDF a lancé en 1991 un vaste programme d'étude centré sur le matériau béton, et visant à prévoir le comportement réel en fluage des enceintes déjà construites. Ici, nous analysons une par une les différentes déformations différées du béton. Une analyse fine des phénomènes physico chimiques qui en sont à l'origine, nous conduit à proposer une modélisation d'ingénieur, sur la base d'un modèle du type milieu continu équivalent. Enfin, nous présentons les résultats de nos simulations et nous concluons en donnant les enseignements théoriques et pratiques de cette étude ainsi que quelques propositions réglementaires.

On présente ici deux types de modélisations, ainsi que de nombreuses données expérimentales, portant sur les déformations des bétons et en particulier celles des bétons à hautes performances. La première a pour but de prévoir les amplitudes des déformations instantanées, de retrait endogène et de fluage propre, en fonction de la formulation. On utilise pour cela une approche par homogénéisation, en adaptant le modèle bisphère de Hashin au calcul du module du béton. La démarche de modélisation, fondée sur la prise en compte de la compacité réelle du squelette granulaire, aboutit à une géométrie trisphère. On fait le choix, après l'avoir justifié, de fixer les coefficients de Poisson à un cinquième. Puis, on applique le modèle trisphère au calcul des amplitudes de fluage propre et de retrait endogène. Les paramètres principaux de formulation (rapports eau sur ciment, fumée de silice sur ciment, proportion granulaire, classe du ciment) sont pris en compte. L'approche appliquée au fluage est confortée par un calcul viscoélastique linéaire. Dans la deuxième approche, on propose un modèle simplifié, de type réglementaire, pour le calcul des déformations des bétons à hautes performances en fonction du temps, et d'un nombre Limité de paramètres. Les apports principaux de la modélisation concernent la prise en compte de l'autodessiccation, de la maturité et de la résistance du matériau au moment du chargement. Le modèle a été validé d'après de nombreuses données, et pour des bétons de maturités très variables au chargement, variant de seize heures à plus de deux ans. Elaboré dans le cadre de l'Association Française pour le Recherche et l'Etude des Matériaux et des Structures (AFREM), le modèle est proposé comme annexe au règlement français de béton précontraint (BPEL), pour les bétons couvrant la gamme B40-B80.

Ce document présente les résultats obtenus pour différents bétons et pâtes de ciment, en ce qui concerne la microstructure, les propriétés hydriques et les déformations de retrait. En particulier, des matériaux à très hautes performances ont été étudiés (formulations avec un faible rapport eau sur ciment, contenant fumée de silice et fluidifiant). De nombreuses techniques expérimentales ont été utilisées à cet effet. Ainsi, les paramètres texturaux ont été déterminés par des méthodes classiques d'investigation de la structure poreuse (intrusion de mercure, adsorption d'azote, ...). Cependant, ces méthodes sont apparues suivant les cas incomplètes ou inadéquates pour fournir les caractéristiques exactes du réseau poreux des pâtes de ciment ou des bétons durcis, notamment dans le cas des matériaux à très hautes performances. La difficulté d'accès à la porosité fine des hydrates (gel C-S-H) est essentiellement responsable des insuffisances de ces méthodes. De plus, le séchage et le dégazage nécessaires, préalablement à la mesure, rendent l'interprétation des résultats délicate. Aussi, un autre type d'expérimentations plus pertinent a été mis en oeuvre : les expériences de désorption et d'adsorption de vapeur d'eau, à partir desquelles ont été obtenues les isothermes, c'est à dire les courbes "Teneur en Eau en fonction de l'Humidité Relative", à 'équilibre, à température constante (T = 23 °C) et pour chaque chemin de désorption/adsorption. L'analyse de ces courbes n'a pas seulement fourni les paramètres de texture (porosité, surface spécifique B.E.T., distribution poreuse,...) dans le domaine des mésopores, non accessible par beaucoup d'autres méthodes, elle a aussi permis l'étude des interactions eau/solide (épaisseur du film liquide adsorbé sur la surface solide, état d'énergie de l'eau des pores, hystérésis, ...). Les cinétiques de transfert d'humidité sont très lentes dans un tel milieu microporeux et elles sont fonction de l'épaisseur de l'échantillon testé ; l'équilibre hydrique a été établi après plusieurs mois avec nos échantillons de quelques millimètres d'épaisseur. La diffusivité a été calculée à partir des courbes de perte (ou gain) de masse en fonction du temps, pour chaque palier d'humidité relative. Au cours de ces expériences, les déformations hygrométriques ont aussi été mesurées, afin d'évaluer le retrait de dessiccation en fonction de l'humidité relative, et de le comparer au retrait endogène. La caractérisation hydrique a été complétée par la détermination de la répartition de la teneur en eau dans des éprouvettes en cours de séchage. A partir d'une large panoplie expérimentale, nous avons obtenu une description complète de la structure poreuse, depuis les pores relatifs au gel C-S-H, jusqu'à l'échelle des macropores. Nous avons mis en évidence les similarités et les différences, du point de vue microstructural et hydrique, entre les différents matériaux. Nous avons ainsi quantifié l'influence des paramètres de formulation tels que, par exemple, le rapport eau sur ciment ou la présence de fumée de silice.

<p>Les travaux de cette thèse s'inscrivent dans le cadre de la validation formelle et le contrôle réactif de calculs à haute performance sur systèmes-sur-puce (SoC). </p> <p>Dans ce contexte, la première contribution est la modélisation synchrone accompagnée d'une transformation d'applications en équations synchrones. Les modéles synchrones permettent de résoudre plusieurs questions liées à la validation formelle via l'usage des outils et techniques formels offerts par la technologie synchrone. Les transformations sont développées selon l'approche d'Ingénierie Dirigé par les Modèles (IDM). </p> <p>La deuxième contribution est une extension et amélioration des mécanismes de contrôle pour les calculs à haute performance, sous forme de constructeurs de langage de haut-niveau et de leur sémantique. Ils ont été défini afin de permettre la vérification, synthèse et génération de code. Il s'agit de déterminer un niveau d'abstraction de représentation des systèmes où soit extraite la partie contrôle, et de la modéliser sous forme d'automates à états finis. Ceci permet de spécifier et implémenter des changements de modes de calculs, qui se distinguent par exemple par les ressources utilisées, la qualité de service fournie, ou le choix d'algorithme remplissant une fonctionnalité. </p> <p>Ces contributions permettent l'utilisation d'outils d'analyse et vérification, tels que la vérification de propriétés d'assignement unique et dépendance acyclique, model checking. L'utilisation de techniques de synthèse de contrôleurs discrets est également traitée. Elles peuvent assurer la correction de faˆ on constructive: à partir d'une spécification partielle du contrôle, la partie manquante pour que les propriétés soient satisfaites est calculée. Grâce à ces techniques, lors du développement de la partie contrôle, la spécification est simplifiée, et le résultat est assuré d'être correct par construction. </p> <p>Les modélisations synchrone et de contrôle reposes sur MARTE et UML. Les travaux de cette thèse sont été partiellement implémentés dans le cadre de Gaspard, dédié aux applications de traitement de données intensives. Une étude de cas est présentée, dans laquelle nous nous intéressont à une application de système embarqué pour téléphone portable multimédia.</p>