Les centrales nucléaires jouent un rôle important au niveau mondial dans la production d'électricité aujourd'hui. Dans l'histoire humaine, deux accidents tragiques ont déjà eu lieu, à Tchernobyl et Fukushima, provoquant des fuites des matières radioactives. Pour éviter ce genre d'accident dans la future, la construction en béton prend un rôle important en tant que matériau de construction de confinement qui est la dernière barrière pour isoler la matière radioactive de l'environnement. La compréhension des transferts dans le béton dans des circonstances graves, en particulier à étudier le comportement du béton sous haute température jusqu'à 200°C. Dans la première partie, un modèle de déshydratation a été appliqué dans le modèle thermo-hygro. Le modèle thermo-hygro a ensuite été mis en oeuvre pour simuler une maquette dont les résultats expérimentaux sont disponibles dans la littérature, sans tenir compte des comportements mécaniques. Enfin, des études paramétrique ont été réalisées pour étudier l'influence de certains paramètres. Le modèle à ensuite été appliqué à des structures sous vitesse de chauffage différente, et le modèle à ensuite été appliqué à des structure sous rythme de chauffage différents et celles avec des épaisseurs différentes pour comparer les profiles de température et de pression de gaz à travers de la mur. / Nuclear power plant now takes an important part in the world's electricity generation. In human history, there have already been two tragic accidents in Chernobyl and Fukushima causing severe radioactive matter leakage. To pre- vent this kind of accident from happening again, concrete structure plays an important role as material of containment building, which is the last barrier to isolate the radioactive matter from the environment. Hence, the transfer properties of concrete under severe circumstances, especially high tempera- ture, are important for this usage. This PhD thesis aims to investigate the behavior of the concrete under high temperature up to 200°C. In the first part, a dehydration model was proposed. The model consists of different dehydration sub-models for main hydrates in the cement paste. In the second part, the dehydration model was implemented in a thermo-hygral model. The thermo- hygral model was then used to simulate a mock-up for which experimental results are available in the literature, without considering the mechanical behaviors. At last, parametric studies were performed to investigate the influ- ence of some parameters, and the model was then applied to structures under different heating rates, and structures with different thicknesses to compare the temperature and gas pressure profiles across a wall.

O estudo das propriedades de transporte do meio poroso é um tema importante para muitas áreas como a Ciência do Solo, onde o conhecimento dos processos que envolvem o movimento da água é de fundamental importância para o manejo e a conservação do mesmo. Nas últimas décadas, as técnicas não-invasivas e o processamento de imagens têm ajudado muito na modelagem e visualização da estrutura do meio, e têm sido aplicadas no estudo da distribuição dos fluidos com diferentes abordagens. Nesse trabalho é apresentado um novo método para análise do movimento da água, baseado na descrição da estrutura tridimensional da frente de molhamento no processo de infiltração vertical na região não-saturada do solo, considerando que a frente de molhamento tem informação relevante sobre os primeiros estágios da infiltração na interface líquido-sólido. O método proposto para a descrição da estrutura tridimensional da frente de molhamento compõe-se de dois principais módulos. O primeiro é responsável pela segmentação de imagens tomográficas para a detecção da frente de molhamento e cujo resultado é crucial para a análise da superfície resultante. O segundo módulo efetua a determinação de descritores da superfície obtida baseada na computação da variabilidade morfológica e a identificação de zonas de máxima adsorção através da análise da curvatura. As imagens usadas nos experimentos foram obtidas usando um tomógrafo de campo específico para estudos de solos, permitindo o processamento sem alterar a estrutura do solo. Os resultados preliminares são encorajadores e indicam que a abordagem utilizada consegue descrever o movimento da água usando informação da frente de molhamento no espaço tridimensional e no tempo. / The study of the transport properties in porous media is an important issue for many areas such as soil science, where knowledge about processes that involve the movement of water in the soil has fundamental importance to soil management and soil conservation. In recent decades noninvasive techniques and image processing algorithms have been very helpful in modeling and visualization of the structure medium and have been applied to study of the distribution of fluid with different approaches. This work present a new method to analysis of the movement of water based on the description of the three-dimensional wetting front structure in vertical infiltration process in unsaturated soil, whereas the wetting front structure has relevant information in the earliest stages of infiltration in liquid-solid interface. The proposed method for the description of the three-dimensional wetting front structure is comprised of two main modules. The first module is responsible of the three-dimensional image segmentation for the wetting front detection and its result is a very crucial step to analysis of the surface obtained. The second module compute features of the surface obtained to analysis based on morphological variability and maximal adsorption zones identification through the curvature. The image used in the experimental test were obtained using a tomograph of field specific to soil study, allowing the processing without changing of the soil structure. Our preliminary results are encouraging and indicate that our approach can describe the movement of water using information from the wetting front in three-dimensional space and time.

Este trabalho propõe um modelo de transferência de calor e massa para prever a formação de gelo em um container preenchido com legumes congelados. O problema físico é modelado como um meio poroso composto pelo próprio produto e o ar em seu entorno. O regime de convecção natural é assumido dentro do container, o qual promove o transporte de massa. Como uma primeira validação, o modelo é simulado considerando diferentes temperaturas de ar externo, causadas por flutuações da vizinhança. Resultados para quatro ciclos de temperaturas foram comparados, variando separadamente a temperatura média do ar, amplitude e frequência de oscilação. De modo geral, é observado que a temperatura do produto se comporta assim como era esperado e este resultado é diretamente associado à formação de gelo dentro do container. A formação de gelo cresce com uma maior amplitude de oscilação, porém decresce com um aumento na frequência e na temperatura média. Os parâmetros do modelo foram obtidos para dois diferentes produtos: fatias de cenouras congeladas e vagens congeladas, ambos em meio ao ar. As definições de parâmetros são oriundas de revisão bibliográfica, medições experimentais e simulações numéricas. Os parâmetros encontrados para a caracterização desses meios porosos foram similares para ambos os produtos, mesmo eles possuindo diferentes geometrias. A validação experimental foi feita para as fatias de cenoura considerando dois ciclos de temperatura O modelo numérico é capaz de prever o campo de velocidades do ar, as temperaturas do produto e a formação de gelo local. Os resultados foram validados em relação a um grupo independente de resultados numéricos, tal comparação apresentou uma boa concordância. A circulação de ar encontrada é, de fato, devido à convecção natural. O comportamento da temperatura dos produtos simulados concorda com os valores medidos e os valores de temperaturas diferem por menos de 12%. Com respeito à formação de gelo, o modelo é capaz de prevê-la corretamente nas regiões mais suscetíveis a este fenômeno. Porém, a quantidade de gelo formado prevista pelo modelo (1,56 g/semana) é menor do que a experimental (4,67 g/semana), apesar de serem de mesma ordem de magnitude. O efeito de cada parâmetro no modelo é estudado visando detectar maneiras de aprimorar o modelo. Foi encontrado que os parâmetros mais importantes para a formação de gelo total são a difusividade de massa efetiva e o coeficiente de transferência de calor convectivo dentro do container. Ajustando estes parâmetros duas vezes foi possível encontrar resultados melhores com respeito à formação de gelo (3,09 g/semana). / A model of heat and mass transfer is proposed in order to predict frost formation into a closed container filled with frozen vegetables. The physical problem is modeled as a macroporous media composed by the product itself and the surrounding air. Natural convection air flow is assumed into the container, who promotes water mass transport. As a first validation, the model is simulated for several exterior air temperatures, under environmental fluctuations (boundary conditions). Results of four temperature cycles were compared, varying average air temperature, amplitude and frequency of oscillation, one by one. As a general result, it is observed that the product temperature behavior is as expected, and it is directly associated with frost formation into the container. Frost formation increases with large amplitude of oscillation, but decreases with higher frequencies and higher mean temperatures. Model parameters were obtained for two assembling: frozen slices of carrots and air, and frozen extra thin green beans and air. Parameter definition and evaluation combines literature review, measurements and numerical simulation. In general, parameters which characterize these porous media were similar for both products, even though they display different geometries. The experimental validation is performed for carrot slices with two temperature cycles The numerical model is able to predict air velocity field, air and product temperatures, and local frost formation. Results are validated in respect to a set of independent experimental results that shown a good agreement. Air flow circulation is as expected due to natural convection. Product temperature simulated behavior agrees with measurements, and temperature values differ by less than 12%. Respect to frost formation predictions, the model predicts correctly the most susceptible regions to frost formation. However, the quantity of frost formed predicted by the model (1.56 g/ week)is lower than the experimental one (4.67g/week), despite being of the same order of magnitude. The effect of each parameter in the model is study in order to detect how to improve the model. The most important parameters affecting total frost formation are effective mass diffusivity and convective heat coefficient into the storage container. Adjusting these parameters to twice, better results in terms of frost formation could be obtained (3.09 g/ week).

Conselho Nacional de Desenvolvimento Científico e Tecnológico / A relação entre porosidade e permeabilidade desperta o interesse de pesquisadores e engenheiros por causa de suas diversas aplicações. Tais como na utilização de filtros, materiais pouco permeáveis, reservatórios naturais, etc. Ao longo do século XX, diversos trabalhos propondo tal relação foram apresentados na literatura e grande parte desses trabalhos desenvolvem modelos baseados na equação clássica de Kozeny-Carman. Nesta dissertação, propomos um modelo mais robusto que a formulação clássica de Kozeny-Carman, ou seja, que não apresenta as limitações dessa equação clássica. Além disso, um estudo baseado na Teoria dos Meios Fractais indica que o modelo estudado, nesta dissertação, generaliza diversas equações que fornecem a relação entre porosidade e permeabilidade. Por fim, será mostrado que o modelo proposto é capaz de descrever a relação entre porosidade e permeabilidade de diversos materiais porosos de natureza fractal. / The relationship between porosity and permeability attracts the attention of researchers and engineers because of their various applications. Such as in utilization of filters, waterproof materials, natural reservoirs, for example.Throughout the twentieth century, several works proposed in the literature they study the relation porosity-permeability, and much of this works they develop models based on the classical equation of Kozeny-Carman. In this dissertation, we propose a model more robust than the classical formulation of Kozeny-Carman, ie, that does not have the limitations of the equation classical. Furthermore, a study based on the Theory of the Media Fractals indicates that the model studied in this dissertation provide the relationship between porosity and permeability of several models presented in the literature. Finally, it will shown that the model proposed is able to describe the relationship between porosity and permeability of porous materials of various fractal nature.

abstract: The understanding of multiphase fluid flow in porous media is of great importance in many fields such as enhanced oil recovery, hydrology, CO2 sequestration, contaminants cleanup, and natural gas production from hydrate bearing sediments. In this study, first, the water retention curve (WRC) and relative permeability in hydrate bearing sediments are explored to obtain fitting parameters for semi-empirical equations. Second, immiscible fluid invasion into porous media is investigated to identify fluid displacement pattern and displacement efficiency that are affected by pore size distribution and connectivity. Finally, fluid flow through granular media is studied to obtain fluid-particle interaction. This study utilizes the combined techniques of discrete element method simulation, micro-focus X-ray computed tomography (CT), pore-network model simulation algorithms for gas invasion, gas expansion, and relative permeability calculation, transparent micromodels, and water retention curve measurement equipment modified for hydrate-bearing sediments. In addition, a photoelastic disk set-up is fabricated and the image processing technique to correlate the force chain to the applied contact forces is developed. The results show that the gas entry pressure and the capillary pressure increase with increasing hydrate saturation. Fitting parameters are suggested for different hydrate saturation conditions and morphologies. And, a new model for immiscible fluid invasion and displacement is suggested in which the boundaries of displacement patterns depend on the pore size distribution and connectivity. Finally, the fluid-particle interaction study shows that the fluid flow increases the contact forces between photoelastic disks in parallel direction with the fluid flow. / Dissertation/Thesis / Doctoral Dissertation Civil and Environmental Engineering 2016

Conselho Nacional de Desenvolvimento Científico e Tecnológico / A relação entre porosidade e permeabilidade desperta o interesse de pesquisadores e engenheiros por causa de suas diversas aplicações. Tais como na utilização de filtros, materiais pouco permeáveis, reservatórios naturais, etc. Ao longo do século XX, diversos trabalhos propondo tal relação foram apresentados na literatura e grande parte desses trabalhos desenvolvem modelos baseados na equação clássica de Kozeny-Carman. Nesta dissertação, propomos um modelo mais robusto que a formulação clássica de Kozeny-Carman, ou seja, que não apresenta as limitações dessa equação clássica. Além disso, um estudo baseado na Teoria dos Meios Fractais indica que o modelo estudado, nesta dissertação, generaliza diversas equações que fornecem a relação entre porosidade e permeabilidade. Por fim, será mostrado que o modelo proposto é capaz de descrever a relação entre porosidade e permeabilidade de diversos materiais porosos de natureza fractal. / The relationship between porosity and permeability attracts the attention of researchers and engineers because of their various applications. Such as in utilization of filters, waterproof materials, natural reservoirs, for example.Throughout the twentieth century, several works proposed in the literature they study the relation porosity-permeability, and much of this works they develop models based on the classical equation of Kozeny-Carman. In this dissertation, we propose a model more robust than the classical formulation of Kozeny-Carman, ie, that does not have the limitations of the equation classical. Furthermore, a study based on the Theory of the Media Fractals indicates that the model studied in this dissertation provide the relationship between porosity and permeability of several models presented in the literature. Finally, it will shown that the model proposed is able to describe the relationship between porosity and permeability of porous materials of various fractal nature.

O estudo das propriedades de transporte do meio poroso é um tema importante para muitas áreas como a Ciência do Solo, onde o conhecimento dos processos que envolvem o movimento da água é de fundamental importância para o manejo e a conservação do mesmo. Nas últimas décadas, as técnicas não-invasivas e o processamento de imagens têm ajudado muito na modelagem e visualização da estrutura do meio, e têm sido aplicadas no estudo da distribuição dos fluidos com diferentes abordagens. Nesse trabalho é apresentado um novo método para análise do movimento da água, baseado na descrição da estrutura tridimensional da frente de molhamento no processo de infiltração vertical na região não-saturada do solo, considerando que a frente de molhamento tem informação relevante sobre os primeiros estágios da infiltração na interface líquido-sólido. O método proposto para a descrição da estrutura tridimensional da frente de molhamento compõe-se de dois principais módulos. O primeiro é responsável pela segmentação de imagens tomográficas para a detecção da frente de molhamento e cujo resultado é crucial para a análise da superfície resultante. O segundo módulo efetua a determinação de descritores da superfície obtida baseada na computação da variabilidade morfológica e a identificação de zonas de máxima adsorção através da análise da curvatura. As imagens usadas nos experimentos foram obtidas usando um tomógrafo de campo específico para estudos de solos, permitindo o processamento sem alterar a estrutura do solo. Os resultados preliminares são encorajadores e indicam que a abordagem utilizada consegue descrever o movimento da água usando informação da frente de molhamento no espaço tridimensional e no tempo. / The study of the transport properties in porous media is an important issue for many areas such as soil science, where knowledge about processes that involve the movement of water in the soil has fundamental importance to soil management and soil conservation. In recent decades noninvasive techniques and image processing algorithms have been very helpful in modeling and visualization of the structure medium and have been applied to study of the distribution of fluid with different approaches. This work present a new method to analysis of the movement of water based on the description of the three-dimensional wetting front structure in vertical infiltration process in unsaturated soil, whereas the wetting front structure has relevant information in the earliest stages of infiltration in liquid-solid interface. The proposed method for the description of the three-dimensional wetting front structure is comprised of two main modules. The first module is responsible of the three-dimensional image segmentation for the wetting front detection and its result is a very crucial step to analysis of the surface obtained. The second module compute features of the surface obtained to analysis based on morphological variability and maximal adsorption zones identification through the curvature. The image used in the experimental test were obtained using a tomograph of field specific to soil study, allowing the processing without changing of the soil structure. Our preliminary results are encouraging and indicate that our approach can describe the movement of water using information from the wetting front in three-dimensional space and time.

Este trabalho propõe um modelo de transferência de calor e massa para prever a formação de gelo em um container preenchido com legumes congelados. O problema físico é modelado como um meio poroso composto pelo próprio produto e o ar em seu entorno. O regime de convecção natural é assumido dentro do container, o qual promove o transporte de massa. Como uma primeira validação, o modelo é simulado considerando diferentes temperaturas de ar externo, causadas por flutuações da vizinhança. Resultados para quatro ciclos de temperaturas foram comparados, variando separadamente a temperatura média do ar, amplitude e frequência de oscilação. De modo geral, é observado que a temperatura do produto se comporta assim como era esperado e este resultado é diretamente associado à formação de gelo dentro do container. A formação de gelo cresce com uma maior amplitude de oscilação, porém decresce com um aumento na frequência e na temperatura média. Os parâmetros do modelo foram obtidos para dois diferentes produtos: fatias de cenouras congeladas e vagens congeladas, ambos em meio ao ar. As definições de parâmetros são oriundas de revisão bibliográfica, medições experimentais e simulações numéricas. Os parâmetros encontrados para a caracterização desses meios porosos foram similares para ambos os produtos, mesmo eles possuindo diferentes geometrias. A validação experimental foi feita para as fatias de cenoura considerando dois ciclos de temperatura O modelo numérico é capaz de prever o campo de velocidades do ar, as temperaturas do produto e a formação de gelo local. Os resultados foram validados em relação a um grupo independente de resultados numéricos, tal comparação apresentou uma boa concordância. A circulação de ar encontrada é, de fato, devido à convecção natural. O comportamento da temperatura dos produtos simulados concorda com os valores medidos e os valores de temperaturas diferem por menos de 12%. Com respeito à formação de gelo, o modelo é capaz de prevê-la corretamente nas regiões mais suscetíveis a este fenômeno. Porém, a quantidade de gelo formado prevista pelo modelo (1,56 g/semana) é menor do que a experimental (4,67 g/semana), apesar de serem de mesma ordem de magnitude. O efeito de cada parâmetro no modelo é estudado visando detectar maneiras de aprimorar o modelo. Foi encontrado que os parâmetros mais importantes para a formação de gelo total são a difusividade de massa efetiva e o coeficiente de transferência de calor convectivo dentro do container. Ajustando estes parâmetros duas vezes foi possível encontrar resultados melhores com respeito à formação de gelo (3,09 g/semana). / A model of heat and mass transfer is proposed in order to predict frost formation into a closed container filled with frozen vegetables. The physical problem is modeled as a macroporous media composed by the product itself and the surrounding air. Natural convection air flow is assumed into the container, who promotes water mass transport. As a first validation, the model is simulated for several exterior air temperatures, under environmental fluctuations (boundary conditions). Results of four temperature cycles were compared, varying average air temperature, amplitude and frequency of oscillation, one by one. As a general result, it is observed that the product temperature behavior is as expected, and it is directly associated with frost formation into the container. Frost formation increases with large amplitude of oscillation, but decreases with higher frequencies and higher mean temperatures. Model parameters were obtained for two assembling: frozen slices of carrots and air, and frozen extra thin green beans and air. Parameter definition and evaluation combines literature review, measurements and numerical simulation. In general, parameters which characterize these porous media were similar for both products, even though they display different geometries. The experimental validation is performed for carrot slices with two temperature cycles The numerical model is able to predict air velocity field, air and product temperatures, and local frost formation. Results are validated in respect to a set of independent experimental results that shown a good agreement. Air flow circulation is as expected due to natural convection. Product temperature simulated behavior agrees with measurements, and temperature values differ by less than 12%. Respect to frost formation predictions, the model predicts correctly the most susceptible regions to frost formation. However, the quantity of frost formed predicted by the model (1.56 g/ week)is lower than the experimental one (4.67g/week), despite being of the same order of magnitude. The effect of each parameter in the model is study in order to detect how to improve the model. The most important parameters affecting total frost formation are effective mass diffusivity and convective heat coefficient into the storage container. Adjusting these parameters to twice, better results in terms of frost formation could be obtained (3.09 g/ week).

L'ingénierie tissulaire représente une solution prometteuse pour la production de substituts osseux. L'utilisation de bioréacteurs à perfusion pour cultiver des cellules ostéo-compétentes sur des matrices poreuses, permet de résoudre les limitations dues au transfert de masse, et d'apporter des stimuli physiques améliorant la prolifération et la différenciation cellulaire. Malgré les récents et importants développements des bioréacteurs en ingénierie tissulaire, les mécanismes menant à la production de substituts osseux en bioréacteurs restent mal compris.Le but de cette thèse est d'améliorer la compréhension de l'influence des phénomènes de transport, sur la croissance cellulaire et tissulaire dans un bioréacteur à perfusion. Dans cet objectif, une approche combinant modélisation et expérimentation est proposée.Dans un premier temps, un cadre théorique rigoureux est développé afin d'étudier les propriétés de transport du bioréacteur. Etant donné la nature hiérarchique du système, l'aspect multi-échelle du problème doit être pris en compte. En se basant sur la méthode de prise de moyenne volumique avec fermeture, les processus de transport d'espèce et de quantité de mouvement sont homogénéisés à partir de l'échelle de la matrice extracellulaire, jusqu'à l'échelle du bioréacteur. Les propriétés effectives des différentes structures rencontrées sont évaluées, et l'influence des dépendances inter-échelles sont mises en valeur. Le model macroscopique obtenu inclus des termes non-conventionnels, dont les contributions sont évaluées pour les conditions de fonctionnement du bioréacteur.Dans un second temps, la prolifération cellulaire et la production de tissu sont étudiées d'un point de vue expérimentale et théorique. Premièrement, des cellules de type fibroblaste, sont cultivées jusqu'à trois semaines sur des billes de verre, dans un bioréacteur perfusé à 10mL/min. Un protocole combinant des techniques d'histologie et d'analyse d'image, permet de quantifier la croissance de cellules et de tissu en fonction du temps et de l'espace. Deuxièmement, une cinétique de production de tissu est introduite dans le modèle de transport multiéchelle développé plus tôt. Finalement, la résolution à l'échelle du bioréacteur permet de discuter les résultats expérimentaux et théoriques au regard des phénomènes de transport ayant lieux dans le bioréacteur à perfusion. / Tissue engineering represents a promising approach for the production of bone substitutes. The use of perfusion bioreactors for the culture of bone-forming cells on a three-dimensional porous scaffold material, resolves mass transport limitations and provides physical stimuli, increasing the overall proliferation and differentiation of cells. Despite the recent and important development of bioreactors for tissue engineering, the underlying mechanisms leading to the production of bone substitutes remain poorly understood. The aim of this thesis is to gain insight on the influence of transport phenomena, on cell and tissue growth within a perfusion bioreactor. To this purpose, a combined modeling and experimental approach is followed.To start with, a rigorous theoretical framework is developed in order to study the transport properties of the bioreactor. Given the hierarchical nature of the system, the multiscale aspect of the problem must be taken into account. Based on the volume averaging theory with closure, mass and momentum transport processes are upscaled from the extracellular matrix scale, to the bioreactor scale. The effective properties of the encountered structures are evaluated, and the influence of the interscale dependencies are emphasized. The resulting macroscopic model includes non-conventional terms, which contributions are evaluated in the case of the bioreactor culture conditions.Then, cell proliferation and tissue growth are studied both, from an experimental and modeling point of view. First, fibroblast cells are cultured on glass beads in a bioreactor, perfused with culture medium at 10mL/min, for up to three weeks. A protocol combining histological techniques and image analysis allows the quantification of cell and tissue growth as a function of space and time. Second, a theoretical tissue production kinetic is introduced in the multiscale transport model previously developed. Finally, the resolution at the bioreactor scale allows to discuss the theoretical and experimental results in regard to the transport phenomena taking place in the perfusion bioreactor.

A Ressonância Magnética Nuclear (RMN) é uma importante técnica empregada nas principais áreas de conhecimento, tais como, Física, Química e Medicina. Importantes trabalhos da RMN aplicada ao estudo da dinâmica de moléculas em fluidos presentes em meios porosos permitiram que esta técnica ganhasse também notoriedade na indústria do petróleo. O presente projeto é fundamentado em alguns destes trabalhos seminais, reproduzindo, através de modelos físico-computacionais, os principais efeitos físicos da difusão e a consequente relaxação dos spins dos núcleos atômicos presentes nas moléculas dos fluidos imersos nos meios porosos. Os métodos teóricos utilizados para a interpretação dos parâmetros de relaxação transversal (T2) e longitudinal (T1), levam em consideração as propriedades macroscópicas da dinâmica de um ensemble de spins nucleares. Para a compreensão da difusão e relaxação nuclear em meios porosos, é de fundamental importância o emprego de métodos físico-computacionais que levam em consideração tanto a microestrutura quanto a difusão translacional dos spins nucleares através desses meios. Dentro dos modelos propostos nesta dissertação, foram utilizados conceitos básicos de mecânica estatística e de relaxação nuclear, resultando em programas envolvendo consideráveis recursos computacionais. Para a realização das simulações numéricas, foi necessária a construção da microestrutura digital dos meios porosos, usando como base do processo as imagens obtidas por microtomografia de raios-X 3D. Tendo disponível a microestrutura digital, o processo de difusão através do meio poroso e a consequente relaxação dos spins nucleares foram simulados através de modelos Random Walk, um exemplo clássico de Método de Monte Carlo. As simulações numéricas dos experimentos clássicos de RMN, empregados para a determinação de distribuições de tempos de relaxação transversal (conhecidas pela sigla CPMG), e mais avançados, para a determinação de taxas de troca entre diferentes sítios dos meios porosos (experimentos de troca ou Exchange), produziram importantes resultados para a compreensão da dinâmica de fluidos em meios porosos através de difusão translacional. / Nuclear Magnetic Resonance (NMR) is a versatile technique applied in many different areas of science such as Physics, Chemistry and Medicine. Several important NMR works applied in the study of the molecular dynamics of fluids inside porous media allowed this technique to become popular in the oil industry. This project is based on some of these seminal works. Using computational physics models, the main physical effects of diffusion and relaxation for the nuclear spins in an ensemble of molecules in a fluid permeating the porous media are reproduced. The theoretical methods used to interpret the transverse (T2) and longitudinal (T1) relaxation parameters take into consideration the macroscopic dynamical properties of an ensemble of nuclear spins. In order to comprehend the diffusion and relaxation in porous media, the use of computational physics methods that take into account the translational diffusion of the nuclear spins inside the media as well as its microstructure is essential. In the proposed models by this dissertation, basic statistical mechanics and nuclear relaxation concepts were used, resulting in programs that demand considerable computational resources. For the numerical simulations it was necessary to digitally construct the microstructure of the porous media and that was achieved with the use of X-ray 3D microtomography. On this digital microstructure, the diffusion process and the resulting relaxation of the nuclear spins were simulated with random walk models, an example of the Monte Carlo Method. The numerical simulation of classic NMR experiments, such as the one used to determine the distribution of relaxation times (known as CPMG), and more advanced ones, such as the one used to determine exchange rates between different sites (2D exchange experiments), produced fundamental results to the understanding of the fluid dynamics in porous media.