No presente trabalho é apresentada a modelagem e solução numérica da convecção mista em cavidade aquecida por baixo com o topo deslizante, preenchida com meio poroso heterogêneo e homogêneo. Na abordagem heterogênea, o domínio do sólido é representado por blocos condutores de calor igualmente espaçados; a fase fluido circunda os blocos, limitada pelas paredes da cavidade. A abordagem homogênea ou poro-contínua é caracterizada através da porosidade e da permeabilidade da cavidade. As equações de conservação da massa, quantidade de movimento e energia são obtidas, adimensionalizadas e generalizadas de modo a representarem tanto o modelo contínuo quanto o poro-contínuo. A solução numérica é obtida através do método dos volumes finitos. As equações são discretizadas via esquema QUICK e é utilizado o algoritmo SIMPLE para o acoplamento pressão - velocidade. Visando o regime laminar, os parâmetros do escoamento são mantidos no intervalo de 102≤Re≤103 e 103≤Ra≤106 tanto para a abordagem heterogênea, quanto para a homogênea. Nas configurações testadas para o modelo contínuo, 9, 16, 36 e 64 blocos são considerados para cada combinação de Re e Ra e a porosidade microscópica é mantida constante φ=0,64 . No modelo poro-contínuo o número de Darcy (Da) é definido em função do número de blocos da cavidade heterogênea e da porosidade φ. Resultados numéricos do estudo comparativo entre a abordagem microscópica e a macroscópica são apresentados. Como resultado, correlações para o Nusselt médio para os modelos contínuo e poro-contínuo são obtidas em função do Ra modificado para cada Re. / In this work is presented mixed convection heat transfer inside a lid-driven cavity heated from below and filled with heterogeneous and homogeneous porous medium. In the heterogeneous approach, the solid domain is represented by heat conductive equally spaced blocks; the fluid phase surrounds the blocks being limited by the cavity walls. The homogeneous or pore-continuum approach is characterized by the cavity porosity and permeability. Generalized mass, momentum and energy conservation equations are obtained in dimensionless form to represent both the continuum and the pore-continuum models. The numerical solution is obtained via the finite volume method. QUICK interpolation scheme is set for numerical treatment of the advection terms and SIMPLE algorithm is applied for pressure-velocity coupling. Aiming the laminar regime, the flow parameters are kept in the range of 102≤Re≤103 and 103≤Ra≤106 for both the heterogeneous and homogeneous approaches. In the tested configurations for the continuous model, 9, 16, 36, and 64 blocks are considered for each combination of Re and Ra being the microscopic porosity set as constant φ=0,64 . For the pore-continuum model the Darcy number (Da) is set according to the number of blocks in the heterogeneous cavity and the φ. Numerical results of the comparative study between the microscopic and macroscopic approaches are presented. As a result, average Nusselt number equations for the continuum and the pore continuum models as a function of Ra and Re are obtained.

Apesar de sua familiaridade e simplicidade, a água apresenta um conjunto propriedades termodinâmicas, dinâmicas e estruturais que são ainda objeto de intensa pesquisa. O aumento da densidade com a temperatura, da difusão com a densidade, ou ainda do ordenamento com a temperatura são exemplos de alguns de seus comportamentos não usuais. Com a finalidade de melhor compreender tais propriedades inúmeras abordagens têm sido utilizadas, tais como o uso geometrias de confinamento, modelos simplificados ou até mesmo misturas. Dentre as geometrias confinantes frequentemente usadas, encontra-se, nanoporos, placas paralelas e meio porosos. Os meios porosos são formados por obstáculos fixos que impõem efeitos de volume excluído adicionais ao sistema. Já no caso de misturas quando elas ocorrem entre líquidos capazes de formar ligações de hidrogênio, o comportamento não usual da água dá origem a um conjunto ainda maior de propriedades anômalas. A mistura água-metanol por exemplo, é munida de um conjunto propriedades de excesso incapazes de serem descritas pelas teorias usuais. São alguns exemplos, o máximo no calor específico e o mínimo no volume e entalpia de excesso. Neste projeto de doutoramento, nós estudamos por simulações numericas o confinamento por meio poroso (desordem queched) e misturas de água com solutos interagentes. O primeiro estudo é realizado usando um modelo 2D tipo-água que é largamente conhecido na literatura. No segundo estágio, estamos a influência de solutos interagentes nas propriedades de modelos em rede e contínuos. Para o modelo em rede, nós desenvolvemos um modelo de soluto e posteriormente uma técnica capaz de simular misturas de modelos em rede a pressão constante. De posse desta técnica estudamos as propriedades de excesso da mistura. / Although the familiarity and simplicity, the water show a set of thermodynamic, dynamics and structural properties which are still subject to intense research. The increase of density as the temperature, of diffusion as the density, or even of ordering with the temperature are examples of some of its unusual behavior. In order to better understand these properties numerous approaches have been used, such as the use of confinement geometries, simplified models, or ever mixtures. Among the confinement geometries used, are those, nanopores, parallel plates and porous media. The porous media are formed by fixed obstacles that impose the additional excluded volume effects to the system. In the case of mixtures, when they occur between liquids able to form hydrogen-bonds, the unusual behavior of water give rise to a set even higher anomalous properties. The water-methanol mixture, for example, has a set of excess properties unable to be described by usual theories. Some examples are the maximum in the specific heat and minimum in excess volume and enthalpy. In this Ph.D. project, we study by numerical simulations, the confinement of water by porous media(or under quenched disorder) and the mixture of water with interacting solutes. The first study is performed using a 2D lattice model which is widely known in the literature. In a second stage, we study the influence of interacting solutes on the properties of lattice and continuous models. For the lattice model, we develop a solute model and a technique to simulate mixtures of lattice models at constant pressure. Using this technique, we study the excess properties of the mixture. For the continuous model we study the influence of a dimeric solute on the TMD of a water-like model and posteriorly we study the excess properties of this type of mixture.

Submitted by José Henrique Krähenbühl Ambiel (zeambiel@hotmail.com) on 2018-09-28T13:24:06Z No. of bitstreams: 1 Dissertação_JoséAmbiel.pdf: 8234688 bytes, checksum: e8355af378aacfa5b31cc3b2d4f77de7 (MD5) / Approved for entry into archive by Lucilene Cordeiro da Silva Messias null (lubiblio@bauru.unesp.br) on 2018-09-28T18:11:05Z (GMT) No. of bitstreams: 1 ambiel_jhk_me_bauru.pdf: 6828894 bytes, checksum: e1a58fe94084b497e2236056afe8f889 (MD5) / Made available in DSpace on 2018-09-28T18:11:05Z (GMT). No. of bitstreams: 1 ambiel_jhk_me_bauru.pdf: 6828894 bytes, checksum: e1a58fe94084b497e2236056afe8f889 (MD5) Previous issue date: 2018-07-24 / Materiais porosos constituem uma grande gama de materiais que podem ser encontrados na natureza ou em forma artificial. Rochas reservatório é um exemplo importante desse tipo de material, sendo o estudo delas a motivação principal desse trabalho. O estudo de rochas reservatório, de onde são extraídos gases e petróleo, consiste em um problema físico no qual os sistemas mecânico e hidráulico são acoplados. O acoplamento ocorre pois as deformações (no sistema mecânico) inuenciam as pressão (no sistema hidráulico), que por sua vez inuenciam as tensões (sistema mecânico). As equações governantes do sistema mecânico são mostradas e as do hidráulico deduzidas. Para a resolução do problema, o Método dos Elementos Finitos (MEF) foi utilizado para ambos os sistemas físicos, logo, as equações governantes são apresentadas em sua forma fraca e, então, aproximada pelo MEF. Numericamente, o acoplamento pode ser tratado de diferentes maneiras, seja considerando um dos sistemas de maneira bem pobre tal como fórmulas empíricas simplistas, seja considerado os sistemas de maneira individual, ou então de maneira completa. Essa última maneira de considerar um acoplamento, o acoplamento total, é formulada, programada e testada nesse trabalho. Para validar a implementação, dois problemas foram analisados: Problema de Terzaghi e Problema Mandel, ambos com solução analítica conhecidas. Os resultados obtidos numericamente comparados aos analíticos indicam que o método totalmente acoplado foi bem implementado, tanto em 2D quanto em 3D. Nesse trabalho também é mostrada a oscilação numérica que há em problemas de acoplamento hidromecânico e uma das formas de amenizá-la. / Porous materials constitute a wide range of materials that can be found in nature and arti cially. Reservoir rock is an important example of this kind of material, which is the main motivation of this work. The study of reservoir rocks, from which gases and oil are extracted, consists of a physical problem in which mechanical and hydraulic systems are coupled. The coupling occurs because the deformations (in the mechanical system) in uence the pressure (in the hydraulic system), which in turn in uence the stresses (mechanical system). The governing equations of the mechanical system are shown and those of the hydraulic system are deduced. To solve the problem, the Finite Element Method (FEM) is used for both physical systems, so the governing equations are presented in their weak form and then approximated according to the FEM. Numerically, the coupling can be handled in di erent ways, either by considering one of the systems in a very poor way by using simplistic empirical formulas, by considering the systems individually, or in a complete manner. The latter one, the fully-coupled treatment, is formulated, programmed and tested in this work. To validate the implementation, two problems has been analyzed: Terzaghi Problem and Mandel Problem, both with known analytical solutions. The comparison between the results obtained numerically and analytically indicates that the fully coupled method has been well implemented in both 2D and 3D cases. The numerical oscillation existing in hydrmechanical coupled problems is also shown and one of the ways to minimize it is presented.

A major problem of our time is the environmental impact and socio-economic consequences of greenhouse gases, caused mainly by burning fossil fuels. In this context, hybrid porous materials known as MOFs have been investigated both experimentally and computationally for CO2 adsorption, among many other applications. The large number of atoms, however, is an obstacle to computationally expensive methods. This work evaluates the performance of semi-empirical quantum methods (usually applied only to treat organic and biological compounds) in the description of CO2 adsorption process for IRMOF-74 series. The AM1, PM3 PM6 and PM7 methods were used in the description of 72 structures, and the impact of MOZYME algorithm in calculation was also evaluated. Chemical and geometrical properties of the system were estimated considering the presence or absence of water in the structure, variations in CO2 concentration in the unit cell, primary and secondary sites occupancy, and different metal ions in the structure. The results were compared with experimental data and computational estimates obtained by DFT, emphasizing the importance of correction for dispersion forces in energy calculations. The presence of open-shell metal ions affect the calculations, but PM6 and PM7 methods are able to reproduce the geometric structure of MOFs, found that the presence of water hinders CO2 adsorption, detected the primary and secondary adsorption sites, providing estimates for the binding energy comparable to that of the most widespread computational methods and in agreement with experimental data. / Um grande problema da atualidade é o impacto ambiental e as consequências socioeconômicas decorrentes da emissão de gases, causada principalmente pela queima de combustíveis fósseis. Nesse contexto, os materiais porosos híbridos conhecidos como MOFs têm sido estudados tanto experimentalmente quanto computacionalmente para adsorção de CO2, dentre várias outras aplicações. O grande número de átomos, no entanto, é um obstáculo para métodos custosos computacionalmente. Este trabalho avalia o desempenho de métodos semiempíricos (geralmente utilizados apenas para tratar compostos orgânicos e biológicos) na descrição do processo de adsorção de CO2 pela série IRMOF-74. Os métodos AM1, PM3, PM6 e PM7 foram utilizados na descrição de até 72 estruturas, avaliando também o impacto do algoritmo MOZYME nos cálculos. Foram estimadas propriedades químicas e geométricas do sistema considerando a presença ou ausência de água na estrutura, variações na quantidade de CO2 na célula unitária, ocupação de sítios primários e secundários, e diferentes íons metálicos na estrutura. Os resultados foram comparados com dados experimentais e com estimativas computacionais obtidas por DFT, enfatizando a importância de correções para forças de dispersão nos cálculos de energia. A presença de íons metálicos com camada incompleta prejudica os cálculos, porém os métodos PM6 e PM7 reproduziram bem a estrutura cristalográfica das MOFs, identificaram que a presença de água dificulta a adsorção de CO2, detectaram os sítios de adsorção primário e secundário, obtendo estimativas para a energia de ligação comparáveis às dos métodos computacionais mais difundidos e em concordância com dados experimentais.

Orientador: Osvair Vidal Trevisan / Tese (doutorado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecanica / Made available in DSpace on 2018-07-28T11:57:08Z (GMT). No. of bitstreams: 1 Moreno_RosangelaBarrosZanoniLopes_D.pdf: 9810356 bytes, checksum: 23584345c809199b204cd503400d8184 (MD5) Previous issue date: 2000 / Resumo: O estudo é focado na remoção eletrocinética de contaminantes de solos. Fenômenos diretos e acoplados em um meio poroso saturado sujeito a gradientes hidráulico, elétrico e quúnico são modelados de forma completa. O trabalho apresenta os modelos teórico, numérico e experimental desenvolvidos. O modelo teórico proposto tem como ponto de partida a Termodinâmica de NãoEquih'brio. Seu desenvolvimento resulta num grupo de equações diferenciais parciais para os parâmetros de estado do sistema: pressão, potencial elétrico, concentração de espécies quúnicas e temperatura. Equações algébricas descrevendo as reações de eletrólise e de formação de água completam o modelo. As parcelas dos fenômenos acoplados foram adimensionalizadas em relação aos fenômenos diretos, possibilitando a análise da contribuição relativa de cada fenômeno nas respectivas equações de transporte, bem como a análise de transferência de escala. O sistema de equações resultante (difusão hidráulica, difusão elétrica, difusão-convecção-retardação do componente i e difusão-convecção-dissipação de calor) foi expresso na forma discreta pelo método de diferenças finitas e resolvido numericamente por um algoritmo em linguagem FORTRAN. O algoritmo foi testado e um critério para minirnizar erros devidos à discretização foi proposto. O potencial de análise do modelo proposto foi verificado para um conjunto de casos hipotéticos. Três sistemas fisicos (30 rnD, 304 rnD e 3040 rnD) foram sujeitos a combinações de gradientes externos de potencial elétrico (O V/cm, 0.1 V/cm e 0.2 V/cm), de pressão (0.01 atmlm, 0.05 atmlm, 0.1 atmlm) e de concentração (NaCl: 5.12 mol/m3 de injeção, 512 mol/m3 de saturação). Perfis e históricos para os parâmetros de estado do sistema, bem como perfis de condutividade e históricos de volume de fluido produzido acumulado e de corrente total foram determinados. A aplicação da técnica de tratamento foi testada, em laboratório, em uma amostra de arenito Botucatu de 7,66 cm de diâmetro e 35 cm de comprimento e permeabilidade na faixa de 400 rnD. Trinta e uma variáveis foram monitoradas em tempo real com a ajuda de um sistema de aquisição automática de dados. A amostra foi submetida a uma seqüência de testes para comparação do processo de remoção de NaCl (512 mol/m / Abstract: The present study focus on the electrokinetic removal of pollutants from soils. Direct and coupled phenomena in saturated porous media under electric, hydraulic and chemical gradients are comprehensively modeled. The approach is done through the development of a theoretical mo dei, a numerical model and an experimental model. The proposed theoretical model is based on NonEquilibrium Thermodynamics. The development results into a set of partial differential equations for the state parameters of the system: pressure, electrical potential, concentration of chemical species and temperature. Algebraic equations describing the electrolysis and water ionization reactions complete the model. The terms for the coupled phenomena are non-dimensionalized with respect to terms for the direct phenomena. Such procedure facilitates the analysis of the relative contribution of each phenomenon in the transport equations, as well as the analysis of scale transfer. The resulting set of equations (hydraulic diffusion, electric diffusion, diffusionconvection-retardation of components i and diffusion-convection-dissipation of heat) was discretised via the finite difference method and numerically solved by a FORTRAN algorithm. The algorithm was tested and a criterion to minimize discretisation errors was proposed. Some hypothetical cases were simulated to confirm the model potential. Three physical systems (30 rnD, 304 mD and 3040 rnD) were submitted to combinations of electric (OV/cm, 0.1 V/cm and 0.2 V/cm), hydraulic (0.01 atm/m, 0.05 atm/m, 0.1 atm/m) and chemical (NaCI: 5.12 mol/m3 injection, 512 mol/m3 saturation) external gradients. Distribution and history profiles were obtained for the state parameters of the system. Conductivity distribution, accumulated volume of produced fluid and total current history were also ca1culated. The application of the treatment was tested in the laboratory for a sample of Botucatu sandstone of 7,66 cm diameter, 35 cm length and permeability in the range of 400 mD. Thirty-one variables were monitored, in real time, by an automatic data acquisition system. In order to compare NaCI (512 mol/m3) removal, injection of a diluted solution (NaCI 5.12 mol/m3) was conducted under different levels of external electric potential DC (O, 5V and 8V) for a constant hydraulic head (-60 cm H20). The restoration of initial concentration conditions and the simultaneous determination of hydraulic permeability of the sample was carrled out by intermediate tests. The results reveal good repeatability. The numeric and experimental results were qualitatively compared and show good agreement with published data / Doutorado / Termica e Fluidos / Doutor em Engenharia Mecânica

This works presents a study about effective properties of porous solids with nonlinear elastic and elastoplastic matrix. For macroscopic mechanics properties evaluation, micromechanics models are used with effective strain concept relative to the modified second method. The porous are assumed as randomly distributed in the matrix, which presents a constitutive law with linear behavior in dilatation and nonlinear in shear. The results are compared with those provided by finite element methods program ABAQUS, assuming porous with spherical geometry for three dimensional solids. Numerical results from ABAQUS were obtained by an implementation of an external subroutine which incorporates at analysis the nonlinear constitutive law. / CAPES - Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / CNPq - Conselho Nacional de Desenvolvimento Científico e Tecnológico / Este trabalho apresenta um estudo sobre as propriedades efetivas de sólidos porosos com matriz elástica não linear e elastoplástica. Na avaliação das propriedades mecânicas macroscópicas empregam-se modelos micromecânicos lineares em conjunto com o conceito de deformação efetiva correspondente ao método secante modificado. Os poros são admitidos como distribuídos randomicamente na matriz, a qual apresenta uma lei constitutiva caracterizada por um comportamento linear em dilatação e não linear em cisalhamento. Os resultados obtidos são confrontados com aqueles fornecidos pelo programa comercial de elementos finitos ABAQUS, admitindo-se que os poros exibem geometrias esféricas para sólidos tridimensionais. A geração dos resultados numéricos oriundos do programa ABAQUS foi viabilizada mediante a implementação de uma sub-rotina externa que incorpora a relação constitutiva não linear considerada nas análises.

Formulamos um modelo simplificado para o estudo do processo de injeção de solvente em reservatórios de petróleo, onde o fluido injetado (um ácido) tem a capacidade de dissolver parcialmente a matriz sólida. Como hipóteses principais, consideramos que o solvente e o soluto (componente químico que constitui o meio poroso) são espécies totalmente miscíveis, a viscosidade da mistura solvente + soluto não varia com a concentração de soluto, há significativa transferência de massa entre as fases e a permeabilidade do meio poroso varia linearmente com a porosidade. O modelo é formado por duas Equações Diferenciais Parciais, uma do tipo Convecção-Difusão a outra é do tipo Convecção-Reação. Para resolução numérica, desenvolvemos uma metodologia que denominamos de EPEC (Explícita Porosidade e Explícita Concentração). Tal metodologia se baseia em um limitador de fluxo do tipo TVD e em diferenças finitas centradas de segunda ordem. Em adição, o EPEC emprega uma técnica de separação de operadores. Deste modo, em cada passo de tempo, realizamos inicialmente o cálculo explícito da porosidade seguido do cálculo explícito da concentração do solvente. Assim, obtemos um desacoplamento natural das equações que descrevem o problema. Resultados de simulações são apresentados para um meio poroso bidimensional, após sessenta dias de injeção de solvente. / We formulate a simplified Model to study the process of solvent injection in petroleum Reservoir, where the injected fluid (an acid) can partially dissolve a solid matrix. As prime hypotheses, we considered that solvent an soluble component are completely mixed, the viscosity of the fluid does not vary with the concentration of the soluble component, theres significant transfer of mass between the parts and, the permeability of media porous changes linearly with porosity. The model is formed by two Partial Differential Equation, one is convection-diffusion type and another is a convection-reaction type. The Numerical Resolution weve developed a method called EPEC (Explicit Porosity Explicit Concentration). Such methodology is based upon a Limiting of Flow of TVD type and, used Centered Finite Differences of second order. In addition, the EPEC use a operators separation technique. This way, every time, first we clearly calculate the porosity and then the concentration of solvent is calculated. Thus we obtain a natural decoupling of the equations that describe the problem. Simulation results are presented to a two dimensional media porous after sixty days of solvent injection.

Formulamos um modelo simplificado para o estudo do processo de injeção de solvente em reservatórios de petróleo, onde o fluido injetado (um ácido) tem a capacidade de dissolver parcialmente a matriz sólida. Como hipóteses principais, consideramos que o solvente e o soluto (componente químico que constitui o meio poroso) são espécies totalmente miscíveis, a viscosidade da mistura solvente + soluto não varia com a concentração de soluto, há significativa transferência de massa entre as fases e a permeabilidade do meio poroso varia linearmente com a porosidade. O modelo é formado por duas Equações Diferenciais Parciais, uma do tipo Convecção-Difusão a outra é do tipo Convecção-Reação. Para resolução numérica, desenvolvemos uma metodologia que denominamos de EPEC (Explícita Porosidade e Explícita Concentração). Tal metodologia se baseia em um limitador de fluxo do tipo TVD e em diferenças finitas centradas de segunda ordem. Em adição, o EPEC emprega uma técnica de separação de operadores. Deste modo, em cada passo de tempo, realizamos inicialmente o cálculo explícito da porosidade seguido do cálculo explícito da concentração do solvente. Assim, obtemos um desacoplamento natural das equações que descrevem o problema. Resultados de simulações são apresentados para um meio poroso bidimensional, após sessenta dias de injeção de solvente. / We formulate a simplified Model to study the process of solvent injection in petroleum Reservoir, where the injected fluid (an acid) can partially dissolve a solid matrix. As prime hypotheses, we considered that solvent an soluble component are completely mixed, the viscosity of the fluid does not vary with the concentration of the soluble component, theres significant transfer of mass between the parts and, the permeability of media porous changes linearly with porosity. The model is formed by two Partial Differential Equation, one is convection-diffusion type and another is a convection-reaction type. The Numerical Resolution weve developed a method called EPEC (Explicit Porosity Explicit Concentration). Such methodology is based upon a Limiting of Flow of TVD type and, used Centered Finite Differences of second order. In addition, the EPEC use a operators separation technique. This way, every time, first we clearly calculate the porosity and then the concentration of solvent is calculated. Thus we obtain a natural decoupling of the equations that describe the problem. Simulation results are presented to a two dimensional media porous after sixty days of solvent injection.

Submitted by Dilene Paulo (dilene.fatima@ufcg.edu.br) on 2018-02-07T13:41:46Z No. of bitstreams: 1 VERA SOLANGE DE OLIVEIRA FARIAS – TESE PPGEP 2011.pdf: 7179434 bytes, checksum: 4a30c9a95f4a089e00fa550fbf1b42b8 (MD5) / Made available in DSpace on 2018-02-07T13:41:46Z (GMT). No. of bitstreams: 1 VERA SOLANGE DE OLIVEIRA FARIAS – TESE PPGEP 2011.pdf: 7179434 bytes, checksum: 4a30c9a95f4a089e00fa550fbf1b42b8 (MD5) Previous issue date: 2011-04-29 / CNPq / Este trabalho apresenta a solução numérica da equação de difusão tridimensional em regime transiente, para um domínio arbitrário. Para atingir os objetivos, a equação de difusão foi discretizada usando coordenadas generalizadas via método dos volumes finitos com uma formulação totalmente implícita, para condições de contorno de equilíbrio e convectiva. Para cada passo no tempo, o sistema de equações obtido para uma dada malha estruturada foi resolvido pelo método de Gauss-Seidel. O código computacional foi desenvolvido em FORTRAN, usando o estúdio CVF 6.6.0, na plataforma Windows Vista. A solução proposta foi validada usando soluções analíticas e numéricas da equação de difusão para várias geometrias, permitindo validar malhas ortogonais e não-ortogonais. A análise e comparação dos resultados mostraram que a solução proposta forneceu resultados coerentes para todos os casos investigados. O código computacional desenvolvido foi aplicado na simulação, a partir de dados experimentais da secagem de telhas cerâmicas para as seguintes condições experimentais: temperaturas de 55,6 °C, 69,7 °C, 82,7 °C e 98,6 °C e teor de umidade inicial variando de 0,2345 até 0,2405 (b.s.). A simulação tornou possível determinar o coeficiente de difusão efetivo em função da razão de umidade e da temperatura do ar de secagem e também o valor do coeficiente de transferência convectivo de massa correspondente para cada temperatura. / This work presents a three-dimensional numerical solution for the diffusion equation in transient state, in an arbitrary domain. The diffusion equation was discretized using the finite volume method with a fully implicit formulation and generalized coordinates, for the equilibrium and convective boundary condition. For each time step, the system of equations obtained for a given structured mesh was solved by the Gauss-Seidel method. A computational code in FORTRAN, using the CFV 6.6.0 Studio, in a Windows Vista platform was developed. The proposed solution was validated by analytical and numerical solutions of the diffusion equation for several geometries. The geometries tested enabled to validate both orthogonal and non-orthogonal meshes. The analysis and comparison of the results showed that the proposed solution provides correct results for all cases investigated. The developed computational code was applied in the simulation, using experimental data of the drying of ceramic roof tiles, for the following experimental conditions: temperature from 55.6; 69.7; 82.7; 72.8 and 98.7 °C, initial moisture content from 0.2345 up to 0.2405 (d.b.). The simulation makes it possible to determine an expression for the diffusion coefficient as a function of the moisture content and temperature of the drying air, and also the value of the convective mass transfer coefficient corresponding to each temperature.

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.