L'objectif de la présente thèse et de développer un cadre d'étude et de simulation apte à prédire les propriétés effectives de diffusion dans les matériaux cimentaires saturés. Le principe suivi, fondé sur les approches multi-échelles séquencées, permettra de modéliser chaque échelle d'intérêt, depuis le gel des C-S-H, jusqu'au matériau béton, et ainsi via des simulations de prédire les coefficients de diffusion effectifs.Le cadre d'étude est donc fondé sur les points clés suivant :- Des morphologies de type "matrice - inclusion" sont générées via des excursions de champs aléatoires corrélés tri-dimensionnels. Cette approche, combinée avec quelques opérations morphologiques, permet la génération de morphologies inclusionnaires 3D suivant des distributions de tailles imposées. Il est ensuite possible de faire évoluer ces morphologies, depuis un état initial, en fonction d'un modèle prescrit (par exemple en termes de fraction volumiques), via le seuil de l'excursion.- Les cibles morphologiques de chaque échelle considérées tentent de refléter les résultats expérimentaux et modèles les plus communément admis. En supplément, à l'échelle de la pâte de ciment, un modèle "simple" d'hydratation est développé afin de prédire l'évolution de celle-ci dans le temps en termes de fractions volumiques. Le choix a été fait d'utiliser pour fondement un modèle simple (Jennings & Tennis), incorporant des paramètres essentiels, tels que la composition minéralogique du ciment ou le rapport e/c. Des modifications mineures ont été apporté au modèle, afin de prendre en compte de récents résultats expérimentaux, concernant principalement la fin de l'hydratation ainsi que sa cinétique.- Le cadre d'homogénéisation a été construit pour être utilisable à toutes les échelles considérées, et il permet ainsi de prédire les propriétés effectives de diffusion avec des considération énergétiques d'équivalence entre échelles. De façon analogue à la théorie de Hill pour l'élasticité, le cadre d'étude assure une séparation d'échelle cohérente vis-à-vis des propriétés effectives de diffusion.- Enfin, d'un point de vue numérique, les simulations utilisent les méthodes E-FEM, couplées à l'utilisation de maillages non-adaptés sur lesquels les morphologies sont "projetées". De cette façon, les discontinuités géométriques (discontinuité de propriétés matériaux) sont sont assurés par les élément finis enrichis, et ne requièrent pas d'éléments joint, ou de "moyennage" de propriétés.Les résultats de chaque échelle sont utilisés comme paramètres d'entrés pour les échelles suivantes. Depuis l'échelle des hydrates (C-S-H), jusqu'à l'échelle du béton, les coefficients de diffusion effectifs sont ainsi estimés. / The aim of the present Ph.D. is to develop a methodology and a simulation framework able to predict as accurately as possible effective properties regarding diffusion in saturated cement based materials. A sequenced multiscale framework is developed to perform numerical homogenization on diffusive properties, and thus predict effective properties. The general outline is to model every scale in concrete-like materials, from the C-S-H gel, to the concrete scale, and with the help of a proper framework, be able to predict effective properties. The methodology relies on the following key points for every considered scale :- Inclusion-matrix like morphologies are generated through level set methods applied to multi-dimensional correlated random fields. This approach, together with few morphological operations allows to generate 3D morphologies with given particle size distributions (PSD), or pore size distributions. Then, it is also possible, from the initial state, to make those morphologies evolve according to any descriptive model (e.g. in terms of volume fraction), via the level set.- The morphological description of each scale is chosen according to the most acknowledged experimentations and models. In addition, at the cement paste scale, a hydration model was developed to predict the evolution of the cement paste through time, in terms of volume fractions. The choice was made to keep a simple model (Jennings & Tennis), with few but essential input parameters (mineral cement composition, w/c ratio, PSD). It has only been modified to include recent experimental results, such as the maximal hydration degree, or the hydration rates.- The upscaling framework was built to suit every considered scale, and allows to assess effective properties based on energy conservation principles. Similar to Hill's theory for elasticity, the framework is build to ensure to proper scale separation regarding effective diffusion coefficients.- Lastely, the numerical context is based on E-FEM techniques, which allows to use unstructured meshes, on which the morphologies are projected. This way, any material discontinuity in dealt within the FE, and thus does not require averaging properties, or joint elements. Results from every scale will be used as inputs at the higher scales, and effective diffusion coefficients are thus estimated.

Dans ce travail réalisé au Centre de Morphologie Mathématique and IFPEN, on s'intéresse à la microstructure et aux propriétés physiques d'alumines mésoporeuses. Il s'agit d'un supporte de catalyseur utilisés notamment dans les processus industriels de raffinage du pétrole. Fortement poreux, ce matériau est formé de ''plaquettes'' distribuées de manière désordonnée à l'échelle de la dizaine de nanomètres. Les propriétés de transport de masse du support de catalyseur sont fortement influencées par la morphologie de la microstructure poreuse. Ce travail porte sur la modélisation de la microstructure et des propriétés de transport des alumines mésoporeuses, à l'aide d'outils numériques et théoriques dérivés de l'analyse d'image et de la théorie des ensembles aléatoires. D'une part, on met en place des méthodes de caractérisation et de modélisation des microstructures, qui s'appuient sur, entre autre, des images obtenues par microscopie électronique en transmission (MET) et des courbes de porosimétrie azote. D'autre part, on utilise des méthodes d'homogénéisation numérique à champs complets par transformées de Fourier rapide (FFT).Dans un premier temps, le matériau est caractérisé expérimentalement par porosimétrie azote et résonance magnétique nucléaire à gradient de champ pulsé (RMN-GCP). Les images MET sont obtenus sur des échantillons d'épaisseur variable, filtrées et caractérisés par des fonctions de corrélation, notamment. Le bruit à haute fréquence issu de la membrane de carbone est identifié et pris en compte dans la modélisation de l'imagerie MET. À partir des images MET 2D, un modèle aléatoire à deux échelles est proposé pour représenter la microstructure 3D. Il prend en compte la forme des plaquettes d'alumines, leurs tailles, les effets d'alignement locaux et d'agrégation, qui sont identifiés numériquement. La procédure est validée à l'aide de comparaisons entre modèle et images expérimentales, en terme notamment de fonctions de corrélation et de surface spécifique, mesurées par porosimétrie azote.Dans un deuxième temps, une méthode de simulation des courbes d'isothermes de porosimétrie dans des milieux poreux périodiques ou aléatoires est développée. Basée sur des opérations morphologiques simples, elle étend un travail antérieur sur la porosimétrie au mercure. L'adsorption multicouche à basse pression est simulée à l'aide d'une dilatation tandis que les ménisques de l'interface vapeur-liquide intervenant pendant l'adsorption sont simulés à l'aide de fermetures de la phase solide par des éléments structurants sphériques. Pour simuler la désorption, une combinaison de fermetures et de bouchages de trou est utilisée. Le seuil de désorption est obtenu par une analyse de la percolation de la phase gazeuse. La méthode, d'abord validée sur des géométries simples, est comparée à des résultats antérieurs. Elle prédit une hystérésis et les distributions de pores associées à la porosimétrie. Nous l'appliquons aux modèles de microstructures 3D d'alumines mésoporeuses et proposons un modèle à trois échelles afin de rendre compte du seuil de pression pendant la désorption. En plus de la courbe de désorption, ce modèle reproduit les fonctions de corrélation mesurées sur les images MET.Dans un troisième temps, la diffusion de Fick, la perméabilité de Darcy, et les propriétés élastiques sont prédits à l'aide de calculs de champs complets par FFT sur des réalisations des modèles d'alumines mésoporeuses à deux et trois échelles. Les coefficients de diffusion effectifs et les facteurs de tortuosité sont prédits à partir de l'estimation du flux. Sont étudiés les effets de forme, d'alignement et d'agrégation des plaquettes sur les propriétés de diffusion à grande échelle. Les prédictions numériques sont validées au moyen des résultats expérimentaux obtenus par méthode RMN-GCP. / In a work made at Centre de Morphologie Mathématique and IFPEN, we study the microstructure and physical properties of mesoporous alumina. This is a catalyst carrier used in the petroleum refining industry. Highly porous, it contains disordered ''platelets'' at the nanoscale. The mass transport properties of the catalyst carrier are strongly influenced by the morphology of the porous microstructure. We focus on the modeling of the microstructure and of transport properties of mesoporous alumina, using numerical and theoretical tools derived from image analysis and random sets models. On the one hand, methods are developed to characterize and model the microstructure, by extracting and combining information from transmission electron microscope (TEM) images and nitrogen porosimetry curves, among others. On the other hand, the numerical homogenization relies on full-field Fourier transform computations (FFT).The material is first characterized experimentally by nitrogen porosimetry and pulse-field gradient nuclear magnetic resonance (PFG-NMR). TEM images, obtained on samples of various thicknesses are filtered and measured in terms of correlation function. The high-frequency noise caused by carbon membrane support is identified and integrated in the TEM image model. Based on the 2D TEM images, a two-scale random set model of 3D microstructure is developed. It takes into account the platelet shape, platelet size, local alignments and aggregations effects which are numerically identified. The procedure is validated by comparing the model and experimental images in terms of correlation function and specific surface area estimated by nitrogen porosimetry.Next, a procedure is proposed to simulate porosimetry isotherms in general porous media, including random microstructures. Based on simple morphological operations, it extends an earlier approach of mercury porosimetry. Multilayer adsorption at low pressure is simulated by a dilation operation whereas the menisci of the vapor-liquid interface occurring during adsorption are simulated by closing the solid phase with spherical structuring elements. To simulate desorption, a combination of closing and hole-filling operations is used. The desorption threshold is obtained from a percolation analysis of the gaseous phase. The method, validated first on simple geometries, is compared to previous results of the literature, allowing us to predict the hysteresis and pore size distribution associated to porosimetry. It is applied on 3D microstructures of mesoporous alumina. To account for the pressure threshold during desorption, we propose a refined three-scale model for mesoporous alumina, that reproduces the correlation function and the desorption branch of porosimetry isotherms.Finally, Fick diffusion, Darcy permeability, and elastic moduli are numerically predicted using the FFT method and the two-scale and three-scale models of mesoporous alumina. The hindering effects in diffusion are estimated by the Renkin's equation. The effective diffusion coefficients and the tortuosity factors are estimated from the flux field, taking into account hindering effects. The effects of platelet shape, alignment and aggregation on the diffusion property are studied. The numerical estimation is validated from experimental PFG-NMR results.