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.

Many growing southern Ontario communities, including Guelph, rely on fractured bedrock aquifers for drinking water. Contamination and overexploitation pose a threat to these water resources, necessitating characterization of vulnerability, risks, and recharge areas. Quaternary sediments southeast of the City of Guelph, including the Paris Moraine, were investigated in order to delineate hydrostratigraphy. This was achieved through study of 9 cored-holes, as well as existing MOE, GRCA, and University of Guelph data. Falling head permeameter measurements and empirical grain-size distribution measurements and analysis were employed for determination of 𝐾𝑠𝑎𝑡 values, which were used to construct a hydrostratigraphy. Of 19 methods evaluated, The Kozeny-Carman empirical grain-size method for determining 𝐾𝑠𝑎𝑡 was found to be more representative of measured values for the study area. The area is dominated by a conductivity regime of 2.72x10-7 – 1.40x10-6 m/s with local heterogeneity present on the scale of 10’s to 100’s of meters. The Paris Moraine, particularly its backslope, is at higher risk due to its relatively high conductivity, greater occurrence of aquifer units, as well as prevalence of small-scale topographic (hummocky topography), and bedrock topographic lows. / Thesis / Master of Science (MSc)

Ce travail présente deux nouvelles approches pour aider à la caractérisation de la topologie des milieux poreux en termes de distribution de tailles de pores. La première des méthodes proposées repose sur les propriétés rhéologiques des fluides à seuil (type Bingham, Casson....) s'écoulant à travers l'échantillon poreux. Par une simple mesure du débit total de fluide en fonction du gradient de pression imposé, il est possible de déterminer la distribution de tailles de pores moyennant une hypothèse sur leur forme. Dans notre approche, un modèle inspiré par Carman-Kozeny a été employé. Cette technique a été validée analytiquement et numériquement sur différents types de distributions classiques (gaussiennes uni-modales ou multi-modales). Elle a été étendue avec succès aux fluides viscoplastiques réels de type Herschel-Bulkley. La seconde méthode proposée repose sur l'analyse dynamique de l'écoulement oscillant d'un fluide newtonien ou non-newtonien à travers le milieu poreux. Elle consiste en l'utilisation de la fonction de transfert hydrodynamique du milieu poreux et plus particulièrement en la caractérisation de son admittance complexe car l'épaisseur de pénétration et donc la taille d'un pore est fonction de la fréquence d'oscillation du gradient de pression imposé. Le modèle de faisceau de capillaires parallèles (de type Carman-Kozeny) a aussi été utilisé. Cette technique a été testée et validée avec succès avec plusieurs types de distributions, dans le cas où le fluide utilisé est newtonien et pour des fluides en loi de puissance. Pour rendre la résolution de ce dernier problème possible nous avons introduit la notion " d'inconsistance complexe ". Par rapport aux techniques existantes, ces deux nouvelles approches se distinguent par leur simplicité, leur non-toxicité et leur faible coût.

Conselho Nacional de Desenvolvimento Científico e Tecnológico / As simulações computacionais tem sido amplamente empregadas no estudo do escoamento darciano e não-darciano em meios porosos consolidados e não-consolidados. Neste trabalho, através de uma nova formulação para a equação de Forchheimer, foram identificadas duas novas propriedades denominados fator de comportamento do fluido, que atua paralelamente a permeabilidade, e permeabilidade equivalente global, resultado da relação anterior. Este comportamento foi estudado e validado através da implementação de um aparato experimental e um código computacional baseado no modelo de regressão-linear que, além disso, demonstrou que o escoamento, ainda que em regime não darciano, comporta-se linearmente como a equação de Darcy, ainda que o coeficiente angular desta diminuiu de acordo com a faixa do número de Reynolds atingida, sendo esta dependente do tipo de leito empregado. Ainda neste trabalho, foi implementado o método de otimização R2W para estimar os parâmetros da equação de Kozeny-Carman a partir de dados experimentais obtidos por Dias et al, a fim de simular o escoamento darciano em meios porosos. Por fim, foi alcançada excelente concordância entre os dados simulados pelo método R2W / equação de Kozeny-Carman e os dados reais. / Computer simulations have been widely used in the study of Darcys flow and non-Darcy porous media in consolidated and non-consolidated. In this work, through a new formulation for the Forchheimer equation, we have been identified two new called Fluid Factor Behavior properties, which acts parallel to permeability, and overall equivalent permeability result of the previous relationship. This behavior has been studied and validated through implementation of an experimental apparatus and a computer code based on the linear regression model, moreover, it demonstrated that flow, even in non darciano system behaves linearly as the Darcy, however, the slope of this decreased according to the range of Reynolds numbers reached, this being dependent on the type of bed used. Although this work was implemented R2W optimization method to estimate the parameters of Kozeny-Carman equation from experimental data provided in the literature in order to simulate the darciano flow in porous media. Finally, it achieved excellent agreement between the data simulated by R2W method / Kozeny-Carman equation and actual data.

Conselho Nacional de Desenvolvimento Científico e Tecnológico / As simulações computacionais tem sido amplamente empregadas no estudo do escoamento darciano e não-darciano em meios porosos consolidados e não-consolidados. Neste trabalho, através de uma nova formulação para a equação de Forchheimer, foram identificadas duas novas propriedades denominados fator de comportamento do fluido, que atua paralelamente a permeabilidade, e permeabilidade equivalente global, resultado da relação anterior. Este comportamento foi estudado e validado através da implementação de um aparato experimental e um código computacional baseado no modelo de regressão-linear que, além disso, demonstrou que o escoamento, ainda que em regime não darciano, comporta-se linearmente como a equação de Darcy, ainda que o coeficiente angular desta diminuiu de acordo com a faixa do número de Reynolds atingida, sendo esta dependente do tipo de leito empregado. Ainda neste trabalho, foi implementado o método de otimização R2W para estimar os parâmetros da equação de Kozeny-Carman a partir de dados experimentais obtidos por Dias et al, a fim de simular o escoamento darciano em meios porosos. Por fim, foi alcançada excelente concordância entre os dados simulados pelo método R2W / equação de Kozeny-Carman e os dados reais. / Computer simulations have been widely used in the study of Darcys flow and non-Darcy porous media in consolidated and non-consolidated. In this work, through a new formulation for the Forchheimer equation, we have been identified two new called Fluid Factor Behavior properties, which acts parallel to permeability, and overall equivalent permeability result of the previous relationship. This behavior has been studied and validated through implementation of an experimental apparatus and a computer code based on the linear regression model, moreover, it demonstrated that flow, even in non darciano system behaves linearly as the Darcy, however, the slope of this decreased according to the range of Reynolds numbers reached, this being dependent on the type of bed used. Although this work was implemented R2W optimization method to estimate the parameters of Kozeny-Carman equation from experimental data provided in the literature in order to simulate the darciano flow in porous media. Finally, it achieved excellent agreement between the data simulated by R2W method / Kozeny-Carman equation and actual data.

Nous développons dans ce travail de thèse un outil de simulation numérique pour les matériaux à changement de phase (MCP), en tenant compte du phénomène de convection naturelle dans la phase liquide, pour des configurations en deux et trois dimensions. Les équations de Navier-Stokes incompressible avec le modèle de Boussinesq pour la prise en compte des forces de flottabilité liées aux effets thermiques, couplées avec une formulation de l’équation d’énergie suivant la méthode d’enthalpie, sont résolues par une méthode d’éléments finis adaptatifs. Une approche mono-domaine, consistant à résoudre les mêmes systèmes d’équations dans les phases solide et liquide, est utilisée. La vitesse est ramenée à zéro dans la phase solide, en introduisant un terme de pénalisation dans l’équation de quantité de mouvement, suivant le modèle de Carman-Kozeny, consistant à freiner la vitesse à travers un milieu poreux. Une discrétisation spatiale des équations utilisant des éléments finis de Taylor-Hood, éléments finis P2 pour la vitesse et éléments finis P1 pour la pression, est appliquée, avec un schéma d’intégration en temps implicite d’ordre deux (GEAR). Le système d’équations non-linéaires est résolu par un algorithme de Newton. Les méthodes numériques sont implémentées avec le logiciel libre FreeFem++ (www.freefem.org), disponible pour tout système d’exploitation. Les programmes sont distribués sous forme de logiciel libre, sous la forme d’une forme de toolbox simple d’utilisation, permettant à l’utilisateur de rajouter d’autres configurations numériques pour des problèmes avecchangement de phase. Nous présentons dans ce manuscrit des cas de validation du code de calcul, en simulant des cas tests bien connus, présentés par ordre de difficulté croissant : convection naturelle de l’air, fusion d’un MCP, le cycle complet fusion-solidification, chauffage par le bas d’un MCP, et enfin, la solidification de l’eau. / In this thesis we develop a numerical simulation tool for computing two and three-dimensional liquid-solid phase-change systems involving natural convection. It consists of solving the incompressible Navier-Stokes equations with Boussinesq approximation for thermal effects combined with an enthalpy-porosity method for the phase-change modeling, using a finite elements method with mesh adaptivity. A single-domain approach is applied by solving the same set of equations over the whole domain. A Carman-Kozeny-type penalty term is added to the momentum equation to bring to zero the velocity in the solid phase through an artificial mushy region. Model equations are discretized using Galerkin triangular finite elements. Piecewise quadratic (P2) finite-elements are used for the velocity and piecewise linear (P1) for the pressure. The coupled system of equations is integrated in time using a second-order Gear scheme. Non-linearities are treated implicitly and the resulting discrete equations are solved using a Newton algorithm. The numerical method is implemented with the finite elements software FreeFem++ (www.freefem.org), available for all existing operating systems. The programs are written and distributed as an easy-to-use open-source toolbox, allowing the user to code new numerical algorithms for similar problems with phase-change. We present several validations, by simulating classical benchmark cases of increasing difficulty: natural convection of air, melting of a phase-change material, a melting-solidification cycle, a basal melting of a phase-change material, and finally, a water freezing case.