La modélisation de la rupture des géo-matériaux constitue un important défi pour les applications telles que la séquestration du CO2 , le stockage de déchets nucléaires, la production des hydrocarbures ainsi que les projets de génie civil concernant les tunnels ou les excavations. L'objectif de cette thèse est de développer des lois d'évolution macroscopiques d'endommagement à partir des descriptions explicites de la rupture à l'échelle microscopique en vue de la modélisation du comportement d'endommagement à long terme des sites de stockage géologique. L'approche adoptée est basée sur l'homogénéisation par développements asymptotiques et la description énergétique de la propagation des micro-fissures, qui permettent l'obtention des lois d'endommagement et conduisent à une quantification explicite de l'énergie de l'émission acoustique associée à la rupture. Les modèles obtenus sont capables de prédire la dégradation des modules d'élasticité en raison de l'évolution des micro-fissures. Cette représentation permet de modéliser la propagation des ondes dans un milieu à endommagement évolutif. Deux types de modèles d'endommagement seront proposés: indépendants de temps et dépendants de temps. Les modèles dépendants de temps décrivent l'évolution progressive quasi-fragile de la micro-fissuration. Dans les modèles dépendants de temps, l'évolution des micro-fissures est décrite à travers un critère sous-critique et la propagation mixte, par branchement. En utilisant le modèle dépendant de temps, des simulations seront faites à trois niveaux: du laboratoire, du tunnel et du réservoir. / Accurate modeling of failure of geomaterials is the key to the success of a diverse range of engineering challenges including the topic of CO2 sequestration, nuclear waste disposal and hydrocarbon production plus civil engineering projects for tunnels or excavations. The aim of this thesis is to develop macroscopic damage evolution laws based on explicit descriptions of fracture at the micro-scale level which can be successfully employed to describe long term damage behavior of geologic storage sites. The approach taken is based on homogenization through asymptotic developments combined with micro-crack propagation energy analysis which leads to an explicit quantification of the acoustic emission (AE) energy associated with damage. Proposed damage models are capable of modeling the degradation of elastic moduli due to the micro-crack evolution. This representation allows the modeling of wave propagation in a medium with evolving damage. Two types of damage models will be considered: time independent and time dependent. Time independent damage models capable of describing progressive micro-cracking propagation (i.e. quasi-brittle type damage law) are considered. In the case of time-dependent damage models, the evolution of the micro-crack length during propagation is described through a sub-critical criterion and mixed mode propagation by branching. Using the time dependent damage model including rotational micro-cracks, simulations will be made at three levels: laboratory, tunnel and reservoir scales.

Dans un objectif de réduction de l'utilisation des pesticides, des méthodes alternatives doivent assurer la protection des cultures. Les dégâts causés par les insectes ravageurs sont les manifestations de dynamiques de population dont certains déterminants paysagers sont des éléments structuraux potentiellement relocalisables. Nous explorons ici, par la simulation numérique, la modification du paysage comme stratégie de protection des cultures.Les éléments paysagers peuvent impacter les populations d'insectes par leur influence sur les comportements individuels. Nous avons développé un modèle de population tenant compte de l'approvisionnement individuel. Ce processus résulte (1) de la perception des ressources distantes par l'individu, (2) de la localisation de la population sur une dimension additionnelle quantifiant les réserves énergétiques, et (3) d'une procédure d'optimisation qui définit des mouvements dirigés adaptés.Une heuristique évolutionnaire de modification de paysages a été développée. Elle produit des réarrangements sous contraintes de composition et de structure. Les modèles de paysages et de dynamique de population ont été appliqués à un problème théorique d'aménagement du paysage. La description mécaniste de l'approvisionnement optimal offre des leviers d'action efficaces, basés sur la perturbation des comportements d'approvisionnement. Nos travaux montrent la réactivité de la population aux modifications du paysage, ainsi que la capacité des algorithmes évolutionnaires à proposer des paysages résistants, sous contraintes agronomiques. / Environmental and health issues call for a switch in crop protection towards less chemically driven strategies. Pest damages on crops result of population dynamics that are influenced by landscape features. Those features may be relocated or dismissed to improve landscapes resistance to pest damages. Here we explore numerically the landscape modification as a crop protection strategy.Landscape features can influence population dynamics through their impact on individuals. Therefore, we developed a multiscale population dynamics model accounting for individual foraging. The foraging process results of (1) the perception of distant features used as resources by the individuals, (2) the localisation of the population along an additional dimension quantifying individual energy supply, (3) an optimisation procedure defining adapted directed motionsWe developed an evolutionary heuristic for landscape modification. It is able to rearrange landscapes with respect to compositional and structural constraints. Population and landscape models were applied to a theoretical landscape planning problem. The mechanistic description of the optimal foraging process enables new and efficient levers for crop protection, building on the disruption of the foraging behaviours. This application enlightens the responsiveness of the simulated population to landscape modifications, as well as the ability of evolutionary algorithms to produce resilient landscapes under agronomic productive constraints.

En France métropolitaine, les projets de géothermie haute température pour la production d’électricité sont principalement localisés dans le socle des fossés d’effondrement liés à la mise en place du Rift Ouest Européen. Le socle de ces fossés a été étudié sur deux analogues à l’affleurement sur les épaules du rift : les Vosges du Nord pour le fossé Rhénan et la bordure Est du Massif central pour la fosse de Valence. Cette étude a permis de montrer que le réseau de failles s’organise selon trois ordres de grandeurs de longueurs et d’espacements caractéristiques qui individualisent des blocs structuraux. Les orientations et l’espacement des failles formant ces blocs et la présence ou l’absence de certains ordres de grandeurs sont le résultat de l’héritage anté-rift du socle, ainsi que du mécanisme d’ouverture du bassin. Le potentiel réservoir des formations de socle et de la couverture surincombante a été analysé au regard de ces zones de failles et de l’altération supergène qui affecte le toit du socle. Ainsi, les lithologies potentiellement rencontrées en base des fossés ont pu être classées en fonction du potentiel de développement de porosité et de perméabilité matricielle dans les cœurs de failles, les zones endommagées et le réseau pervasif de fractures dans le protolithe. L’évolution de la fracturation dans les zones de failles a également pu être appréhendé, et une méthodologie de modélisation double milieu a été élaborée pour caractériser la porosité et la perméabilité de fractures et modéliser le fonctionnement d’un doublet géothermique dans une faille synthéthique équivalente aux cibles des projets géothermiques / High temperature geothermal projects for electricity production are in France mostly localized in the basement of basins linked to the West European Rifting event. The basement of theses basins have been studied on two outcrop analogues at the shoulders of the rift: the Northern Vosges mountains for the Upper Rhine Graben and the Eastern border of the Massif central for the Valence Graben. This study has shown that the fault network is organized in three orders of size with characteristic length and spacing, and that form characteristic structural blocks. The orientation and spacing of these faults and the presence or absence of some size orders are the result of structural inheritance of the basement and the mechanism of the basin opening. The reservoir potential of the basement rocks and the surrounding sedimentary cover has been analysed in light of the fault zones structure and the weathered layer at the top of the basement. Thus, the basement rocks of the basins has been classified in light of the potential of matrix porosity and permeability development in the fault core, the damaged zone and the fractured protolith. The evolution of the fracture network in the fault zone has been studied too, with the development of a double-porosity model in order to characterize the fracture porosity and permeability, and therefore simulate the working of a geothermal doublet in a synthethic fault zone analogue of the geothermal drilling targets

L’os est un matériau vivant avec une structure hiérarchique complexe qui lui confère des propriétés mécaniques remarquables. L’os subit perpétuellement des contraintes mécaniques et physiologiques, ainsi sa qualité et sa résistance à la fracture évoluent constamment au cours du temps à travers le processus de remodelage osseux. La qualité osseuse est non seulement définie par la densité minérale osseuse mais également par les propriétés mécaniques ainsi que la microarchitecture. Dans le cadre de la présente thèse, on a développé une modélisation multiéchelle unifiée couplant à la fois les activités cellulaires au comportement mécanique de l'os tenant compte des différents niveaux hiérarchiques de l'os: de l’ultrastructure au remodelage osseux. Ce modèle permet d’étudier le comportement mécano-bibliologique de l’os et de prédire ses propriétés mécaniques apparentes à différentes échelles allant du nanoscopique au macroscopique en fonction des constituants élémentaires de l'os. Pour atteindre cet objectif, une démarche en quatre phases a été adoptée. La première phase consiste à décrire les constituants élémentaires de l’os. La deuxième phase avait pour objectif la modélisation multiéchelle de l'ultrastructure osseuse constituée de trois échelles nanoscopiques (microfibrille, fibrille et fibre) par la méthode des éléments finis et des réseaux de neurones. La troisième phase correspond à la modélisation des échelles micro-macroscopiques de l’os cortical (lamelle, ostéon, os cortical) en utilisant comme paramètres d’entrée les propriétés de la fibre déterminées dans la deuxième phase. Enfin, dans la dernière phase, on a développé un modèle mécano-biologique du remodelage osseux permettant de simuler le processus d'adaptation osseuse tenant compte explicitement des activités biologiques des cellules osseuses. Les propriétés mécaniques prédites par nos algorithmes multiéchelles ont servi pour alimenter le modèle de remodelage. Ce modèle a été implémenté au code de calcul d’éléments finis ABAQUS/Standard à travers sa routine utilisateur UMAT. Finalement, le modèle EF mécano-biologique multiéchelle du remodelage osseux a été appliqué pour simuler différents scénarii de remodelage sur des fémurs humains (2D et 3D). Différents facteurs ont été ainsi analysés tels que l'âge, le genre, l'amplitude des activités physiques, etc. Les résultats obtenus sont conformes (qualitativement) avec les observations cliniques et cohérents avec les différentes études expérimentales. En conclusion: (i) Les modèles unifiés ainsi développés (modèle multiéchelle, modèle mécano-biologique de remodelage osseux) contribuent à l'analyse fine du comportement de l'os humain. (ii) L'application des algorithmes a permis d'effectuer des essais virtuels pour analyser les effets combinés de nombreux facteurs caractérisant la qualité osseuse. / Bone is a living material with a complex hierarchical structure which entails exceptional mechanical properties. Bone undergoes permanent mechanical and physiological stresses, thus its quality and fracture toughness are constantly evolving over time through the process of bone remodeling. Bone quality is not only defined by bone mineral density but also by the mechanical properties and microarchitecture. The current thesis offers a multiscale modeling approach unifying the cell activity to the mechanical behavior, taking into consideration the hierarchical levels of bone, from the ultrastructure to bone remodeling. This model permits to study the mechanobiological behavior and to predict the mechanical properties of the bone at different scales from nano to macro depending on the elementary constituents of bone. To achieve the objective of the current work, an approach of four phases was adopted. The first phase is to describe the basic components of the bone. The second phase concerns the multiscale modeling of the three nanoscopic levels of bone ultrastructure (microfibril, fibril and fiber) by the finite element method and neural networks. The third phase aims to model the micro-macroscopic structures of cortical bone (lamella, osteon, cortical bone) using the fiber properties predicted from the second phase as input parameters. In the last phase, a mechano-biological model of bone remodeling was achieved to simulate the process of bone adaptation explicitly considering the biological activities of bone cells. Mechanical properties predicted by our multiscale algorithms were used to feed the remodeling model. This model has been implemented into the ABAQUS/Standard finite elements code as a user subroutine. Finally, the finite element mechano-biological multiscale model of bone remodeling was applied to simulate different scenarios on human femurs (2D and 3D). Hence, different factors such as: age, gender, physical activities, etc were analyzed. The obtained results are conformed (qualitatively) to clinical observations and consistent with the various experimental studies. In summary, (i) the models portrayed here (multiscale model, mechanical-biological model of bone remodeling) contribute by their unified approach to the realistic modeling of the response of human bone. (ii) The application of the algorithms permits to perform virtual experiments to scrutinize the combined effects of numerous factors dictating the bone quality.

Représentations séparées pour la simulation multi-échelle du comportementmécanique et de l’endommagement des matériaux composites.Résumé: Le développement de méthodes numériques performantes pour simuler les structurescomposites est un défi en raison de la nature multi-échelle et de la complexité des mécanismed’endommagement de ce type de matériaux. Les techniques classiques de discrétisationvolumique conduisent à des coûts de calcul importants et sont restreintes en pratique à despetites structures.Dans cette thèse, un nouvelle stratégie basée sur une représentation séparée de la solution estexplorée. L’objectif est de proposer un cadre numérique efficace et fiable pour analyser les endommagementsdans les composites stratifiés sous chargements statiques et dynamiques. Ladécomposition propre généralisée (PGD) est utilisée pour construire la solution.Pour traiter l’endommagement, et plus particulière le délaminage, un modèle de zone cohésivea été implémenté dans la PGD. Une approches multi-échelle innovante est également proposéepour simuler le comportement mécanique des composites à microstructure périodique. L’idéeprincipale est de séparer deux échelles : l’échelle du motif périodique (microstructure) et l’échellemacroscopique. Les résultats de la PGD sont très proches des résultats obtenus par la méthodeéléments fini classique. Finalement, la PGD permet de réduire significativement la complexitédes modèles tout en gardant une précision satisfaisante. / Separated representations for the multiscale simulation of the mechanicalbehavior and damages of composite materials.Abstract: The development of efficient simulations for composite structures is very challengingdue to the multiscale nature and the complex damage process of this materials. When usingstandard 3D discretization techniques with advanced models for large structures, the computationalcosts are generally prohibitive.In this thesis, a new strategy based on a separated represenation of the solution is explored todevelop a computationally efficient and reliable numerical framework for the analysis of damagesin laminated composites subjected to quasi-static and dynamic loading. The PGD (Proper GeneralizedDecomposition) is used to build the solution.To treat damage, and especially delamination, a cohesive zone model has been implemented inthe PGD solver. A novel multiscale approach is also proposed to compute the mechanical behaviorof composites with periodic microstructure. The idea is to separate two scales: the scaleof periodic pattern and the macroscopic scale. The PGD results have been compared with theresults obtained with the classcial finite element method. A close agreement is found between thetwo approach and the PGD has significantly reduced the model complexity.

L’amélioration du caractère prédictif des simulations numériques requiert une meilleure maîtrise des modèles constitutifs. Des procédures d’identification exploitant des mesures de champs ont été développées afin de valider les modèles et d’identifier un jeu de paramètres constitutifs à partir d’un nombre réduit d’essais hétérogènes. Bien évidemment, les résultats de telles procédures dépendent grandement des incertitudes des mesures. Dans cette thèse, une stratégie d’identification de paramètres constitutifs à partir de mesures de champs par corrélation d’images numériques éléments finis (CIN-EF) multi-échelles est proposée et développée. Une variante régularisée du recalage par éléments finis (FEMU-R) est adaptée à cette approche multi-échelles. Elle exploite des champs de déplacements mesurés par CIN-EF d'une part à l’échelle de la structure (fournissant les conditions aux limites pour la simulation EF) et d’autre part à une échelle plus locale (fournissant des champs de déplacements mieux résolus pour la comparaison essais/calculs). Un algorithme de recalage d’images est développé pour faire le lien entre les échelles. D’abord validée à l’aide d’images de synthèse, la procédure est ensuite appliquée à un essai de traction sur plaque trouée réalisé sur un composite stratifié verre/époxy. On montre qu’une telle stratégie permet de diminuer nettement non seulement les incertitudes de mesure mais également les incertitudes d’identification. Finalement, la question d’un mouchetis adapté à cette approche multi-échelles est évoquée et des éléments de solutions sont proposés et testés. / Improving the prediction of numerical simulations requires a better control of constitutive models. Identification methods exploiting full-field measurements have been developed in order to validate models and to identify a set of constitutive parameters from a reduced number of heterogeneous tests. The results of those methods largely depends on measurement uncertainties. In this PhD thesis, a strategy is proposed and developed for the identification of constitutive parameters from full-field measurements obtained by multiscale finite element digital image correlation (FE-DIC). A Regularised Finite Element Model Updating method (FEMU-R) is adapted to this multiscale approach. It exploits displacement fields measured by FE-DIC at a structural scale (in order to obtain the boundary conditions for FE simulation)and at a local scale (giving a better resolution on the displacement field for the comparison between experiment and simulation). An image registration algorithm is developed to bridge these scales. First validated on synthetic images, the multiscale approach is then applied to an open-hole tensile test of a glass/epoxy composite laminate. It is shown that such a strategy allows to reduce not only the measurement uncertainties but also the identification uncertainties. Finally the question of a speckle pattern suitable for this approach is discussed and some first technical solutions are proposed and tested.

Émissions polluantes, les motoristes souhaitent contrôler au mieux l’atomisation du carburant, injecté généralement sous forme de jets ou de nappes liquides. Les essais étant long et coûteux, leur remplacement par un outil numérique capable de simuler le processus d’atomisation permettrait non seulement une réduction des coûts importante mais faciliterait également la phase de conception. Toutefois, en raison du caractère multi-échelle du phénomène, il est difficile de le décrire dans son ensemble avec les approches habituellement utilisées en mécanique des fluides numérique.L’objectif de cette thèse est de concevoir une nouvelle approche qui permettra à terme de simuler l’atomisation pour une configuration industrielle complète. Celle-ci consiste à coupler deux types de modèles. Le premier, dit modèle bifluide, est un modèle à deux fluides compressibles basé sur les équations de Navier-Stokes diphasiques. Celui-ci permet de décrire les grandes échelles du phénomène d’atomisation correspondant à la formation de ligaments et d’amas liquides dans la zone proche de l’injecteur. Le second, dit modèle de spray, est basé sur une équation cinétique. Dans la zone située en aval de l’injecteur, ce dernier permet de décrire de manière statistique l’évolution du brouillard de gouttelettes issues de la fragmentation primaire du jet de carburant. Le point délicat, à la fois sur le plan de la modélisation et sur celui de l’algorithmique, réside dans le couplage des deux modèles. Celui ci a été réalisé grâce à l’introduction de deux modèles auxiliaires permettant de traiter le transfert de liquide entre le modèle bifluide et le modèle de spray par atomisation ou ré-impact.L’approche proposée a été appliquée à la simulation numérique de nappes liquides cisaillées. Les comparaisons entre les résultats numériques et des résultats expérimentaux montrent que le modèle bifluide permet de prévoir l’influence de la géométrie et des conditions d’injection sur l’atomisation primaire de la nappe liquide. Le modèle d’atomisation permet quant à lui, de reproduire le caractère instationnaire des mécanismes de production de gouttes lors du transfert de la phase liquide depuis le modèle bifluide vers celui de spray. Des cas de ré-impact valident également la robustesse et la généralité de la méthodologie de couplage. / In order to improve efficiency of aircraft combustion chambers and reduce polluting emissions,engine manufacturers try to achieve a better control on fuel atomization, which is usually injectedas a jet or liquid sheet. As experiments are expensive and time consuming, a numerical tool able to simulate atomization would be a powerful asset in engine conception design. However, simulation ofthe whole atomization process with commonly used approach in computational fluid dynamics is still prohibitive due to the multi-scale nature of the phenomenon.The objective of this thesis is to develop a new approach allowing the simulation of the spray formation for a industrial configuration in the near future. This involves coupling of two types of models.The first one, called two-fluid model, is based on the Navier-Stokes equations for two immiscible compressible fluids. This one is used to describe the large scales of the atomization mechanism corresponding to the formation of ligaments and liquids blobs in the near-injector area. The second one,called spray model, is based on a kinetic equation. Further downstream from the injector, this model describes statistically the evolution of the droplet cloud produced by the primary fragmentation of liquid jet. The main difficulty, in terms of both modeling and algorithmic, is the coupling of these twomodels.This has been achieved by introducing an atomization and an impact models which ensure liquid transfer between the two-fluid model and the spray model.This new approach was applied to the numerical simulation of sheared liquid sheets. Comparisons between numerical and experimental results show how the two-fluid model predicts the influence of injector geometry and injection conditions on the primary atomization of the liquid sheet. Concerning droplets production, the atomization model is able to reproduce the unsteady nature of this mechanism when transferring liquid phase from the two-fluid model to the spray model. Test cases for the impact model also validate the robustness and generality of the coupling approach.

abstract: This work investigates the effects of non-random sampling on our understanding of species distributions and their niches. In its most general form, bias is systematic error that can obscure interpretation of analytical results by skewing samples away from the average condition of the system they represent. Here I use species distribution modelling (SDM), virtual species, and multiscale geographically weighted regression (MGWR) to explore how sampling bias can alter our perception of broad patterns of biodiversity by distorting spatial predictions of habitat, a key characteristic in biogeographic studies. I use three separate case studies to explore: 1) How methods to account for sampling bias in species distribution modeling may alter estimates of species distributions and species-environment relationships, 2) How accounting for sampling bias in fossil data may change our understanding of paleo-distributions and interpretation of niche stability through time (i.e. niche conservation), and 3) How a novel use of MGWR can account for environmental sampling bias to reveal landscape patterns of local niche differences among proximal, but non-overlapping sister taxa. Broadly, my work shows that sampling bias present in commonly used federated global biodiversity observations is more than enough to degrade model performance of spatial predictions and niche characteristics. Measures commonly used to account for this bias can negate much loss, but only in certain conditions, and did not improve the ability to correctly identify explanatory variables or recreate species-environment relationships. Paleo-distributions calibrated on biased fossil records were improved with the use of a novel method to directly estimate the biased sampling distribution, which can be generalized to finer time slices for further paleontological studies. Finally, I show how a novel coupling of SDM and MGWR can illuminate local differences in niche separation that more closely match landscape genotypic variability in the two North American desert tortoise species than does their current taxonomic delineation. / Dissertation/Thesis / Doctoral Dissertation Geography 2018

L'analyse des systèmes de grande taille est confrontée à des difficultés d'ordre syntaxique et sémantique : comment observer un million d'entités distribuées et asynchrones ? Comment interpréter le désordre résultant de l'observation microscopique de ces entités ? Comment produire et manipuler des abstractions pertinentes pour l'analyse macroscopique des systèmes ? Face à l'échec de l'approche analytique, le concept d'émergence épistémique - relatif à la nature de la connaissance - nous permet de définir une stratégie d'analyse alternative, motivée par le constat suivant : l'activité scientifique repose sur des processus d'abstraction fournissant des éléments de description macroscopique pour aborder la complexité des systèmes. Cette thèse s'intéresse plus particulièrement à la production d'abstractions spatiales et temporelles par agrégation de données. Afin d'engendrer des représentations exploitables lors du passage à l'échelle, il apparaît nécessaire de contrôler deux aspects essentiels du processus d'abstraction. Premièrement, la complexité et le contenu informationnel des représentations macroscopiques doivent être conjointement optimisés afin de préserver les détails pertinents pour l'observateur, tout en minimisant le coût de l'analyse. Nous proposons des mesures de qualité (critères internes) permettant d'évaluer, de comparer et de sélectionner les représentations en fonction du contexte et des objectifs de l'analyse. Deuxièmement, afin de conserver leur pouvoir explicatif, les abstractions engendrées doivent être cohérentes avec les connaissances mobilisées par l'observateur lors de l'analyse. Nous proposons d'utiliser les propriétés organisationnelles, structurelles et topologiques du système (critères externes) pour contraindre le processus d'agrégation et pour engendrer des représentations viables sur les plans syntaxique et sémantique. Par conséquent, l'automatisation du processus d'agrégation nécessite de résoudre un problème d'optimisation sous contraintes. Nous proposons dans cette thèse un algorithme de résolution générique, s'adaptant aux critères formulés par l'observateur. De plus, nous montrons que la complexité de ce problème d'optimisation dépend directement de ces critères. L'approche macroscopique défendue dans cette thèse est évaluée sur deux classes de systèmes. Premièrement, le processus d'agrégation est appliqué à la visualisation d'applications parallèles de grande taille pour l'analyse de performance. Il permet de détecter les anomalies présentes à plusieurs niveaux de granularité dans les traces d'exécution et d'expliquer ces anomalies à partir des propriétés syntaxiques du système. Deuxièmement, le processus est appliqué à l'agrégation de données médiatiques pour l'analyse des relations internationales. L'agrégation géographique et temporelle de l'attention médiatique permet de définir des évènements macroscopiques pertinents sur le plan sémantique pour l'analyse du système international. Pour autant, nous pensons que l'approche et les outils présentés dans cette thèse peuvent être généralisés à de nombreux autres domaines d'application. / The analysis of large-scale systems faces syntactic and semantic difficulties: How to observe millions of distributed and asynchronous entities? How to interpret the disorder that results from the microscopic observation of such entities? How to produce and handle relevant abstractions for the systems' macroscopic analysis? Faced with the failure of the analytic approach, the concept of epistemic emergence - related to the nature of knowledge - allows us to define an alternative strategy. This strategy is motivated by the observation that scientific activity relies on abstraction processes that provide macroscopic descriptions to broach the systems' complexity. This thesis is more specifically interested in the production of spatial and temporal abstractions through data aggregation. In order to generate scalable representations, the control of two essential aspects of the aggregation process is necessary. Firstly, the complexity and the information content of macroscopic representations should be jointly optimized in order to preserve the relevant details for the observer, while minimizing the cost of the analysis. We propose several measures of quality (internal criteria) to evaluate, compare and select the representations depending on the context and the objectives of the analysis. Secondly, in order to preserve their explanatory power, the generated abstractions should be consistent with the background knowledge exploited by the observer for the analysis. We propose to exploit the systems' organisational, structural and topological properties (external criteria) to constrain the aggregation process and to generate syntactically and semantically consistent representations. Consequently, the automation of the aggregation process requires solving a constrained optimization problem. We propose a generic algorithm that adapts to the criteria expressed by the observer. Furthermore, we show that the complexity of this optimization problem directly depend on these criteria. The macroscopic approach supported by this thesis is evaluated on two classes of systems. Firstly, the aggregation process is applied to the visualisation of large-scale distributed applications for performance analysis. It allows the detection of anomalies at several scales in the execution traces and the explanation of these anomalies according to the system syntactic properties. Secondly, the process is applied to the aggregation of news for the analysis of international relations. The geographical and temporal aggregation of media attention allows the definition of semantically consistent macroscopic events for the analysis of the international system. Furthermore, we believe that the approach and the tools presented in this thesis can be extended to a wider class of application domains.

In this dissertation, we study the stationary heat equation in a heterogeneous tridimensional plate, using a "dimension reduction" techinique called hierarchical modelling and we generate model the original problem in a two-dimensional domain. To estimate the error modelling, we develop an asymptotic expansion for the original problem solution and for the aproximate solution. Comparing both solutions with their own asymptotic expansions, we obtain an estimative of the error modelling. We perform some computational experiments, using the Residual Free Bubbles (RFB) Method and the Multiscale Finite Element Method for the diffusion problem and for the diffusion-reaction problem in a two-dimensional domain, with small parameters. Finally, we extend the numeric solutions found the original tridimensional problem. / Neste trabalho, estudamos a equação do calor estacionária em uma placa heterogênea tridimensional. Para a modelagem deste problema, utilizamos uma técnica de redução de dimensão conhecida por Modelagem Hierárquica. Desta forma, geramos um modelo para o problema original em um domínio bidimensional. Com o objetivo de estimar o erro de modelagem, desenvolvemos a expansão assintótica da solução do problema original e da solução aproximada. Comparando as soluções com suas respectivas expansões assintóticas, obtemos uma estimativa para o erro de modelagem. Realizamos alguns experimentos computacionais, desenvolvendo o método Residual Free Bubbles (RFB) e o método de Elementos Finitos Multiescala (MEFM) para o problema de difusão e para o problema de difusão-reação em um domínio bidimensional, com parâmetros pequenos. Com base nestes experimentos, encontramos algumas soluções numéricas para o problema da placa tridimensional.