Morphological Modeling and Transport Properties of mesoporous alumina
Modélisation Morphologique et Propriétés de Transport d’Alumines Mésoporeuses
Résumé
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.
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 support de catalyseur utilisé notamment dans les procédés 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 obtenues 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'alumine, 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èles 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'isotherme 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édites à 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.
Loading...