Comprendre l'altération du verre nucléaire dans un réseau de fracture au sein d'un bloc de verre vitrifié est important pour la sûreté du conditionnement des déchets nucléaires (quantification des risques associés au relâchement des radionucléides). L’évaluation de la performance du stockage géologique des déchets nucléaires passe obligatoirement par la modélisation de l’altération aqueuse d’un bloc de verre nucléaire fracturé, l’échelle de temps envisagée (plusieurs milliers d’années) dépassant toute possibilité d’expérience directe. Cette thèse vise donc à combler le fossé entre les simulations d'écoulement et de transport à l'échelle du réservoir et la modélisation à l'échelle micrométrique des processus interfaciaux verre-eau, en apportant l'évaluation quantitative de la dégradation aqueuse du verre à l'échelle d’un bloc.Pour aborder ce problème, les objectifs principaux de cette thèse ont été fixés comme suit : (i) la reproduction des résultats expérimentaux obtenus précédemment (pour quelques fractures modélisées de manière discrète en mode diffusif), (ii) l’analyse de l'impact des géométries de fractures sur la quantité de verre altéré pour quelques fissures modélisées de manière discrète, (iii) l’étude de la possibilité d'adaptation du modèle géochimique à la modélisation dans le cadre de l’approche milieu équivalent, (iv) la mise au point d'une méthodologie de caractérisation, (v) la modélisation géostatistique et géométrique de réseau de fractures à l’échelle d’un conteneur de verre, (vi) le calcul des paramètres équivalents diffusifs, hydrauliques et les paramètres qui contrôlent la cinétique de dissolution de verre, et au final, (vii) la modélisation de transport réactif à l’échelle d’un conteneur.À titre illustratif, la méthodologie de la caractérisation de réseau fracturé proposée, basée sur le traitement des images, a été appliquée aux images bidimensionnelles (2D) de haute résolution de deux blocs de verre. Cette application a permis de mettre en œuvre à la fois les données directes obtenues par mesures des paramètres d’un réseau fracturé de verre vitrifié et les données indirectes explicatives issues des simulations thermomécaniques. L’application a abouti à la création de multiples réalisations de tessellation de réseaux fracturés équivalents qui ont ensuite été utilisées comme représentations physiques pour les calculs de la perméabilité équivalente, de la diffusion équivalente et des paramètres contrôlant la cinétique de dissolution de verre borosilicaté. L'évolution de la quantité de verre altéré obtenue en effectuant la modélisation de transport réactif appliquée à plusieurs réalisations de la tessellation de réseau fracturé équivalent a été comparée aux données expérimentales d’un essai d'altération aqueuse d'un conteneur non radioactif de verre nucléaire. Les résultats montrent que la méthodologie conçue offre une opportunité pour mieux comprendre l'impact de la fracturation sur l'altération aqueuse du verre vitrifié et constitue un outil fiable permettant de prendre en compte différents scénarios d'évolution du stockage. / Understanding the alteration of nuclear glass in a fracture network of a vitrified glass block is important for the safe conditioning of nuclear waste (quantification of the risks associated with radionuclide release). Performance assessment of geological nuclear waste repositories entails modelling of the long-term evolution of the fractured nuclear glass block aqueous alteration, because the considered time scale, of several thousands of years, is beyond the range of any direct experimental perspectives. This dissertation aims then to bridge the gap between the reservoir-scale flow and transport simulations and the micron-scale modeling of the glass-water interfacial processes, by bringing the quantitative evaluation of the glass aqueous degradation at the block scale.To tackle this issue, the main objectives of this thesis were fixed as follows: (i) reproduction of the experimental results previously obtained (for some fractures modeled in a discrete way in the diffusive mode),(ii) analysis of the impact of fractures geometries on the quantity of altered glass at the scale of some fractures modeled in a discrete way, (iii) investigation of the possibilities of the geochemical model adaptation for the equivalent homogenous modeling, (iv) establishment of a methodology for glass block fracture network characterization, (v) geostatistical and geometric modeling, (vi) calculation of the equivalent diffusive, hydraulic and glass dissolution kinetics controlling properties and (vii) upcoming reactive transport modeling at the scale of one canister.As an illustrative example, the proposed image processing-based fracture network characterization methodology was applied to two-dimensional (2D) high-resolution images of two blocks of vitrified glass. This application brought into service both hard data obtained by direct measurement of the fracture network and soft physics-based explanatory data and resulted in the creation of multiple realizations of fracture network equivalent tessellation that were further used as physical representation for the calculation of the equivalent hydraulic, diffusive, and alteration kinetics - controlling properties. The evolution of the quantity of altered glass obtained by conducting reactive transport modeling applied to several realizations of the equivalent fracture network tessellation was compared with the experimental data of the aqueous alteration test of a non-radioactive full-scale nuclear glass canister. The results show that implementation of the devised procedure presents an opportunity for better understanding the impact of fracturing on aqueous alteration of borosilicate glass and provides a reliable tool enabling different scenarios of repository evolution to be accounted for.

Les campagnes de recherche - comprenant la conception de nouveaux matériaux, la prospection, la réalisation de composites et leur caractérisation - sont longues. Pour l'ingénieur ou le chercheur, relier la microstructure et les propriétés physiques macroscopiques des matériaux le plus en amont, permet de réduire la durée et le coût de ces études. Nous nous intéressons aux propriétés diélectriques d'un composite particulier, formé d'une matrice époxy contenant des inclusions d'agrégats de noir de carbone. Nous détaillons le processus global, de la réalisation du composite à sa modélisation. Des échantillons de composite sont réalisés en faisant varier des paramètres de préparation. Des mesures diélectriques, électriques et d'autres mesures physiques situent le niveau des propriétés macroscopiques des échantillons. L'observation du composite par microscopie électronique à transmission à plusieurs grandissement permet d'appréhender sa morphologie. Les micrographies sont traitées par analyse d'images pour caractériser la géométrie du composite en quelques grandeurs. Les lois statistiques de changement d'échelle permettent, à partir de données sur la microstructure, d'obtenir un encadrement de la permittivité de composites. Les lois statistiques de changement d'échelle sont une méthode de modélisation récente, elles permettent en particulier :

Les traitements d'acidification sont généralement utilisés pour stimuler l'injection et/ou la productivité des puits dans les formations de carbonates. Cependant, de nombreux traitements ne produisent pas les résultats attendus en terme de gain de productivité à cause de la mauvaise modélisation autour du puits des mécanismes physiques intervenant durant le processus d'injection acide. La nature instable du phénomène de dissolution en milieu poreux entraîne la formation de canaux fortement conductifs appelés wormholes, qui sont difficiles à modéliser quantitativement. Un modèle de dissolution à l'échelle de Darcy est proposé comprenant une équation de Darcy-Brinkman pour la partie écoulement couplée avec un modèle de dissolution en non-équilibre local. Un simulateur numérique 3D a été développé pour résoudre ce système d'équations en utilisant une méthode de pas fractionnaire et des schémas TVD. Les résultats sont présentés sur des configurations 2D et 3D aussi bien pour des systèmes homogènes qu'hétérogènes. Les résultats numériques sont discutés d'un point de vue qualitatif et quantitatif par rapport à la littérature et comparés aux résultats expérimentaux. Les expériences ont été réalisées sur un massif de sel quasi 2D dans lequel on a injecté une solution d'eau salée sous-saturée. Les instabilités de dissolution, le développement des canaux et la propagation des wormholes ont été enregistrés à l'aide d'une caméra vidéo. En se basant sur les résultats 2D, la possibilité d'une description du phénomène à l'échelle de la section, c'est-à-dire en effectuant des moyennes sur des sections du domaine, a été explorée. Plusieurs approches sont utilisées, tels que les modèles à une équation, où le milieu considéré incorpore la physique du wormhole et la matrice poreuse environnante, et les modèles à deux équations pour lesquels les wormholes sont traités séparément. Les implications théoriques sont discutées en se basant sur les résultats numériques.

Nous présentons des modèles macroscopiques du comportement mécanique du parenchyme pulmonaire humain obtenus par homogénéisation double--échelle. Dans une première partie consacrée au couplage entre parenchyme et arbre bronchique, nous commençons par proposer un modèle de la déformation du parenchyme. Nous modélisons (i) le parenchyme par un matériau élastique linéaire, (ii) les alvéoles comme des cavités de taille epsilon réparties périodiquement dans le domaine macroscopique et (iii) l'arbre bronchique par un arbre dyadique résistif en supposant la loi de Poiseuille valide pour chaque voie aérienne. Nous obtenons en faisant tendre epsilon vers zéro, par homogénéisation, une description macroscopique du parenchyme comme un matériau viscoélastique, sous certaines conditions sur l'irrigation du domaine par l'arbre que nous étudions ensuite. L'arbre induit une dissipation non--locale en espace. Nous illustrons ces résultats par des résultats numériques. Dans une deuxième partie, nous étudions la propagation d'ondes sonores dans le parenchyme, sans prendre en compte l'arbre bronchique. Nous homogénéisons dans le domaine fréquentiel un premier modèle couplant l'élasticité linéarisée dans le parenchyme avec l'équation acoustique dans l'air. Nous déduisons une loi de comportement macroscopique élastique linéaire en dehors d'un ensemble de résonances. Ensuite, nous homogénéisons un deuxième modèle qui prend en compte le caractère viscoélastique et inhomogène du parenchyme au niveau microscopique. Le matériau macroscopique présente des effets de mémoire nouveaux par rapport aux composants microscopiques que nous étudions numériquement.

Les structures planaires de grandes tailles sont de plus en plus utilisées dans les applications des satellites et des radars. Deux grands types de ces structures à savoir les FSS et les Reflectarrays sont particulièrement les plus intéressants dans les domaines de la conception RF. Mais en raison de leur grande taille et de la complexité des cellules élémentaires, l'analyse complète de ces structures nécessite énormément de mémoire et des temps de calcul excessif. Par conséquent, les techniques classiques basées sur maillage linéaire soit ne parviennent pas à simuler de telles structures soit, exiger des ressources non disponibles à un concepteur d'antenne. Une technique appelée " technique par changement d'échelle " tente de résoudre ce problème par partitionnement de la géométrie du réseau par de nombreux domaines imbriqués définis à différents niveaux d'échelle du réseau. Le multi-pôle par changement d'échelle, appelé " Scale changing Network (SCN) ", modélise le couplage électromagnétique entre deux échelles successives, en résolvant une formulation intégral des équations de Maxwell par une technique basée sur la méthode des moments. La cascade de ces multi-pôles par changement d'échelle, permet le calcul de la matrice d'impédance de surface de la structure complète qui peut à son tour être utilisées pour calculer la diffraction en champ lointain. Comme le calcul des multi-pôles par changement d'échelle est mutuellement indépendant, les temps d'exécution peuvent être réduits de manière significative en parallélisant le calcul. Par ailleurs, la modification de la géométrie de la structure à une échelle donnée nécessite seulement le calcul de deux multi-pôles par changement d'échelle et ne requiert pas la simulation de toute la structure. Cette caractéristique fait de la SCT un outil de conception et d'optimisation très puissant. Des structures planaires uniformes et non uniformes excité par un cornet ont étés modélisés avec succès, avec des temps de calcul dé lais intéressants, employant les ressources normales de l'ordinateur.

Antenne de l'avenir surtout pour les applications spatiales sont de plus en plus complexe en raison de la nécessité de reconfigurabilité en termes de fréquence, la puissance rayonnée et reçues, la consommation d'énergie et de la fiabilité. Dans ce contexte, réflecteurs et des surfaces sélectives en fréquence (FSS) sont particulièrement les plus intéressantes dans les domaines de la conception RF. Mais en raison de leur grande taille et de la complexité des cellules élémentaires, l'analyse complète de ces structures nécessite énormément de mémoire et des temps de calcul excessif spécialement lorsque des éléments de réglage tels que RF-MEMS sont également intégrés au sein des cellules. Par conséquent, les techniques classiques basées sur maillage linéaire soit ne parviennent pas à simuler de telles structures soit, exiger des ressources non disponibles à un concepteur d'antenne. Une technique appelée "technique par changement d'échelle" tente de résoudre ce problème par partitionnement de la géométrie du réseau par de nombreux domaines imbriqués définis à différents niveaux d'échelle du réseau. Le multi-pôle par changement d'échelle, appelé "Scale Changing Network (SCN)", modélise le couplage électromagnétique entre deux échelles successives, en résolvant une formulation intégral des équations de Maxwell par une technique basée sur la méthode des moments. La cascade de ces multi-pôles par changement d'échelle, permet le calcul de la matrice d'impédance de surface de la structure complète qui peut à son tour être utilisées pour calculer la diffraction en champ lointain. Comme le calcul des multi-pôles par changement d'échelle est mutuellement indépendant, les temps d'exécution peuvent être réduits de manière significative en parallélisant le calcul. Par ailleurs, la modification de la géométrie de la structure à une échelle donnée nécessite seulement le calcul de deux multi-pôles par changement d'échelle et ne requiert pas la simulation de toute la structure. Cette car actéristique fait de la SCT un outil de conception et d'optimisation très puissant. Des structures planaires uniformes et non uniformes excité par un cornet ont étés modélisés avec succès, avec des temps de calcul délais intéressants, employant les ressources normales de l'ordinateur. Avec l'intérêt croissant pour des réflecteurs hautement reconfigurables pour les applications spatiales, différentes technologies pour réaliser la reconfigurabilité du diagramme de rayonnement ont été explorées ces dernières années. L'une d'elles consiste à imprimer un déphasage approprié en différent point du front d'onde incident à l'aide de cellule déphaseuse à MEMS-RF. En raison de la grande diversité des échelles dans de telles structures reconfigurables, les outils classiques de simulation électromagnétique se révèlent souvent inefficaces en termes des ressources informatiques et de temps de calcul. Les approches basées sur un circuit électrique équivalent se sont avérées plus utiles pour les concepteurs.

Le béton se présente comme un matériau constitué de granulats, jouant le rôle d'inclusions, enchâssés dans une matrice correspondant à la pâte de ciment hydraté. A l'échelle de la pâte, l'hydratation du ciment génère un milieu multiphasique, constitué d'un squelette solide et de pores remplis ou partiellement remplis d'eau selon leur taille. Le composant principal du squelette solide est le gel de C-S-H, les autres composants étant de nature cristalline. La qualification de gel du composant C-S-H est liée à sa nanostructure dont la schématisation la plus admise consiste en une phase aqueuse adsorbée, en sandwich avec des feuillets solides de nature cristalline. Il est bien admis que la structure du C-S-H est à l'origine de son comportement viscoélastique et donc de celui du béton. Ce comportement viscoélastique peut s'expliquer par un réarrangement de sa nanostructure sous l'effet des contraintes mécaniques appliquées à l'échelle macroscopique. La modélisation macroscopique du fluage du béton ne permet pas d'expliquer la variabilité du fluage d'une formulation de béton à une autre. En effet, les paramètres des modèles macroscopiques ne peuvent être identifiés que par l'analyse de résultats expérimentaux obtenus par des essais réalisés sur des éprouvettes de béton. Ces paramètres ne sont valables que pour une formulation donnée. L'identification de ces paramètres conduit donc à des programmes expérimentaux très coûteux et ne fournit pas suffisamment d'informations sur la sensibilité des paramètres macroscopiques à la variabilité des caractéristiques mécaniques et morphologies des constituants. Dans ce travail, on suppose qu'il existe une échelle microscopique à laquelle les mécanismes moteurs du fluage ne sont pas impactés par la formulation du béton. A cette échelle, celle du C-S-H, les propriétés viscoélastiques peuvent être considérées avoir un caractère intrinsèque. L'influence de la formulation ne concerne alors que les concentrations des différents hydrates. Trois approches, analytiques, semi-analytiques et numériques sont alors proposées pour estimer, par une homogénéisation multi-échelle, les propriétés viscoélastiques macroscopiques du béton à partir des propriétés de ses constituants ainsi qu'à partir de sa microstructure. Ces approches sont basées sur l'extension des schémas d'homogénéisation élastique au cas viscoélastique au moyen du principe de correspondance qui utilise la transformée de Laplace-Carson. Les propriétés effectives sont alors déterminées directement dans l'espace de Carson. Par la suite, celles dans l'espace temporel sont obtenues par la transformée inverse. Les approches proposées apportent des solutions aussi bien dans un cadre général que sous certaines hypothèses restrictives : coefficient de Poisson viscoélastique microscopique ou macroscopique constant, module de compressibilité constant. Sur le plan théorique, deux schémas d'homogénéisation ont été étudiés : le schéma de Mori-Tanaka, basé sur le problème de l'inclusion d'Eshelby, et le schéma auto-cohérente généralisé basé sur la neutralité énergique de l'inclusion. Les résultats obtenus montrent que sous ces hypothèses restrictives, le spectre macroscopique se présente comme une famille de sous ensembles de temps caractéristiques bornés par les temps caractéristiques microscopiques. Par ailleurs, les propriétés thermodynamiques, de croissance monotone et de concavité, des fonctions de retard macroscopiques ne sont préservées par l'homogénéisation que sous certaines conditions de compatibilité des spectres microscopiques. Sur le plan pratique, les méthodes développées ont été appliquées pour construire la complaisance de fluage propre macroscopique du béton en connaissant les données communes de toutes sortes de bétons et celles correspondant à une formulation donnée. Les résultats expérimentaux disponibles sont alors exploités pour analyser le caractère intrinsèque des propriétés viscoélastiques à l'origine du fluage du béton

Dans cette thèse nous nous sommes intéressés à la formation des bancs de poissons et à l'étude des interactions proies-prédateurs en présence de comportements collectifs par une approche de modélisation. Ce phénomène, bien qu'il soit le résultat d'interactions qui ont lieu à l'échelle individuelle, il engendre des conséquences à plus grande échelle, spatiale ou temporelle. L'objectif principal de cette thèse a été celui de comprendre l'influence du processus d'agrégation sur la réponse fonctionnelle. Nous avons élaboré un modèle centré sur l'individu (IBM) qui décrit les interactions intra-spécifiques ainsi que celle entre proies et prédateurs. Les agents peuvent ou pas avoir du comportement collectif.Dans un premier modèle les prédateurs sont attirés par les proies qu'ils voient. Ensuite, nous avons ajouté une composante à la vitesse du prédateur dépendante du bruit émis par les proies. Pour les deux cas, nous avons comparé les réponses fonctionnelles dans quatre configurations différentes. Nos résultats suggèrent que la prédation est plus efficace lorsque les proies s'agrègent. De plus, deux différents types de réponse fonctionnelle émergent : celle de Holling type II si les proies ne forment pas des bancs, celle de Holling type III autrement.Ensuite, nous nous sommes focalisés sur les conséquences à l'échelle de la population d'un phénomène (l'agrégation) se déroulant à l'échelle individuelle. Nous avons cherché un indice capable de détecter les bancs dans l'espace. Ensuite, en suivant une démarche mathématique, nous avons écrit une équation aux dérivées partielles représentant l'évolution spatio-temporelle de la densité des proies. / This thesis deals with the modelling of schooling phenomena and its consequences on predator-prey dynamics. Many marine species exhibit collective behaviour. While this phenomenon depends on individual interactions, it can have important effects at larger spatial and temporal scales.The main goal of this work is to understand the influence of aggregation on functional response, which represents predator feeding behaviour. We elaborated an individual-based model (IBM) describing schooling behaviour as well as predator-prey interactions. Predators can either be attracted towards visible prey, or hear, and so move towards, the noise produced by prey at a larger distance and then attack once they are close enough to see them. We analysed four different configurations, in which prey and/or predators can school or just move randomly. Our results shown an increased predation efficiency when prey school, and also different functional response shapes: Holling type II emerges if prey do not aggregate, while Holling type III emerges when prey school. Then, we focused on schooling consequences at higher scale. In particular, we analysed how to fit classical models to our emergent functional responses. Moreover, we found out a possible way to detect aggregates. Finally, we wrote a model representing the spatio-temporal evolution of prey density.

Le béton se présente comme un matériau constitué de granulats, jouant le rôle d'inclusions, enchâssés dans une matrice correspondant à la pâte de ciment hydraté. A l'échelle de la pâte, l'hydratation du ciment génère un milieu multiphasique, constitué d'un squelette solide et de pores remplis ou partiellement remplis d'eau selon leur taille. Le composant principal du squelette solide est le gel de C-S-H, les autres composants étant de nature cristalline. La qualification de gel du composant C-S-H est liée à sa nanostructure dont la schématisation la plus admise consiste en une phase aqueuse adsorbée, en sandwich avec des feuillets solides de nature cristalline. Il est bien admis que la structure du C-S-H est à l'origine de son comportement viscoélastique et donc de celui du béton. Ce comportement viscoélastique peut s'expliquer par un réarrangement de sa nanostructure sous l'effet des contraintes mécaniques appliquées à l'échelle macroscopique. La modélisation macroscopique du fluage du béton ne permet pas d'expliquer la variabilité du fluage d'une formulation de béton à une autre. En effet, les paramètres des modèles macroscopiques ne peuvent être identifiés que par l'analyse de résultats expérimentaux obtenus par des essais réalisés sur des éprouvettes de béton. Ces paramètres ne sont valables que pour une formulation donnée. L'identification de ces paramètres conduit donc à des programmes expérimentaux très coûteux et ne fournit pas suffisamment d'informations sur la sensibilité des paramètres macroscopiques à la variabilité des caractéristiques mécaniques et morphologies des constituants. Dans ce travail, on suppose qu'il existe une échelle microscopique à laquelle les mécanismes moteurs du fluage ne sont pas impactés par la formulation du béton. A cette échelle, celle du C-S-H, les propriétés viscoélastiques peuvent être considérées avoir un caractère intrinsèque. L'influence de la formulation ne concerne alors que les concentrations des différents hydrates. Trois approches, analytiques, semi-analytiques et numériques sont alors proposées pour estimer, par une homogénéisation multi-échelle, les propriétés viscoélastiques macroscopiques du béton à partir des propriétés de ses constituants ainsi qu'à partir de sa microstructure. Ces approches sont basées sur l'extension des schémas d'homogénéisation élastique au cas viscoélastique au moyen du principe de correspondance qui utilise la transformée de Laplace-Carson. Les propriétés effectives sont alors déterminées directement dans l'espace de Carson. Par la suite, celles dans l'espace temporel sont obtenues par la transformée inverse. Les approches proposées apportent des solutions aussi bien dans un cadre général que sous certaines hypothèses restrictives : coefficient de Poisson viscoélastique microscopique ou macroscopique constant, module de compressibilité constant. Sur le plan théorique, deux schémas d'homogénéisation ont été étudiés : le schéma de Mori-Tanaka, basé sur le problème de l'inclusion d'Eshelby, et le schéma auto-cohérente généralisé basé sur la neutralité énergique de l'inclusion. Les résultats obtenus montrent que sous ces hypothèses restrictives, le spectre macroscopique se présente comme une famille de sous ensembles de temps caractéristiques bornés par les temps caractéristiques microscopiques. Par ailleurs, les propriétés thermodynamiques, de croissance monotone et de concavité, des fonctions de retard macroscopiques ne sont préservées par l'homogénéisation que sous certaines conditions de compatibilité des spectres microscopiques. Sur le plan pratique, les méthodes développées ont été appliquées pour construire la complaisance de fluage propre macroscopique du béton en connaissant les données communes de toutes sortes de bétons et celles correspondant à une formulation donnée. Les résultats expérimentaux disponibles sont alors exploités pour analyser le caractère intrinsèque des propriétés viscoélastiques à l'origine du fluage du béton / Concrete can be considered as a multiscale heterogeneous material in which elastic inclusions, i.e., aggregates, are embedded in a viscoelastic matrix which corresponds to the hardened cement paste. At the cement paste scale, the hydration process leads to a partially saturated porous multiphase medium whose solid skeleton is mainly constituted of the C-S-H gel, the other hydrates being crystalline solids. The gel-like properties of the C-S-H are essentially related to the existence of strongly adsorbed structural water in the interlayer nanospace. It is well admitted that this feature is the microscopic origin of the creep mechanism of C-S-H and therefore of concrete at the macroscale. This viscoelastic behavior can be explained by a rearrangement of the nanostructure of the C-S-H component due to a macroscopically applied loading. Therefore, the macroscopic modeling of concrete creep does not allow explaining the variability of creep with respect to the concrete mix. Indeed, the identi?cation of the viscoelastic parameters of macroscopic models requires performing creep tests at the concrete scale. Thus, the validity of the identi?ed parameters is limited to the concrete mix under consideration and the identi?cation process does not permit to carry out the influence of the microstructure and the mechanical constituent properties on the macroscopic creep. In this work, we assume that there exists a microscopic scale at which the driving creep mechanisms are not affected by the concrete mix. At this scale, the C-S-H one, the viscoelastic properties can be considered to be intrinsic. The concrete mix will essentially influence the volume fraction of the hydrates. Analytical, semi-analytical and numerical approaches are developed in order to derive, by multi-scale homogenization, the macroscopic viscoelastic properties of concrete on the basis of the knowledge of the microscopic properties, of the constituents together with the knowledge of the mix-dependent microstructure. These approaches consist in extending results of elastic homogenization schemes to viscoelasticity with the help of the correspondence principle using the Laplace-Carson transform. Therefore, the effective properties are obtained in a straightforward manner in the Carson transform space, which requires the transform inversion in order to determine the time evolution of these properties. The proposed models permit to answer to this requirement, both under some restrictive assumptions (constant Poisson's ratio at either microscopic or macroscopic scale, constant bulk microscopic modulus) and in a general way. At the theoretical level, two homogenization schemes are investigated : the Mori-Tanaka scheme based on the Eshelby solution and the Generalized Self Consistent scheme, based on the energetic neutrality condition of the inclusion. The obtained results show that under the aforementioned restrictive conditions, the macroscopic creep spectrum presents as a family of sets of retardation times which are bounded by two successive microscopic characteristic times. Furthermore, the monotonically increasing and concavity features of the macroscopic creep compliance, which derive from thermodynamic considerations, are not always preserved, except if some compatibility conditions are ful?lled by the microscopic spectrums. From the practical point of view, these methods developed are used to establish the complaisance creep function of concrete from the knowledge of the common data of all types of concrete, as well as the other parameters that characterize a given formulation. Hence, by a back-analysis of creep tests obtained on cement paste or on concrete, the intrinsic basic creep properties of the gel C-S-H are estimated

Les aquifères fissurés de type karst contiennent d'importantes ressources en eau. Ces aquifères sont complexes et hétérogènes sur une gamme d'échelles importantes. Leur gestion nécessite l'utilisation d'outils et de méthodologies adaptés. Dans le cadre de cette étude, différents outils et méthodologies numériques d'étude ont été développés pour la modélisation des aquifères karstiques, et plus généralement, des milieux poreux très fissurés 2D et 3D - en mettant l'accent sur la morphologie et sur le comportement hydrodynamique du milieu à travers la notion de changement d'échelle ("second changement d'échelle", reposant sur un modèle d'écoulement local de type Darcy et/ou Poiseuille avec quelques généralisations). Plusieurs axes sont explorés concernant la morphologie du milieu poreux fissuré (milieux aléatoires, milieux booléens avec réseaux statistiques de fissures, mais aussi, modèles morphogénétiques). L'étude du changement d'échelle hydrodynamique tourne autour du concept de macro perméabilité. Dans un premier temps, l'étude porte sur un modèle de perte de charge linéaire darcien. Les perméabilités effectives sont calculées numériquement en termes des fractions volumiques de fissures et du contraste de perméabilité matrice/fissures. Elles sont analysées et comparées à des modèles théoriques (analytiques). Une étude particulière des effets de quasi-percolation pour les grands contrastes aboutit à la définition de trois fractions critiques liées à des seuils de percolation. Pour tenir compte des effets inertiels dans les fissures, l'étude est étendue au cas d'une loi locale comprenant un terme quadratique en vitesse (Darcy/Ward-Forchheimer). Une perméabilité macroscopique équivalente non linéaire est définie et analysée à l'aide d'un modèle inertiel généralisé (linéaire/puissance). Enfin, l'anisotropie hydraulique à grande échelle du milieu fissuré est étudiée, en termes de perméabilités directionnelles, à l'aide d'une méthode numérique d'immersion. / Karstic aquifers contain large subsurface water resources. These aquifers are complex and heterogeneous on a large range of scales. Their management requires appropriate numerical tools and approaches. Various tools and numerical methodologies have been developed to characterize andmodel the geometry and hydraulic properties of karstic aquifers, more generally, of highly fissured 2D and 3D porous media. In this study, we emphasize morphological characterization, and we analyze hydrodynamic behavior through the concept of upscaling ("second upscaling"). Concerning the morphology of fissured porous media, several axes are explored : random media, composite random Boolean media with statistical properties, and morphogenetic models. Hydrodynamic upscaling is developed using the macro-permeability concept. This upscaling method is based on either Darcy's linear law, or on a linear/quadratic combination of Darcy's and Ward-Forchheimer's quadratic law (inertial effects). First, the study focuses on Darcy's linear head loss law, and Darcian effective permeabilities are calculated numerically in terms of volume fractions of fissures and "fissure/matrix" permeability contrasts. The results are analysed and compared with analytical results and bounds. A special study of percolation and quasi-percolation effects, for high contrasts, leads to defined three critical fractions. These critical fractions are "connected" to percolation thresholds. Secondly, in order to consider inertial effect in fissures, the study is extended to a local law with a quadratic velocity term (Darcy/Ward-Forchheimer). Then, an equivalent nonlinear macroscopic permeability is defined and analysed using a generalized inertial model (linear/power). Finally, the large scale hydraulic anisotropy of fissured medium is studied, in terms of directional permeabilities, using an "immersion" numerical method.